Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale Continuum 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
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 OSTI2020

ChIMES: A Machine-Learned Interatomic Model Targeting Improved Description of Condensed Phase Chemistry in Energetic Materials

In this report we detail completion of a Physics and Engineering Model Level Two Milestone targeting improved reactive interatomic potentials (IAPs) for energetic materials (EM) through machine learning. The specific goals of this milestone were to develop, validate, and document a new reactive molecular dynamics method for EM, based on machine learning by (1) generating databases of first-principles-derived forces, stresses, and energies for HN3 and 3,4-bis(3-nitrofurazan- 4-yl)furoxan (DNTF) (2) generate atomistic force fields from these databases via ML, and (3) benchmark model performance against first principles calculations. These goals were achieved by (1) further developing a machine learned reactive IAP and generation approach (i.e. the Chebyshev Interaction Model for Efficient Simulation or “ChIMES”), for which resulting IAPs can approach the predictive power of quantum-mechanical approaches at a fraction of the computational expense, and (2) applying the ChIMES framework to develop models for HN3 and DNTF. We find that for simple energetic materials like HN3, high accuracy ChIMES models can be obtained through application of a fitting approach that does not use active machine learning. We demonstrate the suitability of ChIMES models for simulations involving EM by using the HN3 model in multiscale shock technique simulations to predict the HN3 Chapman-Jouguet detonation state and investigate chemical evolution out to 1 ns following shock compression. This model is then used in larger direct shock (DS) simulations for a preliminary investigation of how bubbles (i.e. voids) influence material response under shock compression. We find that more complex EM (i.e. DNTF) necessitate a more sophisticated fitting approach, and develop a new active learning method and python tool to meet this challenge. We demonstrate that this fitting approach yields ChIMES models that out-perform commonly used standard reactive IAPs as well as semi-empirical quantum methods, and discuss the systematic improvability of these actively learned ChIMES models. We also describe challenges related to model development for EM such as DNTF, for which few experimental or previous simulation data are available (e.g. which could otherwise inform generation of training data). To overcome this issue, we establish a semi-empirical quantum ChIMES capability which can be used to efficiently map out relevant thermodynamic and configurational space, and generate ChIMES-IAP training data in a multiscale manner. We also show that these semi-empirical quantum ChIMES models can be used to generate predictions for the shock Hugoniot (the Hugoniot is the locus of thermodynamic states found in a shocked material) equation of state, investigate related thermochemistry, and explore carbon condensation following shock compression. This work represents a substantial advance in our atomistic modeling capability for EM that will provide much needed information on the chemistry of detonation for continued development of continuum models based on the Cheetah thermochemical code.

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

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Finite-element-based simulations of electrodes for CO 2 cascade reduction reactions

The multielectron reduction of CO 2 to liquid fuels could be a path to scalable energy storage, but reaching this goal requires major advances in catalysis and systems engineering. Cascade catalysis, which couples sequential reactions without isolating intermediates, has emerged as a promising route to enhance selectivity and efficiency in CO 2 reduction (CO 2 R). In this review, we examine how finite-element-based simulations of continuum model [finite element method (FEM)] approaches are being used to analyze and guide CO 2 R cascade systems. We first outline the fundamentals of cascade catalysis and recent advances in catalytic materials (metallic, molecular, and hybrid architectures). We then focus on FEM developments at the electrode and device scales, emphasizing how these models capture transport phenomena, local microenvironments, and geometry-dependent effects. To clarify design principles, we present case studies of cascade electrodes organized in systems without and with integrated semiconductors. We further emphasize the integration of FEM with multiscale frameworks (density functional theory, molecular dynamics, kinetic Monte Carlo) and its role in bridging atomic-level insights with device-level performance. Finally, we identify current limitations and future prospects, including improved boundary conditions, coupling with operando experiments, and machine learning-accelerated model development. Together, these insights provide design principles for next-generation CO 2 R cascade systems for efficient solar fuel production.

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

Model for Humidity-Mediated Diffusion on Aluminum Surfaces and Its Role in Accelerating Atmospheric Aluminum Corrosion

Bare aluminum metal surfaces are highly reactive, which leads to the spontaneous formation of a protective oxide surface layer. Because many subsequent corrosive processes are mediated by water, the structure and dynamics of water at the oxide interface are anticipated to influence corrosion kinetics. Here, using molecular dynamics simulations with a reactive force field, we model the behavior of aqueous aluminum metal ions in water adsorbed onto aluminum oxide surfaces across a range of ion concentrations and water film thicknesses corresponding to increasing relative humidity. We find that the structure and diffusivity of both the water and the metal ions depend strongly on the humidity of the environment and the relative height within the adsorbed water film. Aqueous aluminum ion diffusion rates in water films corresponding to a typical indoor relative humidity of 30% are found to be more than 2 orders of magnitude slower than self-diffusion of water in the bulk limit. Connections between metal ion diffusivity and corrosion reaction kinetics are assessed parametrically with a reductionist model based on a 1D continuum reaction–diffusion equation. Our results highlight the importance of incorporating the properties specific to interfacial water in predictive models of aluminum corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Dynamics and lipid membrane coupling of the RAS-RAF complex revealed via multiscale simulations

