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 163 records · Page 9

High-pressure hydrogen decompression in sulfur crosslinked elastomers

High-pressure storage and cyclic (de)pressurization of hydrogen gas is known to result in degradation and failure of gas canisters, hoses, linings, and O-rings as the relatively small hydrogen molecule can readily permeate most materials. Hence, identifying material compositions that are less susceptible to hydrogen-induced damage is of significant importance to the hydrogen energy infrastructure. Here, we use classical atomistic molecular dynamics simulations to study hydrogen exposed ethylene-propylene-diene monomer (EPDM) rubber, an elastomer typically used in O-rings. We make chemical modifications to the model by adjusting the crosslink density and report on gas solubility, diffusivity, and molecular restructuring in response to rapid decompression. Our simulations indicate that increases in crosslink density can reduce volumetric expansion during decompression and result in smaller free volume pore sizes. However, these favorable properties for sealing materials come with a tradeoff. At pressure, crosslinks introduce extra free volume, providing potential sites for gas localization, the precursor to cavitation-induced failure.

08 HYDROGEN↗

3D Printing of Cement-Based Materials Using Seawater for Simulated Marine Environments

Global demand for adaptable and rapidly deployable construction solutions in offshore, coastal, and fluvial environments continues to rise, driven by pressing needs to develop energy platforms, improve coastal resilience, and support emergency response in the face of natural disasters. Increased investment in human-made coastal infrastructure, such as piers, support structures for power lines, offshore wind farms, and seawall protection systems, further underscores this trend. This study investigates the development of printable concrete mixtures for underwater environments using seawater as a replacement for freshwater, using a 3D printing syringe-based extrusion system. The effect of seawater addition and the printing medium (in air vs. underwater) was assessed via rheological and mechanical performance characterization. The results indicate rheological properties are favorable for seawater adoption by producing mixtures with higher yield stress and viscosity with the same levels of admixtures used for freshwater. Seawater-based mixtures demonstrated superior dimensional stability compared to freshwater counterparts, maintaining cross-sectional geometry, while compressive strength results showed no statistical differences between in-air and underwater samples. However, flexural strength was significantly influenced by geometry and printing medium. These findings establish critical rheological parameters for printable underwater mixtures and highlight the need for optimized curing strategies and layer bonding techniques to improve interfacial strength in underwater 3D printing applications.

36 MATERIALS SCIENCE↗

Characterization of boundary precipitation in a heavy ion irradiated tungsten heavy alloy under the simulated fusion environment

In the concerted effort to identify materials capable of surviving the adverse environment of a fusion reactor interior, tungsten heavy alloys have been put forth as leading candidates. Experimental trials and behavioral studies have yielded positive results for their adoption by taking advantage of the alloy’s unique balance of inherently high fracture toughness and low sputter yields; yet due to their relative novelty in the fusion community, there remains a lack of understanding on the response of these materials to the extended high temperature irradiation environment of the reactor interior. Further, to alleviate this issue and provide the necessary data on the behavior of tungsten heavy alloys to the simulated fusion environment, a 90W-7Ni-3Fe alloy has been subjected to elevated temperature sequential Ni + and He + ion irradiations to mimic the expected displacement damage and gas production expected after five years of service as a plasma facing material component. Atomic-scale structural analyses and nanoscale chemical mapping have identified the formation of two distinct precipitation structures, a surface localized η-carbide and an a-W 2 C type tungsten carbide, both of which appear to originate at the bi-phase interface between W and the ductile phase. This irradiation assisted and induced precipitate formation respectively is anticipated to adversely affect these materials by selective embrittlement at bi-phase interfaces leading to a reduction in the material’s overall fracture toughness during prolonged high temperature irradiation. It is asserted that any potential exposure to C during fusion reactor operational service should be minimized as to prevent the formation of these phases.

36 MATERIALS SCIENCE↗

Comparison of simulated and experimentally measured detector impulse responses to pulsed x-rays

