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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

On computing stress in polymer systems involving multi-body potentials from molecular dynamics simulation

Hardy stress definition has been restricted to pair potentials and embedded-atom method potentials due to the basic assumptions in the derivation of a symmetric microscopic stress tensor. Force decomposition required in the Hardy stress expression becomes obscure for multi-body potentials. In this work, we demonstrate the invariance of the Hardy stress expression for a polymer system modeled with multi-body interatomic potentials including up to four atoms interaction, by applying central force decomposition of the atomic force. The balance of momentum has been demonstrated to be valid theoretically and tested under various numerical simulation conditions. The validity of momentum conservation justifies the extension of Hardy stress expression to multi-body potential systems. Computed Hardy stress has been observed to converge to the virial stress of the system with increasing spatial averaging volume. This work provides a feasible and reliable linkage between the atomistic and continuum scales for multi-body potential systems.

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

Towards computational polar-topotronics: Multiscale neural-network quantum molecular dynamics simulations of polar vortex states in SrTiO3/PbTiO3 nanowires

Recent discoveries of polar topological structures ( e.g ., skyrmions and merons) in ferroelectric/paraelectric heterostructures have opened a new field of polar topotronics. However, how complex interplay of photoexcitation, electric field and mechanical strain controls these topological structures remains elusive. To address this challenge, we have developed a computational approach at the nexus of machine learning and first-principles simulations. Our multiscale neural-network quantum molecular dynamics molecular mechanics approach achieves orders-of-magnitude faster computation, while maintaining quantum-mechanical accuracy for atoms within the region of interest. This approach has enabled us to investigate the dynamics of vortex states formed in PbTiO 3 nanowires embedded in SrTiO 3 . We find topological switching of these vortex states to topologically trivial, uniformly polarized states using electric field and trivial domain-wall states using shear strain. These results, along with our earlier results on optical control of polar topology, suggest an exciting new avenue toward opto-electro-mechanical control of ultrafast, ultralow-power polar topotronic devices.

Linker, Thomas↗
DOE OSTI2021

Accelerated molecular dynamics simulations of dislocation climb in nickel

We report the mechanical behavior of materials operating under high temperatures is strongly influenced by creep mechanisms such as dislocation climb, which is controlled by the diffusion of vacancies. However, atomistic simulations of these mechanisms have traditionally been impractical due to the long time scales required. To overcome these time scale challenges, we use Parallel Trajectory Splicing (ParSplice), an accelerated molecular dynamics method, to simulate dislocation climb in nickel. We focus on modeling the activity of a vacancy near a jog on an edge dislocation in order to observe vacancy pipe diffusion and vacancy absorption at the jog. From rigorously constructed trajectories encompassing more than 2000 vacancy absorption events over a simulation time of more than 4μs at 900 K, a comprehensive sampling of available atomistic mechanisms is collated and analyzed further with molecular statics calculations. We estimate average rates for pipe diffusion and vacancy absorption into the jog using data from the dynamic and static calculations, finding very good agreement. Our results strongly suggest that the dominant mechanism for vacancy absorption by jogs is via biased diffusion to the dislocation core followed by fast pipe diffusion to the jog.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Communication: Quantum molecular dynamics simulation of liquid para-hydrogen by nuclear and electron wave packet approach

Liquid para-hydrogen (p-H{sub 2}) is a typical quantum liquid which exhibits strong nuclear quantum effects (NQEs) and thus anomalous static and dynamic properties. We propose a real-time simulation method of wave packet (WP) molecular dynamics (MD) based on non-empirical intra- and inter-molecular interactions of non-spherical hydrogen molecules, and apply it to condensed-phase p-H{sub 2}. The NQEs, such as WP delocalization and zero-point energy, are taken into account without perturbative expansion of prepared model potential functions but with explicit interactions between nuclear and electron WPs. The developed MD simulation for 100 ps with 1200 hydrogen molecules is realized at feasible computational cost, by which basic experimental properties of p-H{sub 2} liquid such as radial distribution functions, self-diffusion coefficients, and shear viscosities are all well reproduced.

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

Modification of a force field for molecular dynamics simulations of silicon etching by chlorine atoms

Here, a modified classical molecular dynamics (MD) force field is presented for the interaction between silicon (Si) and chlorine (Cl). The original version of the force field is shown to significantly overestimate the probability of Si etching by thermal Cl atoms. However, the modified force field corrects this problem and results in generally good agreement with experimental data. Further, it is shown that while the modification of the force field improves the prediction of Si spontaneous etching with Cl atoms, it does not degrade predictions of atomic-layer etching of Si with Cl[Formula: see text] molecules.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1989

Development of efficient computer program for dynamic simulation of telerobotic manipulation

