Engineering Papers⌕ Search

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.

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

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Simulation of Inelastic Neutron Scattering Spectra Directly from Molecular Dynamics Trajectories

Inelastic neutron scattering (INS) is a widely used technique to study atomic and molecular vibrations. With the increasing complexity of materials and thus the INS spectra, being able to simulate the spectra from various atomistic models becomes an essential step and also a major bottleneck for INS data analysis. The conventional approach using density functional theory and lattice dynamics often falls short when the materials of interest are complex (e.g., defective, disordered, heterogeneous, amorphous, large-scale), for which molecular dynamics driven by an interatomic force field is a more common approach. In this paper, we demonstrate a method to directly convert molecular dynamics trajectories into simulated INS spectra, including not only fundamental but also higher order excitations. In this work, the results are compared with data collected on various representative samples from different neutron spectrometers. This development will open great opportunities by providing the key tool to perform in-depth analysis of INS data and to validate and optimize computer models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

The influence of counterion structure identity on conductivity, dynamical correlations, and ion transport mechanisms in polymerized ionic liquids

We used equilibrium and non-equilibrium atomistic simulations to probe the influence of anion chemistry on the true conductivity, dynamical correlations, and ion transport mechanisms in polymeric ionic liquids. An inverse correlation was found between anion self-diffusivities, ionic mobilities, and the anion size for spherical anions. While some larger asymmetric anions had higher diffusivities than smaller spherical anions, their diffusivities and mobilities did not exhibit a direct correlation to the anion volumes. Here, the conductivity and anion dynamical correlations also followed the same trends as displayed by the diffusivity and mobility of anions. All the systems we examined displayed positively correlated motion among anions, suggesting a contribution that enhances the conductivity beyond the ideal Nernst–Einstein value. Analysis of ion transport mechanisms demonstrated very similar hopping characteristics among the spherical anions despite differences in their sizes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigating structure and dynamics of unentangled poly(dimethyl- co -diphenyl)siloxane via molecular dynamics simulation

Polysiloxane is one of the most important polymeric materials in technological use. Polydimethylsiloxane displays glass-like mechanical properties at low temperatures. Incorporation of phenyl siloxane, via copolymerization for example, improves not only the low-temperature elasticity but also enhances its performance over a wide range of temperatures. Copolymerization with the phenyl component can significantly change the microscopic properties of polysiloxanes, such as chain dynamics and relaxation. However, despite much work in the literature, the influence of such changes is still not clearly understood. In this work, we systematically study the structure and dynamics of random poly(dimethyl- co -diphenyl)siloxane via atomistic molecular dynamics simulations. As the molar ratio Φ of the diphenyl component increases, we find that the size of the linear copolymer chain expands. At the same time, the chain-diffusivity slows down by over an order of magnitudes. In conclusion, the reduced diffusivity appears to be a result of a complex interplay between the structural and dynamic changes induced by phenyl substitution.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Triple junction solute segregation in Al-based polycrystals

Solute segregation is a crucial means of stabilizing nanostructured alloys, and at very small grain sizes, this requires consideration of triple junctions (TJs), which attain a meaningful volume fraction and thus become relevant for bulk material behavior. Here, in this study, the solute segregation spectra for grain boundary (GB) and TJ sites are calculated for a large number of dilute Al-based binary systems with available interatomic potentials. A defect-identifying algorithm is applied to quantify the average GB thickness and classify the intergranular site spectra into GB and junction subspectra. The algorithm is also applied to a hybrid electronic-level database for GB segregation of various solutes in Al, yielding polycrystalline TJ solute segregation spectra from first principles. The results suggest that TJ segregation is alloy or interatomic potential dependent and can exhibit either boundary or junction preference. With these spectra as inputs, the spectral GB segregation model gives quantitative predictions of segregation as a function of grain size, temperature, and total solute concentration, suitable for alloy screening and design.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamic data-driven multiscale modeling for predicting the degradation of a 316L stainless steel nuclear cladding material

Here, we have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Permeability simulation of kidney-bean shaped carbon fibers

Cost saving efforts in carbon fiber manufacturing may result in irregularly shaped cross-sections, often as a kidney-bean shape. Here, in this research, flow simulations were performed to determine transverse permeability of these kidney-bean shaped fibers. Three different unit-cell packing arrangements for the kidney-bean shaped fibers and unit-cell orientations of 0°, + 90°, and – 90° were considered over the fiber volume fraction range 0.30 – 0.75. The kidney-bean shaped fibers had two distinct hexagonally packed unit-cell geometries because of the fiber’s asymmetry with different packing efficiency. The maximum fiber packing fraction for kidney-bean shaped fibers was 0.923, which was larger than the 0.907 for circular fibers. When compared to circular fibers with hexagonal packing, the kidney-bean shaped fibers had up to a 74% reduction in permeability. The ±90° orientations for KB-fibers with hexagonal packing had an infiltration time up to 3.86 longer than the hexagonal packed circular fibers for the same volume fraction.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Impact of boundary layer simulation on predicting radioactive pollutant dispersion: a case study for HANARO Research Reactor using the WRF-MMIF-CALPUFF modeling system

