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 271 records · Page 15

MELCOR

MELCOR is a fully integrated, engineering-level computer code for modeling the progression of severe accidents in light water reactors (LWR) at nuclear power plants and nuclear fuel cycle facilities. Originally developed to assess severe accidents following Three Mile Island, MELCOR’s flexible modeling framework has enabled it to be applied to safety assessments of a much broader range of nuclear power reactor designs and other types of nuclear facilities processing radioactive material. Further, MELCOR can model a broad spectrum of severe accident phenomena such as thermal-hydraulic response in a reactor coolant system; core heat-up, degradation, and relocation; and transport behavior in both boiling water and pressurized water reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Library of Advanced Materials for Engineering (LAMÉ) (V.5.24)

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options an

36 MATERIALS SCIENCE↗

Multiscale Thermomechanical Damage Model with Internal State Variables for Ceramic Matrix Composites

A length scale-dependent matrix damage model has been developed to simulate the response of ceramic matrix composites (CMCs). This model captures damage initiation and propagation in the brittle matrix material by employing internal state variable (ISV) theory within a multiscale modeling framework to obtain damaged matrix stress-strain constitutive relationships at each length scale. Material degradation due to matrix cracking is captured by a damage variable determined using fracture mechanics and the self-consistent scheme. An additional state variable is introduced to capture the nucleation and growth of matrix porosity, a key mechanism contributing to matrix nonlinearity. Porosity occurs as a result of material diffusion around grain boundaries and the growth of the porosity state variable is related to the material entropy dissipation and the volumetric strain. The nonlinear predictive capabilities of the material model are demonstrated for monolithic silicon carbide, unidirectional (UD) carbon fiber/silicon carbide matrix (C/SiC) CMC, and five-harness satin (5HS) woven C/SiC CMC.

Skinner, Travis D.↗

High strain-rate strength response of single crystal tantalum through in-situ hole closure imaging experiments

The properties of crystalline materials often depend on directionality and operating conditions. Specifically, the strength of materials can depend anisotropically on crystal direction and the loading condition. To probe these effects, a preliminary series of high strain-rate (> 105/s) strength plate-impact hole closure experiments were performed on high purity single crystal Tantalum cubes. The orientation of the single crystals with respect to impact/loading were varied to provide data to inform crystal plasticity modeling efforts. The experiments consist of in-situ high-resolution X-ray radiographic imaging of the hole collapse under dynamic compression conditions to infer the material strength via its resistance to closure at increasing levels of plastic strain. The experiments are compared against hydrocode simulation predictions. Here, a comparison with simple elastic perfectly plastic strength model predictions is presented to elucidate the response of the different crystal orientations at high strain-rate and large plastic strains.

36 MATERIALS SCIENCE↗

Imposing correct jellium response is key to predict the density response by orbital-free DFT

Orbital-free density functional theory constitutes a computationally highly effective tool for modeling electronic structures of systems ranging from room-temperature materials to warm dense matter. Its accuracy critically depends on the employed kinetic energy (KE) density functional, which has to be supplied as an external input. In this work we consider several nonlocal and Laplacian-level KE functionals and use an external harmonic perturbation to compute the static density response at T = 0 K in the linear and beyond-linear response regimes. We test for the satisfaction of exact conditions in the limit of uniform densities and for how approximate KE functionals reproduce the density response of realistic materials (e.g., Al and Si) against the Kohn-Sham DFT reference, which employs the exact KE. The results illustrate that several functionals violate exact conditions in the uniform electron gas (UEG) limit. We find a strong correlation between the accuracy of the KE functionals in the UEG limit and in the strongly inhomogeneous case. This empirically demonstrates the importance of imposing the limit of UEG response for uniform densities and validates the use of the Lindhard function in the formulation of kernels for nonlocal functionals. This conclusion is substantiated by additional calculations for bulk aluminum (Al) with a face-centered cubic (fcc) lattice and silicon (Si) with an fcc lattice, body-centered cubic (bcc) lattice, and semiconducting crystal diamond state. The analysis of fcc Al, and fcc as well as bcc Si data follows closely the conclusions drawn for the UEG, allowing us to extend our conclusions to realistic systems that are subject to density inhomogeneities induced by ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of nickel on the kinematic stability of retained austenite in carburized bearing steels – In-situ neutron diffraction and crystal plasticity modeling of uniaxial tension tests in AISI 8620, 4320 and 3310 steels

