Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic simulations”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Helium flux effects on bubble growth and surface morphology in plasma-facing tungsten from large-scale molecular dynamics simulations

We investigate helium flux effects on helium transport and surface evolution in plasma-facing tungsten using molecular dynamics. The simulations span two orders of magnitude, from ITER-relevant levels to those more typical of simulations published to date. Simulation times of up to 2.5 µs (corresponding to actual fluences of 1.4 x 10 21 m -2 ) are achieved, revealing concerted bubble-bursting events that are responsible for significant and very sudden changes in surface morphology. The depth distribution of helium depends very strongly on helium flux: helium self-trapping becomes more probable near the surface at high flux, and a layer of near-surface bubbles forms. Helium retention prior to the onset of bubble bursting is also substantially lower at low flux than it is at high flux. Surface features at low fluence are correlated with the positions of bubbles, but at high fluence, bubbles tend to coalesce, venting to the surface at one or more locations and leaving large interconnected cavities below the surface. Ruptured bubbles may serve as pathways deeper into the material, allowing helium to bypass the layer of near-surface bubbles and fill deeper, potentially much larger, bubbles that can produce more substantial surface features. Deeper bubbles also emit prismatic dislocation loops that can fill in cavities closer to the surface. Our results suggest that nearly all molecular dynamics simulations published to date are hampered by finite-size effects, and that helium flux is a very important parameter in determining the behavior of helium in plasma-facing components.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Molecular dynamics simulations of graphene oxide frameworks

We use quantum mechanical calculations to develop a full set of force field parameters in order to perform molecular dynamics simulations to understand and optimize the molecular storage properties inside Graphene Oxide Frameworks (GOFs). A set of boron-related parameters for commonly used empirical force fields is determined to describe the non-bonded and bonded interactions between linear boronic acid linkers and graphene sheets of GOF materials. The transferability of the parameters is discussed and their validity is quantified by comparing quantum mechanical and molecular mechanical structural and vibrational properties. In conclusion, the application of the model to the dynamics of water inside the GOFs reveals significant variations in structural flexibility of GOF depending on the linker density, which is shown to be usable as a tuning parameter for desired diffusion properties.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Molecular dynamics simulations of the structure and single-particle dynamics of mixtures of divalent salts and ionic liquids

We report a molecular dynamics study of the structure and single-particle dynamics of mixtures of a protic (ethylammonium nitrate) and an aprotic (1-butyl-3-methylimidazolium hexaflurophosphate [BMIM][PF{sub 6}]) room-temperature ionic liquids doped with magnesium and calcium salts with a common anion at 298.15 K and 1 atm. The solvation of these divalent cations in dense ionic environments is analyzed by means of apparent molar volumes of the mixtures, radial distribution functions, and coordination numbers. For the protic mixtures, the effect of salt concentration on the network of hydrogen bonds is also considered. Moreover, single-particle dynamics of the salt cations is studied by means of their velocity autocorrelation functions and vibrational densities of states, explicitly analyzing the influence of salt concentration, and cation charge and mass on these magnitudes. The effect of the valency of the salt cation on these properties is considered comparing the results with those for the corresponding mixtures with lithium salts. We found that the main structural and dynamic features of the local solvation of divalent cations in ionic liquids are similar to those of monovalent salts, with cations being localized in the polar nanoregions of the bulk mixture coordinated in monodentate and bidentate coordination modes by the [NO{sub 3}]{sup −} and [PF{sub 6}]{sup −} anions. However, stronger electrostatic correlations of these polar nanoregions than in mixtures with salts with monovalent cations are found. The vibrational modes of the ionic liquid (IL) are seen to be scarcely affected by the addition of the salt, and the effect of mass and charge on the vibrational densities of states of the dissolved cations is reported. Cation mass is seen to exert a deeper influence than charge on the low-frequency vibrational spectra, giving a red shift of the vibrational modes and a virtual suppression of the higher energy vibrational modes for the heavier Ca{sup 2+} cations. No qualitative difference with monovalent cations was found in what solvation is concerned, which suggests that no enhanced reduction of the mobility of these cations and their complexes in ILs respective to those of monovalent cations is to be expected.

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

