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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Consistent and reproducible computation of the glass transition temperature from molecular dynamics simulations

In many fields, from semiconductors for opto-electronic applications to ionic liquids (ILs) for separations, the glass transition temperature (Tg) of a material is a useful gauge for its potential use in practical settings. As a result, there is a great deal of interest in predicting Tg using molecular simulations. However, the uncertainty and variation in the trend shift method, a common approach in simulations to predict Tg, can be high. This is due to the need for human intervention in defining a fitting range for linear fits of density with temperature assumed for the liquid and glass phases across the simulated cooling. The definition of such fitting ranges then defines the estimate for the Tg as the intersection of linear fits. We eliminate this need for human intervention by leveraging the Shapiro–Wilk normality test and proposing an algorithm to define the fitting ranges and, consequently, Tg. Through this integration, we incorporate into our automated methodology that residuals must be normally distributed around zero for any fit, a requirement that must be met for any regression problem. Consequently, fitting ranges for realizing linear fits for each phase are statistically defined rather than visually inferred, obtaining an estimate for Tg without any human intervention. The method is also capable of finding multiple linear regimes across density vs temperature curves. We compare the predictions of our proposed method across multiple IL and semiconductor molecular dynamics simulation results from the literature and compare other proposed methods for automatically detecting Tg from density–temperature data. We believe that our proposed method would allow for more consistent predictions of Tg. We make this methodology available and open source through GitHub.

Chemistry↗
DOE OSTI2019

Molecular Dynamics Simulation of Hydration and Swelling of Mixed-Layer Clays in the Presence of Carbon Dioxide

Swelling of clay minerals is a prime anomaly that is relevant to many problems in geoscience, oil and gas reservoirs, geophysics, and engineering. The phenomenon is even more complex when carbon dioxide is also present. Since sequestration of CO 2 in sedimentary rock has been under active consideration as a way of alleviating the increasing concerns over climate change for which CO 2 is a main culprit, studying swelling of clays minerals in the presence of CO 2 has taken on added urgency. The swelling changes the permeability and porosity of porous formations, which in turn may lead to a reactivation of dormant fractures, triggering seismic activities, and opening up new pathways for CO 2 to return to the porous formations’ surface. Although the large majority of sedimentary rocks contain various types of mixed-layer clays (MLCs) with intermixed stacking sequences of several types of distinct layers within a single crystal, the vast majority of the past experimental and computational studies of the swelling was focused on pure clays. We present the results of extensive molecular dynamics simulation of hydration and swelling of illite-montmorillonite (I-MMT) MLCs, the most common type of mixed clays, in the presence of both water and CO 2 and various combinations of interlayer cations, Na + and K + . To understand the differences with pure clays, we also report on the same phenomenon in the MMT only. At low CO 2 concentrations in Na-MMT, which has its layers’ charge concentrated in its octahedral sheet, weak ion–surface interactions result in fully hydrated ions and, therefore, more extensive swelling than in I-MMT. Without CO 2 , adsorption of the cations at the illite surface increases the hydrophobicity of its surface. Thus, in the asymmetric interlayer of I-MMT, illite with stronger surface–ion interactions causes accumulation of cations near its surface, which limits their hydration and, therefore, controls swelling of the MLCs. Further inhibition of swelling of Na-MLC can be achieved by increasing the concentration of K + in the interlayer. At higher CO 2 concentrations, however, intercalation of water and CO 2 results in a completely different behavior, since in the 2W and 3W hydration states, the hydrophobicity of the MMT surface is stronger than that of illite. Thus, the density distributions of the intercalated molecules are considerably different from that with water only, which together with disruption of the water network by CO 2 reduces the difference between the extent of swelling of the MLC and pure MMT.

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

Dynamic Simulation with on the fly Doppler Broadening treatment in RMC

