Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Liquid dynamics”

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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Role of Surface Features on the Initial Dissolution of CH 3 NH 3 PbI 3 Perovskite in Liquid Water: An Ab Initio Molecular Dynamics Study

In this study, the degradation of CH 3 NH 3 PbI 3 (MAPbI 3 ) hybrid organic inorganic perovskite (HOIP) by water has been the major issue hampering its use in commercial perovskites solar cells (PSCs) as MAPbI 3 HOIP have been known to easily degrade in the presence of water molecules. However, even though there have been numerous studies investigating this phenomenon, there is still no consensus on the mechanisms of initial stages of dissolution. Here, we attempt to consolidate differing mechanistic interpretations previously reported in the literature through the use of the first-principles constrained ab-initio molecular dynamics (AIMD) to study the mechanisms, kinetics, and thermodynamics that accompany the degradation of MAPbI 3 HOIP in liquid water. We consider not only the dissolution of the species found on the pristine MAPbI 3 HOIP surfaces, but also the dissolution from defect sites to imitate the successive degradation steps, and propose a sequence of events in the initial phase of MAPbI 3 HOIP dissolution. By comparing the dissolution free energy barrier between surface species of different surficial types, we find that the dominant dissolution mechanisms of surface species varies widely based on the specific surface features adjacent to the dissolving ion. The high sensitivity of dissolution mechanism to surface features has contributed to the many dissolution mechanisms proposed in the literature. In contrast, the dissolution free energy barriers are mainly determined by the dissolving species rather than the type of surfaces, and the type of surfaces the ions are dissolving from are inconsequential toward the dissolution free energy barrier. However, the presence of surface defects such as vacancy sites are found to significantly lower the dissolution free energy barriers. Based on the estimated dissolution free energy barriers from different types of surfaces that we investigated in this study, we proposed that the dissolution of MAPbI 3 HOIP in liquid water originates from surface defect sites that propagate laterally along the surface layer of the MAPbI 3 HOIP crystal.

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

Geometric and energetic considerations of surface fluctuations during ion transfer across the water-immiscible organic liquid interface

Molecular dynamics simulations and umbrella sampling free energy calculations are used to examine the thermodynamics, energetics, and structural fluctuations that accompany the transfer of a small hydrophilic ion (Cl{sup −}) across the water/nitrobenzene interface. By examining several constrained interface structures, we isolate the energetic costs of interfacial deformation and co-transfer of hydration waters during the ion transfer. The process is monitored using both energy-based solvation coordinates and a geometric coordinate recently introduced by Morita and co-workers to describe surface fluctuations. Our simulations show that these coordinates provide a complimentary description of the water surface fluctuations during the transfer and are necessary for elucidating the mechanism of the ion transfer.

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

Slosh wave and geyser excitations due to liquid hydrogen shut-off during draining in microgravity

The dynamical behavior of liquid hydrogen shut-off during draining, and shut-off at the moment of the incipience of a suction dip have been investigated. It shows that a large amplitude surge is observed for liquid in the container at the moment of liquid hydrogen shut-off in reduced gravity. It also shows that slosh waves accompanied by a strong geyser are developed for surge-related flow fields induced by liquid hydrogen shut-off at the incipience of a suction dip. In the slosh wave excitation, both a lower gravity environment and higher flow rate before the shut-off of liquid draining are resonsible for the initiation of greater amplitude slosh waves. Slosh wave excitation, due to shut-off during liquid hydrogen draining, shift the fluid mass distribution in the container which imposes time-dependent variation in the spacecraft moment of inertia. This provides important information necessary for on-orbit guidance and attitude control of spacecraft.

Hung, R. J.↗
DOE OSTI2013

Velocity measurements and free-shear-layer instabilities in a rotating liquid metal

