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

Reactive Molecular Dynamics Simulation of Kerogen Thermal Maturation and Cross-Linking Pathways

For this study, molecular dynamics simulations were performed with a ReaxFF reactive force field to investigate bond breaking and bond formation mechanisms during the early thermal maturation of kerogen and potential crosslinking pathways towards a three-dimensional (3D) quasi-infinite molecular network (crosslinked kerogen). Starting with small ensembles of high molecular mass models for immature type I Green River Shale kerogen (Kerogen 1-I), top of the oil window type II kerogen (Kerogen 2-L), and low maturity type III kerogen (Kerogen 3-L), low molecular mass species including H 2 O, H 2 S, C 2 H 4 , C 3 H 6 were produced as the maturities of the remaining kerogens increased. Highly reactive fragments, which are not detected in pyrolysis experiments, were also produced. Further, the crosslinking mechanism in the newly develop polymeric kerogen networks appears to be highly complex and it is influenced by the nature of the kerogen. Covalent —C—S—, —C—O— and —C—C— bonds are the primary crosslinks that structurally binds kerogen monomers together. The trends observed in the thermochemical transformation of kerogen and the kerogen crosslinking pathways are consistent with the theoretical and experimental studies reported in the scientific literature. The reactive force field molecular dynamics simulation provides a potentially valuable alternative approach to the development of a realistic 3D molecular network of kerogen. However, the conversion of kerogen molecules into a 3D molecular network was low.

02 PETROLEUM↗
DOE OSTI2023

The selective heating effect of microwave irradiation on a binary mixture of water and polyethylene oxide: a molecular dynamics simulation approach

In this study, we investigate the molecular mechanisms of a microwave-driven selective heating process by performing molecular dynamics simulations for three different systems including pure water, pure polyethylene oxide (PEO), and water-PEO mixed systems in the presence of a microwave with two different intensities of electric field such as 0.001 V/Å -1 and 0.01 V/Å -1 at a frequency of 100 GHz. First, from performing molecular dynamics simulations of CO and CO 2 in the presence of the microwave, it is confirmed that the molecular dipole moment is responsible for the rotational motion induced by the oscillating electric field. Second, by analyzing the MD simulations of the pure water system, we discover that the dipole moment of water exhibits a time lag with respect to the microwave. During the heating process, however, the temperature, kinetic, and potential energies increase synchronously with the oscillating electric field of the microwave, showing that the heating of the water system is caused by the molecular reaction of water to the microwave. Comparing the water-PEO mixed system to the pure water and pure PEO systems, the water-PEO mixed system has a higher heating rate than the pure PEO system but a lower heating rate than the pure water system. Therefore, we conclude that heating the water-PEO mixed system is driven by water molecules selectively activated by microwave irradiation. Furthermore, we also calculate the diffusion coefficients of water molecules and PEO chains by describing their mean square displacements, demonstrating that the diffusion coefficients are increased in the presence of microwaves for both water and PEO in pure and mixed systems. Lastly, during the microwave heating process, the structures of the water-PEO mixed system are altered as a function of the intensity of electric field, which is mainly driven by the response of water molecules.

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

Ensemble Docking in Drug Discovery: How Many Protein Configurations from Molecular Dynamics Simulations are Needed To Reproduce Known Ligand Binding?

Ensemble docking in drug discovery or chemical biology uses dynamical simulations of target proteins to generate binding site conformations for docking campaigns. We demonstrate that 600 ns molecular dynamics simulations of four G-protein-coupled receptors in their membrane environments generate ensembles of protein configurations that, collectively, are selected by 70–99% of the known ligands of these proteins. Thus, the process of ligand recognition by conformational selection can be reproduced by combining molecular dynamics and docking calculations. Clustering of the molecular dynamics trajectories, however, does not necessarily identify the protein conformations that are most often selected by the ligands.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Recovering position-dependent diffusion from biased molecular dynamics simulations

All atom molecular dynamics (MD) models provide valuable insight into the dynamics of biophysical systems, but are limited in size or length by the high computational demands. The latter can be reduced by simulating long term diffusive dynamics (also known as Langevin dynamics or Brownian motion) of the most interesting and important user-defined parts of the studied system, termed collective variables (colvars). A few hundred nanosecond-long biased MD trajectory can therefore be extended to millisecond lengths in the colvars subspace at a very small additional computational cost. In this work, we develop a method for determining multidimensional anisotropic position- and timescale-dependent diffusion coefficients (D) by analysing the changes of colvars in an existing MD trajectory. As a test case, we obtained D for dihedral angles of the alanine dipeptide. An open source Mathematica{sup ®} package, capable of determining and visualizing D in one or two dimensions, is available at https://github.com/lbf-ijs/DiffusiveDynamics . Given known free energy and D, the package can also generate diffusive trajectories.

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

