Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale 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.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Difference in neoclassical edge flows between strongly negative and positive triangularities in the XGC gyrokinetic simulation

The neoclassical baseline study of a strongly negative triangularity (NT) plasma and the corresponding positive triangularity (PT) plasma is performed using the edge-specialized, total-f gyrokinetic code (X-point Gyrokinetic Code, XGC). A DIII-D-like plasma is used, based on the negative triangularity discharge of DIII-D #193 793. An artificial positive triangularity (PT) equilibrium has been constructed to compare the edge rotation physics at the same triangularity strength but with opposite sign, while keeping the same elongation and other geometric parameters. Carbon(6+) ions are added to the deuterium plasma at an experimentally relevant level. By using the experimental profile of carbon toroidal rotation profile as an input, XGC finds that the deuteron rotation is significantly different from the carbon rotation at the inboard and outboard midplanes, mostly caused by the difference in the Pfirsch–Schlüter rotation. More importantly, a significant difference in the X-point orbit loss physics, thus the rotation source, is found between the positive and negative triangularity equilibrium models. However, it is also found that the agreement between the present neoclassical simulation and the experimental NT data is validated only within the middle of the pedestal slope, indicating the importance of edge turbulence. Finally, this study could establish a baseline for the multiphysics, multiscale studies that include turbulence of negative triangularity plasmas.

Deuterium↗
DOE OSTI2021

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

Computer Simulation of Proton Transport in Fuel Cell Membranes (Final Report)

This DOE-supported research grant focused on understanding the nature of proton transport in complex systems such as proton exchange membranes (PEMs). The most unique aspect of the research was the development and implementation of a novel multiscale reactive molecular dynamics (MS-RMD) methodology. In this approach, covalent bonds can dynamically break and form, allowing one to accurately treat the proton hopping process essential to capturing the physics of proton transport. The Voth group applied this method to proton exchange membrane systems, providing insight into their proton transport mechanism. They found that protons can diffuse most rapidly in the water-rich regions, but that protons actually spend so little time in such regions that transport along the hydrophobic – hydrophilic interface controls the membrane performance. In addition, the group worked to increase understanding of acidic solutions, developing methods for simulating and interpreting experimental infrared vibrational spectroscopy for excess protons (acidic solutions). The group also implemented novel tools for developing proton transport reactive MD models using a relative entropy minimization scheme.

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

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

Ab initio generalized Langevin equation

We introduce a machine learning–based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables (CVs) in materials and molecules. In this scheme, the parameters are learned from atomistic simulations based on ab initio quantum mechanical models. Force field, memory kernel, and noise generator are constructed in the context of the Mori–Zwanzig formalism, under the constraint of the fluctuation–dissipation theorem. Combined with deep potential molecular dynamics and electronic density functional theory, this approach opens the way to multiscale modeling in a variety of situations. Here, we demonstrate this capability with a study of two mesoscale processes in crystalline lead titanate, namely the field-driven dynamics of a planar ferroelectric domain wall, and the dynamics of an extensive lattice of coarse-grained electric dipoles. In the first case, AIGLE extends the reach of ab initio simulations to a regime of noise-driven motions not accessible to molecular dynamics. In the second case, AIGLE deals with an extensive set of CVs by adopting a local approximation for the memory kernel and retaining only short-range noise correlations. The scheme is computationally more efficient than molecular dynamics by several orders of magnitude and mimics the microscopic dynamics at low frequencies where it reproduces accurately the dominant far-infrared absorption frequency.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Manifold learning for coarse-graining atomistic simulations: Application to amorphous solids

In this work, we introduce a generalized machine learning framework to probabilistically parameterize upper-scale models in the form of nonlinear PDEs consistent with a continuum theory, based on coarse-grained atomistic simulation data of mechanical deformation and flow processes. The proposed framework utilizes a hypothesized coarse-graining methodology with manifold learning and surrogate-based optimization techniques. Coarse-grained high-dimensional data describing quantities of interest of the multiscale models are projected onto a nonlinear manifold whose geometric and topological structure is exploited for measuring behavioral discrepancies in the form of manifold distances. A surrogate model is constructed using Gaussian process regression to identify a mapping between stochastic parameters and distances. Derivative-free optimization is employed to adaptively identify a unique set of parameters of the upper-scale model capable of rapidly reproducing the system's behavior while maintaining consistency with coarse-grained atomic-level simulations. The proposed method is applied to learn the parameters of the shear transformation zone (STZ) theory of plasticity that describes plastic deformation in amorphous solids as well as coarse-graining parameters needed to translate between atomistic and continuum representations. We show that the methodology is able to successfully link coarse-grained microscale simulations to macroscale observables and achieve a high-level of parity between the models across scales.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Investigation of flow regimes in arc plasma–gas interactions using a two-temperature arc in crossflow model

