Search NASA⌕ Search

SEARCH · Search NASA

Results for “material response 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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Dynamics of deuterium retention and desorption from plasma-facing materials in fusion reactor-relevant conditions

Hydrogen isotopes retention and desorption during and after discharges in fusion devices are still not well understood due to the complex device conditions and limitations of in-situ diagnostics and measurements. Here, we simulated well-diagnosed recent experiments at the DIII-D facility to benchmark our ITMC-DYN integrated package of modeling deuterium diffusion, retention, and desorption during and after D discharge irradiation. Modeling results were compared with detail experimental data of D desorption fluxes for various irradiation conditions. We predicted the temporal evolution of free and trapped D distribution in tungsten (W) plasma-facing material (PFM). Effects of key parameters namely diffusion coefficient, recombination rate, trapping energies against different defect types, were examined in these simulations. Existing experimental data of these parameters in literature varies significantly which makes it harder to identify key mechanisms and physics responsible for hydrogen isotope retention and desorption. The purpose of this work is to accurately simulate recent well-diagnosed reactor experiments given the uncertainties in such parameters and identify mechanisms responsible for the retention and desorption. We implemented the best identified diffusion, recombination, and trapping parameters in ITMC-DYN package that integrate both various collisional and thermal processes. We predicted, for example, that sample cooling between discharges in DIII-D operations can significantly affect the spatial distribution of trapped D in W under reactor irradiation conditions. Correct prediction of desorption spectra from samples irradiated during 10 DIII-D discharges showed that up to 35% of D can be retained in high binding energy defects such as vacancy clusters or voids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Grain2Mesh: Mesh Generation for Grain-Scale Nonlinear Elasticity Modeling

The nonlinear hysteretic behavior of rocks under cyclic loading is a crucial area of study in geomechanics. The macroscopic response of a variety of materials has been found to be contingent upon the behavior of the micro-scale structure. This project aims to develop a functional and maintainable software package for generating a multi-phase numerical mesh and accompanying simulation files for finite element modeling used in computational mechanics solvers. Meshes generated from images often lack key preprocessing that reduces noise and prevents mesh element distortion that can increase computational cost. By incorporating user feedback throughout, grain2mesh ensures a high-fidelity mesh that can be used to model grain-scale interactions such as shearing, crack propagation, and interfacial material contrast. Scientific applications of this software include material fracturing, stress-strain analysis for natural and engineered materials, and nonlinear meso-scale analysis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗
DOE OSTI2025

Nuclear Astrophysics through simulations of neutron star mergers using Monte‐Carlo neutrino radiation transport (DE-SC0020435 Final Technical Report)

Neutron star mergers are an important source of information for nuclear physics. Except for black holes, neutron stars are the densest macroscopic objects known to exist in the Universe. They provide us with a remarkable laboratory to study the poorly understood high-density states of matter, complementing the study of heavy nuclei performed by nuclear physics experiments on Earth. Collisions of neutron stars are of particular interest. They are among the most energetic events observable in the Universe, powering a broad range of signals across the electromagnetic spectrum as well as gravitational wave signals. These signals contain important information about the properties of extremely dense matter. Colliding neutron stars additionally eject large amounts of neutron-rich material into the surrounding interstellar medium -- material that then undergoes rapid neutron-capture (r-process) nucleosynthesis, the mechanism thought to be responsible for the production of about half of the heavy nuclei. In order to understand current and future observations of neutron star mergers, we need reliable models for the signals that they power. An important component in the construction of such models are numerical simulations of colliding neutron stars. These simulations are costly, running for multiple months on supercomputers, and require the inclusion of complex physics (general relativity, magnetohydrodynamics, neutrino physics, nuclear reactions). In this document, we describe the development of new methods for the treatment of neutrinos in merger simulations, as well as the first simulations capable of evolving the equations of neutrino radiation transport directly and their impact on our understanding of neutron star mergers and more broadly in nuclear astrophysics. We also discuss additional work performed in the study of nucleosynthesis and neutrino physics in neutron star mergers as part of Early Career Award DE-SC0020435.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A Continuum Dislocation Dynamics Crystal Plasticity Approach to Irradiated Body-Centered Cubic α-Iron

