Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical studies”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Molecular dynamics studies of the size and temperature dependence of the kinetics of freezing of Fe nanoparticles

Molecular dynamics (MD) computer simulations have been carried out and a novel modified technique of Voronoi polyhedra has been performed to identify solid-like particles in a molten nanoparticle. This technique works quite well in analyzing the effects of particle size on nucleation rates of iron nanoparticles in the temperature range of 750–1160 K. Nanoparticles with 1436 and 2133 Fe atoms have been examined and the results are compared with those obtained earlier with Fe{sub 331} nanoparticles. Nucleation rates for freezing obtained from MD simulations for Fe{sub 2133} vary from 8.8×10{sup 34} m{sup 3}/s to 4.1×10{sup 35} m{sup 3}/s at over a temperature range from 1160 K to 900 K, Rates for. Fe{sub 1436} and Fe{sub 331} are somewhat higher. Nucleation rates increase as supercooling deepens until the viscosity of the liquid increases sharply enough to slow down the rate. Bt applying classical nucleation theory, the interfacial free energy between solid and liquid cab be estimated From this and other thermodynamic information can be derived a theoretical expression for the size-dependence of the heat of fusion of nanoparticles. Results agreed quite well with those observed in our MD observations. An earlier expression in the literature for this size-dependence was shown to be incorrect. The size dependence of melting point is discussed. - Graphical abstract: Critical nuclei of crystallization. Display Omitted - Highlights: • Solid state material synthesis. • Material structure. • Experimental study of nucleation in condensed materials. • Computation study of nucleation in condensed materials.

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

All-Atom Molecular Dynamics Study of Water–Dodecane Interface in the Presence of Octanol

In this work, classical atomistic molecular dynamics (MD) simulations were performed to investigate the structural properties of octanol at water/dodecane interface. Multiple combinations of simulation conditions and force fields were employed in order to :obtain high fidelity interfacial structure. The analysis of noncovalent interactions shows that the convergence of interfacial octanol distribution requires a simulation time scale larger than 100 ns. Finite size effect diminishes with increasing dimension of simulation box both laterally and vertically. In the largest simulated systems with vanishing finite size effect, the additive CHARMM force field shows similar interfacial aggregation behavior with the DRUDE polarizable model. Self-assembly of octanol molecules at water/oil interface is observed, expected to reduce;the free energy barrier between water phase and oil phase. This work provides an MD simulation benchmark for structural study of complex liquid/liquid interface with cosurfactant reducing interfacial tension between two immiscible phases.

Qiao, Baofu↗
DOE OSTI2019

Ligand-Dependent Sodium Ion Dynamics within the A 2A Adenosine Receptor: A Molecular Dynamics Study

Sodium ions have long been known to reduce the binding of agonists in many class-A GPCRs while having little effect on antagonist binding. In this study, using long-time scale classical all-atom molecular dynamics simulations, we explore, in atomic detail, the motion of sodium ions within the ligand-binding pocket of the A 2A adenosine receptor (A2A-AR) both in the presence and absence of ligands and in the active and inactive state. We identify novel secondary ion binding sites within the pocket and find that the types of ion motions within the pocket are highly dependent on the presence and type of ligand within the pocket. Finally, our results provide a first step toward developing a molecular understanding of the impact of sodium ions on class-A GPCRs.

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

Molecular dynamics studies of helium bubble effects on grain boundary fracture vulnerabilities in an Fe 70 Ni 11 Cr 19 –1%H austenitic stainless steel

Comprehensive molecular dynamics tensile test simulations have been performed to study the delamination processes of seven different grain boundaries / cleavage planes (Σ1{111}, Σ3{111}, Σ5{100}, Σ7{111}, Σ9{411}, Σ11{311}, and R{100}/{411}) containing a helium bubble. Combinations of a variety of conditions are explored including different strain rates, system dimensions, bubble density, bubble radius, bubble pressure, and temperature. We found that in general, grain boundaries absorb less energies with decreasing strain rate but increasing bubble areal density, bubble pressure, bubble radius, and temperature. The propensity of grain boundary delamination is sensitive to grain boundary type: The random grain boundary R{100}/{411} is one of the most brittle boundaries whereas the Σ1{111} cleavage plane and the Σ3{111} twin boundary are two of the toughest boundaries. The sorted list of grain boundary fracture vulnerability obtained from our dynamic tensile test simulations differs from the one obtained from our decohesion energy calculations, confirming the important role of plastic deformation during fracture. Finally, detailed mechanistic analyses are performed to interpret the simulated results.

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