To gain molecular and mechanistic insights into initiation of the RAS-RAF signaling cascade, we developed and used a combination of multiscale simulation and experimental approaches. The influence and impact of the membrane on RAS and RAF proteins is a factor we are just beginning to understand and appreciate in more detail. Molecular simulation is an ideal methodology to further study this complicated relationship between the membrane and associated proteins. Our previous work using Multiscale Machine-learned Modeling Infrastructure investigated different lipid compositions solely around the KRAS4b protein and the interplay between protein behavior and these membrane environments. Multiscale Machine-learned Modeling Infrastructure uses machine learning to couple adjacent simulation scales and has been efficiently scaled across some of the world’s largest high-performance computers. Recently, we have expanded this multiresolution framework to include the all-atom simulation scale and to incorporate the RAF RBDCRD domains. Here, we present the overall analysis results from this new simulation campaign comprising a mixture of RAS and RAF RBDCRD proteins. Approximately 35,000 coarse-grained and 10,000 all-atom molecular dynamics simulations were completed, sampled from a variety of protein/lipid composition configurations that were generated from a micron-scale continuum simulation containing hundreds of copies of the proteins. Our studies suggest that orientations of the RAS-RBDCRD complex on the membrane occupy distinct configurational states, and the spatial patterns of lipid arrangements around these different protein states are unique to each state. The extent and size of lipid “fingerprints” imposed on the membrane by the RAS-RBDCRD protein complex are significantly larger than observed for just the RAS protein on its own. These protein complexes strongly associate, but we do not observe statistically significant preferred protein-protein orientations. These observations indicate that spatial colocalization of RAS-RBDCRD proteins in the same vicinity may be assisted by specific membrane environments, acting to increase the probability of signaling complex formation.

