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

Dissipative particle dynamics simulation of dilute polymer solutions—Inertial effects and hydrodynamic interactions

We examine the accuracy of dissipative particle dynamics (DPD) simulations of polymers in dilute solutions with hydrodynamic interaction (HI), at the theta point, modeled by setting the DPD conservative interaction between beads to zero. We compare the first normal-mode relaxation time extracted from the DPD simulations with theoretical predictions from a normal-mode analysis for theta chains. We characterize the influence of bead inertia within the coil by a ratio L{sub m}/R{sub g}, where L{sub m} is the ballistic distance over which bead inertia is lost, and R{sub g} is the radius of gyration of the polymer coil, while the HI strength per bead h* is determined by the ratio of bead hydrodynamic radius (r{sub H}) to the equilibrium spring length. We show how to adjust h* through the spring length and monomer mass, and how to optimize the accuracy of DPD for fixed h* by increasing the friction coefficient (γ ≥ 9) and by incorporating a nonlinear distance dependence into the frictional interaction. Even with this optimization, DPD simulations exhibit deviations of over 20% from the theoretical normal-mode predictions for high HI strength with h* ≥ 0.20, for chains with as many as 100 beads, which is a larger deviation than is found for Stochastic rotation dynamics simulations for similar chains lengths and values of h*.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Constant-potential molecular dynamics simulations of molten salt double layers for FLiBe and FLiNaK

We report the results of constant-potential molecular dynamics simulations of the double layer interface between molten 2LiF–BeF 2 (FLiBe) and 23LiF–6NaF–21KF (FLiNaK) fluoride mixtures and idealized solid electrodes. Employing methods similar to those used in studies of chloride double layers, we compute the structure and differential capacitance of molten fluoride electric double layers as a function of applied voltage. The role of molten salt structure is probed through comparisons between FLiBe and FLiNaK, which serve as models for strong and weak associate-forming salts, respectively. In FLiBe, screening involves changes in Be–F–Be angles and alignment of the oligomers parallel to the electrode, while in FLiNaK, the electric field is screened mainly by rearrangement of individual ions, predominantly the polarizable potassium cation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Fracture Prediction of Epoxy Resin using Morse Bond Potential Embedded in GAFF using Molecular Dynamics Simulations

A generic method is developed to investigate the fracture behavior of an epoxy resin using atomistic molecular (MD) dynamics simulations. The epoxy system consisted of the stoichiometric mixture of tetra-/tri-functional epoxies and di-functional hardener molecules and was cured with a step-growth crosslinking algorithm. A hybrid force field (hFF) in which Morse bond function were added in second generation-general Amber Force Field (GAFF2) was used to capture bond breaking of the crosslinked epoxy system under uniaxial deformation. The Morse bond parameters for covalent bonds in the crosslinked backbone were fit to dissociation curves computed by CASPT2/6-311+G**. After the systems in various crosslink densities were obtained and equilibrated, the full length of stress-strain (𝜎-𝜀) curve was produced including initial elastic regime, yielding, plastic flow, strain hardening, and progressive failure. We found that the hFF was effective in studying the dissociation of crosslinked polymers and understanding the failure mechanism. In the results, the hFF delivered the same elastic property obtained from the unmodified force field GAFF. The stress-strain curve, however, began to be distinct in the plastic regime due to broken covalent bonds capable of altering the tensile behavior. Once the stress reached maximum, it progressively decreased until the ultimate failure (𝜎=0). Since hFF employed the existing GAFF functions with including a minimal change of its bond potential, this approach is easier to perform and computationally more efficient for unveiling the fracture behavior of polymer materials at the molecular level than other approaches with reactive force fields. We expect this approach will be utilized to accelerate the material-by-design process for thermosets by incorporate data from molecular models.

Changwoon Jang↗
DOE OSTI2025

Accelerating charge estimation in molecular dynamics simulations using physics-informed neural networks: corrosion applications