The dynamics of rotating flows in magnetohydrodynamics (MHD) are expected to play an important role in astrophysical accretion disks, specifically via the generation of the magnetorotational instability (MRI). The Princeton MRI experiment, a Taylor-Couette device with split axial endcap rings and with a gallium-indium-tin working fluid, was constructed to study rotating MHD flows in the laboratory. This work uses an ultrasound Doppler velocimetry (UDV) diagnostic to measure internal fluid velocities of the mean background flow and also of fluctuations driven by instabilities in this experiment. Mean azimuthal velocities have been measured to show the effect of an applied axial magnetic field. For moderate magnetic field strengths, with the Elsasser number less than 1, the azimuthal velocity at a point, normalized to the inner cylinder velocity, is constant for a given Alfven speed normalized to the inner cylinder velocity. In the Elsasser much greater than 1 regime, the shear at the split in the endcaps extends axially as a free Shercliff layer through the full height of the experiment. When the inner endcap ring rotates faster than the outer endcap ring, the free shear layer produces a Kelvin-Helmholtz instability with an eigenmode that fills the experimental volume when Elsasser is greater than 1, or when there is sufficient background rotation so that the Rossby number is greater than 2.35. When the inner endcap ring rotates slower than the outer endcap ring, the Kelvin-Helmholtz instability is present in the absence of an applied magnetic field or background rotation. Using a linear eigenmode analysis of a free shear layer, the threshold for the appearance of the Kelvin-Helmholtz instability is presented as a competition between Rayleigh's centrifugal instability, which is driven by a negative radial gradient in the specific angular momentum, and the Kelvin-Helmholtz instability, which is driven by the velocity shear. Both instabilities act to smooth out the shear layer. While the eigenmodes of the centrifugal instability are confined to the free shear layer, allowing quiet flow in the bulk of the fluid when that instability dominates, the eigenmodes of the Kelvin-Helmholtz instability fill the fluid volume, creating large-amplitude coherent velocity fluctuations throughout.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Structure and lifetimes in ionic liquids and their mixtures

With the aid of molecular dynamics simulations, we study the structure and dynamics of different ionic liquid systems, with focus on hydrogen bond, ion pair and ion cage formation. To do so, we report here radial distribution functions, their number integrals, and various time–correlation functions, from which we extract well-defined lifetimes by means of the reactive flux formalism. We explore the influence of polarizable force fields vs. non-polarizable ones with downscaled charges (±0.8) for the example of 1-butyl-3-methylimidazolium bromide. Furthermore, we use 1-butyl-3-methylimidazolium trifluoromethanesulfonate to investigate the impact of temperature and mixing with water as well as with the chloride ionic liquid. Smaller coordination numbers, larger distances, and tremendously accelerated dynamics are observed when the polarizable force field is applied. The same trends are found with increasing temperature. Adding water decreases the ion–ion coordination numbers whereas the water–ion and water–water coordination is enhanced. A domain analysis reveals that the nonpolar parts of the ions are dispersed and when more water is added the water clusters increase in size. The dynamics accelerate in general upon addition of water. In the ionic liquid mixture, the coordination number around the cation changes between the two anions, but the number integrals of the cation around the anions remain constant and the dynamics slow down with increasing content of the chloride ionic liquid.

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

Level crossing analysis of chemically induced dynamic nuclear polarization: Towards a common description of liquid-state and solid-state cases

