Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale simulations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

Multiscale Multifunctional Progressive Fracture of Composite Structures

A new approach is described for evaluating fracture in composite structures. This approach is independent of classical fracture mechanics parameters like fracture toughness. It relies on computational simulation and is programmed in a stand-alone integrated computer code. It is multiscale, multifunctional because it includes composite mechanics for the composite behavior and finite element analysis for predicting the structural response. It contains seven modules; layered composite mechanics (micro, macro, laminate), finite element, updating scheme, local fracture, global fracture, stress based failure modes, and fracture progression. The computer code is called CODSTRAN (Composite Durability Structural ANalysis). It is used in the present paper to evaluate the global fracture of four composite shell problems and one composite built-up structure. Results show that the composite shells. Global fracture is enhanced when internal pressure is combined with shear loads. The old reference denotes that nothing has been added to this comprehensive report since then.

Chamis, C. C.↗

Quantifying the Agyrotropy of Proton and Electron Heating in Turbulent Plasmas

An important aspect of energy dissipation in weakly collisional plasmas is that of energy partitioning between different species (e.g., protons and electrons) and between different energy channels. Here we analyse pressure–strain interaction to quantify the fractions of isotropic compressive, gyrotropic, and nongyrotropic heating for each species. An analysis of kinetic turbulence simulations is compared and contrasted with corresponding observational results from Magnetospheric Multiscale Mission data in the magnetosheath. In assessing how protons and electrons respond to different ingredients of the pressure–strain interaction, we find that compressive heating is stronger than incompressive heating in the magnetosheath for both electrons and protons, while incompressive heating is stronger in kinetic plasma turbulence simulations. Concerning incompressive heating, the gyrotropic contribution for electrons is dominant over the nongyrotropic contribution, while for protons nongyrotropic heating is enhanced in both simulations and observations. Variations with plasma β are also discussed, and protons tend to gain more heating with increasing β.

Yan Yang↗

Maneuver Recovery Analysis for the Magnetospheric Multiscale Mission

The use of spacecraft formations creates new and more demanding requirements for orbit determination accuracy. In addition to absolute navigation requirements, there are typically relative navigation requirements that are based on the size or shape of the formation. The difficulty in meeting these requirements is related to the relative dynamics of the spacecraft orbits and the frequency of the formation maintenance maneuvers. This paper examines the effects of bi-weekly formation maintenance maneuvers on the absolute and relative orbit determination accuracy for the four-spacecraft Magnetospheric Multiscale (MMS) formation. Results are presented from high fidelity simulations that include the effects of realistic orbit determination errors in the maneuver planning process. Solutions are determined using a high accuracy extended Kalman filter designed for onboard navigation. Three different solutions are examined, considering the effects of process noise and measurement rate on the solutions.

Volle, Michael↗

Maneuver Recovery Analysis for the Magnetospheric Multiscale Mission

The use of spacecraft formations creates new and more demanding requirements for orbit determination accuracy. In addition to absolute navigation requirements, there are typically relative navigation requirements that are based on the size or shape of the formation. The difficulty in meeting these requirements is related to the relative dynamics of the spacecraft orbits and the frequency of the formation maintenance maneuvers. This paper examines the effects of bi-weekly formation maintenance maneuvers on the absolute and relative orbit determination accuracy for the four-spacecraft Magnetospheric Multiscale (MMS) formation. Results are presented from high fidelity simulations that include the effects of realistic orbit determination errors in the maneuver planning process. Solutions are determined using a high accuracy extended Kalman filter designed for onboard navigation. Three different solutions are examined, considering the effects of process noise and measurement rate on the solutions.

Gramling, Cheryl↗

Impacts of Estimated Plume Rise on PM 2.5 Exceedance Prediction During Extreme Wildfire Events: A Comparison of Three Schemes (Briggs, Freitas, and Sofiev)

