Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural 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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Summary of Structural Material Modeling Development for the NEAMS Program in Fiscal Year 2020

This report summarizes work performed during Fiscal Year (FY) 2020 at Idaho National Laboratory (INL) for the U.S. Department of Energy?s Nuclear Engineering Advanced Modeling and Simulation (NEAMS) program for the Structural Materials and Chemistry Technical Area in the work package entitled "MS- 20IN050104 - Structural Materials - INL." The Structural Materials and Chemistry Technical Area is a relatively new component of the NEAMS program, and is currently focusing on developing simulation capabilities to support the deployment of nuclear energy in the areas of molten salt reactor chemistry, light water reactor (LWR) structural material degradation, and structural material behavior for advanced reactor applications. INL performed work for to advance capabilities for simulation of structural material behavior in both LWR and advanced reactor applications in the work described here.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Charge transfer effect on local lattice distortion in a HfNbTiZr high entropy alloy

We report it is often assumed that atoms are hard spheres in the estimation of local lattice distortion (LLD) in high-entropy alloys (HEAs). However, our study demonstrates that the hard sphere model misses the key effect, charge transfer among atoms with different electronegativities, in the understanding of the severely distorted local structure in the body-centered cubic (BCC) HfNbTiZr HEA. Through the characterization and simulations of the local structure of the HfNbTiZr HEA, we found that the LLD regions in HfNbTiZr have larger average size compared with face-centered cubic (FCC) HEAs, but the magnitude of LLD significantly varies from one atomic site to another because of electron transfer from larger atoms, Hf and Zr, to smaller ones, Nb and Ti. Our finding may form the basis for the design of high performance HEAs through tuning their site-to-site chemical heterogeneity.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An Integrated Simulation of Multiple-Pass U-10Mo Alloy Hot Rolling and Static Recrystallization

To achieve a desired microstructure and minimize the thickness variation in rolled foils, researchers must understand the effects of foil fabrication process variables on microstructure evolution. We developed an integrated simulation of deformation and recrystallization that employs the finite element method (FEM) and the kinetic Monte Carlo (KMC) Potts model, respectively, to investigate microstructure evolution during multiple-pass hot rolling and heat treatment in polycrystalline U-10Mo fuel. Scanning electron microscopy and electron backscatter diffraction images of microstructures were directly used as input in FEM calculation of deformation, and the calculated strains were used to determine the driving force of nucleation and growth of recrystallized grains in the Potts model. Grain structures predicted by the Potts model were used to update the grain structure and material properties for FEM. Simulation alternated between FEM and the Potts model to simulate grain structure evolution during multiple rolling and heat treatments. The initial model parameters were determined by benchmarking the recrystallization kinetics against experimental data. Then, the model was applied to predict the grain structure evolution. Results showed that our model can capture the coupling between deformation and recrystallization and can quantitatively reproduce the observed U-10Mo recrystallization and grain growth kinetics. The simulation results demonstrated that the developed model can predict U-10Mo grain structures as a function of initial microstructure and foil fabrication parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Time-resolved photoelectron spectroscopy via trajectory surface hopping

Time-resolved photoelectron spectroscopy is a powerful pump-probe technique which can probe nonadiabatic dynamics in molecules. Interpretation of the experimental signals however requires input from theoretical simulations. Advances in electronic structure theory, nonadiabatic dynamics, and theory to calculate the ionization yields, have enabled accurate simulation of time-resolved photoelectron spectra leading to successful applications of the technique. Here, we review the basic theory and steps involved in calculating time-resolved photoelectron spectra, and highlight successful applications.

Chakraborty, Pratip↗
DOE OSTI2022

Finding Electronic Structure Machine Learning Surrogates without Training

A myriad of phenomena in materials science and chemistry rely on quantum-level simulations of the electronic structure in matter. While moving to larger length and time scales has been a pressing issue for decades, such large-scale electronic structure calculations are still challenging despite modern software approaches and advances in high-performance computing. The silver lining in this regard is the use of machine learning to accelerate electronic structure calculations – this line of research has recently gained growing attention. The grand challenge therein is finding a suitable machine-learning model during a process called hyperparameter optimization. This, however, causes a massive computational overhead in addition to that of data generation. We accelerate the construction of machine-learning surrogate models by roughly two orders of magnitude by circumventing excessive training during the hyperparameter optimization phase. We demonstrate our workflow for Kohn-Sham density functional theory, the most popular computational method in materials science and chemistry.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phase-field simulation of domain size effect on dielectric and piezoelectric responses in K 0.5 Na 0.5 NbO 3 epitaxial thin films with superdomain structures