Synthesis of formamide and related organic species in the interstellar medium via chemical dynamics simulations

We show, by means of direct dynamics simulations, how it is possible to define possible reactants and mechanisms leading to the formation of formamide in the interstellar medium. In particular, different ion–molecule reactions in the gas phase were considered: NH{sub 3}OH{sup +}, NH{sub 2}OH{sub 2}{sup +}, H{sub 2}COH{sup +}, and NH{sub 4} {sup +} for the ions and NH{sub 2}OH, H{sub 2}CO, and NH{sub 3} for the partner neutrals. These calculations were combined with high level ab initio calculations to investigate possible further evolution of the products observed. In particular, for formamide, we propose that the NH{sub 2}OH{sub 2}{sup +} + H{sub 2}CO reaction can produce an isomer, NH{sub 2}OCH{sub 2}{sup +}, that, after dissociative recombination, can produce neutral formamide, which was observed in space. The direct dynamics do not pre-impose any reaction pathways and in other reactions, we did not observe the formation of formamide or any possible precursor. On the other hand, we obtained other interesting reactions, like the formation of NH{sub 2}CH{sub 2}{sup +}. Finally, some radiative association processes are proposed. All of the results obtained are discussed in light of the species observed in radioastronomy.

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

Traction–separation relationships for hydrogen induced grain boundary embrittlement in nickel via molecular dynamics simulations

A statistical approach combined with molecular dynamics simulations is used to study the influence of hydrogen on intergranular decohesion. This methodology is applied to a Ni Σ3(112)[11¯0] symmetric tilt grain boundary. Hydrogenated grain boundaries with different H concentrations are constructed using an energy minimization technique with initial H atom positions guided by Monte Carlo simulation results. Decohesion behavior is assessed through extraction of a traction–separation relationship during steady-state crack propagation in a statistically meaningful approach, building upon prior work employing atomistic cohesive zone volume elements (CZVEs). A sensitivity analysis is performed on the numerical approach used to extract the traction–separation relationships, clarifying the role of CZVE size, threshold parameters necessary to differentiate elastic and decohesion responses, and the numerical averaging technique. Results show that increasing H coverage at the Ni Σ3(112)[11¯0] grain boundary asymmetrically influences the crack tip velocity during propagation, leads to a general decrease in the work of separation required for crack propagation, and provides a reduction in the peak stress in the extracted traction–separation relationship. Furthermore the present framework offers a meaningful vehicle to pass atomistically derived interfacial behavior to higher length scale formulations for intergranular fracture.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Staggered solution procedures for multibody dynamics simulation

The numerical solution procedure for multibody dynamics (MBD) systems is termed a staggered MBD solution procedure that solves the generalized coordinates in a separate module from that for the constraint force. This requires a reformulation of the constraint conditions so that the constraint forces can also be integrated in time. A major advantage of such a partitioned solution procedure is that additional analysis capabilities such as active controller and design optimization modules can be easily interfaced without embedding them into a monolithic program. After introducing the basic equations of motion for MBD system in the second section, Section 3 briefly reviews some constraint handling techniques and introduces the staggered stabilized technique for the solution of the constraint forces as independent variables. The numerical direct time integration of the equations of motion is described in Section 4. As accurate damping treatment is important for the dynamics of space structures, we have employed the central difference method and the mid-point form of the trapezoidal rule since they engender no numerical damping. This is in contrast to the current practice in dynamic simulations of ground vehicles by employing a set of backward difference formulas. First, the equations of motion are partitioned according to the translational and the rotational coordinates. This sets the stage for an efficient treatment of the rotational motions via the singularity-free Euler parameters. The resulting partitioned equations of motion are then integrated via a two-stage explicit stabilized algorithm for updating both the translational coordinates and angular velocities. Once the angular velocities are obtained, the angular orientations are updated via the mid-point implicit formula employing the Euler parameters. When the two algorithms, namely, the two-stage explicit algorithm for the generalized coordinates and the implicit staggered procedure for the constraint Lagrange multipliers, are brought together in a staggered manner, they constitute a staggered explicit-implicit procedure which is summarized in Section 5. Section 6 presents some example problems and discussions concerning several salient features of the staggered MBD solution procedure are offered in Section 7.

