Search NASA⌕ Search

SEARCH · Search NASA

Results for “Material Response 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 55 records · Page 3

Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal

Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.

36 MATERIALS SCIENCE↗

X-ray-induced piezoresponse during X-ray photon correlation spectroscopy of PbMg 1/3 Nb 2/3 O 3

X-ray photon correlation spectroscopy (XPCS) holds strong promise for observing atomic-scale dynamics in materials, both at equilibrium and during non-equilibrium transitions. Here an in situ XPCS study of the relaxor ferroelectric PbMg 1/3 Nb 2/3 O 3 (PMN) is reported. A weak applied AC electric field generates strong response in the speckle of the diffuse scattering from the polar nanodomains, which is captured using the two-time correlation function. Correlated motions of the Bragg peak are also observed, which indicate dynamic tilting of the illuminated volume. This tilting quantitatively accounts for the observed two-time speckle correlations. The magnitude of the tilting would not be expected solely from the modest applied field, since PMN is an electrostrictive material with no linear strain response to the field. A model is developed based on non-uniform static charging of the illuminated surface spot by the incident micrometre-scale X-ray beam and the electrostrictive material response to the combination of static and dynamic fields. The model qualitatively explains the direction and magnitude of the observed tilting, and predicts that X-ray-induced piezoresponse could be an important factor in correctly interpreting results from XPCS and nanodiffraction studies of other insulating materials under applied AC field or varying X-ray illumination.

36 MATERIALS SCIENCE↗

Modification of the MTS model for high strain-rate behavior of TI-6AL-4V

The Mechanical Threshold Stress (MTS) model provides excellent predictive capabilities for the material constitutive response for a wide range of temperatures and strain rates. However, the MTS model fails to capture the rapidly increasing yield stress at high strain rate behavior as the deformation controlling mechanism transitions from thermal activation to drag mechanisms, only capturing the linear behavior. Further, the model typically over predicts the flow stress behavior at yield and post yield due to its use of a constant work hardening rate parameter derived from the stress–strain response at constant saturation stress. An alternative approach to fitting portions of the MTS model is investigated and mathematical models are developed to address these issues. The results show that with appropriate experimental data, the mechanical threshold stress and work hardening rate parameters within the MTS model can quite easily and accurately be modified to extend applicability to high strain rate behavior and more accurately model the initial flow stress behavior at early work hardening rates without modification of the functions core to the MTS model itself.

36 MATERIALS SCIENCE↗

Identify the influence of microstructure on mesoscale creep and fatigue damage

This report describes the development of a microstructural model that can quantify the uncertainty in the observed rupture life of Grade 91 steel. The model is microstructural, meaning it relates microstructural characteristics of the material to the resulting material response. As such, one of the uses of this model is to identify the key microstructural parameters controlling the development of damage in Grade 91 operating at elevated temperatures. The report describes two veins of work: improvements to the crystal plasticity model required to run the uncertainty quantification analysis and the results of that UQ analysis. For creep, the model identifies the grain boundary diffusivity as the critical parameter controlling the rupture life of the material. The report demonstrates that a reasonable microstructural distribution of grain boundary diffusivity can account for the observed macroscale variation in rupture life at fixed temperature and load.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Direct Imaging of Hydrogen‐Driven Dislocation and Strain Field Evolution in a Stainless Steel Grain