Molecular Dynamics (MD) simulations are used to understand the effects of corrosion on metallic materials in salt brine. Reactive force fields in classical MD enable accurate modeling of bond formation and breakage in the aqueous medium and at the metal-electrolyte interface, while also facilitating dynamic partial charge equilibration. However, MD simulations are computationally intensive and unsuitable for modeling the long time scales characteristic of corrosive phenomena. To address this, we develop reduced-order machine learning models that provide accurate and efficient predictions of charge density in corrosive environments. Specifically, we use Long Short-Term Memory (LSTM) networks to forecast charge density evolution based on atomic environments represented by Smooth Overlap of Atomic Positions (SOAP) descriptors. A physics-informed loss function enforces charge neutrality and electronegativity equivalence. The atomic charges predicted by the deep learning model trained on this work were obtained two orders of magnitude faster than those from molecular dynamics (MD) simulations, with an error of less than 3% compared to the MD-obtained charges, even in extrapolative scenarios, while adhering to physical constraints. This demonstrates the excellent accuracy, computational efficiency, and validity of the developed model. Lastly, even though developed for corrosion, these protocols are formulated in a phenomenon-agnostic manner, allowing application to various variable-charge interatomic potentials and related fields.

Atomistic models↗
DOE OSTI2018

Molecular Dynamics Simulation of the Structure and Ion Transport in the Ce{sub 1–x}Gd{sub x}O{sub 2–δ}|YSZ Heterosystem

Molecular dynamics simulation has been used to develop a realistic atomistic model of two-layer Ce{sub 1–x}Gd{sub x}O{sub 2–δ}|YSZ heterosystem. It is shown that Ce{sub 1–x}Gd{sub x}O{sub 2–δ} and YSZ layers (about 15 and 16 Å thick, respectively) retain their crystal structure on the whole. The main structural distortions are found to occur near the Ce{sub 1–x}Gd{sub x}O{sub 2–δ}|YSZ geometric interface, within a narrow interfacial region of few angstroms thick. Both the generalized diffusion characteristics of the system as a whole and the oxygen diffusion coefficients in the layers are calculated, and the diffusion activation energies are determined.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Molecular dynamics simulation studies of the wild type and E92Q/N155H mutant of Elvitegravir-resistance HIV-1 integrase

Although Elvitegravir (EVG) is a newly developed antiretrovirals drug to treat the acquired immunodeficiency syndrome (AIDS), drug resistance has already been found in clinic, such as E92Q/N155H and Q148H/G140S. Several structural investigations have already been reported to reveal the molecular mechanism of the drug resistance. As full length crystal structure for HIV-1 integrase is still unsolved, we use in this paper the crystal structure of the full length prototype foamy virus (PFV) in complex with virus DNA and inhibitor Elvitegravir as a template to construct the wild type and E92Q/N155H mutant system of HIV-1 integrase. Molecular dynamic simulations was used to revel the binding mode and the drug resistance of the EVG ligand in E92Q/N155H. Several important interactions were discovered between the mutated residues and the residues in the active site of the E92Q/N155H double mutant pattern, and cross correlation and clustering methods were used for detailed analysis. The results from the MD simulation studies will be used to guide the experimental efforts of developing novel inhibitors against drug-resistant HIV integrase mutants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Alumina(0001)/Water Interface: Structural Properties and Infrared Spectra from First-Principles Molecular Dynamics Simulations

We investigated the atomic structure and infrared spectra of the alumina(0001)/water interface, using first-principles molecular dynamics simulations based on density functional theory within the generalized gradient approximation. The computed structural properties of the interface are in good agreement with the results of synchrotron X-ray experiments. In this work detailed analyses of the computed infrared spectra revealed two types of water molecules at the hydrophilic oxide/water interface: molecules participating in strong “ice-like” hydrogen bonding with the oxide surface and molecules involved in weaker “liquid-like” hydrogen bonding. Our results provide a molecular interpretation of the “ice-like” and “liquid-like” bands observed in sum-frequency vibrational spectroscopy experiments and underscore the significance of strong hydrogen-bonding interactions in determining the orientation of interfacial water molecules.

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

Spectral shapes of Ar-broadened HCl lines in the fundamental band by classical molecular dynamics simulations and comparison with experiments

