Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Real-time dynamics simulation of the Cassini spacecraft using DARTS. Part 1: Functional capabilities and the spatial algebra algorithm

This paper describes the Dynamics Algorithms for Real-Time Simulation (DARTS) real-time hardware-in-the-loop dynamics simulator for the National Aeronautics and Space Administration's Cassini spacecraft. The spacecraft model consists of a central flexible body with a number of articulated rigid-body appendages. The demanding performance requirements from the spacecraft control system require the use of a high fidelity simulator for control system design and testing. The DARTS algorithm provides a new algorithmic and hardware approach to the solution of this hardware-in-the-loop simulation problem. It is based upon the efficient spatial algebra dynamics for flexible multibody systems. A parallel and vectorized version of this algorithm is implemented on a low-cost, multiprocessor computer to meet the simulation timing requirements.

Jain, A.↗
DOE OSTI2021

Prediction of Bi 2 Te 3 -Sb 2 Te 3 Interfacial Conductance and Superlattice Thermal Conductivity Using Molecular Dynamics Simulations

Bismuth telluride (Bi 2 Te 3 ) and its alloys with antimony telluride (Sb 2 Te 3 ) have long been considered to be the best room-temperature bulk thermoelectric (TE) materials. In recent decades, proof-of-concept demonstrations on Bi 2 Te 3 -Sb 2 Te 3 nanostructures have shown high TE performance due to reduction in lattice thermal conductivities. Particularly, ultra-low thermal conductivities have been observed in Bi 2 Te 3 -Sb 2 Te 3 1D superlattices, leading to thermoelectric figures of merit (ZT) as high as 2.4. In contrast, very few computational studies have been performed to provide insight into the phonon transport across these nanostructures. In this work, we use non-equilibrium molecular dynamics simulations with previously developed force fields to simulate thermal transport across Bi 2 Te 3 -Sb 2 Te 3 interfaces and superlattices. We first calculate the thermal conductance associated with a Bi 2 Te 3 -Sb 2 Te 3 interface across a temperature range of 200–400 K. Furthermore, the values are also compared with thermal conductances calculated by a modified Landauer transport formalism using phonon transmission coefficients obtained from the diffuse mismatch model. Our results show that inelastic scattering processes contribute to an increase in interfacial thermal conductance at higher temperatures. Finally, we calculate the thermal conductivities of Bi 2 Te 3 -Sb 2 Te 3 superlattices with varying period lengths from 2 to 18 nm. A minimum thermal conductivity of 0.27 W/mK is observed at a period length of 4 nm, which is attributed to the competition between incoherent and coherent phonon transport regimes. In comparison with previous experimental measurements in the literature, our results show good agreement with respect to the range of thermal conductivity values and the period length corresponding to the minimum superlattice thermal conductivity.

42 ENGINEERING↗
DOE OSTI2018

Molecular Dynamics Simulation of Hydration and Swelling of Mixed-Layer Clays

Swelling of clay minerals is important to a broad class of problems in science and engineering. While the problem has been extensively studied experimentally, past molecular modeling of the phenomenon was focused on pure clays of one type or another. In practice, however, there is a diverse class of mixed-layer clays (MLCs) in sedimentary rock with intermixed stacking sequence of two or more types of distinct layers within a single crystal. In fact, more than 60% of sedimentary rocks in the U.S. contain various types of MLCs. In this work, we present the results, to our knowledge, of the first molecular dynamics simulation of hydration energetics and swelling of illite–montmorillonite (I–MMT) MLCs, the most common type of mixed clays. The swelling is studied as a function of the water concentration with three combinations of interlayer cations, namely, Na + and K + . The hydration energies, the radial distribution functions, and the density profiles in the interlayer region are computed. For regular Na–MMT with layer charge concentrated in the octahedral sheet, weak cation–surface interaction results in fully hydrated ions and significant swelling. In the asymmetric interlayer of the MLC, however, the illite sheet with stronger interaction of surface and ions causes adsorption of the cations deep in the ditrigonal cavities of the siloxane surface. Given that the hydration enthalpy of K + is smaller than that of Na + , its hydration shell is quite unstable compared with that of Na+. Therefore, swelling is inhibited as the ratio K + /Na + increases. Finally, the results demonstrate the significant differences between the hydration and swelling properties of pure clays and the mixed ones, which have important implications in practice, particularly for sequestration of CO 2 in sedimentary rock.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Equation of state calculations for two-dimensional dust coulomb crystal at near zero temperature by molecular dynamics simulations