Hydrogen embrittlement (HE) poses a significant challenge to the durability of materials used in hydrogen production and utilization. Disentangling the competing nanoscale mechanisms driving HE often relies on simulations and electron-transparent sample techniques, limiting experimental insights into hydrogen-induced dislocation behavior in bulk materials. This study employs in situ Bragg coherent X-ray diffraction imaging to track three-dimensional (3D) dislocation and strain field evolution during hydrogen charging in a bulk grain of austenitic 316 stainless steel. Tracking a single dislocation reveals hydrogen-enhanced mobility and relaxation, consistent with dislocation dynamics simulations. Subsequent observations reveal dislocation unpinning and climb processes, likely driven by osmotic forces. Additionally, nanoscale strain analysis around the dislocation core directly measures hydrogen-induced elastic shielding. These findings experimentally validate theoretical predictions and offer mechanistic insights into hydrogen-driven dislocation behavior. The quantified nanoscale phenomena serve as critical inputs for multiscale modeling frameworks to predict bulk material responses and accelerate the development of HE-resistant alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Distilling physical origins of hardness in multi-principal element alloys directly from ensemble neural network models

Abstract Despite a plethora of data being generated on the mechanical behavior of multi-principal element alloys, a systematic assessment remains inaccessible via Edisonian approaches. We approach this challenge by considering the specific case of alloy hardness, and present a machine-learning framework that captures the essential physical features contributing to hardness and allows high-throughput exploration of multi-dimensional compositional space. The model, tested on diverse datasets, was used to explore and successfully predict hardness in Al x Ti y (CrFeNi) 1- x - y , Hf x Co y (CrFeNi) 1- x - y and Al x (TiZrHf) 1- x systems supported by data from density-functional theory predicted phase stability and ordering behavior. The experimental validation of hardness was done on TiZrHfAl x . The selected systems pose diverse challenges due to the presence of ordering and clustering pairs, as well as vacancy-stabilized novel structures. We also present a detailed model analysis that integrates local partial-dependencies with a compositional-stimulus and model-response study to derive material-specific insights from the decision-making process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Real-time characterization of dislocation slip and twinning of shock-compressed molybdenum single crystals

Characterizing the fundamental micromechanisms activated during plastic deformation is critical to explain the macroscopic shock response of materials and develop accurate material models. In this paper, we investigate the orientation dependence, and the mediated slip and twin systems on [100] and [111] bcc molybdenum single crystals shock compressed up to 18 GPa with real-time Laue x-ray diffraction measurements. We report that dislocation slip along the {110} < 111 > and {112} < 111 > systems is the governing deformation mechanism during compression with negligible anisotropy observed at the Hugoniot state. Here, we provide real-time evidence that molybdenum undergoes deformation twinning along {11¯2} < 111 > during shock release.

36 MATERIALS SCIENCE↗

Material Model Parameters Optimization in Liquid Mercury Target Dynamics Simulation With Machine Learning Surrogates

A pulsed spallation target is subjected to very short (∼0.7μs) but intense loads (23.3 kJ) from repeated proton pulses, which knock away neutrons from the mercury atoms’ nuclei for a wide range application in physics, engineering, medicine, petroleum exploration, biology, chemistry, etc. The effect of this pulsed loading on the stainless-steel target module which contains the flowing mercury target material is difficult to predict not only due to its short but intense explosive-like physical reaction, but also the nonlinear material behavior of the liquid mercury in the structure. Injecting small helium bubbles in the mercury has been an efficient method of mitigating the pressure wave at high power level stage. However, prediction of the resultant loading on the target is more difficult when helium gas is intentionally injected into the mercury. A 2-phase material model that incorporates the Rayleigh-Plesset (R-P) model is expected to address this complex multi-physics dynamics problem by including the bubble dynamics in the liquid mercury. A parameter sensitivity study was firstly employed to understand their impact on the simulation strains. The investigated parameters included E, μ, γ, σ, n, VFgas, and gas cumulative volume curve control parameters a and b. Verification and validation results from sparse polynomial expansions (SPE) method and directional Gaussian smoothing (DGS) optimization show that the surrogate model had training error of ∼7% and validation error of ∼15%, indicating that machine learning methods and surrogate models can help optimize the uncertain parameters in the complex 2-phase material model. This approach is expected to fill the knowledge gap between unknown liquid-gas mixture material model and measured vessel strain responses.

Lin, Lianshan↗

U-net architected deep material network training with microstructure local field information