Radiation-induced embrittlement of reactor pressure vessel (RPV) steels can potentially limit the operating life of nuclear power plants. Over extended exposure to radiation doses, these body-centered cubic (BCC) irons demonstrate irradiation damage. Here, we present a continuum dislocation density (CDD) crystal plasticity model to capture the interaction among dislocations and self-interstitial atom (SIA) loops in α-iron. We demonstrate the importance of modeling cross slip using a combined stochastic Monte Carlo approach and the role of slip system strength anisotropy in capturing stochastic cross slip interactions. Through these captured interactions, the CDD crystal plasticity model can capture both the stress response and the physical evolution of dislocations on different slip system planes. Single-crystal verification experiments are used to calibrate the CDD crystal plasticity model, and a set of simplified polycrystalline simulations demonstrates the model’s ability to capture the stress response from tensile experiments on α-iron.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fast and accurate reduced-order modeling of a MOOSE-based additive manufacturing model with operator learning

One predominant challenge in additive manufacturing (AM) is to achieve specific material properties by manipulating manufacturing process parameters during the runtime. Such manipulation tends to increase the computational load imposed on existing simulation tools employed in AM. The goal of the present work is to construct a fast and accurate reduced-order model (ROM) for an AM model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, ultimately reducing the time/cost of AM control and optimization processes. Our adoption of the operator learning (OL) approach enabled us to learn a family of differential equations produced by altering process variables in the laser’s Gaussian point heat source. More specifically, we used the Fourier neural operator (FNO) and deep operator network (DeepONet) to develop ROMs for time-dependent responses. Furthermore, we benchmarked the performance of these OL methods against a conventional deep neural network (DNN)-based ROM. Ultimately, we found that OL methods offer comparable performance and, in terms of accuracy and generalizability, even outperform DNN at predicting scalar model responses. The DNN-based ROM afforded the fastest training time. Furthermore, all the ROMs were faster than the original MOOSE model yet still provided accurate predictions. FNO had a smaller mean prediction error than DeepONet, with a larger variance for time-dependent responses. Unlike DNN, both FNO and DeepONet were able to simulate time series data without the need for dimensionality reduction techniques. Finally, the present work can help facilitate the AM optimization process by enabling faster execution of simulation tools while still preserving evaluation accuracy.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Quasi-Distributed Fiber Sensor-Based Approach for Pipeline Health Monitoring: Generating and Analyzing Physics-Based Simulation Datasets for Classification

This study presents a framework for detecting mechanical damage in pipelines, focusing on generating simulated data and sampling to emulate distributed acoustic sensing (DAS) system responses. The workflow transforms simulated ultrasonic guided wave (UGW) responses into DAS or quasi-DAS system responses to create a physically robust dataset for pipeline event classification, including welds, clips, and corrosion defects. This investigation examines the effects of sensing systems and noise on classification performance, emphasizing the importance of selecting the appropriate sensing system for a specific application. The framework shows the robustness of different sensor number deployments to experimentally relevant noise levels, demonstrating its applicability in real-world scenarios where noise is present. Overall, this study contributes to the development of a more reliable and effective method for detecting mechanical damage to pipelines by emphasizing the generation and utilization of simulated DAS system responses for pipeline classification efforts. The results on the effects of sensing systems and noise on classification performance further enhance the robustness and reliability of the framework.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High Energy Density Physics of Inertial Confinement Fusion Ablator Materials