Chemically Induced Dynamic Nuclear Polarization (CIDNP) is an efficient method of creating non-equilibrium polarization of nuclear spins by using chemical reactions, which have radical pairs as intermediates. The CIDNP effect originates from (i) electron spin-selective recombination of radical pairs and (ii) the dependence of the inter-system crossing rate in radical pairs on the state of magnetic nuclei. The CIDNP effect can be investigated by using Nuclear Magnetic Resonance (NMR) methods. The gain from CIDNP is then two-fold: it allows one to obtain considerable amplification of NMR signals; in addition, it provides a very useful tool for investigating elusive radicals and radical pairs. While the mechanisms of the CIDNP effect in liquids are well established and understood, detailed analysis of solid-state CIDNP mechanisms still remains challenging; likewise a common theoretical frame for the description of CIDNP in both solids and liquids is missing. Difficulties in understanding the spin dynamics that lead to the CIDNP effect in the solid-state case are caused by the anisotropy of spin interactions, which increase the complexity of spin evolution. In this work, we propose to analyze CIDNP in terms of level crossing phenomena, namely, to attribute features in the CIDNP magnetic field dependence to Level Crossings (LCs) and Level Anti-Crossings (LACs) in a radical pair. This approach allows one to describe liquid-state CIDNP; the same holds for the solid-state case where anisotropic interactions play a significant role in CIDNP formation. In solids, features arise predominantly from LACs, since in most cases anisotropic couplings result in perturbations, which turn LCs into LACs. We have interpreted the CIDNP mechanisms in terms of the LC/LAC concept. This consideration allows one to find analytical expressions for a wide magnetic field range, where several different mechanisms are operative; furthermore, the LAC description gives a way to determine CIDNP sign rules. Thus, LCs/LACs provide a consistent description of CIDNP in both liquids and solids with the prospect of exploiting it for the analysis of short-lived radicals and for optimizing the polarization level.

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

Stockmayer fluid simulations for viscosity and glass transition temperature of ionic liquids

We develop a Stockmayer fluid model for molecular dynamics simulations of ionic liquids that captures molecular polarization, ionic conductivity, viscosity, and glass transition temperature, using ethylammonium nitrate (EAN) as an example. The ions in EAN are treated as spheres interacting via the Lennard-Jones potential with an embedded point charge and a permanent dipole moment. We show that our simulation results for EAN are consistent with experimental data and then explore the effects of the molecular parameters on the viscosity of ionic liquids. Our results indicate that viscosity monotonically increases with ionic charge and dipole moment but non-monotonically changes with ionic diameter (or molar volume). This non-monotonic trend arises from the competition among the electrostatic interactions, molecular packing, and size asymmetry between the cation and anion. In conclusion, our model also shows that long-lived ion pairs result in higher viscosities.

Coarse-grained simulations↗
DOE OSTI2017

Quantum Dynamics and Spectroscopy of Ab Initio Liquid Water: The Interplay of Nuclear and Electronic Quantum Effects

Understanding the reactivity and spectroscopy of aqueous solutions at the atomistic level is crucial for the elucidation and design of chemical processes. However, the simulation of these systems requires addressing the formidable challenges of treating the quantum nature of both the electrons and nuclei. Exploiting our recently developed methods that provide acceleration by up to 2 orders of magnitude, we combine path integral simulations with on-the-fly evaluation of the electronic structure at the hybrid density functional theory level to capture the interplay between nuclear quantum effects and the electronic surface. Here we show that this combination provides accurate structure and dynamics, including the full infrared and Raman spectra of liquid water. This allows us to demonstrate and explain the failings of lower-level density functionals for dynamics and vibrational spectroscopy when the nuclei are treated quantum mechanically. Finally, these insights thus provide a foundation for the reliable investigation of spectroscopy and reactivity in aqueous environments.

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

Pairing mechanism among ionic liquid ions in aqueous solutions. A molecular dynamics study

In this study, we carried out molecular dynamics simulations to examine the molecular mechanism for ionic liquid pair association in aqueous solutions. We chose the commonly studied imidazolium-based ionic liquid pairs. We computed potentials of mean force (PMF) for four systems—1,3-dimethlylimidazoliumchloride; 1,3-dimethlylimidazolium iodide; 1-methly-3-octylimidazolium chloride; and 1-methly-3-octylimidazolium iodide. Our PMF studies show a stronger interaction for the ion pairs of systems involving dimethlylimidazolium as the cation species compared to that of the systems containing octylimidazolium. This result indicates a decrease in ion-pair association as the cation alkyl tail length increases. We also studied the kinetics of ion-pair dissociation using different rate theories such as the Grote-Hynes and Kramer’s theories. As expected, the computed rate results significantly deviated from results obtained from transition state theory because it does not account for dynamical solvent effects. Dissociative barrier curvatures are found to be very small for the systems investigated because the transmission coefficients computed using Grote-Hynes theory and Kramer’s theory are approximately equal. Our analysis of the rotational dynamics of cations revealed that the time scales for molecular reorientation are longer for cations with longer alkyl tails. This work was supported by the U.S. Department of Energy (DOE), Office of Basic Energy Sciences (BES), Division of Chemical Sciences, Geosciences, and Biosciences. Pacific Northwest National Laboratory is a multiprogram national laboratory operated for DOE by Battelle. The calculations were carried out using computer resources provided by BES.

