Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic 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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Thermal conductivity of MoS 2 monolayers from molecular dynamics simulations

Quantification of lattice thermal conductivity of two-dimensional semiconductors like MoS 2 is necessary for the design of electronic and thermoelectric devices, but direct experimental measurements on free-standing samples is challenging. Molecular dynamics simulations, with appropriate corrections, can provide a reference value for thermal conductivity for these material systems. Here, we construct a new empirical forcefield of the Stillinger-Weber form, parameterized to phonon dispersion relations, lattice constants and elastic moduli and we use it to compute a material-intrinsic thermal conductivity of 38.1 W/m-K at room temperature and estimate a maximum thermal conductivity of 85.4 W/m-K at T = 200 K. We also identify that phonon scattering by the large isotopic mass distribution of Mo and S contributes a significant correction (>45%) to the thermal conductivity at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Molecular dynamics simulation of phase competition in terbium

The competition among multiple solid phases determines the final microstructures of a material. Such competition can originate at the very beginning of the solidification process. We report the results of molecular dynamics simulation of the phase competition between the hexagonal close-packed (hcp), face-centered cubic (fcc), and body-centered cubic (bcc) phases during the solidification of pure Tb. We found that the liquid supercooled below the hcp melting temperature has both bcc and hcp/fcc nuclei, but only the bcc nuclei grow such that the liquid always solidifies into the bcc phase, even at temperatures where the hcp phase is more stable. The hcp phase can only form in the last liquid droplet or at the bcc grain boundaries. Depending on the bcc grain orientations, the hcp phase jammed between the bcc grains either completely disappears or slowly grows via a solid-state massive transformation mechanism. Once the hcp phase becomes large enough, the stresses associated with its appearance can trigger a martensitic transformation. Yet, not the entire bcc phase is consumed by the martensitic transformation and the remaining bcc phase is transformed into the hcp phase via the solid-state massive transformation mechanism. If the supercooling is too large, the nucleation becomes almost barrier free and the liquid solidifies into a structure consisting of ultra-fine hcp and bcc grains after which the bcc phase quickly disappears.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Anion pairs in room temperature ionic liquids predicted by molecular dynamics simulation, verified by spectroscopic characterization

Molecular-level spectroscopic analyses of an aprotic and a protic room-temperature ionic liquid, BMIM OTf and BMIM HSO4, respectively, have been carried out with the aim of verifying molecular dynamics simulations that predict anion pair formation in these fluid structures. Fourier-transform infrared spectroscopy, Raman spectroscopy and nuclear magnetic resonance spectroscopy of various nuclei support the theoretically-determined average molecular arrangements.

ionic liquids↗
DOE OSTI2021

Ion Pairing, Clustering and Transport in a LiFSI-TMP Electrolyte as Functions of Salt Concentration using Molecular Dynamics Simulations

Battery capacity is highly related to ion-pairing mechanisms in electrolytes, since a cluster formation can lead to dead Li formation, reducing the number of charge carriers and leading to capacity fading. We use molecular dynamics simulations to model an electrolyte comprising trimethyl phosphate (TMP) solvent and a lithium bis(fluorosulfonyl)imide (LiFSI) salt, exploring effects of salt concentration on solvation and ion-transport. We simulate the LiFSI-TMP electrolyte for salt concentrations of 0.7, 1.43 and 3.82 molar. A statistical analysis was performed to study ion-pairing, clustering, diffusivity, conductivity, and coordination of Li-ions, providing insights into relations between molecular structures and transport properties. Molecular structure of ionic components changes as concentration increases, from a predominant solvent separated ion pair (SSIP) and contact ion pair (CIP) to aggregate salt (AGG) and ionic cluster formation. Given the formation of the ionic cluster, the diffusion mechanism followed by Li-ions changes from a hopping/exchange to a vehicular mechanism as concentration increases; this is reflected in a decrease of ionic conductivities. Ionicity was also calculated to reveal how the ionic motion changes from an uncorrelated to a correlated one as the salt concentration increases. Furthermore, we also compared our results with experimental calculations performed for similar electrolyte systems