Park, K. C.↗
NASA NTRS1997

Molecular Dynamics Simulations of Laser Powered Carbon Nanotube Gears

Dynamics of laser powered carbon nanotube gears is investigated by molecular dynamics simulations with Brenner's hydrocarbon potential. We find that when the frequency of the laser electric field is much less than the intrinsic frequency of the carbon nanotube, the tube exhibits an oscillatory pendulam behavior. However, a unidirectional rotation of the gear with oscillating frequency is observed under conditions of resonance between the laser field and intrinsic gear frequencies. The operating conditions for stable rotations of the nanotube gears, powered by laser electric fields are explored, in these simulations.

Srivastava, Deepak↗
DOE OSTI2015

Dynamic Simulation over Long Time Periods with 100% Solar Generation.

This project aimed to identify the path forward for dynamic simulation tools to accommodate these needs by characterizing the properties of power systems (with high PV penetration), analyzing how these properties affect dynamic simulation software, and offering solutions for potential problems.

14 SOLAR ENERGY↗
DOE OSTI2018

ReaxFF molecular dynamics simulation of intermolecular structure formation in acetic acid-water mixtures at elevated temperatures and pressures

The intermolecular structure formation in liquid and supercritical acetic acid-water mixtures was explored using ReaxFF-based molecular dynamics simulations. The microscopic structures of acetic acid-water mixtures with different acetic acid mole fractions (1.0 ≥ x HAc ≥ 0.2) at ambient and critical conditions were examined. The potential energy surface associated with the dissociation of acetic acid molecules was calculated using a metadynamics procedure to optimize the dissociation energy of ReaxFF potential. At ambient conditions, depending on the acetic acid concentration, either acetic acid clusters or water clusters are dominant in the liquid mixture. When acetic acid is dominant (0.4 ≤ x HAc ), cyclic dimers and chain structures between acetic acid molecules are present in the mixture. Both structures disappear at increased water content of the mixture. It was found by simulations that the acetic acid molecules released from these dimer and chain structures tend to stay in a dipole-dipole interaction. These structural changes are in agreement with the experimental findings. When switched to critical conditions, the long-range interactions (e.g., second or fourth neighbor) disappear and the water-water and acetic acid-acetic acid structural formations become disordered. The simulated radial distribution function for water-water interactions is in agreement with experimental and computational studies. The first neighbor interactions between acetic acid and water molecules are preserved at relatively lower temperatures of the critical region. As higher temperatures are reached in the critical region, these interactions were observed to weaken. These simulations suggest that ReaxFF molecular dynamics simulations are an appropriate tool for studying supercritical water/organic acid mixtures.

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

Characterization of three phases of liquid carbon by tight-binding molecular dynamics simulations

In this work, we have performed systematic molecular dynamics simulations to study the structures of liquid carbon at 5000 K with the weight density ranging from 1.4 to 3.5 g cm –3 , using a three-center tight-binding potential of carbon. The simulation results show that the bonding characteristics of the liquid changes predominately from twofold to threefold, and then to fourfold coordination as the density increases. Signals of two structural changes at the densities of about 1.9 and 3.0 g cm –3 respectively are revealed by the slope changes in the density dependence of structural, electronic and dynamical properties. Our simulation results suggest that there are three distinct liquid carbon phases in this density range. However, further thermodynamics calculations, e.g. , free energy calculations, would be required to clarify the possible liquid–liquid transitions.

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

Equilibration and analysis of first-principles molecular dynamics simulations of water

First-principles molecular dynamics (FPMD) simulations based on density functional theory are becoming increasingly popular for the description of liquids. In view of the high computational cost of these simulations, the choice of an appropriate equilibration protocol is critical. In this study, we assess two methods of estimation of equilibration times using a large dataset of first-principles molecular dynamics simulations of water. The Gelman-Rubin potential scale reduction factor [A. Gelman and D. B. Rubin, Stat. Sci. 7, 457 (1992)] and the marginal standard error rule heuristic proposed by White [Simulation 69, 323 (1997)] are evaluated on a set of 32 independent 64-molecule simulations of 58 ps each, amounting to a combined cumulative time of 1.85 ns. The availability of multiple independent simulations also allows for an estimation of the variance of averaged quantities, both within MD runs and between runs. We analyze atomic trajectories, focusing on correlations of the Kohn-Sham energy, pair correlation functions, number of hydrogen bonds, and diffusion coefficient. The observed variability across samples provides a measure of the uncertainty associated with these quantities, thus facilitating meaningful comparisons of different approximations used in the simulations. We find that the computed diffusion coefficient and average number of hydrogen bonds are affected by a significant uncertainty in spite of the large size of the dataset used. Lastly, a comparison with classical simulations using the TIP4P/2005 model confirms that the variability of the diffusivity is also observed after long equilibration times. Complete atomic trajectories and simulation output files are available online for further analysis.

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