Six-dimensional quantum dynamics study for the dissociative adsorption of HCl on Au(111) surface

The six-dimensional quantum dynamics calculations for the dissociative chemisorption of HCl on Au(111) are carried out using the time-dependent wave-packet approach, based on an accurate PES which was recently developed by neural network fitting to density functional theory energy points. The influence of vibrational excitation and rotational orientation of HCl on the reactivity is investigated by calculating the exact six-dimensional dissociation probabilities, as well as the four-dimensional fixed-site dissociation probabilities. The vibrational excitation of HCl enhances the reactivity and the helicopter orientation yields higher dissociation probability than the cartwheel orientation. A new interesting site-averaged effect is found for the title molecule-surface system that one can essentially reproduce the six-dimensional dissociation probability by averaging the four-dimensional dissociation probabilities over 25 fixed sites.

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

Interfacial thermal conductance between multi-layer graphene sheets and solid/liquid octadecane: A molecular dynamics study

Mixtures of paraffin and carbon nanofillers have promising potential for thermal storage, as paraffin (the matrix) possesses high latent heat and the nanofillers compensate for the low thermal conductivity (TC) of paraffin. Understanding thermal transport in these materials is essential for practical applications, as weak thermal transport hinders fast charge/discharge of thermal energy. Here, we use non-equilibrium molecular dynamics (NEMD) simulations to study the interfacial thermal conductance (ITC) between graphene sheets and octadecane (C18H38) matrix under the limiting conditions of the sheets being parallel or perpendicular to the direction of the imposed heat flux. The findings show that the systems containing thin graphene layers exhibit higher values of ITC. This study captures the asymptotic saturation of thermal conductance for the liquid phase of the perpendicular structure. Besides, given the greater number of structured layers of paraffin upon phase change, the ITC for the solid paraffin-graphene system is higher than the conductance of the liquid paraffin-graphene interface. We use the Pearson correlation coefficients of the vibrational power spectrum (VPS) of interfacing materials to explain the orders of magnitude variations of the observed ITC.

25 ENERGY STORAGE↗
DOE OSTI2014

Six-dimensional quantum dynamics study for the dissociative adsorption of DCl on Au(111) surface

We carried out six-dimensional quantum dynamics calculations for the dissociative adsorption of deuterium chloride (DCl) on Au(111) surface using the initial state-selected time-dependent wave packet approach. The four-dimensional dissociation probabilities are also obtained with the center of mass of DCl fixed at various sites. These calculations were all performed based on an accurate potential energy surface recently constructed by neural network fitting to density function theory energy points. The origin of the extremely small dissociation probability for DCl/HCl (v = 0, j = 0) fixed at the top site compared to other fixed sites is elucidated in this study. The influence of vibrational excitation and rotational orientation of DCl on the reactivity was investigated by calculating six-dimensional dissociation probabilities. The vibrational excitation of DCl enhances the reactivity substantially and the helicopter orientation yields higher dissociation probability than the cartwheel orientation. The site-averaged dissociation probability over 25 fixed sites obtained from four-dimensional quantum dynamics calculations can accurately reproduce the six-dimensional dissociation probability.

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

Ab Initio Molecular Dynamics Study of Pd Nucleation on a γ-Al 2 O 3 Surface

