Search NASA⌕ Search

SEARCH · Search NASA

Results for “dynamics 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 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

High-resolution neutron time-of-flight measurements for light water at the Spallation Neutron Source (SNS), Oak Ridge National Laboratory

Series of light water inelastic neutron scattering experiments have been made at the Oak Ridge National Laboratory (ORNL), Spallation Neutron Source (SNS) covering temperatures ranging from 295 K to 600 K and pressures of 1 bar and 150 bar. The temperatures and pressures ranges correspond to that of pressurized light water reactors. The inelastic scattering measurements will help the development of light water thermal scattering kernels, also known as S (α,β) thermal scattering law (TSL), in a consistent fashion given the amount and the quality of the measured data. Light water thermal scattering evaluations available in existing nuclear data libraries have certain limitations and pitfalls. This paper introduces the state of the art of the light water thermal scattering cross-section data not only for room temperature but as well as for reactor operating temperatures, i.e. 550 - 600 K. During the past few years there has been a renewed interest in re-investigating the existing TSL models and utilize the recent experimental data or perform molecular dynamics simulations. It should be pointed out that no single TSL evaluation is based entirely on experimental data and one has to rely on TSL models or a combination of both. New TOF measurement of light water at the SNS, with a detailed description of the experimental setup, measurement conditions, and the associated foreseen results is presented in this paper. The analysis of the experimental data would help in validating the existing approach based on old experimental data or based on molecular dynamic simulations using classical water models, knowledge of which is very important to generate TSL libraries at reactor operating conditions.

Leal, Luiz C.↗
DOE OSTI2016

Structure of ionic liquids with cationic silicon-substitutions

Significantly lower viscosities result when a single alkyl carbon is replaced by a silicon atom on the side chain of an ionic liquid cation. To further explore this effect, we compare liquid structure factors measured using high-energy X-ray scattering and calculated using molecular dynamics simulations. In this work, four ionic liquids are studied that each has a common anion, bis(trifluoromethylsulfonyl)amide (NTf${^–_2}$). The four cations for this series of NTf${^–_2}$-anion ionic liquids are 1-methyl-3-trimethylsilylmethylimidazolium (Si-mim + ), 1-methyl-3-neopentylimidazolium (C-mim + ), 1-methyl-3-pentamethyldisiloxymethylimidazolium (SiOSi-mim + ), and 1-methyl-1- trimethylsilylmethylpyrrolidinium (Si-pyrr + ). To achieve quantitative agreement between the structure factors measured using high-energy X-ray scattering and molecular dynamics simulations, new transferable parameters for silicon were calibrated and added to the existing force fields.

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

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

Dynamic models for simulation of the 70-M antenna axis servos

Dynamic models for the various functional elements of the 70 m antenna axis servos are described. The model representing the digital position controller, the linear and nonlinear properties of the physical hardware, and the dynamics of the flexible antenna structure are encoded in six major function blocks. The general modular structure of the function blocks facilitates their adaptation to a variety of dynamic simulation studies. Model parameter values were calculated from component specifications and design data. A simulation using the models to predict limit cycle behavior produced results in excellent agreement with field test data from the DSS 14 70-m antenna.

Hill, R. E.↗
DOE OSTI2015

Propensity of bond exchange as a window into the mechanical properties of metallic glasses

We investigated the mechanical properties of Zr-Cu-Al bulk metallic glasses, by compression experiment and molecular dynamics simulations. From the simulation, we found that the large, solvent atom, Zr, has high propensity of bond exchange compared to those of the smaller solute atoms. The difference in bond exchange is consistent with the observed disparity in mechanical behaviors: Zr-rich metallic glass exhibits low elastic modulus and large plastic strain. Furthermore, X-ray photoelectron spectroscopy measurements suggest that the increased propensity in bond exchange is related to the softening of Zr bonds with increasing Zr content.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Online Static Load Model Estimation in Distribution Systems: Preprint