Carpenter, Timothy S. [Lawrence Livermore National↗
DOE OSTI2022

Continuum Modeling of Porous Electrodes for Electrochemical Synthesis

Electrochemical synthesis possesses substantial promise to utilize renewable energy sources to power the conversion of abundant feedstocks to value-added commodity chemicals and fuels. Of the potential system architectures for these processes, only systems employing 3-D structured porous electrodes have the capacity to achieve the high rates of conversion necessary for industrial scale. However, the phenomena and environments in these systems are not well understood and are challenging to probe experimentally. Fortunately, continuum modeling is well-suited to rationalize the observed behavior in electrochemical synthesis, as well as to ultimately provide recommendations for guiding the design of next-generation devices and components. In this review, we begin by presenting an historical review of modeling of porous electrode systems, with the aim of showing how past knowledge of macroscale modeling can contribute to the rising challenge of electrochemical synthesis. We then present a detailed overview of the governing physics and assumptions required to simulate porous electrode systems for electrochemical synthesis. Here, leveraging the developed understanding of porous-electrode theory, we survey and discuss the present literature reports on simulating multiscale phenomena in porous electrodes in order to demonstrate their relevance to understanding and improving the performance of devices for electrochemical synthesis. Lastly, we provide our perspectives regarding future directions in the development of models that can most accurately describe and predict the performance of such devices and discuss the best potential applications of future models.

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

Micro-Mechanically Guided High-Throughput Alloy Design Exploration Towards Metastability-Induced H Embrittlement Resistance

We develop a high-throughput approach for studying H embrittlement (HE)-resistance in alloys, which is based on combinatorial compositional screening of metastability effects by in situ scanning electron microscopy H-analyses. The project objective included: (i) technique development of high-throughput screening (HTS) for HE-resistance, (ii) discovery of new metallic materials with superior HE-resistance, (iii) Multiscale verification of HE-resistance of the new alloys and H-barrier layers, from atomic scale to an engineering scale. The investigation focused on a model system that enables exploration of different metastable states in Fe-based complex-concentrated alloys. Composition spread islands with hundreds of varying compositions are fabricated on a single substrate using a combinatorial co-sputtering technique. To screen these alloys, we developed and employed a home-built SEM-based integrated analysis system capable of characterizing H-permeability, H-trapping, H-influence on mechanical properties, and other material properties, as well as atomistic to continuum simulations to study the underlying physics.

08 HYDROGEN↗
DOE OSTI2024

Multiscale modeling-enabled design of multifunctional composites

This study aims to create a comprehensive model that considers multiple scales and physics for predicting the electromechanical behavior of fiber-reinforced composites enhanced with barium titanate (BaTiO3). In our earlier work, we have demonstrated that depositing BaTiO3 microparticles of 200-nm-diameter, on fiber surfaces during fiber-reinforced composite fabrication enhances mechanical strength, passive self-sensing, and energy harvesting properties. The key is to carefully control the microparticle concentration to prevent agglomeration. Since the particles are micron-sized, understanding how agglomeration affects the composites' electromechanical properties is crucial for guiding such multifunctional materials’ design. This study introduces a micromechanics-based approach to explore the impact of microparticle dispersion on the bulk composites' electromechanical properties. Insights gained from this investigation are applied in experiments, enabling accurate predictions of mechanical and self-sensing responses in BaTiO3-enhanced fiber-reinforced composites. Micro-level findings from this computational approach can be integrated into larger continuum models to comprehensively capture the electromechanical behavior of the composite structures at bulk scale. The proposed model is validated by comparing predictions with experimental results, accounting for the nonlinear mechanical and electromechanical behaviors of constituent materials. Consequently, this computational model serves as a digital platform for efficiently designing multifunctional composites.

Gupta, Sumit↗
DOE OSTI2021

Modeling the Nonlinear Rheology of Polymer Additive Manufacturing

This report summarizes molecular and continuum simulation studies focused on developing physics - based predictive models for the evolution of polymer molecular order during the nonlinear processing flows of additive manufacturing. Our molecular simulations of polymer elongation flows identified novel mechanisms of fluid dissipation for various polymer architectures that might be harnessed to enhance material processability. In order to predict the complex thermal and flow history of polymer realistic additive manufacturing processes, we have developed and deployed a high - performance mesh - free hydrodynamics module in Sandia's LAMMPS software. This module called RHEO – short for Reproducing Hydrodynamics and Elastic Objects – hybridizes an updated - Lagrange reproducing - kernel method for complex fluids with a bonded particle method (BPM) to capture solidification and solid objects in multiphase flows. In combination, our two methods allow rapid, multiscale characterization of the hydrodynamics and molecular evolution of polymers in realistic processing geometries.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Novel Multiphysics Multiscale Multiporosity Shale Gas Transport Model for Geomechanics/Flow Coupling in Steady and Transient States

Summary A novel multiphysics multiscale multiporosity shale gas transport (M3ST) model was developed to investigate shale gas transport in both transient and steady states. The microscale model component contains a kerogen domain and an inorganic matrix domain, and each domain has its own geomechanical and gas transport properties. Permeabilities of various shale cores were measured in the laboratory using a pulse decay permeameter (PDP) with different pore pressure and confining stress combinations. The PDP-measured apparent permeability as a function of pore pressure under two effective stresses was fitted using the microscale M3ST model component based on nonlinear least squares fitting (NLSF), and the fitted model parameters were able to provide accurate model predictions for another effective stress. The parameters and petrophysical properties determined in the steady state were then used in the transient-state, continuum-scale M3ST model component, which performed history matching of the evolutions of the upstream and downstream gas pressures. In addition, a double-exponential empirical model was developed as a powerful alternative to the M3ST model to fit laboratory-measured apparent permeability under various effective stresses and pore pressures. The developed M3ST model and the research findings in this study provided critical insights into the role of the multiphysics mechanisms, including geomechanics, fluid dynamics and transport, and the Klinkenberg effect on shale gas transport across different spatial scales in both steady and transient states.

Engineering↗
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 OSTI2020

Coordinated Multiscale Modeling and Synthesis of Novel Nanostructured Composite Membranes for Solar Fuels Generation (Final Report)

Generation of carbon-based liquid fuels from reduction of CO 2 using solar energy requires the development of improved membranes that provide good ionic conductivity and mechanical properties while minimizing the crossover of gases and reduction products. Using closely coordinated multiscale modeling, synthesis and characterization studies we propose to investigate and develop novel nanostructured composite membranes with applications to conversion of carbon dioxide to storable chemical fuels. The envisioned proton exchange membranes (PEMs) will be based on a support of self- assembling polymer-modified nanoparticles with controllable porosity and tortuosity. The use of a nanoparticle support will provide control of transport and mechanical properties of the membranes. Tuning the copolymer modification of nanoparticles (polymer architectures, brush grafting density) will allow to optimize the ion (proton) transport via the Grotthuss mechanism through narrow hydrated channels while minimizing crossover of gases and reaction products due to the high density of the grafted brush and use of a non-ionic conducting glassy polymers to fill spaces between polymer-grafted nanoparticles. In Phase I synthesis of the copolymers and their attachment to the nanoparticle support will be guided by multiscale simulations that include atomistic molecular dynamics and continuum-level transport simulations. The former will include reactive and non-reactive atomistic simulations that will provide key insight into nanoscale polymer morphology and ion/molecular transport mechanisms need to optimize the copolymer structure. The latter will provide insight into the role of nanoparticle self-assembly on the global transport of ionic and molecular species. The proton conduction and CO 2 /methanol permeability of the membranes will be experimentally characterized as well as simulated and will be correlated with morphology and polymer structure. Prototypes of optimal membranes will be fabricated and characterized at the end of the project. Successful development of improved nanostructured composite PEMs and demonstration of the simulation-guided materials-by-design paradigm will result in our ability to design and synthesize membranes for a wide variety of solar fuels generation and related electrochemical applications.

14 SOLAR ENERGY↗
DOE OSTI2024

Development of multi-scale computational frameworks to solve fusion materials science challenges

Over the past two decades, the US-DOE has funded multiple projects that rely on high-performance computing and exascale computing platforms to accelerate scientific discoveries and address grand scientific challenges, such as harnessing fusion energy. In this article, we review in detail one of these efforts aimed at enhancing our capability to model plasma-facing materials subject to plasma and high-energy ion/neutron irradiation. The plasma surface interactions project has built a multi-scale modeling framework where many of the plasma- and high-energy ion/neutron irradiation-induced effects occurring in tungsten are explored. Here, this knowledge is used to develop atomistically-informed, high-fidelity continuum and meso-scale models that can be validated against experiments. We review the developments within this project, with attention to experimental validation efforts, and specifically highlight activities associated with: helium bubble bursting and equation of state, and hydrogen-helium interactions in tungsten; atomistically-informed model development for beryllium-tungsten material mixing; coupling of scrape-of-layer plasma, sheath and material models; and coupling of stochastic cluster-dynamics and crystal plasticity models to address radiation effects in tungsten under stress. Finally, we present how the project is preparing for future computational architectures, for instance through efforts to adapt atomistic methods to exascale computing.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multiscale Reactive Model for 1,3,5-Triamino-2,4,6-trinitrobenzene Inferred by Reactive MD Simulations and Unsupervised Learning

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

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

Long-range hydrogen-binding effects of carbide interfaces in iron

A micromechanics model was developed to evaluate the elastic binding energy between carbide precipitates and hydrogen interstitials using Eshelby's equivalent inclusion method. In this work, density functional theory (DFT) simulations were performed to obtain the material-specific quantities, e.g., lattice constants and the elastic constants, for the continuum model. Using this model, we find that for coherent carbide precipitates, hydrogen atoms are more likely to bind on the broad surfaces of the disk-like precipitates, which is consistent with experimental observations. For semicoherent and incoherent precipitates, our model suggests that it is possible for semicoherent precipitates to have significant hydrogen binding capability while there is no hydrogen-binding capability of incoherent precipitates, which also agrees with experimental findings. In addition, several factors that influence the binding energies between hydrogen atoms and carbide precipitates were quantitatively analyzed, including the precipitate size, morphology, orientation, and interface. These collective results include both the position and the value of the strongest hydrogen-binding interaction for a wide range of carbide stoichiometries, which contributes to our understanding of hydrogen trapping in steel-based materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Accelerating multiscale electronic stopping power predictions with time-dependent density functional theory and machine learning

Knowing the rate at which particle radiation releases energy in a material, the “stopping power,” is key to designing nuclear reactors, medical treatments, semiconductor and quantum materials, and many other technologies. While the nuclear contribution to stopping power, i.e., elastic scattering between atoms, is well understood in the literature, the route for gathering data on the electronic contribution has for decades remained costly and reliant on many simplifying assumptions, including that materials are isotropic. We establish a method that combines time-dependent density functional theory (TDDFT) and machine learning to reduce the time to assess new materials to hours on a supercomputer and provide valuable data on how atomic details influence electronic stopping. Our approach uses TDDFT to compute the electronic stopping from first principles in several directions and then machine learning to interpolate to other directions at a cost of 10 million times fewer core-hours. We demonstrate the combined approach in a study of proton irradiation in aluminum and employ it to predict how the depth of maximum energy deposition, the “Bragg Peak,” varies depending on the incident angle—a quantity otherwise inaccessible to modelers and far outside the scales of quantum mechanical simulations. The lack of any experimental information requirement makes our method applicable to most materials, and its speed makes it a prime candidate for enabling quantum-to-continuum models of radiation damage. The prospect of reusing valuable TDDFT data for training the model makes our approach appealing for applications in the age of materials data science.

36 MATERIALS SCIENCE↗