In conventional transient analysis, thermal hydraulics calculations are generally coupled to deterministic neutronics solvers. Due to the limited accuracy of deterministic methods and the rapid development of computer power, Monte Carlo neutronics calculations have received increasing attention nowadays. One of challenges for Monte Carlo methods coupled to thermal hydraulics calculations is the temperature dependent cross sections. Taking the detailed temperature distribution feedback into account, each nuclide may be in a large range of temperatures, pre-generated libraries at every temperature point are no longer applicable owing to unaffordable memory footprint. Different solutions have been proposed to deal with this problem. TMS, accounting for Doppler broadening effect by the explicit treatment of thermal motion, was used to model fast transients in Serpent 2. Stochastic mixing, using two sets of libraries for each nuclide and mixing them to get the correct temperature, was used in transient analysis by Sjenitzer and Hoogenboom. Besides, there are other on the fly Doppler broadening methods being developed, including series expansion of cross sections, multipole representation, different integration algorithms based on Doppler broadening equation. In this paper, on the fly Doppler broadening treatment based on integration algorithms was applied to transient calculation in Reactor Monte Carlo code RMC . In the following section, the integration algorithms implemented and dynamic simulation method in RMC is briefly introduced. Finally, a pure physical transient fuel assembly model was chosen for numerical validation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Molecular dynamics simulations of Zika virus NS3 helicase: Insights into RNA binding site activity

America is still suffering with the outbreak of Zika virus (ZIKV) infection. Congenital ZIKV syndrome has already caused a public health emergency of international concern. However, there are still no vaccines to prevent or drugs to treat the infection caused by ZIKV. The ZIKV NS3 helicase (NS3h) protein is a promising target for drug discovery due to its essential role in viral genome replication. NS3h unwinds the viral RNA to enable the replication of the viral genome by the NS5 protein. NS3h contains two important binding sites: the NTPase binding site and the RNA binding site. Here, we used molecular dynamics (MD) simulations to study the molecular behavior of ZIKV NS3h in the presence and absence of ssRNA and the potential implications for NS3h activity and inhibition. Although there is conformational variability and poor electron densities of the RNA binding loop in various apo flaviviruses NS3h crystallographic structures, the MD trajectories of NS3h-ssRNA demonstrated that the RNA binding loop becomes more stable when NS3h is occupied by RNA. Our results suggest that the presence of RNA generates important interactions with the RNA binding loop, and these interactions stabilize the loop sufficiently that it remains in a closed conformation. This closed conformation likely keeps the ssRNA bound to the protein for a sufficient duration to enable the unwinding/replication activities of NS3h to occur. In addition, conformational changes of this RNA binding loop can change the nature and location of the optimal ligand binding site, according to ligand binding site prediction results. These are important findings to help guide the design and discovery of new inhibitors of NS3h as promising compounds to treat the ZIKV infection. - Highlights: • Ligand binding site predictions could provide useful insights into design of compounds to block ZIKV helicase activities. • Insights into dynamic features of NS3h and its interactions with ssRNA are provided through molecular dynamics simulations. • Presence of RNA generates stabilizing interactions with RNA binding loop of NS3h, which holds it in a closed conformation. • Closed conformation of the RNA binding loop keeps the ssRNA bound to the protein to enable the activity of NS3h.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1982

Digital system for structural dynamics simulation

State-of-the-art digital hardware and software for the simulation of complex structural dynamic interactions, such as those which occur in rotating structures (engine systems). System were incorporated in a designed to use an array of processors in which the computation for each physical subelement or functional subsystem would be assigned to a single specific processor in the simulator. These node processors are microprogrammed bit-slice microcomputers which function autonomously and can communicate with each other and a central control minicomputer over parallel digital lines. Inter-processor nearest neighbor communications busses pass the constants which represent physical constraints and boundary conditions. The node processors are connected to the six nearest neighbor node processors to simulate the actual physical interface of real substructures. Computer generated finite element mesh and force models can be developed with the aid of the central control minicomputer. The control computer also oversees the animation of a graphics display system, disk-based mass storage along with the individual processing elements.