Spectral shapes of isolated lines of HCl perturbed by Ar are investigated for the first time using classical molecular dynamics simulations (CMDS). Using reliable intermolecular potentials taken from the literature, these CMDS provide the time evolution of the auto-correlation function of the dipole moment, whose Fourier-Laplace transform leads to the absorption spectrum. In order to test these calculations, room temperature spectra of various lines in the fundamental band of HCl diluted in Ar are measured, in a large pressure range, with a difference-frequency laser spectrometer. Comparisons between measured and calculated spectra show that the CMDS are able to predict the large Dicke narrowing effect on the shape of HCl lines and to satisfactorily reproduce the shapes of HCl spectra at different pressures and for various rotational quantum numbers.

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

Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations

Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.

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

Insights on the structure and properties of sodium iron phosphate glasses from molecular dynamics simulations

Iron phosphate glasses are promising nuclear waste forms while more detailed understanding of their structures and structure-property relations are still needed to better design waste glass compositions. In this work we report studies of three series of sodium iron phosphate (NFP) glasses: 60P 2 O 5 -(40-x)Fe 2 O 3 -xNa 2 O (x = 0→40), (100–2x)P 2 O 5 -xFe 2 O 3 -xNa 2 O (x = 5→17.5) and one with different iron redox ratio, to understand the composition as well as the iron redox effects on the structure and properties of these glasses using molecular dynamics simulations with effective two-body and three-body potentials. Structural analyses, including pair distribution function, bond angle distribution, Q n distribution, and polyhedral connectivity, were performed to obtain in-depth information on short-range and medium-range structural features. The P-O pair distributions showed a first peak splitting with phosphorus-bridging and non-bridging oxygen contributions. This and the average P-O and other cation-oxygen bond distances are in excellent agreement with experiments. The coordination number of P 5+ remained four while that of Fe 3+ increased from 4.30 to 4.72 with decreasing Fe/Na ratio. Polyhedral linkage analysis showed [PO 4 ] units linked with [PO 4 ] and [FeO x ] through corner-sharing while the [PO 4 ]-[FeO x ] linkages become dominant for compositions with Fe 2 O 3 larger than 15 mol%. The effect of iron redox ratio on the structure of NFP glasses was also studied and it was found that bond lengths and coordination numbers were not strongly affected, while the reduction of iron introduced higher network distortions, as evident by O-P-O bond angle and Q n distribution. The glass transition temperature (T g ) showed a monotonic increase with Fe 2 O 3 in the first series, in good agreement with experiments, while those of the second series showed a maximum at P 2 O 5 = 82 mol%. Here, calculated elastic moduli were found to increase with Fe 2 O 3 in the first glass series, which was be explained by the increase of network connectivity, while those of the second series decrease with Fe 2 O 3 due to decrease of P 2 O 5 .

36 MATERIALS SCIENCE↗
DOE OSTI2021

Molecular Dynamics Simulations of a Catalytic Multivalent Peptide–Nanoparticle Complex

Molecular modeling of a supramolecular catalytic system is conducted resulting from the assembling between a small peptide and the surface of cationic self-assembled monolayers on gold nanoparticles, through a multiscale iterative approach including atomistic force field development, flexible docking with Brownian Dynamics and µs-long Molecular Dynamics simulations. Self-assembly is a prerequisite for the catalysis, since the catalytic peptides do not display any activity in the absence of the gold nanocluster. Atomistic simulations reveal details of the association dynamics as regulated by defined conformational changes of the peptide due to peptide length and sequence. Our results show the importance of a rational design of the peptide to enhance the catalytic activity of peptide–nanoparticle conjugates and present a viable computational approach toward the design of enzyme mimics having a complex structure–function relationship, for technological and nanomedical applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

MixPI: Mixed-time slicing path integral software for quantized molecular dynamics simulations

We introduce the MixPI software to implement path integral molecular dynamics (PIMD) simulations for the study of condensed phase systems where nuclear quantum effects (NQEs) are important. In contrast to existing PIMD simulation software, MixPI enables the implementation of mixed quantum–classical path integral simulations where only a subset of system degrees of freedom (dofs) are treated quantum mechanically in an extended phase space while the remaining dofs are described classically. We expect this software to be particularly useful for simulations of electron and proton transfer in condensed phase systems, as well as for the study of biological and material systems where only a handful of dofs contribute significantly to the observed NQEs. We demonstrate the use of MixPI in two different systems. The first is a simple water model where we implement a set of mixed quantum–classical simulations to compute average energy and radial distribution functions. We use these simulations to benchmark the effectiveness of MixPI and to demonstrate how it enables systematic investigation into the origin of observed NQEs. We then compute radial distribution functions for a system where MixPI is essential: a solvated metal (M 2+ ) cation described using an explicit quantized electron localized on an M 3+ ion in water.

