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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 955 records · Page 53

All sources NASA NTRS DOE OSTI

DOE OSTI2017

Numerical Analysis of 2-D and 3-D MHD Flows Relevant to Fusion Applications

Here, the analysis of many fusion applications such as liquid-metal blankets requires application of computational fluid dynamics (CFD) methods for electrically conductive liquids in geometrically complex regions and in the presence of a strong magnetic field. A current state of the art general purpose CFD code allows modeling of the flow in complex geometric regions, with simultaneous conjugated heat transfer analysis in liquid and surrounding solid parts. Together with a magnetohydrodynamics (MHD) capability, the general purpose CFD code will be a valuable tool for the design and optimization of fusion devices. This paper describes an introduction of MHD capability into the general purpose CFD code CFX, part of the ANSYS Workbench. The code was adapted for MHD problems using a magnetic induction approach. CFX allows introduction of user-defined variables using transport or Poisson equations. For MHD adaptation of the code three additional transport equations were introduced for the components of the magnetic field, in addition to the Poisson equation for electric potential. The Lorentz force is included in the momentum transport equation as a source term. Fusion applications usually involve very strong magnetic fields, with values of the Hartmann number of up to tens of thousands. In this situation a system of MHD equations become very rigid with very large source terms and very strong variable gradients. To increase system robustness, special measures were introduced during the iterative convergence process, such as linearization using source coefficient for momentum equations. The MHD implementation in general purpose CFD code was tested against benchmarks, specifically selected for liquid-metal blanket applications. Results of numerical simulations using present implementation closely match analytical solutions for a Hartmann number of up to 1500 for a 2-D laminar flow in the duct of square cross section, with conducting and nonconducting walls. Results for a 3-D test case are also included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Highly tail-asymmetric lipids interdigitate and cause bidirectional ordering

Phospholipids form structurally and compositionally diverse membranes. A less studied type of compositional diversity involves phospholipid tail variety. Some phospholipids contain two acyl tails which differ in length. These tail-asymmetric lipids are shown to contribute to temperature sensitivity, oxygen adaptability, and membrane fluidity. Membranes of a highly virulent intracellular bacterium, Francisella tularensis, contain highly tail-asymmetric 1-lignoceroyl-2-decanoyl-sn-glycero-3-phosphatidylethanolamine (XJPE) lipids which were previously shown to inhibit inflammatory responses in host cells. XJPE tails have unusually high asymmetry, and how they contribute to membrane properties on a molecular level is unknown. Here, we use small angle X-ray scattering and molecular dynamics simulations to investigate how varying XJPE ratios alters properties of simple membranes. Our results demonstrate that at high concentration they promote liquid-to-gel transition in otherwise liquid membranes, while at low concentration they are tolerated well, minimally altering membrane properties. In liquid membranes, XJPE lipids dynamically adopt two main conformations; with the long tail extended into the opposing leaflet or bent-back residing in its own leaflet. When added to both leaflets XJPE primarily adopts an extended confirmation, while asymmetric addition results in more bent-back orientations. The former increases tail ordering and the latter decreases it. XJPE tails adopt different conformations that induce composition- and leaflet-dependent bidirectional effect on membrane fluidity and this suggests that Francisella tularensis could use tail asymmetry to facilitate vesicle fusion and destabilize host cells. The effect of tail-asymmetric lipids on complex membranes should be further investigated to reveal the regulatory roles of high tail asymmetry.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

Tuning the melting point of selected ionic liquids through adjustment of the cation’s dipole moment.

In prior research with thermally robust salts [Cassity et al., PCCP, 2017, 19, 31560] it was noted that an increase in the dipole moment of the cation generally led to a decrease in the melting point. Molecular dynamics simulations of the liquid state revealed that an increased dipole moment reduces cation-cation repulsions through dipole-dipole alignment. This was believed to reduce the liquid phase enthalpy, which would tend to lower the melting point of the IL. Here, we also test this principle by replacing hydrogen atoms with fluorine atoms at selected positions within the cation. This allows us to alter the electrostatics of the cation without substantially affecting the sterics. Furthermore, the strength of the dipole moment can be controlled by choosing different positions within the cation for replacement. We studied variants of four different parent cations paired with bistriflimide and determined their melting points, and enthalpies and entropies of fusion through DSC experiments. The decreases in the melting point were determined to be enthalpically driven. We found that the dipole moment of the cation, as determined by quantum chemical calculations, is inversely correlated with the melting point of the given compound. Molecular dynamics simulations of the crystalline and solid states of two isomers showed differences in their enthalpies of fusion that closely matched those seen experimentally. Moreover, this reduction in the enthalpy of fusion was determined to be caused by an increase in the enthalpy of the crystalline state. We provide evidence that dipole-dipole interactions between cations leads to the formation of cationic domains in the crystalline state. These cationic associations partially block favourable cation–anion interactions, which are recovered upon melting. If, however, the dipole- dipole interactions between cations is too strong they have a tendency to form glasses. This study provides a design rule for lowering the melting point of structurally similar ILs by altering their dipole moment.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2016

