Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale modeling”

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 721 records · Page 40

Monte Carlo Simulations of the Formation Flying Dynamics for the Magnetospheric Multiscale (MMS) Mission

The MMS mission is an ambitious space physics mission that will fly 4 spacecraft in a tetrahedron formation in a series of highly elliptical orbits in order to study magnetic reconnection in the Earth's magnetosphere. The mission design is comprised of a combination of deterministic orbit adjust and random maintenance maneuvers distributed over the 2.5 year mission life. Formal verification of the requirements is achieved by analysis through the use of the End-to-End (ETE) code, which is a modular simulation of the maneuver operations over the entire mission duration. Error models for navigation accuracy (knowledge) and maneuver execution (control) are incorporated to realistically simulate the possible maneuver scenarios that might be realized These error models, coupled with the complex formation flying physics, lead to non-trivial effects that must be taken into account by the ETE automation. Using the ETE code, the MMS Flight Dynamics team was able to demonstrate that the current mission design satisfies the mission requirements.

Schiff, Conrad↗

Multiscale Simulations of Magnetic Island Coalescence

We describe a new interactive parallel Adaptive Mesh Refinement (AMR) framework written in the Python programming language. This new framework, PyAMR, hides the details of parallel AMR data structures and algorithms (e.g., domain decomposition, grid partition, and inter-process communication), allowing the user to focus on the development of algorithms for advancing the solution of a systems of partial differential equations on a single uniform mesh. We demonstrate the use of PyAMR by simulating the pairwise coalescence of magnetic islands using the resistive Hall MHD equations. Techniques for coupling different physics models on different levels of the AMR grid hierarchy are discussed.

Dorelli, John C.↗

Hydration-Controlled Proton Transport in Respiratory Complex I

Proton pumping by respiratory complex I is one essential element for generating the proton motive force that drives ATP synthesis in mitochondria. Although it is understood that electrons from NADH reduce ubiquinone at the peripheral arm and that four protons are transferred in the membrane domain, the mechanism by which this redox reaction initiates proton translocation remains unclear. A lateral pathway linking the quinone binding site to the membrane domain via ND1, ND3, and ND4L subunits has been proposed as a possible initial path of an excess proton. However, experimental structures indicate that the hydration connectivity between D66 ND3 and E34 ND4L is comparatively weaker than in neighboring segments, suggesting a potential regulatory point for proton transfer. Using multiscale reactive molecular dynamics (MS-RMD) and a water wire connectivity metric, we directly simulate proton transport through this region as coupled to the hydration by water molecules. Our results reveal that proton transfer is thermodynamically feasible when transient hydration aligns with the presence of an excess proton, revealing the strong coupling between hydration and proton transfer (PT) in this region of Complex I. These findings support a model where proton injection enhances local hydration, dynamically opening the pathway for proton transfer and regulating the onset of proton pumping in Complex I.

Chemical structure↗

Large-eddy simulation of an atmospheric bore and associated gravity wave effects on wind farm performance in the southern Great Plains

Gravity waves are a common occurrence in the atmosphere, with a variety of generation mechanisms. Their impact on wind farms has only recently gained attention, with most studies focused on wind farm-induced gravity waves. In this study, the interaction between a wind farm and gravity waves generated by an atmospheric bore event is assessed using multiscale large-eddy simulations. The atmospheric bore is created by a thunderstorm downdraft from a nocturnal mesoscale convective system (MCS). The associated gravity waves impact the wind resource and power production at a nearby wind farm during the American Wake Experiment (AWAKEN) in the US southern Great Plains. A two-domain nested setup (Δx=300 and 20 m) is used in the Weather Research and Forecasting (WRF) model, forced with data from the High-Resolution Rapid Refresh model, to capture both the formation of the bore and its interaction with individual wind turbines. The MCS is resolved on the large outer domain, where the structure of the bore and the associated gravity waves are found to be especially sensitive to parameterized microphysics processes. On the finer inner domain, gravity wave interactions with individual wind turbines are resolved; wake dynamics are captured using a generalized actuator disk parameterization in WRF. The gravity waves are found to have a strong effect on the atmosphere above the wind farm; however, the effect of the waves is more nuanced closer to the surface where there is additional turbulence, both ambient and wake-generated. Notably, the gravity waves modulate the mesoscale environment by weakening and dissipating the preexisting low-level jet, which reduces hub-height wind speed and hence the simulated power output, which is confirmed by the observed supervisory control and data acquisition (SCADA) power data. Additionally, the gravity waves induce local wind direction variations correlated with fluctuations in pressure, which lead to fluctuations in the simulated power output as various turbines within the farm are subjected to waking from nearby turbines.