Size effects of mesoscale ferroelectric domains on the macroscopic dielectric and piezoelectric responses in domain-engineered bulk piezocrystals have been extensively studied for more than two decades. However, less is known about the domain size effects in ferroelectric epitaxial thin films, especially for films consisting of low-symmetry ferroelectric phases and exhibits hierarchical superdomain structures. Herein, using phase-field simulations, we systemically evaluate the effective out-of-plane dielectric and piezoelectric coefficients, κ$^{*}_{33}$ and d$^{*}_{33}$, as a function of the domain periods for two types of superdomain structures in ferroelectric K 0.5 Na 0.5 NbO 3 epitaxial thin films. In one type of the superdomain structures, we find that more than 70% increase of κ$^{*}_{33}$ and nearly 20% increase of d$^{*}_{33}$ can be achieved by tuning the domain period by a few tens of nanometers. Dissimilar behaviors are found when the domain period varies along different lateral directions, suggesting anisotropic domain size effects in thin films. By analyzing the local dielectric and piezoelectric responses from each domain variants and domain walls, we reveal that the domain size effect is governed by the variation of out-of-plane polarization inside the domains. Moreover, we also demonstrate enhanced domain size effects by modulating the misfit strains and temperature to approach the polymorphic phase boundaries, suggesting tunability of the size effect by doping and strain engineering. Furthermore, our results reveal multimodal domain size dependence of dielectric and piezoelectric responses in low-symmetry ferroelectric epitaxial thin films, implying that domain size engineering can be used to tune macroscopic properties of thin-film ferroelectrics, similar to their bulk counterparts.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Atomistic Characterization of Gramicidin Channel Formation

We investigated gramicidin A (gA) subunit dimerization in lipid bilayers using microsecond-long replica-exchange umbrella sampling simulations, millisecond-long unbiased molecular dynamics simulations, and machine learning. Our simulations led to a dimer structure that is indistinguishable from the experimentally determined gA channel structures, with the two gA subunits joined by six hydrogen bonds (6HB). The simulations also uncovered two additional dimer structures, with different gA–gA stacking orientations that were stabilized by four or two hydrogen bonds (4HB or 2HB). When examining the temporal evolution of the dimerization, we found that two bilayer-inserted gA subunits can form the 6HB dimer directly, with no discernible intermediate states, as well as through paths that involve the 2HB and 4HB dimers.

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

Investigations of irradiation effects in crystalline and amorphous SiC

The results of irradiation on 3C-silicon carbide (SiC) and amorphous SiC (a-SiC) are explored using both in situ transmission electron microscopy (TEM) and complementary molecular dynamics (MD) simulations. The single ion strikes identified in the in situ TEM irradiation experiments, utilizing a 1.7 MeV Au 3+ ion beam with nanosecond resolution, are contrasted to MD simulation results of the defect cascades produced by 10–100 keV Si primary knock-on atoms (PKAs). The MD simulations also investigated defect structures that could possibly be responsible for the observed strain fields produced by single ion strikes in the TEM ion beam irradiation experiments. Both MD simulations and in situ TEM experiments show evidence of radiation damage in 3C-SiC but none in a-SiC. Selected area electron diffraction patterns, due to the results of MD simulations and in situ TEM irradiation experiments, show no evidence of structural changes in either 3C-SiC or a-SiC.

3C-SiC↗
DOE OSTI2024

Structural and Spectroscopic Properties of Butanediol-Modified Boehmite Materials

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, in this study, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O–H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

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

Comparative CHARMM and AMOEBA Simulations of Lanthanide Hydration Energetics and Experimental Aqueous-Solution Structures

The accurate understanding of metal ion hydration in solutions is a prerequisite for predicting stability, reactivity, and solubility. Herein, additive CHARMM force field parameters were developed to enable molecular dynamics simulations of lanthanide (Ln) speciation in water. Quantitatively similar to the much more resource-intensive polarizable AMOEBA potential, the CHARMM simulations reproduce the experimental hydration free energies and correlations in the first shell (Ln-oxygen distance and hydration number). Comparisons of difference pair-distribution functions obtained from the two simulation approaches with those from high-energy X-ray scattering experiments reveal good agreement of first coordination sphere correlations for the Lu3+ ion (CHARMM only), but further improvement to both approaches is required to reproduce the broad, non-Gaussian distribution seen from the La3+ experiment. Second-coordination sphere comparisons demonstrate the importance of explicitly including an anion in the simulation. This work describes the usefulness of less resource-intensive additive potentials in some complex chemical systems such as solution environments where multiple interactions have similar energetics. In addition, 3-dimensional descriptions of the La3+ and Lu3+ coordination geometries are extracted from the CHARMM simulations and generally discussed in terms of potential improvements to solute-structure modeling within solution environments.