Wind plays an important role in cases of unexpected radioactive pollutant dispersion, deciding distribution and concentration of the leaked substance. The accurate prediction of wind has been challenging in numerical weather prediction models, especially near the surface because of the complex interaction between turbulent flow and topographic effect. As such, in this study, we investigated the characteristics of atmospheric dispersion of radioactive material (i.e. 137 Cs) according to the simulated boundary layer around the HANARO research nuclear reactor in Korea using the Weather Research and Forecasting (WRF)-Mesoscale Model Interface (MMIF)-California Puff (CALPUFF) model system. We examined the impacts of orographic drag on wind field, stability calculation methods, and planetary boundary layer parameterizations on the dispersion of radioactive material under a radioactive leaking scenario. We found that inclusion of the orographic drag effect in the WRF model improved the wind prediction most significantly over the complex terrain area, leading the model system to estimate the radioactive concentration near the reactor more conservatively. We also emphasized the importance of the stability calculation method and employing the skillful boundary layer parameterization to ensure more accurate low atmospheric conditions, in order to simulate more feasible spatial distribution of the radioactive dispersion in leaking scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Atomistic thermodynamics and kinetics of dicalcium silicate dissolution

Low-temperature cement manufacturing has garnered academic and industrial attention for its low environmental footprints. However, the sluggish hydration kinetics of the resultant cement affects their early-age strength development. This motivates fundamental studies to unravel the mechanistic picture of the dissolution process and discover science-informed pathways to accelerate hydration. Standard atomistic simulations seldomly exceed a microsecond making them impractical to study slow dissolution processes. We are using rare event sampling techniques, we provide the mechanistic picture of Ca 2+ ion dissolution from a kink site on the dicalcium silicate surface. The Ca 2+ ion dissolution is comprised of two sequential stages: breaking restraints from the kink sites to form a ledge adatom and detaching from the ledge/terrace adatom sites into the solution. The first and second stages feature free energy barriers of ~63 kJ/mol and ~ 29 kJ/mol respectively, making the first stage the rate-limiting step of the entire Ca 2+ dissolution kinetics. Using the reactive flux method, the rate and equilibrium constants for each reaction step are calculated, which yield the Ca 2+ ion activity of ~1.03 × 10 -5 . The diffusion calculations indicate that the surface effects lower the self-diffusion coefficient of Ca 2+ ions at the solid-water interface.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Exponential Escape Rate of Filamentary Incubation in Mott Spiking Neurons

Mott materials such as vanadium oxides, when subject to a strong applied voltage, present an inhomogeneous insulator-to-metal transition with formation of metallic filaments within the insulating bulk. This property is enabling the development of compact and power-efficient neuromorphic devices known as Mott neurons. However, the nature of the transition has not been fully understood yet, as it may be attributed to different effects, including Joule self-heating and hot-carrier injection. Moreover, the experimental determination of the threshold voltage needed to induce the transition has proven to be challenging, as the transition becomes increasingly unpredictable when the threshold is approached. The physical understanding of these issues would not only deepen our understanding of Mott insulators, but would also be an important step toward the realization of neuromorphic devices based on such materials. In this work we use numerical simulations based on the Mott resistor network model to study the nature of the filament incubation and formation process. Here, we show that both electronic and thermal effects, in the form of current density focusing and Joule self-heating, respectively, contribute to the filamentary incubation and growth. Remarkably, we find that the percolation of the metallic filaments near the threshold is intrinsically stochastic, qualitatively similar to the familiar Arrhenius activated behavior and to the stochastic firing of biological neurons. More precisely, we characterize the filament percolation as a Poisson point process, which has the same probability distribution as mathematical models of neuronal firing with an exponential escape rate. Finally, we support the numerical simulation results by performing experiments in VO 2 that are in agreement with the exponential escape rate behavior. Thus, we establish a functionality of Mott insulators that opens a path toward implementing neuromorphic hardware with quantum materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Strain-driven oxygen vacancy ordering in LaNiO 3 thin films revealed by integrated differential phase contrast imaging in scanning transmission electron microscopy