The perpendicular impingement of a gas stream on an electric arc, a configuration known as the arc in crossflow, is of primary relevance in the study of plasma-gas interactions, as well as in industrial applications such as circuit breakers and wire-arc spraying. The flow dynamics in the arc in crossflow are the result of coupled fluid-thermal-electromagnetic phenomena accompanied by large property gradients, which can produce significant deviations from Local Thermodynamic Equilibrium (LTE) among electrons and gas species. These characteristics can lead to the establishment of distinct flow regimes depending on the relative values of the controlling parameters of the system, such as inflow velocity, arc current, and inter-electrode spacing. A two-temperature non-LTE (NLTE) model is used to investigate the arc dynamics and the establishment of flow regimes in the arc in crossflow. The plasma flow model is implemented within a nonlinear Variational Multiscale (VMS) numerical discretization approach that is less dissipative, and hence better suited to capture unstable behavior, than traditional VMS methods commonly used in computational fluid dynamics simulations. The Reynolds and the Enthalpy dimensionless numbers, characterizing the relative flow strength and arc strength, respectively, are chosen as the controlling parameters of the system. Simulation results reveal the onset of dynamic behavior and the establishment of steady, periodic, quasi-periodic, and chaotic or potentially turbulent, regimes, as identified by distinct spatiotemporal fluctuations. The computational findings reveal the role of increasing the relative arc strength on enhancing flow stability by delaying the growth of fluctuating and unstable flow behavior.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Theory and simulations of linear and nonlinear two-dimensional Rayleigh–Taylor dynamics with variable acceleration

Interfacial Rayleigh–Taylor mixing is crucial to describing important natural and engineering processes, such as exploding supernovae, laser micromachining, hot spots in inertial confinement fusion, and optical telecommunications. These require the characterization of the time dependence of the driving acceleration. We compare our theoretical formulation based on group theory foundations with interface-capturing numerical simulations for linear and nonlinear two-dimensional Rayleigh–Taylor instabilities in a finite-sized domain with time-varying acceleration over broad ranges of Atwood numbers and acceleration exponents. Detailed corroboration between theory and simulations is provided for this foundational case. Both demonstrate the strong interfacial nature of Rayleigh–Taylor instabilities, which suggests that practical flow fields can be reconstructed from the derived fluid potential using the proposed theory. A robust agreement is also obtained for the early and late-time evolution of the amplitudes of the bubble and spike, which demonstrate that the Rayleigh–Taylor flow can transition to the mixing regime even for a single-mode initial perturbation. Corroboration with experiments of high energy density plasmas motivated by studies of supernovae is also achieved. In addition, a long-standing puzzle in Rayleigh–Taylor dynamics on the interplay between the acceleration, the shear, and the interface morphology in the theory and simulations is resolved by accounting for finite viscosity of the fluids. The characterization of Rayleigh–Taylor instabilities as a highly interfacial phenomenon provides valuable insight into its multiscale nature, which enhances the design and understanding of numerous processes of practical interest.

42 ENGINEERING↗
DOE OSTI2022

Machine learning predictions for local electronic properties of disordered correlated electron systems

We present a scalable machine learning (ML) model to predict local electronic properties such as on-site electron number and double occupation for disordered correlated electron systems. Our approach is based on the locality principle, or the nearsightedness nature, of many-electron systems, which means local electronic properties depend mainly on the immediate environment. A ML model is developed to encode this complex dependence of local quantities on the neighborhood. We demonstrate our approach using the square-lattice Anderson-Hubbard model, which is a paradigmatic system for studying the interplay between Mott transition and Anderson localization. We develop a lattice descriptor based on the group-theoretical method to represent the on-site random potentials within a finite region. The resultant feature variables are used as input to a multilayer fully connected neural network, which is trained from data sets of variational Monte Carlo (VMC) simulations on small systems. We show that the ML predictions agree reasonably well with the VMC data. Our work underscores the promising potential of ML methods for multiscale modeling of correlated electron systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Coalescence dynamics of prismatic dislocation loops due to vacancy supersaturation