25 ENERGY STORAGE↗
NASA NTRS1975

Dispersion analysis and linear error analysis capabilities of the space vehicle dynamics simulation program

Previous error analyses conducted by the Guidance and Dynamics Branch of NASA have used the Guidance Analysis Program (GAP) as the trajectory simulation tool. Plans are made to conduct all future error analyses using the Space Vehicle Dynamics Simulation (SVDS) program. A study was conducted to compare the inertial measurement unit (IMU) error simulations of the two programs. Results of the GAP/SVDS comparison are presented and problem areas encountered while attempting to simulate IMU errors, vehicle performance uncertainties and environmental uncertainties using SVDS are defined. An evaluation of the SVDS linear error analysis capability is also included.

Snow, L. S.↗
DOE OSTI2023

Quantum mechanics based non-bonded force field functions for use in molecular dynamics simulations of materials and systems: The nitrogen and oxygen columns

Accurate Force Fields (FFs) are essential for Molecular Dynamics (MD) simulations of the dynamics of realistic materials in terms of atomic-level interactions. The FF parameters of short-range valence interactions can be derived through Quantum Mechanical (QM) calculations on model systems practical for QM (<300 atoms). Similarly, the dynamic electrostatic interactions can be described with methods such as QEq or PQEq that allow charges and polarization to adjust dynamically. However, accurately extracting long-range van der Waals (vdW) interactions from QM calculations poses challenges due to the absence of a definitive method to distinguish between the different energetic components of electrostatics, polarization, vdW, hydrogen bonding, and valence interactions. To do this we use the Perdew–Burke–Ernzerhof flavor of Density Functional Theory, including empirical D3 vdW corrections, to predict the Equation of State for each element (keeping any covalent bonds fixed), from which we obtain the two-body vdW nonbond potential. Here, we extend these calculations to include non-bonded parameters for the N and O columns of the periodic table so that we now describe columns 15 (N), 16 (O), 17 (F), and 18 (Ne) of the periodic table. For these 20 elements, we find that the two-body vdW potentials can all be mapped to a single universal two-body curve, with just three scaling parameters: Re, De, and L. We refer to this as the Universal NonBond (UNB) potential. We expect this to be useful for new MD simulations and a helpful starting point to obtain UNB parameters for the remainder of the periodic table.

Chemistry↗
DOE OSTI2020

Molecular dynamics simulations of alkaline earth metal ions binding to DNA reveal ion size and hydration effects

The identity of metal ions surrounding DNA is key to its biological function and materials applications. In this work, we compare atomistic molecular dynamics simulations of double strand DNA (dsDNA) with four alkaline earth metal ions (Mg 2+ , Ca 2+ , Sr 2+ , and Ba 2+ ) to elucidate the physical interactions that govern DNA–ion binding. Simulations accurately model the ion–phosphate distance of Mg 2+ and reproduce ion counting experiments for Ca 2+ , Sr 2+ , and Ba 2+ . Our analysis shows that alkaline earth metal ions prefer to bind at the phosphate backbone compared to the major groove and negligible binding occurs in the minor groove. Larger alkaline earth metal ions with variable first solvation shells (Ca 2+ , Sr 2+ , and Ba 2+ ) show both direct and indirect binding, where indirect binding increases with ion size. Mg 2+ does not fit this trend because the strength of its first solvation shell predicts indirect binding only. Ions bound to the phosphate backbone form fewer contacts per ion compared to the major groove. Within the major groove, metal ions preferentially bind to guanine–cystosine base pairs and form simultaneous contacts with the N7 and O6 atoms of guanine. Overall, we find that the interplay among ion size, DNA–ion interaction, and the size and flexibility of the first solvation shell are key to predicting how alkaline earth metal ions interact with DNA.

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

Molecular dynamics simulation of hollow thick-walled cylinder collapse