Rare-earth nickelates, such as LaNiO 3 (LNO), exhibit complex electronic properties, with ordered oxygen vacancies (OOV) influencing conductivity and magnetic behavior. We investigate the structural stability of strain-induced OOV phases in LNO thin films grown on SrTiO 3 substrates and the impact of Ruddlesden–Popper (RP) faults. Using high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) and integrated differential phase contrast (iDPC) STEM imaging, we conducted atomic-scale structural and compositional analyses of OOV. Geometric phase analysis (GPA) was employed to measure the strain in fault-free and RP fault regions, while density functional theory (DFT) calculations explored different OOV arrangements in the LNO phase. Simulated iDPC-STEM imaging of energy-stabilized structures was performed to correlate with experimental results. Here, our findings reveal superstructure modulation in the chemical composition and atomic-scale lattice structure in LNO, primarily due to the formation of the OOV in Ni–O layers of the LaNiO 2.5 phase. The out-of-plane compressive strain of about 2% stabilizes this phase, reducing the strain, diminishing OOV, and transforming them into LNO.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A cycle-jump acceleration method for the crystal plasticity simulation of high cycle fatigue of the metallic microstructure

For metallic materials, high-cycle fatigue life is sensitive to underlying microstructure features including secondary phases, textures, grains morphology, etc. The traditional, data-based safe-life approaches for modeling fatigue don’t explicitly consider the microstructure and can’t guide study in microstructure modification for improved fatigue property. Crystal plasticity-based simulation provides increased model fidelity at the expense of immense computation time, making it inapplicable for high cycle fatigue. In this work, an acceleration method based on cycle-jump approach (Lesne and Savalle, 1989) was developed for microstructure-based high-cycle FE simulation using crystal plasticity constitutive-model. We report this method demonstrated high efficiency in benchmark tests of various conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Grain incompatibility determines the local structure of amorphous grain boundary complexions

Amorphous grain boundary complexions lack long-range crystalline order but are not featureless, as distinct gradients in structural short-range order have been reported through their thickness. In this work, we test the hypothesis that the distribution of short-range order is determined by the confining crystals using atomistic simulations of both Cu-Zr bicrystals and a random polycrystal. Voronoi polyhedra with structures similar to that of perfect face-centered cubic serve as signatures of high structural order and are only found at the amorphous-crystalline transition regions. The density of the ordered structural motifs within a specific amorphous-crystalline transition region is found to not be directly determined by the orientation and symmetry of the grain which touches it, but rather by the incompatibility between the two confining grains. Further, ordered polyhedra density is found to be inversely related to grain incompatibility, meaning that large incompatibilities between the confining crystals lead to less order in the amorphous-crystalline transition region. The finding that the entire grain-film-grain system must be considered to understand local structure unequivocally demonstrates that amorphous complexions are not simply a collection of independent phases which happen to nucleate at a grain boundary. Rather, an amorphous grain boundary complexion is a single entity that finds a local equilibrium configuration.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Excess carrier concentration in silicon devices and wafers: How bulk properties are expected to accelerate light and elevated temperature degradation

Light and elevated temperature induced degradation (LeTID) is accelerated nearly linearly by the presence of excess carriers. It is therefore important to understand how excess carrier concentration (Δn) changes as a function of exposure conditions, materials properties, and sample structure. We simulate Δn as a function of wafer thickness and bulk minority carrier lifetime (τ) in solar cells and wafers using SCAPS and Quokka3. We also derive closed-form analytic expressions. For wafers, there is a near-linear relationship between Δn and τ or thickness, whereas for solar cells, Δn in the bulk may become limited by rear surface recombination. Thus, LeTID may progress more quickly in wafers than in cells, with a stronger dependence on τ. When comparing experiments, observed degradation rates must be corrected between samples or conditions to account for differences in Δn. This study demonstrates three tools to estimate the magnitude of such corrections, which can aid in the quantitative interpretation of LeTID data and performance predictions. Finally, while each tool yields similar results, there are advantages to each approach that must be weighed in terms of simplicity of inputs versus sophistication of treatment. Incomplete specification of back contact characteristics in commercial products is identified as an important contributor to uncertainty in expected LeTID rates.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A finite viscoelastic constitutive model for low to high strain rate response of elastomers with application of strain rate-induced glass transition

Amorphous elastomers exhibit significant rate-stiffening and unique viscous flow characteristics across a wide range of strain rates, often undergoing glass transition above a strain rate threshold. We have developed a thermodynamically-consistent and micromechanically-inspired constitutive model for soft elastomeric materials to capture the rate-dependent stress-strain behavior and hysteresis when subjected to low to high strain rates. Here, our proposed constitutive model encapsulates the viscous flow of materials through molecular motion at low strain rates and local rearrangement and alignment of the molecules trying to overcome the intermolecular resistance at high strain rates, essentially covering the glass transition. We applied our constitutive model to uniaxial compression experiments performed at low and high strain rates for polyborosiloxane (PBS) to identify the material parameters, and subsequently, performed numerical simulations of single and multi-cycle compression, stress relaxation, and small amplitude oscillatory tension-compression. Our analyses indicate that the model predicts higher total energy dissipation with increasing strain rate; however, dissipation associated with molecular relaxation decreases (forming a cusp) because, beyond a crossover strain rate, intermolecular rearrangement and alignment become dominant, which is consistent with the onset of the glass transition. For cyclic loading-unloading, we observed that dissipation over a cycle remains constant at low strain rates but decreases non-monotonically at high strain rates before becoming constant, with the peak stress over the cycle becoming higher, which can be interpreted as more loading being carried elastically by the polymer network as the intermolecular rearrangement process occurs. Additionally, our model was able to predict the qualitative nature of the storage modulus and loss modulus in the limit of small strain over a wide range of frequency sweeps.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Understanding Asteroid 16 Psyche’s composition through 3D impact crater modeling