Collective Vibrational Strong Coupling Effects on Molecular Vibrational Relaxation and Energy Transfer: Numerical Insights via Cavity Molecular Dynamics Simulations**

Abstract For a small fraction of hot CO 2 molecules immersed in a liquid‐phase CO 2 thermal bath, classical cavity molecular dynamics simulations show that forming collective vibrational strong coupling (VSC) between the C=O asymmetric stretch of CO 2 molecules and a cavity mode accelerates hot‐molecule relaxation. This acceleration stems from the fact that polaritons can be transiently excited during the nonequilibrium process, which facilitates intermolecular vibrational energy transfer. The VSC effects on these rates 1) resonantly depend on the cavity mode detuning, 2) cooperatively depend on Rabi splitting, and 3) collectively scale with the number of hot molecules. For larger cavity volumes, the average VSC effect per molecule can remain meaningful for up to N ≈10 4 molecules forming VSC. Moreover, the transiently excited lower polariton prefers to relax by transferring its energy to the tail of the molecular energy distribution rather than distributing it equally to all thermal molecules. As far as the parameter dependence is concerned, the vibrational relaxation data presented here appear analogous to VSC catalysis in Fabry–Pérot microcavities.

Li, Tao E.↗
DOE OSTI2021

Collective Vibrational Strong Coupling Effects on Molecular Vibrational Relaxation and Energy Transfer: Numerical Insights via Cavity Molecular Dynamics Simulations**

Abstract For a small fraction of hot CO 2 molecules immersed in a liquid‐phase CO 2 thermal bath, classical cavity molecular dynamics simulations show that forming collective vibrational strong coupling (VSC) between the C=O asymmetric stretch of CO 2 molecules and a cavity mode accelerates hot‐molecule relaxation. This acceleration stems from the fact that polaritons can be transiently excited during the nonequilibrium process, which facilitates intermolecular vibrational energy transfer. The VSC effects on these rates 1) resonantly depend on the cavity mode detuning, 2) cooperatively depend on Rabi splitting, and 3) collectively scale with the number of hot molecules. For larger cavity volumes, the average VSC effect per molecule can remain meaningful for up to N ≈10 4 molecules forming VSC. Moreover, the transiently excited lower polariton prefers to relax by transferring its energy to the tail of the molecular energy distribution rather than distributing it equally to all thermal molecules. As far as the parameter dependence is concerned, the vibrational relaxation data presented here appear analogous to VSC catalysis in Fabry–Pérot microcavities.

Li, Tao E.↗
DOE OSTI2019

Transfer of an erbium ion across the water/dodecane interface: Structure and thermodynamics via molecular dynamics simulations

The thermodynamics and structural changes involved in the transfer of the Er 3+ ion across the water/dodecane interface are investigated by molecular dynamics simulations. We show that the Er 3+ ion is transferred as a highly conserved 8-water coordinated species and that the transfer involves significant perturbation of the interfacial water structure. Furthermore, several structural properties are used to quantify this process. Implications for the ion extraction process are discussed.

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

Distance Restraints from Crosslinking Mass Spectrometry: Mining a Molecular Dynamics Simulation Database to Evaluate Lysine-Lysine Distances

Integrative structural biology models the structures of protein complexes that are intractable by classical structural methods (because of extreme size, dynamics, or heterogeneity) by combining computational structural modeling with data from experimental methods. One such method is chemical cross-linking mass spectrometry (XL-MS), in which cross-linked peptides, derived from a covalently cross-linked protein complex and identified by liquid chromatography-mass spectrometry, pinpoint protein residues close in three-dimensional space. The commonly used lysine-reactive N-hydroxysuccinimide ester reagents disuccinimidylsuberate (DSS) and bis(sulfosuccinimidyl)suberate (BS3) have a linker arm that is 11.4 Å long when fully extended. However, XL-MS studies on proteins of known structure frequently report cross-links that exceed this distance. Typically, a tolerance of ~3 Å is added to the theoretical maximum to account for this observation, with little justification for the value chosen. We used the Dynameomics database, a repository of high-quality molecular dynamics simulations of 807 proteins representative of all protein folds, to investigate the change in lysine-lysine distances resulting from native-state dynamics on the time-scale of tens of nanoseconds. We conclude that observed cross-links are consistent with a protein structure if the distance between cross-linked lysine Nζ atoms is less than the cross-linker length plus 11.3 Å. For DSS or BS3, this corresponds to a Cα to Cα distance of 30.4 Å. This analysis provides a theoretical basis for the widespread practice of adding a tolerance to the crosslinker length when comparing XL-MS results to structures, and indicates the appropriate values of an XLMS derived distance constraint to use in structural modeling.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2023

