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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Ab initio implementation of quantum trajectory mean-field approach and dynamical simulation of the N{sub 2}CO photodissociation

In this work, the recently introduced quantum trajectory mean-field (QTMF) approach is implemented and employed to explore photodissociation dynamics of diazirinone (N{sub 2}CO), which are based on the high-level ab initio calculation. For comparison, the photodissociation process has been simulated as well with the fewest-switches surface hopping (FSSH) and the ab initio multiple spawning (AIMS) methods. Overall, the dynamical behavior predicted by the three methods is consistent. The N{sub 2}CO photodissociation at λ > 335 nm is an ultrafast process and the two C—N bonds are broken in a stepwise way, giving birth to CO and N{sub 2} as the final products in the ground state. Meanwhile, some noticeable differences were found in the QTMF, FSSH, and AIMS simulated time constants for fission of the C—N bonds, excited-state lifetime, and nonadiabatic transition ratios in different intersection regions. These have been discussed in detail. The present study provides a clear evidence that direct ab initio QTMF approach is one of the reliable tools for simulating nonadiabatic dynamics processes.

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

System dynamic simulation: A new method in social impact assessment (SIA)

Many complex social questions are difficult to address adequately with conventional methods and techniques, due to the complicated dynamics, and hard to quantify social processes. Despite these difficulties researchers and practitioners have attempted to use conventional methods not only in evaluative modes but also in predictive modes to inform decision making. The effectiveness of SIAs would be increased if they were used to support the project design processes. This requires deliberate use of lessons from retrospective assessments to inform predictive assessments. Social simulations may be a useful tool for developing a predictive SIA method. There have been limited attempts to develop computer simulations that allow social impacts to be explored and understood before implementing development projects. In light of this argument, this paper aims to introduce system dynamic (SD) simulation as a new predictive SIA method in large development projects. We propose the potential value of the SD approach to simulate social impacts of development projects. We use data from the SIA of Gareh-Bygone floodwater spreading project to illustrate the potential of SD simulation in SIA. It was concluded that in comparison to traditional SIA methods SD simulation can integrate quantitative and qualitative inputs from different sources and methods and provides a more effective and dynamic assessment of social impacts for development projects. We recommend future research to investigate the full potential of SD in SIA in comparing different situations and scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Quantifying Dimer and Trimer Formation by Tri- n -butyl Phosphates in n -Dodecane: Molecular Dynamics Simulations

Tri-n-butyl phosphate (TBP), a representative of neutral organophosphorous ligands, is an important extractant used in solvent extraction process for the recovery of uranium and plutonium from spent nuclear fuel. Microscopic pictures of TBP isomerism and its behavior in n-dodecane diluent were investigated utilizing MD simulations with previously optimized force field parameters for TBP and n-dodecane. Potential Mean Force (PMF) calculations on a single TBP molecule show seven probable TBP isomers. Radial Distribution Functions (RDF) of TBP suggests the existence of TBP trimers at high TBP concentrations in addition to dimers. 2D PMF calculations were performed to determine the angle and distance criteria for TBP trimers. The dimerization and trimerization constants of TBP in n-dodecane were obtained and match our own experimental values using FTIR technique. The new insights into the conformational behaviors of TBP molecule as a monomer and as part of an aggregate could greatly aid the understanding of the complexation between TBP and metal ions in solvent extraction system. The U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences and Biosciences funded the work performed by LXD.

equilibrium constants↗
DOE OSTI2021

Local structure analysis of low-temperature neutron pair distribution function coupled with molecular dynamics simulations of CH 4 and CO 2 hydrates from 2 to 210 K

CH 4 hydrates occur naturally and are an abundant potential fuel source with a corresponding risk of potent greenhouse gas release, due to their stability conditions at low temperature and high pressure. Byproduct CO 2 can be exchanged with CH 4 in these natural deposits and can potentially stabilize the hydrates at higher temperatures. CH 4 , CO 2 , and mixed CH 4 -CO 2 hydrates are studied with in situ neutron pair distribution function experiments from 2 to 210 K to investigate the impact of varying the CH 4 -CO 2 guest composition in the gas hydrate structure. These experiments combined with Reverse Monte Carlo analysis allow for the characterization of intermolecular CH 4 and CO 2 interactions with the water molecules which form the hydrate lattice and how they impact the local structure of the lattice itself. Results indicate that when CH 4 and CO 2 co-occupy the hydrate, the host is more strongly distorted than in the pure CH 4 and CO 2 hydrates, but this becomes less defined with increasing temperature. The presence of CO 2 in mixed hydrate increases the stability range and creates a barrier for CH 4 to completely leave the structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