17 WIND ENERGY↗

Multiscale Simulation of Deployable Composite Structures

In this paper, a multiscale simulation method for analyzing deployable composite structures is presented. Effective shell properties of the composites are obtained based on Mechanics of Structure Genome (MSG) homogenization, and then implemented into a user-subroutine UGENS for structural simulation with shell elements in Abaqus. The column bending test(CBT) of a flat thin flexure and lenticular composite boom in a simplified deployer structure are studied for demonstration. A viscoelastic material model with direct integration is adopted in this paper. The CBT simulation shows good agreement with experiments during relaxation, while errors are observed when comparing residual deformation. It is shown that this CBT model can be calibrated to CBT test results. For the lenticular boom analysis, the complete process of flattening, coiling, stowage, deployment and recovery is simulated with the viscoelastic shell model.

Finite element analysis↗

Canopy reflectance modelling of semiarid vegetation

Three different types of remote sensing algorithms for estimating vegetation amount and other land surface biophysical parameters were tested for semiarid environments. These included statistical linear models, the Li-Strahler geometric-optical canopy model, and linear spectral mixture analysis. The two study areas were the National Science Foundation's Jornada Long Term Ecological Research site near Las Cruces, NM, in the northern Chihuahuan desert, and the HAPEX-Sahel site near Niamey, Niger, in West Africa, comprising semiarid rangeland and subtropical crop land. The statistical approach (simple and multiple regression) resulted in high correlations between SPOT satellite spectral reflectance and shrub and grass cover, although these correlations varied with the spatial scale of aggregation of the measurements. The Li-Strahler model produced estimated of shrub size and density for both study sites with large standard errors. In the Jornada, the estimates were accurate enough to be useful for characterizing structural differences among three shrub strata. In Niger, the range of shrub cover and size in short-fallow shrublands is so low that the necessity of spatially distributed estimation of shrub size and density is questionable. Spectral mixture analysis of multiscale, multitemporal, multispectral radiometer data and imagery for Niger showed a positive relationship between fractions of spectral endmembers and surface parameters of interest including soil cover, vegetation cover, and leaf area index.

Franklin, Janet↗

Stochastic Simulation of Mudcrack Damage Formation in an Environmental Barrier Coating

The FEAMAC/CARES program, which integrates finite element analysis (FEA) with the MAC/GMC (Micromechanics Analysis Code with Generalized Method of Cells) and the CARES/Life (Ceramics Analysis and Reliability Evaluation of Structures / Life Prediction) programs, was used to simulate the formation of mudcracks during the cooling of a multilayered environmental barrier coating (EBC) deposited on a silicon carbide substrate. FEAMAC/CARES combines the MAC/GMC multiscale micromechanics analysis capability (primarily developed for composite materials) with the CARES/Life probabilistic multiaxial failure criteria (developed for brittle ceramic materials) and Abaqus (Dassault Systèmes) FEA. In this report, elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio elements enabled the formation of crack boundaries, leading to development of mudcrack-patterned damage. Finite element models of a disk-shaped three-dimensional specimen and a twodimensional model of a through-the-thickness cross section subjected to progressive cooling from 1,300 °C to an ambient temperature of 23 °C were made. Mudcrack damage in the coating resulted from the buildup of residual tensile stresses between the individual material constituents because of thermal expansion mismatches between coating layers and the substrate. A two-parameter Weibull distribution characterized the coating layer stochastic strength response and allowed the effect of the Weibull modulus on the formation of damage and crack segmentation lengths to be studied. The spontaneous initiation of cracking and crack coalescence resulted in progressively smaller mudcrack cells as cooling progressed, consistent with a fractal-behaved fracture pattern. Other failure modes such as delamination, and possibly spallation, could also be reproduced. The physical basis assumed and the heuristic approach employed, which involves a simple stochastic cellular automaton methodology to approximate the crack growth process, are described. The results ultimately show that a selforganizing mudcrack formation can derive from a Weibull distribution that is used to describe the stochastic strength response of the bulk brittle ceramic material layers of an EBC.

