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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

ReaxFF molecular dynamics simulations of electrolyte–water systems at supercritical temperature

We have performed ReaxFF molecular dynamics simulations of alkali metal–chlorine pairs in different water densities at supercritical temperature (700 K) to elucidate the structural and dynamical properties of the system. The radial distribution function and the angular distribution function explain the inter-ionic structural and orientational arrangements of atoms during the simulation. The coordination number of water molecules in the solvation shell of ions increases with an increase in the radius of ions. We find that the self-diffusion coefficient of metal ions increases with a decrease in density under supercritical conditions due to the formation of voids within the system. The hydrogen bond dynamics has been interpreted by the residence time distribution of various ions, which shows Li+ having the highest water retaining capability. The void distribution within the system has been analyzed by using the Voronoi polyhedra algorithm providing an estimation of void formation within the system at high temperatures. We observe the formation of salt clusters of Na+ and K+ at low densities due to the loss of dielectric constants of ions. The diffusion of ions gets altered dramatically due to the formation of voids and nucleation of ions in the system.

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

RDX Compression, α→ γ Phase Transition, and Shock Hugoniot Calculations from Density-Functional-Theory-Based Molecular Dynamics Simulations

Prediction of the density and lattice compression properties of the α and γ phases of the hexahydro-1,3,5-trinitro-1,3,5-s-triazine (RDX) crystal and of the low-pressure α → γ phase transition upon pressure increase are general tests used to assess the accuracy of density-functional-theory- (DFT-) based computational methods and to identify the essential parameters that govern the behavior of this high-energy-density material under extreme conditions. The majority of previous DFT studies have analyzed such issues under static optimization conditions by neglecting the corresponding temperature effects. In this paper, we extend previous investigations and analyze the performance of dispersion-corrected density functional theory to predict the compression of RDX in the pressure range of 0–9 GPa and the corresponding α → γ phase transition under realistic temperature and pressure conditions. We demonstrate that, by using static dispersion-corrected density functional theory calculations, direct interconversion between the α and γ phases upon compression is not observed. This limitation can be addressed by using isobaric–isothermal molecular dynamic simulations in conjunction with DFT-D2-calculated potentials, an approach that is shown to provide an accurate description of both the crystallographic RDX lattice parameters and the dynamical effects associated with the α→ γ phase transformation. An even more comprehensive and demanding analysis was done by predicting the corresponding shock Hugoniot curve of RDX in the pressure range of 0–9 GPa. It was found that the theoretical results reproduce reasonably well the available experimental Hugoniot shock data for both the α and γ phases. Finally, the results obtained demonstrate that a satisfactory prediction of the shock properties in high-energy-density materials undergoing crystallographic and configurational transformations is possible through the combined use of molecular dynamics simulations in the isobaric–isothermal ensemble with dispersion-corrected density functional theory methods.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS2020

Multibody Based Digital Astronaut Dynamics Simulation

BACKGROUND: This study provides the Software, Robotics, & Simulation Division at the NASA Johnson Space Center with a verification tool for multibody dynamics simulation requiring human motion. The motivation stems from current studies of several Vibration Isolation & Stabilization(VIS)system designs that attenuate the moments and forces which would be transmitted to a spacecraft during an exercise. A multibody dynamics model for a proposed VIS was available previously[1], therefore modeling of the VIS was not needed for this work. The interest here is in creating the multibody dynamics model of an astronaut in motion which may be utilized independently or while attached to a mechanism. An existing simulation [2] that utilizes OpenSim [3,4] and an in-house multibody dynamics package (MBDyn) [5] is used in order to verify the astronaut model. The main advantage this model will have over the existing simulation is that everything will be processed in one tool. METHODS AND RESULTS: Creating the simulation required; estimation of Body Segment Inertial Parameters (BSIP),a multibody model of the human-VIS system, joint acceleration profiles, and input files for MBDyn, which is used for this analysis. The scaling factors provided by Dumas et al. [6] are utilized in estimating the BSIP. Anthropometric data are used for estimating these parameters, the Anthropometric Survey of US Army Personnel (ANSUR II) [7] was the source. The astronaut model consists of 15 bodies, 14 joints and 32 degrees of freedom, with the dynamics topology generated using the center of mass locations and anthropometric data. MBDyn has an option for prescribed joint motion (PJM), which requires joint acceleration data as input. The joint angle data is first obtained from a motion capture system and then processed through code that has been created to generate approximate joint acceleration profiles. The topology tree for the astronaut model begins at the right foot up to the pelvis where there is one branch for going down the left leg and another for the torso. The torso branch leads to branches for the arms and a leaf body for the head/neck segment. For attachment to the VIS, the heel of the right foot is connected to the VIS platform through a fixed joint, resembling a foot restraint. The left foot does not attach to the platform in order to prevent a system with a closed loop. Topology and symmetry of the astronaut model were verified through kinematic analysis. Further verification of the forces and moments transmitted to the VIS were verified against the existing simulation. There was a satisfactory level of agreement when testing a simple motion, for example, rocking back and forth. Full exercise motions are to be tested soon. The main outcome has been a novel application of MBDyn for biomechanics modeling that is now available for dynamic simulations involving human motion. The estimation of BSIP was another useful result of this study, requiring only 15 inputs for generating mass properties of a theoretical astronaut model. Expansion on this work is possible by going through an alternative route in obtaining the joint motion data. Instead of high-tech and often expensive motion capture systems, an individual may watch videos with high focus and at a slow motion for each individual segment in order to determine the initial and final time and angle for that specific degree of freedom. Synthetic trajectories may also be created if there is no video reference available.