Transport properties and Stokes-Einstein relation in Al-rich liquid alloys

We use ab initio molecular dynamics simulations to study the transport properties and the validity of the Stokes-Einstein relation in Al-rich liquid alloys with Ni, Cu, and Zn as alloying elements. First, we show that the composition and temperature dependence of their transport properties present different behaviors, which can be related to their local structural ordering. Then, we evidence that the competition between the local icosahedral ordering and the local chemical ordering may cause the breakdown of the Stokes–Einstein relation even in the liquid phase. We demonstrate that this breakdown can be captured by entropy-scaling relationships developed by Rosenfeld and using the two-body excess entropy. Our findings provide a unique framework to study the relation between structure, thermodynamics, and dynamics in metallic melts and pave the way towards the explanation of various complex transport properties in metallic melts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Atomic dynamics in fluids: Normal mode analysis revisited

Developing microscopic understanding of the thermal properties of liquids is challenging due to their strong dynamic disorder, which prevents characterization of the atomic degrees of freedom. There have been significant research interests in the past few decades to extend the normal mode analysis for solids to instantaneous structures of liquids. However, the nature of normal modes that arise from these unstable structures is still elusive. Here, in this paper, we explore the instantaneous eigenmodes of dynamical matrices of various Lennard-Jones argon liquid and gas systems at high temperatures and show that the normal modes can be interpreted as an interpolation of T →∞ (gas) and T=0 (solid) mode descriptions. We find that normal modes become increasingly collisional and translational, recovering atomistic gaslike behavior rather than vibrational with increase in temperature, suggesting that normal modes in liquids may be described by both solidlike and gaslike modes.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2019

Effect of an external electric field on the dynamics and intramolecular structures of ions in an ionic liquid

Simulations of the ionic liquid 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide in an external electric field have been performed using a Drude particle polarizable force field. The structure of the ions has been analyzed, with close attention paid to the configurations of the ions. The “charge arm” concept is used to explain certain changes of these configurations that would be difficult to rationalize otherwise, e.g., trans → cis isomerization of the bis(trifluoromethylsulfonyl)imide anion and extension of the alkyl chain of the cation. It has also been shown that the ions orient themselves so that their charge arms align with and stretch out along the field, and these effects occur at lower external electric field strengths than cause a change in the inherent diffusion of the ions. Here, the dynamics of the system parallel and perpendicular to the field were analyzed, and it was found that the applied field affected the diffusion normal to the field. This is explained as a secondary effect of a change in the ion cage formed by the surrounding counterions of a given ion in the ionic liquid. The breakdown of the ion cages was rationalized by correlating changes to the inherent diffusion of the ions with other changes to the diffusion and bulk structure of the liquid, as well as considering the average forces on the ions compared to the force the ions would be expected to experience in an electric field. Parallel to the field, a drift was observed at every electric field studied. In electric fields with no changes to the ion cage structure, the relationship between the drift and electric field was found to be linear, becoming nonlinear as the ion cage structure breaks down.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRSDate not supplied

Principles of Fluid Dynamic Similarity Analysis for Slosh Experiments

Most tank slosh test programs are not conducted with flight hardware in a flight-like environment due to cost and schedule constraints. Instead, subscale tanks are typically tested on the ground. If the “real”, flight liquid in the tank is a cryogenic propellant, sometimes a cryogen is used during testing, but more often a simulant liquid is used. Fluid dynamic similarity analysis, or “scaling analysis” as it is sometimes called, for slosh examines nondimensional numbers to determine which fluid dynamic regime(s) a tank is operating in. If the relevant nondimensional numbers of a scale tank test are similar enough to the real tank’s nondimensional numbers, then it is likely that the scale test data will be similar to what would be obtained with the real tank. This type of analysis is useful for guiding experiment design, test liquid selection, and for understanding the physical relevance of prior test data to a specific application. The authors have performed scaling analyses for many slosh test programs conducted in ground, parabolic flight, suborbital, and in-orbit environments. Basics of how to perform the analysis, along with select examples from ground and microgravity environments, are presented. Suggestions for various assumptions and compromises to make are discussed. Limitations and challenges are also discussed.

slosh↗

DOE OSTI2015

A Comparison of the Predictive Capabilities of the Embedded-Atom Method and Modified Embedded-Atom Method Potentials for Lithium