The paper focuses on parameter identification for time varying voltage dependent polynomial load models, namely the constant impedance (Z), constant current (I), and constant power (P), in the distribution systems. Typically, parameter identification has been done for aggregated loads at substations for the transmission planning studies. Given the higher voltage volatility faced by active distribution systems, this paper takes the load model estimation to the distribution networks. The paper assumes the presence of measurements at the distribution nodes from various sources, each at varying temporal resolutions. Least squares estimation (LSE) method with bounded variables is used to convert the measurements into ZIP parameter estimations for loads at various locations of the feeder. The results and discussions focus on some of the estimation issues faced in a distribution system due to varying voltage sensitivities of loads along a feeder, measurement resolution, impact of window size and sampling rate of measurement data, and the use of dynamic simulations to create synthetic measurement data in the absence of real ones.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Understanding the Influence of Chain Architecture on the Transport Quantities of Polymer Electrolytes with Covalently Bonded Anions

Here, we use a combination of experiments and coarse-grained molecular dynamics simulations to elucidate the structure–property relationships in polymer electrolytes obtained by the copolymerization of poly(vinyl ethylene carbonate─lithium styrene bis(trifluoromethanesulfonyl)imide) or p(VEC-LiSTFSI). Experiments show that the conductivity reduces with increasing anion (i.e., STFSI) fraction on the chain, and the cation transference number (t + ) is found to be dependent on the anion fraction. Furthermore, a significant fraction of unpolymerized VEC monomers are observed. Since it is inherently difficult to experimentally control the chain architecture and the amount of unpolymerized VEC in these systems, we perform coarse-grained molecular dynamics simulations on model polymer systems with different chain architectures to mimic the plausible experimental systems. Specifically, we look at the differences in transference numbers arising from (i) a random copolymer of VEC and STFSI monomers; (ii) a blend of VEC-STFSI copolymer with VEC monomers; and (iii) a ternary blend of the VEC homopolymer, STFSI homopolymers, and VEC monomers. The ternary blend model demonstrates the closest resemblance with the experimental transference numbers and diffusivities. The lithium diffusivity obtained from the coarse-grained models with VEC monomers (plasticizers) is about 1.5 times that of the model without VEC monomers, showing that the plasticizing effect of VEC monomers is modest. We rationalize the experimental observations based on aggregate and cluster analyses obtained from molecular simulations. This work reveals that polymer electrolyte chain architecture and plasticizers can critically influence the transport properties, and these parameters should be considered when designing single ion conducting polymeric electrolytes.

cluster distribution↗
DOE OSTI2023

Impact of Soret effect on hydrogen and helium retention in PFC tungsten under ELM-like conditions

Abstract In our previous work, we have demonstrated using nonequilibrium molecular-dynamics simulations that the fluxes of helium and self-interstitial atoms in the presence of a thermal gradient in tungsten are directed opposite to the heat flux, indicating that species transport is governed by a Soret effect, namely, thermal-gradient-driven diffusion, characterized by a negative heat of transport that drives species transport uphill, i.e. from the cooler to the hot regions of the tungsten sample. In this work, the findings of our thermal and species transport analysis have been implemented in our cluster-dynamics code, Xolotl, which has been used to compute temperature and species profiles over spatiotemporal scales representative of plasma-facing component (PFC) tungsten under typical reactor operating conditions, including extreme heat loads at the plasma-facing surface characteristic of plasma instabilities that induce edge localized modes (ELMs). We demonstrate that the steady-state species profiles, when properly accounting for the Soret effect, vary significantly from those where temperature-gradient-driven transport is not accounted for and discuss the implications of such a Soret effect on the response to plasma exposure of plasma-facing tungsten. Although our cluster-dynamics simulations do not yet include self-clustering of helium or hydrogen blister formation, our simulation results show that the Soret effect substantially reduces helium and hydrogenic species retention inside PFC tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2022

E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials

Abstract This work presents Neural Equivariant Interatomic Potentials (NequIP), an E(3)-equivariant neural network approach for learning interatomic potentials from ab-initio calculations for molecular dynamics simulations. While most contemporary symmetry-aware models use invariant convolutions and only act on scalars, NequIP employs E(3)-equivariant convolutions for interactions of geometric tensors, resulting in a more information-rich and faithful representation of atomic environments. The method achieves state-of-the-art accuracy on a challenging and diverse set of molecules and materials while exhibiting remarkable data efficiency. NequIP outperforms existing models with up to three orders of magnitude fewer training data, challenging the widely held belief that deep neural networks require massive training sets. The high data efficiency of the method allows for the construction of accurate potentials using high-order quantum chemical level of theory as reference and enables high-fidelity molecular dynamics simulations over long time scales.

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