Recovery processes in metals with a high defect concentration are at the heart of a structural alloy's usability in engineering applications. However, such processes take place via nonconservative mechanisms mediated by point defect diffusion and are governed by their interactions with the underlying microstructure. Indeed, one of the main challenges in studying recovery processes is their intrinsic multiscale complexity, which often negates the direct use of experimental techniques. Here, in this Research Letter, we study the mechanism of dislocation loop coalescence in metals by simulating the interaction dynamics between vacancies and prismatic loops using a stochastic framework that seamlessly bridges point defect kinetics with dislocation mechanics. Our simulations reveal two fundamental discoveries: first, that climb processes are governed by local fluctuations that are not captured by continuum elastodiffusion models, and, second, that prismatic loops create internal traps for vacancies, leading to the formation of rings around the perimeter of the loops that affect how coalescence takes place. Our results show an excellent agreement with in situ experimental observations of loop coalescence in Fe under irradiation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

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

On the Representation of Hyporheic Exchange in Models for Reactive Transport in Stream and River Corridors

Efforts to include more detailed representations of biogeochemical processes in basin-scale water quality simulation tools face the challenge of how to tractably represent mass exchange between the flowing channels of streams and rivers and biogeochemical hotspots in the hyporheic zones. Multiscale models that use relatively coarse representations of the channel network with subgrid models for mass exchange and reactions in the hyporheic zone have started to emerge to address that challenge. Two such multiscale models are considered here, one based on a stochastic Lagrangian travel time representation of advective pumping and one on multirate diffusive exchange. The two models are formally equivalent to well-established integrodifferential representations for transport of non-reacting tracers in steady stream flow, which have been very successful in reproducing stream tracer tests. Despite that equivalence, the two models are based on very different model structures and produce significantly different results in reactive transport. In a simple denitrification example, denitrification is two to three times greater for the advection-based model because the multirate diffusive model has direct connections between the stream channel and transient storage zones and an assumption of mixing in the transient storage zones that prevent oxygen levels from dropping to the point where denitrification can progress uninhibited. By contrast, the advection-based model produces distinct redox zonation, allowing for denitrification to proceed uninhibited on part of the hyporheic flowpaths. These results demonstrate that conservative tracer tests alone are inadequate for constraining representation of mass transfer in models for reactive transport in streams and rivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

SQMS Quantum R&D in Machine Learning, Optimization and Sensing beyond Fundamental Physics Applications

This newly formed team at SQMS under the Ecosystem Thrust is looking to develop capabilities impacting societal advances outside the core domain of HEP and condensed matter physics. We explicitly leverage the experimental and algorithmic innovations developed across all groups as well as connect to broad-scope external projects of the diverse team of PIs. As the inaugural set of projects, we are studying numerically quantum machine learning models inspired by efficiently trainable echo-state and orthogonal neural networks and developing designs for related experiments to be performed on quantum processors based on SQMS SRF cQED technology and Rigetti s transmon arrays. Investigated models exploit ideas and lessons learned from multiple prior work by SQMS team members in a variety of internal and external activities [R1]. Target initial applications include noisy signal processing, potentially captured by quantum sensors or noisy QPUs, as well as simulation and classification of healthcare data. For instance, image reconstruction of the brain s electrical properties by solving the inverse Maxwell equation problem with uncertainty [R2] through a hybrid quantum-classical physics-informed architecture for time-dependent processes [R3]. The group is also investigating the application and development of novel quantum sensors based on magnetic levitation of a superconducting sphere coupled to a superconducting qubit. This coupling enables high-precision measurements of the position of the sphere, which can be used for sensitive detection of forces, enabling practical applications such as gravimetry for geophysics analysis, or accelerometry for GPS-denied navigation [R4] [R1] Rieffel, Eleanor G., Ata Akbari Asanjan, M. Sohaib Alam, Namit Anand, David E. Bernal Neira, Sophie Block, Lucas T. Brady et al. "Assessing and advancing the potential of quantum computing: A NASA case study." Future Generation Computer Systems (2024). [R2] Yu, X., Serrall s, J.E., Giannakopoulos, I.I., Liu, Z., Daniel, L., Lattanzi, R. and Zhang, Z., 2023. Pifon-ept: Mr-based electrical property tomography using physics-informed fourier networks. IEEE Journal on Multiscale and Multiphysics Computational Techniques. [R3] Wudarski, Filip, Daniel OConnor, Shaun Geaney, Ata Akbari Asanjan, Max Wilson, Elena Strbac, P. Aaron Lott, and Davide Venturelli. "Hybrid quantum-classical reservoir computing for simulating chaotic systems." arXiv preprint arXiv:2311.14105 (2023). [R4] Higgins, Gerard, Saarik Kalia, and Zhen Liu. "Maglev for dark matter: Dark-photon and axion dark matter sensing with levitated superconductors." Physical Review D 109.5 (2024): 055024.

Venturelli, Davide↗
DOE OSTI2025