The Deep Material Network (DMN) has recently emerged as a powerful reduced-order modeling framework for simulating the mechanical response of heterogeneous materials such as composites. Unlike most data-driven approaches that directly learn a material’s response under prescribed loading, the DMN acts as a homogenization operator, learning the kinematic constraints and mechanical interactions of the underlying microstructure. However, traditional DMN training relies exclusively on homogenized effective properties derived from Direct Numerical Simulations (DNS), discarding the rich local field data that govern microstructural interactions. In this work, we extend the DMN framework to incorporate such local field information into the offline training process. Utilizing a U-Net architecture, we augment the DMN training objective to include the first and second statistical moments of the local stress fields obtained from linear DNS. This ensures that the learned network topology not only fits the effective stiffness but also accurately reflects the internal local stress and strain partitioning of the microstructure. The results confirm that supervising the localization process during training yields a superior surrogate model, reducing local prediction errors by an order of magnitude and significantly improving generalization to unseen nonlinear constitutive behaviors compared to traditional DMNs.

36 MATERIALS SCIENCE↗

Interfacial thermal transport in spin caloritronic material systems

Interfaces often govern the thermal performance of nanoscale devices and nanostructured materials. As a result, accurate knowledge of thermal interface conductance is necessary to model the temperature response of nanoscale devices or nanostructured materials to heating. Here, we report the thermal boundary conductance between metals and insulators that are commonly used in spin-caloritronic experiments. We use time-domain thermoreflectance to measure the interface conductance between metals such as Au, Pt, Ta, Cu, and Al with garnet and oxide substrates, e.g., NiO, yttrium iron garnet (YIG), thulium iron garnet (TmIG), Cr 2 O 3 , and sapphire. We find that, at room temperature, the interface conductance in these types of material systems range from 50 to 300 MW m –2 K –1 . Furthermore, we also measure the interface conductance between Pt and YIG at temperatures between 80 and 350 K. At room temperature, the interface conductance of Pt/YIG is 170 MW m –2 K –1 and the Kapitza length is ~40 nm. A Kapitza length of 40 nm means that, in the presence of a steady-state heat current, the temperature drop at the Pt/YIG interface is equal to the temperature drop across a 40-nm-thick layer of YIG. At 80 K, the interface conductance of Pt/YIG is 60 MW m –2 K –1 , corresponding to a Kapitza length of ~300 nm.

36 MATERIALS SCIENCE↗

Viscoelastic Material Characterization and Modeling of Photovoltaic Module Packaging Materials for Direct Finite-Element Method Input

Numerical tools, such as the finite-element method, are increasingly used to design and evaluate the photovoltaic (PV) modules, providing for the reduction of development time and improved performance and reliability. However, high-fidelity material models are necessary to accurately model the complex structural behavior of the involved packaging materials. A common simplification used in recent years is to model the polymer materials (i.e., encapsulant and backsheet) as linear elastic, which will lead to inaccurate results. Therefore, in this work, we present a thorough characterization of the time- and temperature-dependent mechanical response of predominant PV module encapsulant and backsheet materials. Based on this material characterization, we developed and experimentally validated generalized Maxwell models to describe each material's viscoelastic response. In addition, we included measurements of the coefficient of thermal expansion and presented all material models in such a fashion for direct input into commercial finite-element method modeling software.

FEM↗

Experimental verification of a crystal plasticity-based simulation framework for predicting microstructure and geometric shape changes: Application to bending and Taylor impact testing of Zr