The Debye-Waller factor of disordered solids: The case for amorphous silicon

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (α-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in α-Si networks. It has been observed that the MSDs associated with tetrahedrally bonded Si atoms are considerably lower than their dangling-bond counterparts—originating from isolated and vacancy-induced clustered defects—and those on the surface of microvoids, leading to an asymmetric non-Gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in α-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of α-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of α-Si at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Ab initio studies of the impact of the Debye-Waller factor on the structural and dynamical properties of amorphous semiconductors: The case of a-Si

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (a-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in a-Si networks. It has been observed that the MSDs associated with tetrahedrally-bonded Si atoms are considerably lower than their dangling-bond counterparts – originating from isolated and vacancy-induced clustered defects – and those on the surface of microvoids, leading to an asymmetric non-gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in a-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of a-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of a-Si at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2018

GhoT of the GhoT/GhoS toxin/antitoxin system damages lipid membranes by forming transient pores

Highlights: • Determined that toxin GhoT disrupts membranes by forming transient pores. • Determined that membrane charge affects the activity of GhoT. • Determined through molecular dynamic simulations why GhoT F38R is inactive. GhoT is a bacterial toxin of the type V toxin/antitoxin system that allows Escherichia coli to reduce its metabolism in response to oxidative and bile stress. GhoT functions by increasing membrane permeability and reducing both ATP levels and the proton motive force. However, how GhoT damages the inner membrane has not been elucidated. Here we investigated how GhoT damages membranes by studying its interaction with lipid bilayers and determined that GhoT does not cause macroscopic disruption of the lipid bilayer to increase membrane permeability to the dye carboxyfluorescein. Using circular dichroism, we found that GhoT forms an alpha helical structure in lipid bilayers that agrees with the structure predicted by the I-TASSER protein structure prediction program. The structure generated using I-TASSER was used to conduct coarse-grained molecular dynamics simulations, which indicate that GhoT damages the cell membrane, as a multimer, by forming transient transmembrane pores.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Dynamic Structure of Magnetized High Energy Density Plasmas (Final Report)

This report summarizes results of project DE-SC0022202 ``Dynamic Structure of Magnetized High Energy Density Plasmas,'' which was originally funded for the period 08/01/2021 - 7/31/2022, but was no-cost extended to 7/31/2023. The primary objectives of the work, as stated in the original proposal, were ``develop a theory to describe a novel state of magnetized high energy density plasmas, benchmark the theory using first-principles molecular dynamics simulations, and apply the results to develop a means to interpret x-ray scattering measurements in future experiments.'' These were the goals of the original project proposal, which was for a 3 year project. The project was ultimately funded for 1 year, and so the work completed focused on the goals described for only the first year of the proposed work. This included further developing a kinetic theory for strongly magnetized plasmas, applying it to compute fluid-scale (MHD) transport coefficients, and to test the predictions using molecular dynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A ReaxFF molecular dynamics study of molecular-level interactions during binder jetting 3D-printing

In the present work, we study one of the major additive manufacturing processes, i.e., the binder jetting printing (BJP) process, at the molecular level through atomistic-scale level representations of powders and binder solutions with chromium-oxide (Cr-oxide) nanoparticles and water-based diethylene glycol solutions, respectively. The results show that both diethylene glycol and water contribute to the bonding of Cr-oxide particles during the print and curing stages by forming a hydrogen bond network. Heating the system to the burn-out temperature results in the oxidation of diethylene glycol and the decomposition of the hydrogen bond network. Subsequently, Cr-oxide particles are partially sintered by forming Cr–O bonds. The final sintering facilitates further Cr–O bond formation. Additionally, the influence of the chemical composition of the binder solution is investigated by performing ReaxFF molecular dynamics simulations on two sets of systems, which control the number of water and diethylene glycol molecules, respectively. Our results demonstrate that adding both diethylene glycol and water to the binder solution can raise the number of “useful” hydrogen bonds, resulting in a higher breaking strength at the print and curing stages. During the burn-out and sintering stages, the influence of water on the breaking strength is not obvious. In contrast, an optimal quantity of binder species exists for the breaking strength after sintering. A comparison of the ReaxFF molecular dynamics simulations using 2-ethoxyethanol, diethylene glycol and 1-(2,2,2-trihydroxyethoxy)ethane-2,2,2-triol as the binder phase indicates that an increasing number of hydroxyl groups leads to higher breaking strength at the print and curing stages. The findings from this study can be extended to identify the optimal binder chemistry, curing and sintering conditions for different material systems.

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

An Integrated Study on a Novel High Temperature High Entropy Alloy

This report summarizes our recent works of theoretical modeling, simulation, and experimental validation of the simulation results on the new refractory high entropy alloy (HEA) design and oxide doped refractory HEA research. The simulation of the stability and thermal dynamics simulation on potential thermal stable candidates were performed and related HEA with oxide doped samples were synthesized and characterized. The HEA ab initio density functional theory and molecular dynamics physical property simulation methods and experimental texture validation techniques development, achievements already reached, course work development, students and postdoc training, and future improvement research directions are briefly introduced.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Transferable density functional tight binding for carbon, hydrogen, nitrogen, and oxygen: Application to shock compression

A new parameterization for density functional tight binding (DFTB) theory, lanl31, has been developed for molecules containing carbon, hydrogen, nitrogen, and oxygen. Optimal values for the Hubbard Us, on-site energies, and the radial dependences of the bond integrals and repulsive potentials were determined by numerical optimization using simulated annealing to a modest database of ab initio-calculated atomization energies and interatomic forces. The transferability of the optimized DFTB parameterization has been assessed using the CHNO subset of the QM-9 database [R. Ramakrishnan et al., Sci. Data, 1, 140022 (2014)]. These analyses showed that the errors in the atomization energies and interatomic forces predicted by our model are small and in the vicinity of the di erences between density functional theory calculations with di erent basis sets and exchange-correlation functionals. Good correlations between the molecular dipole moments and HOMO-LUMO gaps predicted by lanl31 and the QM-9 data set are also found. Furthermore, the errors in the atomization energies and forces derived from lanl31 are signi cantly smaller than those obtained from the ReaxFF-lg reactive force eld for organic materials [L. Liu et al., J. Phys. Chem. A, 115, 11016 (2011)]. The lanl31 DFTB parameterization for C, H, N, and O has been applied the to the molecular dynamics simulation of the principal Hugoniot of liquid nitromethane, liquid benzene, liquid nitrogen, pentaerythritol tetranitrate, trinitrotoluene, and cyclotetramethylene tetranitramine. The computed and measured Hugoniot loci are in excellent agreement with experiment and we discuss the sensitivity of the loci to the underestimated shock heating that is a characteristic of classical molecular dynamics simulations.

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

Machine learning predictions of diffusion in bulk and confined ionic liquids using simple descriptors

Ionic liquids have many intriguing properties and widespread applications such as separations and energy storage. However, ionic liquids are complex fluids and predicting their behavior is difficult, particularly in confined environments. We introduce fast and computationally efficient machine learning (ML) models that can predict diffusion coefficients and ionic conductivity of bulk and nanoconfined ionic liquids over a wide temperature range (350–500 K). The ML models are trained on molecular dynamics simulation data for 29 unique ionic liquids as bulk fluids and confined in graphite slit pores. This model is based on simple physical descriptors of the cations and anions such as molecular weight and surface area. Here, we also demonstrate that accurate results can be obtained using only descriptors derived from SMILES (simplified molecular-input line-entry system) codes for the ions with minimal computational effort. This offers a fast and efficient method for estimating diffusion and conductivity of nanoconfined ionic liquids at various temperatures without the need for expensive molecular dynamics simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Structural connection between gallium crystals and near- T m liquids under ambient pressure

We study the short-range structural order in liquid gallium prepared by ab initio molecular dynamics simulations. We found that at 400 K, which is close to the melting point (Tm) of Ga, the dominant motif in the liquid phase is identical to that in the stable solid phase Ga-I. Strong directional bonding in this motif prevents it from nucleation in the liquid phase. Meanwhile, a newly identified motif, topologically distinct from yet related to those in β-Ga and Ga-III phases, is also abundant in Ga liquid. This new motif could serve as precursors for β-Ga or Ga-III during crystallization.

Zhang, F.↗