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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Parallel point defect identification in molecular dynamics simulations without post-processing: A compute and dump style for LAMMPS

Two classes for the molecular dynamics program LAMMPS – one compute style and one dump style – are presented that are designed to identify, count, and output point defects in cascade damage and related molecular dynamics simulations. The calculations are done in parallel across multiple MPI processes, and can be done at the same time as the original simulation. This drastically reduces storage requirements by eliminating the need to post-process every atom in the system. Non-cubic lattices, free surfaces, and large voids can be eliminated from the output by suitable choices of the reference lattice. The classes are derived from LAMMPS’s Dump and Compute classes, and pose no additional overhead to LAMMPS if they are not used.

42 ENGINEERING↗
DOE OSTI2019

Understanding the release of helium atoms from nanochannel tungsten: a molecular dynamics simulation

The design of highly radiation-tolerant plasma-facing materials (PFMs) is of great importance for fusion reactors. Our recent experiments have shown that nanochannel tungsten (W) films have clearly superior radiation tolerance properties. In the present work, helium clustering and release from nanochannel tungsten were studied by molecular dynamics simulations. The effects of temperature and vacancy concentration on the helium release from a tungsten cylinder were investigated. Our results show that nanochannel W that consists of thin W cylinders releases He atoms more quickly than bulk W with flat surfaces, thus greatly reducing the He concentration and suppressing the formation and growth of He bubbles, which leads to increased radiation tolerance. Moreover, the micro-structural changes due to increasing He fluence are smaller in nanochannel W than those in bulk W. Although vacancies in nanochannel W will trap He atoms, the nanochannel W also has a stronger tendency to stabilize helium retention than bulk W. The mechanism of helium release from nanochannel W was also examined. The results reported here are beneficial for guiding future work in the design of radiation resistant PFMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Variational control forces for enhanced sampling of nonequilibrium molecular dynamics simulations

We introduce a variational algorithm to estimate the likelihood of a rare event within a nonequilibrium molecular dynamics simulation through the evaluation of an optimal control force. Optimization of a control force within a chosen basis is made possible by explicit forms for the gradients of a cost function in terms of the susceptibility of driven trajectories to changes in variational parameters. We consider probabilities of time-integrated dynamical observables as characterized by their large deviation functions and find that in many cases, the variational estimate is quantitatively accurate. Additionally, we provide expressions to exactly correct the variational estimate that can be evaluated directly. We benchmark this algorithm against the numerically exact solution of a model of a driven particle in a periodic potential, where the control force can be represented with a complete basis. Additionally, we then demonstrate the utility of the algorithm in a model of repulsive particles on a line, which undergo a dynamical phase transition, resulting in singular changes to the form of the optimal control force. In both systems, we find fast convergence and are able to evaluate large deviation functions with significant increases in statistical efficiency over alternative Monte Carlo approaches.

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

Molecular dynamics simulations of concentration-dependent defect production in Fe-Cr and Fe-Cu alloys

Molecular dynamics simulations are conducted to study the effects of alloying elements on the primary damage behaviors in three Fe-based ferritic alloy systems: (1) a Fe-Cr system in which the heat of mixing changes its sign with the Cr concentration; (2) a Fe-Cu system that has a positive heat of mixing; and (3) an ideal but artificial Fe-Cr system that has a zero heat of mixing, which is used as a reference system to investigate solute interstitial formation based on probability. It is found that in these alloys, the solute type and concentration do not have a significant effect on the total number of surviving Frenkel pairs. However, the fraction of solute interstitials has distinct behaviors. In Fe-Cr, the Cr interstitial fraction is much higher than the Cr solute concentration and the Cr interstitial production efficiency decreases with the increasing Cr concentration. By contrast, in Fe-Cu, Cu interstitials are barely produced. In the ideal alloy, the solute interstitial fraction is close to the solute concentration. The defect formation energies in both dilute and concentrated alloys, interstitial binding energies, liquid diffusivities of Fe and solute atoms, and heat of mixing have been calculated for both Fe-Cr and Fe-Cu alloys. Among them, we find that the relative thermodynamic stability between Fe self-interstitials and solute interstitials plays the most important role in the solute interstitial production behaviors. The decrease of Cr interstitial production efficiency with increasing Cr concentration can be explained by the probability distribution functions of solute interstitial formation energy in concentrated alloys.

Zhang, Yaxuan↗
DOE OSTI2017

Thermal Hydraulic Computational Fluid Dynamics Simulations and Experimental Investigation of Deformed Fuel Assemblies