Asteroid 16 Psyche is the largest M-type (metallic) Main Belt Asteroid (MBA). Radar albedo data indicate Psyche’s surface is rich in metallic content, but estimates for Psyche’s bulk structure vary widely, with bulk density estimates ranging from 1.4 ± 0.3 g/cm 3 to 4.5 ± 1.4 g/cm 3 and porosity estimates ranging from 30% to 70%. A forthcoming NASA discovery mission will visit Psyche. Narrowing the range of estimated compositions of Psyche prior to mission launch can help to determine the necessary tools for asteroid analysis. Psyche has two large impact structures in its Southern hemisphere, with estimated diameters from 50 km to 70 km and estimated depths up to 6.4 km. In this work, we present results from 2D and 3D hydrocode simulations of the largest of these impact structures. Through quantitative comparison of the simulated crater dimensions with measured values, our models suggest that Psyche is largely composed of porous, metallic material. In addition, to match the crater dimensions, our simulations indicate that the angle of impact was significantly altered from vertical by 45–60 degrees.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Integrated simulation, machine learning, and experimental approach to characterizing fracture instability in indentation pillar-splitting of materials

Measuring fracture toughness of materials at small scales remains challenging due to limited experimental testing configurations. A recently developed indentation pillar-splitting method has shown promise of improved flexibility in fracture toughness measurements at the microscale, partly due to the occurrence of an unusual fracture instability, i.e., a transition from stable to unstable crack propagation. In spite of growing interest in this method, the underlying mechanism of this phenomenon is yet to be elucidated. Furthermore, we provide a comprehensive description of fracture instability in indentation pillar-splitting by combining in situ experiments with high-fidelity simulations based on cohesive zone and J-integral methods. In addition, a machine-learning-based solution for predicting the critical indentation load of fracture instability is established through Gaussian processes regression for broad use of this method by the community.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An ab initio molecular dynamics investigation of the thermophysical properties of molten NaCl-MgCl 2

Molten salts have many applications in the nuclear and solar energy industries for thermal storage and heat transfer applications. However, there is a knowledge gap in molten salt thermophysical properties which hinders the technical readiness level of molten salt applications, especially in the nuclear industry. A common method of investigating new materials is through ab initio Molecular Dynamics (AIMD) simulations which is an effective tool to investigate structural and thermophysical properties at realistic temperatures. NaCl-MgCl 2 is an inexpensive salt that is a good candidate for use as a heat transfer medium in solar power applications or in the secondary loop of a nuclear reactor. In this article, the thermophysical properties of NaCl-MgCl 2 are computed via AIMD calculations to supplement the limited experimental studies in the literature. Here a wide range of compositions and temperatures for the pseudo-binary NaCl-MgCl 2 were used to calculate the density, heat capacity, compressibility, enthalpy of mixing, and volumetric thermal expansion coefficient. AIMD is shown to accurately model the densities of molten NaCl-MgCl 2 as there is good agreement with the available literature. This work observed a transition to a monotonic increase of the density with respect to MgCl 2 composition occurring above 1100 K. The heat capacity values increase uniformly with respect to concentration of MgCl 2 at a rate of 2.85 J/mol-K per 10 mol% of MgCl 2 . Select thermophysical properties are fit to a Redlich-Kister expansion for utilization in multiphysics simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Understanding the surface wave characteristics using 2D particle-in-cell simulation and deep neural network

Here, the characteristics of the surface waves along the interface between a plasma and a dielectric material have been investigated using kinetic particle-in-cell simulations. A microwave source of GHz frequency has been used to trigger the surface wave in the system. The outcome indicates that the surface wave gets excited along the interface of plasma and the dielectric tube and appears as light and dark patterns in the electric field profiles. The dependency of radiation pressure on the dielectric permittivity and supplied input frequency has been investigated. Further, we assessed the capabilities of neural networks to predict the radiation pressure for a given system. The proposed deep neural network model is aimed at developing accurate and efficient data-driven plasma surface wave devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