The generation and evolution of plastic deformation in a hollow single-crystal cylinder under high-rate axisymmetric loading were studied. An advantage of the proposed loading scheme is that all loading modes are applied simultaneously within the chosen crystallographic plane of the cylinder base and different strain degrees are achieved along the specimen cross section. Molecular dynamics simulation was performed to show that the achievement of a certain strain causes the formation of structural defects on the inner surface of the specimen. The obtained results can be used to explain the main plastic deformation mechanisms of crystalline solids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Molecular Dynamics Simulation of the Structure and Properties of CaO-SiO{sub 2}-CaF{sub 2} Slag Systems

As a fundamental slag system in metallurgical processes, the CaO-SiO{sub 2}-CaF{sub 2} system, its properties, and related microstructures are of critical importance. The effect of the binary basicity and addition of CaF{sub 2} on the structure and viscosity of CaO-SiO{sub 2}-CaF{sub 2} slag were investigated using molecular dynamics simulations and Fourier transform infrared spectrometry. The average bond lengths of Si-O, Ca-O, and O-O were 1.61, 2.31, and 2.61 Å, respectively, and the coordinate number of Si remained constant. O-Si-O maintained an angle of 109.2, which was considered to be unaffected by the basicity changes. With increasing basicity, the network connectivity degrees Q{sup 3} and Q{sup 4} decreased, and the microstructure of the melt was significantly simpler. The effect of CaF{sub 2} on the depolymerization of the network structure was not obvious, and CaF{sub 2} acted mainly as a diluent. The calculations showed that the viscosity decreased upon setting up a ratio of nonbridging oxygen over tetrahedrally coordinated atoms (NBO/T). A good linear relationship was confirmed between the viscosity of the molten slag and the NBO/T.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Dielectric relaxation of ethylene carbonate and propylene carbonate from molecular dynamics simulations

Ethylene carbonate (EC) and propylene carbonate (PC) are widely used solvents in lithium (Li)-ion batteries and supercapacitors. Ion dissolution and diffusion in those media are correlated with solvent dielectric responses. Here, we use all-atom molecular dynamics simulations of the pure solvents to calculate dielectric constants and relaxation times, and molecular mobilities. The computed results are compared with limited available experiments to assist more exhaustive studies of these important characteristics. As a result, the observed agreement is encouraging and provides guidance for further validation of force-field simulation models for EC and PC solvents.

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

Structures of fluoride containing aluminosilicate low activity nuclear waste glasses: A molecular dynamics simulations study

Fluorine and other halides commonly exist in nuclear waste forms, and due to their volatile nature, halide retention poses an issue affecting waste loading during vitrification. The compositional effect on fluorine incorporation in aluminosilicate glasses is investigated through molecular dynamics simulations. Oxygen and fluorine coordination numbers around glass former and modifier cations, bond angle distributions, and medium range structure features such as Q n distributions, ring size distributions and neutron diffraction structure factors were calculated. It was found that fluorine has higher preference to bond to Ca 2+ than to Na + , both in the melt and the glass, and there is no Si-F bond formation in the glass but they do exist in the melt. Consequently, CaO for Na 2 O substitution can be an effective way to help fluorine retention without significantly changing the glass chemistry. Furthermore, these results thus provide insights on fluorine incorporation in the aluminosilicate nuclear waste glasses and the strategy on how to improve fluorine retention both in the glass and the melt.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Testing Dynamic Simulation Models for Different Types of Advanced Pumped Storage Hydro Units

A major task of the U.S. Department of Energy (DOE)-sponsored project “Modeling and Analysis of Value of Advanced Pumped Storage Hydropower in the U.S.” was to develop new dynamic simulation models to represent advanced pumped storage hydro (PSH) technologies. This report describes the testing that was performed to demonstrate the performance of these simulation models and illustrates how these models can now be used in analyses required for investigations into applications of these technologies.

Koritarov, Vladimir↗
DOE OSTI2018

Validating Molecular Dynamics Simulations against Experimental Observables in Light of Underlying Conformational Ensembles

