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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Investigating the time dependence of neutron-proton equilibration using molecular dynamics simulations

Neutron-proton equilibration has previously been studied experimentally in dynamically deformed nuclei in heavy ion collisions [Phys. Rev. Lett. 118, 062501 (2017); Phys. Rev. C 95, 044604 (2017)] using the NIMROD detector array at the Cyclotron Institute at Texas A&M University. Results indicated the composition of the two heaviest fragments originating from the excited projectile-like fragment evolved exponential with respect to their orientation angle to be more similar. Constrained Molecular Dynamics and Antisymmetrized Molecular Dynamics simulations were performed for different formulations of the density dependence of the asymmetry energy term of the nuclear equation of state. The simulations are compared to experimental results, which indicate a better agreement with a softer interaction in the nuclear equation of state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2011

Molecular Dynamics Simulation of the Antiamoebin Ion Channel: Linking Structure and Conductance

Molecular dynamics simulations were carried out in order to ascertain which of the potential multimeric forms of the transmembrane peptaibol channel, antiamoebin, is consistant with its measured conductance. Estimates of the conductance obtained through counting ions that cross the channel and by solving the Nernst-Planck equation yield consistent results, indicating that the motion of ions inside the channel can be satisfactorily described as diffusive.The calculated conductance of octameric channels is markedly higher than the conductance measured in single channel recordings, whereas the tetramer appears to be non-conducting. The conductance of the hexamer was estimated to be 115+/-34 pS and 74+/-20 pS, at 150 mV and 75 mV, respectively, in satisfactory agreement with the value of 90 pS measured at 75 mV. On this basis we propose that the antiamoebin channel consists of six monomers. Its pore is large enough to accommodate K(+) and Cl(-) with their first solvation shells intact. The free energy barrier encountered by K(+) is only 2.2 kcal/mol whereas Cl(-) encounters a substantially higher barrier of nearly 5 kcal/mol. This difference makes the channel selective for cations. Ion crossing events are shown to be uncorrelated and follow Poisson statistics. keywords: ion channels, peptaibols, channel conductance, molecular dynamics

Wilson, Michael A.↗
DOE OSTI2015

Photoinduced reactions of both 2-formyl-2H-azirine and isoxazole: A theoretical study based on electronic structure calculations and nonadiabatic dynamics simulations

In the present work, the combined electronic structure calculations and dynamics simulations have been performed to explore photocleavages of 2-formyl-2H-azirine and isoxazole in the gas phase and the subsequent rearrangement reactions. The carbonyl n → π{sup *} transition induces a cleavage of the C—N single bond of 2-formyl-2H-azirine to yield β-formylvinylnitrene in open-shell singlet state. However, the n → π{sup *} excitation of the imine chromophore results in a cleavage of the C—C single bond, producing a nitrile ylide intermediate through an internal conversion to the ground state. β-formylvinylnitrene and nitrile ylide with the carbonyl group are easily transformed into 2-formyl-2H-azirine and oxazole, respectively. The N—O bond cleavages on both S{sub 1}({sup 1}ππ{sup *}) and S{sub 2}({sup 1}n{sub N}π{sup *}) of isoxazole are ultrafast processes, and they give products of 2-formyl-2H-azirine, 3-formylketenimine, HCN + CHCHO, and HCO + CHCHN. Both 2H-azirines and ketenimines were suggested to be formed from the triplet vinylnitrenes by intersystem crossing in the previous studies. However, our calculations show that the singlet β-formylvinylnitrene is responsible for the formation of 2-formyl-2H-azirine and 3-formylketenimine, and the singlet vinylnitrenes can play a key role in the photoinduced reactions of both 2H-azirines and isoxazoles.

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

Sensitivity Analysis and Uncertainty Quantification for the LAMMPS Molecular Dynamics Simulation Code

We examine sensitivity analysis and uncertainty quantification for molecular dynamics simulation. Extreme (large or small) output values for the LAMMPS code often occur at the boundaries of input regions, and uncertainties in those boundary values are overlooked by common SA methods. Similarly, input values for which code outputs are consistent with calibration data can also occur near boundaries. Upon applying approaches in the literature for imprecise probabilities (IPs), much more realistic results are obtained than for the complacent application of standard SA and code calibration.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1989

A model for including thermal conduction in molecular dynamics simulations

A technique is introduced for including thermal conduction in molecular dynamics simulations for solids. A model is developed to allow energy flow between the computational cell and the bulk of the solid when periodic boundary conditions cannot be used. Thermal conduction is achieved by scaling the velocities of atoms in a transitional boundary layer. The scaling factor is obtained from the thermal diffusivity, and the results show good agreement with the solution for a continuous medium at long times. The effects of different temperature and size of the system, and of variations in strength parameter, atomic mass, and thermal diffusivity were investigated. In all cases, no significant change in simulation results has been found.

