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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

VOTCA: multiscale frameworks for quantum and classical simulations in soft matter

Many physical phenomena in liquids and soft matter are multiscale by nature and can involve processes with quantum and classical degrees of freedom occurring over a vast range of length- and timescales. Examples range from structure formation processes of complex polymers or even polymer blends (Svaneborg & Everaers, 2023) on the classical side to charge and energy transport and conversion processes (Lee et al., 2019) involving explicit electronic and, therefore, quantum information. The Versatile Object-oriented Toolkit for Coarse-graining Applications (VOTCA) provides multiscale frameworks built on a comprehensive set of methods for the development of classical coarse-grained potentials (VOTCA-CSG) as well as state-of-the art excited state electronic structure methods based on density-functional and many-body Green’s function theories, coupled in mixed quantum-classical models and used in kinetic network models (VOTCA-XTP).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Videos, photos, and AI-derived grain size data associated with “High-throughput AI Video Surveys Enable Reproducible Multiscale Sediment Size Mapping, with Implications for Hydrobiogeochemical Parameterization”

NOTE: The manuscript associated with this data package is currently in review. The data may be revised based on reviewer feedback. Upon manuscript acceptance, this data package will be updated with the final dataset and additional metadata. This data package is associated with the manuscript “High-throughput AI Video Surveys Enable Reproducible Multiscale Sediment Size Mapping, with Implications for Hydrobiogeochemical Parameterization” under review. This data package includes five data types: 1) raw photos and videos from drone survey and walking smartphone surveys; 2) images derived from raw videos; 3) manual labeling of reference scales; 4) metadata for all images and photo resolution derived from artificial intelligence (AI) models or manual labels, 5) grain size data obtained from AI models for all photos, 6) metadata and grain size data after quality control, 7) summaries of sample efficiency for all data, and 8) computational fluid dynamics (CFD) data used to support hydro-biogeochemical (HBGC) parameter estimation. Such data is used to 1) demonstrate significant improvements in accuracy, efficiency, and quality control for grain size data collection with the help of AI models, 2) study the spatial heterogeneity of grain size and observation reproducibility based on tens of thousands of data points generated by the AI models, and 3) evaluate the impacts of grain size heterogeneity on key HBGC parameters across sediment-to-reach and hourly-to-yearly scales. In particular, the data package contains 116 folders and 179696 files. The files include 41 videos in .mov format, 64047 photos in .jpg format, 13541 video-derived photos in .png format, 12747 segmentation mask data in .tif format, 12747 segmentation data in .json format, 24771 .csv files that with metadata and grain size for each individual photo as well as water depth and velocity data from CFD and observation, 51791 .txt files of raw AI predicted labels, and 11 flight record data in .srt format. The summary for all metadata and grain size statistics information is included in “Scales_V3_NG.csv” and “Statistics_V3_NG.csv”. The summary for data that pass data quality control (QC) level 0-2 is included in “QCStatistics_V3_NG.csv”. The QC level 0 represents photos whose photo resolution is positive, excluding photos that miss reference scale. The QC level 1 means reference scale circularity uncertainty is less than 5% for smartphone images while representing photo resolution is larger than 0.44 mm/pixel for drone images. The QC level 2 means excluding photos whose grain number is less than 100, a minimum number of grains recommended by classic literature. The summary for each video’s name, length, frame rates, survey area, grain number, survey efficiency, etc. can be found in “QCSummary_V3_NG.csv”. The summary for site name, GPS coordinates, and number of images at each site can be found in “SitesSummary_V3_*.csv” files. Overall computational efficiency summary is reported in Table 4 of accompanying manuscript. Additionally, the nitrate concentration data used in this work was downloaded from an existing dataset published on ESS-DIVE (Boat-Dragged Sensor Hanford Reach.csv; Conner A. et al., 2020). We thank the United States Forest Service, Washington Department of Fish and Wildlife, Washington Department of Natural Resources, Cowiche Canyon Conservatory, Port of Benton, and the Confederated Tribes and Bands of the Yakama Nation for access to field locations where the data were collected. We also thank the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate data collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Dispatch Manager for NEML2 Constitutive Model Calculations Embedded in MOOSE

This report describes the extended capabilities of the NEML2 constitutive modeling library, including a flexible and efficient work dispatching system designed to leverage both CPU and GPU resources. This enhancement addresses one of the primary computational challenges in large-scale simulations: the ability to distribute and execute batches of material model evaluations across heterogeneous computing devices. The new dispatch system introduces a modular set of dispatcher and scheduler classes that coordinate the flow of data and execution between devices. The dispatcher is responsible for efficiently packaging work, managing device-specific memory operations, and synchronizing results. This modularity allows for extensibility, making it straightforward to integrate additional computing backends in the future. From an implementation standpoint, the dispatcher system interfaces seamlessly with NEML2's existing models. They handle device-aware tensor operations, optimize memory transfers, and support asynchronous execution when applicable. This design ensures that batches of material points can be evaluated concurrently, substantially improving throughput compared to previous single-device or serial implementations. These improvements not only enhance the raw performance of NEML2 but also improve its usability in multiscale and high-fidelity simulations, where the simultaneous evaluation of large material point batches is critical. Benchmarks included in the report demonstrate the system’s scalability, highlighting its effectiveness when leveraging modern GPU architectures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Ionic Interdiffusion at Cathode|Solid-Electrolyte Interface: A Machine Learning–Assisted Multiscale Investigation and Mitigation Strategies

