SEARCH · Engineering Papers

Results for “materials simulations”

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.

SourceResearch connectionFrom yearThrough year

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

At least 451 records · Page 25

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Predicting initial dissolution rates using structural features from molecular dynamics simulations

Predicting chemical durability of glass materials is important for various applications from daily life such as drink glass and kitchen ware to advanced technologies such as nuclear waste disposal and biomedicine. In this work, we explored prediction of initial dissolution rate through structural features from molecular dynamics (MD) simulations for a wide range of glass compositions (total 28) including borosilicate and aluminosilicate glasses, ZrO 2 -containing and V 2 O 5 -containing boroaluminosilicate glasses. The initial dissolution rates (r 0 ) measured experimentally at 90 °C with varying solution conditions were correlated with structural features (e.g., polyhedral linkages and non-bridging oxygen species) obtained from MD simulations, either from this study or from literature. Since hydrolysis of the glass network through breaking of the network former linkages (e.g., Si-O-Si, Si-O-Al, etc.) is a critical step of network glass dissolution, the statistics of these linkages obtained from MD were correlated to r0 through linear regression, where the coefficient of determination (R 2 ) and root mean square error are found to be 0.949 and 0.681, respectively. This model was compared and discussed with existing models developed by various approaches including machine learning, the kinetic rate equation, topological constraint theory, and other descriptors from MD simulations. The discussion provides insights on future model improvements to predict glass dissolution. In addition, the impact of V 2 O 5 on the glass dissolution was examined in detail, implicating that the impact is not the same across all glass compositions and test conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Simulants for Testing the Blend Can Loading System

Eleven simulants have been provided to the Savannah River National Laboratory (SRNL) Mechanical Systems & Custom Equipment Development (MSCED) organization to further test the capabilities of the Blend Can Loading System (BCLS). The basis for the materials used in these simulants are based on SRNL-STI-2020-00503, “Simulant Recommendation for the Blend Can Loading System” (Ref. i). The 11 different simulants provided to MSCED are summarized in Table ES-1. This table provides a basis of why they were selected and what properties were specifically targeted, given the composition of the procured materials.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

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↗

DOE OSTI2021

Dynamics of laser-generated magnetic fields using long laser pulses

In this work, we report on the experimental investigation of magnetic field generation with a half-loop gold sheet coil driven by long-duration (10 ns) and high-power (0.5 TW) laser pulses. The amplitude of the magnetic field was characterized experimentally using proton deflectometry. The field rises rapidly in the first 1 ns of laser irradiation, and then increases slowly and continuously up to 10 ns during further laser irradiation. The transient dynamics of current shape were investigated with a two-dimensional (2D) numerical simulation that included Ohmic heating of the coil and the resultant change of electrical resistivity determined by the coil material temperature. The numerical simulations show rapid heating at the coil edges by current initially localized at the edges. This current density then diffuses to the central part of the sheet coil in a way that depends both on normal current diffusion as well as temporal changes of the coil resistance induced by the Ohmic heating. The measured temporal evolution of the magnetic field is compared with a model that determines a solution to the coil current and voltage that is consistent with a plasma diode model of the drive region and a 2D simulation of current diffusion and dynamic resistance due to Ohmic heating in the laser coil.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Radioactive Waste Sludge Washing and Demonstration of the Nitric-Glycolic Acid Flowsheet for Sludge Batch 10 Qualification

For each sludge batch that is processed in the Defense Waste Processing Facility (DWPF), the Savannah River National Laboratory (SRNL) performs qualification testing to demonstrate that the sludge batch (SB) is processible. During processing of SB9, DWPF will be transitioning from the Nitric-Formic Acid (NFA) flowsheet to the Nitric-Glycolic Acid (NGA) flowsheet. Thus, the qualification of SB10 was requested to only be performed using the NGA flowsheet. In order to qualify the batch for the NGA flowsheet, Sludge Receipt and Adjustment Tank (SRAT) and Slurry Mix Evaporator (SME) cycles, designated SC-19, were performed using SB10 Tank 51 sample material. SRNL received Tank 51 material in the midst of Tank Farm washing. SRNL continued the washing in the SRNL Shielded Cells. The SRNL process included the addition of Sodium Reactor Experiment (SRE) material from H Canyon, simulating the transfer of SRE from H Canyon to Tank 51 during Tank Farm washing. The washed SB10 Tank 51 material, with SRE, was characterized prior to flowsheet qualification testing.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Material Interactions in Severe Accidents – Benchmarking the MELCOR V2.2 Eutectics Model for a BWR-3 MARK-I Station Blackout: Part I – Single Case Analysis

