Search NASA⌕ Search

SEARCH · Search NASA

Results for “material modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Dispatch Manager for NEML2 Constitutive Model Calculations Embedded in MOOSE

This report describes the extended capabilities of the NEML2 constitutive modeling library, including a flexible and efficient work dispatching system designed to leverage both CPU and GPU resources. This enhancement addresses one of the primary computational challenges in large-scale simulations: the ability to distribute and execute batches of material model evaluations across heterogeneous computing devices. The new dispatch system introduces a modular set of dispatcher and scheduler classes that coordinate the flow of data and execution between devices. The dispatcher is responsible for efficiently packaging work, managing device-specific memory operations, and synchronizing results. This modularity allows for extensibility, making it straightforward to integrate additional computing backends in the future. From an implementation standpoint, the dispatcher system interfaces seamlessly with NEML2's existing models. They handle device-aware tensor operations, optimize memory transfers, and support asynchronous execution when applicable. This design ensures that batches of material points can be evaluated concurrently, substantially improving throughput compared to previous single-device or serial implementations. These improvements not only enhance the raw performance of NEML2 but also improve its usability in multiscale and high-fidelity simulations, where the simultaneous evaluation of large material point batches is critical. Benchmarks included in the report demonstrate the system’s scalability, highlighting its effectiveness when leveraging modern GPU architectures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Using Hydrodynamic Similarity as a Verification Method for Impact Cratering Simulations in the FLAG Hydrocode

Hydrodynamic codes (hydrocodes) are common tools for modeling hypervelocity impacts to provide insight into the physical phenomenon. Hydrocodes can simulate impacts from micrometer to kilometer spatial scales and reach impact velocities difficult to achieve in experimental settings. However, numerical models are approximations, and demonstrating that a numerical method is capable of providing physical results for these models is essential. In this work, we employ a hydrocode verification technique that leverages hydrodynamic similarity, a mathematical property of the conservation equations of fluid mechanics that form the basis for hydrocode models. Using the FLAG hydrocode, we simulate aluminum (Al) and basalt projectiles and targets at spatial scales spanning 7 orders of magnitude (hundreds of micrometers to kilometers). These materials were chosen because Al-6061 is a common material in spacecraft and satellites and basalt is a useful approximation of rocky astronomical bodies. Our results show that hydrodynamic similarity holds for each material model used and across spatial scales. We show that under certain conditions hydrodynamic similarity can apply in the presence of gravity and that similarity does not hold in the presence of strength models. We conclude that the FLAG hydrocode preserves important mathematical properties of fluid dynamics in hypervelocity impacts of Al-6061 and basalt.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Report on DIC Cyclic Dryout Testing of ATF Cladding Materials

Boiling water reactors can experience anticipated operational occurrences, where an expected deviation from normal operating conditions is witnessed for short time scales. One example of such transients arises from cladding dryout, where stable boiling becomes unstable, resulting in cladding heat up before re-wetting occurs. To understand the impact of such transients on fuel cladding lifetime and performance, cyclic dryout conditions were imposed upon fresh fuel cladding in the absence of irradiation. Pressurized cladding tube segments were exposed to rapid temperature increases to intermediate temperatures below a cladding rupture threshold and deformation was monitored in situ via 3-D digital image correlation. The impact of cycle frequency on cladding deformation was assessed. A clear understanding of elastic and plastic strain responses were developed for Zr alloy deformation under temperature transients, and these results can be used to support material model development and validate fuel performance codes.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Anisotropic Effects of Loading on Zircaloy Cladding Performance

This presentation provides recent results related to the mechanical behavior of Zircaloy cladding for the Light Water Reactor Fuels Workshop/CRAFT Initiative Update presented by the Advanced Fuels Campaign. This presentation focuses on Zircaloy cladding response to pellet-cladding interactions and how mechanical anisotropy makes the cladding sensitive to the testing mode. Results are compared to existing material models used by BISON, and implications for reactor performance are briefly discussed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Generalizable machine learning potentials for quantum-accurate predictions of non-equilibrium behavior in 2D materials