Dust particles observed in universe as well as in laboratory and technological plasma devices are still under investigation. At low temperature, these particles are strongly negatively charged and are able to form a 2D or 3D coulomb crystal. In this work, our aim was to check the ideal gas law validity for a 2D single-layer dust crystal recently reported in the literature. For this purpose, we have simulated, using the molecular dynamics method, its thermodynamic properties for different values of dust particles number and confinement parameters. The obtained results have allowed us to invalidate the ideal gas behaviour and to propose an effective equation of state which assumes a near zero dust temperature. Furthermore, the value of the calculated sound velocity was found to be in a good agreement with experimental data published elsewhere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Fluid-phase helium: Shock-compression experiments, quantum molecular dynamics simulations, and development of an equation of state

Helium (He) plays a critical role in numerous areas ranging from the study of celestial objects like brown dwarfs and gas giants to modern-day technologies like nuclear energy and rocket propulsion. For many of these applications, it is essential to have a reliable equation of state (EOS) for He that yields an accurate representation of its thermodynamic behavior. To help constrain and develop such EOS models, we have performed a series of shock-compression experiments on cryogenic liquid He to pressures exceeding 100 GPa using a magnetically accelerated flyer plate on Sandia National Laboratories' Z-machine. We have also performed quantum molecular dynamics simulations that are consistent with our shock measurements. None of the previously available EOSs agree with our experimental and simulation results, motivating the development of a fluid-phase He EOS that we present in this study. Here, we show that our EOS yields good agreement with published data that span temperatures and pressures encountered across a diverse array of applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Fluid dynamic simulation and analysis of water-cooling systems for the Electron-Ion Collider

The Electron-Ion Collider is the newest large-scale project at Brookhaven National Laboratory. The collider’s purpose is to provide further advancements in the knowledge of the universe’s origin by accelerating particles near the speed of light. Our project for this 3.8 km ring was to create a thermal hydraulic steady-state simulation design of the water-cooling system to be cost-effective and energy efficient, as envisioned by Charlie Foltz, the EIC Infrastructure Division Director. The system would include a supply and return header, which cools several thousand components of the ring. The water would then be returned and cooled down using a system of cooling towers and plate and frame heat exchangers. Due to the size of the system and the complexity of the network analysis, a fluid dynamic simulation software, AFT Fathom, was used. Since previous methods of maintaining systems relied on building upon smaller real-life models and implementing empirical data, this flow model was unique and first of a kind in the domain of accelerator design, construction and operation. Therefore, our hydraulic team piloted a new method to perform network analysis on a large scale cooling system. We successfully created several test scenarios for system behavior in a shorter time compared to the method of performing hand calculations. Cooling specifications for heat rejection, pressure drop, flow rate, and pipe sizing were changed based on the individual systems of the vacuum, radio frequency (RF), magnet and power supply, and cryogenics sections. Finally, we used DOE guidelines to perform life-cycle cost analysis with net present value and carbon saving analysis on the systems where pipe size could be optimized.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Microscale mechanical property variations of Al-substituted LLZO: insights from compression testing and molecular dynamics simulations

Ceramic solid electrolytes based on LLZO (Li 7 La 3 Zr 2 O 12 ) are promising candidates for all-solid-state batteries due to their high ionic conductivity and good apparent stability vs. lithium metal, however they are prone to mechanical failure. Lithium metal intrusions, alongside cell stack pressure, transition polycrystalline solid electrolyte grains into a compressed state that promotes crack propagation and fracture. Here this work examines the mechanical response of Al-substituted LLZO to compressive forces by measuring ultimate strength under pillar compression with a flat punch tip. Failure modes characterized by in situ scanning electron microscopy show diverse splitting patterns arising from internal porosity, grain boundaries, and slip planes. Large correlated variations in compressive strength (0.93–2.63 GPa) and Young's modulus (72.1–150.97 GPa) are observed across microscale regions of the solid electrolyte. Molecular dynamics simulations of LLZO with different porosities describe the variation of compressive strength and Young's modulus, and enable a microscale porosity model to be fit accounting for Young's modulus reduction across the solid electrolyte. Overall, the results indicate the importance of microscale mechanical testing of ceramic solid electrolytes to identify preferential sites for mechanical degradation and Li intrusion, and ensure the robust design of all-solid-state lithium metal batteries.

25 ENERGY STORAGE↗
DOE OSTI2016

Li 3 Y(PS 4 ) 2 and Li 5 PS 4 Cl 2 : New Lithium Superionic Conductors Predicted from Silver Thiophosphates using Efficiently Tiered Ab Initio Molecular Dynamics Simulations

