Accelerating Large Multiscale Composite Simulations with a GNN/LSTM Microscale Surrogate
Explore the source record for details and available documents.
Multiscale modeling↗
SEARCH · Engineering Papers
Results for “multiscale modeling and simulation”
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.
Predicting Unreinforced Fabric Mechanical Behavior with Recurrent Neural Networks
Unreinforced woven fabrics are widely employed in various high-performance applications, including parachute deployment systems, airbags, and ballistic armor. The analysis of such materials is inherently complex due to the multiscale structure of these materials, and the dependence of macroscale behavior on changes that occur at lower scales. Previously, NASA’s Multiscale Analysis Tool (NASMAT) showed its capability in predicting unreinforced fabric behavior at the macroscale by capturing finite rotations that occur at the mesoscale. Though effective, the tool can face high computational cost for large, complex problems, motivating the need for the development of a surrogate model that can capture the same behavior. A recurrent neural network (RNN) was developed and trained on virtual NASMAT data to mimic the physics-based solutions while improving the computational runtime. The architecture of the RNN to best simulate the fabric behavior was carefully crafted based on heuristic knowledge of predicting physics-based temporal data, manual hyperparameter case studies, and Hyperband optimization.. The resultant model was able to predict a variety of stress-strain curves for fabrics with different mesoscale geometries, and was further validated by comparing to experimental data for the K706 style Kevlar plain-weave fabric, demonstrating the ability of the model to effectively capture the geometric changes in the fabric without explicitly calculating them, as is done in NASMAT. Furthermore, the tool showed its ability to improve on the runtime by a factor of 10 for fabric solutions compared to the multiscale tool, which would further enable the simulation of complex loading scenarios on unreinforced fabrics.
Fabric↗
Machine learning-enabled multiscale modeling of mechanical deformation of aluminum and Al-SiC nanocomposites
A machine learning-enabled multiscale framework is developed for modeling the mechanical response of both pure metal and nanoparticle-reinforced metal matrix nanocomposites (MMNCs). Using aluminum–silicon carbide (Al-SiC) as an example MMNC, atomistic simulations reveal three distinct deformation mechanisms (i.e., defect-free, dislocation-based, and interface separation) governed by the interfaces between the Al matrix and SiC nanoparticles. As compared with single crystal Al, the lattice undergoes a more abrupt failure once the dislocation network becomes extensive and void nucleation initiates, whereas in Al-SiC, nanoparticle interfaces enable a more gradual progression of damage. These mechanisms are captured through a combined classification-regression neural network surrogate model that bridges atomic-scale insights with continuum-scale finite element analysis. Machine learning-enabled multiscale modeling of pure Al accurately predicted strain localization and confirmed by in-situ scanning electron microscopic tensile testing on perforated Al specimens. This study underscores the promise of integrating physics-informed machine learning with hierarchical modeling to capture the interface dominated phenomena and guide the design of advanced MMNCs.
Al-SiC↗
Enabling AI-Accelerated Multiscale Modeling of Thrombogenesis at Millisecond and Molecular Resolutions on Supercomputers
We report the first congruent integration of HPC, AI, and multiscale modeling (MSM) for solving a mainstream biomechanical problem of thrombogenesis involving 6 million particles at record molecular-scale resolutions in space and at simulation rates of milliseconds per day. The two supercomputers, the IBM Summit-like AiMOS and our University’s SeaWulf, are used for scalability analysis of, and production runs with, the LAMMPS with our customization and AI augmentation and they attained optimal simulation speeds of 3,077 µs/day and 266 µs/day respectively. The long-time and large scales simulations enable the first study of the integrated platelet flowing, flipping, aggregating dynamics in one dynamically-coupled production run. The platelets’ angular and translational speeds, membrane particles’ speeds, and the membrane stress distributions are presented for the analysis of platelets’ aggregations.
Zhu, Yicong↗
A Multiscale Approach to Simulate Non‐Isothermal Multiphase Flow in Deformable Porous Materials
Coupled thermal, hydraulic, and mechanical processes in porous materials play important roles in several energy and environmental technologies. The Darcy-Brinkman-Biot (DBB) framework has proven effective in modeling multiphase fluid flow in deformable porous solids across both pore and Darcy scales, including in systems where fractures coexist with a porous matrix. In this study, we extend the DBB framework, originally designed for isothermal conditions, to address non-isothermal problems by incorporating an energy conservation equation. The resulting solver, hybridBiotThermalInterFoam, enables simulations of coupled multiphase fluid flow, heat transfer, and solid deformation in hybrid-scale systems containing both solid-free regions and ductile porous domains. The new solver is validated through comparisons with analytical solutions and, also, against established heat transfer solvers chtMultiRegionFoam and compressibleInterFoam. Further, a series of 2D and 3D case studies, including two-phase heat transfer in solid-free, static, or deformable porous media, highlights the solver's capacity to simulate complex flow dynamics and heat transport in systems involving high mobility ratios, viscous fingering, and fracture propagation. Our results establish the feasibility of incorporating thermal effects in simulations of a wide variety of energy geotechnics and environmental applications, including enhanced hydrocarbon recovery, soil remediation, and enhanced geothermal energy systems.
04 OIL SHALES AND TAR SANDS↗
Neural operator transformers capture bifurcating drift-wave turbulence in fusion plasma simulations
Self-consistent modeling of turbulence-driven transport is critical for optimizing confinement in magnetically confined fusion plasmas, such as tokamaks and stellarators. In particular, capturing the long-term co-evolution of turbulence, flow, and background plasma profiles remains computationally challenging. Direct numerical simulation of these multiscale, highly nonlinear processes is often demanding and impractical for real-time control or design optimization. To address this bottleneck, we investigate transformer-based neural operator partial differential equation surrogates for emulating the dynamics of drift-wave turbulence bifurcation mediated by zonal flows, using the modified Hasegawa–Wakatani (MHW) model as a prototypical system. We find that the finetuned neural operator model has excellent performance in capturing the multi-spatiotemporal-scales of MHW turbulence bifurcation and is robust to testing on rare and out-of-distribution dynamics. Specifically, we demonstrate that a single unified model accurately predicts both quasi-steady-state turbulence and a wide range of dynamical transition processes, such as nonlinear saturation, spontaneous suppression of turbulence, and the emergence of macroscopic zonal flows, over time horizons vastly exceeding the local turbulence correlation time. This computationally efficient approach establishes a strong foundation for fast, AI-based modeling of complex, multiscale phenomena in magnetized fusion plasmas.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Investigating the Effect of Water on the Mechanical Properties of Cellulose from Multiscale Molecular Dynamics Simulations
Classical molecular dynamics (MD) simulations provide insight into the structure and physicochemical properties of materials with atomic resolution. However, the length and time scales accessible to atomistic MD are orders of magnitude smaller than many relevant processes such as the response of a bulk material to experimentally accessible strain rates, which presents challenges when comparing models to experimental measurements. Bottom-up coarse-graining provides a means for systematically mapping atomistic information to lower resolution models to increase the length and time scales achievable by simulation. Cellulose is an abundant carbohydrate biopolymer with applications to many fields of research, such as materials science and renewable energy, due to its desirable mechanical properties and viability for conversion into biofuel. The effect of moisture content on the Young's modulus of cellulose is of special interest due to its native environment often being in the hydrated secondary plant cell wall and the grinding energy requirements for biomass feedstock preprocessing. The current work investigates the effects of water solvent on the Young's modulus of cellulose calculated from coarse-grained MD mechanical stress simulations. The coarse-grained model was parametrized from atomistic MD calculations of cellulose-cellulose potentials of mean force using umbrella sampling techniques under vacuum and solvated conditions. The Young's moduli of the coarse-grained cellulose assemblies parametrized from cellulose in vacuum or solvated in water were computed via mechanical stress simulations to highlight the importance of capturing solvent interactions for modeling the mechanical behavior of cellulose.
BASIC BIOLOGICAL SCIENCES,RADIATION PROTECTION AND↗
Wavelet flow for extragalactic foreground simulations
Extragalactic foregrounds in cosmic microwave background (CMB) observations are both a source of cosmological and astrophysical information and a nuisance to the CMB. Effective field-level modeling that captures their non-Gaussian statistical distributions is increasingly important for optimal information extraction, particularly given the low-noise observations from current and upcoming experiments. Here, we explore the use of Wavelet Flow (WF) models to tackle the novel task of modeling the field-level probability distributions of multi-component CMB secondaries and foregrounds. Specifically, we jointly train correlated CMB lensing convergence (κ) and cosmic infrared background (CIB) maps with a WF model and obtain a network that statistically recovers the input to high accuracy — the trained network generates samples of κ and CIB fields whose average power spectra are within a few percent of the inputs across all scales, and whose Minkowski functionals are similarly accurate compared to the inputs. Leveraging the multiscale architecture of these models, we fine-tune both the model parameters and the priors at each scale independently, optimizing performance across different resolutions. These results demonstrate that WF models can accurately simulate correlated components of CMB secondaries, supporting improved analysis of cosmological data. Our code and trained models can be found on this GitHub repo.
cosmological simulations↗
A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping
A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.
Feng, Min↗
Nonlinear ionization dynamics of hot dense plasma observed in a laser-plasma amplifier
Understanding the behaviour of matter under conditions of extreme temperature, pressure, density and electromagnetic fields has profound effects on our understanding of cosmologic objects and the formation of the universe. Lacking direct access to such objects, our interpretation of observed data mainly relies on theoretical models. However, such models, which need to encompass nuclear physics, atomic physics and plasma physics over a huge dynamic range in the dimensions of energy and time, can only provide reliable information if we can benchmark them to experiments under well-defined laboratory conditions. Due to the plethora of effects occurring in this kind of highly excited matter, characterizing isolated dynamics or obtaining direct insight remains challenging. High-density plasmas are turbulent and opaque for radiation below the plasma frequency and allow only near-surface insight into ionization processes with visible wavelengths. Here, the output of a high-harmonic seeded laser-plasma amplifier using eight-fold ionized krypton as the gain medium operating at a 32.8 nm wavelength is ptychographically imaged. A complex-valued wavefront is observed in the extreme ultraviolet (XUV) beam with high resolution. Ab initio spatio-temporal Maxwell–Bloch simulations show excellent agreement with the experimental observations, revealing overionization of krypton in the plasma channel due to nonlinear laser-plasma interactions, successfully validating this four-dimensional multiscale model. This constitutes the first experimental observation of the laser ion abundance reshaping a laser-plasma amplifier. The presented approach shows the possibility of directly modelling light-plasma interactions in extreme conditions, such as those present during the early times of the universe, with direct experimental verification.
36 MATERIALS SCIENCE↗
Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model: Modeling Archive
This archive contains the data and Python scripts required to reproduce the analyses and figures in the study: Gomez-Velez, J. D., Rathore, S. S., Cohen, M. J., & Painter, S. L. (2025). Hyporheic-zone Processes and Stream Oxygen Dynamics: Insights from a Multiscale Reactive Transport Model. Submitted to Water Resources Research. The analysis utilizes the subgrid model Advection Dispersion Equation with Lagrangian Subgrids (ADELS) implemented in the Advanced Terrestrial Simulator (ATS; https://amanzi.github.io/ats/stable/). In this case, the ATS and Amanzi versions are (1) ATS version 1.5.1_f5ba18f8 and (2) Amanzi version 1.6-dev_53444cca4. The repository includes a Jupyter Notebook and the necessary data (Pandas DataFrames stored as pickle files) to generate the figures for the manuscript. Additionally, it contains Python scripts to create ATS input files, run the ATS simulations, and post-process the results. Finally, it provides routines for parameter estimation using the Single-Station Metabolism (SSM) model with the Differential Evolution Adaptive Metropolis (DREAM) Markov Chain Monte Carlo (MCMC) algorithm with ZS enhancements (DREAM-ZS).
54 ENVIRONMENTAL SCIENCES↗
Techniques and resources for storm-scale numerical weather prediction
The topics discussed include the following: multiscale application of the 5th-generation PSU/NCAR mesoscale model, the coupling of nonhydrostatic atmospheric and hydrostatic ocean models for air-sea interaction studies; a numerical simulation of cloud formation over complex topography; adaptive grid simulations of convection; an unstructured grid, nonhydrostatic meso/cloud scale model; efficient mesoscale modeling for multiple scales using variable resolution; initialization of cloud-scale models with Doppler radar data; and making effective use of future computing architectures, networks, and visualization software.
Droegemeier, Kelvin↗
SULI Intern Final Report: Computationally Investigating Hydrogen Thermo-Diffusion in Yttrium Hydride Using Multiscale Methods
The renaissance of nuclear energy has arrived, heralding an age of abundant inexpensive clean energy, and renewed space exploration. In nuclear-powered spacecraft and microreactors, safety and size are of utmost importance. Yttrium Hydride (YHx) is being researched for its utility as a neutron moderator in nuclear reactors; the hydrogen in YHx slows down neutrons, enabling a continuous nuclear reaction in the reactor. This has the benefit of allowing reactors to be more safe, compact, and efficient. The goal of this effort is to computationally predict the coefficient of temperature-dependent hydrogen diffusion within YHx, the Soret coefficient. This parameter is essential for determining the safe operating modes of YHx moderators. Zirconium Hydride (ZrHx) is used in the Training, Research, Isotopes, General Atomics (TRIGA) reactor, is the reference material for these calculations. In this work, nanoscale atomic modeling in the Vienna Ab initio Simulation Package (VASP) is combined with the mesoscale finite element phase-field module in the Multiphysics Object-Oriented Simulation Environment (MOOSE); this culminates in a new multiscale computational method to simulate Soret diffusion of hydrogen in YHx. This data is useful for predicting experimental outcomes. This workflow involves convergence testing followed by static, Nudged Elastic Band (NEB), Quasi-Harmonic Approximation (QHA), and Molecular Dynamics (MD) calculations - linked with phase field simulation. NEB simulates hydrogen migration, while QHA and MD predict temperature-dependent properties. The static calculations align with literature, and preliminary NEB and QHA calculations yield accurate results. Once the atomic calculations are complete, we will incorporate Electron Backscatter Diffraction (EBSD) images and VASP-generated parameters into the phase field module to simulate intra- and intergranular transport of hydrogen in ZrHx and YHx. Future research will extend our approach to fuel-moderator materials systems such as Uranium-Yttrium Hydride (U-YHx). This work contributes to the development of advanced nuclear energy solutions for space travel.
36 - MATERIALS SCIENCE↗
Sea Ice Rheology Experiment (SIREx): 2. Evaluating Linear Kinematic Features in High-Resolution Sea Ice Simulations
Simulating sea ice drift and deformation in the Arctic Ocean is still a challenge because of the multiscale interaction of sea ice floes that compose the Arctic Sea ice cover. The Sea Ice Rheology Experiment (SIREx) is a model intercomparison project of the Forum of Arctic Modeling and Observational Synthesis (FAMOS). In SIREx, skill metrics are designed to evaluate different recently suggested approaches for modeling linear kinematic features (LKFs) to provide guidance for modeling small-scale deformation. These LKFs are narrow bands of localized deformation that can be observed in satellite images and also form in high resolution sea ice simulations. In this contribution, spatial and temporal properties of LKFs are assessed in 36 simulations of state-of-the-art sea ice models and compared to deformation features derived from the RADARSAT Geophysical Processor System. All simulations produce LKFs, but only very few models realistically simulate at least some statistics of LKF properties such as densities, lengths, or growth rates. All SIREx models overestimate the angle of fracture between conjugate pairs of LKFs and LKF lifetimes pointing to inaccurate model physics. The temporal and spatial resolution of a simulation and the spatial resolution of atmospheric boundary condition affect simulated LKFs as much as the model's sea ice rheology and numerics. Only in very high resolution simulations (≤2 km) the concentration and thickness anomalies along LKFs are large enough to affect air-ice-ocean interaction processes.
54 ENVIRONMENTAL SCIENCES↗
Advances in Multiscale Modeling of Lignocellulosic Biomass
Applications and associated processing technologies of lignocellulosic biomass are becoming more important as we endeavor to meet societal demand for fuels, chemicals, and materials from renewable resources. Meanwhile, the rapidly expanding availability and capabilities of high-performance computing present an unprecedented opportunity to accelerate development of technologies surrounding lignocellulose utilization. In order to realize this potential, suitable modeling frameworks must be constructed that effectively capture the multiscale complexity and tremendous variety exhibited by lignocellulosic materials. In our assessment of previous endeavors towards this goal, several important shortcomings have been identified: (1) the lack of multiscale integration strategies that capture emergent properties and behaviors spanning different length scales, and (2) the inability of many modeling approaches to effectively capture the variability and diversity of lignocellulose that arise from both natural and process-induced sources. In this perspective, we survey previous modeling approaches for lignocellulose and simulations processes involving its chemical and mechanical transformation, and suggest opportunities for future development to enhance the utility of computational tools to address barriers to widespread adoption of a renewable bioeconomy.
09 BIOMASS FUELS↗
The development of a high-resolution Eulerian radiation-hydrodynamics simulation capability for laser-driven Hohlraums
Hohlraums are hollow cylindrical cavities with high-Z material walls used to convert laser energy into uniform x-ray radiation drives for inertial confinement fusion capsule implosions and high energy density physics experiments. Credible computational modeling of hohlraums requires detailed modeling and coupling of laser physics, hydrodynamics, radiation transport, heat transport, and atomic physics. We report on improvements to Los Alamos National Laboratory's xRAGE radiation-hydrodynamics code in order to enable hohlraum modeling. xRAGE's Eulerian hydrodynamics and adaptive mesh refinement make it uniquely well suited to study the impacts of multiscale features in hohlraums. In order to provide confidence in this new modeling capability, we demonstrate xRAGE's ability to produce reasonable agreement with data from several benchmark hohlraum experiments. We also use xRAGE to perform integrated simulations of a recent layered high density carbon capsule implosion on the National Ignition Facility in order to evaluate the potential impacts of the capsule support tent, mixed cell conductivity methodologies, plasma transport, and cross-beam energy transfer (XBT). We find that XBT, seeded by plasma flows in the laser entrance hole (LEH), causes a slight decrease in energy coupling to the capsule and that all of these impact the symmetry of the x-ray drive such that they have an appreciable impact on the capsule implosion shape.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Optimal renormalization of multiscale systems
While model order reduction is a promising approach in dealing with multi-scale time-dependent systems that are too large or too expensive to simulate for long times, the resulting reduced order models can suffer from instabilities. We have recently developed a time-dependent renormalization approach to stabilize such reduced models. In the current work, we extend this framework by introducing a parameter that controls the time-decay of the memory of such models and optimally selecting this parameter based on limited fully resolved simulations. First, we demonstrate our framework on the inviscid Burgers equation whose solution develops a finite-time singularity. Our renormalized reduced order models are stable and accurate for long times while using for their calibration only data from a full order simulation before the occurrence of the singularity. Furthermore, we apply this framework to the 3D Euler equations of incompressible fluid flow, where the problem of finite-time singularity formation is still open and where brute force simulation is only feasible for short times. Our approach allows us to obtain for the first time a perturbatively renormalizable model which is stable for long times and includes all the complex effects present in the 3D Euler dynamics. We find that, in each application, the renormalization coefficients display algebraic decay with increasing resolution, and that the parameter which controls the time-decay of the memory is problem-dependent.
Price, Jacob↗
Microstructure Modeling of 3rd Generation Disk Alloys
The objective of this program is to model, validate, and predict the precipitation microstructure evolution, using PrecipiCalc (QuesTek Innovations LLC) software, for 3rd generation Ni-based gas turbine disc superalloys during processing and service, with a set of logical and consistent experiments and characterizations. Furthermore, within this program, the originally research-oriented microstructure simulation tool will be further improved and implemented to be a useful and user-friendly engineering tool. In this report, the key accomplishment achieved during the second year (2008) of the program is summarized. The activities of this year include final selection of multicomponent thermodynamics and mobility databases, precipitate surface energy determination from nucleation experiment, multiscale comparison of predicted versus measured intragrain precipitation microstructure in quench samples showing good agreement, isothermal coarsening experiment and interaction of grain boundary and intergrain precipitates, primary microstructure of subsolvus treatment, and finally the software implementation plan for the third year of the project. In the following year, the calibrated models and simulation tools will be validated against an independently developed experimental data set, with actual disc heat treatment process conditions. Furthermore, software integration and implementation will be developed to provide material engineers valuable information in order to optimize the processing of the 3rd generation gas turbine disc alloys.
Jou, Herng-Jeng↗