Machine learning interatomic potentials (ML-IAPs) are emerging as transformative tools in materials modeling, promising quantum-level accuracy at a fraction of the computational cost. However, their ability to generalize beyond equilibrium configurations and to reliably capture defect- and temperature-driven behavior remains underexplored. Here, we develop and benchmark two state-of-the-art ML-IAPs, Spectral Neighbor Analysis Potential (SNAP) and Allegro, on a comprehensive dataset for monolayer MoSe₂. Using density functional theory (DFT) as the reference, we evaluate their performance in capturing stress–strain behavior, phase transition energetics, defect evolution, edge stability, and fracture toughness. Allegro, a deep equivariant neural network potential, surpasses both SNAP and the classical Tersoff potential in accuracy, efficiency, and transferability. Importantly, both ML potentials accurately reproduce experimental fracture measurements and ab initio predictions of inversion domain formation—phenomena well beyond their training sets. Our findings establish ML-IAPs as viable replacements for traditional force fields in the study of non-equilibrium mechanical phenomena, enabling large-scale, high-fidelity simulations in 2D materials and beyond. In conclusion, this work provides a broadly applicable framework for data-driven modeling of structural and functional transformations under extreme conditions.

2D materials↗
DOE OSTI2025

Failure Strain and Related Triaxiality of Aluminum 6061-T6, A36 Carbon Steel, 304 Stainless Steel, and Nitronic 60 Metals, Part I: Experimental Investigation

The objective of this study is to develop failure-limit material models for Aluminum 6061-T6, A36 Carbon Steel, 304 Stainless Steel, and Nitronic 60 metals, based on parameters of plastic equivalent strain (failure strain) and stress triaxiality. The research is conducted in two parts. This paper presents Part One of the study. In Part One, custom-designed test specimens undergo controlled uniaxial tension and compression testing at ambient temperature. These tests are performed at quasi-static speeds using Universal Testing Machines (UTMs) in accordance with ASTM E8 and ASTM E9 standards. Experimental data, specifically engineering stress–strain and force–displacement curves, are recorded from the onset of loading until specimen fracture, or in the case of compression tests, until the capacity of the testing machine is reached. In Part Two, the emphasis shifts to the calibration of Finite Element Analysis (FEA) models of the custom-designed test specimens. Plastic equivalent strain and the corresponding stress triaxiality values at failure are extracted from each test specimen for the given metal. These values are then systematically plotted onto a single graph to construct the failure-limit curve, which delineates the boundary conditions for material failure. This approach will facilitate the development of a comprehensive material property definition that correlates plastic equivalent strain with stress triaxiality at failure for Aluminum 6061-T6, A36 Carbon Steel, 304 Stainless Steel, and Nitronic 60 metals.

Harwell, Ron↗
DOE OSTI2026

Fast and accurate calculation of EXAFS Debye-Waller factors in U⁢O2 using the dynamical matrix method

Theoretical modeling of bonding dynamics in metal oxides is required for predicting their thermal conductivity, catalytic activity, and mechanical properties. A primary challenge is the scarcity of experimental methods for validating theoretical predictions of these atomic-scale dynamics. This work presents a workflow that uses experimental extended x-ray absorption fine structure (EXAFS) data collected at high temperatures to validate an interatomic force field for uranium dioxide (UO2), an important model material. The validated force field is then used to drive computationally intensive molecular dynamics (MD) simulations and as input for the much faster dynamical matrix Debye-Waller (DMDW) method. The predicted values of the Debye-Waller factors from the DMDW calculations are in good agreement with those obtained from the MD simulations, with residual pair-specific differences attributable to quantum zero-point motion at low temperatures and lattice anharmonicity at high temperatures. We further show that theoretical EXAFS spectra constructed directly from DMDW-derived Debye-Waller factors reproduce the experimental data (at relatively low temperatures) with accuracy comparable to full MD-EXAFS, providing an additional validation of the choice of the potential. This study establishes a validated, rapid computational pathway for modeling bond dynamics, naturally incorporating quantum nuclear\\\\r\\\\nstatistics absent in classical simulations, which are essential for the mechanistic understanding of complex oxide materials.

58 GEOSCIENCES↗
DOE OSTI2025

Automated Redox Titrations via Interdigitated Electrode Arrays: Application to the Mediated Electron Transfer Interrogation of Charge and Rate on Electrodeposited Polymers

Mediated electron transfer (MET) plays a crucial role in energy storage and conversion technologies such as redox targeting flow batteries (RTFBs), yet its experimental investigation often requires labor-intensive and low-throughput setups. To address this, we developed a microfabricated interdigitated electrode array (IDA) platform that enables automated, high-throughput electrochemical redox titration measurement to be performed to study the MET process. Our redox titration method enables simultaneous measurement of the charge capacity and rate of MET processes on a material or surface. Automated redox titration (ART) facilitates systematic investigation of the MET process across a broad parameter space, exemplified through the study of polypyrrole (PPy) and a pyrene-4,5,9,10-tetrone azo group-based polymer (PTAP), both redox-active polymers relevant to various energy storage applications. Using PPy as a model material, 500 redox titration measurements were conducted within 50 h, varying the electrode gap widths, polymer charging potentials, voltammetric scan rates, and electrolyte concentrations. Finite-element simulations confirmed the electrochemical responses and elucidated the kinetics of the MET reactions. Our automated methodology was further tested with PTAP, revealing a surprising charging potential dependence on the rate of MET. The automation, flexibility, and scalability of our redox titration platform pave the way not only for advanced studies of MET processes relevant to RTFBs, but also with implications in the understanding of next-generation energy storage materials, molecular electrocatalysis, and biosensing.