Qiao, Baofu↗
DOE OSTI2023

Learning local equivariant representations for large-scale atomistic dynamics

Abstract A simultaneously accurate and computationally efficient parametrization of the potential energy surface of molecules and materials is a long-standing goal in the natural sciences. While atom-centered message passing neural networks (MPNNs) have shown remarkable accuracy, their information propagation has limited the accessible length-scales. Local methods, conversely, scale to large simulations but have suffered from inferior accuracy. This work introduces Allegro, a strictly local equivariant deep neural network interatomic potential architecture that simultaneously exhibits excellent accuracy and scalability. Allegro represents a many-body potential using iterated tensor products of learned equivariant representations without atom-centered message passing. Allegro obtains improvements over state-of-the-art methods on QM9 and revMD17. A single tensor product layer outperforms existing deep MPNNs and transformers on QM9. Furthermore, Allegro displays remarkable generalization to out-of-distribution data. Molecular simulations using Allegro recover structural and kinetic properties of an amorphous electrolyte in excellent agreement with ab-initio simulations. Finally, we demonstrate parallelization with a simulation of 100 million atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Aeroelastic Modeling and Full-Scale Loads Measurements for Investigation of Single-Axis PV Tracker Wind-Driven Dynamic Instabilities

While wind tunnel testing and proprietary industry modeling tools have been used for years to design and develop PV tracker systems, wind induced dynamic failures are becoming more prevalent and high visibility. Designers and manufactures have reacted by developing new systems or add-on products to improve system dynamics adding to overall system costs. To better understand the physics and aeroelastic behavior that leads to dynamic instabilities NREL researches have taken up a first of a kind study to both develop open source aeroelastic modeling tools and measure tracker loads on a single axis full-scale tracker system at a high wind site. These modeling tools can simulate the fluid-structure interaction driving torsional instabilities under a wide range of turbulent inflow conditions and stow angles. The results of these simulations reveal the flow features associated with panel rotation and the induced loads in the structure. The experimental load measurements can be used to validate models and provide a quantifiable understanding of tracker dynamics, critical loads paths, identify instability markers, inform resilient design, and suggest the most favorable stow approach.

14 SOLAR ENERGY↗
DOE OSTI2018

Effect of atomic order/disorder on vacancy clustering in concentrated NiFe alloys

Using molecular dynamics simulations, we elucidate in this paper the effect of atomic structure on vacancy clustering in ordered (L1 0 ) and random NiFe. Based on our simulations, we predict the vacancy evolution to be in complete contrast between the two systems. While large vacancy clusters, i.e., stacking fault tetrahedra (SFT) are formed in the random structure, no clustering is observed in the ordered-L1 0 structure. Similar simulations are performed on L1 0 -CuAu and L1 0 -TiAl to understand whether SFT formation is generic in L1 0 structures, or is specific to NiFe. Both materials show SFT formation, thereby highlighting specific defect energetics in L1 0 NiFe that lead to the lack of vacancy clustering. We elucidate that L1 0 -NiFe has unique thermodynamic and kinetic defect energetics, i.e., antisite energies, vacancy sublattice preference, and directional migration energy barriers that collectively lead to the lack of vacancy clustering. Finally, understanding such defect energetics could open avenues to prevent defect clustering in the vision towards development of radiation-tolerant concentrated alloys for nuclear reactor applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Quantitative Structure–Property Relationship (QSPR) Analysis of ZrO 2 -Containing Soda-Lime Borosilicate Glasses