This paper is concerned with experimental verification of a recently developed multi-scale simulation framework for plastic deformation of metallic materials from quasi-static to impact deformation conditions. The framework is a visco-plastic self-consistent (VPSC) polycrystalline model embedded in an implicit finite element method (FE-VPSC) to provide a microstructure-sensitive constitutive response at each material point. Each material point of the FEM model is a polycrystalline aggregate with crystallographic deformation mechanisms operating at the single crystal scale with their evolving activity based on a dislocation density-based hardening law and texture. Four beams and three cylinders machined in different orientations from a textured plate of high-purity zirconium are tested quasi-statically in 4-point bending and at speeds of 100 m/s, 170 m/s and 243 m/s during Taylor impact tests, respectively. The variation in dimensional changes resulting from different sample orientations in the plate with respect to loading directions is measured for each sample. Moreover, texture and twinning characterization is performed using electron backscattered diffraction (EBSD). The deformation processes and underlying evolution of microstructure are successfully simulated using the FE-VPSC framework. In doing so, the model parameters are optimized and validated across a broad range of strain rates and temperatures. Simulation results in terms of geometrical changes and microstructural evolution are compared with the experimental measurements. The model predicts anisotropic material flow resulting from the hard-to-deform crystallographic directions, the development of gradients in texture and twinning through the geometries, tension–compression asymmetry, as well as the extent of plasticity under impact.

42 ENGINEERING↗

A Machine Learning Framework for Modeling Ensemble Properties of Atomically Disordered Materials

Atomic disorder can strongly influence material properties such as charge transport, optical response, and catalytic activity. However, efficiently modeling these disorder effects remains challenging for first-principles methods due to the cost of sampling large configurational spaces and computing complex physical quantities. Recent advances of machine learning techniques, particularly graph neural networks (GNNs), has enabled the efficient and accurate predictions of complex material properties, offering promising tools for studying disordered systems. In this work, we present a general machine-learning-assisted computational framework that integrates equivariant GNNs with Monte Carlo simulations to compute the thermodynamic and ensemble-averaged functional properties of disordered materials. Using the surface-termination-disordered MXene monolayer Ti 3 C 2 T 2–x as a representative system, we find that electrical conductivity exhibits an emergent peak near the order–disorder phase transition temperature due to the interplay between electron scattering and doping. In contrast, optical conductivity remains largely insensitive to local atomic disorder and reflects the global surface chemical composition. These results highlight the role of atomic disorder in affecting material properties and demonstrate the potential of our approach for statistically modeling disorder effects in a wide range of materials such as high-entropy alloys and spin liquids.

MXene↗

Direct imaging of shock wave splitting in diamond at Mbar pressure

Understanding the behavior of matter at extreme pressures of the order of a megabar (Mbar) is essential to gain insight into various physical phenomena at macroscales—the formation of planets, young stars, and the cores of super-Earths, and at microscales—damage to ceramic materials and high-pressure plastic transformation and phase transitions in solids. Under dynamic compression of solids up to Mbar pressures, even a solid with high strength exhibits plastic properties, causing the induced shock wave to split in two: an elastic precursor and a plastic shock wave. This phenomenon is described by theoretical models based on indirect measurements of material response. The advent of x-ray free-electron lasers (XFELs) has made it possible to use their ultrashort pulses for direct observations of the propagation of shock waves in solid materials by the method of phase-contrast radiography. However, there is still a lack of comprehensive data for verification of theoretical models of different solids. Here, we present the results of an experiment in which the evolution of the coupled elastic–plastic wave structure in diamond was directly observed and studied with submicrometer spatial resolution, using the unique capabilities of the x-ray free-electron laser (XFEL). The direct measurements allowed, for the first time, the fitting and validation of the 2D failure model for diamond in the range of several Mbar. Our experimental approach opens new possibilities for the direct verification and construction of equations of state of matter in the ultra-high-stress range, which are relevant to solving a variety of problems in high-energy-density physics.

36 MATERIALS SCIENCE↗

Integrating Advanced Modeling and Accelerated Testing for a Modernized Fuel Qualification Paradigm

