Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale simulations”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Solitary Magnetic Structures at Quasi–Parallel Collisionless Shocks: Formation

Solitary magnetic structures known as SLAMS (short large-amplitude magnetic structures) have been considered as essential elements of collisionless shocks with quasi-parallel geometries. Yet the physics underlying their formation remains an open question. In this paper, we use measurements from the magnetospheric multiscale mission combined with fully kinetic simulations to study the formation of SLAMS. We find that gyro-resonance between solar wind ions and right-hand circularly polarized electromagnetic waves results in magnetic field amplification. Gyro-trapping by the growing magnetic field builds up the plasma density that further enhances the current and field. The solitary nature of SLAMS stems from a beat-like magnetic field envelope where the maximum sets the initial location for nonlinear growth. Our results present a conceptual advance on SLAMS, and may shed new light on the open question of magnetic field amplification at astrophysical shocks.

58 GEOSCIENCES↗
DOE OSTI2017

Strong attractions and repulsions mediated by monovalent salts

Significance The physical properties of proteins and polyions in electrolyte solutions are of crucial importance for understanding biological processes and for controlling biotechnological processes. The concentration of salt controls the range and strength of the interactions in these complex electrolytes. Here, multiscale coarse-grained molecular dynamics simulations reveal long-range interactions between nanoparticles in aqueous solutions at high monovalent salt concentrations. The simulations explain the crystallization of nanoparticles even when they are strongly charged and explain the reentrant attractive to repulsive to attractive interactions as the nanoparticle charge density increases in high monovalent salt-concentration conditions.

Li, Yaohua↗
DOE OSTI2019

Ultrafast X-Ray Scattering Measurements of Coherent Structural Dynamics on the Ground-State Potential Energy Surface of a Diplatinum Molecule

We present x-ray free electron laser experiments addressing ground-state structural dynamics of the diplatinum anion Pt 2 POP 4 following photoexcitation. The structural dynamics are tracked with < 100 fs time resolution by x-ray scattering, utilizing the anisotropic component to suppress contributions from the bulk solvent. The x-ray data exhibit a strong oscillatory component with period 0.28 ps and decay time 2.2 ps, and structural analysis of the difference signal directly shows this as arising from ground-state dynamics along the PtPt coordinate. Here, these results are compared with multiscale Born-Oppenheimer molecular dynamics simulations and demonstrate how off-resonance excitation can be used to prepare a vibrationally cold excited-state population complemented by a structure-dependent depletion of the ground-state population which subsequently evolves in time, allowing direct tracking of ground-state structural dynamics.

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

Enhanced Planar Antenna Efficiency Through Magnetic Thin-Films

This work proposes to use magnetic material as the substrate of planar antennas to overcome the platform effect caused by the conducting ground plane. The upper bound of the radiation efficiency of an electric-current-driven low-profile antenna is theoretically derived, which is inversely proportional to the Gilbert damping factor of the magnetic material. Meanwhile, the improvement of radiation due to the use of magnetic material is demonstrated by a three-dimensional (3D) multiphysics and multiscale time-domain model. The simulation results match the theoretical derivation, showing 25% radiation efficiency from a planar antenna backed by a FeGaB thin film with 2.56 m thickness. Furthermore, for conductive ferromagnetic materials, it is shown that the eddy current loss can be well suppressed by laminating the thin film into multiple layers. Here, the radiation efficiency of the modeled antenna with a conductive ferromagnetic substrate is improved from 2.2% to 11.8% by dividing the substrate into 10 layers, with a ferromagnetic material fill factor of 93%.

42 ENGINEERING↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2016

3D Representative Volume Element Reconstruction of Fiber Composites via Orientation Tensor and Substructure Features

To provide a seamless integration of manufacturing processing simulation and fiber microstructure modeling, two new stochastic 3D microstructure reconstruction methods are proposed for two types of random fiber composites: random short fiber composites, and Sheet Molding Compounds (SMC) chopped fiber composites. A Random Sequential Adsorption (RSA) algorithm is first developed to embed statistical orientation information into 3D RVE reconstruction of random short fiber composites. For the SMC composites, an optimized Voronoi diagram based approach is developed for capturing the substructure features of SMC chopped fiber composites. The proposed methods are distinguished from other reconstruction works by providing a way of integrating statistical information (fiber orientation tensor) obtained from material processing simulation, as well as capturing the multiscale substructures of the SMC composites.