Here in this analysis, the two material interaction models available in the MELCOR code are benchmarked for a severe accident at a BWR under representative Fukushima Daiichi boundary conditions. This part of the benchmark investigates the impact of each material interaction model on accident progression through a detailed single case analysis. It is found that the eutectics model simulation exhibits more rapid accident progression for the duration of the accident. The slower accident progression exhibited by the interactive materials model simulation, however, allows for a greater degree of core material oxidation and hydrogen generation to occur, as well as elevated core temperatures during the ex-vessel accident phase. The eutectics model simulation exhibits more significant degradation of core components during the late in-vessel accident phase – more debris forms and relocates to the lower plenum before lower head failure. The larger debris bed observed in the eutectics model simulation also reaches higher temperatures, presenting a more significant thermal challenge to the lower head until its failure. At the end of the simulated accident scenario, however, core damage is comparable between both simulations due to significant core degradation that occurs during the ex-vessel phase in the interactive materials model simulation. A key difference between the two models’ performance is the maximum temperatures that can be reached in the core and therefore the maximum ΔT between any two components. When implementing the interactive materials model, users have the option to modify the liquefaction temperature of the ZrO 2 -interactive and UO 2 -interactive materials as a way to mimic early fuel rod failure due to material interactions. Through modification of the liquefaction of high melting point materials with significant mass, users may inadvertently limit maximum core temperatures for fuel, cladding, and debris components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Towards exact finite temperature electronic structure in solids and molecules (Final Technical Report)

This report describes the University of Iowa portion of a project that is now continuing at Michigan State University. We are developing novel methods, algorithms, and software to enable simulations of molecules and materials at high temperature. This is a key challenge in chemistry and materials science. By refining an approach called Density Matrix Quantum Monte Carlo (DMQMC), we developed faster and more accurate ways to conduct these simulations. These advances will help us understand how temperature affects the behavior of electrons, chemical bonds, and phase transitions in solids and molecules. These breakthroughs are especially important for applications where light and heat drive chemical reactions, superconductivity, and materials used in energy and sensing. In addition, this project involved the development of the open-source HANDE-QMC software package, supporting the broader community in benchmarking and developing finite-temperature electronic structure methods.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Equivalent Fundamental-Mode Source Simulations for Spherical Uranium and Plutonium Systems [Abstract]

The equivalent fundamental-mode source refers to a “source” which is identically distributed in space, energy, and angle as that of a fundamental-mode fission source distribution. It has previously been shown that the system fixed source multiplication ($M_{fs}$) relates to the effective multiplication factor ($\kappa_{eff}$) as $M_{fs} = \frac{g^*}{1-k_{eff}}$ The term $\mathcal{g}$* therefore relates the system source distribution to the system equivalent fundamental-mode source. As previously shown, the $\mathcal{g}$* term is approximately 1 for deeply subcritical systems and then diverges as the system multiplication increases. For a system with a point source located in the center of an assembly (such a sphere), $\mathcal{g}$* will be greater than 1; this happens because there is a higher probability of induced fission when all of the starter neutrons originate from the center of the sphere. For a uniform source, $\mathcal{g}$* is less than 1, but does not diverge from 1 as dramatically as seen by a point source. One reason that $\mathcal{g}$* is an important parameter is that it is used in measurement methods to infer reactor kinetics parameters ($β_{eff}$ being one example). This work will discuss multiple ways to simulate $\mathcal{g}$* using MCNP®6.2, which include use of traditional methods versus a single input file method that was recently published. A discussion related to methods to simulate $M_{fs}$ will be presented. Spherical systems of Highly Enriched Uranium (HEU) will be compared with previous works. In addition, results of plutonium including multiple reflector materials will be investigated. Simulated systems will include both parametric studies as well as configurations used in recent experiments performed at the National Criticality Experiments Research Center (NCERC). These studies will include subtleties associated with system geometry that have not been previously shown. Last, this work will investigate the relationship between $\mathcal{g}$* and reactor kinetics parameters (such as $β_{eff}$).

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Disorder effects in nitride semiconductors: impact on fundamental and device properties