electrochemical analysis↗
DOE OSTI2026

Finite deformation implementation of a mixed-mode single-integral type cohesive zone with reorienting surfaces of separation

To model material ductile failure and crack propagation, cohesive zone elements can be embedded along potential fracture paths in a finite element simulation. When damage criteria are met, elements in the mesh decohere, simulating the formation and propagation of a crack. In this paper, we present a novel computational algorithm based on finite deformation theory, essential to modeling crack initiation and growth in solids undergoing large deformations. This new algorithm was formulated within a Lagrangian frame of reference to extend previous cohesive zone algorithms to include modeling crack growth in finite deformation contexts. The local coordinate system, necessary for defining an embedded cohesive zone, is constructed based upon the current configuration and is updated within the nonlinear iteration process, thereby resulting in the convergence of the solution for a growing crack in a large deformation quasi-static setting. The model’s accuracy was demonstrated by comparing finite element model simulation results with the analytic case of a constant surface separation, as shown in the verification examples. The power and efficacy of the algorithm to capture large deformations during crack growth were then demonstrated with a double cantilever beam example case. It indicates that the model can be applied to a variety of physical circumstances for predicting crack initiation and growth with delamination and fracture.

42 ENGINEERING↗
DOE OSTI2025

Electron Microscopy Transfer System to Protect Atmosphere‐Sensitive Materials for Scanning Electron Microscopy Characterization

Atmosphere- and/or moisture-sensitive materials can be challenging to characterize using electron microscopy techniques due to sample preparation workflows that generally require exposure to ambient conditions. Here, we describe a novel preparation method that uses aluminum foil in combination with a commercial cryo-EM transfer system to circumvent undesired exposure to the atmosphere. First, hygroscopic MgCl 2 was used as a model material, and prepared samples (both protected and unprotected) were placed in a controlled-humidity environment (> 80% relative humidity) for various exposure lengths (circa seconds to hours). Following this, the effectiveness of the sample preparation method was determined by comparing qualitative photos and quantitative X-ray diffraction patterns between the two sample subsets. The combined results of these experiments suggest that the outlined preparation method effectively protects MgCl 2 from atmospheric contamination compared to MgCl 2 samples that had no protective measures taken. Finally, the preparation method was utilized to protect a highly hygroscopic crystalline BaO thin film for characterization via scanning electron microscopy, thereby demonstrating a functional application of the outlined preparation technique and an additional use for the commercial cryo-EM transfer system beyond its intended application.

atmosphere-sensitive materials↗
DOE OSTI2025

Exploring a New Regime of Molecular Orbital Physics in 4⁢d Cluster Magnets with Resonant Inelastic X-Ray Scattering

Molecular orbital systems with clusters of heavy transition metal (TM) ions are one of the most important classes of model materials for studying the interplay between local physics and effects of itinerancy. Despite a large number of candidates identified in the family of 4⁢d TM materials, an understanding of their physics from competing microscopic energy scales is still missing. Here, we bridge this gap by reporting the first resonant inelastic x-ray scattering (RIXS) measurement on a well-known series of Ru cluster magnets with a 6⁢H-perovskite structure Ba 3 ⁢M⁢Ru 2 ⁢O 9 (M 3+ = In 3+ , Y 3+ , La 3+ ) comprised of Ru dimers. In addition to providing a microscopic explanation for their anomalous magnetic properties, our RIXS measurements combined with theoretical modeling uncover a new regime of molecular orbital physics where the combined effect of large hopping and small spin-orbit coupling results in highly fragile electronic states in the Ru-dimer compounds directly manifested as an abrupt change in the RIXS spectrum accompanying a tiny change in the local structure tuned by the M-site ion. This unique combination of energy scales found only in the 4⁢d but not the 5⁢d cluster magnets highlights the value of these materials as a new platform for studying quantum phase transition involving molecular orbitals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

WHA Sphere Impacting RHA Plate: Modeling Brittle Behavior in Metal Alloys