chemical physics↗
DOE OSTI2022

A high-dimensional neural network potential for molecular dynamics simulations of condensed phase nickel and phase transitions

A high-dimensional neural network interatomic potential was developed and used in molecular dynamics simulations of condensed phase Ni and Ni systems with liquid–solid phase coexistence. The reference data set was generated by sampling the potential energy surface over a broad temperature-pressure domain using ab initio MD simulations to train a unified potential. Excellent agreement was achieved between bulk face-centered cubic nickel thermal expansion simulations and relevant experimental data. The same potential also yields accurate structures and diffusivities in the liquid state. The phase transition between liquid and solid phases was simulated using the two-phase interface method. The predicted melting point temperature is within a few kelvins of the literature value. Here, the general methodology could be applied to describe crystals with much more complex phase behaviors.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

The Potential of Computational Fluid Dynamics Simulation on Serial Monitoring of Hemodynamic Change in Type B Aortic Dissection

PurposeWe aimed to assess the potential of computational fluid dynamics simulation (CFD) in detecting changes in pressure and flow velocity in response to morphological changes in type B aortic dissection.Materials and MethodsPressure and velocity in four morphological models of type B aortic dissection before and after closure of the entry tear were calculated with CFD and analyzed for changes among the different scenarios. The control model (Model 1) was patient specific and built from the DICOM data of CTA, which bore one entry tear and three re-entry tears. Models 2–4 were modifications of Model 1, with two re-entry tears less in Model 2, one re-entry tear more in Model 3, and a larger entry tear in Model 4.ResultsThe pressure and velocity pertaining to each of the morphological models were unique. Changes in pressure and velocity findings were accountable by the changes in morphological features of the different models. There was no blood flow in the false lumen across the entry tear after its closure, the blood flow direction across the re-entry tears was reversed after closure of the entry tear.ConclusionCFD simulation is probably useful to detect hemodynamic changes in the true and false lumens of type B aortic dissection in response to morphological changes, it may potentially be developed into a non-invasive and patient-specific tool for serial monitoring of hemodynamic changes of type B aortic dissection before and after treatment.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2016

Modelling the local atomic structure of molybdenum in nuclear waste glasses with ab initio molecular dynamics simulations

The nature of chemical bonding of molybdenum in high level nuclear waste glasses has been elucidated by ab initio molecular dynamics simulations. Two compositions, (SiO2)57.5 – (B2O3)10 – (Na2O)15 – (CaO)15 – (MoO3)2.5 and (SiO2)57.3 – (B2O3)20 – (Na2O)6.8 – (Li2O)13.4 – (MoO3)2.5 , were considered in order to investigate the effect of ionic and covalent components on the glass structure and the formation of the crystallisation precursors (Na2MoO4 and CaMoO4). The coordination environments of Mo cations and the corresponding bond lengths calculated from our model are in excellent agreement with experimental observations. The analysis of the first coordination shell reveals two different types of molybdenum host matrix bonds in the lithium sodium borosilicate glass. Based on the structural data and the bond valence model, we demonstrate that the Mo cation can be found in a redox state and the molybdate tetrahedron can be connected with the borosilicate network in a way that inhibits the formation of crystalline molybdates. These results significantly extend our understanding of bonding in Mo-containing nuclear waste glasses and demonstrate that tailoring the glass composition to specific heavy metal constituents can facilitate incorporation of heavy metals at high concentrations. K.K. was supported through the Impact Studentship scheme at UCL co-funded by the IHI Corporation and UCL. P.V.S. thanks the Royal Society, which supported preliminary work on this project, and the Laboratory Directed Research and Development program at PNNL, a multiprogram national laboratory operated by Battelle for the U.S. Department of Energy. Via our membership of the UK's HEC Materials Chemistry Consortium, which is funded by EPSRC (EP/L000202), this work used the ARCHER UK National Supercomputing Service (http://www.archer.ac.uk).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