Nemeth, Noel N.↗

Alumina – Stabilized SEI and CEI in Potassium Metal Batteries

Aluminum oxide (Al 2 O 3 ) nanopowder is spin-coated onto both sides of commercial polypropene separator to create artificial solid-electrolyte interphase (SEI) and artificial cathode electrolyte interface (CEI) in potassium metal batteries (KMBs). Here, this significantly enhances the stability, including of KMBs with Prussian Blue (PB) cathodes. For example, symmetric cells are stable after 1,000 cycles at 0.5 mA/cm 2 –0.5 mAh/cm 2 and 3.0 mA/cm 2 –0.5 mAh/cm 2 . Alumina modified separators promote electrolyte wetting and increase ionic conductivity (0.59 vs. 0.2 mS/cm) and transference number (0.81 vs. 0.23). Cryo-stage focused ion beam (cryo-FIB) analysis of cycled modified anode demonstrates dense and planar electrodeposits, versus unmodified baseline consisting of metal filaments (dendrites) interspersed with pores and SEI. Alumina-modified CEI also suppresses elemental Fe crossover and reduces cathode cracking. Mesoscale modeling of metal – SEI interactions captures crucial role of intrinsic heterogeneities, illustrating how artificial SEI affects reaction current distribution, conductivity and morphological stability.

25 ENERGY STORAGE↗

Data mining and computational screening of Rashba-Dresselhaus splitting and optoelectronic properties in two-dimensional perovskite materials

Recent developments highlighting the promise of two-dimensional perovskites have vastly increased the compositional search space in the perovskite family. This presents a great opportunity for the realization of highly performant devices and practical challenges associated with the identification of candidate materials. High-fidelity computational screening offers great value in this regard. In this study, we carry out a multiscale computational workflow, generating a dataset of two-dimensional perovskites in the Dion-Jacobson and Ruddlesden-Popper phases. Our dataset comprises ten B-site cations, four halogens, and over 20 organic cations across over 2000 materials. We compute electronic properties, thermoelectric performance, and numerous geometric characteristics. Furthermore, we introduce a framework for the high-throughput computation of Rashba-Dresselhaus splitting. Finally, we use this dataset to train machine learning models for the accurate prediction of band gaps, candidate Rashba-Dresselhaus materials, and partial charges. The work presented herein can aid future investigations of two-dimensional perovskites with targeted applications in mind.

14 SOLAR ENERGY↗

Multiscale Computer Simulation of Failure in Aerogels

Aerogels have been of interest to the aerospace community primarily for their thermal properties, notably their low thermal conductivities. While such gels are typically fragile, recent advances in the application of conformal polymer layers to these gels has made them potentially useful as lightweight structural materials as well. We have previously performed computer simulations of aerogel thermal conductivity and tensile and compressive failure, with results that are in qualitative, and sometimes quantitative, agreement with experiment. However, recent experiments in our laboratory suggest that gels having similar densities may exhibit substantially different properties. In this work, we extend our original diffusion limited cluster aggregation (DLCA) model for gel structure to incorporate additional variation in DLCA simulation parameters, with the aim of producing DLCA clusters of similar densities that nevertheless have different fractal dimension and secondary particle coordination. We perform particle statics simulations of gel strain on these clusters, and consider the effects of differing DLCA simulation conditions, and the resultant differences in fractal dimension and coordination, on gel strain properties.

Good, Brian S.↗

Characterizing multiscale variability of zero intermittency in spatial rainfall