With the increasing interest in sodium fast reactor technology, as seen by applications to the U.S. Nuclear Regulatory Commission for the OKLO Aurora plant, fuel testing for the TerraPower Traveling Wave Reactor, and the impending construction and startup of the versatile test reactor (VTR), a modernized, accelerated approach to fuel qualification is needed. To guide this effort, a Phenomena Identification Ranking Table–styled analysis was performed for a U-Pu-Zr sodium-free annular fuel system. This analysis evaluated a series of fuel design properties and parameters against their contributions to key fuel performance phenomena. The resulting priority parameters were then reviewed against existing modeling and experimental capabilities to support investigation of the highest-priority parameters. A pathway for qualification was then established using highthroughput, high-volume experiments from MiniFuel and FAST in parallel with advanced physics-based model development. This effort outlines how the first stages of qualification can be reduced from the typical 20+-year development cycle to 5 to 7 years by deploying accelerated irradiation testing platforms. As with any accelerated test, these methods are prototypic in some aspects and less so in others; however, by coupling with advanced fuel performance modeling and simulation capabilities, the larger space of irradiation parameters and material response provided offers advantages for the validation of physics-based models supporting the deployment of novel fuel designs. As a test case, this paper utilizes a proposed Mark II fuel system for the upcoming VTR. Thus, an accelerated qualification method can be tested for the development of MARK II driver fuel so that by the time of VTR startup, lead test assemblies for a Mark II fuel can be initiated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Development and validation of a gamma-ray detection response model for neutron-diagnosed subcritical experiments

In a neutron-diagnosed subcritical experiment, a highly sensitive and fast detector is required to accurately determine the true fission γ-ray emission rate from an object composed of special nuclear material under interrogation, in addition to a proper detector response model for use in forward modeling. We present the development of a detector response model that is based in GEANT4 and uses post-processing techniques to include photomultiplier tube time and pulse height response effects. We also define a set of detector performance evaluation metrics that emphasized our experimental objectives. A suite of time and pulse height response measurements for several diverse detector assemblies are described. These assemblies were simulated with this response model and benchmarked against the laboratory measurements. This response model reproduces the pulse height curve shapes and performance hierarchy of the measured detectors well. Additionally, the simulated time response curves exhibit reasonable agreement with counterpart measurements. These measurements and simulation results together suggest that a large volume EJ-299-49 plastic scintillator that is mated to an ultraviolet transmitting plastic light guide and relies on total internal reflection for light transport to the photocathode will provide significantly faster timing, increase pulse height by a factor of 3, and require 25% less γ rays to meet experimental needs, compared to the candidate detector design developed prior to the availability of this response model. Application of this response model and associated performance metrics in the context of large-scale design sensitivity studies is also discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions

Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.

36 MATERIALS SCIENCE↗

Micrometer: Micromechanics transformer for predicting full field mechanical responses of heterogeneous materials

Predicting mechanical responses of heterogeneous materials across scales remains a significant challenge. Traditional computational methods often struggle with complex and multiscale nature of these materials, limiting their effectiveness in real-world applications. Here, in this paper, we introduce Micrometer, a vision transformer based deep learning model designed to predict full field mechanical responses of heterogeneous materials, bridging the gap between computer vision and solid mechanics problems. We show that Micrometer, trained on a large-scale high-resolution dataset of 2D fiber-reinforced composites, can achieve state-of-the-art performance in predicting microscale strain fields across a wide range of material properties and loading conditions. Our model demonstrates accuracy and computational efficiency in applications such as computational homogenization and multiscale modeling, reducing computational time by up to two orders of magnitude compared to conventional numerical solvers while maintaining less than 1 % errors in predicting macroscale stress fields. Furthermore, we showcase Micrometer’s adaptability through transfer learning experiments on new materials with limited data, highlighting its potential to tackle diverse scenarios in computational solid mechanics. These results represent a significant step towards AI-driven innovation in materials science, addressing the limitations of traditional numerical methods and paving the way for more efficient simulations of heterogeneous materials across various industrial applications.

Composite materials↗