F N Matari↗
NASA NTRS1995

Molecular dynamics simulation of a piston driven shock wave in a hard sphere gas

Molecular dynamics simulation is used to study the piston driven shock wave at Mach 1.5, 3, and 10. A shock tube, whose shape is a circular cylinder, is filled with hard sphere molecules having a Maxwellian thermal velocity distribution and zero mean velocity. The piston moves and a shock wave is generated. All collisions are specular, including those between the molecules and the computational boundaries, so that the shock development is entirely causal, with no imposed statistics. The structure of the generated shock is examined in detail, and the wave speed; profiles of density, velocity, and temperature; and shock thickness are determined. The results are compared with published results of other methods, especially the direct simulation Monte-Carlo method. Property profiles are similar to those generated by direct simulation Monte-Carlo method. The shock wave thicknesses are smaller than the direct simulation Monte-Carlo results, but larger than those of the other methods. Simulation of a shock wave, which is one-dimensional, is a severe test of the molecular dynamics method, which is always three-dimensional. A major challenge of the thesis is to examine the capability of the molecular dynamics methods by choosing a difficult task.

Woo, Myeung-Jouh↗
NASA NTRS1989

Efficient dynamic simulation for multiple chain robotic mechanisms

An efficient O(mN) algorithm for dynamic simulation of simple closed-chain robotic mechanisms is presented, where m is the number of chains, and N is the number of degrees of freedom for each chain. It is based on computation of the operational space inertia matrix (6 x 6) for each chain as seen by the body, load, or object. Also, computation of the chain dynamics, when opened at one end, is required, and the most efficient algorithm is used for this purpose. Parallel implementation of the dynamics for each chain results in an O(N) + O(log sub 2 m+1) algorithm.

Lilly, Kathryn W.↗
DOE OSTI2017

Ab initio molecular dynamics simulations of low energy recoil events in MgO

In this paper, low-energy recoil events in MgO are studied using ab initio molecular dynamics simulations to reveal the dynamic displacement processes and final defect configurations. Threshold displacement energies, E d , are obtained for Mg and O along three low-index crystallographic directions, [100], [110], and [111]. The minimum values for E d are found along the [110] direction consisting of the same element, either Mg or O atoms. Minimum threshold values of 29.5 eV for Mg and 25.5 eV for O, respectively, are suggested from the calculations. For other directions, the threshold energies are considerably higher, 65.5 and 150.0 eV for O along [111] and [100], and 122.5 eV for Mg along both [111] and [100] directions, respectively. These results show that the recoil events in MgO are partial-charge transfer assisted processes where the charge transfer plays an important role. Finally, there is a similar trend found in other oxide materials, where the threshold displacement energy correlates linearly with the peak partial-charge transfer, suggesting this behavior might be universal in ceramic oxides.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Effects of Surface Orientation and Termination Plane on Glass‐to‐Crystal Transformation of Lithium Disilicate by Molecular Dynamics Simulations

Glass‐to‐crystal transformation of lithium disilicate is studied using molecular dynamics simulations using an effective partial charge potential. The structural evolution of the interface between glassy and crystalline lithium disilicate is analyzed to simulate crystallization of glass on pre‐existing crystal seeds. Besides previously used atomic number density, the distribution of Q n species (Si tetrahedra with n bridging oxygen) is shown to be an effective parameter for following this transformation quantitatively. The early stages of crystal growth are significantly affected by the orientation and termination of the surface of adjacent crystal, as indicated by calculated atomic density, partial ordering, atomic segregation, and an increase in Q 3 concentration. In particular, under‐coordinated Si within the outer crystal layer is found to be most effective in transforming the amorphous structure toward crystallinity. The increase in Q 3 in the glass close to interface region most clearly shows the initial stage of lithium disilicate crystal growth.