Here, we compare six lithium potentials by examining their ability to predict coexistence properties and liquid structure using molecular dynamics. All potentials are of the embedded-atom-method (EAM) type. The coexistence properties we focus on are the melting curve, vapor pressure, saturated liquid density, and vapor-liquid surface tension. For each property studied, the simulation results are compared to available experimental data in order to properly assess the accuracy of each potential. We find that the Cui 2NN MEAM is the most robust potential, giving adequate agreement with most of the properties examined. For example, the zero-pressure melting point of this potential is shown to be around 443 K, while experimentally is it about 454 K. This potential also gives excellent agreement with saturated liquid densities, even though no liquid properties were used in the fitting procedure. Our study allows us to conclude that the Cui 2NN MEAM should be used for further simulations of lithiums.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

A polarizable QM/MM approach to the molecular dynamics of amide groups solvated in water

The infrared (IR) spectra of polypeptides are dominated by the so-called amide bands. Because they originate from the strongly polar and polarizable amide groups (AGs) making up the backbone, their spectral positions sensitively depend on the local electric fields. Aiming at accurate computations of these IR spectra by molecular dynamics (MD) simulations, which derive atomic forces from a hybrid quantum and molecular mechanics (QM/MM) Hamiltonian, here we consider the effects of solvation in bulk liquid water on the amide bands of the AG model compound N-methyl-acetamide (NMA). As QM approach to NMA we choose grid-based density functional theory (DFT). For the surrounding MM water, we develop, largely based on computations, a polarizable molecular mechanics (PMM) model potential called GP6P, which features six Gaussian electrostatic sources (one induced dipole, five static partial charge distributions) and, therefore, avoids spurious distortions of the DFT electron density in hybrid DFT/PMM simulations. Bulk liquid GP6P is shown to have favorable properties at the thermodynamic conditions of the parameterization and beyond. Lennard-Jones (LJ) parameters of the DFT fragment NMA are optimized by comparing radial distribution functions in the surrounding GP6P liquid with reference data obtained from a “first-principles” DFT-MD simulation. Finally, IR spectra of NMA in GP6P water are calculated from extended DFT/PMM-MD trajectories, in which the NMA is treated by three different DFT functionals (BP, BLYP, B3LYP). Method-specific frequency scaling factors are derived from DFT-MD simulations of isolated NMA. The DFT/PMM-MD simulations with GP6P and with the optimized LJ parameters then excellently predict the effects of aqueous solvation and deuteration observed in the IR spectra of NMA. As a result, the methods required to accurately compute such spectra by DFT/PMM-MD also for larger peptides in aqueous solution are now at hand.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Activated relaxation in supercooled monodisperse atomic and polymeric WCA fluids: Simulation and ECNLE theory

Here, we combine simulation and Elastically Collective Nonlinear Langevin Equation (ECNLE) theory to study the activated relaxation in monodisperse atomic and polymeric Weeks–Chandler–Andersen (WCA) liquids over a wide range of temperatures and densities in the supercooled regime under isochoric conditions. By employing novel crystal-avoiding simulations, metastable equilibrium dynamics is probed in the absence of complications associated with size polydispersity. Based on a highly accurate structural input from integral equation theory, ECNLE theory is found to describe well the simulated density and temperature dependences of the alpha relaxation time of atomic fluids using a single system-specific parameter, a c , that reflects the nonuniversal relative importance of local cage and collective elastic barriers. For polymer fluids, the explicit dynamical effect of local chain connectivity is modeled at the fundamental dynamic free energy trajectory level based on a different parameter, N c , that quantifies the degree of intramolecular correlation of bonded segment activated barrier hopping. For the flexible chain model studied, a physically intuitive value of N c ≈ 2 results in good agreement between simulation and theory. A direct comparison between atomic and polymeric systems reveals that chain connectivity can speed up activated segmental relaxation due to weakening of equilibrium packing correlations but can slow down relaxation due to local bonding constraints. The empirical thermodynamic scaling idea for the alpha time is found to work well at high densities or temperatures but fails when both density and temperature are low. The rich and subtle behaviors revealed from simulation for atomic and polymeric WCA fluids are all well captured by ECNLE theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Fluid Flow Observations of The Spray Near-Field Using High-Speed X-Ray Imaging

Sprays are important to may industrial processes, but the near-field region where spray formation takes place is difficult to visualize and characterize through visible light imaging because this region is typically optically dense. X-ray imaging is one technique that can penetrate this optically dense region and provide unique observations of the spray formation process. This study uses high-speed white beam X-ray imaging conducted at the Advanced Photon Source to capture spray formation dynamics from a canonical airblast atomizer. In this work, a range of momentum flux ratios are visualized by varying the gas flow rate while the liquid flow rate is held constant. The use of X-rays show overlapping flow structures including bags, ligaments, mushrooms, and webs. Air bubbles are captured within the liquid regions, and a very dynamic crown forms at the liquid needle exit at higher momentum flux ratios. Finally, at the highest momentum flux ratios, air bubbles and liquid droplets are observed migrating upstream in the center of the crown, and this motion is attributed to a toroidal recirculation cell at the nozzle exit.