Far from the static, idealized conformations deposited into structural databases, proteins are highly dynamic molecules that undergo conformational changes on temporal and spatial scales that may span several orders of magnitude. These conformational changes, often intimately connected to the functional roles that proteins play, may be obscured by traditional biophysical techniques. Over the past 40 years, molecular dynamics (MD) simulations have complemented these techniques by providing the “hidden” atomistic details that underlie protein dynamics. However, there are limitations of the degree to which molecular simulations accurately and quantitatively describe protein motions. Here we show that although four molecular dynamics simulation packages (AMBER, GROMACS, NAMD, and il mm) reproduced a variety of experimental observables for two different proteins (engrailed homeodomain and RNase H) equally well overall at room temperature, there were subtle differences in the underlying conformational distributions and the extent of conformational sampling obtained. This leads to ambiguity about which results are correct, as experiment cannot always provide the necessary detailed information to distinguish between the underlying conformational ensembles. Furthermore, the results with different packages diverged more when considering larger amplitude motion, for example, the thermal unfolding process and conformational states sampled, with some packages failing to allow the protein to unfold at high temperature or providing results at odds with experiment. While most differences between MD simulations performed with different packages are attributed to the force fields themselves, there are many other factors that influence the outcome, including the water model, algorithms that constrain motion, how atomic interactions are handled, and the simulation ensemble employed. Here four different MD packages were tested each using best practices as established by the developers, utilizing three different protein force fields and three different water models. Differences between the simulated protein behavior using two different packages but the same force field, as well as two different packages with different force fields but the same water models and approaches to restraining motion, show how other factors can influence the behavior, and it is incorrect to place all the blame for deviations and errors on force fields or to expect improvements in force fields alone to solve such problems.

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

Molecular dynamics simulations of stratum corneum lipid mixtures: A multiscale perspective

The lipid matrix of the stratum corneum (SC) layer of skin is essential for human survival; it acts as a barrier to prevent rapid dehydration while keeping potentially hazardous material outside the body. While the composition of the SC lipid matrix is known, the molecular-level details of its organization are difficult to infer experimentally, hindering the discovery of structure-property relationships. To this end, molecular dynamics simulations, which give molecular-level resolution, have begun to play an increasingly important role in understanding these relationships. Furthermore, most simulation studies of SC lipids have focused on preassembled bilayer configurations, which, owing to the slow dynamics of the lipids, may influence the final structure and hence the calculated properties. Self-assembled structures would avoid this dependence on the initial configuration, however, the size and length scales involved make self-assembly impractical to study with atomistic models. Here, we report on the development of coarse-grained models of SC lipids designed to study self-assembly. Building on previous work, we present the interactions between the headgroups of ceramide and free fatty acid developed using the multistate iterative Boltzmann inversion method. Validation of the new interactions is performed with simulations of preassembled bilayers and good agreement between the atomistic and coarse-grained models is found for structural properties. The self-assembly of mixtures of ceramide and free fatty acid is investigated and both bilayer and multilayer structures are found to form. This work therefore represents a necessary step in studying SC lipid systems on multiple time and length scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Thermal conductivity of tungsten: Effects of plasma-related structural defects from molecular-dynamics simulations

We report results on the lattice thermal conductivities of tungsten single crystals containing nanoscale-sized pores or voids and helium (He) nanobubbles as a function of void/bubble size and gas pressure in the He bubbles based on molecular-dynamics simulations. For reference, we calculated lattice thermal conductivities of perfect tungsten single crystals along different crystallographic directions at room temperature and found them to be about 10% of the overall thermal conductivity of tungsten with a weak dependence on the heat flux direction. The presence of nanoscale voids in the crystal causes a significant reduction in its lattice thermal conductivity, which decreases with increasing void size. Filling the voids with He to form He nanobubbles and increasing the bubble pressure leads to further significant reduction of the tungsten lattice thermal conductivity, down to ~20% of that of the perfect crystal. The anisotropy in heat conduction remains weak for tungsten single crystals containing nanoscale-sized voids and He nanobubbles throughout the pressure range examined. Analysis of the pressure and atomic displacement fields in the crystalline region that surrounds the He nanobubbles reveals that the significant reduction of tungsten lattice thermal conductivity in this region is due to phonon scattering from the nanobubbles, as well as lattice deformation around the nanobubbles and formation of lattice imperfections at higher bubble pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1991