In this paper the authors study how zero intermittency in spatial rainfall, as described by the fraction of area covered by rainfall, changes with spatial scale of rainfall measurement or representation. A statistical measure of intermittency that describes the size distribution of 'voids' (nonrainy areas imbedded inside rainy areas) as a function of scale is also introduced. Morphological algorithms are proposed for reconstructing rainfall intermittency at fine scales given the intermittency at coarser scales. These algorithms are envisioned to be useful in hydroclimatological studies where the rainfall spatial variability at the subgrid scale needs to be reconstructed from the results of synoptic- or mesoscale meteorological numerical models. The developed methodologies are demsonstrated and tested using data from a severe springtime midlatitude squall line and a mild midlatitude winter storm monitored by a meteorological radar in Norman, Oklahoma.

Kumar, Praveen↗

Multiscale investigation of thermomechanical and compositional developments in Ni alloy 718 under laser processing

Laser processing has been widely employed in various applications due to its exceptional spatial resolution. However, the rapid temperature gradients generated in localized areas present significant challenges for experimental characterization using conventional instruments. To characterize Ni alloy 718 during laser processing, we employed in-situ synchrotron X-ray diffraction with a high-speed detector, a method particularly well-suited for probing processes with high temporal and spatial resolution. Through a series of in-situ experiments, we investigated the local variations in the evolution of microstructures and thermomechanical behaviors within a keyhole mode melt pool. The in situ macroscopic thermomechanical behaviors were quantified using an empirical model derived from diffraction patterns, with experimental results showing reasonable agreement with finite element analysis. Various laser parameters were tested to assess their influences on the residual strains in the melt pools. The results revealed that the residual strain in the keyhole mode melt pool is relatively insensitive to variations in the parameters and is smaller than that in the melt pool created under conduction mode laser scanning. Additionally, we analyzed the shapes of individual diffraction spots, providing insights into the plastic behaviors and compositional developments in the resolidified alloy. The analysis confirmed that compositional variations in a dendritic microstructure manifest as asymmetric broadening of the diffraction spots.

Ni alloy 718↗

MMS Observations of Energized He+ Pickup Ions at Quasiperpendicular Shocks

Accelerated He+ PUIs, downstream of quasi-perpendicular shocks, are studied as 10a function of the fast mode Mach number (Mf) and shock obliquity (𝜃Bn).We analyze quasi-perpendicular shocks with Mach numbers in the range [1, 7] observed by the Magnetospheric MultiScale (MMS) mission, and compare upstream and downstream He+ velocity distribution functions (VDFs). For each shock event, we characterize the upstream PUI distribution and derive reduced 1D velocity distributions for the selected upstream and downstream intervals. We also compare the upstream to downstream ratio of spectral indices, computed from the He+ perpendicular distributions, to Mf and 𝜃Bn. We find a positive correlation between this spectral index ratio and Mf, which suggests that perpendicular energization of He+ PUIs is enhanced as the shock becomes stronger. These results inform modeling efforts of PUIs and shock acceleration processes, particularly those taking place at the termination shock.

Michael J. Starkey↗

Structural characterization of high-protein system through ultra-small and small-angle X-ray scattering

High-protein systems exhibit a hierarchical structure consisting of interconnected multiscale length assemblies. Recent investigations are directed toward identifying these structural units under various processing and environmental conditions to establish structure–function correlations. Ultra-small and small-angle X-ray scattering (USAXS/SAXS) has become a crucial tool for characterizing the structures of proteins and their clusters/aggregates, ranging from nanometers to micrometers, with minimal disruption to their original state. Here, this review first describes the facilities and principles of X-ray scattering, followed by discussions on the analysis of scattering data, including the interpretation of fitting models. It then delves into the main applications of USAXS/SAXS in plant and dairy protein-rich systems. Future strategies to enhance the utilization of scattering techniques for elucidating the structure of high-protein systems are also included.