Wu, Yue↗
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 OSTI2015

A model of lipid-free Apolipoprotein A-I revealed by iterative molecular dynamics simulation

Apolipoprotein A-I (apo A-I), the major protein component of high-density lipoprotein, has been proven inversely correlated to cardiovascular risk in past decades. The lipid-free state of apo A-I is the initial stage which binds to lipids forming high-density lipoprotein. Molecular models of lipid-free apo A-I have been reported by methods like X-ray crystallography and chemical cross-linking/mass spectrometry (CCL/MS). Through structural analysis we found that those current models had limited consistency with other experimental results, such as those from hydrogen exchange with mass spectrometry. Through molecular dynamics simulations, we also found those models could not reach a stable equilibrium state. Therefore, by integrating various experimental results, we proposed a new structural model for lipidfree apo A-I, which contains a bundled four-helix N-terminal domain (1–192) that forms a variable hydrophobic groove and a mobile short hairpin C-terminal domain (193–243). This model exhibits an equilibrium state through molecular dynamics simulation and is consistent with most of the experimental results known from CCL/MS on lysine pairs, fluorescence resonance energy transfer and hydrogen exchange. This solution-state lipid-free apo A-I model may elucidate the possible conformational transitions of apo A-I binding with lipids in high-density lipoprotein formation.

59 BASIC BIOLOGICAL SCIENCES↗
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↗
DOE OSTI2014

Molecular dynamics simulations of intergranular fracture in UO2 with nine empirical interatomic potentials

The intergranular fracture behavior of UO2 was studied using molecular dynamics simulations with a bicrystal model. The anisotropic fracture behavior due to the different grain boundary characters was investigated with the View the MathML source symmetrical tilt S5 and the View the MathML source symmetrical tilt S3 ({1 1 1} twin) grain boundaries. Nine interatomic potentials, seven rigid-ion plus two core–shell ones, were utilized to elucidate possible potential dependence. Initiating from a notch, crack propagation along grain boundaries was observed for most potentials. The S3 boundary was found to be more prone to fracture than the S5 one, indicated by a lower energy release rate associated with the former. However, some potential dependence was identified on the existence of transient plastic deformation at crack tips, and the results were discussed regarding the relevant material properties including the excess energies of metastable phases and the critical energy release rate for intergranular fracture. In general, local plasticity at crack tips was observed in fracture simulations with potentials that predict low excess energies for metastable phases and high critical energy release rates for intergranular fracture.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

Peptide dynamics by molecular dynamics simulation and diffusion theory method with improved basis sets

Improved basis sets for the study of polymer dynamics by means of the diffusion theory, and tests on a melt of cis-1,4-polyisoprene decamers, and a toluene solution of a 71-mer syndiotactic trans-1,2-polypentadiene were presented recently [R. Gaspari and A. Rapallo, J. Chem. Phys. 128, 244109 (2008)]. The proposed hybrid basis approach (HBA) combined two techniques, the long time sorting procedure and the maximum correlation approximation. The HBA takes advantage of the strength of these two techniques, and its basis sets proved to be very effective and computationally convenient in describing both local and global dynamics in cases of flexible synthetic polymers where the repeating unit is a unique type of monomer. The question then arises if the same efficacy continues when the HBA is applied to polymers of different monomers, variable local stiffness along the chain and with longer persistence length, which have different local and global dynamical properties against the above-mentioned systems. Important examples of this kind of molecular chains are the proteins, so that a fragment of the protein transthyretin is chosen as the system of the present study. This peptide corresponds to a sequence that is structured in β-sheets of the protein and is located on the surface of the channel with thyroxin. The protein transthyretin forms amyloid fibrils in vivo, whereas the peptide fragment has been shown [C. P. Jaroniec, C. E. MacPhee, N. S. Astrof, C. M. Dobson, and R. G. Griffin, Proc. Natl. Acad. Sci. U.S.A. 99, 16748 (2002)] to form amyloid fibrils in vitro in extended β-sheet conformations. For these reasons the latter is given considerable attention in the literature and studied also as an isolated fragment in water solution where both experimental and theoretical efforts have indicated the propensity of the system to form β turns or α helices, but is otherwise predominantly unstructured. Differing from previous computational studies that employed implicit solvent, we performed in this work the classical molecular dynamics simulation on a realistic model solution with the peptide embedded in an explicit water environment, and calculated its dynamic properties both as an outcome of the simulations, and by the diffusion theory in reduced statistical-mechanical approach within HBA on the premise that the mode-coupling approach to the diffusion theory can give both the long-range and local dynamics starting from equilibrium averages which were obtained from detailed atomistic simulations.

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