Sun, Wei↗
DOE OSTI2018

Model Evaluation Report for LNG Dispersion on Fire Dynamics Simulator Version 6.5.3

This report is an assessment of the computational fluid dynamics (CFD) code Fire Dynamics Simulator (FDS), version 6.5.3, using the Model Evaluation Protocol (MEP) for Liquefied Natural Gas (LNG). The MEP consists of model validation and a scientific assessment that encompasses model information and model verification activities. Model validation entails comparison against various field and wind-tunnel trials for LNG and other dense gases with and without obstructions. The results indicate that FDS is generally over-predictive for most of the field scale dense gas releases and cases involving obstructions without requiring a safety factor of 2, however, since there are uncertainties in extrapolating to accident scales a safety factor of 2 (1/2 LFL criteria) may be appropriate. The results of this scientific assessment and prior verification efforts demonstrate the soundness of the numerical techniques and robustness of FDS. The validation results indicate that FDS is suitable for modeling dense gas dispersion with and without obstructions.

03 NATURAL GAS↗
NASA NTRS1983

A large-signal dynamic simulation for the series resonant converter

A simple nonlinear discrete-time dynamic model for the series resonant dc-dc converter is derived using approximations appropriate to most power converters. This model is useful for the dynamic simulation of a series resonant converter using only a desktop calculator. The model is compared with a laboratory converter for a large transient event.

King, R. J.↗
NASA NTRS1995

Program For Simulating Dynamics Of Aerospace Vehicles

SORT (Simulation and Optimization of Rocket Trajectories) is general-purpose three-degree-of-freedom with three axis static moment balance simulation of flight dynamics of arbitrary aerospace vehicle. Modular structure facilitates application to variety of trajectory-analysis problems. Contains math model of aerodynamics completely generalized. Computes both longitudinal and lateral forces and moments. In addition to fore-body coefficients, computes longitudinal base effect aerodynamic forces and moments. Simplified ballistic-coefficient model also available for analysis of ballistic entry. Written using ANSI FORTRAN 77.

Berning, M. J.↗
DOE OSTI2017

Diffusion of Selective Penetrants in Interfacially Modified Block Copolymers from Molecular Dynamics Simulations

In order to show the influence of the interface on structure and dynamics of microphase separated polymer systems, we study interfacially modified AB block copolymers with small molecule penetrants. The polymers have a random midblock or tapered midblock whose composition varies from pure A to pure B (or from pure B to pure A for an inverse taper) between two pure blocks of A and B. We perform simple coarse-grained molecular dynamics simulations of symmetric polymers that form lamellae. With normal tapering, both polymer and penetrant diffusion parallel to the lamellae increases as taper length increases. Inverse tapered polymers exist in different conformational states (e.g. stretched versus folded back and forth across the interface) with different dynamic behavior, leading to non-monotonic trends in their diffusion. However, the local mixing of monomers (rather than polymer conformation) appears to be the most important factor in determining penetrant diffusion.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Semi-Empirical Interatomic Potential for Large-Scale Molecular Dynamics Simulation of Metal-Oxide Systems

Classical molecular dynamics (MD) simulation can be applied to systems containing billions of atoms during times up to microseconds. However, utilization of a large-scale MD simulation requires reliable but computationally cheap interatomic potentials. In the case of metallic systems, embedded atom method (EAM) and Finnis-Sinclair (FS) potentials are the reasonable choices, but their development for multicomponent alloys is a challenge. Addition of oxygen atoms makes the problem of developing semi-empirical potentials even more difficult because the long-range Coulomb interaction cannot be neglected, and different atoms have different charges which vary during MD relaxation such that a charge equilibration procedure should be applied every MD step. As results researchers tend to avoid simulating metal-oxide systems. In this talk, a new Finnis-Sinclair (FS) Ni-Co-Cr potential will be presented. Special attention will be paid to reproducing of the element melting temperatures and clustering. Next, a charge transfer ionic potential (CTIP) developed to simulate the NiO properties will be presented. This potential utilizes the Ni FS potential and accounts for the Coulomb interaction in the vicinity of oxygen atoms. Fitting of the CTIP parameters and FS potential functions describing the interaction with oxygen atoms to the available experimental and ab initio data will be discussed. Results of MD simulation of interaction of NiO particles with dislocations in Ni will be shown.

molecular dynamics simulation↗
NASA NTRS1985

Software life cycle dynamic simulation model: The organizational performance submodel