Gain in computational efficiency by vectorization in the dynamic simulation of multi-body systems

An improved technique for the identification and extraction of the exact quantities associated with the degrees of freedom at the element as well as the flexible body level is presented. It is implemented in the dynamic equations of motions based on the recursive formulation of Kane et al. (1987) and presented in a matrix form, integrating the concepts of strain energy, the finite-element approach, modal analysis, and reduction of equations. This technique eliminates the CPU intensive matrix multiplication operations in the code's hot spots for the dynamic simulation of the interconnected rigid and flexible bodies. A study of a simple robot with flexible links is presented by comparing the execution times on a scalar machine and a vector-processor with and without vector options. Performance figures demonstrating the substantial gains achieved by the technique are plotted.

Amirouche, F. M. L.↗
DOE OSTI2016

Homogeneous hydride formation path in α-Zr: Molecular dynamics simulations with the charge-optimized many-body potential

A formation path for homogeneous γ hydride formation in hcp α-Zr, from solid solution to the ζ and then the γ hydride, was demonstrated using molecular static calculations and molecular dynamic simulations with the charge-optimized many-body (COMB) potential. Hydrogen has limited solubility in α-Zr. Once the solubility limit is exceeded, the stability of solid solution gives way to that of coherent hydride phases such as the ζ hydride by planar precipitation of hydrogen. At finite temperatures, the ζ hydride goes through a partial hcp-fcc transformation via 1/3 <1¯100> slip on the basal plane, and transforms into a mixture of γ hydride and α-Zr. In the ζ hydride, slip on the basal plane is favored thermodynamically with negligible barrier, and is therefore feasible at finite temperatures without mechanical loading. The transformation process involves slips of three equivalent shear partials, in contrast to that proposed in the literature where only a single shear partial was involved. The adoption of multiple slip partials minimizes the macroscopic shape change of embedded hydride clusters and the shear strain accumulation in the matrix, and thus reduces the overall barrier needed for homogeneous γ hydride formation. In conclusion, this formation path requires finite temperatures for hydrogen diffusion without mechanical loading. Therefore, it should be effective at the cladding operating conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2016

Experimental and molecular dynamics simulation study of structure of liquid and amorphous Ni 62 Nb 38 alloy

The state-of-the-art experimental and atomistic simulation techniques were utilized to study the structure of the liquid and amorphous Ni 62 Nb 38 alloy. First, the ab initio molecular dynamics (AIMD) simulation was performed at rather high temperature where the time limitations of the AIMD do not prevent to reach the equilibrium liquid structure. A semi-empirical potential of the Finnis-Sinclair (FS) type was developed to almost exactly reproduce the AIMD partial pair correlation functions (PPCFs) in a classical molecular dynamics simulation. This simulation also showed that the FS potential well reproduces the bond angle distributions. The FS potential was then employed to elongate the AIMD PPCFs and determine the total structure factor (TSF) which was found to be in excellent agreement with X-ray TSF obtained within the present study demonstrating the reliability of the AIMD for the simulation of the structure of the liquid Ni–Nb alloys as well as the reliability of the developed FS potential. The glass structure obtained with the developed potential was also found to be in excellent agreement with the X-ray data. We report the analysis of the structure revealed that a network of the icosahedra clusters centered on Ni atoms is forming during cooling the liquid alloy down to T g and the Nb Z14, Z15, and Z16 clusters are attached to this network. This network is the main feature of the Ni 62 Nb 38 alloy and further investigations of the properties of this alloy should be based on study of the behavior of this network.

36 MATERIALS SCIENCE↗