The project, Toward a Longer Life Core: Thermal Hydraulic CFD Simulations and Experimental Investigation of Deformed Fuel Assemblies, DOE Project code DE-NE0008321, was a verification and validation project for flow and heat transfer through wire wrapped simulated liquid metal fuel assemblies that included both experiments and computational fluid dynamics simulations of those experiments. This project was a two year collaboration between AREVA, TerraPower, Argonne National Laboratory and Texas A&M University. Experiments were performed by AREVA and Texas A&M University. Numerical simulations of these experiments were performed by TerraPower and Argonne National Lab. Project management was performed by AREVA Federal Services. The first of a kind project resulted in the production of both local point temperature measurements and local flow mixing experiment data paired with numerical simulation benchmarking of the experiments. The project experiments included the largest wire-wrapped pin assembly Mass Index of Refraction (MIR) experiment in the world, the first known wire-wrapped assembly experiment with deformed duct geometries and the largest numerical simulations ever produced for wire-wrapped bundles.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1996

Applying Parallel Processing Techniques to Tether Dynamics Simulation

The focus of this research has been to determine the effectiveness of applying parallel processing techniques to a sizable real-world problem, the simulation of the dynamics associated with a tether which connects two objects in low earth orbit, and to explore the degree to which the parallelization process can be automated through the creation of new software tools. The goal has been to utilize this specific application problem as a base to develop more generally applicable techniques.

Wells, B. Earl↗
DOE OSTI2025

Combining Reactive Quantum-Mechanical Molecular-Dynamics Simulations with Mutagenesis, Crystallography, and Enzyme Kinetics to Reveal Plausible Steps of Isocyanide Hydratase Catalysis

A complete understanding of enzyme mechanisms requires atomistic details of chemical reactions. Quantum-based molecular dynamics simulations (QMD) are a potential source of this information, but trade-offs between accuracy and computational cost have limited their use. We previously developed extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) methods that leverage a negligible compromise in accuracy to substantially decrease the cost of QMD simulations. Here, we develop a reactive QMD approach using the latest XL-BOMD formulation, which enables efficient simulations of highly reactive systems, and use it to investigate mechanisms of intermediate formation in isocyanide hydratase (ICH) catalysis. In QMD simulations, molecular analogs of ICH active site residues reacted with para-nitrophenyl isocyanide, forming a thioimidate. Analysis of simulated atomic configurational and charge dynamics revealed a pathway where protonation of the isocyanide carbon occurs prior to thioimidate formation and suggested a possible role of Asp17 as a proton donor in the early phase of ICH catalysis. To test whether the pathway seen using the reactive QMD approach might be relevant to ICH catalysis, we performed X-ray crystallography and pre-steady-state enzyme kinetics studies of wild-type and D17N mutant ICH. Both the structure and kinetics are sensitive to the D17N mutation in a manner that is consistent with the order of the reaction steps seen in the simulations. Mobile protons play essential roles in many enzymes, yet they are difficult to observe experimentally, making the ordering of proton-dependent steps ambiguous in many enzyme mechanisms. The ability to directly simulate model reactions for the design of experiments that provide information about enzyme mechanisms involving mobile protons demonstrates the significance of our reactive QMD approach and motivates further biological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Internal protein motions in molecular-dynamics simulations of Bragg and diffuse X-ray scattering

Molecular-dynamics (MD) simulations of Bragg and diffuse X-ray scattering provide a means of obtaining experimentally validated models of protein conformational ensembles. This paper shows that compared with a single periodic unit-cell model, the accuracy of simulating diffuse scattering is increased when the crystal is modeled as a periodic supercell consisting of a 2 × 2 × 2 layout of eight unit cells. The MD simulations capture the general dependence of correlations on the separation of atoms. There is substantial agreement between the simulated Bragg reflections and the crystal structure; there are local deviations, however, indicating both the limitation of using a single structure to model disordered regions of the protein and local deviations of the average structure away from the crystal structure. Although it was anticipated that a simulation of longer duration might be required to achieve maximal agreement of the diffuse scattering calculation with the data using the supercell model, only a microsecond is required, the same as for the unit cell. Rigid protein motions only account for a minority fraction of the variation in atom positions from the simulation. The results indicate that protein crystal dynamics may be dominated by internal motions rather than packing interactions, and that MD simulations can be combined with Bragg and diffuse X-ray scattering to model the protein conformational ensemble.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Molecular Dynamics Simulations of Complexation of Am(III) with a Preorganized Dicationic Ligand in an Ionic Liquid