Research in robot control has generated interest in computationally efficient forms of dynamic equations for multi-body systems. For a simply connected open-loop linkage, dynamic equations arranged in recursive form were found to be particularly efficient. A general computer program capable of simulating an open-loop manipulator with arbitrary number of links has been developed based on an efficient recursive form of Kane's dynamic equations. Also included in the program is some of the important dynamics of the joint drive system, i.e., the rotational effect of the motor rotors. Further efficiency is achieved by the use of symbolic manipulation program to generate the FORTRAN simulation program tailored for a specific manipulator based on the parameter values given. The formulations and the validation of the program are described, and some results are shown.

Chen, J.↗
NASA NTRS1988

Obstacles to high fidelity multibody dynamics simulation

The authors take a wide ranging look at issues involved in the dynamic modeling of complex, multibodied orbiting space systems. The nature of the multibody problem is reviewed and capabilities and limitations of two major codes (DISCOS, TREETOPS) are assessed. Problem areas that limit accuracy or contribute to uncertainty in the modeling, simulation process are identified. As a specific case, the significance of including, or of not including, second-degree geometric nonlinearity in the elastic displacement field (foreshortening) is evaluated. Conclusions and recommendations are made concerning the direction future development should take to achieve higher-fidelity and more computationally efficient multibody software solutions.

Lips, K. W.↗
DOE OSTI2017

Cluster analysis of accelerated molecular dynamics simulations: A case study of the decahedron to icosahedron transition in Pt nanoparticles

Modern molecular-dynamics-based techniques are extremely powerful to investigate the dynamical evolution of materials. With the increase in sophistication of the simulation techniques and the ubiquity of massively parallel computing platforms, atomistic simulations now generate very large amounts of data, which has to be carefully analyzed in order to reveal key features of the underlying trajectories, including the nature and characteristics of the relevant reaction pathways. We show that clustering algorithms, such as the Perron Cluster Cluster Analysis, can provide reduced representations that greatly facilitate the interpretation of complex trajectories. To illustrate this point, clustering tools are used to identify the key kinetic steps in complex Accelerated Molecular Dynamics trajectories exhibiting shape uctuations in Pt nanoclusters. This analysis provides an easily-interpretable coarse representation of the reaction pathways in terms of a handful of clusters, in contrast to the raw trajectory that contains thousands of unique states and tens of thousands of transitions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Water Dynamics in Nanoporous Silica: Ultrafast Vibrational Spectroscopy and Molecular Dynamics Simulations

Nanoporous silica materials are important in catalysis, energy, and materials applications in which water is an essential component. System performance is intimately connected to the water dynamics occurring in the confined environment. However, the dynamics and associated structures of water in nanoporous silica have proven challenging to measure and predict. Here, confined water dynamics are examined via the ultrafast infrared spectroscopy of selenocyanate (SeCN – ) dissolved in the hydrated ~2.4 nm silica mesopores of MCM41. Polarization selective pump–probe and two-dimensional infrared measurements on the CN stretching mode of SeCN – are used to probe the effect of confinement on orientational relaxation and spectral diffusion dynamics. Here, the dynamics of SeCN – provide information on the surrounding water hydrogen bond dynamics. The long CN stretch lifetime (~36 ps), relative to the water hydroxyl stretch (<2 ps), significantly extends the time scales that can be accessed. Complete orientational relaxation (C 2 (t), orientational correlation function) and spectral diffusion (C ω (t), frequency–frequency correlation function) dynamics are presented and compared to the simulated time correlation functions in a model silica pore of the same size. A slow decay component not present in the bulk liquid is observed in both experiments, indicating that the hydrogen bond dynamics are significantly altered by confinement. The simulations reveal a qualitative difference in the functional dependence of C 2 (t;d) and C ω (t;d) on d, the distance from the interface. The former becomes exponentially faster with distance while the latter makes an abrupt transition from slower to faster dynamics midway between the surface and pore center, d ≅ 6 Å.

14 SOLAR ENERGY↗
DOE OSTI2016

Inward open characterization of EmrD transporter with molecular dynamics simulation