We perform ab initio molecular dynamics (AIMD) free energy simulations to investigate the stability and dynamics of dispersed Pd atoms and subnanometer clusters on γ-Al 2 O 3 . We consider nucleation for varying Pd loadings on dry γ-Al 2 O 3 (100), dry γ-Al 2 O 3 (110) and wet γ-Al 2 O 3 (110) and examine temperature, entropy and water coverage effects on the stability and dynamics of Pd clusters. At low Pd loading (1.07 Pd atoms/nm 2 on the (100) and 0.75 Pd/nm 2 on the (110)), association of the Pd atoms is thermodynamically unfavorable on the dry surfaces and almost thermoneutral on the wet (110) owing to water blocking of the AlIII defect sites. In agreement with past studies, the Pd atoms are quite mobile but the diffusion free energy barriers do not increase with the water coverage, as previously suggested, but rather remain similar in magnitude owing to translational entropy at finite temperatures. At these Pd loadings, oxidized Pd(II) atoms on the wet (110) surface are dispersed and immobile. At higher Pd loading (1.12 and 1.5 Pd atoms/nm 2 ) on the (110) surface, Pd 3 and Pd 4 clusters are thermodynamically very stable and form fast. At the 1.12 Pd atoms/nm 2 loading, the formation of the Pd 3 cluster is two orders of magnitude slower when the surface is hydroxylated as the system becomes somewhat trapped before the dinuclear cluster Pd 2 collides with a Pd atom.

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

First-Principles Molecular Dynamics Study of a Deep Eutectic Solvent: Choline Chloride/Urea and Its Mixture with Water

First-principles molecular dynamics simulations in the canonical ensemble at temperatures of 333 and 363 K and at the corresponding experimental densities are carried out to investigate the behavior of the 1:2 choline chloride/urea (reline) deep eutectic solvent and its equimolar mixture with water. Analysis of atom–atom radial and spatial distribution functions and of the H-bond network reveals the microheterogeneous structure of these complex liquid mixtures. In neat reline, the structure is governed by strong H-bonds of the trans- and cis-H atoms of urea to the chloride ion. In hydrous reline, water competes for the anions, and the H atoms of urea have similar propensities to bond to the chloride ions and the O atoms of urea and water. Finally, the vibrational spectra exhibit relatively broad peaks reflecting the heterogeneity of the environment. Although the 100 ps trajectories allow only for a qualitative assessment of transport properties, the simulations indicate that water is more mobile than the other species and its addition also fosters faster motion of urea.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Dynamical Study of Adsorbate-Induced Restructuring Kinetics in Bimetallic Catalysts Using the PdAu(111) Model System

The dynamic restructuring of bimetallic catalysts plays a crucial role in their catalytic activity and selectivity. In particular, catalyst pretreatment with species such as carbon monoxide and oxygen has been shown to be an effective strategy for tuning the surface composition and morphology. Mechanistic and kinetic understanding of such restructuring are fundamental to the chemistry and engineering of surface active sites but have remained challenging due to the large structural, chemical, and temporal degrees of freedom. Here, we combine time-resolved temperature-programmed infrared reflection absorption spectroscopy, ab initio thermodynamics, and machine-learning molecular dynamics to uncover previously unidentified timescale and kinetic parameters of in situ restructuring in Pd/Au(111), a highly relevant model system for dilute Pd-in-Au nanoparticle catalysts. The key innovation lies in utilizing CO not only as a chemically sensitive probe of surface Pd, but also as an agent that induces restructuring of the surface. Upon annealing in vacuum, as-deposited Pd islands became encapsulated by Au and partially dissolved into the subsurface, leaving behind isolated Pd monomers on the surface. Subsequent exposure to 0.1 mbar CO enabled Pd monomers to repopulate the surface up to 373 K, above which complete Pd dissolution occurred by 473 K, with apparent activation energies of 0.14 and 0.48 eV, respectively. Furthermore, these restructuring processes occurred over the span of ~1000 s at a given temperature. Such a minute-timescale dynamics not only elucidates the fluxional nature of alloy catalysts but also presents an opportunity to fine-tune the surface at moderate temperature and pressure conditions.

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

Structural dynamics studies of rotating bladed-disk assemblies coupled with flexible shaft motions

In order to analyze the dynamic behavior of the first stage compressor/fan of the 'E3' turbofan engine, a classical structural dynamics approach is employed to couple the motions of a flexible bladed disk to a rotating flexible shaft. The analysis accounts for flexible disk displacements which are transverse to the plane of rotation, and radial as well as tangential, and also accounts for rigid disk translations along, and rotations about, axes normal to the undeformed shaft axes. In the case of a wide range of E3 engine shaft flexibilities and speeds, some of the one-diametral node frequencies are shown to be affected by shaft degrees of freedom whose stiffness values are in general range of design practice. Coriolis forces are also found to significantly affect natural frequencies where strong coupling between certain modes is present.

