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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2023

NASA GRC ICME Schema for Materials Data Management: An Executive Summary

Integrated Computational Materials Engineering (ICME) has received a growing emphasis in attention due its potential impact on rapid material design, reduction in cost and time to market for new applications, and the promise of ‘fit-for-purpose’ materials coupled with recent advances in high performance computing and material characterization tools. However, for an organization to implement ICME practices for material discovery and design, a series of both technical and cultural challenges must be overcome to foster an environment that enables efficient, traceable, and predictive multiscale simulations of material behavior to enable virtual design of materials. In 2016, NASA sponsored a 2040 Vision study to define the potential 25-year future state required for integrated multiscale modeling of materials and systems to improve both the associated time and cost for aerospace and aeronautical innovation. The study envisions a cyber-physical-social ecosystem of experimentally validated computational models, tools, and techniques, along with the associated digital tapestry, that can enable rapid, optimized, ‘fit-for-purpose’ design of materials, components, and systems. A key requirement for such an ecosystem is the development of a robust information management system for materials across their full lifecycle, including material pedigree, experimental (real) and virtual (simulation) data, developed material models, and the implementation of models in engineering applications, such that process-structure-property-performance relationships can be established, thereby enabling the virtual design and optimization of materials. Such an information management system must be able to effectively capture: i) material information at each length scale; ii) test data and analysis; iii) associated material models; and iv) material and model deployment in engineering applications. These systems must also provide traceability between experimental and virtual representations of the material to ensure, when appropriate, the material digital twin is maintained. Additionally, this robust material information management system must be able to seamlessly connect with both commercial and an organization’s in-house software tools, be they analysis tools, other material databases, product lifecycle management (PLM) or simulation data management (SDM) tools, etc., such that automation of the design and analysis of a material across multiple length scales is possible. In this paper, an executive summary of the NASA GRC ICME Schema for materials information management is presented. The database best practices and schema design philosophy specifically for ICME materials data management and an overview description of each element in the schema is given, along with its associated role in an ICME workflow. Additionally, auxiliary tools that interact with the database and provide judicious automation with regards to importing, exporting, and analyzing materials data are presented. Such tools are critical to an ICME ecosystem, not only for their role in enabling optimization, but also in relieving users of tedious manual tasks, thus helping to promote adoption and combat the cultural challenges organizations face in enabling ICME.

Materials↗
DOE OSTI2024

Limitations of Hydrogen Detection After 150 Years of Research on Hydrogen Embrittlement

Hydrogen's significance in contemporary society lies in its remarkable energy density, yet its integration into the worldwide energy grid presents a substantial challenge. Exposing materials to hydrogen environments leads to degradation of mechanical properties, damage, and failure. While the current approach for assessing hydrogen's impact on materials involves mainly multiscale modeling and mechanical testing, there exists a significant deficiency in detecting the intricate interactions between hydrogen and materials at the nanoatomic scales and under in situ conditions. This perspective review highlights the experimental endeavors aimed at bridging this gap, pointing toward the imminent need for new experimental techniques that can detect and map hydrogen in materials’ microstructures and their site‐specific dependencies.

Tunes, Matheus A.↗
DOE OSTI2025

Direct Imaging of Hydrogen‐Driven Dislocation and Strain Field Evolution in a Stainless Steel Grain

Hydrogen embrittlement (HE) poses a significant challenge to the durability of materials used in hydrogen production and utilization. Disentangling the competing nanoscale mechanisms driving HE often relies on simulations and electron-transparent sample techniques, limiting experimental insights into hydrogen-induced dislocation behavior in bulk materials. This study employs in situ Bragg coherent X-ray diffraction imaging to track three-dimensional (3D) dislocation and strain field evolution during hydrogen charging in a bulk grain of austenitic 316 stainless steel. Tracking a single dislocation reveals hydrogen-enhanced mobility and relaxation, consistent with dislocation dynamics simulations. Subsequent observations reveal dislocation unpinning and climb processes, likely driven by osmotic forces. Additionally, nanoscale strain analysis around the dislocation core directly measures hydrogen-induced elastic shielding. These findings experimentally validate theoretical predictions and offer mechanistic insights into hydrogen-driven dislocation behavior. The quantified nanoscale phenomena serve as critical inputs for multiscale modeling frameworks to predict bulk material responses and accelerate the development of HE-resistant alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Evaluating grid stress and reliability in future electricity grids across a range of demand, generation mix, and weather trends