Abstract Semiconductor structures used for fundamental or device applications most often incorporate alloy materials. In “usual” or “common” III–V alloys, based on the InGaAsP or InGaAlAs material systems, the effects of compositional disorder on the electronic properties can be treated in a perturbative approach. This is not the case in the more recent nitride-based GaInAlN alloys, where the potential changes associated with the various atoms induce strong localization effects, which cannot be described perturbatively. Since the early studies of these materials and devices, disorder effects have indeed been identified to play a major role in their properties. Although many studies have been performed on the structural characterization of materials, on intrinsic electronic localization properties, and on the impact of disorder on device operation, there are still many open questions on all these topics. Taking disorder into account also leads to unmanageable problems in simulations. As a prerequisite to address material and device simulations, a critical examination of experiments must be considered to ensure that one measures intrinsic parameters as these materials are difficult to grow with low defect densities. A specific property of nitride semiconductors that can obscure intrinsic properties is the strong spontaneous and piezoelectric fields. We outline in this review the remaining challenges faced when attempting to fully describe nitride-based material systems, taking the examples of LEDs. The objectives of a better understanding of disorder phenomena are to explain the hidden phenomena often forcing one to use ad hoc parameters, or additional poorly defined concepts, to make simulations agree with experiments. Finally, we describe a novel simulation tool based on a mathematical breakthrough to solve the Schrödinger equation in disordered potentials that facilitates 3D simulations that include alloy disorder.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Atomic scale etching of diamond: insights from molecular dynamics simulations

Diamond is a promising material for multiple applications in quantum information processing and sensing as well as applications in microelectronics. However, diamond devices can be limited by surface defects that compromise charge stability and spin coherence, among others. Improved strategies in plasma etching of diamond could play an important role in minimizing or eliminating these defects. In this work, we explore plasma-assisted atomic scale etching of diamond using argon ions (Ar + ), hydrogen ions (H + ) and hydrogen atoms (H). We employ classical molecular dynamics (MD) simulations and test several interatomic potentials based on the Reactive Empirical Bond Order (REBO) form with comparisons to a variety of published experimental results. We performed MD simulations of low-energy hydrogen ($\leqslant$50 eV) and argon ( $\leqslant$200 eV) ion bombardment of diamond surfaces. Ar + bombardment can be used to locally smooth initially rough diamond surfaces via the formation of an amorphous C layer, the thickness of which increases with argon ion energy. Subsequent exposure with hydrogen ions (or fast neutrals) will selectively etch this amorphous C layer, leaving the underlying diamond layer mostly intact if the H energy is maintained below about 10 eV. The simulations suggest that combining Ar + smoothing with selective, near threshold energy H removal of amorphous C can be an effective strategy for diamond surface engineering, leading to more reliable and sensitive diamond color center devices.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2023

Frustration- and doping-induced magnetism in a Fermi–Hubbard simulator

Geometrical frustration in strongly correlated systems can give rise to a plethora of novel ordered states and intriguing magnetic phases, such as quantum spin liquids. Promising candidate materials for such phases can be described by the Hubbard model on an anisotropic triangular lattice, a paradigmatic model capturing the interplay between strong correlations and magnetic frustration. However, the fate of frustrated magnetism in the presence of itinerant dopants remains unclear, as well as its connection to the doped phases of the square Hubbard model. Here we investigate the local spin order of a Hubbard model with controllable frustration and doping, using ultracold fermions in anisotropic optical lattices continuously tunable from a square to a triangular geometry. At half-filling and strong interactions U/t ≈ 9, we observe at the single-site level how frustration reduces the range of magnetic correlations and drives a transition from a collinear Néel antiferromagnet to a short-range correlated 120° spiral phase. Away from half-filling, the triangular limit shows enhanced antiferromagnetic correlations on the hole-doped side and a reversal to ferromagnetic correlations at particle dopings above 20%, hinting at the role of kinetic magnetism in frustrated systems. Furthermore, this work paves the way towards exploring possible chiral ordered or superconducting phases in triangular lattices and realizing t–t' square lattice Hubbard models that may be essential to describe superconductivity in cuprate materials.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2022