EmrD is a member of the multidrug resistance exporter family. Up to now, little is known about the structural dynamics that underline the function of the EmrD protein in inward-facing open state and how the EmrD transits from an occluded state to an inward open state. For the first time the article applied the AT simulation to investigate the membrane transporter protein EmrD, and described the dynamic features of the whole protein, the domain, the helices, and the amino acid residues during an inward-open process from its occluded state. The gradual inward-open process is different from the current model of rigid-body domain motion in alternating-access mechanism. Simulation results show that the EmrD inward-open conformational fluctuation propagates from a C-terminal domain to an N-terminal domain via the linker region during the transition from its occluded state. The conformational fluctuation of the C-terminal domain is larger than that of the N-terminal domain. In addition, it is observed that the helices exposed to the surrounding membrane show a higher level of flexibility than the other regions, and the protonated E227 plays a key role in the transition from the occluded to the open state. -- Highlights: •This study described the dynamic features of the whole EmrD protein, during an inward-open process from its occluded state. •The EmrD inward-open conformational fluctuation propagates from a C-terminal domain to an N-terminal domain via the linker region during the transition from its occluded state. •The conformational fluctuation of the C-terminal domain is larger than that of the N-terminal domain. •The protonated E227 plays a key role in the transition from the occluded to the open state.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1974

The effect of simulator dynamics on pilot response

The effects of visual display dynamics on the altitude tracking performance of a subject in a fixed base flight simulator are considered. The subject, flying the linearized longitudinal equations of motion, attempts to maintain the same altitude as two airplanes positioned three hundred feet ahead, as in level formation flying. The horizon together with the two leading aircraft are represented symbolically on a CRT display. The subject's aircraft is disturbed by atmospheric turbulence. The data indicate a relationship between the bandwidth of the display dynamics and the short period characteristics of the simulated airplane. For an airplane with a relatively fast pitch response the presence of altitude display dynamics, with a bandwidth as high as five times the short period natural frequency, causes significant degradation of altitude tracking performance.

Weener, E. F.↗
NASA NTRSDate not supplied

Validation of Cryogenic Propellant Tank Filling using Computational Fluid Dynamics Simulation

Validation of cryogenic propellant tank filling was performed using the Computational Fluid Dynamics (CFD) solver Loci/STREAM-VoF. The validation effort helped identify modeling methodologies that enable NASA to best support its partners in both launch pad and on-orbit filling operations. Data from liquid hydrogen ground tests filled via jet injection were used for the validation effort which include a sensitivity to initial tank wall temperature. Initially hot walls are expected to yield rapid evaporation and possibly boiling. A Volume of Fluid (VoF) methodology was used to capture the gas-liquid interface. Rapid breakup of the liquid jet was observed in simulation results as liquid evaporated and expanded. Inflowing liquid transitioned to a contiguous jet as tank temperatures decreased, gas pressure increased, and saturation conditions at the incoming liquid temperature were approached. The final phase of filling was distinguished by rapid gas pressure rise due to a higher rate of gas volume compression than condensation at the liquid surface. Key physics of propellant tank filling were captured in the computational predictions, and opportunities for added simulation robustness and efficiency in future modeling efforts were identified.

CFM↗
DOE OSTI2022

Beam Dynamics Simulations of Transient Beam Loading Effects in the 10GeV EIC Electron Storage Ring

We report on beam dynamics studies of transient beam loading effects in the 10 GeV EIC electron storage ring [1]. The studies are carried out with time-dependent Vlasov-Fokker-Planck simulations performed with the parallel, particle tracking code SPACE [2], which allows to follow self-consistently the dynamics of h bunches, where h in the number of RF buckets, in arbitrary multi-bunch configurations. The specific goal of the numerical simulations is to determine stable RF cavity settings under heavy beam loading. We also study the option to operate with a passive, third-harmonic cavity (3HC) system for bunch lengthening, addressing both stability and the performance limitation due to a gap in the uniform filling pattern for ion clearing.

43 PARTICLE ACCELERATORS↗
NASA NTRS1997

Molecular Dynamics Simulation of a Multi-Walled Carbon Nanotube Based Gear

We used molecular dynamics to investigate the properties of a multi-walled carbon nanotube based gear. Previous work computationally suggested that molecular gears fashioned from (14,0) single-walled carbon nanotubes operate well at 50-100 gigahertz. The gears were formed from nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. The gear in this study was based on the smallest multi-walled nanotube supported by some experimental evidence. Each gear was a (52,0) nanotube surrounding a (37,10) nanotube with approximate 20.4 and 16,8 A radii respectively. These sizes were chosen to be consistent with inter-tube spacing observed by and were slightly larger than graphite inter-layer spacings. The benzyne teeth were attached via 2+4 cycloaddition to exterior of the (52,0) tube. 2+4 bonds were used rather than the 2+2 bonds observed by Hoke since 2+4 bonds are preferred by naphthalene and quantum calculations by Jaffe suggest that 2+4 bonds are preferred on carbon nanotubes of sufficient diameter. One gear was 'powered' by forcing the atoms near the end of the outside buckytube to rotate to simulate a motor. A second gear was allowed to rotate by keeping the atoms near the end of its outside buckytube on a cylinder. The ends of both gears were constrained to stay in an approximately constant position relative to each other, simulating a casing, to insure that the gear teeth meshed. The stiff meshing aromatic gear teeth transferred angular momentum from the powered gear to the driven gear. The simulation was performed in a vacuum and with a software thermostat. Preliminary results suggest that the powered gear had trouble turning the driven gear without slip. The larger radius and greater mass of these gears relative to the (14,0) gears previously studied requires a smaller rotation rate and multiple rows of teeth to avoid excessive force on the gear teeth resulting, in slip and failure of the driven gear to turn. We hope that studies such as these will eventually lead to synthesis of components that can be assembled into atomically precise fullerene machines. These machines, in turn, may someday be used in machine-phase fullerene materials with remarkable properties.