The Gamma Array Simulation Toolkit (GAST) is a multi-physics software stack comprised of Geant4, MATLAB, and experimentally measured photomultiplier tube (PMT) impulse response waveforms that is used to model and inform the end-toend performance of radiation detector systems. In order to validate the ability of GAST to predict the impulse response of a modeled detector configuration operated in current mode, experimental measurements using the Nevada National Security Site Transformational Diagnostics and Imaging 2 MeV endpoint energy electron linear accelerator (linac) with a tungsten target were performed. Seven different detector configurations consisting of unique combinations of scintillator material, size, reflectivity conditions, presence or absence of light guide, and PMT were evaluated. The resultant impulse response of each detector configuration was recorded and compared to a GAST simulation that emulated the respective linac experiment conditions. The comparison of the experimental and simulated detector impulse responses indicated that GAST can reliably predict detector impulse responses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Comparison of Simulated and Experimentally Measured Detector Impulse Responses to Pulsed X-rays

The Gamma Array Simulation Toolkit (GAST) is a multi-physics software stack comprised of Geant4, MATLAB, and experimentally measured photomultiplier tube (PMT) impulse response waveforms that is used to model and inform the end-to-end performance of radiation detector systems. In order to validate the ability of GAST to predict the impulse response of a modeled detector configuration operated in current mode, experimental measurements using the Nevada National Security Site Transformational Diagnostics and Imaging 2 MeV endpoint energy electron linear accelerator (linac) with a tungsten target were performed. Seven different detector configurations consisting of unique combinations of scintillator material, size, reflectivity conditions, presence or absence of light guide, and PMT were evaluated. The resultant impulse response of each detector configuration was recorded and compared to a GAST simulation that emulated the respective linac experiment conditions. The comparison of the experimental and simulated detector impulse responses indicated that GAST can reliably predict detector impulse responses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Unstructured mesh-based multi-physics reactor transient simulation with iMC code

In this study, we present a time-dependent Monte Carlo simulation result about a complex fuel-loaded light water reactor system with thermal-hydraulic feedback using the KAIST iMC code. The intra-pin temperature distribution is calculated based on Monte Carlo-tallied detailed power distribution and the finite element heat transfer method. The temperature effect of fuel, absorber, and coolant are adequately reflected into the dynamic response by the on-the-fly Doppler broadened cross-section correction. For the time-dependent Monte Carlo transport simulation, the predictor-corrector quasi-static Monte Carlo method is used for this study. We tested the coupled analysis framework with a 17-by-17 fuel assembly loaded with the centrally-shielded burnable absorber (CSBA) fuel element. The temperature feedback effect on the material cross-section and the moderator density incurred negative reactivity in response to the external positive reactivity insertion, showing the inherent stability of the reactor system. The pin-wise power and temperature change during the system transient are presented and discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Humidity Disrupts Structural and Chiroptical Properties of Chiral 2D Perovskites

Chiral two-dimensional (2D) hybrid organic–inorganic metal halide perovskite semiconductors have emerged as an exceptional material platform with many design opportunities for spintronic applications. However, a comprehensive understanding of changes to the crystal structure and chiroptical properties upon exposure to atmospheric humidity has not been established. We demonstrate phase degradation to the 1D (MBA)PbI 3 (MBA = methylbenzylammonium) and the hypothetical (MBA) 3 PbI 5 ·H 2 O hydrate phases, accompanied by a reduction and disappearance of the chiroptical response. Firstprinciple simulations show that water molecules preferentially locate at the interface between the organic cations and the inorganic framework, thereby disrupting the hydrogen bonding, impacting both the structural chirality and stability of the material. These findings provide critical insights into phase degradation mechanisms and their impact on chiroptical activity in chiral 2D perovskites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shock compression response of an insensitive high explosive single crystal: 1,1-diamino-2,2-dinitroethene (FOX-7)

Insensitive high explosives (IHEs) – that do not compromise performance – are of immense interest as a safer alternative to conventional high explosives, such as pentaerythritol tetranitrate (PETN) and cyclotrimethylenetrinitramine (RDX). Despite the strong interest in using IHEs, shock compression experiments on IHE single crystals have not been reported. To address this need, plate impact experiments were conducted to measure wave profiles in 1,1-diamino-2,2-dinitroethene (FOX-7) single crystals – a representative IHE crystal – shocked to 21 GPa longitudinal stress. Particle velocity histories, measured using laser interferometry, show a clear two-wave structure (elastic-inelastic response) at modest stresses (<3.8 GPa). Wave profiles at higher stresses show a single (overdriven) wave. Measured shock velocities and wave profiles provide accurate Hugoniot data to 21 GPa. The measured wave profiles to 21 GPa show no sign of energy release due to chemical decomposition and constitute the first demonstration of an IHE single crystal insensitivity under plane shock compression. Numerical simulations using a phenomenological material model developed for FOX-7 showed good agreement with the measured wave profiles. The experimental results and continuum simulations presented here constitute a significant first step in gaining insight into the shock compression response of IHE single crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A modular methodology for time-domain stochastic seismic wave propagation

