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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

New Programs Utilizing Light Scattering and Flow Imaging Techniques for Macromolecular Crystal Growth and Fluid Dynamics Studies

Dr. Phil Segre, a physicist by training, is a recent addition to the Biotech group, SD46, having joined NASA in August of 2000. Over the past two years he has been developing a laboratory for the study of macromolecular and protein crystal growth. The main apparatus for this work is a Dynamic Light Scattering apparatus, DLS, which is capable of making highly precise measurements of size distributions of both protein solutions and protein crystals. With Drs. Chernov and Thomas (USRA), he has begun a collaboration studying the affects of protein impurities on protein crystal growth and subsequent crystal quality. One of the hypotheses behind the differences between Earth and space grown protein crystals is that the absorption of harmful impurities is reduced in space due to the absence of convective flows. Using DLS measurements we are examining crystal growth with varying amounts of impurities and testing whether there is a strong physical basis behind this hypothesis. With Dr. Joe Ng of UAH he has been collaborating on a project to examine the folding/unfolding dynamics of large RNA complexes. A detailed understanding of this process is necessary for the handling of RNA in biotech applications, and the DLS instrument gives details and results beyond that of other instruments. With Prof. Jim McClymer of the University of Maine (summer faculty visitor to NASA in 2001, 2002), we have been studying the crystallization process in model colloidal suspensions whose behavior in some cases can mimic that of much smaller protein solutions. An understanding of the self-assembly of colloids is the first step in the process of engineering novel materials for photonic and light switching applications. Finally, he has begun an investigation into the physics of particle sedimentation. In addition to the DLS instrument he also has an instrument (called PIV) that can measure flow fields of fluids. The applications are to the dynamics of protein crystal motions both on earth and in low-gravity.

Source record↗
DOE OSTI2014

Fragmentation of displacement cascades into subcascades: A molecular dynamics study

The fragmentation of displacement cascades into subcascades in copper and iron has been investigated through the molecular dynamics technique. A two-point density correlation function has been used to analyze the cascades as a function of the primary knock-on (PKA) energy. This approach is used as a tool for detecting subcascade formation. The fragmentation can already be identified at the end of the ballistic phase. Its resulting evolution in the peak damage state discriminates between unconnected and connected subcascades. The damage zone at the end of the ballistic phase is the precursor of the extended regions that contain the surviving defects. A fractal analysis of the cascade exhibits a dependence on both the stage of the cascade development and the PKA energy. This type of analysis enables the minimum and maximum displacement spike energies together with the subcascade formation threshold energy to be determined. (C) 2014 Elsevier B.V. All rights reserved.

Antoshchenkova, E.↗
DOE OSTI2021

Effective Fragment Potential-Based Molecular Dynamics Studies of Diffusion in Acetone and Hexane

To facilitate more reliable descriptions of transport properties in liquids, molecular dynamics (MD) simulations are performed based on the effective fragment potential (EFP) method derived from first-principles quantum mechanics (in contrast to MD based upon empirically fitted potentials). The EFP method describes molecular interactions in terms of Coulomb, polarization/induction, dispersion, exchange-repulsion, and charge-transfer interactions. The EFP MD simulations described in this paper, performed on hexane and acetone, are able to track the mean-square displacement of molecules for sufficient time to reliably extract translational diffusion coefficients. Finally, the results reported here are in reasonable agreement with experiment.

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

Role of electronic thermal transport in amorphous metal recrystallization: A molecular dynamics study

Recrystallization of glasses is important in a wide range of applications including electronics and reactive materials. Molecular dynamics (MD) has been used to provide an atomic picture of this process, but prior work has neglected the thermal transport role of electrons, the dominant thermal carrier in metallic systems. We characterize the role of electronic thermal conductivity on the velocity of recrystallization in Ni using MD coupled to a continuum description of electronic thermal transport via a two-temperature model. Our simulations show that for strong enough coupling between electrons and ions, the increased thermal conductivity removes the heat from the exothermic recrystallization process more efficiently, leading to a lower effective temperature at the recrystallization front and, consequently, lower propagation velocity. We characterize how electron-phonon coupling strength and system size affect front propagation velocity. Interestingly, we find that initial recrystallization velocity increases with decreasing system size due to higher overall temperatures. Overall, we show that a more accurate description of thermal transport due to the incorporation of electrons results in better agreement with experiments.

McClure, Zachary D. (ORCID:0000000183300217)↗
DOE OSTI2020

Molecular dynamics study on interface formation and bond strength of impact-welded Mg-steel joints