Li, Yi↗
DOE OSTI2015

Summary of BISON Development and Validation Activities - NEAMS FY15 Report

The objective of the NEAMS ToolKit is to develop a “pellet-to-plant” simulation capability useful for predicting performance and safety for a broad range of nuclear reactor power systems. The NEAMS ToolKit has been organized into a Fuels Product Line (FPL) and a Reactor Product Line (RPL) and is modular in design. Within the FPL, a multiscale approach has been adopted in which simulations of fuel performance at the engineering scale are informed by material property and irradiation performance models developed from mesoscale simulations of microstructural evolution. The focus in this report is on development and validation of the engineering-scale fuel performance analysis tool within the FPL, which is BISON. This summary report contains an overview of work performed under the work package enti- tled “FY2015 NEAMS INL-Engineering Scale Fuel Performance & Interface with RPL Tools.” A first chapter identifies the specific FY-15 milestones, providing a basic description of the asso- ciated work and references to related detailed documentation. Where applicable, a representative technical result is provided. A second chapter summarizes major additional accomplishments including 1) a BISON validation paper, 2) a multiphysics coupling paper, 3) application of BI- SON to Reactivity Insertion Accident analysis and 4) use of BISON for uncertainity analysis in fuel modeling. A final chapter outlines FY-16 future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Micropolar deep material network

This study extends the Deep Material Network (DMN), a physics-informed machine learning framework, to predict the homogenized mechanical response of composite materials with micropolar (Cosserat-type) constitutive behavior. This extension incorporates microstructure-dependent size effects, enabling accurate, efficient, and size-aware predictions for composites with complex internal architectures. While traditional, direct numerical simulation micropolar models effectively capture size effects by introducing extra local degrees of freedom, they bring significant computational challenges, particularly for multiscale analyses relevant to engineering applications. The micropolar DMN developed in this paper achieves high accuracy while significantly reducing computation time compared to micropolar direct numerical simulations. This advancement enables multiscale analyses and parameter studies that were previously impractical, such as high-cycle fatigue simulations and comprehensive investigations of internal length scale effects notably in size-dependent plastic response and the optimization of lattice structures. By uniting microstructure-sensitive modeling, physics-driven learning, and scalable surrogate modeling, the micropolar DMN paves the way for accelerated material design, large-scale parametric studies, and the reliable incorporation of size-dependent effects across a wide range of engineering applications, including optimization and next-generation composite design.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Mesoscale Physicochemical Interactions in Lithium–Sulfur Batteries: Progress and Perspective

The shuttle effect and poor conductivity of the discharge products are among the primary impediments and scientific challenges for lithium–sulfur batteries. The lithium–sulfur battery is a complex energy storage system, which involves multistep electrochemical reactions, insoluble polysulfide precipitation in the cathode, soluble polysulfide transport, and self-discharge caused by chemical reactions between polysulfides and Li metal anode. These phenomena happen at different length and time-scales and are difficult to be entirely gauged by experimental techniques. In this paper, we reviewed the multiscale modeling studies on lithium–sulfur batteries: (1) the atomistic simulations were employed to seek alternative materials for mitigating the shuttle effect; (2) the growth kinetics of Li2S film and corresponding surface passivation were investigated by the interfacial model based on findings from atomistic simulations; (3) the nature of Li2S2, which is the only solid intermediate product, was revealed by the density functional theory simulation; and (4) macroscale models were developed to analyze the effect of reaction kinetics, sulfur loading, and transport properties on the cell performance. The challenge for the multiscale modeling approach is translating the microscopic information from atomistic simulations and interfacial model into the meso-/macroscale model for accurately predicting the cell performance.

Electrochemistry↗
DOE OSTI2018

Rapid freezing of water under dynamic compression