Plume height plays a vital role in wildfire smoke dispersion and the subsequent effects on air quality and human health. In this study, we assess the impact of different plume rise schemes on predicting the dispersion of wildfire air pollution and the exceedances of the National Ambient Air Quality Standards (NAAQS) for fine particulate matter (PM 2.5 ) during the 2020 western United States wildfire season. Three widely used plume rise schemes (Briggs, 1969; Freitas et al., 2007; Sofiev et al., 2012) are compared within the Community Multiscale Air Quality (CMAQ) modeling framework. The plume heights simulated by these schemes are comparable to the aerosol height observed by the Multi-angle Imaging SpectroRadiometer (MISR) and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO). The performance of the simulations with these schemes varies by fire case and weather conditions. On average, simulations with higher plume injection heights predict lower aerosol optical depth (AOD) and surface PM 2.5 concentrations near the source region but higher AOD and PM 2.5 in downwind regions due to the faster spread of the smoke plume once ejected. The 2-month mean AOD difference caused by different plume rise schemes is approximately 20 %–30 % near the source regions and 5 %–10 % in the downwind regions. Thick smoke blocks sunlight and suppresses photochemical reactions in areas with high AOD. The surface PM 2.5 difference reaches 70 % on the West Coast of the USA, and the difference is lower than 15 % in the downwind regions. Moreover, the plume injection height affects pollution exceedance (>35 µg m−3) predictions. Higher plume heights generally produce larger downwind PM 2.5 exceedance areas. The PM 2.5 exceedance areas predicted by the three schemes largely overlap, suggesting that all schemes perform similarly during large wildfire events when the predicted concentrations are well above the exceedance threshold. At the edges of the smoke plumes, however, there are noticeable differences in the PM 2.5 concentration and predicted PM 2.5 exceedance region. For the whole period of study, the difference in the total number of exceedance days could be as large as 20 d in northern California and 4 d in the downwind regions. This disagreement among the PM 2.5 exceedance forecasts may affect key decision-making regarding early warning of extreme air pollution episodes at local levels during large wildfire events.

Yunyao Li↗

Magnetospheric Multiscale Observations of Earth’s Oblique Bow Shock Reformation by Foreshock Ultra-Low Frequency Waves

Collisionless shocks can be nonstationary with periodic reformation shown in many simulation results, but direct observations are still tenuous and difficult to conclusively interpret. In this study, using Magnetospheric Multiscale (MMS) observations, we report direct observational evidence of Earth’s oblique bow shock reformation driven by the foreshock ultra-low frequency (ULF)waves. When the four MMS spacecraft were in a string-of-pearls formation roughly along the bow shock normal, they observed that when each period of foreshock ULF waves encountered the bow shock, a new shock ramp formed. Meanwhile, in the magnetosheath, the old bow shock’s remnants were observed periodically convecting downstream. We propose that the reformation mechanism of the oblique bow shock is the variation of the upstream conditions by the periodic ULF waves as they encounter the bow shock. We also examine the nature of reflected ions during the reformation process.

Terry Z Liu↗

Multiscale Modeling of Damage Processes in fcc Aluminum: From Atoms to Grains

Molecular dynamics (MD) methods are opening new opportunities for simulating the fundamental processes of material behavior at the atomistic level. However, current analysis is limited to small domains and increasing the size of the MD domain quickly presents intractable computational demands. A preferred approach to surmount this computational limitation has been to combine continuum mechanics-based modeling procedures, such as the finite element method (FEM), with MD analyses thereby reducing the region of atomic scale refinement. Such multiscale modeling strategies can be divided into two broad classifications: concurrent multiscale methods that directly incorporate an atomistic domain within a continuum domain and sequential multiscale methods that extract an averaged response from the atomistic simulation for later use as a constitutive model in a continuum analysis.

Glaessgen, E. H.↗

Reconnection in a Pinch

A recently published analysis of current sheets has updated the classic Harris 1D static solution by considering multiple classes of charged particle trajectories in a generalized and dynamic current sheet. It uses a 1D PIC simulation to describe dynamic pinching and bifurcation of the current sheet. These 1D results strongly suggest that properties of the inflowing plasma, including the plasma beta, have an important effect on the equilibrium thickness of the pinched current sheet. Since 1D studies cannot describe magnetic reconnection, the time appears right to carry such 1D studies over to 2D or 3D simulations to explore current sheet thickness effects on reconnection. The Magnetospheric Multiscale Mission (MMS), with its well-resolved multipoint measurements of collisionless plasma and fields, has found that collisionless reconnection is accompanied by non-adiabatic motions of electrons that only occur in magnetic structures with a narrow scale comparable to electron inertial lengths (d e ). The recent 1D studies suggest that a plasma pinch to such scales may only occur for inflowing magnetized plasmas with relatively low plasma beta. We conclude that a parametric exploration of simulated and observed reconnection inflow conditions, particularly plasma beta, should shed light on the enablement of reconnection in collisionless plasmas.

reconnection↗

Thermoplastic Matrix Composite Design for Cryotanks Using Multiscale Modeling and Bayesian Optimization