CRADA NO. NFE-22-00001 Between CNS and Northrop Grumman Systems Corporation

Consolidated Nuclear Security, LLC (CNS) and Northrop Grumman Systems Corporation (NGC) have established a partnership through the subject Cooperative Research and Development Agreement (CRADA) in three areas of collaboration: (1) Additive Manufacturing (AM) of manufactured carbon-carbon composite materials, (2) modeling and simulation, and (3) specialized chip and circuit board manufacturing for supply chain assurance.

36 MATERIALS SCIENCE↗

DOE OSTI2024

GeSn-on-GaAs with photoconductive carrier lifetime >400 ns: role of substrate orientation and atomistic simulation

Group IV GeSn quantum material finds application in electronics and silicon-compatible photonics. Synthesizing these materials with low defect density and high carrier lifetime is a potential challenge due to lattice mismatch induced defects and tin segregation at higher growth temperature. Recent advancements in the growth of these GeSn materials on Si, Ge, GaAs, and with substrate orientations, demonstrated different properties using epitaxial and chemical deposition methods. This article addresses the effect of GaAs substrate orientation and misorientation on the materials’ properties and carrier lifetimes in epitaxial Ge 0.94 Sn 0.06 layers. With starting GaAs substrates of (100)/2°, (100)/6°, (110) and (111)A orientations, Ge 0.94 Sn 0.06 epitaxial layers were grown with an intermediate Ge buffer layer by molecular beam epitaxy and analyzed by several analytical tools. X-ray analysis displayed good crystalline quality, and Raman spectroscopy measurements showed blue shifts in phonon wavenumber due to biaxial compressive strain in Ge 0.94 Sn 0.06 epilayers. Cross-sectional transmission electron microscopy analysis confirmed the defect-free heterointerface of Ge 0.94 Sn 0.06 /Ge/GaAs heterostructure. Minority carrier lifetimes of the unintentionally doped n-type Ge 0.94 Sn 0.06 epilayers displayed photoconductive carrier lifetimes of >400 ns on (100)/6°, 319 ns on (100)/2°, and 434 ns on (110) GaAs substrate at 1500 nm excitation wavelength. On the other hand, Ge 0.94 Sn 0.06 layer showed poor carrier lifetime on (111)A GaAs substrate. The observed differences in carrier lifetimes were correlated with the formation energy of the Ge on (100)/6° and (100)/2° GaAs heterointerface using Stillinger–Weber interatomic potential model-based atomistic simulation with different heterointerfacial bonding by Synopsys QuantumATK tool. Total energy computation of 6280-atom Ge/GaAs supercell on (100)/6° leads to lower formation energy than (100)/2°, making it more thermodynamically stable. Hence, the growth of the GeSn/III–V material system using misoriented (100) substrates that are more thermodynamically stable will enhance the performances of optoelectronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Effects of Strain Rate and Temperature on the Mechanical Properties of Simulated Silica Ionogels

Ionogels are hybrid materials formed by impregnating the pore space of a solid matrix with a conducting ionic liquid. By combining the properties of both component materials, ionogels can act as self-supporting electrolytes in Li batteries. In this study, molecular dynamics simulations are used to investigate the dependence of mechanical properties of silica ionogels on solid fraction, temperature, and pore width. Comparisons are made with corresponding aerogels. We find that the solid matrix fraction increases the moduli and strength of the ionogel. This varies nonlinearly with temperature and strain rate, according to the contribution of the viscous ionic liquid to resisting deformation. Additionally, owing to the temperature and strain sensitivity of the ionic liquid viscosity, the mechanical properties approach a linear mixing law at high temperature and low strain rates. The median pore width of the solid matrix plays a complex role, with its influence varying qualitatively with deformation mode. Narrower pores increase the relevant elastic modulus under shear and uniaxial compression but reduce the modulus obtained under uniaxial tension. Conversely, shear and tensile strength are increased by narrowing the pore width. All of these pore size effects become more pronounced as the silica fraction increases. Pore size effects, similar to the effects of temperature and strain rate, are linked to the ease of fluid redistribution within the pore space during deformation-induced changes in the geometry of the pores.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Viscometric flow of dense granular materials under controlled pressure and shear stress