Krauter, A. I.↗
NASA NTRS1980

Analysis and improvement of contact dynamics simulation using the six degree of freedom motion system

The 6 DOF motion system was employed for closed loop simulation of the contact dynamics involved in the docking of two satellites. Under most initial conditions the system gives inaccurate results. It is well known that the presence of a time lag can generate instability in closed loop systems. In order to see if the time lag (about 200 ms.) in the 6 DOF system can can satisfactorily account for its errors, a simple computer model of the system incorporating this lag and other revelant parameters was implemented. The model in fact does predict all the observed rebound errors with reasonable precision. Also, since it is probable that the time delay can only be reduced but not eliminated, various compensation schemes were devised. Two give adequate correction for delays under 100 ms for a non preloaded probe. More powerful schemes may be adequate for the present 200 ms. environment and the presence of a preload.

Brown, R. C.↗
DOE OSTI2022

Characterization of Heparin’s Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin’s solution conformational dynamics and structure–function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule’s underlying conformational ensemble, which then provides important information necessary for understanding structure–function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin’s conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. In this work, the anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

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

Molecular Dynamic Simulations of Interaction of an AFM Probe with the Surface of an SCN Sample

Molecular dynamic (MD) simulations is conducted in order to estimate forces of probe-substrate interaction in the Atomic Force Microscope (AFM). First a review of available molecular dynamic techniques is given. Implementation of MD simulation is based on an object-oriented code developed at the University of Delft. Modeling of the sample material - succinonitrile (SCN) - is based on the Lennard-Jones potentials. For the polystyrene probe an atomic interaction potential is used. Due to object-oriented structure of the code modification of an atomic interaction potential is straight forward. Calculation of melting temperature is used for validation of the code and of the interaction potentials. Various fitting parameters of the probe-substrate interaction potentials are considered, as potentials fitted to certain properties and temperature ranges may not be reliable for the others. This research provides theoretical foundation for an interpretation of actual measurements of an interaction forces using AFM.

Bune, Adris↗
DOE OSTI2015

Utilizing fast multipole expansions for efficient and accurate quantum-classical molecular dynamics simulations

Recently, a novel approach to hybrid quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) simulations has been suggested [Schwörer et al., J. Chem. Phys. 138, 244103 (2013)]. Here, the forces acting on the atoms are calculated by grid-based density functional theory (DFT) for a solute molecule and by a polarizable molecular mechanics (PMM) force field for a large solvent environment composed of several 10{sup 3}-10{sup 5} molecules as negative gradients of a DFT/PMM hybrid Hamiltonian. The electrostatic interactions are efficiently described by a hierarchical fast multipole method (FMM). Adopting recent progress of this FMM technique [Lorenzen et al., J. Chem. Theory Comput. 10, 3244 (2014)], which particularly entails a strictly linear scaling of the computational effort with the system size, and adapting this revised FMM approach to the computation of the interactions between the DFT and PMM fragments of a simulation system, here, we show how one can further enhance the efficiency and accuracy of such DFT/PMM-MD simulations. The resulting gain of total performance, as measured for alanine dipeptide (DFT) embedded in water (PMM) by the product of the gains in efficiency and accuracy, amounts to about one order of magnitude. We also demonstrate that the jointly parallelized implementation of the DFT and PMM-MD parts of the computation enables the efficient use of high-performance computing systems. The associated software is available online.

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

Review of the fundamental theories behind small angle X-ray scattering, molecular dynamics simulations, and relevant integrated application