Future lithium batteries are expected to use solid electrolytes to achieve higher energy density and fast charge capabilities. However, most solid electrolytes are thermodynamically unstable against layered oxide cathodes. In this study, the stability of LiCoO2 (LCO) cathode with Li10GeP2S12 (LGPS) solid electrolyte is investigated using ab initio molecular dynamics (AIMD) and machine learning molecular dynamics (MLMD). The propensity of ionic interdiffusion, formation of a passivating interphase layer, and corresponding decay in cell performance is addressed using a continuum model. Large-scale MLMD simulations confirm that the LCO|LGPS interface permits interdiffusion of cobalt (Co) and other ionic species, leading to the formation and growth of a resistive interphase and to dramatic capacity fade even in the first cycle. We examine the literature evidence that incorporating a thin layer of LiNb0.5Ta0.5O3 (LNTO) between LCO and LGPS prevents the interdiffusion of ions. Atomistic simulations suggest that substituting lithium (Li) in LNTO with Co is thermodynamically unfavorable, thereby inhibiting ionic interdiffusion. The stable Nb5+/Ta5+ states form a rigid metal-oxide framework, which consequently also prevents the substitution of niobium (Nb) or tantalum (Ta). However, continuum-level analysis suggests that the higher mechanical stiffness of LNTO can lead to interfacial delamination between the LCO and LNTO. This phenomenon reduces the effectiveness of the protective layer. This paper, therefore, highlights the need to develop novel interlayers that balance low ionic interdiffusion with low mechanical stiffness.

Ncube, Musawenkosi K.↗
DOE OSTI2025

The role of moisture in MgCl 2 salt: A multiscale approach to TES performance

This study presents a comprehensive multiscale analysis to evaluate the influence of moisture on the thermal performance of thermal energy storage systems using magnesium chloride (MgCl 2 ) as the phase-changing material. The system uses graphite foam with 90% relative density to enhance thermal conductivity. The analysis includes thermal conductivity calculations and specific heat capacity for different hydrate phases of MgCl 2 : Anhydrous, Mono, Di, Tetra, and Hexa. These properties were evaluated for the first time using the phonon density of states from Density Functional Theory simulations. Results showed that thermal conductivity decreased, while specific heat capacity increased by a factor of two as the phase changed from anhydrous to hexahydrate. Meso-scale models were created to account for the anisotropy of graphite foam and property variations of the MgCl 2 hydrate phases. Asymptotic Expansion homogenization simulations determined the anisotropic thermal conductivity for all phases. This unique methodology improved simulation accuracy, which matched experimental data for anhydrous MgCl 2 . A parametric study examined various operating conditions and their effect on TES performance. It revealed that higher charging temperatures did not enhance exergy efficiency, but increased discharge mass flow rates improved it due to better heat transfer. Thermal performance evaluated by the exergy efficiency remained consistent across hydrate systems under the tested conditions. Furthermore, the study suggests that this uniformity is linked to missing key data, particularly latent heat and melting point for different hydrates. Overall, it highlights the importance of multiscale effects and accurate material properties in designing and optimizing TES systems.

Molten salt degradation↗
DOE OSTI2025

Multiscale maps of Active Layer Depth for Teller site Mile Marker 27 and Kougarok Mile Marker 80, Seward Peninsula, AK

Remote sensing maps of active layer depth derived from Unmanned Areal System (UAS) data. The UAS datasets were stepwise scaled until matching the AVIRIS-NG (Airborne Visible / Infrared Imaging Spectrometer - Next Generation) and Sentinel-2 spatial resolutions. Using the field observed Active Layer Depth (ALD) measurement in combination with spectral and topographic predictors derivatives from DJI UAS imagery, we used a spatially explicit RF regression model to predict and map ALD across our study landscapes. This package includes maps for Next-Generation Ecosystem Experiment Arctic (NGEE Arctic)’s Teller Mile Marker (MM) 27, and Kougarok MM80 (aka Mile 80) watersheds. The field, map data, and metadata are provided as geoTIF and text (*.csv) formats. These datasets are provided in support of Hantson et al., 2024 (accepted) “Scaling Arctic landscape and permafrost features improves active layer depth modeling”

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Dynamic density functional theory of polymers with salt in electric fields