It was recently demonstrated that the vaporizing foil actuator welding (VFAW) method can directly join immiscible magnesium and steel alloys without coating or a third chemical element based intermetallic compound layer. The VFAW Mg/steel joint exhibits a mixed interface layer of up to 200μm thickness consisting of Mg matrix and Fe particles. Computer simulations have suggested the formation of the interlayer is from the high-velocity frictional shearing between the Mg/steel substrates during the oblique impact in the VFAW process. This paper investigates the formation of Mg-Fe interlayer under VFAW condition with different shearing velocities using molecular dynamics (MD) model, and studies the bonding strength under different scenarios. Finally, the results elucidate the critical role of shearing velocity and surface roughness in achieving Mg/Fe joint.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Contribution of the Hydroxyl Group on Interfacial Heat Conduction of Monohydric Alcohols: A Molecular Dynamics Study

Monohydric alcohols have been used as promising phase change materials (PCMs) for low-temperature latent heat storage. However, the heat storage/retrieval rates are limited due to the low thermal conductivity of such alcohols. In this work, nonequilibrium molecular dynamics (NEMD) simulations were performed to study the microscopic heat conduction in example monohydric alcohols, i.e., 1-dodecanol (C12H26O), 1-tetradecanol (C14H30O), and 1-hexadecanol (C16H34O). A simplified ideal crystal model was proposed to exploit the potential for improving the thermal conductivity of monohydric alcohols. The effect of ideal crystalline structures, especially the contribution of the hydroxyl group, on the microscopic heat conduction process was analyzed. The thermal conductivity of the ideal crystals of the various monohydric alcohols was predicted to be more than twice as compared to that of their respective solids. The major thermal resistance in the ideal crystals was found around the molecular interfaces, as a result of the excellent heat conduction performance along the linear molecular chains. The calculated vibrational density of states (VDOS) and interfacial heat transfer were then investigated. When the interfaces are surrounded by hydroxyl groups as walls, strong hydrogen bond (HB) interactions were observed. The interfacial heat transfer coefficient of the ideal crystalline structures of 1-tetradecanol was found to reach up to ~735.6 MW/m2 W. It was elucidated that the high interfacial heat transfer rate is clearly related to the stronger intermolecular interactions.

Engineering↗
DOE OSTI2024

Molecular Dynamics Study of the Effect of Grafting Density on Ion Diffusivity in a MARTINI Coarse-Grained Strong Polyelectrolyte Brush

Because surface-grafted polyelectrolyte brushes (PEBs) are responsive to external stimuli, such as electric fields and ionic strength, PEBs are attractive for applications ranging from drug delivery to separation technologies. Essential to PEB utilization is understanding how critical parameters like grafting density (σ) impact the PEB structure and the dynamics of the PEB and counterions. To study the effect of σ on PEB and the counterion structure and dynamics, we fine-tune a coarse-grained model that retains the chemical specificity of a strong polyelectrolyte, poly[(2-(methacryloyloxy)ethyl) trimethylammonium chloride] (PMETAC), using the MARTINI force field. Using “salt-free” conditions where the counterion concentration balances the charge on the brush, we build coarse-grained (CG) molecular dynamics simulations for MARTINI PMETAC brushes (N = 150 monomers; M W = 31.2 kg/mol) at experimentally relevant values of σ = 0.05, 0.10, 0.20, and 0.40 chains/nm 2 . Using 5 μs simulations, we investigate the effects of grafting density on the PEB structure, ion dissociation dynamics, polymer mobility, and counterion diffusivity. Results show that competition between electrostatic interactions, steric hindrance, and polymer mobility controls counterion diffusivity. Finally, the interplay of these factors leads to diffusivity that depends non-monotonically on σ, with counterion diffusivity peaking at an intermediate σ = 0.10 chains/nm 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2014

Origin of the reverse optical-contrast change of Ga-Sb phase-change materials—An ab initio molecular-dynamics study

A large number of phase-change materials (PCMs) have been developed experimentally; however, only Ge{sub 2}Sb{sub 2}Te{sub 5}-based PCMs have been significantly explored using ab initio molecular-dynamics (AIMD) simulations. We present an AIMD study of the full melt/quench/anneal PC cycle for Ga-Sb materials, namely, the stoichiometric composition, GaSb, and the near-eutectic alloy, Ga{sub 16}Sb{sub 84}. The calculated electronic densities of states and optical reflectivities are compared between the amorphous and crystalline phases for both compositions, and it is shown that the contrasting opto-electronic properties of each crystalline material can be attributed to different structural transformations of Ga and Sb on crystallization from the amorphous state.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Ab initio molecular dynamics study of the properties of cerium in liquid sodium at 1000 K temperature