Designing lightweight, robust cryogenic storage tanks is critical for future launch vehicles, in-space propellant storage, and hydrogen powered aircraft. This work presents a multiscale modeling and Bayesian optimization framework for the design of thermoplastic matrix composite cryotanks. Molecular dynamics simulations are first used to determine temperature-dependent constituent properties for candidate thermoplastic matrices, which are homogenized to the lamina scale using NASA’s Multiscale Analysis Tool (NASMAT). These lamina properties, in combination with laminate family generation rules, are evaluated in HyperX structural optimization software to identify stacking sequences that meet all cryogenic load requirements. A Bayesian optimization framework is applied, with HyperX in the loop (via the HyperX API) to efficiently search across material and laminate design variables, yielding an optimized cryotank configuration with significant reductions in design cycle time compared to exhaustive search approaches.

thermoplastics↗

Multiscale Studies of the Three-Dimensional Dayside X-Line

We review recent experience from the Cluster, Double Star, and THEMIS missions for lessons that apply to the upcoming Magnetospheric Multiscale Mission (MMS) being developed for launch in 2014. On global scales, simulation and statistical studies lead to mean congurations of dayside reconnection, implying specific relative alignments of the inflow magnetic fields and X-line, with implications for MMS operations designed to maximize the number of close encounters with the diffusion region. At intermediate MHD-to-ion scales, reconstruction of features created by one or two X-lines have developed to the point where data from a cluster of spacecraft can determine their temporal trends and the approximate three-dimensional X-line structure. Recent petascale particle-in-cell (PIC) simulations of reconnection encompass three spatial dimensions with excellent resolution, and make striking predictions of electron scale physics that creates complex interacting flux ropes under component reconnection. High time resolution measurements from MMS will determine the detailed electron scale kinetics embedded within the global and MHDion scale contexts. These developments will lead to the refinement of our three-dimensional multiscale picture of reconnection, yielding improved understanding of the global, MHD, and local physics controlling the onset or quenching, variability, and mean rate of reconnection. This in turn will enable improved predictability of the structural features created by transient reconnection, and their space weather consequences.

Reconnection↗

Magnetospheric Multiscale Observations of an Expanding Oxygen Wave in Magnetic Reconnection

Heavier plasma species such as oxygen ions can have a large impact on the magnetic reconnection process. It has been hypothesized that the acceleration of demagnetized oxygen ions by the Hall electric field will lead to the formation of an oxygen wave that expands into the exhaust. By comparing data from NASA's Magnetospheric Multiscale mission to a fully kinetic particle-in-cell simulation, we can for the first time provide observational evidence of such an expanding oxygen wave. The wave is characterized by an oxygen jet consisting of cold ions directed toward the neutral sheet associated with a density cavity. This density cavity forms as the EO are subject to collective acceleration by the Hall electric field leaving behind a region of low-density oxygen ions. Our results are important for the understanding of the role and effect of oxygen ions in magnetic reconnection.

Hakon Midthun Kolsto↗

Modeling the Stochastic Response of Fiber Reinforced Composites with Varied Representative Volume Element Sizes

Fiber reinforced composites are desirable in applications where low weight and high strength are needed, but are susceptible to variability and flaws during manufacturing, making failure predictions difficult. These flaws may occur at the microscale where mechanical properties vary locally due to regions of fiber clusters and matrix pockets. In this study, a multiscale approach was taken to model 3-point bend, 4-point bend, and tensile experiments of a unidirectional composite from only having microstructure scans of these samples and constituent properties from literature. These scans were sampled with different sized windows, and statistically equivalent microstructures were generated, then simulated for stiffness, strength, and fracture toughness using a reduced order micromechanical model and NASA’s Multiscale Analysis Tool (NASMAT). Mesoscale models were created with equivalent element sizes to microstructures and properties sampled from microscale simulation results. Results showed how microscale size affects certain mechanical properties. Also shown is how well mesoscale models agree to experiments when using stochastic element properties and varying element size.

microstructure↗

Unstructured LES of Reacting Multiphase Flows in Realistic Gas Turbine Combustors

As part of the Accelerated Strategic Computing Initiative (ASCI) program, an accurate and robust simulation tool is being developed to perform high-fidelity LES studies of multiphase, multiscale turbulent reacting flows in aircraft gas turbine combustor configurations using hybrid unstructured grids. In the combustor, pressurized gas from the upstream compressor is reacted with atomized liquid fuel to produce the combustion products that drive the downstream turbine. The Large Eddy Simulation (LES) approach is used to simulate the combustor because of its demonstrated superiority over RANS in predicting turbulent mixing, which is central to combustion. This paper summarizes the accomplishments of the combustor group over the past year, concentrating mainly on the two major milestones achieved this year: 1) Large scale simulation: A major rewrite and redesign of the flagship unstructured LES code has allowed the group to perform large eddy simulations of the complete combustor geometry (all 18 injectors) with over 100 million control volumes; 2) Multi-physics simulation in complex geometry: The first multi-physics simulations including fuel spray breakup, coalescence, evaporation, and combustion are now being performed in a single periodic sector (1/18th) of an actual Pratt & Whitney combustor geometry.