The presence of kinematically metastable retained austenite in the microstructure of bearing components can significantly affect the macro and micro-mechanical material response. In the present work, the influence of Ni on the stability of the retained austenite within three different grades of high carbon bearing steel using in-situ neutron diffraction is investigated. For the first time, the results show that presence of Ni increases the stability of the austenite in the elastic regime whereas the transformation rate remains unaffected. Crystal plasticity finite element (CPFE) modeling was used to study the deformation in these three steels and shows that the predominant factor causing the difference in mechanical behavior of these steels is the austenite stability. Finally, the elastic and plastic response of the matrix martensite was found to be identical among all specimens while the austenite demonstrates similar elastic behavior but remarkably different stabilities under monotonic loading.

36 MATERIALS SCIENCE↗

Computational Model for Microballistic Perforation of Multilayer Graphene

The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. The computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.

36 MATERIALS SCIENCE↗

Physics augmented machine learning discovery of composition-dependent constitutive laws for 3D printed digital materials

Multi-material 3D printing, particularly through polymer jetting, enables the fabrication of digital materials by mixing distinct photopolymers at the micron scale within a single build to create a composite with tunable mechanical properties. Here, this work presents an integrated experimental and computational investigation into the composition-dependent mechanical behavior of 3D printed digital materials. We experimentally characterize five formulations, combining soft and rigid UV-cured polymers under uniaxial tension and torsion across three strain and twist rates. The results reveal nonlinear and rate-dependent responses that strongly depend on composition. To model this behavior, we develop a physics-augmented neural network (PANN) that combines a partially input convex neural network (pICNN) for learning the composition-dependent hyperelastic strain energy function with a quasi-linear viscoelastic (QLV) formulation for time-dependent response. The pICNN ensures convexity with respect to strain invariants while allowing non-convex dependence on composition. To enhance interpretability, we apply $L_0$ sparsification. For the time-dependent response, we introduce a multilayer perceptron (MLP) to predict viscoelastic relaxation parameters from composition. The proposed model accurately captures the nonlinear, rate-dependent behavior of 3D printed digital materials in both uniaxial tension and torsion, achieving high predictive accuracy for interpolated material compositions. This approach provides a scalable framework for automated, composition-aware constitutive model discovery for multi-material 3D printing.

Constitutive modeling↗

Complex Oxides under Simulated Electric Field: Determinants of Defect Polarization in AB O 3 Perovskites

Abstract Polarization of ionic and electronic defects in response to high electric fields plays an essential role in determining properties of materials in applications such as memristive devices. However, isolating the polarization response of individual defects has been challenging for both models and measurements. Here the authors quantify the nonlinear dielectric response of neutral oxygen vacancies, comprised of strongly localized electrons at an oxygen vacancy site, in perovskite oxides of the form AB O 3 . Their approach implements a computationally efficient local Hubbard U correction in density functional theory simulations. These calculations indicate that the electric dipole moment of this defect is correlated positively with the lattice volume, which they varied by elastic strain and by A‐site cation species. In addition, the dipole of the neutral oxygen vacancy under electric field increases with increasing reducibility of the B‐site cation. The predicted relationship among point defect polarization, mechanical strain, and transition metal chemistry provides insights for the properties of memristive materials and devices under high electric fields.

36 MATERIALS SCIENCE↗

RADAI: A Large-Scale Realistic Dataset for Radiation Detection Algorithm Development

Open, realistic datasets are essential for developing and benchmarking radiation detection algorithms, yet they remain scarce. The Radiological Anomaly Detection and Identification (RADAI) project was develop to create datasets that meet the training and testing needs for sophisticated radiation detection algorithms. The RADAI dataset is a large-scale synthetic resource that integrates high-fidelity Monte Carlo simulations with realistic urban scenarios to capture both background variability and source signatures. RADAI models construction-material NORM, people and vehicles, urban clutter, and dynamic environmental effects such as cosmic-ray and rain-induced transients, and they provide list-mode detector data with motion and response modeling suitable for algorithm training and evaluation. The RADAI project resulted in three publicly-released complementary datasets together with an online scoring portal for standardized performance assessment and an open software toolkit that supports data access, augmentation, model development, and evaluation. These resources enable reproducible comparisons across methods and promote rigorous studies at the scale required by contemporary machine learning. By grounding algorithm development in realistic, well-documented conditions, RADAI supports progress toward more robust detection, identification, and localization in complex urban environments.