Here we present a dynamic density functional theory for modeling the effects of applied electric fields on the local structure of polymers with added salt (polymer electrolytes). Time-dependent equations for the local electrostatic potential and volume fractions of polymer, cation, and anion of added salt are developed using the principles of linear irreversible thermodynamics. For such a development, a field theoretic description of the free energy of polymer melts doped with salts is used, which captures the effects of local variations in the dielectric function. Connections of the dynamic density functional theory with experiments are established by relating the three phenomenological Onsager’s transport coefficients of the theory to the mutual diffusion of electrolyte, ionic conductivity, and transference number of one of the ions. The theory is connected with a statistical mechanical model developed by Bearman and Kirkwood [J. Chem. Phys. 28, 136 (1958)] after relating the three transport coefficients to friction coefficients. The steady-state limit of the dynamic density functional theory is used to understand the effects of dielectric inhomogeneity on the phase separation in polymer electrolytes. The theory developed here provides not only a way to connect with experiments but also to develop multi-scale models for studying connections between local structure and ion transport in polymer electrolytes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Development of Composite Photocatalyst Materials that are Highly Selective for Solar Hydrogen Production and their Evaluation in Z-Scheme Reactor Designs

The key technology gap preventing a vertically stacked dual-bed particle suspension reactor from achieving the DOE MYRD&D ultimate cost target for H 2 production remains the lack of materials in particle form factor that exhibit ≥10% solar-to-H 2 energy conversion (STH) efficiency as a suspension. Therefore, our project goals centered around strategies to increase the STH efficiency by enhancing photophysical properties of perovskite oxide particles including increased visible-light absorption, increased selectivity for electrocatalysis of the H 2 evolution reaction (HER) and the O 2 evolution reaction (OER) through development of ultrathin oxide coatings, correlating composition and structure to function, and improving understanding of multiscale transport and kinetic processes.

08 HYDROGEN↗
DOE OSTI2026

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↗
DOE OSTI2025

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↗
DOE OSTI2024

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↗
DOE OSTI2026

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↗
DOE OSTI2025

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↗
DOE OSTI2024

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↗
DOE OSTI2024

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↗
DOE OSTI2025

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↗
DOE OSTI2025

Atomic-scale understanding of oxide growth and dissolution kinetics of Ni-Cr alloys

Aqueous corrosion of metals is governed by formation and dissolution of a passivating, multi-component surface oxide. Unfortunately, a detailed atomistic description is challenging due to the compositional complexity and the need to consider multiple kinetic factors simultaneously. To this end, we combine experiments with a first-principles-derived, multiscale computational framework that transcends thermodynamic descriptions to explicitly simulate the kinetic evolution of surface oxides of Ni-Cr alloys as a function of composition, temperature, pH, and applied voltage. In the absence of pitting, we identify three distinct voltage regimes, which are kinetically dominated by oxide growth, dissolution, and competitive dissolution and reprecipitation. Evolving compositional gradients and oxide thickness are revealed, including a transition between a metastable Ni-Cr mixed oxide and a thick, porous Ni-dominated oxide. Beyond elucidating the underlying physics, we highlight the need for competing kinetics in models to properly predict the transition from passivation to corrosion. Our results provide a key step towards co-design of alloy composition alongside environmental conditions for sustainable use across a variety of critical energy and infrastructure applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

pH-Dependent Vibrational Dynamics Drives Excited-State Quenching in the Phycobiliprotein Complex PC645

Phycocyanin 645 (PC645) is a closed-form lightharvesting complex found in the lumen of the photosynthetic membrane of cryptophyte algae. These peripheral antenna complexes contain bilin chromophores that absorb sunlight and transfer excitation energy to the core antenna complexes embedded in the thylakoid membrane. The location of cryptophyte antenna complex on the luminal side of the membrane is unusual. During photosynthetic activity, the pH of the lumen drops, by up to two pH units. There is little known about how this pH-change affects the light-harvesting complexes. In this study, we report multiscale simulations using a computationally efficient density functional tight-binding framework to investigate the spectroscopy and excitation energy transfer in the PC645 complex. Complementary experiments were conducted using both steady-state and time-resolved spectroscopic measurements at low, neutral, and high pH values. Our study shows that (de)protonation of specific bilin pigments, namely, the mesobiliverdins (MBVs), modulates the excitation energies, excitonic couplings, and spectral densities. These changes cause excitation transfer rates to increase by up to a factor of two to three, leading to pH-dependent energy transfer pathways in the complex. Using this model, we calculated the pH-dependent fluorescence quantum yield of the system, obtaining quantitative agreement with the experimental results. These computational simulations, supported by experiments, identify MBVs as a more prominent excitation sink than previously realized, and that this role is tuned by pH.

Maity, Sayan [Constructor Univ., Bremen (Germany);↗