The reliability of power grids in the future will depend on how system planners account for the integration of new technologies, extreme weather events, and uncertainties in demand growth from increased electrification and data centers. This study introduces an open-source, multisectoral, multiscale modeling framework that projects grid stress and reliability trends between 2020 and 2055 in the Western Interconnection of the United States. The framework integrates global to national energy-water-land dynamics with power plant siting and hourly grid operations modeling. We analyze future wholesale electricity price shocks and unserved energy events across eight scenarios spanning a range of population growth and economic change, generation mixes, and weather conditions. Our results show future grids with high percentage of non-renewable generation and strong economic growth are characterized by higher reliability and lower wholesale electricity prices than lower growth scenarios because of larger reliance on dispatchable generators and lower fossil fuel extraction costs. Scenarios with high percentage of renewable resources have lower median but more volatile wholesale electricity prices as well as more frequent and severe unserved energy events compared to scenarios relying more on dispatchable generators. These events occur because higher proportion of solar and wind energy causes net demand curves to deepen during midday (duck curves get progressively severe), exacerbating the challenge of meeting demand during summer evening peaks. This study suggests that robust and co-optimized transmission and energy storage planning could help maintain low wholesale electricity prices and high reliability levels in future electricity grids across uncertainties in generation mixes.

Electric grid reliability↗
DOE OSTI2025

Mesh-based super-resolution of fluid flows with multiscale graph neural networks

A graph neural network (GNN) approach is introduced in this work which enables mesh-based three-dimensional super-resolution of fluid flows. In this framework, the GNN is designed to operate not on the full mesh-based field at once, but on localized meshes of elements (or cells) directly. To facilitate mesh-based GNN representations in a manner similar to spectral (or finite) element discretizations, a baseline GNN layer (termed a message passing layer, which updates local node properties) is modified to account for synchronization of coincident graph nodes, rendering compatibility with commonly used element-based mesh connectivities. Furthermore, the architecture is multiscale in nature, and is comprised of a combination of coarse-scale and fine-scale message passing layer sequences (termed processors) separated by a graph unpooling layer. The coarse-scale processor embeds a query element (alongside a set number of neighboring coarse elements) into a single latent graph representation using coarse-scale synchronized message passing over the element neighborhood, and the fine-scale processor leverages additional message passing operations on this latent graph to correct for interpolation errors. Demonstration studies are performed using hexahedral mesh-based data from Taylor–Green Vortex and backward-facing step flow simulations at Reynolds numbers of 1600 and 3200. Through analysis of both global and local errors, the results ultimately show how the GNN is able to produce accurate super-resolved fields compared to targets in both coarse-scale and multiscale model configurations. Reconstruction errors for fixed architectures were found to increase in proportion to the Reynolds number. Geometry extrapolation studies on a separate cavity flow configuration show promising cross-mesh capabilities of the super-resolution strategy.

Backward-facing step↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Grain2mesh: A Python and cubit mesh generator from unprocessed mesoscale images

Predicting bulk behavior from microscale features constitutes a key objective in multiscale modeling research, often involving numerical models composed of finite elements that capture the diversity of constituent phases, shapes, and orientations within the material. The Grain2mesh toolbox allows the user to input unprocessed mesoscopic images for automatic segmentation, pre-processing, quality control, and numerical mesh generation. The numerical mesh generation incorporates Cubit routines to generate robust multi-phase mesh structure for use in computational mechanics solvers. The python classes developed contain detailed documentation and examples to support standard usage and case-specific alternative options.

58 GEOSCIENCES↗
DOE OSTI2026

Pathways to High-Performance Salt Hydrate Thermochemical Energy Storage Materials and Systems

Thermochemical materials (TCMs) based on salt hydrates are promising for thermal energy storage as they combine high energy densities with low reaction temperatures. However, their adoption is hindered by poor structural integrity and degradation under hygrothermal cycling. Storage performance is governed not only by the chemical reaction, but also by the coupled thermo-chemo-mechanical behavior that evolves with cycling. Understanding and controlling this coupling across length scales (material-to-reactor) is necessary to improve TCM stability and lifetime. In this perspective, we discuss the shortcomings of current characterization approaches and emphasize the need for measuring transport properties and structural transformations using in situ techniques that capture the dynamic evolution of these materials. We also outline opportunities for multiscale modeling frameworks that link thermodynamics and mechanics, enabling predictive evaluation of composite architectures designed for cycling stability. We conclude by identifying research questions that must be addressed to transform TCMs into viable energy storage technologies.

Composites↗
DOE OSTI2026

How Does Water Dissociation Work in Bipolar Membranes?