The submodel structure of a software life cycle dynamic simulation model is described. The software process is divided into seven phases, each with product, staff, and funding flows. The model is subdivided into an organizational response submodel, a management submodel, a management influence interface, and a model analyst interface. The concentration here is on the organizational response model, which simulates the performance characteristics of a software development subject to external and internal influences. These influences emanate from two sources: the model analyst interface, which configures the model to simulate the response of an implementing organization subject to its own internal influences, and the management submodel that exerts external dynamic control over the production process. A complete characterization is given of the organizational response submodel in the form of parameterized differential equations governing product, staffing, and funding levels. The parameter values and functions are allocated to the two interfaces.

Tausworthe, Robert C.↗
DOE OSTI2019

Scalable Reactive Molecular Dynamics Simulations for Computational Synthesis

Reactive molecular dynamics (MD) simulation is a powerful research tool for describing chemical reactions. We eliminate the speed-limiting charge iteration in MD with a novel extended-Lagrangian scheme. The extended-Lagrangian reactive MD (XRMD) code drastically improves energy conservation while substantially reducing time-to-solution. Furthermore, we introduce a new polarizable charge equilibration (PQEq) model to accurately predict atomic charges and polarization. The XRMD code based on hybrid message passing+multithreading achieves a weak-scaling parallel efficiency of 0.977 on 786 432 IBM Blue Gene/Q cores for a 67.6 billion-atom system. The performance is portable to the second-generation Intel Xeon Phi, Knights Landing. Blue Gene/Q simulations for the computational synthesis of materials via novel exfoliation mechanisms for synthesizing atomically thin transition metal dichalcogenide layers will dominate nanomaterials science in this century.

Simulation, Modeling, Performance, and Visualizati↗
DOE OSTI2025

Machine learning-accelerated path integral molecular dynamics simulations of reactive organic electrolytes

Hydrogen bonded electrolytes that exhibit accelerated proton transport via sequential reactive hops have drawn interest for their promise in clean energy applications. Molecular dynamics simulations of these electrolytes offer the opportunity to uncover microscopic mechanistic details that could be used to design and tune the properties of candidate electrolyte technologies. However, accurately modeling the proton transfer reactions and transport properties that give rise to high charge conductivites in these electrolytes proves computationally challenging because of the need to perform lengthy condensed phase simulations, treating both the electronic and nuclear degrees of freedom quantum mechanically. In this paper, we demonstrate that such a modeling task can be efficiently achieved with the use of density functional theory (DFT)-trained machine learning potentials (MLP) to accelerate path integral molecular dynamics (PIMD) simulations. We highlight the practical utility of this approach by using it to benchmark how closely PIMD simulations employing different DFT exchange–correlation functionals reproduce the composition-dependent densities, diffusion coefficients, and electrical conductivities of mixtures consisting of imidazole and levulinic acid. Even with the speedup afforded by our MLPs, PIMD simulations remain quite expensive. Furthermore, in order to render PIMD more computationally tractable, we introduce and benchmark the accuracy of a ring polymer contraction approach that leverages a computationally efficient short-range MLP to accelerate our PIMD simulations by an additional factor of four.

Chemical bonding↗
DOE OSTI2025

Acceleration of Power System Dynamic Simulations Using a Deep Equilibrium Layer and Neural ODE Surrogate

The dominant paradigm for power system dynamic simulation is to build system-level simulations by combining physics-based models of individual components. The sheer size of the system along with the rapid integration of inverter-based resources exacerbates the computational burden of running time domain simulations. Here, in this paper, we propose a data-driven surrogate model based on implicit machine learningspecifically deep equilibrium layers and neural ordinary differential equationsto learn a reduced order model of a portion of the full underlying system. The data-driven surrogate achieves similar accuracy and reduction in simulation time compared to a physics-based surrogate, without the constraint of requiring detailed knowledge of the underlying dynamic models. This work also establishes key requirements needed to integrate the surrogate into existing simulation workflows; the proposed surrogate is initialized to a steady state operating point that matches the power flow solution by design.

Neural ordinary differential equations↗
DOE OSTI2018

Thermodynamic analysis of the Livermore molecular-dynamics simulations of dislocation-mediated plasticity

Results of recent large-scale molecular dynamics simulations of dislocation-mediated solid plasticity are compared with predictions of the statistical thermodynamic theory of these phenomena. These computational and theoretical analyses are in substantial agreement with each other in both their descriptions of strain-rate-dependent steady plastic flows and of transient stress peaks associated with initially small densities of dislocations. Lastly, the comparisons between the numerical simulations and basic theory reveal inconsistencies in some conventional phenomenological descriptions of solid plasticity.

97 MATHEMATICS AND COMPUTING↗