Annapureddy, Harsha V.↗
NASA NTRSDate not supplied

Time Resolved Visualization of Liquid Jet Interaction With H2-Air Rotating Detonations Using MHz Rate Diesel PLIF

The characterization of the dynamic response of liquid jets to transient detonation wave passage is critical for optimization and modeling of liquid fueled rotating detonation combustors. In this work, a rotating detonation combustor (RDC) is operated on hydrogen and air to sustain stable detonation waves that acts as a detonation driver and interact in a one-way coupled manner with a single liquid fuel jet that propagates into the combustion chamber with cycle periods of ~ 250 μs. Diesel is utilized as a realistic fuel surrogate with higher aromatic compounds to enable fluorescence excitation, using the 355 nm thirdharmonic output of a burst-mode Nd:YAG laser, imaged at rates up to 1 MHz. By optimizing the technique to accommodate orders of magnitude variations in the fuel density throughout the injection process, the PLIF data enable quantitative measurements including the refill time, the relative recovery between liquid and gaseous jets and jet trajectory across various momentum flux ratio. As the passage of the detonation wave imparts significant changes in the momentum flux ratio, the qualitative liquid break-up process and spatial distribution varies significantly in time. As the injection system recovers ~ 70% of the cycle period and return to a quasi-steady position and allow comparisons with theoretical jet trajectories. These data, enabled by ultra-high-speed PLIF imaging, represent some of the first detailed measurements for quantifying the dynamic response and recovery of liquid jets exposed to periodic detonations in an operating RDC.

Propulsion↗
DOE OSTI2020

Tuning the dynamics of imidazolium-based ionic liquids via hydrogen bonding. I. The viscous regime

Combining results from impedance spectroscopy and oscillatory shear rheology, this study focuses on the relation between the mass and charge flows and on how these are affected by the H-bonding in viscous ionic liquids (ILs). In particular, we compare the relaxational behaviors of the paradigmatic IL 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (EMIM-TFSI) and its OH-functionalized counterpart 1-(2-hydroxyethyl)-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (OHEMIM-TFSI). Our results and their analysis demonstrate that the presence of cationic OH-groups bears a strong impact on the overall dynamics of OHEMIM-TFSI, although no signatures of suprastructural relaxation modes could be identified in their dielectric and mechanical responses. To check whether at the origin of this strong variation is the H-bonding or merely the difference between the corresponding cation sizes (controlling both the hydrodynamic volume and the inter-charge distance), the present study includes 1-propyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide (PMIM-TFSI), mixtures of EMIM-TFSI and PMIM-TFSI with lithium bis(trifluoromethylsulfonyl)imide (Li-TFSI), and mixtures of OHEMIM-TFSI with PMIM-TFSI. Their investigation clearly reveals that the dynamical changes induced by H-bonding are significantly larger than those that can be attributed to the change in the ion size. Moreover, in the mixtures of OHEMIM-TFSI with PMIM-TFSI, a dilution of the OH-groups leads to strong deviations from ideal mixing behavior, thus highlighting the common phenomenological ground of hydroxy-functionalized ILs and other H-bonded liquids.

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

Molecular simulation study of dynamical properties of room temperature ionic liquids with carbon pieces