Brittle behavior of metal alloys is often critical to modeling ballistic impact and penetration. The ALEGRA multiphysics finite element software incorporates calibrated models for the equation of state, elasticity, yield stress, plasticity and fracture, but simulations do not always capture expected metal fracture. Here we report concerted efforts to do so for one important case where experiments clearly show shear fractures: a tungsten sphere impacting a steel plate at various angles. Our best simulations show fractures that are qualitatively similar to experiments, but there are significant differences in quantitative metrics. Specifically, velocities of tracers used to quantify simulated plug parameters consistently fall short of measured plug velocities. Also, simulated plugs break apart more than expected from experimental evidence. We attribute these shortfalls to the lack of an explicit shear fracture mechanism in the material models, leading to over-estimated resistance to plug formation and movement.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Synthetic Band Structure Engineering of Graphene Using Block Copolymer-Templated Dielectric Superlattices

Engineering the electronic band structure of two-dimensional (2D) materials by imposing spatially periodic superlattice (SL) potentials opens a pathway to unconventional electronics. Nanopatterning the gate electrode or surface dielectric near 2D crystals provides a powerful strategy for realizing electrostatically tunable “remote” SLs with flexibility in lattice design. Here, we demonstrate the effectiveness of block copolymer (BCP)-templated dielectric nanopatterns for fabricating etch-free high-grade metal oxide SLs. Alumina (AlO x ) nanopatterns with hexagonal symmetry and a 38 nm SL wavelength are produced as a model material by directly converting a self-assembled BCP film via block-selective vapor phase infiltration. Despite micrometer-scale rotational disorder inherent to BCP self-assembly, electronic transport measurements of graphene reveal replica Dirac points at zero field and Hofstadter mini-gaps under finite magnetic fields. These results indicate the successful formation of remote SL potentials in graphene resulting from optimized AlO x nanopattern fabrication to achieve consistent lattice symmetry and periodicity at a macroscopic scale. The findings of this study, combined with the versatile, scalable, and cost-effective nature of BCP nanopatterning, highlight the potential of BCP-templated nanostructures for remote SL engineering in 2D crystals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Autocorrelation and Multifractal Detrended Fluctuation Analyses Reveal Superdiffusive Mass Transport in Solvent-Filled Nanoporous Media

Fluorescence fluctuation spectroscopy experiments were conducted to better understand the complex mass transport dynamics of organic molecules in liquid-filled nanoporous media. Anodic aluminum oxide (AAO) membranes incorporating 10 and 20 nm diameter cylindrical pores were employed as model materials. Nile red (NR) dye was used as a fluorescent tracer. The dye was dissolved separately in ethanol and toluene at a concentration of 20 nM and used to fill the membrane nanopores. Confocal fluorescence microscopy was employed to capture photon intensity time series data reflecting apparent diffusion of the dye within the pores. Autocorrelation of these data revealed that NR diffusion within the membranes occurred over a broad range of time scales. The autocorrelation decays were fit to a model for one-dimensional diffusion incorporating both fast and slow components having apparent diffusion coefficients, D f and D s , differing by a factor of ∼100. The fast mechanism was attributed to hindered bulk-like diffusion in the central pore cavity, while slow diffusion likely involved absorption of the dye to the pore surfaces. Unfortunately, important evidence of diffusion anomalies is lost in the broad autocorrelation decays obtained. The method of multifractal detrended fluctuation analysis (MF-DFA) was applied to the same data as a means to overcome this limitation. MF-DFA revealed that time series acquired from within the nanopores were multifractal and exhibited evidence of anomalous superdiffusion, likely resulting from the participation of a desorption-mediated diffusion mechanism. Monte Carlo simulations of time series modeling desorption-mediated diffusion in cylindrical nanopores provided support for this assignment. Here, the new knowledge gained affords an improved understanding of hydrocarbon dynamics within nanoporous oil and gas shales.

Diffusion↗
DOE OSTI2026

Modeling thermal radiation waves in silica plasmas for the Mooncat NIF experiment

The Mooncat experiment on the National Ignition Facility uses a laser-driven hohlraum to create a thermal radiation wave in a titanium-doped silica plasma. The titanium dopant enables absorption spectroscopy measurements to infer the temperature of the wave as it propagates. This measurement can be used to constrain multi-physics simulation codes to better understand when simulations do not match an experiment. In this paper, we present radiation-hydrodynamics simulations of the thermal radiation wave in the first full-platform shots of the Mooncat experiment. We examine the important parameters of the simulation, focusing on the radiation temperature source, the material model of the silica plasma as it pertains to radiation transport, and lateral leakage through a beryllium tube enclosing the silica. We compare different simulation modeling strategies to an analytic model of diffusive radiation transport and find that the simulation agrees with the analytic model when it is sufficiently simplified. These simulations show how radiation energy couples to matter to develop a shock wave in a radiative heat wave, an important topic in astrophysics and nuclear fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Structural and Mechanical Analysis of Individual Mineralized Collagen Fibrils Using In Situ Transmission Electron Microscopy