Loewy, R. G.↗
NASA NTRS1995

Molecular dynamics studies of simple membrane-water interfaces: Structure and functions in the beginnings of cellular life

Molecular dynamics computer simulations of the structure and functions of a simple membrane are performed in order to examine whether membranes provide an environment capable of promoting protobiological evolution. Our model membrane is composed of glycerol 1-monooleate. It is found that the bilayer surface fluctuates in time and space, occasionally creating thinning defects in the membrane. These defects are essential for passive transport of simple ions across membranes because they reduce the Born barrier to this process by approximately 40%. Negative ions are transferred across the bilayer more readily than positive ions due to favorable interactions with the electric field at the membrane-water interface. Passive transport of neutral molecules is, in general, more complex than predicted by the solubility-diffusion model. In particular, molecules which exhibit sufficient hydrophilicity and lipophilicity concentrate near membrane surfaces and experience 'interfacial resistance' to transport. The membrane-water interface forms an environment suitable for heterogeneous catalysis. Several possible mechanisms leading to an increase of reaction rates at the interface are discussed. We conclude that vesicles have many properties that make them very good candidates for earliest protocells. Some potentially fruitful directions of experimental and theoretical research on this subject are proposed.

Pohorille, Andrew↗
NASA NTRS1999

Early Events in the Folding of an Amphipathic Peptide A Multi- Nanosecond Molecular Dynamics Study

Folding of the capped LQQLLQQLLQL peptide is investigated at the water-hexane interface by molecular dynamics simulations over 161.5 nanoseconds. Initially placed in the aqueous phase as a beta-strand, the peptide rapidly adsorbs to the interface, where it adopts an amphipathic conformation. The marginal presence of non-amphipathic structures throughout the complete trajectory indicate- that the corresponding conformations are strongly disfavored at the interface. It is further suggestive that folding in an interfacial environment proceeds through a pathway of successive amphipathic intermediates. The energetic and entropic penalties involved in the conformational changes along this pathway markedly increase the folding time-scales of LQQLLQQLLQL, explaining why the alpha-helix, the hypothesized lowest free energy structure for a sequence with a hydrophobic periodicity of 3.6, has not been reached yet. The formation of a type I beta-turn at the end of the simulation confirms the importance of such motifs as initiation sites allowing the peptide to coalesce towards a secondary structure.

Chipot, Christophe↗
NASA NTRS1999

Early events in the folding of an amphipathic peptide: A multinanosecond molecular dynamics study

Folding of the capped LQQLLQQLLQL peptide is investigated at the water-hexane interface by molecular dynamics simulations for 161.5 ns. Initially placed in the aqueous phase as a beta-strand, the peptide rapidly adsorbs to the interface, where it adopts an amphipathic conformation. The marginal presence of nonamphipathic structures throughout the complete trajectory indicates that the corresponding conformations are strongly disfavored at the interface. It is further suggestive that folding in an interfacial environment proceeds through a pathway of successive amphipathic intermediates. The energetic and entropic penalties involved in the conformational changes along this pathway markedly increase the folding time scales of LQQLLQQLLQL, explaining why the alpha-helix, the hypothesized lowest free energy structure for a sequence with a hydrophobic periodicity of 3.6, has not been reached yet. The formation of a type I beta-turn at the end of the simulation confirms the importance of such motifs as initiation sites allowing the peptide to coalesce towards a secondary structure. Proteins 1999;36:383-399. Copyright 1999 Wiley-Liss, Inc.

NASA Center ARC↗
DOE OSTI2023

Molecular dynamics study of carbon dioxide and nitrogen selectivity through poly[bis((methoxyethoxy)ethoxy)phosphazene] (MEEP) membrane