60 APPLIED LIFE SCIENCES↗

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. Ab initio methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system’s bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision. RAPTOR may also be interfaced with enhanced sampling methods to drive simulations toward the analysis of reactive rare events, and a number of collective variables (CVs) have been developed to facilitate this. Key advances to this methodology, including GPU acceleration efforts and novel CVs to model water wire formation are reviewed, along with recent applications of the method which demonstrate its versatility and robustness.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Arctic Sea Ice Motion from Wavelet Analysis of Satellite Data

Wavelet analysis of DMSP SSM/I (Special Sensor Microwave/Imager) 85 GHz and 37 GHz radiance data, SMMR (Scanning Multichannel Microwave Radiometer) 37 GHz, and NSCAT (NASA Scatterometer) 13.9 GHZ data can be used to obtain daily sea ice drift information for both the northern and southern polar regions. The derived maps of sea ice drift provide both improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing data from satellite synthetic aperture radars (SAR). Examples of derived ice-drift maps in the Arctic illustrate large-scale circulation reversals within a period of a couple weeks. Comparisons with ice displacements derived from buoys show good quantitative agreement. NASA Scatterometer (NSCAT) 13.9 GHZ data have been also used for wavelet analysis to derive sea-ice drift. First, the 40' incidence-angle, sigma-zero (surface roughness) daily map of whole Arctic region with 25 km of pixel size from satellite's 600 km swath has been constructed. Then, the similar wavelet transform procedure to SSM/I data can be applied. Various scales of wavelet transform and threshold have been tested. By overlaying , neighbor filtering, and block-averaging the results of multiscale wavelet transforms, the final sea ice drift vectors are much smooth and representative to the sea ice motion. This wavelet analysis procedure is robust and can make a major contribution to the understanding of ice motion over large areas at relatively high temporal resolutions. The results of wavelet analysis of SSM/I and NSCAT images and buoy data can be merged by some data fusion techniques and will help to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model.

Liu, Antony K.↗

NavCube3-mini Lunar GNSS Receiver

This paper describes development and testing of NASA Goddard Space Flight Center’s new NavCube3-mini (NC3m) spaceborne, weak-signal GNSS receiver, which targets all Earth orbit regimes with special focus on lunar applications. NC3m derives from the ground-breaking Magnetospheric Multiscale (MMS) mission Navigator GPS receiver. The MMS-Navigator (launched 2015) holds a Guinness World Record for highest altitude GPS fix and is currently in a highly elliptic orbit with a 29 Earth radii apogee, nearly half lunar distance. NC3m has reduced size, weight, and power compared to the MMS-Navigator, making it suitable for smallsat applications, and adds multi-frequency and multi-GNSS capabilities, among other improvements. A NC3m engineering test unit was subjected to and successfully completed a comprehensive Technology Readiness Level 6 (system/subsystem model or prototype demonstration in a relevant environment) testing campaign in the second half of 2022. This involved high-fidelity simulations in LEO, GEO, and lunar regimes and full environmental testing, including vibration, thermal vacuum, and electromagnetic compatibility. An overview of the NC3m receiver, the test setup, and a sample of results is presented. The results include predicted performance at high-altitude and in the lunar regime.

NavCube3↗

Towards a multiscale approach for understanding irradiation induced swelling and creep in 316 stainless steels - A coupled cluster dynamics and crystal plasticity approach

Structural materials undergo mechanical degradation, in part due to irradiation-induced swelling and creep, under nuclear reactor conditions. While swelling results from the migration and clustering of irradiation-induced atomic-scale mobile defects, the interaction of mesoscale dislocations with these defects causes creep deformation. A coupled crystal plasticity (CP) and mean-field cluster dynamics (CD) approach is presented to investigate the effect of irradiation on the long-term mechanical behavior of 316 stainless steel, which are under consideration for use in nuclear reactors. The temporal evolution of Frenkel pairs and extended defect population, under a chosen irradiation flux and temperature, is predicted using the CD model. The impact of the irradiation defects on the stress state, and the resulting dislocation-mediated inelastic deformation, is modeled concurrently with the CP model. The inelastic deformation is irradiation flux dependent, and early-stage defect evolution determines the later-stage mechanical behavior in 316 stainless steel.

36 - MATERIALS SCIENCE↗