47 OTHER INSTRUMENTATION↗

DOE OSTI2014

Ab initio study of electron-ion structure factors in binary liquids with different types of chemical bonding

Electron-ion structure factors, calculated in ab initio molecular dynamics simulations, are reported for several binary liquids with different kinds of chemical bonding: metallic liquid alloy Bi–Pb, molten salt RbF, and liquid water. We derive analytical expressions for the long-wavelength asymptotes of the partial electron-ion structure factors of binary systems and show that the analytical results are in good agreement with the ab initio simulation data. The long-wavelength behaviour of the total charge structure factors for the three binary liquids is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1990

A high-speed photography study of cavitation in a dynamically loaded journal bearing

The earlier study made by Jacobson and Hamrock on the cavitation of liquid lubricant films in a dynamically loaded journal bearing was repeated with a quartz sleeve, which was more rigid than the Polymethylmethacrylate (PMMA) sleeve used previously. Various improvements of the test rig were made concomitantly so that the experimental errors could be better controlled and assessed. The updated speed photography experiment and its results are described.

Sun, D. C.↗

NASA NTRS1991

A high speed photography study of cavitation in a dynamically loaded journal bearing

Sun, D. C.↗

DOE OSTI2014

A molecular dynamics study of nuclear quantum effect on the diffusion of hydrogen in condensed phase

In this paper, the quantum effect of hydrogen molecule on its diffusivity is analyzed using Molecular Dynamics (MD) method. The path integral centroid MD (CMD) method is applied for the reproduction method of time evolution of the molecules. The diffusion coefficient of liquid hydrogen is calculated using the Green-Kubo method. The simulation is performed at wide temperature region and the temperature dependence of the quantum effect of hydrogen molecule is addressed. The calculation results are compared with those of classical MD results. As a result, it is confirmed that the diffusivity of hydrogen molecule is changed depending on temperature by the quantum effect. It is clarified that this result can be explained that the dominant factor by quantum effect on the diffusivity of hydrogen changes from the swollening the potential to the shallowing the potential well around 30 K. Moreover, it is found that this tendency is related to the temperature dependency of the ratio of the quantum kinetic energy and classical kinetic energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

The Proton Density of States in Confined Water (H2O)

The hydrogen density of states (DOS) in confined water has been probed by inelastic neutron scattering spectra in a wide range of its P–T phase diagram. The liquid–liquid transition and the dynamical crossover from the fragile (super-Arrhenius) to strong (Arrhenius) glass forming behavior have been studied, by taking into account the system polymorphism in both the liquid and amorphous solid phases. The interest is focused in the low energy region of the DOS ( E < 10 meV) and the data are discussed in terms of the energy landscape (local minima of the potential energy) approach. In this latest research, we consider a unit scale energy (EC) linked to the water local order governed by the hydrogen bonding (HB). All the measured spectra, scaled according to such energy, evidence a universal power law behavior with different exponents ( γ ) in the strong and fragile glass forming regions, respectively. In the first case, the DOS data obey the Debye squared-frequency law, whereas, in the second one, we obtain a value predicted in terms of the mode-coupling theory (MCT) ( γ ≃ 1.6 ).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Insertion of liquid crystal molecules into hydrocarbon monolayers

Atomistic molecular dynamics simulations were carried out to investigate the molecular mechanisms of vertical surface alignment of liquid crystals. We study the insertion of nCB (4-Cyano-4{sup ′}-n-biphenyl) molecules with n = 0,…,6 into a bent-core liquid crystal monolayer that was recently found to provide good vertical alignment for liquid crystals. The results suggest a complex-free energy landscape for the liquid crystal within the layer. The preferred insertion direction of the nCB molecules (core or tail first) varies with n, which can be explained by entropic considerations. The role of the dipole moments was found to be negligible. As vertical alignment is the leading form of present day liquid crystal displays (LCD), these results will help guide improvement of the LCD technology, as well as lend insight into the more general problem of insertion of biological and other molecules into lipid and surfactant layers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Computing the viscous effect in early-time drop impact dynamics

The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. ( Phys. Rev. Lett. , vol. 112, no. 13, 2014 b , p. 134501) show that the liquid–air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier–Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square-root scaling relationship. We also predict the effect of interfacial tension at the liquid–gas interface on the drop impact, showing that it mediates the lift-off behaviour.

42 ENGINEERING↗