For liquid-sodium-cooled fast nuclear reactor systems, it is crucial to understand the behavior of lanthanides and other potential fission products in liquid sodium or other liquid metal solutions such as liquid cesium-sodium. In this study, we focus on lanthanide behavior in liquid sodium. Using ab initio molecular dynamics, we found that the solubility of cerium in liquid sodium at 1000 K was less than 0.78 at. %, and the diffusion coefficient of cerium in liquid sodium was calculated to be 5.57 × 10{sup −9} m{sup 2}/s. Furthermore, it was found that cerium in small amounts may significantly alter the heat capacity of the liquid sodium system. Our results are consistent with the experimental results for similar materials under similar conditions.

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

Hydration of divalent lanthanides, Sm 2+ and Eu 2+ : A molecular dynamics study with polarizable AMOEBA force field

The chemistry of divalent lanthanides, Ln 2+ , is a growing sub-field of heavy element chemistry owing to new synthetic approaches. However, some theoretical aspects of these unusual cations are currently underdeveloped, especially as they relate to their dynamic properties in solution. In this study, we address the hydration of two of the classical Ln 2+ cations, Sm 2+ and Eu 2+ , using atomic multipole optimized energetic for biomolecular applications (AMOEBA) force fields. These cations have not been parameterized to date with AMOEBA, and few studies are available because of their instability with respect to oxidation in aqueous media. Coordination numbers (CN's) of 8.2 and 8.1 respectively for Sm 2+ and Eu 2+ , and 8.8 for both Sm 3+ and Eu 3+ have been obtained and are in good agreement with the few available AIMD and X-ray absorption fine structures studies. The decreased CN of Ln 2+ compared with Ln 3+ arises from progressive water exchange events that indicates the gradual stabilization of 8-coordinate structures with respect to 9-coordinate geometries. Moreover, the effects of the chloride counter anions on the coordination of Ln 2+ cations have been studied at different chloride concentrations in this work. Lastly, water exchange times of Ln 2+ cations have been calculated to provide a comprehensive understanding of the behavior of Eu 2+ and Sm 2+ in aqueous chloride media.

AMOEBA force field parameterization↗
DOE OSTI2023

Phase transition in bilayer MoS 2 under tensile loading: a molecular dynamics study

Abstract Molybdenum disulfide (MoS 2 ), especially single-layer MoS 2 , has been experimentally and computationally discovered to exist in several different polymorphs exhibiting various electronic and mechanical properties. The morphology of MoS 2 can be tuned through strain engineering. Molecular dynamics simulations are conducted to systematically study the phase transition of single-layer MoS 2 and bilayer MoS 2 under the uniaxial tensile condition at room temperature. The roles of edge and S-line vacancy are investigated. Phase transitions are always triggered near the edge and vacancy sites. The initiation of the metastable T″ phase can release the tensile stress in the lattice, followed by I4/mmm phase initiation, regardless of the edge conditions. The growth of the I4/mmm phase can cause the local buckling of the MoS 2 plane. With a tilted S-line vacancy, I4/mmm phase is first initiated to reduce the local shear stress accumulated near the vacancy line. Overall, the phase transition mechanism of single layer and bilayer MoS 2 under the uniaxial tensile loading is provided, which guides the future strain engineering of MoS 2 in nanoelectronics applications.

Materials Science↗
DOE OSTI2025

Molecular-dynamics study of diffusional creep in uranium mononitride

Uranium mononitride (UN) is a promising advanced nuclear fuel due to its high thermal conductivity and high fissile density. Yet, many aspects of its mechanical behavior and microstructural features are currently unknown. In this paper, molecular dynamics (MD) simulations are used to study UN's diffusional creep. Nanometer-sized polycrystals are used to simulate diffusional creep and to calculate an effective GB width. It is found that Nabarro-Herring creep is not dominant in the temperature range of 1700–2000 K and that the dominant diffusional creep mechanism is Coble creep with an activation energy of 2.28 ± 0.09 eV. A method is proposed to calculate the diffusional GB width and its temperature dependence in polycrystals. The effective GB width of UN is calculated as 2.69 ± 0.08 nm. This value fits very well with the prefactor of the phenomenological Coble creep formula. It is demonstrated that the most comprehensive thermal creep model for UN can be represented as the combination of our Coble creep model and the dislocation creep model proposed by Hayes et al.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2013

A non-adiabatic dynamics study of octatetraene: The radiationless conversion from S{sub 2} to S{sub 1}

Simulation of the excited state dynamics of all-trans-1,3,5,7-octatetraene has been performed to investigate the ultrafast radiationless S{sub 2} → S{sub 1} internal conversion process. Multireference configuration interaction with single excitation method has been employed to optimize the equilibrium structure of the excited states, as well as the S{sub 2}/S{sub 1} conical intersection, and to investigate the non-adiabatic molecular dynamics of the S{sub 2}/S{sub 1} state transition. At the conical intersection, the molecule is found to be distorted from the original planar trans structure to a nearly perpendicular conformation around C{sub 3}−C{sub 4} bond, with the torsion angle being about 107°. Such structural change can result in mutual approaching of states S{sub 2} and S{sub 1} in energy, and drastically increase the nonadiabatic coupling between the two states by destroying the inter-state symmetry prohibition in the electronic wavefunctions. Surface-hopping molecular dynamics simulations are performed to describe the non-adiabatic process. Upon the Franck-Condon excitation to the S{sub 2} state, the molecule quickly twists its C{sub 3}−C{sub 4} bond and approaches the conical intersection region, where it can undergo efficient internal conversion to S{sub 1}. The decay time constant (τ) of S{sub 2} state is estimated to be around 251 fs by fitting the occupation number of average fraction of trajectories using an exponential damping function. This value is reasonably consistent with previous experimental measurements of around 300–400 fs.

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

Dynamic Studies of Photo- and Electron-Induced Reactions on Nanostructured Surfaces

Our current research program examines the photon- and electron-initiated reaction mechanisms, half-collision dynamics, and other nonequilibrium-excited dynamics effects, occurring with excitation of adsorbates on well-characterized metal-oxide of metalchalcogenide surfaces, both in single-crystal and in nanocrystalline forms. The program has made extensive use of STM methods, collaborations with the Nanoprobing SPLEEM system at Brookhaven National Lab, and XPS methods. As our current program terminates, it brings to an end a sustained and major effort to understand reactions and electronic structure on materials of interest for photochemistry on precision surfaces. The last two close-out years of our program allowed us to complete our STM surface studies, interact with scientists at BNL on programs useful for two dimensional (2D) materials, and initiate a new effort in ultrafast two-photon photoemission spectroscopy. Our work has then shown how nanoscale stressed surface regions on sold crystal surfaces (typically TiO 2 ) lead to enhanced reactivity of adsorbate molecules – a result of much importance for nanocatalysis. In addition we have developed several new surface and near-surface probing methods – these, also, are of general use to the community. Finally due to the pervasive interest in DOE labs in these materials, we have exploited our earlier synchrotron-based understanding of the electronic structure in (2D) materials to begin a study of surface chemical reactions in these same materials and their corresponding molecular dynamics. These efforts are described briefly below and in the references or DOE publications presented in the publications section of this report. We are pleased to have been able to contribute so many publications in the last two years of this many year DOE program. We also wish to thank DOE for the sustained funding that has made this program such a success.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2014

A molecular dynamics study of intramolecular proton transfer reaction of malonaldehyde in solutions based upon mixed quantum-classical approximation. I. Proton transfer reaction in water

In order to investigate proton transfer reaction in solution, mixed quantum-classical molecular dynamics calculations have been carried out based on our previously proposed quantum equation of motion for the reacting system [A. Yamada and S. Okazaki, J. Chem. Phys. 128, 044507 (2008)]. Surface hopping method was applied to describe forces acting on the solvent classical degrees of freedom. In a series of our studies, quantum and solvent effects on the reaction dynamics in solutions have been analysed in detail. Here, we report our mixed quantum-classical molecular dynamics calculations for intramolecular proton transfer of malonaldehyde in water. Thermally activated proton transfer process, i.e., vibrational excitation in the reactant state followed by transition to the product state and vibrational relaxation in the product state, as well as tunneling reaction can be described by solving the equation of motion. Zero point energy is, of course, included, too. The quantum simulation in water has been compared with the fully classical one and the wave packet calculation in vacuum. The calculated quantum reaction rate in water was 0.70 ps{sup −1}, which is about 2.5 times faster than that in vacuum, 0.27 ps{sup −1}. This indicates that the solvent water accelerates the reaction. Further, the quantum calculation resulted in the reaction rate about 2 times faster than the fully classical calculation, which indicates that quantum effect enhances the reaction rate, too. Contribution from three reaction mechanisms, i.e., tunneling, thermal activation, and barrier vanishing reactions, is 33:46:21 in the mixed quantum-classical calculations. This clearly shows that the tunneling effect is important in the reaction.

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

A seven-degree-of-freedom, time-dependent quantum dynamics study on the energy efficiency in surmounting the central energy barrier of the OH + CH{sub 3} → O + CH{sub 4} reaction

A time-dependent, quantum reaction dynamics calculation with seven degrees of freedom was carried out to study the energy efficiency in surmounting the approximate center energy barrier of OH + CH{sub 3}. The calculation shows the OH vibration excitations greatly enhance the reactivity, whereas the vibrational excitations of CH{sub 3} and the rotational excitations hinder the reactivity. On the basis of equal amount of total energy, although this reaction has a slight early barrier, it is the OH vibrational energy that is the dominate force in promoting the reactivity, not the translational energy. The studies on both the forward O + CH{sub 4} and reverse OH + CH{sub 3} reactions demonstrate, for these central barrier reactions, a small change of the barrier location can significantly change the energy efficacy roles on the reactivity. The calculated rate constants agree with the experimental data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A discrete dislocation dynamics study of precipitate bypass mechanisms in nickel-based superalloys

Order strengthening in nickel-based superalloys is associated with the extra stress required for dislocations to bypass the $\gamma'$ precipitates distributed in the $\gamma$ matrix. Depending on the operating conditions and microstructure, a rich variety of bypass mechanism has been identified, with various shearing and Orowan looping processes gradually giving way to climb bypass as the operating conditions change from the low/intermediate temperatures and high stress regime, to the high temperature and low stress regime. When anti phase boundary (APB) shearing and Orowan looping mechanisms operate, the classical picture is that, at for a given volume fraction, the bypass mechanism changes from shearing to looping with increased particle size and within a broad coexistence size window. Another possibility, which is supported by indirect experimental evidence, is that a third “hybrid” transition mechanism may operate. Here, in this paper we use discrete dislocation dynamics (DDD) simulations to study dislocation bypass mechanisms in Ni-based superalloys. We develop a new method to compute generalized stacking fault forces in DDD simulations, based on a concept borrowed from complex analysis and known as the winding number of a closed curve about a point. We use this method to study the mechanisms of bypass of a square lattice of spherical $\gamma'$ precipitates by $a/2\langle110\rangle$ {111} edge dislocations, as a function of the precipitates volume fraction and size. We show that not only the hybrid mechanism is possible, but also that it operates as the transition mechanism between the shearing and looping regimes over a wide range of precipitates volume fraction and radii. Based on our simulation results, we propose a simple model for the strength of this mechanism. We also consider the effects of a $\gamma$/$\gamma'$ lattice misfit on the bypass mechanisms, which we approximate by an additional precipitate stress computed according to Eshelby’s inclusion theory. We show that in the shearing and hybrid looping-shearing regimes, a lattice misfit generally results in an increased bypass stress. For sufficiently high lattice misfit, the critical bypass configuration in attractive dislocation-precipitates interactions changes dramatically, and the bypass stress is controlled by the pinning of the trailing dislocation on the exit side of the precipitates, similar to what has been reported in the high-temperature creep literature.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Molecular dynamics study of Si(100)-oxidation: SiO and Si emissions from Si/SiO{sub 2} interfaces and their incorporation into SiO{sub 2}

Dynamics of Si(100)-oxidation processes at the Si/SiO{sub 2} interface and in the SiO{sub 2} region are investigated focusing on SiO and Si emissions from the interface and the following incorporation into the SiO{sub 2} and/or substrate. Classical molecular dynamics (MD) simulations with variable charge interatomic potentials are performed to clarify these atomic processes. By incorporating oxygen atoms, two-folded Si atoms are formed after structural relaxation at the interface and are emitted as SiO molecules into SiO{sub 2}. The energy barrier of the SiO emission is estimated to be 1.20 eV on the basis of the enthalpy change in an MD simulation. The emitted SiO molecule is incorporated into the SiO{sub 2} network through a Si-O rebonding process with generating an oxygen vacancy. The energy barrier of the SiO incorporation is estimated to be 0.79–0.81 eV. The elementary process of oxygen vacancy diffusion leading to the complete SiO incorporation is also simulated, and the energy barriers are found to be relatively small, 0.71–0.79 eV. The energy changes of Si emissions into the substrate and SiO{sub 2} are estimated to be 2.97–7.81 eV, which are larger than the energy barrier of the SiO emission. This result suggests that, at the ideally flat Si/SiO{sub 2} interface, the SiO emission into the SiO{sub 2} region occurs prior to the Si emission, which is consistent with previous theoretical and experimental studies. The above mentioned typical atomic processes are successfully extracted from some (or one) of MD simulations among many trials in which a statistical procedure is partly employed. Our results give a unified understanding of Si oxidation processes from an atomistic point of view.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