Bone serves as an example of nature’s architectured material with its characteristic blend of strength and toughness, all at a lightweight design. Given the hierarchical nature of these materials, it is essential to understand the governing mechanisms and organization of their constituents across length scales for bioinspired structural design. Despite recent advances in transmission electron microscopy (TEM) that have allowed us to witness the hierarchical arrangement of bone at micro-down to the nanoscale, we are still missing the details about the structural organization and mechanical properties of the main building blocks of bone─mineralized collagen fibrils (MCFs). Here, we present a method to extract individual MCFs from nature’s model material, mineralized turkey leg tendon, using a dropcasting procedure. By isolating the MCFs onto TEM supporting grids, we visualized the arrangement of organic and mineral phases within individual MCFs at the nanoscale. Using a four-dimensional scanning transmission electron microscopy (4D-STEM) approach, the orientation of individual mineral crystals within the MCFs was examined. Furthermore, we conducted in situ tensile experiments, revealing exceptional tensile strains of at least 8%, demonstrating the intricate relationship between structural organization and the mechanical behavior of MCFs. These insights into the ultrastructure of mineralized tissue building blocks, as well as the proposed sample-extraction method compatible with in situ mechanical testing, provide a strong basis for research into nature-inspired material design.

4D-STEM↗
DOE OSTI2025

A hereditary integral, transient network approach to modeling permanent set and viscoelastic response in polymers

An efficient numerical framework is presented for modeling viscoelasticity and permanent set of polymers. It is based on the hereditary integral form of transient network theory, in which polymer chains belong to distinct networks each with different natural equilibrium states. Chains continually detach from previously formed networks and reattach to new networks in a state of zero stress. The free energy of these networks is given in terms of the deformation gradient relative to the configuration at which the network was born . A decomposition of the kernel for various free energies allows for a recurrence relationship to be established, bypassing the need to integrate over all time history. The technique is established for both highly compressible and nearly incompressible materials through the use of neo-Hookean, Blatz–Ko, Yeoh, and Ogden-Hill material models. Multiple examples are presented showing the ability to handle rate-dependent response and residual strains under complex loading histories.

Finite element method↗
DOE OSTI2025

A Standardized Analysis Process Using Digital Image Correlation to Calculate In Situ Cladding Strain from Modified Burst Tests for Fuel Performance Code Validation

Historical data collection on nuclear fuel cladding materials has focused on generating a statistically significant amount of data to assess the material and its failure behavior. Furthermore, data generated to support material model and failure criteria development were previously posttest evaluations, so a large number of tests was required to gain new understanding. A way to expedite this process is to develop techniques capable of generating large, high-fidelity data sets from a single test with lower uncertainty or quantified uncertainty. One such example of this approach is Oak Ridge National Laboratory’s use of modified burst tests (MBTs) to analyze the mechanical behavior and failure conditions of cladding during a simulated reactivity-initiated accident (RIA). Each test incorporates digital image correlation (DIC) analysis techniques that are used to assess the accumulated strain in situ, as well as eventual cladding failure. This work has been fruitful in defining strain-to-failure conditions for materials like silicon carbide (SiC) fiber–reinforced/SiC matrix composite tubes (SiC/SiC), iron-chromium-aluminum (FeCrAl) alloy tubes, and chromium-coated Zircaloy-4 tubes. However, there are numerous DIC software available, including open-source and proprietary software. The different DIC software use various algorithms to process images and calculate displacement values. Using these different software and algorithms can lead to varying results, and perhaps larger-than-expected uncertainties. In the present study, previously published MBT data encompassing a variety of test conditions were reanalyzed with two different DIC software to assess the variance in the calculated strain results. The data consisted of SiC/SiC, FeCrAl, and chromium-coated Zircaloy-4 tubes. Plots of the calculated strains during the transient revealed good agreement between the two DIC software. The average root-mean-square errors between the two software was 0.20% strain, which is slightly larger than a previously reported error value for these tests. In conclusion, this variance in results is low enough that this analysis method can be used for code validation.

Reactivity-initiated accident↗