We report two novel, earth-abundant lithium superionic conductors, Li3Y(PS4)2 and Li5PS4Cl2, that are predicted to satisfy the necessary combination of good phase stability, high Li+ conductivity, wide band gap and good electrochemical stability for solid electrolyte applications in all-solid-state rechargeable lithium-ion batteries. These candidates were identified from a high-throughput first-principles screening of the Li–P–S ternary and Li–M–P–S (where M is a non-redox-active element) quaternary chemical spaces, including candidates obtained by replacing Ag with Li in the Ag–P–S and Ag–M–P–S chemical spaces. An efficient tiered screening strategy was developed that combines topological analysis with ab initio molecular dynamics simulations to exclude rapidly candidates unlikely to satisfy the stringent conductivity requirements of lithium superionic conductors. In particular, we find Li3Y(PS4)2 to be an extremely promising candidate exhibiting a room-temperature Li+ conductivity of 2.16 mS/cm, which can be increased multifold to 7.14 and 5.25 mS/cm via aliovalent doping with Ca2+ and Zr4+, respectively. More critically, we show that the phase and electrochemical stability of Li3Y(PS4)2 is expected to be better than current state-of-the-art lithium superionic conductors.

Zhu, Zhuoying↗
DOE OSTI2014

Conformational dynamics of a crystalline protein from microsecond-scale molecular dynamics simulations and diffuse X-ray scattering

Significance A major challenge of protein crystallography is to accurately describe protein structure variations. Whereas Bragg peaks only yield a picture in which the variations are superimposed, diffuse X-ray scattering (diffraction between the peaks) reports on the way different atoms move together and can be used to increase the accuracy of models. We have performed a molecular-dynamics (MD) simulation of crystalline staphylococcal nuclease, yielding a rich picture of functionally important motions. This simulation is 100-fold longer than previous studies and yields reproducible calculations of diffuse intensity, overcoming a major challenge in the field. Experimental data support the MD models and indicate that diffuse scattering can be used both to validate and to improve MD simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Molecular Dynamics Simulation and Experimental Studies of Gold Nanoparticle Templated HDL-like Nanoparticles for Cholesterol Metabolism Therapeutics

High-density lipoprotein (HDL) plays an important role in the transport and metabolism of cholesterol. Mimics of HDL are being explored as potentially powerful therapeutic agents for removing excess cholesterol from arterial plaques. Gold nanoparticles (AuNPs) functionalized with apolipoprotein A-I and with the lipids 1,2-dipalmitoyl-sn-glycero-3-phosphocholine and 1,2-dipalmitoyl-sn-glycero-3-phosphoethanolamine-N-[3-(2-pyridyldithio)propionate] have been demonstrated to be robust acceptors of cellular cholesterol. However, detailed structural information about this functionalized HDL AuNP is still lacking. Here in this study, we have used X-ray photoelectron spectroscopy and lecithin/cholesterol acyltransferase activation experiments together with coarse-grained and all-atom molecular dynamics simulations to model the structure and cholesterol uptake properties of the HDL AuNP construct. By simulating different apolipoprotein-loaded AuNPs, we find that lipids are oriented differently in regions with and without apoA-I. We also show that in this functionalized HDL AuNP, the distribution of cholesteryl ester maintains a reverse concentration gradient that is similar to the gradient found in native HDL.

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

Effects of polarizability on the structural and thermodynamics properties of [C{sub n}mim][Gly] ionic liquids (n = 1–4) using EEM/MM molecular dynamic simulations

An extended electronegativity equalization method/molecular mechanics (EEM/MM) model for ionic liquids is used to investigate the structures and properties of 1-alkyl-3-methylimidazolium glycine ionic liquids [C{sub n}mim][Gly] (n = 1–4) with alkyl substituents of different lengths. The EEM/MM model describes the electrostatic interactions of atoms and their changes in different ambient environments. This property is the most outstanding characteristic of the model. EEM parameters (i.e., valence electronegativities and valence hardness parameters) are calibrated using linear regression and least-squares methods, which can accurately predict the gas-phase properties of [C{sub n}mim]{sup +}, [Gly]{sup −}, and [C{sub n}mim][Gly] ion pairs. We utilize the EEM/MM force field to systematically investigate the effects of polarizability on the accuracy of [C{sub n}mim][Gly] properties predicted through the molecular dynamic simulations. EEM/MM explicitly describes the atom-based polarizability of [C{sub n}mim][Gly]; thus, the densities, enthalpies of vaporization, self-diffusion coefficients, and conductivities of the [C{sub n}mim][Gly] are consistent with the experimental values. The calculated radial distribution functions provide a mechanistic understanding of the effects of polarizability on ionic aggregations in amino acid ionic liquids. The effects of alkyl chain length on the diffusion coefficient and conductivity are also discussed.

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