Presented here is a modular methodology for time-domain stochastic seismic wave propagation analysis. Presented methodology is designed to analyse uncertain seismic motions as an input, propagating through uncertain material. Traditional approach for uncertain wave propagation relies on models that include deep bedrock, local soil site, and their random process and random field information. Such models can become quite large and computationally intractable. The modular approach proposed herein features two step approach that allows separate consideration of the deep bedrock and local site along with corresponding random field information. In this work, the first step considers an auxiliary stochastic motions problem in the bedrock. Stochastic local site response can then be simulated in a reduced domain within certain depth from the surface. Application of uncertain seismic motions at depth, for local uncertain site response is done using stochastic effective forces developed through the Domain Reduction Method. By using Hermite polynomial chaos expansion to represent the non-Gaussian random field of material parameters and non-stationary random process of seismic motion, the proposed modular methodology is formulated using intrusive stochastic Galerkin approach, as seen in the Stochastic Elastic–Plastic Finite Element Method (SEPFEM). Developed modular methodology is illustrated using a 1-D stochastic seismic wave propagation analysis with three cases, and simulation results are also verified with results from conventional approach.

58 GEOSCIENCES↗

Generation of Synthetic Data for a Radiation Detection Algorithm Competition

This paper details the generation of synthetic radiation data using large-scale Monte Carlo transport models to evaluate radiation search detection algorithms. Modular 3-D Monte Carlo models spanning multiple city blocks were constructed, loosely based on downtown Knoxville, TN, containing buildings composed of multiple materials (brick, granite, and concrete), sidewalks, a four-lane road, side streets, parking lots, and grassy fields. Background and simulated source detector response calculations from these models were used to create synthetic list mode data sets for a 2"×4"×16" NaI(Tl) detector moving through a city street at a constant speed. For the background simulations, major isotopes were computed individually so that background composition and variability could be computed efficiently outside Monte Carlo. The source detector response included six simulated sources placed at 15 source locations. Detector response was developed to be periodic through the city street so that a detector path could begin at one end of the model and wrap around to the other. This framework allowed for the creation of diverse data sets, each with its own unique background and simulated source detector response. Synthetic data allows for high-quality labels, which are useful in developing data-driven radiation detection algorithms. This methodology was used to create synthetic data sets which were released as part of a public data competition to spur the development of new radiation detection algorithms for radiological search applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Experimental Study of Mechanical Response of Polyurethane Rubber at Varying Strain Rates for W80-4 System

This report describes an experimental study to determine the mechanical behavior of the polyurethane rubber material that was used in the W80-4 MCC Shock/Breach Phase 1 test series. Compression experiments were conducted on cylindrical specimens over a wide range of loading rates to characterize the material over the range of strain rates that were experienced in the shock/breach testing. Additionally, specimen diameter was varied to determine the effect of confinement on the material response and was found to be significant. This data is used to populate a material model to enable accurate analyses and finite element simulations of the shock/breach test series.

36 MATERIALS SCIENCE↗

PRISMS-Fatigue computational framework for fatigue analysis in polycrystalline metals and alloys

Abstract The PRISMS-Fatigue open-source framework for simulation-based analysis of microstructural influences on fatigue resistance for polycrystalline metals and alloys is presented here. The framework uses the crystal plasticity finite element method as its microstructure analysis tool and provides a highly efficient, scalable, flexible, and easy-to-use ICME community platform. The PRISMS-Fatigue framework is linked to different open-source software to instantiate microstructures, compute the material response, and assess fatigue indicator parameters. The performance of PRISMS-Fatigue is benchmarked against a similar framework implemented using ABAQUS. Results indicate that the multilevel parallelism scheme of PRISMS-Fatigue is more efficient and scalable than ABAQUS for large-scale fatigue simulations. The performance and flexibility of this framework is demonstrated with various examples that assess the driving force for fatigue crack formation of microstructures with different crystallographic textures, grain morphologies, and grain numbers, and under different multiaxial strain states, strain magnitudes, and boundary conditions.

Chemistry↗

Parallel Simulation of Beam Dynamics in Particle Accelerators [Slides]