Preorganized ligands with imidazolium arms have been found to be highly selective in extracting Am(III) into ionic liquids (ILs), but the detailed structure and mechanism of the complexation process in the ionic solvation environment are unclear. Here, we carry out molecular dynamics simulation of the complexation of Am(III) with a preorganized 1,10-phenanthroline-2,9-dicarboxamide complexant (L) functionalized with alkyl chains and imidazolium cations in the butylmethylimidazolium bistriflimide ([BMIM][NTf 2 ]) IL. Both Am:L (1:1) and Am:L 2 (1:2) complexes are examined. In the absence of the ligand, Am(III) is found to be coordinated by six NTf 2 anions via nine O donors in the first solvation shell. In the Am:L complex, Am(III) is coordinated to the ligand via two O donors and four NTf 2 anions via seven O donors in the first coordination shell. In the Am:L 2 complex, Am(III) is coordinated to the two ligands via four O donors and four NTf 2 anions via five O donors. The imidazolium arms of the ligands play an important role in the secondary solvation environment by attracting NTf 2 anions closer to the metal center. As a result, we find that the binding free energy for the second L 2+ ligand is twice that for the first L 2+ ligand, making the Am:L 2 complex significantly more stable than the Am:L complex. This study highlights the multiple factors and tunability in using preorganized ligands with charged functional groups in an ionic solvation environment, which could hold the key to achieving desired selectivity in ion extraction efficiency.

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

Dynamics of the EAG1 K{sup +} channel selectivity filter assessed by molecular dynamics simulations

EAG1 channels belong to the KCNH family of voltage gated potassium channels. They are expressed in several brain regions and increased expression is linked to certain cancer types. Recent cryo-EM structure determination finally revealed the structure of these channels in atomic detail, allowing computational investigations. In this study, we performed molecular dynamics simulations to investigate the ion binding sites and the dynamical behavior of the selectivity filter. Our simulations suggest that sites S2 and S4 form stable ion binding sites, while ions placed at sites S1 and S3 rapidly switched to sites S2 and S4. Further, ions tended to dissociate away from S0 within less than 20 ns, due to increased filter flexibility. This was followed by water influx from the extracellular side, leading to a widening of the filter in this region, and likely non-conductive filter configurations. Simulations with the inactivation-enhancing mutant Y464A or Na{sup +} ions lead to trapped water molecules behind the SF, suggesting that these simulations captured early conformational changes linked to C-type inactivation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Stochastic symplectic reduced-order modeling for model-form uncertainty quantification in molecular dynamics simulations in various statistical ensembles

Here, this work focuses on the representation of model-form uncertainties in molecular dynamics simulations in various statistical ensembles. In prior contributions, the modeling of such uncertainties was formalized and applied to quantify the impact of, and the error generated by, pair-potential selection in the microcanonical ensemble (NVE). In this work, we extend this formulation and present a linear-subspace reduced-order model for the canonical (NVT) and isobaric (NPT) ensembles. The symplectic reduced-order basis is randomized on the tangent space of the Stiefel manifold to provide topological relationships and capture model-form uncertainty. Using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), we assess the relevance of these stochastic reduced-order atomistic models on canonical problems involving a Lennard-Jones fluid and an argon crystal melt.

42 ENGINEERING↗
NASA NTRSDate not supplied

Molecular dynamics simulation of effects of solutes on dislocation propagation in Ni-based superalloys

Ni-based superalloys are used in the hot sections of jet turbine engines because of their high strength, stability and resistance to oxidation. The properties of these alloys can be further optimized by adding alloying elements. Therefore, a fundamental knowledge on the effect of different elements on properties of Ni-based superalloys is required. Molecular dynamics simulation could shed light here but its application is hindered by the absence of reliable and computationally cheap semi-empirical potential of the interatomic interaction for 4 and more element alloys. We will present a new Ni-Al-Cr-Nb Finnis-Sinclair (FS) potential specially designed to simulate the dislocation propagation from to  phase. In order to construct this potential, we designed a special algorithm to incorporate the data on element partitioning in the potential development procedure. For example, it is known from experiment, that Cr is mostly present in the gamma phase. Figure 1 shows a snapshot obtained after equilibration of the model of the Ni68Al17Cr15 alloy at T=1000 K using the hybrid Monte-Carlo (MC)/molecular dynamics (MD) simulation with the developed semi-empirical potential. One can clearly see that the Cr partitioning is in agreement with the experimental data. We will discuss the developed algorithm to incorporate the solute partition data in details. Using the developed semi-empirical potential, we first investigated the effect of anti-site defects in the  phase on the single dislocation propagation. It was found that the dislocation velocity increases with the increasing of the anti-site defect concentration. This effect was attributed to smaller number of Al-Al pairs forming during the dislocation migration in the presence of the anti-site defects. Next, we investigated the effect of Nb on the dislocation pair propagation in the Kolbe mechanism. It was found that the addition of Nb leads to considerable decrease in the dislocation propagation rate. This is in agreement with the experimental data on the effect of Nb on the creep resistance of the Ni-based superalloys. We will discuss the origin of this effect.