Han, Jie↗
NASA NTRS1990

Staggered solution procedures for multibody dynamics simulation

A computational procedure is presented for the direct integration of the multibody dynamical (MBD) equations with constraints. Numerical experiments conducted with this MBD procedure indicate that it yields robust solutions in cases where the step size furnishes more than 20 samples for the period of highest apparent response frequency in a given multibody system. The MBD solution procedure is implemented in two separate modules: the generalized coordinated solver, 'CINT', and the constraint Lagrange multiplier solver, 'LINT'.

Park, K. C.↗
DOE OSTI2025

First Principles Molecular Dynamics Simulations of Ammonia Adsorption onto MFI Zeolite Nanosheets

MFI zeolite nanosheets membranes are promising candidates for ammonia separation from nitrogen and hydrogen, yet questions remain on the origin of their high selectivity. Silanols, Si-OH, are present in high concentration at the surface of zeolite nanosheets, and force-field-based simulations indicate that surface adsorption at the silanols contributes to selectivity. Silanols can chemically react with ammonia, which may further contribute to the ability of zeolite nanosheet membranes to separate it from other gases. In this work, we use first-principles molecular dynamics techniques to simulate ammonia’s behavior within stacked MFI zeolite nanosheets. We find that at 523 K and a loading corresponding to 35 bar, conditions desired for the ammonia separation, about 30% of ammonia reacts with surface silanols. Our work explores H-bonding and proton transfer within this system.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Orientational Pair Correlations and Local Structure of Benzonitrile from Molecular Dynamics Simulations with Comparisons to Experiments [plus supplemental information]

Here, we present an experimentally parametrized molecular dynamics study of single-molecule and collective orientational relaxation in neat benzonitrile through the analysis of the reorientational anisotropy and polarizability anisotropy time-correlation function (PA-TCF). The simulations show that the PA-TCF is dominated by collective reorientation after 20 ps. Collective reorientation is found to be slower than single-molecule reorientation by a factor of 1.67, consistent with recent experiments. The simulations provide direct evidence of local antiparallel benzonitrile configurations. These structures, which have been the center of some debate, are responsible for the slower rate of collective vs. single-molecule reorientation in the liquid. Further structural analysis indicates that significant Coulombic interactions between the nitrile group and hydrogen atoms on adjacent molecules play a role in the formation of the antiparallel structures. The single-molecule dynamics reflected in the anisotropy are complex and consist of a ballistic regime, restricted angular diffusion, and spatially anisotropic free diffusion. The principal components of the rotational diffusion tensor are independently obtained and shown to reproduce the free diffusion regime of the anisotropy for each principal axis according to the predictions of a previous theory.

14 SOLAR ENERGY↗
DOE OSTI2018

The Kinetics of the Liquid Phase Nucleation in a Stretched FCC Crystal: A Molecular Dynamics Simulation

The kinetics of a spontaneous formation of liquid phase in a stretched (superheated) Lennard-Jones fcc crystal is studied. Molecular dynamics experiments are used to determine the main parameters of the nucleation process: nucleation frequency J, diffusion coefficient of nuclei D{sub *}, nonequilibrium Zel’dovich factor Z, and critical nucleus size R{sub *}. The calculations are performed at negative pressures from the endpoint of the melting line and at positive pressures that are higher by a factor of eight than the critical pressure. The simulation results are compared to the classical homogeneous nucleation theory. It is found that the theory qualitatively correctly reproduces the dynamics of developing the process. The theory and the simulation demonstrate good quantitative agreement for the transition rate of the liquid phase nucleus through the critical size, but there is large difference in the numbers of critical nuclei in the unit volume of the metastable phase. In the case of significant superheatings and negative pressures, the contribution of the energy of elastic stresses to the moving force of the phase transformation is small and it can be neglected in a first approximation. The mismatch between the theory and the simulation results can be eliminated taking that the surface free energy of a curved “crystal–liquid droplet” interface is smaller than that of a plane interface by 30–35%.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