Hypercrosslinked polystyrene networks: An atomistic molecular dynamics simulation combined with a mapping/reverse mapping procedure

An algorithm involving classical molecular dynamics simulations with mapping and reverse mapping procedure is here suggested to simulate the crosslinking of the polystyrene dissolved in dichloroethane by monochlorodimethyl ether. The algorithm comprises consecutive stages: molecular dynamics atomistic simulation of a polystyrene solution, the mapping of atomistic structure onto coarse-grained model, the crosslink formation, the reverse mapping, and finally relaxation of the structure dissolved in dichloroethane and in dry state. The calculated values of the specific volume and the elastic modulus are in reasonable quantitative correspondence with experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1995

Determination of Quantum Chemistry Based Force Fields for Molecular Dynamics Simulations of Aromatic Polymers

Ab initio quantum chemistry calculations for model molecules can be used to parameterize force fields for molecular dynamics simulations of polymers. Emphasis in our research group is on using quantum chemistry-based force fields for molecular dynamics simulations of organic polymers in the melt and glassy states, but the methodology is applicable to simulations of small molecules, multicomponent systems and solutions. Special attention is paid to deriving reliable descriptions of the non-bonded and electrostatic interactions. Several procedures have been developed for deriving and calibrating these parameters. Our force fields for aromatic polyimide simulations will be described. In this application, the intermolecular interactions are the critical factor in determining many properties of the polymer (including its color).

Jaffe, Richard↗
DOE OSTI2016

Molecular dynamics simulation of the structural, elastic, and thermal properties of pyrochlores

We present a comprehensive simulation study of the effect of composition on the structural, elastic and thermal properties of 25 different compounds from the pyrochlore family. We joined a repulsive potential to an existing interatomic potential to enable molecular dynamics simulations of conditions away from equilibrium. We systematically varied the chemistry of the pyrochlore by substituting different cations in the A and B sites of the A2B2O7 formula unit. The A cations varied from Lu3+ to La3+, and the B cations from Ti4+ to Ce4+. The lattice parameter increased steadily with increasing the radius of A or B cations, but the bulk modulus showed a decreasing trend with increasing cation radius. However, the specific heat capacity and thermal expansion coefficient remained almost unchanged with increasing the radii of A and B cations. It is of interest to note that Ce on the B site significantly reduces the specific heat capacity and thermal expansion coefficient, which could have implications for annealing of radiation damage in cerate pyrochlores. The present results are consistent with the experimental measurements, which validates these potentials for simulation of dynamical processes, such as radiation damage, in pyrochlores.

computer simulation↗
DOE OSTI2015

Molecular dynamic simulation of Ar-Kr mixture across a rough walled nanochannel: Velocity and temperature profiles

This paper presents the results from a molecular dynamics simulation of mixture of argon and krypton in the Poiseuille flow across a rough walled nanochannel. The roughness effect on liquid nanoflows has recently drawn attention The computational software used for carrying out the molecular dynamics simulations is LAMMPS. The fluid flow takes place between two parallel plates and is bounded by horizontal rough walls in one direction and periodic boundary conditions are imposed in the other two directions. Each fluid atom interacts with other fluid atoms and wall atoms through Leenard-Jones (LJ) potential with a cut off distance of 5.0. To derive the flow a constant force is applied whose value is varied from 0.1 to 0.3 and velocity profiles and temperature profiles are noted for these values of forces. The velocity profile and temperature profiles are also looked at different channel widths of nanochannel and at different densities of mixture. The velocity profile and temperature profile of rough walled nanochannel are compared with that of smooth walled nanochannel and it is concluded that mean velocity increases with increase in channel width, force applied and decrease in density also with introduction of roughness in the walls of nanochannel mean velocity again increases and results also agree with the analytical solution of a Poiseuille flow.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