Room temperature ionic liquids (RTILs) with dispersed carbon pieces exhibit distinctive physiochemical properties. In order to explore the molecular mechanism, RTILs/carbon pieces mixture we investigated it by molecular dynamics (MD) simulation in this work. Rigid and flexible carbon pieces in the form of graphene with different thicknesses and carbon nanotubes in different sizes were dispersed in a representative RTIL 1-butyl-3-methyl-imidazolium dicyanamide ([Bmim][DCA]). Our study demonstrated that the diffusion coefficients of RTILs in the presence of flexible carbons are similar to those of bulk RTILs at varying temperatures, which is in contrast to the decreased diffusion of RTILs in the presence of rigid carbons. In addition, interfacial ion number density at rigid carbon surfaces was higher than that at flexible ones, which is correlated with the accessible external surface area of carbon pieces. The life time of cation-anion pair in the presence of carbon pieces also exhibited a dependence on carbon flexibility. RTILs with dispersed rigid carbon pieces showed longer ion pair life time than those with flexible ones, in consistence with the observation in diffusion coefficients. Furthermore, this work highlights the necessity of including the carbon flexibility when performing MD simulation of RTILs in the presence of dispersed carbon pieces in order to obtain the reliable dynamical and interfacial structural properties.

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

In Situ Synchrotron X-ray Study of Ultrasound Cavitation and Its Effect on Solidification Microstructures

Considerable progress has been made in studying the mechanism and effectiveness of using ultrasound waves to manipulate the solidification microstructures of metallic alloys. However, uncertainties remain in both the underlying physics of how microstructures evolve under ultrasonic waves, and the best technological approach to control the final microstructures and properties. We used the ultrafast synchrotron X-ray phase contrast imaging facility housed at the Advanced Photon Source, Argonne National Laboratory, US to study in situ the highly transient and dynamic interactions between the liquid metal and ultrasonic waves/bubbles. The dynamics of ultrasonic bubbles in liquid metal and their interactions with the solidifying phases in a transparent alloy were captured in situ. The experiments were complemented by the simulations of the acoustic pressure field, the pulsing of the bubbles, and the associated forces acting onto the solidifying dendrites. The study provides more quantitative understanding on how ultrasonic waves/bubbles influence the growth of dendritic grains and promote the grain multiplication effect for grain refinement.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Liquid methanol under a static electric field

We report on an ab initio molecular dynamics study of liquid methanol under the effect of a static electric field. We found that the hydrogen-bond structure of methanol is more robust and persistent for field intensities below the molecular dissociation threshold whose value (≈0.31 V/Å) turns out to be moderately larger than the corresponding estimate obtained for liquid water. A sustained ionic current, with ohmic current-voltage behavior, flows in this material for field intensities above 0.36 V/Å, as is also the case of water, but the resulting ionic conductivity (≈0.40 S cm{sup −1}) is at least one order of magnitude lower than that of water, a circumstance that evidences a lower efficiency of proton transfer processes. We surmise that this study may be relevant for the understanding of the properties and functioning of technological materials which exploit ionic conduction, such as direct-methanol fuel cells and Nafion membranes.

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

Non-equilibrium dynamics in disordered materials: Ab initio molecular dynamics simulations

The dynamic properties of liquid B{sub 2}O{sub 3} under pressure and highly-charged bromophenol molecule are studied by using molecular dynamics (MD) simulations based on density functional theory (DFT). Diffusion properties of covalent liquids under high pressure are very interesting in the sense that they show unexpected pressure dependence. It is found from our simulation that the magnitude relation of diffusion coefficients for boron and oxygen in liquid B{sub 2}O{sub 3} shows the anomalous pressure dependence. The simulation clarified the microscopic origin of the anomalous diffusion properties. Our simulation also reveals the dissociation mechanism in the coulomb explosion of the highly-charged bromophenol molecule. When the charge state n is 6, hydrogen atom in the hydroxyl group dissociates at times shorter than 20 fs while all hydrogen atoms dissociate when n is 8. After the hydrogen dissociation, the carbon ring breaks at about 100 fs. There is also a difference on the mechanism of the ring breaking depending on charge states, in which the ring breaks with expanding (n = 6) or shrink (n = 8)

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