Direct Dynamics Simulation of Dissociation of the [CH3--I--OH]- Ion-Molecule Complex

Direct dynamics simulations were used to study dissociation of the [CH3--I--OH]- complex ion, which was observed in a previous study of the OH- + CH3I gas phase reaction (J. Phys. Chem. A 2013, 117, 7162). Restricted B97-1 simulations were performed to study dissociation at 65, 75 and 100 kcal/mol and the [CH3--I--OH]- ion dissociated exponentially, in accord with RRKM theory. For these energies the major dissociation products are CH3I + OH-, CH2I- + H2O, and CH3OH + I-. Unrestricted B97-1 and restricted and unrestricted CAM-B3LYP simulations were also performed at 100 kcal/mol to compare with the restricted B97-1 results. The {CH3I + OH-}:{CH2I- + H2O}:{CH3OH + I-} product ratio is 0.72 : 0.15 : 0.13, 0.81 : 0.05 : 0.14, 0.71 : 0.19 : 0.10 , and 0.83 : 0.13 : 0.04 for the restricted B97-1, unrestricted B97-1, restricted CAM-B3LYP, and unrestricted CAM-B3LYP simulations, respectively. Other product channels found are CH2 + I- + H2O, CH2 + I-(H2O), CH4 + IO-, CH3 - + IOH, and CH3 + IOH-. The CH3 - + IOH singlet products are only given by the restricted B97-1 simulation and the lower energy CH3 + IOH- doublet products are only formed by the unrestricted B97-1 simulation. Also studied were the direct and indirect atomic-level mechanisms for forming CH3I + OH-, CH2I- + H2O, and CH3OH + I-. The majority of CH3I + OH- were formed through a direct mechanism. For both CH2I- + H2O and CH3OH + I-, the direct mechanism is overall more important than the indirect mechanisms, with the round-about like mechanism the most important indirect mechanism at high excitation energies. Mechanism comparisons between the B97-1 and CAM-B3LYP simulations showed that formation of the CH3OH---I- complex is favored for the B97-1 simulations, while formation of the HO----HCH2I complex is favored for the CAM-B3LYP simulations. The unrestricted simulations give a higher percentage of indirect mechanisms than the restricted simulations. The possible role of the self-interaction error in the simulations is also discussed. The work presented here gives a detailed picture of the [CH3--I--OH]- dissociation dynamics, and is very important for unraveling the role of [CH3--I--OH]- in the dynamics of the OH-(H2O)n=1,2 + CH3I reactions.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2016

Molecular dynamics simulation of diffusion and electrical conductivity in montmorillonite interlayers

In this study, the diffusion of water and ions in the interlayer region of smectite clay minerals represents a direct probe of the type and strength of clay–fluid interactions. Interlayer diffusion also represents an important link between molecular simulation and macroscopic experiments. Here we use molecular dynamics simulation to investigate trends in cation and water diffusion in montmorillonite interlayers, looking specifically at the effects of layer charge, interlayer cation and cation charge (sodium or calcium), water content, and temperature. For Na-montmorillonite, the largest increase in ion and water diffusion coefficients occurs between the one-layer and two-layer hydrates, corresponding to the transition from inner-sphere to outer-sphere surface complexes. Calculated activation energies for ion and water diffusion in Na-montmorillonite are similar to each other and to the water hydrogen bond energy, suggesting the breaking of water–water and water–clay hydrogen bonds as a likely mechanism for interlayer diffusion. A comparison of interlayer diffusion with that of bulk electrolyte solutions reveals a clear trend of decreasing diffusion coefficient with increasing electrolyte concentration, and in most cases the interlayer diffusion results are nearly coincident with the corresponding bulk solutions. Trends in electrical conductivities computed from the ion diffusion coefficients are also compared.

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

Ab initio molecular dynamics simulation study of successive hydrogenation reactions of carbon monoxide producing methanol

Doing ab initio molecular dynamics simulations, we demonstrate a possibility of hydrogenation of carbon monoxide producing methanol step by step. At first, the hydrogen atom reacts with the carbon monoxide molecule at the excited state forming the formyl radical. Formaldehyde was formed after adding one more hydrogen atom to the system. Finally, absorption of two hydrogen atoms to formaldehyde produces methanol molecule. This study is performed by using the all-electron mixed basis approach based on the time dependent density functional theory within the adiabatic local density approximation for an electronic ground-state configuration and the one-shot GW approximation for an electronic excited state configuration.

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