Understanding the behavior of materials at extreme pressures is a central issue in fields like aerodynamics, astronomy, and geology, as well as for advancing technological grand challenges such as inertial confinement fusion. Dynamic compression experiments to probe high-pressure states often encounter rapid phase transitions that may cause the materials to behave in unexpected ways, and understanding the kinetics of these phase transitions remains an area of great interest. We examine experimental and theoretical/computational efforts to study the freezing kinetics of water to a high-pressure solid phase known as ice VII. We first present a detailed analysis of dynamic compression experiments in which water has been observed to freeze on sub-microsecond time scales to ice VII. This is followed by a discussion of the limitations of currently available molecular and continuum simulation methods in modeling these experiments. We then describe how our phase transition kinetics models, which are based on classical nucleation theory, provide a more physics-based framework that overcomes some of these limitations. We give suggestions on future experimental and modeling work on the liquid–ice VII transition, including an outline of the development of a predictive multiscale model in which molecular and continuum simulations are intimately coupled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Develop a new integrated macro→micro←nano (MMN) multiscale modeling framework to optimize high strength aluminum alloys and processes for vehicle light-weighting​

Bending tests provide a means to study plane strain fracture performance of 6000 series aluminum alloys. Metrics from bending tests have been correlated with self-pierce riveting (SPR) performance of a high strength AA6111 automotive aluminum alloy in previous works. Using the ORNL HPC resources, this project developed an innovative macro→micro←nano (MMN) multiscale microstructure-based finite element (FE) code to further understanding of the relationship between microstructure and fracture properties of high-strength 6000-series alloys. This work started with microstructural characterization in both mesoscale and nanoscale and bending performance characterization of AA6111 HS2-T6 alloy at Ford, the MMN framework was applied to this alloy to simulate 3-point VDA bending. From the results of the macro-modeling of 3-point VDA bending, the critical region of fracture was identified, and the region geometry was used to construct the micro-model. The fracture criterion of micron-scale precipitates and aluminum matrix which contains submicron and nano particles (AL-SMP-NP) in the micro-model was calibrated and validated by comparing simulated and measured bending results. With the AL-SMP-NP fracture strain obtained, the fracture strain of Al-matrix containing nano particles (ALNP) will similarly be determined by a submicron scale-model using an edge-constrained FE modeling approach developed by Hu et al. With the ALNP fracture strain obtained, the fracture strain of Al-matrix containing no particles will similarly be determined by a nano scale-model using an edge-constrained FE modeling approach. After the MMN framework is built and fracture criterion calibrated, nano-model FE simulations with virtual microstructures was performed to obtain a reduced order model (ROM) of the fracture criterion of the ALNP as a function of volume fraction, size, and distribution of the nanoparticle. This nano→submicron→micro modeling part allows exploration of the influence of different material nanostructures from different process conditions on the bending properties within the multiscale bending simulation framework and the ROM of material bendability as a function of nanoparticle size and shape was established. This obtained reduced order model (ROM) could help guide the design and selection of materials to improve existing SPR process models that could replace trial-and-error rivet/die selection and help to design new rivet/die combinations capable of robustly joining new higher strength 6000 5 series alloys in automotive body structures. This would enable lightweighting of Ford vehicles leading to greater fuel efficiency and reduce manufacturing time and energy.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Multiscale electro-chemo-mechanical model of high-capacity amorphous silicon anode material in sodium-ion batteries

The growing interest in sodium-ion batteries (SIBs) is fueled by their abundant resources and environmentally friendly nature, with amorphous silicon (a-Si) emerging as a promising anode material for enhancing capacity. However, the key challenge lies in sustaining reversible capacity during cycling. Here, in this work, we developed a multiscale electrochemical model incorporating an a-Si anode to elucidate the performance parameters of SIBs. Additionally, we integrated an electro-chemo-mechanical model at the particle level to capture stress generation, an essential factor in the degradation of high-capacity electrodes. Unlike existing models, our approach accounts for large-deformation chemo-mechanics at the particle scale and includes simulations under varying charge rates to explore multiscale behavior. The results reveal that coupled sodiation significantly prolongs complete cycling times and influences discharge dynamics, indicating that neglecting this coupling leads to an underestimation of actual capacity. Furthermore, we observed pronounced polarization effects at higher charge rates, resulting in heterogeneous stress distributions across the anode. With the identification of critical failure parameters for both active particles and binder materials, offering novel insights for mitigating degradation in high-capacity electrode systems.