In this paper, the fundamental concepts and equations necessary for performing small angle X-ray scattering (SAXS) experiments, molecular dynamics (MD) simulations, and MDSAXS analyses were reviewed. Furthermore, several key biological and non-biological applications for SAXS, MD, and MD-SAXS are presented in this review; however, this article does not cover all possible applications. SAXS is an experimental technique used for the analysis of a wide variety of biological and non-biological structures. SAXS utilizes spherical averaging to produce one- or two-dimensional intensity profiles, from which structural data may be extracted. MD simulation is a computer simulation technique that is used to model complex biological and non-biological systems at the atomic level. MD simulations apply classical Newtonian mechanics’ equations of motion to perform force calculations and to predict the theoretical physical properties of the system. This review presents several applications that highlight the ability of both SAXS and MD to study protein folding and function in addition to non-biological applications, such as the study of mechanical, electrical, and structural properties of non-biological nanoparticles. Lastly, the potential benefits of combining SAXS and MD simulations for the study of both biological and non-biological systems are demonstrated through the presentation of several examples that combine the two techniques.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Molecular dynamics simulations of NMR relaxation and diffusion of bulk hydrocarbons and water

Molecular dynamics (MD) simulations are used to investigate 1 H nuclear magnetic resonance (NMR) relaxation and diffusion of bulk n-C 5 H 12 to n-C 17 H 36 hydrocarbons and bulk water. The MD simulations of the 1H NMR relaxation times T 1,2 in the fast motion regime where T 1 = T 2 agree with measured (de-oxygenated) T2 data at ambient conditions, without any adjustable parameters in the interpretation of the simulation data. Likewise, the translational diffusion DT coefficients calculated using simulation con gurations are well-correlated with measured diffusion data at ambient conditions. The agreement between the predicted and experimentally measured NMR relaxation times and diffusion coefficient also validate the force elds used in the simulation. The molecular simulations naturally separate intramolecular from intermolecular dipole-dipole interactions helping bring new insight into the two NMR relaxation mechanisms as a function of molecular chain-length (i.e. carbon number). Comparison of the MD simulation results of the two relaxation mechanisms with traditional hard-sphere models used in interpreting NMR data reveals important limitations in the latter. With increasing chain length, there is substantial deviation in the molecular size inferred on the basis of the radius of gyration from simulation and the tted hard-sphere radii required to rationalize the relaxation times. This deviation is characteristic of the local nature of the NMR measurement, one that is well-captured by molecular simulations.

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

Molecular dynamics simulations of grain boundary thermal resistance in UO2

By means of molecular dynamics (MD) simulations, we have calculated Kaptiza resistance of UO2 with or without radiation damage. For coincident site lattice boundaries of different configurations, the boundary thermal resistance of unirradiated UO2 can be well described by a parameter-reduced formula by using boundary energies as variables. We extended the study to defect-loaded UO2 by introducing damage cascades in close vicinity to the boundaries. Following cascade annealing and defect migrations towards grain boundaries, the boundary energy increases and so does Kaptiza resistance. The correlations between these two still follow the same formula extracted from the unirradiated UO2. The finding will benefit multi-scale modeling of UO2 thermal properties under extreme radiation conditions by combining effects from boundary configurations and damage levels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

Biomolecular Solvation Structure Revealed by Molecular Dynamics Simulations

To compare ordered water positions from experiment with those from molecular dynamics (MD) simulations, a number of MD models of water structure in crystalline endoglucanase were calculated. The starting MD model was derived from a joint X-ray and neutron diffraction crystal structure, enabling the use of experimentally assigned protonation states. Simulations were performed in the crystalline state, using a periodic 2 × 2 × 2 supercell with explicit solvent. Water X-ray and neutron scattering density maps were computed from MD trajectories using standard macromolecular crystallography methods. In one set of simulations, harmonic restraints were applied to bias the protein structure toward the crystal structure. For these simulations, the recall of crystallographic waters using strong peaks in the MD water electron density was very good, and there also was substantial visual agreement between the boomerang-like wings of the neutron scattering density and the crystalline water hydrogen positions. An unrestrained simulation also was performed. For this simulation, the recall of crystallographic waters was much lower. For both restrained and unrestrained simulations, the strongest water density peaks were associated with crystallographic waters. Overall, the results demonstrate that it is now possible to recover crystallographic water structure using restrained MD simulations but that it is not yet reasonable to expect unrestrained MD simulations to do the same. Further development and generalization of MD water models for force-field development, macromolecular crystallography, and medicinal chemistry applications is now warranted. In particular, the combination of room-temperature crystallography, neutron diffraction, and crystalline MD simulations promises to substantially advance modeling of biomolecular solvation.

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