Cluster dynamics simulations of tritium and helium diffusion in lithium ceramics

Tritium (T) and He diffusion in LiAlO 2 and LiAl 5 O 8 phases influences the performance of tritium producing burnable absorber rods (TPBARs) by affecting the gas release, swelling and thermal conductivity of Li-bearing ceramic pellets. Frenkel pair defects and clusters created by irradiation can attract T and He interstitials and form clusters of the type He i x Li, He i x Al, He i x O, T i x Li, T i x Al, and T i x O, 1 ≤ x ≤ 4 in a Li, Al or O vacancy site (notation denotes x He or T atoms in a 1 Li, 1 Al or 1 O vacant site). The concentration and mobility of each of these clusters collectively contribute to the diffusion of the He and T gases in LiAlO 2 and LiAl 5 O 8 . Here, in this work, free energy cluster dynamics simulations implemented in the Centipede code, are used to obtain the concentration and diffusivities of these clusters which are then used to calculate the total diffusivity of T and He gases in LiAlO 2 and LiAl 5 O 8 . The results show that diffusivity of T is at least one order of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 whereas He diffusion is 2–13 orders of magnitude higher in LiAlO 2 as compared to that in LiAl 5 O 8 . There is a higher concentration of highly diffusive species (T interstitials and T i 03 Li for the case of tritium and He i 01 Li, He i 02 Li, and He i 03 Li for the case of He) in LiAlO 2 than in LiAl 5 O 8 which increase the total diffusion of T and He in LiAlO 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2020

Molecular Dynamics Simulation of Atomic Interactions in the Vancomycin Binding Site

Vancomycin is a glycopeptide antibiotic produced by Amycolaptopsis orientalis used to treat serious infections by Grampositive pathogens including methicillin-resistant Staphylococcus aureus. Vancomycin inhibits cell wall biosynthesis by targeting lipid II, which is the membrane-bound peptidoglycan precursor. The heptapeptide aglycon structure of vancomycin binds to the D-Ala-D-Ala of the pentapeptide stem structure in lipid II. The third residue of vancomycin aglycon is asparagine, which is not directly involved in the dipeptide binding. Nonetheless, asparagine plays a crucial role in substrate recognition, as the vancomycin analogue with asparagine substituted by aspartic acid (V D ) shows a reduction in antibacterial activities. To characterize the function of asparagine, binding of vancomycin and its aspartic-acid-substituted analogue V D to L-Lys-D-Ala-D-Ala and L-Lys-D-Ala-D-Lac was investigated using molecular dynamic simulations. Binding interactions were analyzed using root-mean-square deviation (RMSD), two-dimensional (2D) contour plots, hydrogen bond analysis, and free energy calculations of the complexes. The analysis shows that the aspartate substitution introduced a negative charge to the binding cleft of V D , which altered the aglycon conformation that minimized the repulsive lone pair interaction in the binding of a depsipeptide. Our findings provide new insight for the development of novel glycopeptide antibiotics against the emerging vancomycin-resistant pathogens by chemical modification at the third residue in vancomycin to improve its binding affinity to the D-Ala-D-Lac-terminated peptidoglycan in lipid II found in vancomycin-resistant enterococci and vancomycinresistant S. aureus.

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

EXPLORING THE STRUCTURE OF GLASS-FORMING LIQUIDS USING HIGH ENERGY X-RAY DIFFRACTION, CONTAINERLESS METHODOLOGY AND STATE-OF-THE-ART MOLECULAR DYNAMICS SIMULATION

High energy X-ray diffraction can be combined with containerless techniques to provide information on the atomic arrangements in glass-forming liquids in stable and metastable regimes. The high incident energies provide bulk diffraction data to high values of scattering vector which enables significantly more robust analysis of the local and medium-range order that influences important physical properties such as viscosity and crystal nucleation. These combined techniques have been applied to a range of oxide liquids. In this contribution we illustrate addition of further dimensions to phase space by controlling the partial pressure of oxygen that permits the study liquids containing iron. The advantages of rapid data acquisition are also demonstrated in a study of tellurite glass-forming systems where a transition from ergodic to non-ergodic regimes in the deeply supercooled liquid is shown. Finally we demonstrate how descriptions of the liquid structure can be developed by combining HEXRD with molecular dynamics simulations.

containerless techniques↗
DOE OSTI2020

Efficient perturbation-tracking method for directly probing the spectral phonon properties from molecular dynamics simulations