Amorphous silicon anode↗
DOE OSTI2025

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗
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 OSTI2019

Data-driven material models for atomistic simulation

The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parametrizations for potentials using traditional approaches. Machine learning has emerged as an effective alternative approach to develop accurate and robust interatomic potentials. Beginning with a very general model form, the potential is learned directly from a database of electronic structure calculations and therefore can be viewed as a multiscale link between quantum and classical atomistic simulations. Risk of inaccurate extrapolation exists outside the narrow range of time and length scales where the two methods can be directly compared. Here, we use the spectral neighbor analysis potential (SNAP) and show how a fit can be produced with minimal interpolation errors which is also robust in extrapolating beyond training. To demonstrate the approach, we have developed a tungsten-beryllium potential suitable for the full range of binary compositions. Subsequently, large-scale molecular dynamics simulations were performed of high energy Be atom implantation onto the (001) surface of solid tungsten. The machine learned W-Be potential generates a population of implantation structures consistent with quantum calculations of defect formation energies. A very shallow ( < 2 nm ) average Be implantation depth is predicted which may explain ITER diverter degradation in the presence of beryllium.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Constraint energy minimizing generalized multiscale finite element method for multi-continuum Richards equations

In fluid flow simulation, the multi-continuum model is a useful strategy. When the heterogeneity and contrast of coefficients are high, the system becomes multiscale, and some kinds of reduced order methods are demanded. Combining these techniques with nonlinearity, we will consider in this paper a dual-continuum model which is generalized as a multi-continuum model for a coupled system of nonlinear Richards equations as unsaturated flows, in complex heterogeneous fractured porous media; and we will solve it by a novel multiscale approach utilizing the constraint energy minimizing generalized multiscale finite element method (CEM-GMsFEM). In particular, such a nonlinear system will be discretized in time and then linearized by Picard iteration (whose global convergence is proved theoretically). Subsequently, we tackle the resulting linearized equations by the CEM-GMsFEM and obtain proper offline multiscale basis functions to span the multiscale space (which contains the pressure solution). More specifically, we first introduce two new sources of samples, and the GMsFEM is used over each coarse block to build local auxiliary multiscale basis functions via solving local spectral problems, that are crucial for detecting high-contrast channels. Second, per oversampled coarse region, local multiscale basis functions are created through the CEM as constrainedly minimizing an energy functional. Various numerical tests for our approach reveal that the error converges with the coarse-grid size and that only few oversampling layers as well as basis functions are needed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Hierarchical top-down bottom-up calibration with consideration for uncertainty and inter-scale discrepancy of Peierls stress of bcc Fe

A hierarchical multiscale model of plasticity in single crystal bcc Fe, framed at the mesoscale pertaining to a statistically representative set of dislocations for each slip system, relies on specification of a set of parameters which are informed using both bottom-up (BU) (atomistic) and top-down (TD) (experimental) information. We term this specification process a 'connection' between BU and TD pathways to inform the mesoscale model. The connection is considered in the presence of error, uncertainty, and discrepancy between the models. Furthermore, we expand upon a previously developed reconciled TD and BU calibration method to account for anticipated discrepancy between information from different length scales. The results of a previously formulated likelihood-based connection test of the multiscale model suggest a 'missing link' may be responsible for part of the inter-scale discrepancy. In this case, this link is assumed to be a relation between the unit-process (single dislocation line) of coordinated kink-pair nucleation on a screw dislocation segment and the many-body dislocation process which manifests in the onset of experimentally observed plastic deformation in a single crystal. A physics-informed discrepancy layer is formulated to improve the connection in the presence of additional persistent uncertainties. This physics-informed hypothesis testing is demonstrated as an alternative to fully data-driven search methods, particularly applicable to datasets of limited size that are often encountered in such multiscale material modeling applications. Discrete dislocation dynamics simulations are used to investigate the effect of the Peierls stress on straight screw dislocations and in the operation of Frank–Read sources. Discussion concerns the critical role of uncertainty quantification in providing a basis for this form of incremental hypothesis testing.

42 ENGINEERING↗