Dynamic Simulation of a Wave Rotor Topped Turboshaft Engine

The dynamic behavior of a wave rotor topped turboshaft engine is examined using a numerical simulation. The simulation utilizes an explicit, one-dimensional, multi-passage, CFD based wave rotor code in combination with an implicit, one-dimensional, component level dynamic engine simulation code. Transient responses to rapid fuel flow rate changes and compressor inlet pressure changes are simulated and compared with those of a similarly sized, untopped, turboshaft engine. Results indicate that the wave rotor topped engine responds in a stable, and rapid manner. Furthermore, during certain transient operations, the wave rotor actually tends to enhance engine stability. In particular, there is no tendency toward surge in the compressor of the wave rotor topped engine during rapid acceleration. In fact, the compressor actually moves slightly away from the surge line during this transient. This behavior is precisely the opposite to that of an untopped engine. The simulation is described. Issues associated with integrating CFD and component level codes are discussed. Results from several transient simulations are presented and discussed.

Greendyke, R. B.↗
NASA NTRS1989

An order (n) algorithm for the dynamics simulation of robotic systems

The formulation of an Order (n) algorithm for DISCOS (Dynamics Interaction Simulation of Controls and Structures), which is an industry-standard software package for simulation and analysis of flexible multibody systems is presented. For systems involving many bodies, the new Order (n) version of DISCOS is much faster than the current version. Results of the experimental validation of the dynamics software are also presented. The experiment is carried out on a seven-joint robot arm at NASA's Goddard Space Flight Center. The algorithm used in the current version of DISCOS requires the inverse of a matrix whose dimension is equal to the number of constraints in the system. Generally, the number of constraints in a system is roughly proportional to the number of bodies in the system, and matrix inversion requires O(p exp 3) operations, where p is the dimension of the matrix. The current version of DISCOS is therefore considered an Order (n exp 3) algorithm. In contrast, the Order (n) algorithm requires inversion of matrices which are small, and the number of matrices to be inverted increases only linearly with the number of bodies. The newly-developed Order (n) DISCOS is currently capable of handling chain and tree topologies as well as multiple closed loops. Continuing development will extend the capability of the software to deal with typical robotics applications such as put-and-place, multi-arm hand-off and surface sliding.

Chun, H. M.↗
DOE OSTI2014

Accelerating dissipative particle dynamics simulations for soft matter systems

Dissipative particle dynamics (DPD) is a coarse-grained particle-based simulation method that offers microscopic-scale insights into soft matter systems. We present an efficient implementation of a DPD model for graphical processing units (GPUs). As implemented in the LAMMPS molecular dynamics package, it can run effectively on current-generation supercomputers which often have hybrid nodes containing multi-core CPUs and (one or more) GPUs. Using efficient communication of information between the CPUs and GPUs, DPD interactions can be computed on the GPU while other portions of a full simulation model (boundary conditions, constraints, bonded interactions, diagnostic calculations, etc.) can be performed on the CPU. Our GPU-enhanced runs show a speedup of up to 9.5× versus many-core CPU simulations, and can run scalably across thousands of compute nodes. We briefly discuss how the new GPU implementation was validated against the CPU version for thermodynamics, diffusion, and hydrodynamic behavior. We also highlight large-scale models which the faster DPD implementation has enabled, for studies of monolayer self-assembly and thin-film instabilities.

36 MATERIALS SCIENCE↗
NASA NTRS2017

Modeling Chemical Reactions in Classical Molecular Dynamics Simulations

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

Jacob R Gissinger↗