Existing methods for directly extracting the spectral phonon properties from molecular dynamics (MD) simulations, like the normal mode analysis (NMA) and spectral energy density analysis, all require a very long simulation time to produce reliable results with good convergence. So far, these methods are mainly applied in studies using small systems and with empirical potentials, as the heavy computational load has greatly hindered their further applications. Here we propose a perturbation-tracking (PT) method for directly probing the mode-wise phonon anharmonic frequencies and lifetimes. We show that results obtained from our method are in excellent agreement with those from the conventional NMA approach, using Si as the model material system. Comparing with the NMA approach, the PT method offers a greater accuracy and significant improvement of efficiency. It takes an average of two orders of magnitude and up to three orders of magnitude less simulation time to obtain the same lifetime result of a phonon mode with intermediate to high accuracy. Meanwhile, our method preserves all the dynamics of probed phonon mode from a particular state, which means it is capable of studying the transient thermal transport processes in a nonequilibrium system. Besides the exceptional efficiency, our method also comes with freedom to choose to probe only those modes of interest. This makes it ideal for use with large systems and in computationally demanding applications, such as ab initio MD simulations. Moreover, the PT method we propose here is very straightforward and easy to implement.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Development of a machine-learning-based ionic-force correction model for quantum molecular dynamic simulations of warm dense matter

In this work Δ learning is used to map orbital-free density functional theory (OF-DFT) ionic forces to the corresponding Kohn-Sham (KS) DFT ionic forces. The development of the approximate force difference in terms of the ion positions is constructed and serves as a stand in for the ground truth force difference. Descriptor vectors for ion configurations are constructed using all distance between ions in conjunction with an indexing based on a nearest neighbor ranking. It is demonstrated that such a scheme of descriptors can uniquely describe an ionic configuration up to a rotation and reflection when no ambiguity in the nearest neighbor ranking exists. How to handle the case when an ambiguity exists in the nearest neighbor ranking is discussed. As a proof of principle, the model is trained and tested on warm dense hydrogen at temperatures between 1 and 15 eV. Once tested, the model was used to perform molecular dynamic simulations of warm dense hydrogen. Furthermore, the resulting energies and pressures are within 1% and 2% of their respective target KS values.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Self-Assembly of Silica Nanoparticles at Water–Hydrocarbon Interfaces: Insights from In Operando Small-Angle X-ray Scattering Measurements and Molecular Dynamics Simulations

To achieve tunable controls on the interactions of siliceous materials in subsurface environments bearing water and hydrocarbons, it is essential to determine the influence of hydrocarbon–water interfaces on the self-assembly of the silica particles. The hydrophilic silica particle have a tendency to aggregate on the water front of the interface. Self-assembly of silica nanoparticles proceeds via the migration of these nanoparticles to the water–hydrocarbon interface followed by aggregation at the interface. Fractal-like morphologies of assembled silica nanoparticles at water–hydrocarbon interfaces are observed. Rapid assembly of the hydrophilic silica nanoparticles at water–hydrocarbon interfaces corresponds to an overall reduction in the surface tension of water–toluene and water–heptane systems. These studies demonstrate the silica aggregation that is the precursor of silica polymerization, or nucleation and growth is influenced by the presence of hydrocarbons in subsurface geologic environments. These insights were derived from in operando ultrasmall and small-angle X-ray scattering (USAXS/SAXS) measurements, cryo-scanning electron microscopy (Cryo-SEM) imaging, and classical molecular dynamics simulations. As a result, these studies are intended to inform current and future efforts aimed at tuning silica reactivity in subsurface geologic environments to enhance permeability, the design and use of silica-based proppants to enhance fractures, and the development of effective strategies to control silica-based scaling behavior in subsurface reservoirs.

42 ENGINEERING↗
NASA NTRS1991

System dynamic simulation of precision segmented reflector

A joint effort was undertaken on a Precision Segmented Reflector (PSR) Project. The missions in which the PSR is to be used will use large (up to 20 m in diameter) telescopes. The essential requirement for the telescopes is that the reflector surface of the primary mirror must be made extremely precise to allow no more than a few microns of errors and, additionally, this high surface precision must be maintained when the telescope is subjected to on-orbital mechanical and thermal disturbances. Based on the mass, size, and stability considerations, reflector surface formed by segmented, probably actively or passively controlled, composite panels are regarded as most suitable for future space based astronomical telescope applications. In addition to the design and fabrication of composite panels with a surface error of less than 3 microns RMS, PSR also develops related reflector structures, materials, control, and sensing technologies. As part of the planning effort for PSR Technology Demonstration, a system model which couples the reflector, consisting of panels, support truss and actuators, and the optical bench was assembled for dynamic simulations. Random vibration analyses using seismic data obtained from actual measurements at the test site designated for PSR Technology Demonstration are described.

Shih, Choon-Foo↗