The historic December 5, 2022 experiment at Lawrence Livermore National Lab’s (LLNL) National Ignition Facility (NIF) reached fusion energy ignition for the first time. This is the most important scientific breakthrough of the 21st century paves the way to future clean inertial fusion energy (IFE). The diamond (high density carbon (HDC)) ablator material used in this experiment displays detrimental effects due to the development of hydrodynamic instabilities at the diamond/fuel interface under shock compression. New alternatives to diamond ablators are required to step up the energy yield in ICF experiments. The unique combination of mechanical strength (approaching that of diamond), the ability to accommodate high-Z dopants (in contrast to diamond), and the tunability of the properties (through synthesis material with varying sp 3 content) make amorphous carbon (a-C) a promising material for next-generation IFE ablative capsules. However, despite its critical importance to the IFE program, the behavior of a-C carbon at extreme temperatures and pressures remains largely unexplored. The primary goals of this project were to perform groundbreaking dynamic compression experiments and predictive simulations to uncover the fundamental high-energy-density physics of amorphous carbon. Our goals were (1) to uncover the metastability range of amorphous carbon and probe phase transitions to diamond or metastable supercooled liquid carbon; (2) to acquire high-quality equation of state (EOS) data and develop an experimentally validated EOS from machine-learning MD simulations of the complex states of carbon; and (3) to uncover the complex behavior of carbon liquid in both thermodynamically stable and metastable supercooled states by accessing large areas of carbon phase diagram with amorphous samples with variable sp 3 content. Our proposed experimental program included measurements of equation of state and diffraction measurements using the Omega EP laser at the Laboratory of Laser Energetics at the University of Rochester. The theoretical/simulation program involved the development of machine-learning models of the complex response of amorphous carbon under dynamic compression by performing molecular dynamics simulations at experimental time and length scales using leadership class DOE supercomputers. Simulations guided experiments to observe predicted phenomena and acquire critical experimental data in specific pressure-temperature domains to validate theoretical models. This research delivered fundamental properties of novel amorphous carbon IFE ablator material, including phase diagram and EOS. These results will aid in IFE target design and implosion experiments. A unique combination of predictive simulations and dynamic and static experiments provided a highly inspirational intellectual environment for graduate students and postdocs involved in this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Electrodeposition of ligand-free copper nanoparticles from aqueous nanodroplets

Copper nanoparticles have emerged as promising electrocatalysts for energy storage and conversion. Generally, homogenous nanoparticle synthesis requires a stabilizing ligand, which may influence the electrocatalysis. Ligands can be avoided by direct nanoparticle electrodeposition. Here, we extend nanodroplet-mediated electrodeposition to the electrodeposition of copper nanoparticles from aqueous nanodroplets suspended in 0.1 M tetrabutylammonium hexafluorophosphate ([TBA][PF 6 ]) and 1,2-dichloroethane. Stochastic electrochemistry on microelectrodes was used to elucidate nanoparticle growth kinetics for various solution conditions and substrate materials. A study on the effect of surfactants on nanoparticle morphology and growth kinetics demonstrated the surfactant's role as an avenue for morphological control. Nanoparticle morphology was studied by Scanning Electron Microscopy (SEM) and Transmission Electron Microscopy (TEM). Energy-dispersive X-ray Spectroscopy (EDX) confirmed the presence of copper nanoparticles. TEM and Selected Area Electron Diffraction (SAED) confirmed nanoparticle size, morphology, and polycrystallinity. We demonstrate tunable nanoparticle size using cyclic voltammetry on single aqueous nanodroplets by altering the voltammetric sweeping potentials of copper electroreduction. Finite element simulations validate the voltammetric response and ability to control the nanoparticle size with nanometer resolution during a voltammetric sweep. These results inform the electrodeposition of copper nanoparticles from aqueous nanodroplets for important applications, such as the conversion of CO 2 to hydrocarbon fuels.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Field-Dependent Magnetic Domain Behavior in van der Waals Fe 3 GeTe 2

Two-dimensional magnetic van der Waals (vdW) materials can show a variety of topological nontrivial spin textures, such as Bloch- or Néel-type stripe, skyrmion, or bubble domains under certain external stimuli. It is critical to understand the magnetic domain behavior in vdW materials in order to control their size and density in response to external stimuli, such as electric and magnetic fields. In this study, we examine the magnetic field dependence of topologically non-trivial magnetization spin textures in vdW Fe 3 GeTe 2 . Néel-type stripe domains and skyrmions are formed depending on the magnetic field-cooling protocol used during in situ Lorentz transmission electron microscopy (LTEM) experiments. Use of quantitative reconstruction of magnetic induction maps and micromagnetic simulations allow for the understanding of the LTEM results of Néel-type stripe domains as well as skyrmions. In addition, the deformation of skyrmion contrast is observed as a result of the introduction of an in-plane magnetic field. We demonstrate the stability of the stripe domains and skyrmions in response to an externally applied magnetic field due to an energy barrier for domain wall annihilation. Our results establish an understanding of the energy landscape that governs the behavior of the topologically non-trivial spin textures in vdW materials which can be harnessed for spintronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Predicting non-linear stress–strain response of mesostructured cellular materials using supervised autoencoder