Here, a molecular dynamics simulation model was developed to comprehensively understand carbon dioxide over nitrogen (CO 2 /N 2 ) selectivity through polyphosphazene-based membrane comprising of poly[bis((methoxyethoxy)ethoxy)phosphazene] (MEEP) selective layer at the molecular level. The effects of temperature, pressure, and initial feed gas composition on the CO 2 transport on a polymer membrane were studied. The computed free energy and density profile of the permeating gas mixture exhibited that CO 2 molecules express higher interactions with the membrane than N 2 molecules, resulting in higher CO 2 /N 2 selectivity. Statistical analysis of gas molecules (CO 2 , water (H 2 O), and N 2 ) transportation suggested that hydro- and CO 2 -philic functional groups in the membrane significantly impact CO 2 permeability and CO 2 /N 2 selectivity. This study suggested that Lewis acid–base and hydrogen bonding combinations contribute to CO 2 permeation and CO 2 /N 2 selectivity. An equal CO 2 /N 2 selectivity was observed with and without water vapor in the feed gas suggesting that water does not hinder CO 2 transport through the membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Molecular dynamics study of the effect of alkyl chain length on melting points of [CnMIM][PF6] ionic liquids

Based on molecular dynamics simulations, the melting points T-m of a series of 1-alkyl-3-methylimidazolium hexafluorophosphate ionic liquids [CnMIM][PF6] with n = 2, 4, 10, 12, and 14 were studied using the free energy-based pseudosupercritical path (PSCP) method. The experimental trend that the Tm decreases with increasing alkyl chain length for ILs with short alkyl chains and increases for the ones with long alkyl chains was correctly captured. Further analysis revealed that the different trends are the results of the balance between fusion enthalpy and fusion entropy. For the ILs with short alkyl chains (ethyl and butyl groups), fusion entropy plays the dominant role so that [C4MIM][PF6], which has a larger fusion entropy due to its higher liquid phase entropy has the lower melting temperature. As for the ILs with long alkyl chains, due to the enhanced van der Waals interactions brought about by the long non-polar alkyl chains, enthalpy becomes the deciding factor and the melting points increase when the alkyl chain goes from C10 to C14. While the melting points for [C2MIM][PF6] and [C4MIM][PF6] were quantitatively predicted and the trends for the long chain ILs were captured correctly, the absolute melting points for [C10MIM][PF6], [C12MIM][PF6] and [C14MIM][PF6] were systematically overestimated in the simulations. Three possible reasons for the overestimation were studied but all ruled out. Further simulation or experimental studies are needed to explain the difference.

Zhang, Y↗
DOE OSTI2013

Femtosecond two-photon photoassociation of hot magnesium atoms: A quantum dynamical study using thermal random phase wavefunctions

Two-photon photoassociation of hot magnesium atoms by femtosecond laser pulses, creating electronically excited magnesium dimer molecules, is studied from first principles, combining ab initio quantum chemistry and molecular quantum dynamics. This theoretical framework allows for rationalizing the generation of molecular rovibrational coherence from thermally hot atoms [L. Rybak, S. Amaran, L. Levin, M. Tomza, R. Moszynski, R. Kosloff, C. P. Koch, and Z. Amitay, Phys. Rev. Lett. 107, 273001 (2011)]. Random phase thermal wavefunctions are employed to model the thermal ensemble of hot colliding atoms. Comparing two different choices of basis functions, random phase wavefunctions built from eigenstates are found to have the fastest convergence for the photoassociation yield. The interaction of the colliding atoms with a femtosecond laser pulse is modeled non-perturbatively to account for strong-field effects.

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

Hybrid Modeling of Three-Phase Grid-Supporting Inverters for Dynamic Studies

Grid technologies connected by power electronic converter (PEC) interfaces continually implement grid support functions mandated by grid codes and standards. The transition to converter-based generation demands precise PEC models to assess system dynamics, which have been previously overlooked in conventional power systems. This study proposes a hybrid method for analyzing grid-connected three-phase PEC dynamics with the IEEE standard 1547-2018 Volt-VAr mode that combines physics and data-driven techniques. The physics model reflects the PEC’s internal behavior, whereas the data-driven modeling technique evaluates the grid-supporting capabilities of the smart PEC. The system identification approach is used to generate dynamic PEC models based on changing grid voltage and measured current injected into the grid by the PEC. In the Volt-VAr support mode, a detailed topological model including switches is utilized to compare the goodness-of-fit of the extracted hybrid dynamic model. The results demonstrate that the hybrid PEC model in the Volt-VAr mode accurately matches the dynamics with the topological model.

Subedi, Sunil↗