Mikhail I. Mendelev↗
NASA NTRS1990

Gamma Ray Observatory (GRO) dynamics simulator requirements and mathematical specifications, revision 1

The requirements and mathematical specifications for the Gamma Ray Observatory (GRO) Dynamics Simulator are presented. The complete simulator system, which consists of the profie subsystem, simulation control and input/output subsystem, truth model subsystem, onboard computer model subsystem, and postprocessor, is described. The simulator will be used to evaluate and test the attitude determination and control models to be used on board GRO under conditions that simulate the expected in-flight environment.

Harman, R.↗
NASA NTRS1988

The use of the articulated total body model as a robot dynamics simulation tool

The Articulated Total Body (ATB) model is a computer sumulation program which was originally developed for the study of aircrew member dynamics during ejection from high-speed aircraft. This model is totally three-dimensional and is based on the rigid body dynamics of coupled systems which use Euler's equations of motion with constraint relations of the type employed in the Lagrange method. In this paper the use of the ATB model as a robot dynamics simulation tool is discussed and various simulations are demonstrated. For this purpose the ATB model has been modified to allow for the application of torques at the joints as functions of state variables of the system. Specifically, the motion of a robotic arm with six revolute articulations with joint torques prescribed as functions of angular displacement and angular velocity are demonstrated. The simulation procedures developed in this work may serve as valuable tools for analyzing robotic mechanisms, dynamic effects, joint load transmissions, feed-back control algorithms employed in the actuator control and end-effector trajectories.

Obergfell, Louise A.↗
DOE OSTI2025

Data from "Deep Potential Molecular Dynamics Simulations of Low-Temperature Plasma-Surface Interactions"

Data and input files related to the paper "Deep Potential Molecular Dynamics Simulations of Low-Temperature Plasma-Surface Interactions" (https://doi.org/10.1116/6.0004027). This includes the final DP model used in all simulations, training data set, example input files to run DeepMD (with LAMMPS), and data tables summarizing the results obtained from the simulations.

machine learning models↗
DOE OSTI2017

Solvation and Dynamics of Sodium and Potassium in Ethylene Carbonate from ab Initio Molecular Dynamics Simulations

The development of sodium and potassium batteries offers a promising way to meet the scaling and cost challenges of energy storage. However, compared to Li + , several intrinsic properties of Na + and K + , including their solvation and dynamics in typical organic electrolytes utilized in battery applications, are less well-understood. Here in this paper, we report a systematic investigation of Na + and K + in ethylene carbonate (EC) using first-principles molecular dynamics simulations. Our simulations reveal significant differences in the solvation structure and dynamical properties of Na + and K + compared to Li + . We find that, in contrast to Li + which exhibits a well-defined first solvation shell, the larger Na+ and K+ ions show more disordered and flexible solvation structures. These differences in solvation were found to significantly influence the ion dynamics, leading to larger diffusion coefficients of Na + and K + compared to Li + . Our simulations also reveal a clear and interesting analog in the behavior of the ions in EC and aqueous environments, particularly in the specific ion effects on the solvent dynamics. Lastly, this work provides fundamental understanding of the intrinsic properties of Na + and K + in organic electrolytes, which may ultimately influence the intercalation mechanism at the electrode–electrolyte interface and therefore battery performance, lifetime, and safety.

25 ENERGY STORAGE↗
DOE OSTI2022

Constant-potential molecular dynamics simulations of molten salt double layers for FLiBe and FLiNaK

We report the results of constant-potential molecular dynamics simulations of the double layer interface between molten 2LiF–BeF 2 (FLiBe) and 23LiF–6NaF–21KF (FLiNaK) fluoride mixtures and idealized solid electrodes. Employing methods similar to those used in studies of chloride double layers, we compute the structure and differential capacitance of molten fluoride electric double layers as a function of applied voltage. The role of molten salt structure is probed through comparisons between FLiBe and FLiNaK, which serve as models for strong and weak associate-forming salts, respectively. In FLiBe, screening involves changes in Be–F–Be angles and alignment of the oligomers parallel to the electrode, while in FLiNaK, the electric field is screened mainly by rearrangement of individual ions, predominantly the polarizable potassium cation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