Ghawaly, James M. [Division of Computer Science an↗

Thermal expansion of plasma-exposed tungsten

We report results from a systematic analysis of thermal expansion of plasma-exposed tungsten based on molecular-dynamics simulations using models of tungsten with distributions of helium (He) bubbles in the tungsten matrix. We distinguish between two approaches of filling the bubbles with He, where the amount of He in the bubble can or cannot vary with temperature. In the former case, the thermal expansion coefficient decreases monotonically with the porosity and He content of the tungsten matrix, while in the latter case, the thermal expansivity increases monotonically with increasing porosity and He content. The latter condition, where the He content in the bubble is determined at the implantation temperature and remains constant with varying temperature in the tungsten matrix, is consistent with He species transport in tungsten used as a plasma-facing component (PFC) in nuclear fusion reactors and implies the development of biaxial compressive thermal strains in the PFC material that contribute to accelerating the growth of a nanostructure on PFC tungsten surfaces. Furthermore our analysis advances the fundamental understanding of thermal expansion in PFC tungsten and contributes to the development of a thermophysical property database for properly incorporating effects of realistic heat loads into modeling the dynamical response of PFC tungsten under fusion reactor operating conditions.

36 MATERIALS SCIENCE↗

Towards predictive control of reversible nanoparticle assembly with solid-binding proteins

Although a broad range of ligand-functionalized nanoparticles and physico-chemical triggers have been exploited to create stimuli-responsive colloidal systems, little attention has been paid to the reversible assembly of unmodified nanoparticles with non-covalently bound proteins. Previously, we reported that a derivative of green fluorescent protein engineered with oppositely located silica-binding peptides mediates the repeated assembly and disassembly of 10-nm silica nanoparticles when pH is toggled between 7.5 and 8.5. We captured the subtle interplay between interparticle electrostatic repulsion and their protein-mediated short-range attraction with a multiscale model energetically benchmarked to collective system behavior captured by scattering experiments. Here, in this work, we show that both solution conditions (pH and ionic strength) and protein engineering (sequence and position of engineered silica-binding peptides) provide pathways for reversible control over growth and fragmentation, leading to clusters ranging in size from 25 nm protein-coated particles to micrometer-size aggregate. We further find that the higher electrolyte environment associated with successive cycles of base addition eventually eliminates reversibility. Our model accurately predicts these multiple length scales phenomena. The underpinning concepts provide design principles for the dynamic control of other protein- and particle-based nanocomposites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A microstructure-based modeling approach to predict the mechanical properties of Zr alloy with hydride precipitates

In nuclear reactors, hydrides can form in fuel cladding due to hydrogen absorption in ZIrcaloy and cause embrittlement. This work presents a microstructure-based finite element model to predict the stress-strain response of Zircaloy containing hydrides. Quantitative microstructural details extracted from scanning electron microscopy (SEM) images were used to generate heterogeneous microstructures including the morphology and spatial distribution of hydrides. The constitutive material model for zircaloy in this study is based on crystal plasticity theory which considers the hexagonal close-packed (HCP) atomic structure of Zircaloy material. The hydrides were modeled as brittle material along with a damage model. Hydride formation inside the zircoloy matrix results in residual stress. This phenomenon is also captured in this model. A parametric study has been conducted to understand the effect of volume fraction, orientation, and lamellae thickness of the hydride phase on the mechanical properties of the overall material.

Kulkarni, Shank S.↗

BISON capability and validation for U3Si2, Cr2O3-doped UO2, FeCrAl, and Cr-coated Zircaloy ATF concepts

Development, implementation, and validation of material and behavior models for accident tolerant fuel (ATF) concepts in the BISON fuel performance code began in 2014 in response to the events that occurred at the Fukushima Daichii nuclear power plant in March 2011. Early on, the focus was on U3Si2 fuel and FeCrAl cladding as part of a high-impact problem through the Nuclear Energy and Advanced Model- ing Simulation program [1]. Then, developments for Cr2O3- doped UO2 fuel, as well as SiC-SiC and Cr-coated Zircaloy cladding began based upon industry interests. Near the end of fiscal year 2018 the Consortium for Advanced Simulation of Light Water Reactors took over further Accident Tolerant Fuel (ATF) work in BISON in support of Nuclear Regulatory Commission (NRC) engagement. Discussions with the NRC clarified their list of priority fuel and cladding concepts, which included Cr2O3-doped UO2 and U3Si2 fuels, as well as Cr-coated Zircaloy and FeCrAl claddings. Further, reports similar in form to NUREG/CR-7024 [2] that was developed for traditional light water reactor (LWR) materials, UO2 fuel with Zircaloy claddings were created for the priority fuel and cladding concepts. In this paper select details from the reports entitled "ATF material model development and validation for priority fuel concepts" [3] and "ATF material model development and validation for priority cladding concepts" [4] are presented. First, a brief overview of the available models and their associated uncertainty are provided for the priority fuel and cladding concepts. Then, for the concepts with experimental data, comparisons between simulation and experiment are shown taking into account forward uncertainty quantification (UQ) techniques. For the concepts without publicly available experimental data, parametric studies incorporating forward UQ are completed to verify claims made about the materials’ behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Numerical and experimental evaluation of mechanical performance of the multifunctional energy storage composites

This work presents numerical simulation methods to model the mechanical behavior of the multifunctional energy storage composites (MESCs), which consist of a stack of multiple thin battery layers reinforced with through-the-hole polymer rivets and embedded inside carbon fiber composite laminates. MESC has been demonstrated through earlier experiments on its exceptional behavior as a structural element as well as a battery. However, the inherent complex infrastructure of the MESC design has created significant challenges in simulation and modeling. A novel homogenization technique was adopted to characterize the multi-layer properties of battery material using physics-based constitutive equations combined with nonlinear deformation theories to handle the interface between the battery layers. Second, mechanical damage and failure modes among battery materials, polymer reinforcements, and carbon fiber-polymer interfaces were characterized through appropriate models and experiments. The model of MESCs has been implemented in a commercial finite element code in ABAQUS. A comparison of structural response and failure modes from numerical simulations and experimental tests are presented. The results of the study showed that the predictions of elastic and damage responses of MESCs at various loading conditions agreed well with the experimental data. © 2021

Materials Science↗

Temperature-Assisted Piezoresponse Force Microscopy: Probing Local Temperature-Induced Phase Transitions in Ferroics

The combination of local heating and biasing at the tip-surface junction in temperature-assisted piezoresponse force microscopy (TPFM) opens a pathway for probing local temperature-induced phase transitions in ferroics, exploring the temperature dependence of polarization dynamics in ferroelectrics and potentially discovering coupled phenomena driven by strong temperature and electric field gradients. In this study we analyze the signal-formation mechanism in TPFM and explore the interplay between thermal- and bias-induced switching in model ferroelectric materials. Furthermore, we explore the contributions of the flexoelectric and thermopolarization effects to the local electromechanical response and demonstrate that the latter can be significant for “soft” ferroelectrics. These results establish a framework for the quantitative interpretation of TPFM observations, predict the emergence of nontrivial switching and relaxation phenomena driven by nonlocal thermal-gradient-induced polarization switching, and open a pathway for exploring the physics of thermopolarization effects in various noncentrosymmetric and centrosymmetric materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Formal Experimentation and Analysis of Handheld RFID Readers as a Tool for Nuclear Material Accounting

Commercial, off-the-shelf Radio Frequency Identification (RFID) systems have been successfully deployed for inventory tracking in numerous industries, but their viability in the tracking of complex environments containing nuclear material is less understood. Of primary interest in this setting is a RFID tracking system for nuclear material accounting which can reliably identify as many tags as possible with high accuracy, while also reducing an operator’s exposure to radiation. In this work, we develop a formal statistical approach to identify relevant handheld RFID reader settings which optimize tagging performance. To achieve this goal, we design a full factorial split-plot experiment for a static shelf configuration scene with 50 randomly placed nuclear material containers affixed with RFID tags. We use Bayesian inference to fit a second-order response surface model which expresses the probability of a successful match for each container as a function of the experimental factors. Such effects are allowed to vary across individual containers and the containers’ population in its entirety to estimate overall effects. Uncertainties of estimates and predictions are quantified via their corresponding posterior distributions. Following extensive model checking and validation, the fitted model is used to identify experimental factors which maximize matching probabilities at both the container-level and for the full shelf configuration scene. We also analyze sensitivity of performance to relevant factors.

25 ENERGY STORAGE↗

Modeling detonation with CartaBlanca simulations

The accurate modeling of high explosive (HE) detonation and the con- sequent large solid deformation, failure, plastic flow, porosity growth, and shock wave propagation is important because simulations can capture spatial and temporal features that experimental diagnostics cannot capture. However, the simulation of the explosive event poses challenges to a computational scientist. These include the accurate modeling of ductile damage, crack formation, plastic deformation, as well as physical and nu- merical instabilities. The material response can be history-dependent and subject to large material deformation. Our research simulates the impact of a high explosive (Detasheet) onto a tantulum metal plate. We performed the simulations using CartaBlanca at different mesh resolutions. We decided the study would be impactful if we perform the simulations with the Material Point Method. Differences were observed at the different mesh resolutions in velocity and nodal stress magnitude, so increased mesh resolutions may be required. In addition, the Discontinuous Galerkin method would be needed to account for the large gas expansion.

97 MATHEMATICS AND COMPUTING↗