Molecular-dynamics simulation methods for macromolecular crystallography

It is investigated whether molecular-dynamics (MD) simulations can be used to enhance macromolecular crystallography (MX) studies. Historically, protein crystal structures have been described using a single set of atomic coordinates. Because conformational variation is important for protein function, researchers now often build models that contain multiple structures. Methods for building such models can fail, however, in regions where the crystallographic density is difficult to interpret, for example at the protein–solvent interface. To address this limitation, a set of MD–MX methods that combine MD simulations of protein crystals with conventional modeling and refinement tools have been developed. In an application to a cyclic adenosine monophosphate-dependent protein kinase at room temperature, the procedure improved the interpretation of ambiguous density, yielding an alternative water model and a revised protein model including multiple conformations. The revised model provides mechanistic insights into the catalytic and regulatory interactions of the enzyme. The same methods may be used in other MX studies to seek mechanistic insights.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Direct comparison of elastic incoherent neutron scattering experiments with molecular dynamics simulations of DMPC phase transitions

Neutron scattering techniques have been employed to investigate 1,2-dimyristoyl-sn-glycero-3- phosphocholine (DMPC) membranes in the form of multilamellar vesicles (MLVs) and deposited, stacked multilamellar-bilayers (MLBs), covering transitions from the gel to the liquid phase. Neutron diffraction was used to characterise the samples in terms of transition temperatures, whereas elastic incoherent neutron scattering (EINS) demonstrates that the dynamics on the sub-macromolecular length-scale and pico- to nano-second time-scale are correlated with the structural transitions through a discontinuity in the observed elastic intensities and the derived mean square displacements. Molecular dynamics simulations have been performed in parallel focussing on the length-, time- and temperature-scales of the neutron experiments. They correctly reproduce the structural features of the main gel-liquid phase transition. Particular emphasis is placed on the dynamical amplitudes derived from experiment and simulations. Two methods are used to analyse the experimental data and mean square displacements. They agree within a factor of 2 irrespective of the probed time-scale, i.e. the instrument utilized. Mean square displacements computed from simulations show a comparable level of agreement with the experimental values, albeit, the best match with the two methods varies for the two instruments. Consequently, experiments and simulations together give a consistent picture of the structural and dynamical aspects of the main lipid transition and provide a basis for future, theoretical modelling of dynamics and phase behaviour in membranes. The need for more detailed analytical models is pointed out by the remaining variation of the dynamical amplitudes derived in two different ways from experiments on the one hand and simulations on the other.

Aoun, Bachir↗
DOE OSTI2017

Reference Computational Meshing Strategy for Computational Fluid Dynamics Simulation of Departure from Nucleate Boiling

The objective of this effort is to establish a strategy and process for generation of suitable computational mesh for computational fluid dynamics simulations of departure from nucleate boiling in a 5 by 5 fuel rod assembly held in place by PWR mixing vane spacer grids. This mesh generation process will support ongoing efforts to develop, demonstrate and validate advanced multi-phase computational fluid dynamics methods that enable more robust identification of dryout conditions and DNB occurrence.Building upon prior efforts and experience, multiple computational meshes were developed using the native mesh generation capabilities of the commercial CFD code STAR-CCM+. These meshes were used to simulate two test cases from the Westinghouse 5 by 5 rod bundle facility. The sensitivity of predicted quantities of interest to the mesh resolution was then established using two evaluation methods, the Grid Convergence Index method and the Least Squares method. This evaluation suggests that the Least Squares method can reliably establish the uncertainty associated with local parameters such as vector velocity components at a point in the domain or surface averaged quantities such as outlet velocity magnitude. However, neither method is suitable for characterization of uncertainty in global extrema such as peak fuel surface temperature, primarily because such parameters are not necessarily associated with a fixed point in space. This shortcoming is significant because the current generation algorithm for identification of DNB event conditions relies on identification of such global extrema. Ongoing efforts to identify DNB based on local surface conditions will address this challenge.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Molecular dynamics simulation of seeded crystal growth in glass