This study examines the flow of dense granular materials under external shear stress and pressure using discrete element method simulations. In this method, the material is allowed to strain along all periodic directions and adapt its solid volume fraction in response to an imbalance between the internal state of stress and the external applied stress. By systematically varying the external shear stress and pressure, the steady rheological response is simulated for: (1) rate-independent quasi-static flow; and (2) rate-dependent inertial flow. The simulated flow is viscometric with non-negligible first and second normal stress differences. While both normal stress differences are negative in inertial flows, the first normal stress difference switches from negative to slightly positive, and second normal stress difference tends to zero in quasi-static flows. The first normal stress difference emerges from a lack of coaxiality between a second-rank contact fabric tensor and strain rate tensor in the flow plane, while the second normal stress difference is linked to an excess of contacts in the shear plane compared with the vorticity direction. A general rheological model of second order (in terms of strain rate tensor) is proposed to describe the two types of flow, and the model is calibrated for various values of interparticle friction from simulations on nearly monodisperse spheres. The model incorporates normal stress differences in both regimes of flow and provides a complete viscometric description of steady dense granular flows.

42 ENGINEERING↗

DOE OSTI2022

Thermal transport modeling of laser-irradiated spheres

Thermal transport of uniformly laser-irradiated spheres of various materials is investigated computationally. One-dimensional simulations of low- to mid-Z materials (Be, Al, and Cu) are performed to evaluate the impact of nonlocal electron transport on experimental observables under laser intensities of relevance to direct-drive inertial confinement fusion. Here we compare thermal transport models of different levels of fidelity: flux-limited Spitzer–Harm diffusion, the Schurtz–Nicolai–Busquet (SNB) reduced-order nonlocal model, and a Fokker–Planck description. Spitzer–Harm diffusion with different flux-limiter factors are compared with different implementations of the SNB model in the HYDRA radiation hydrodynamics code. Under the conditions of interest, the peak heat flux in the thermal front with the SNB model shows good agreement with Fokker–Planck calculations, with the largest errors below 10% at 10 15 W/cm 2 laser intensity. From HYDRA-SNB simulations, two experimentally relevant effects are observed from nonlocal heat transport when compared to flux-limited Spitzer–Harm modeling: coronal temperatures are cooler due to reduced heat fluxes in the expanding plasma and (for mid-Z materials) x-ray emissions are enhanced due to preheating in the dense plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Efficient mapping between void shapes and stress fields using Deep Convolutional Neural Networks with sparse data

Establishing fast and accurate structure-to-property relationships is an important component in the design and discovery of advanced materials. Physics-based simulation models like the finite element method (FEM) are often used to predict deformation, stress, and strain fields as a function of material microstructure in material and structural systems. Such models may be computationally expensive and time intensive if the underlying physics of the system is complex. This limits their application to solve inverse design problems and identify structures that maximize performance. In such scenarios, surrogate models are employed to make the forward mapping computationally efficient to evaluate. However, the high dimensionality of the input microstructure and the output field of interest often renders such surrogate models inefficient, especially when dealing with sparse data. Deep convolutional neural network (CNN) based surrogate models have shown great promise in handling such high-dimensional problems. In this paper, a single ellipsoidal void structure under a uniaxial tensile load represented by a linear elastic, high-dimensional and expensive-to-query, FEM model. We consider two deep CNN architectures, a modified convolutional autoencoder framework with a fully connected bottleneck and a UNet CNN, and compare their accuracy in predicting the von Mises stress field for any given input void shape in the FEM model. Additionally, a sensitivity analysis study is performed using the two approaches, where the variation in the prediction accuracy on unseen test data is studied through numerical experiments by varying the number of training samples from 20 to 100.

surrogate modeling; convolutional neural networks;↗