Recent breakthroughs in advanced manufacturing capabilities have made it possible to design and print sophisticated topologies of cellular structures using diverse engineering materials such as metals, polymers, and ceramics. In these architectured materials, it is often desirable to tailor the mechanical properties by altering the unit cell topology. This necessitates an in-depth understanding of how the topology of the unit cell structure affects the macroscopic behavior of the material in both the linear and the non-linear regimes encountered under large compression. Here, we have developed a machine learning (ML) approach capable of accelerating the prediction of the stress–strain response of a polymer-based cellular structure under uniaxial confined compression. As part of generating the training data for ML, 60,000 mesostructures were generated using a relatively novel approach based on cellular automata, and their corresponding stress–strain responses were obtained from the finite element simulations. Principal component analysis (PCA) was used to reduce the dimensionality of the stress–strain curves. With only 20 principal components, PCA captured 99.89% of the variance in the stress–strain curves while reducing the dimensionality by 5X. ML using supervised autoencoder was able to successfully speed up the prediction of the non-linear stress–strain response of a unit cell by up to 4600X. The proposed method can serve as an efficient data generation tool and a rapid means for predicting the structure–property relationship through accelerated forward modeling of cellular materials under compaction, in cases where the macroscopic stress–strain response is governed by the unit-cell topology.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Dynamic compression and spallation behavior of NbTaTiVZr high-entropy alloy

High-entropy alloys are a new class of materials with promising properties for aerospace and defense applications. The unique behavior of these materials is driven by the complex interactions between dissimilar atoms in a crystal and requires combined theoretical and experimental efforts to unlock their full potential. Here, we evaluate the relationship between the microstructure and the dynamic response of the equiatomic NbTaTiVZr alloy. Specifically, the shock Hugoniot, or equation of state, was measured up to a particle velocity of 0.5 mm μs −1 using gas gun plate impacts. Shock wave profile and incipient spallation experiments were used to characterize wave propagation and damage formation with post-mortem recovery experiments. Molecular dynamics simulations confirm the experimental findings and extend them up to 2.1 mm μs −1 particle velocity. Computational thermodynamic calculation of phase diagram simulations explain details of the material microstructure, which explains the measured strength and experimental damage patterns. Overall, this work provides detailed high-strain-rate characterization of a refractory high-entropy alloy, and more importantly, a framework and demonstration of the utility and necessity of a combined theoretical and experimental approach, outlining the importance of considering processing and manufacturing conditions when evaluating the performance of new materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Size-dependent errors in real-time electron density propagation

Real-time (RT) electron density propagation with time-dependent density functional theory (TDDFT) or Hartree–Fock (TDHF) is one of the most popular methods to model the charge transfer in molecules and materials. However, both RT-TDHF and RT-TDDFT within the adiabatic approximation are known to produce inaccurate evolution of the electron density away from the ground state in model systems, leading to large errors in charge transfer and erroneous shifting of peaks in absorption spectra. Given the poor performance of these methods with small model systems and the widespread use of the methods with larger molecular and material systems, here we bridge the gap in our understanding of these methods and examine the size-dependence of errors in RT density propagation. We analyze the performance of RT density propagation for systems of increasing size during the application of a continuous resonant field to induce Rabi-like oscillations, during charge-transfer dynamics, and for peak shifting in simulated absorption spectra. We find that the errors in the electron dynamics are indeed size dependent for these phenomena, with the largest system producing the results most aligned with those expected from linear response theory. The results suggest that although the RT-TDHF and RT-TDDFT methods may produce severe errors for model systems, the errors in charge transfer and resonantly driven electron dynamics may be much less significant for more realistic, large-scale molecules and materials.

Chemistry↗
DOE OSTI2020

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗
DOE OSTI2023

Development and testing of coarse-grained models for ultrasonic simulations of cast austenitic stainless steel