The mechanism of non-congruent growth of a crystal from glass has been sought using molecular dynamics simulations. Specifically, as a model of this process, the growth of a lithium niobate (LiNbO 3 ) crystal seed sandwiched between two lithium niobosilicate (LNS) glass slabs has been simulated as a function of time and temperature. The growth of pre-existing crystal is strongly affected by the orientation of crystal seed, temperature, and the SiO 2 concentration in the surrounding LNS glass matrix. The orientation of LiNbO 3 seed surface that has inherently larger interplanar distance results in a relatively slower crystal growth. The addition of SiO 2 to LNS system significantly decreases the crystal growth, which primarily occurs in the region devoid of Si. The suppressive effect of SiO 2 on growth rate can be traced to the existence of defect complex comprising of Si substituted at the Nb site and a nearby Nb vacancy.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Parallel processing for nonlinear dynamics simulations of structures including rotating bladed-disk assemblies

The principal objective of this research is to develop, test, and implement coarse-grained, parallel-processing strategies for nonlinear dynamic simulations of practical structural problems. There are contributions to four main areas: finite element modeling and analysis of rotational dynamics, numerical algorithms for parallel nonlinear solutions, automatic partitioning techniques to effect load-balancing among processors, and an integrated parallel analysis system.

Hsieh, Shang-Hsien↗
DOE OSTI2014

The FADE mass-stat: A technique for inserting or deleting particles in molecular dynamics simulations

The emergence of new applications of molecular dynamics (MD) simulation calls for the development of mass-statting procedures that insert or delete particles on-the-fly. In this paper we present a new mass-stat which we term FADE, because it gradually “fades-in” (inserts) or “fades-out” (deletes) molecules over a short relaxation period within a MD simulation. FADE applies a time-weighted relaxation to the intermolecular pair forces between the inserting/deleting molecule and any neighbouring molecules. The weighting function we propose in this paper is a piece-wise polynomial that can be described entirely by two parameters: the relaxation time scale and the order of the polynomial. FADE inherently conserves overall system momentum independent of the form of the weighting function. We demonstrate various simulations of insertions of atomic argon, polyatomic TIP4P water, polymer strands, and C{sub 60} Buckminsterfullerene molecules. We propose FADE parameters and a maximum density variation per insertion-instance that restricts spurious potential energy changes entering the system within desired tolerances. We also demonstrate in this paper that FADE compares very well to an existing insertion algorithm called USHER, in terms of accuracy, insertion rate (in dense fluids), and computational efficiency. The USHER algorithm is applicable to monatomic and water molecules only, but we demonstrate that FADE can be generally applied to various forms and sizes of molecules, such as polymeric molecules of long aspect ratio, and spherical carbon fullerenes with hollow interiors.

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

Molecular dynamics simulation of a binary mixture near the lower critical point

2,6-lutidine molecules mix with water at high and low temperatures but in a wide intermediate temperature range a 2,6-lutidine/water mixture exhibits a miscibility gap. We constructed and validated an atomistic model for 2,6-lutidine and performed molecular dynamics simulations of 2,6-lutidine/water mixture at different temperatures. We determined the part of demixing curve with the lower critical point. The lower critical point extracted from our data is located close to the experimental one. The estimates for critical exponents obtained from our simulations are in a good agreement with the values corresponding to the 3D Ising universality class.

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

Evaporative cooling of microscopic water droplets in vacuo : Molecular dynamics simulations and kinetic gas theory

In the present study, we investigate the process of evaporative cooling of nanometer-sized droplets in vacuum using molecular dynamics simulations with the TIP4P/2005 water model. The results are compared to the temperature evolution calculated from the Knudsen theory of evaporation which is derived from kinetic gas theory. The calculated and simulation results are found to be in very good agreement for an evaporation coefficient equal to unity. Lastly, our results are of interest to experiments utilizing droplet dispensers as well as to cloud micro-physics.

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