Effective Li-Ion Transport Quantification in Composite Cathodes for All-Solid-State Batteries via Multiscale Modeling and Experiments

The tortuosity factor of composite cathodes significantly affects the rate performance of all-solid-state batteries (ASSBs) and has significant differences from systems with liquid electrolytes. Here, in this work, we report a simulation-experiment combined approach that quantifies the effective Li-ion transport in an ASSB composite cathode, which links tortuosity factor on ∼ μm scale to terminal voltage during cycling at the cell level (on ∼ cm scale). Two independent approaches of tortuosity factor quantification are considered: fitting electrochemical cycling data and verifying at different cycling rates and calculating from segmented tomography images, with the tortuosity factor quantified from both methods reaching self-consistency. The simulated terminal voltage using the quantified tortuosity factor has a small relative error of <3% compared to the experimental measurements. We find a significantly reduced value of the Bruggeman exponent of the catholyte phase (1.75), and using shape analysis, we show that rod-shaped catholyte particles play an important role in lowering the tortuosity factor.

Yao, Archie Mingze [Univ. of Michigan, Ann Arbor, ↗
DOE OSTI2023

Direct comparison between experiments and dislocation dynamics simulations of high rate deformation of single crystal copper

A long standing challenge in computational materials science is to establish a quantitative connection between the macroscopic properties of plastic deformation with the microscopic mechanisms of dislocations in crystalline materials. Although the discrete dislocation dynamics (DDD) simulation method has been developed for several decades with the goal of addressing this challenge, a one-to-one comparison between the DDD predictions on single crystal stress–strain curves and experimental measurements under identical conditions has not been possible to date. Such a comparison is an essential step towards establishing a dislocation-physics based theory of plasticity and a multiscale framework of the plastic behaviors of crystalline materials. Here we provide direct comparisons between the stress–strain curves of Cu single crystals under high strain rate loading in the [0 0 1] and [0 1 1] directions obtained from miniaturized desktop Kolsky bar experiments and those from DDD simulations under identical loading conditions. With an appropriate set of parameters, DDD simulations can produce stress–strain curves that are in reasonable agreement with the experimental results. However, the dislocation mobility values needed to achieve this agreement are an order of magnitude lower than expected based on previous measurements and atomistic simulations. We hypothesize that this discrepancy could be caused by drag forces from jogs and point defects produced during the plastic deformation. Cross-slip of screw dislocations is also found to be necessary to capture the experimental stress–strain behavior, especially for the [0 1 1] loading direction. Finally, this work provides an example of how direct comparisons between DDD simulations and experimental measurements can provide new insight into the fundamental mechanisms of plastic deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Permeation of CO 2 and N 2 through glassy poly(dimethyl phenylene) oxide under steady- and presteady-state conditions

Glassy polymers are often used for gas separations because of their high selectivity. Although the dual-mode permeation model correctly fits their sorption and permeation isotherms, its physical interpretation is disputed, and it does not describe permeation far from steady state, a condition expected when separations involve intermittent renewable energy sources. To develop a more comprehensive permeation model, we combine experiment, molecular dynamics, and multiscale reaction–diffusion modeling to characterize the time-dependent permeation of N 2 and CO 2 through a glassy poly(dimethyl phenylene oxide) membrane, a model system. Simulations of experimental time-dependent permeation data for both gases in the presteady-state and steady-state regimes show that both single- and dual-mode reaction–diffusion models reproduce the experimental observations, and that sorbed gas concentrations lag the external pressure rise. The results point to environment-sensitive diffusion coefficients as a vital characteristic of transport in glassy polymers.

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

Deterioration of concrete due to ASR: Experiments and multiscale modeling

The process of ASR (Alkali-Silica Reaction) induced expansion and damage in pavement concrete specimens is investigated using laboratory experiments and computational modeling. In the experimental program, the concrete specimens are subject to CS-CPT (climate simulation concrete prism test) to obtain ASR induced expansion with and without external supply of alkali. The dissolution rates of the granodiorite used in the concrete mix and the gel formation rates are determined under concrete-like conditions (pH 13.8, with/without Ca(OH){sub 2} and NaCl) at different temperatures. A micromechanics based computational model with aggregate-scale diffusion and reaction kinetics coupled to an Eigenstrain based micromechanics damage model is developed for the simulation of ASR induced expansion and damage. Data from the experimental program are used to calibrate and validate the computational model. Model predictions show that for the given concrete mixture, ASR induced expansion in the specimen exposed to water is predominantly governed by microcracking in the aggregate, while the expansion in the specimen subjected to external alkali supply is governed by microcracking in both the aggregates and the cement paste.

36 MATERIALS SCIENCE↗