Bipolar membranes (BPMs) create counteracting spatial gradients of pH and electrostatic potential in electrochemical systems, enabling applications in pH regulation, electrocatalysis, and separations. At the polarized junction of a BPM the water dissociation (WD, 2H2O ⇌ H3O+ + OH-) reaction can be driven, but it remains poorly understood. In this Perspective, we integrate molecular insights from bulk-water autoionization and the associated field effects with continuum descriptions of BPM electrostatics and experimental WD kinetic analyses to describe possible mechanisms of voltage-driven WD. Pristine BPM junctions highlight both the limits of primarily electric-field-driven WD and the practical challenges of junction stability at extreme reverse bias. Introducing heterogeneous catalyst layers, commonly metal oxides and graphene oxides, accelerates WD by orders of magnitude through hypothesized coupled effects in which surface acid-base functionality and high-density hydroxyl sites mediate proton-transfer steps, and catalyst mobile electronic/ionic charges redistribute the junction electric potential drop to shape the local electric fields and reactive microenvironments. Kinetic analyses suggest two regimes of heterogeneous WD mechanism, including field-driven ordering of interfacial water and a Second-Wien-Effect dissociation-barrier lowering. We conclude by defining the key unknown variables (local pH, electrostatic potential, catalyst charge state and relationships among mechanisms) and outlining experimental and multiscale modeling strategies needed for predictive WD catalysis and for controlling related ion-transfer reactions.

Wu, Yifan↗
DOE OSTI2024

Vertically Resolved Analysis of the Madden-Julian Oscillation Highlights the Role of Convective Transport of Moist Static Energy