Quantitative structure–property relationship (QSPR) analysis is a promising approach to correlate structural features with properties of glass materials that lack long-range order and usually have complex structures. By using carefully chosen descriptors based on structural models generated from molecular dynamics (MD) simulations, correlations with properties and insights on glass behaviors can be obtained. Zirconia can significantly alter glass properties including chemical durability, even in a small amount, and hence plays an important role in vitrification of nuclear waste where long-term chemical durability is desired. In this study, borosilicate glasses with the composition of xZrO 2 –(61 – x)SiO 2 –17B 2 O 3 –18Na 2 O–4CaO with x = 0, 1, 2, 4, 6, and 8 were simulated using classical MD simulations with the recently developed composition-dependent potentials. Short-range (e.g., bond distance and coordination numbers) and medium-range (e.g., Qn distribution, network connectivity, and ring-size distribution) structural features altered by ZrO 2 were obtained and analyzed. The use of a descriptor (Fnet descriptor) that combines short-range structural characteristics, from MD simulations, and the cation-oxygen single bond strength was found to provide excellent linear correlations with the density and initial dissolution rate of these glasses. Furthermore, the results show that by combining MD simulations and QSPR analysis the composition and structural effect on the properties of complex multicomponent glasses can be elucidated, thus suggesting that this is a promising approach for future glass research and new composition design.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Structural Anomalies and Electronic Properties of an Ionic Liquid under Nanoscale Confinement

Ionic liquids (ILs) promise far greater electrochemical performance compared to aqueous systems, yet key physicochemical properties governing their assembly at interfaces within commonly used graphitic nanopores remain poorly understood. Here, we combine synchrotron X-ray scattering with first-principles molecular dynamics simulations to unravel key structural characteristics of 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([TFSI] – ) ionic liquids confined in carbon slit pores. X-ray scattering reveals selective pore filling due to size exclusion, while filled pores exhibit disruption in the IL intermolecular structure, the extent of which increases for narrower slit pores. First-principles simulations corroborate this finding and quantitatively describe how perturbations in the local IL structure, particularly the hydrogen-bond network, depend strongly on the degree of confinement. Despite significant deviations in structure under confinement, electrochemical stability remains intact, which is important for energy storage based on nanoporous carbon electrodes (e.g., supercapacitors).

36 MATERIALS SCIENCE↗
DOE OSTI2022

An incompressible Eulerian method for fluid–structure interaction with mixed soft and rigid solids

Here, we present a general simulation approach for incompressible fluid–structure interactions in a fully Eulerian framework using the reference map technique. The approach is suitable for modeling one or more rigid or finitely deformable objects or soft objects with rigid components interacting with the fluid and with each other. It is also extended to control the kinematics of structures in fluids. The model is based on our previous Eulerian fluid–soft solver [Rycroft et al., “Reference map technique for incompressible fluid–structure interaction,” J. Fluid Mech. 898, A9 (2020)] and generalized to rigid structures by constraining the deformation-rate tensor in a projection framework. Several numerical examples are presented to illustrate the capability of the method.

42 ENGINEERING↗
DOE OSTI2023

The Sluggish Diffusion of Cations in CeO 2 Probed through Molecular Dynamics and Metadynamics Simulations

Cation diffusion in fluorite‐structured CeO 2 , though far slower than anion diffusion, is an important, high‐temperature process because it governs diverse fabrication and degradation phenomena. Herein, cation diffusion is studied by means of classical molecular dynamics and metadynamics simulations. Three different mechanisms are examined: migration involving an isolated cerium vacancy, migration involving a cerium vacancy in a defect associate with an oxygen vacancy, and migration involving a cation divacancy. For each mechanism, defect diffusion coefficients are calculated as a function of temperature, from which the respective activation enthalpy of defect migration is obtained. Through comparisons with experimental cation diffusion data (specifically, of the absolute magnitude of the cation diffusivity as well as its activation enthalpy), it is concluded that cation diffusion takes place predominantly neither by isolated vacancies nor by cation vacancy–oxygen vacancy associates but by cation divacancies.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Exploration of Tertiary Structure in Sequence-Defined Polymers Using Molecular Dynamics Simulations

Peptoids are a class of sequence-defined biomimetic polymers with peptide-like backbones and side chains located on backbone nitrogens rather than alpha carbons. These materials demonstrate a strong ability for precise control of single-chain structure, multiunit self-assembly, and macromolecular assembly through careful tuning of sequence due to the diversity of available side chains, although the driving forces behind these assemblies are often not understood. Prior experimental work has shown that linked 15mer peptoids can mimic the protein helical hairpin structure by leveraging the chirality-inducing nature of bulky side chains and hydrophobicity, but there are still gaps in our understanding of the relationship between sequence, stability, and particular secondary or tertiary structure. Here, we present a molecular dynamics (MD) study on the folding behavior of these polymers into hairpins, discussing the differences in structure from sequences with various characteristics in water and acetonitrile, and then compare the handedness preference of common helical motifs between solvents.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