Ultrasonic inspection of cast austenitic stainless steel (CASS) in the nuclear industry is particularly challenging because of sound field scatter and attenuation caused by the coarse-grained microstructure. Modeling and simulation are important tools in ultrasonic testing, as they can be used to help address key aspects of inspections, such as developing new probe designs, predicting inspection reliability, and testing phased-array focal laws. However, developing a useful and reliable CASS model is challenging due to the many grain interfaces and crystalline orientations that must be captured. We demonstrate a method of creating a realistic CASS model that is usable in CIVA, a commercially available modeling and simulation software platform. Using polished and chemically etched sections, we generate models of a coarse-grained equiaxed specimen and a columnar specimen. We also test an alternative method of generating a coarse-grained model using Voronoi regions. We qualitatively compare sound field scatter and quantitatively compare sound field attenuation and beam partitioning in simulated sound fields to those of laboratory-measured sound fields. Results show that the Voronoi models perform as well as or better than the models based on actual grain morphology. Here we also show that model-to-model randomness in Voronoi grain structure can impact the magnitude of a simulated echo response by a factor of two or more. Although CASS models are potentially a good depiction of reality for a given scenario, they should not be considered representative since CASS morphology can change significantly from specimen to specimen or within the same specimen.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The impact of chromium coatings on Zircaloy cladding deformation behavior under reactivity-initiated accident-like mechanical loading conditions

A reactivity-initiated accident (RIA) occurs when a control rod ejection or control blade drop causes an increase in the fission rate. The injection of energy results in an increase in fuel temperature which in turn causes rapid thermal expansion of the fuel pellet. This thermal expansion may result in pellet-cladding mechanical interaction (PCMI) in which the fuel imparts a mechanical strain to the cladding. PCMI may cause the cladding to fail, and thus, the mechanical response of cladding due to PCMI must be investigated when characterizing new cladding materials. Chromium-coated Zircaloy-4 is a near-term accident-tolerant fuel cladding that exhibits improved high-temperature oxidation resistance. Modified burst testing was utilized to experimentally simulate the effects of PCMI on both uncoated and chromium-coated Zircaloy cladding samples at hot zero power conditions. Samples were coated using either cold spraying or physical vapor deposition to understand the differences in behavior that the coating application method may cause. Digital image correlation was used to analyze images of the deforming specimens to extract the in-situ strain behavior of the cladding. The uncoated specimens burst at hoop strains ranging from 8.8 % to 17.2 %. The cold-spray chromium-coated Zircaloy specimens burst at hoop strains of 7.0 % to 11.0 %. The physical vapor deposition coated tubes burst at hoop strains of 9.1 % to 11.5 %. These results indicate that the chromium coating causes a loss in the ductility of the cladding. Furthermore, the higher burst hoop strains of the physical vapor deposition-coated samples relative to the cold-spray samples indicate that the cold-spraying technique causes more of a loss in ductility than physical vapor deposition. All samples burst at higher hoop strains than those expected to occur in an RIA for fresh cladding.

Accident-tolerant fuel↗
DOE OSTI2020

Battery Abuse Testing

Abuse tests are designed to determine the safe operating limits of HEV\PHEV energy storage devices. Testing is intended to achieve certain worst-case scenarios to yield quantitative data on cell\module\pack response, allowing for failure mode determination and guiding developers toward improved materials and designs. Standard abuse tests with defined start and end conditions are performed on all devices to provide comparison between technologies. New tests and protocols are developed and evaluated to more closely simulate real-world failure conditions. While robust mechanical models for vehicles and vehicle components exist, there is a gap for mechanical modeling of EV batteries. The challenge with developing a mechanical model for a battery is the heterogeneous nature of the materials and components (polymers, metals, metal oxides, liquids). This year saw the stand up of a new drop tower tester capable of providing dynamic mechanical test results.

25 ENERGY STORAGE↗
DOE OSTI2021

I.6.C Battery Safety Testing

Abuse tests are designed to determine the safe operating limits of HEV\PHEV energy storage devices. Testing is intended to achieve certain worst-case scenarios to yield quantitative data on cell\module\pack response, allowing for failure mode determination and guiding developers toward improved materials and designs. Standard abuse tests with defined start and end conditions are performed on all devices to provide comparison between technologies. New tests and protocols are developed and evaluated to more closely simulate real world failure conditions. While robust mechanical models for vehicles and vehicle components exist, there is a gap for mechanical modeling of EV batteries. The challenge with developing a mechanical model for a battery is the heterogeneous nature of the materials and components (polymers, metals, metal oxides, liquids).

25 ENERGY STORAGE↗