Ham, Frank↗

Design of a Variational Multiscale Method for Turbulent Compressible Flows

A spectral-element framework is presented for the simulation of subsonic compressible high-Reynolds-number flows. The focus of the work is maximizing the efficiency of the computational schemes to enable unsteady simulations with a large number of spatial and temporal degrees of freedom. A collocation scheme is combined with optimized computational kernels to provide a residual evaluation with computational cost independent of order of accuracy up to 16th order. The optimized residual routines are used to develop a low-memory implicit scheme based on a matrix-free Newton-Krylov method. A preconditioner based on the finite-difference diagonalized ADI scheme is developed which maintains the low memory of the matrix-free implicit solver, while providing improved convergence properties. Emphasis on low memory usage throughout the solver development is leveraged to implement a coupled space-time DG solver which may offer further efficiency gains through adaptivity in both space and time.

Design↗

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

A multiscale repeating unit cell model of a single spherulite containing four disparate length scales was developed to predict the thermoelastic behavior of semicrystalline thermoplastic materials for composite aerospace applications. The continuum level scales were fully coupled and modeled using the generalized method of cells and the high fidelity generalized method of cells micromechanics theories. Data from molecular dynamics simulations were used as inputs for the amorphous and crystalline constituents in the multiscale continuum models. Effective Young’s modulus, shear modulus, Poisson’s ratio, coefficient of thermal expansion, and thermal conductivity were predicted for polyether ether ketone and polyether ketone ketone, showing good agreement with the available experimental data from the open literature. Moreover, it is shown that predicted properties are fairly insensitive to the fidelity of the micromechanics model used at the highest continuum scale or the assumed shape of the spherulite.

thermoplastics↗

2018-2019 Mid-Term Credibility Plan Review

This presentation summarizes the Interagency Modeling and Analysis Group (IMAG) / Multiscale Modeling (MSM)'s Committee on Credible Practice in Modeling and Simulation in Healthcare's (CPMS) mid-term review of each MSM U01 project's progress toward fulfilling their Model Credibility plan. Generally, the model credibility reporting and update process utilized several pieces of information formatted to communicate to the uninformed reviewer/user an assessment of steps taken toward model credibility with respect to the intended use of each U01 model, as well as communicate where the more interested reviewer/user can go to find more information. The reviewers, consisting of the members of the CPMS Executive Committee, focused on the documentation and communication of each PI's credibility efforts and did not perform an assessment of the implemented credible practice given the evolving nature of the projects. In this inaugural activity, the CPMS members sought insight on how PIs may view the importance of communicating credibility, as well as how PIs interpret establishing evidence for credibility within the scope of CPMS Tens simple rules. To that end, reviewers provided feedback on documentation and communication adequacy, and, only in few rare cases, on the credible practice itself.

Myers, Jerry↗

Expanding Hardware-in-the-Loop Formation Navigation and Control with Radio Frequency Crosslink Ranging

The Formation Flying Testbed (FFTB) at the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC) provides a hardware-in-the-loop test environment for formation navigation and control. The facility continues to evolve as a modular, hybrid, dynamic simulation facility for end-to-end guidance, navigation, and control (GN&C) design and analysis of formation flying spacecraft. The core capabilities of the FFTB, as a platform for testing critical hardware and software algorithms in-the-loop, are reviewed with a focus on recent improvements. With the most recent improvement, in support of Technology Readiness Level (TRL) 6 testing of the Inter-spacecraft Ranging and Alarm System (IRAS) for the Magnetospheric Multiscale (MMS) mission, the FFTB has significantly expanded its ability to perform realistic simulations that require Radio Frequency (RF) ranging sensors for relative navigation with the Path Emulator for RF Signals (PERFS). The PERFS, currently under development at NASA GSFC, modulates RF signals exchanged between spacecraft. The RF signals are modified to accurately reflect the dynamic environment through which they travel, including the effects of medium, moving platforms, and radiated power.

Mitchell, Jason W.↗

A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Yamakov, Vesselin↗