Particle accelerators are among the most versatile and important tools of scientific discovery. The Nation's accelerators are responsible for a wealth of advances in materials science, chemistry, the biosciences, particle physics, and nuclear physics. They also have important applications to national security, the environment, energy, medicine, and on the quality of people's lives. LANL has a long history of making pioneering contributions to Accelerator Science including key contributions to the field of Computational Accelerator Physics. These include the development of early beam dynamics codes with space charge (such as PARMILA and PARMELA), the development of rf cavity codes and magnet codes (including Poisson and Superfish), and the development and distribution of codes to the accelerator community through the Los Alamos Accelerator Code Group. LANL researchers also helped pioneer the development of massively parallel space-charge codes. In project t22_accelsim we have moved beyond electrostatic models of collective effects (i.e., solving the Poisson equation in the bunch frame) to fully electromagnetic models based on the Lienard-Wiechert formalism. This approach enables the large-scale simulation of radiation production and collective effects in high brightness electron beams. This is highly relevant to LANL given its future goal of developing an X-ray Free Electron Laser (XFEL). It also directly impacts a LANL LDRD project to develop an undulator-based non-invasive beam profile monitor for beams created in laser-plasma accelerator systems.

43 PARTICLE ACCELERATORS↗

Antiferromagnet thickness dependence and rotatable spins in exchange biased CoO/Fe films

The emergence of exchange bias and coercivity enhancement has been investigated in epitaxial CoO/Fe films with varied antiferromagnet (AF) thicknesses, even smaller than the critical value where the frozen CoO spins are detectable. Vector magnetometry and first-order reversal curve (FORC) measurements reveal different CoO thickness dependence of the exchange bias and coercivity enhancement, including the evolution of magnetization reversal from a high coercivity, low bias phase due to rotatable CoO moments to a high bias, low coercivity phase due to frozen CoO moments. Furthermore, the AF domain state is found to be metastable, which can be reoriented by external and exchange fields prior to the appearance of frozen spins, pointing to a generic origin of the training effect. Monte Carlo simulations show that the AF anisotropy energy barrier and the rotatable spins induced by magnetic field and exchange interaction at the interface are responsible for the observed effects.

36 MATERIALS SCIENCE↗

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗

A sensing material-free and simple readout MEMS sensor for detecting Helium

Abstract In this work, we report a method that enables a standard electrostatic MEMS device to perform complex sensing functionalities, such as detecting the presence of helium without a sensing material or a conditioning circuit. Helium is a noble, odorless, non-reactive gas that is very challenging to detect. It is used in critical applications such as storing nuclear fuel waste inside a dry cask. In these applications, its leakage from the dry cask may indicate the cask's safe operation's degradation. A departure from the common practice of exciting the MEMS around its mechanical resonance, the method is based on exciting the MEMS around its electrical resonance circuit. This method shows that the tiny difference between the air dielectric constant (1.00059) and helium (1.000067) corresponding to only a few Femtofarad level capacitances produces a 25 mV difference without a conditioning circuit. Simulation results confirmed those findings and explored the sensor response at different operation conditions. This method eliminates the need for a heated microstructure and the need for absorption material. This method is not limited to gas sensing. It can be applied to other sensing mechanisms, such as acceleration and pressure measurements, and eliminate the complex circuit to read small capacitance in these applications.

Mohaidat, Sulaiman↗

A multiscale anisotropic polymer network model coupled with phase field fracture

The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.

42 ENGINEERING↗

What dictates soft clay-like lithium superionic conductor formation from rigid salts mixture

Abstract Soft clay-like Li-superionic conductors, integral to realizing all-solid-state batteries, have been recently synthesized by mixing rigid-salts. Here, through computational and experimental analysis, we clarify how a soft clay-like material can be created from a mixture of rigid-salts. Using molecular dynamics simulations with a deep learning-based interatomic potential energy model, we uncover the microscopic features responsible for soft clay-formation from ionic solid mixtures. We find that salt mixtures capable of forming molecular solid units on anion exchange, along with the slow kinetics of such reactions, are key to soft-clay formation. Molecular solid units serve as sites for shear transformation zones, and their inherent softness enables plasticity at low stress. Extended X-ray absorption fine structure spectroscopy confirms the formation of molecular solid units. A general strategy for creating soft clay-like materials from ionic solid mixtures is formulated.

25 ENERGY STORAGE↗