We simulate the Madden-Julian oscillation (MJO) over an aquaplanet with uniform surface temperature using the multiscale modeling framework (MMF) configuration of the Energy Exascale Earth System Model (E3SM-MMF). The model produces MJO-like features that have a similar spatial structure and propagation behavior to the observed MJO. To explore the processes involved in the propagation and maintenance of these MJO-like features, we perform a vertically resolved moist static energy (MSE) analysis for the MJO (Yao et al., 2022, https://doi.org/10.1175/jas-d-20-0254.1). Unlike the column-integrated MSE analysis, our method emphasizes the local production of MSE variance and quantifies how individual physical processes amplify and propagate the MJO's characteristic vertical structure. We find that radiation, convection, and boundary layer (BL) processes all contribute to maintaining the MJO, balanced by the large-scale MSE transport. Furthermore, large-scale dynamics, convection, and BL processes all contribute to the propagation of the MJO, while radiation slows the propagation. Additionally, we perform mechanism-denial experiments to examine the role of radiation and associated feedbacks in simulating the MJO. We find that the MJO can still self-emerge and maintain its characteristic structures without radiative feedbacks. This study highlights the role of convective MSE transport in the MJO dynamics, which was overlooked in the column-integrated MSE analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Energy Infrastructure Futures: A Multiscale Evaluation of Projected Power Plant Siting Across the Western Interconnection

The US Western Interconnection is facing unprecedented challenges in the form of less predictable peak demand, increasingly diverse generating resources, and fast-growing loads due to the onset of artificial intelligence, hyperscale computing, and electrification. Projecting where future generation may be developed is critical to maintaining a robust and resilient electric grid under this mounting uncertainty and variability. Using an integrated multisectoral, multiscale modeling framework that links a human-Earth systems model, an hourly load model, a geospatial power plant siting model, and an hourly grid operations model, we evaluate the power plant landscape evolution under eight alternative futures between 2020 and 2055. These futures represent a wide but plausible range of atmospheric conditions, emissions constraints, and economic, technological, and population growth assumptions. We find that local-level development can vary substantially both by generation type and capacity buildout across these futures. Specific regions of the Western Interconnection are projected to see large amounts of capacity development regardless of the future scenario. We additionally determine that projected power plant locations are more heavily influenced by the cost to interconnect to the electric grid than the locational energy value.

Mongird, Kendall↗
DOE OSTI2025

Resolving Low Cloud Feedbacks Globally With E3SM High‐Res MMF: Agreement With LES but Stronger Shortwave Effects

This study investigates low cloud feedback in a warmer climate using global simulations from the High-Resolution Multi-scale Modeling Framework (HR-MMF), which explicitly simulates small-scale eddies globally. Two 5-year simulations—one with present-day sea surface temperatures (SSTs) and a second with SSTs warmed uniformly by 4 K—reveal a positive global shortwave cloud radiative effect (SWCRE = 0.3 W/m 2 /K), comparable to estimates from CMIP models. As the climate warms, significant reductions in low cloud cover occur over stratocumulus regions. This study is the first attempt to compare HR-MMF results with predictions from idealized large-eddy simulations from the CGILS intercomparison. Despite different underlying assumptions, we find qualitative agreement in SWCRE and inversion height changes between HR-MMF and CGILS predictions. This suggests reasonable credibility for the CGILS framework in predicting cloud responses under the out-of-sample conditions found in HR-MMF. However, the HR-MMF exhibits stronger SWCRE changes than predicted by CGILS. We explore potential causes for this discrepancy, examining variations in cloud-controlling factors (CCFs) and cloud conditions. Our results show a fairly homogeneous SWCRE response, with little systematic variation tied to the variations in CCFs. This reveals a dominant role for SST forcing in modulating SWCRE.

boundary-layer clouds↗
DOE OSTI2025

Advancing simulations of coupled electron and phonon nonequilibrium dynamics using adaptive and multirate time integration

Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.

Yao, Jia [California Institute of Technology (CalT↗
DOE OSTI2025

Gas-mediated defect engineering in earth-abundant Mn-rich layered oxides for non-aqueous sodium-based batteries

Gases are often by-products of battery materials during cell formation and degradation, affecting the cycle life and safety of rechargeable batteries. However, understanding gas-mediated (electro)-chemical reactions and nanoscale structural transformations during the synthesis of battery electrode materials remains challenging because of the lack of suitable characterization routes and the complexity of the interplay between thermodynamics and kinetics. Here, in this study, we use operando synchrotron X-ray diffraction, in situ transmission X-ray microscopy and multiscale modelling to elucidate the reaction pathways and microstructural defect development of earth-abundant Mn-rich layered oxides as positive electrode materials for sodium-based batteries. In particular, we demonstrate the dominant role of CO 2 over O 2 and H 2 O (g) in modulating the competition between entropy and enthalpy during solid-state synthesis. Using Ni 0.25 Mn 0.75 CO 3 as a model precursor, we reveal that CO 2 generation favours the formation of entropy-driven metastable intermediates, suppresses closed pore/nanovoids formation and decreases chemical heterogeneity and residual lattice strain of Mn-rich layered oxides during the synthesis. This result motivates a fast-sintering strategy to promote CO 2 release, which ultimately leads to improved chemo-mechanical and electrochemical stability of the Mn-rich positive electrodes when tested in non-aqueous Na metal coin cells.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Toward durable stacks: glass-ceramic sealants for intermediate-temperature protonic ceramic electrochemical systems

Protonic ceramic electrochemical cells (PCECs) are emerging as promising technologies for efficient energy conversion and hydrogen production because they operate at intermediate temperatures with improved efficiency and durability compared with conventional solid oxide electrochemical cells. However, the long-term reliability and commercialization of PCEC stacks remain strongly limited by the performance of sealants, which are required to maintain gas tightness, electrical insulation, and mechanical integrity under harsh thermal and chemical environments. Among various sealing approaches, glass-ceramic sealants are considered the most practical and scalable due to their excellent wettability, chemical tunability, and strong interfacial adhesion. This review provides a comprehensive overview of recent advances in glass-ceramic sealants for intermediate-temperature protonic ceramic electrochemical systems. The fundamental design principles of sealant compositions are first discussed, followed by recent developments in deposition methods, sintering strategies, surface treatments, and degradation monitoring techniques. Particular attention is given to the unique challenges associated with PCEC operating conditions, including hydrothermal degradation, interfacial reactions with barium-containing electrolytes, and thermal mismatch. Finally, future opportunities involving sustainable materials, multiscale modeling, additive manufacturing, and artificial intelligence-assisted sealant optimization are highlighted.

glass–ceramic sealants↗
DOE OSTI2025

Temperature dependence of vacancy/self-interstitial recombination volumes in copper

We use molecular dynamics to calculate rate coefficients for recombination of vacancies with self-interstitial atoms (SIAs) in Cu at temperatures from 300 to 700 K. From these results, we calculate vacancy/SIA recombination volumes and find that they decrease from around 290 (where is one atomic volume) at 300 K to 160 at 500 K and above. By counting the number of distinct pathways by which a stable SIA may migrate to a site of spontaneous recombination with a nearby vacancy, we find a lower bound estimate of 168 for the recombination volume. We furthermore rationalize its temperature dependence based on differences between the activation energies for recombination and SIA migration. Furthermore, our work sheds light on the fundamental nature of the recombination volume and provides information that may be incorporated into multiscale models of radiation response in solids.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Method of Finite Averages

The Method of Finite Averages (MoFA) is a rigorous multiscale modeling methodology for efficiently modeling multi-physical phenomena in heterogeneous porous media. The code developed in this project aims to perform the numerical calculations required to formulate, implement, and verify MoFA models for Earth and Energy systems (i.e., model verification refers to performing fully-resolved simulations of the systems and comparing their results to those of the models). In general, MoFA transforms partial differential equations (PDEs) describing the fine-scale physics of a system into coupled ordinary differential equations (ODEs)---in time---that describe the coarse-scale---or "average"---physical behaviors of the system. This transformation significantly expedites system simulation, as the coarse-scale ODEs involve vastly fewer degrees of freedom than the fine-scale PDEs. The code developed under this project will allow users to 1.) generate system geometries and numerical meshes, 2.) solve the PDE and ODE systems required for MoFA model formulation and implementation, 3.) solve the PDE systems required to obtain fully-resolved simulation results for model verification, and 4.) compare and plot results (e.g., the model and fully-resolved simulation solutions, the error between the solutions, etc.).

Pietrzyk, KyleM [Lawrence Livermore National Labor↗