Rigidity and soft percolation in the glass transition of an atomistic model of ionic liquid, 1-ethyl-3-methyl imidazolium nitrate, from molecular dynamics simulations—Existence of infinite overlapping networks in a fragile ionic liquid

The typical ionic liquid, 1-ethyl-3-methyl imidazolium nitrate (EMIM-NO{sub 3}), was examined by molecular dynamics simulations of an all-atomistic model to show the characteristics of networks of cages and/or bonds in the course of vitrification of this fragile glass-former. The system shows changes of dynamics at two characteristic temperatures, T{sub B} (or T{sub c}) and the glass transition temperature T{sub g}, found in other fragile glass forming liquids [K. L. Ngai and J. Habasaki, J. Chem. Phys. 141, 114502 (2014)]. On decreasing temperature, the number of neighboring cation-anion pairs, N{sub B}, within the first minimum of the pair correlation function, g(r){sub min}, increases. On crossing T{sub B} (>T{sub g}), the system volume and diffusion coefficient both show changes in temperature dependence, and as usual at T{sub g}. The glass transition temperature, T{sub g}, is characterized by the saturation of the total number of “bonds,” N{sub B} and the corresponding decrease in degree of freedom, F = [(3N − 6) − N{sub B}], of the system consisting of N particles. Similar behavior holds for the other ion-ion pairs. Therefore, as an alternative, the dynamics of glass transition can be interpreted conceptually by rigidity percolation. Before saturation occurring at T{sub g}, the number of bonds shows a remarkable change at around T{sub B}. This temperature is associated with the disappearance of the loosely packed coordination polyhedra of anions around cation (or vice versa), related to the loss of geometrical freedom of the polyhedra, f{sub g}, of each coordination polyhedron, which can be defined by f{sub g} = [(3N{sub V} − 6) − N{sub b}]. Here, 3N{sub v} is the degree of freedom of N{sub V} vertices of the polyhedron, and N{sub b} is number of fictive bonds. The packing of polyhedra is characterized by the soft percolation of cages, which allows further changes with decreasing temperature. The power spectrum of displacement of the central ion in the cage is found to be correlated with the fluctuation of N{sub b} of cation-cation (or anion-anion) pairs in the polyhedron, although the effect from the coordination shells beyond the neighboring ions is not negligible.

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

Chemical dynamics simulations of the monohydrated OH{sup −}(H{sub 2}O) + CH{sub 3}I reaction. Atomic-level mechanisms and comparison with experiment

Direct dynamics simulations, with B97-1/ECP/d theory, were performed to study the role of microsolvation for the OH{sup −}(H{sub 2}O) + CH{sub 3}I reaction. The S{sub N}2 reaction dominates at all reactant collision energies, but at higher collision energies proton transfer to form CH{sub 2}I{sup −}, and to a lesser extent CH{sub 2}I{sup −} (H{sub 2}O), becomes important. The S{sub N}2 reaction occurs by direct rebound and stripping mechanisms, and 28 different indirect atomistic mechanisms, with the latter dominating. Important components of the indirect mechanisms are the roundabout and formation of S{sub N}2 and proton transfer pre-reaction complexes and intermediates, including [CH{sub 3}--I--OH]{sup −}. In contrast, for the unsolvated OH{sup −} + CH{sub 3}I S{sub N}2 reaction, there are only seven indirect atomistic mechanisms and the direct mechanisms dominate. Overall, the simulation results for the OH{sup −}(H{sub 2}O) + CH{sub 3}I S{sub N}2 reaction are in good agreement with experiment with respect to reaction rate constant, product branching ratio, etc. Differences between simulation and experiment are present for the S{sub N}2 velocity scattering angle at high collision energies and the proton transfer probability at low collision energies. Equilibrium solvation by the H{sub 2}O molecule is unimportant. The S{sub N}2 reaction is dominated by events in which H{sub 2}O leaves the reactive system as CH{sub 3}OH is formed or before CH{sub 3}OH formation. Formation of solvated products is unimportant and participation of the (H{sub 2}O)CH{sub 3}OH---I{sup −} post-reaction complex for the S{sub N}2 reaction is negligible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