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 1,063 records · Page 59

Elucidating the Mechanism of Helium Evaporation from Liquid Water

Here, we investigate the evaporation of trace amounts of helium solvated in liquid water using molecular dynamics simulations and theory. Consistent with experimental observations, we find a super-Maxwellian distribution of kinetic energies of evaporated helium. This excess of kinetic energy over typical thermal expectations is explained by an effective continuum theory of evaporation based on a Fokker–Planck equation, parametrized molecularly by a potential of mean force and position-dependent friction. Using this description, we find that helium evaporation is strongly influenced by the friction near the interface, which is anomalously small near the Gibbs dividing surface due to the ability of the liquid–vapor interface to deform around the gas particle. Our reduced description provides a mechanistic interpretation of trace gas evaporation as the motion of an underdamped particle in a potential subject to a viscous environment that varies rapidly across the air–water interface. From it we predict the temperature dependence of the excess kinetic energy of evaporation, which is yet to be measured.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural connection between gallium crystals and near- T m liquids under ambient pressure

We study the short-range structural order in liquid gallium prepared by ab initio molecular dynamics simulations. We found that at 400 K, which is close to the melting point (Tm) of Ga, the dominant motif in the liquid phase is identical to that in the stable solid phase Ga-I. Strong directional bonding in this motif prevents it from nucleation in the liquid phase. Meanwhile, a newly identified motif, topologically distinct from yet related to those in β-Ga and Ga-III phases, is also abundant in Ga liquid. This new motif could serve as precursors for β-Ga or Ga-III during crystallization.

Zhang, F.↗

Hybrid nanocapsules for in situ TEM imaging of gas evolution reactions in confined liquids

Liquid cell transmission electron microscopy (TEM) enables the direct observation of dynamic physical and chemical processes in liquids at the nanoscale. Quantitative investigations into reactions with fast kinetics and/or multiple reagents will benefit from further advances in liquid cell design that facilitate rapid in situ mixing and precise control over reagent volumes and concentrations. This work reports the development of inorganic–organic nanocapsules for high-resolution TEM imaging of nanoscale reactions in liquids with well-defined zeptoliter volumes. These hybrid nanocapsules, with 48 nm average diameter, consist of a thin layer of gold coating a lipid vesicle. As a model reaction, the nucleation, growth, and diffusion of nanobubbles generated by the radiolysis of water is investigated inside the nanocapsules. When the nanobubbles are sufficiently small (10–25 nm diameter), they are mobile in the nanocapsules, but their movement deviates from Brownian motion, which may result from geometric confinement by the nanocapsules. Gases and fluids can be transported between two nanocapsules when they fuse, demonstrating in situ mixing without using complex microfluidic schemes. The ability to synthesize nanocapsules with controlled sizes and to monitor dynamics simultaneously inside multiple nanocapsules provides opportunities to investigate nanoscale processes such as single nanoparticle synthesis in confined volumes and biological processes such as biomineralization and membrane dynamics.

36 MATERIALS SCIENCE↗

Error correction in multi-fidelity molecular dynamics simulations using functional uncertainty quantification

We use functional, Fréchet, derivatives to quantify how thermodynamic outputs of a molecular dynamics (MD) simulation depend on the potential used to compute atomic interactions. Our approach quantifies the sensitivity of the quantities of interest with respect to the input functions as opposed to its parameters as is done in typical uncertainty quantification methods. We show that the functional sensitivity of the average potential energy and pressure in isothermal, isochoric MD simulations using Lennard–Jones two-body interactions can be used to accurately predict those properties for other interatomic potentials (with different functional forms) without re-running the simulations. This is demonstrated under three different thermodynamic conditions, namely a crystal at room temperature, a liquid at ambient pressure, and a high pressure liquid. The method provides accurate predictions as long as the change in potential can be reasonably described to first order and does not significantly affect the region in phase space explored by the simulation. The functional uncertainty quantification approach can be used to estimate the uncertainties associated with constitutive models used in the simulation and to correct predictions if a more accurate representation becomes available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Gravity jitter excited cryogenic liquid slosh waves in microgravity environment

The dynamical behavior of fluids, in particular the effect of surface tension on partially-filled rotating fluids (cryogenic liquid helium and helium vapor) in a full-scale Gravity Probe-B Spacecraft liquid helium container tank without probe imposed by various frequencies of gravity jitter have been investigated. Results disclose the conditions for the excitation of large amplitude slosh waves which should be avoided in the design of cryogenic liquid propellant system.

Hung, R. J.↗

Estimation of the viscosities of silicate liquids at high pressure from measurements of oxygen diffusivity

The dynamics and evolution of a magma ocean depend significantly on the structure and physical properties of silicate liquids at high pressure. Viscosity is particularly important because of its effect on the dynamics of convection and crystal setting. Traditionally, the viscosities of silicate liquids have been measured at high pressure by falling sphere viscometry but, due to some limitations of this technique, few measurements have been made at pressures exceeding 2.5 GPa. An alternative approach, which has previously been proven up to 2 GPa, is to use high pressure measurements of oxygen self-diffusivity (D) to estimate viscosity (eta) from the Eyring relationship eta = kT/D lambda (where k is the Boltzmann constant, T is absolute temperature, and lambda is the diffusive jump distance). We have performed preliminary experiments on Na2Si4O9 liquid in order to explore the feasibility of this method up to pressures in excess of 10 GPa. Diffusion couples were prepared from glass starting materials with one half of the sample enriched in O-18 and the other half containing the natural abundance (0.2 wt percent). Diffusion experiments were performed at 1600 to 1800 C at pressures in the range of 2.5 to 10 GPa for times up to 6 minutes using a 1200 ton multianvil apparatus. Oxygen diffusivities were derived from the resulting O-18 concentration profiles, which were analyzed using an ion microprobe. At 1800 C, oxygen diffusivities increase continuously from 1 x 10(exp -10) m(exp -2) s(exp -1) at 2.5 Gpa to 5 x 10(exp -10) m(exp -2) s(exp -1) at 10 GPa. From the Eyring relationship, the viscosity of Na2SiO9 liquid is predicted to decrease by approx. 0.7 log units over this pressure range at 1800 C. These trends, which can be related to changes in the structure of the liquid at high pressure, are in agreement with results of molecular dynamics calculations. These results suggest that it should be possible to estimate the viscosities of silicate liquids up to at least 15 GPa and 2200 C by this method.

Rubie, D. C.↗

A new dipolar potential for numerical simulations of polar fluids on the 4D hypersphere

We present a new method for Monte Carlo or Molecular Dynamics numerical simulations of three-dimensional polar fluids. The simulation cell is defined to be the surface of the northern hemisphere of a four-dimensional (hyper)sphere. The point dipoles are constrained to remain tangent to the sphere and their interactions are derived from the basic laws of electrostatics in this geometry. The dipole-dipole potential has two singularities which correspond to the following boundary conditions: when a dipole leaves the northern hemisphere at some point of the equator, it reappears at the antipodal point bearing the same dipole moment. We derive all the formal expressions needed to obtain the thermodynamic and structural properties of a polar liquid at thermal equilibrium in actual numerical simulation. We notably establish the expression of the static dielectric constant of the fluid as well as the behavior of the pair correlation at large distances. We report and discuss the results of extensive numerical Monte Carlo simulations for two reference states of a fluid of dipolar hard spheres and compare these results with previous methods with a special emphasis on finite size effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exchangeable Liquid Crystalline Elastomers: Enabling Rapid Processing and Enhanced Actuation Stability through On-Demand Deactivation

Exchangeable liquid crystalline elastomers (xLCEs) bearing dynamic covalent bonds are promising candidates for soft actuators due to their unique capability to adjust both network structure and liquid crystalline (LC) alignment after polymerization. While current xLCEs with low exchange temperatures are convenient for processing, they suffer from issues such as creep and loss of LC alignment during repeated thermal actuation. Herein, we present an effective solution using dynamic anhydride chemistry within a thiol-ene-based xLCE. This approach enables a catalyst-free, low-temperature bond exchange of the xLCE after polymerization, allowing for the adjustment of LC orientation under mild conditions. More importantly, it enables on-demand deactivation of the bond exchange via anhydride hydrolysis, effectively eliminating creep and enhancing actuation stability for over 100 cycles. Furthermore, the hydrolysis process results in the formation of carboxylic acid groups, which can be converted into carboxylates via alkali treatment, thereby providing the xLCE with humidity responsiveness. Finally, these findings highlight the use of dynamic anhydride bonds in the fabrication and optimization of xLCEs with enhanced durability and functionality, which is expected to facilitate significant advancements in their applications in soft actuators and robotics.

36 MATERIALS SCIENCE↗

Rapid sample delivery for megahertz serial crystallography at X-ray FELs

Liquid microjets are a common means of delivering protein crystals to the focus of X-ray free-electron lasers (FELs) for serial femtosecond crystallography measurements. The high X-ray intensity in the focus initiates an explosion of the microjet and sample. With the advent of X-ray FELs with megahertz rates, the typical velocities of these jets must be increased significantly in order to replenish the damaged material in time for the subsequent measurement with the next X-ray pulse. Here in this paper, this work reports the results of a megahertz serial diffraction experiment at the FLASH FEL facility using 4.3 nm radiation. The operation of gas-dynamic nozzles that produce liquid microjets with velocities greater than 80 m s -1 was demonstrated. Furthermore, this article provides optical images of X-ray-induced explosions together with Bragg diffraction from protein microcrystals exposed to trains of X-ray pulses repeating at rates of up to 4.5 MHz. Finally, the results indicate the feasibility for megahertz serial crystallography measurements with hard X-rays and give guidance for the design of such experiments.

36 MATERIALS SCIENCE↗

On the Statistical Mechanics of Mass Accommodation at Liquid–Vapor Interfaces

Here we propose a framework for describing the dynamics associated with the adsorption of small molecules to liquid-vapor interfaces using an intermediate resolution between traditional continuum theories that are bereft of molecular detail and molecular dynamics simulations that are replete with them. In particular, we develop an effective single particle equation of motion capable of describing the physical processes that determine thermal and mass accommodation probabilities. The effective equation is parametrized with quantities that vary through space away from the liquid-vapor interface. Of particular importance in describing the early time dynamics is the spatially dependent friction, for which we propose a numerical scheme to evaluate from molecular simulation. Taken together with potentials of mean force computable with importance sampling methods, we illustrate how to compute the mass accommodation coefficient and residence time distribution. Throughout, we highlight the case of ozone adsorption in aqueous solutions and its dependence on electrolyte composition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tuning the Functionalities of Porous Liquids for Emergent Gas-Capture Properties

Type 3 Porous Liquids (PLs) are a class of materials with the potential to revolutionize gas capture, storage, and utilization. These PLs are formed by suspending sorbent nanoparticles (e.g., metal−organic frameworks) in sterically excluded solvents, creating permanent porosity for gas capture in a processable, low-viscosity phase. Herein, a computational study revealed sorbent surface functionalization strategies to enhance CO 2 sorption, and the molecular structural signatures underpinning the enhancements in gas uptake. PLs composed of a ZIF-8 surface functionalized with 3-amino-1,2,4-triazole (Atz) in glyceryl triacetate were targeted for emergent CO 2 capture, exceeding that of the unfunctionalized ZIF-8 PL. ZIF-8 was surface functionalized with Atz at various surface coverage fractions (f), and classical molecular dynamics simulations predicted an increase in CO 2 sorption capacity with increasing f, up to f = 0.75. Additionally, detailed structural analyses revealed that solvent orientational order, derived from the solvent triplet-angle distribution, can identify the gas-capture potential of a PL without requiring computationally expensive direct modeling of the CO 2 sorption. Combined with experimental validation, initial computational screening of PL compositions promises to accelerate the discovery of PL compositions for novel gas separation materials platforms.

gas capture↗

Fragile-to-strong transition in liquid As 2 S 3 under pressure: The effect of melt metallization

The well-known classification of glass-forming melts into fragile and strong liquids has several notable exceptions, including water, silica, and certain phase-change materials (PCMs). These exceptional fluid systems exhibit a fragile-to-strong transition (FST) upon cooling: a transformation from a high-temperature liquid with fast atomic dynamics, low viscosity, and low flow activation energy, to a viscous supercooled melt with high energy barriers near the glass transition temperature T g . This behavior is critically important for non-volatile memories, photonic tensor cores, reconfigurable metamaterials, and other devices, that use PCMs, enabling nanosecond-scale crystallization in the fragile regime and long data retention in the strong regime near or below T g . A significant structural transformation is expected between these two viscosity regimes, along with a semiconductor-metal (SC-M) transition upon heating, driven by high internal pressure and associated density increase. By applying high external pressure to the canonical low-conducting chalcogenide melt As 2 S 3 , we observed both the FST and the SC-M transition, occurring simultaneously within the same domain of the P, T−phase space. These findings suggest that the FST is not limited to a few exceptional liquids but is a common phenomenon, at least in systems that exhibit melt metallization within specific regions of their P, T−phase diagrams.

first-principles molecular dynamics↗

Liquid-liquid-solid transition in viscoelastic liquids

Liquid-liquid-solid transitions (LLST) are known to occur in confined liquids, exist in supercooled liquids and emerge in liquids driven from equilibrium. Molecular dynamics (MD) simulations claim many successes in forecasting the phenomena. The transitions are also studied in the framework of thermodynamics based methods and minimalistic models. In here, the proposed approach is derived in the framework of continuum and includes spatial and temporal dynamic heterogeneities; the approach is meant to capture the material behavior at small scales. We conjecture that the liquid-like and solid-like behaviors are dissimilar enough for the two to be governed by different constitutive relations. In this way, we gain additional degree of freedom, which is found essential when predicting the transitional phenomena. As a result, we derive the LLST criteria for liquids in equilibrium, during steady flow and at transient conditions. Lastly, we forecast short-lived LLSTs in human blood during cardiac cycle.

SELF-ASSEMBLY↗

Capacitance enhancement via electrode patterning

The necessity of increasing the energy density in electric double layer capacitors to meet current demand is fueling fundamental and applied research alike. We report here molecular dynamics simulation results for aqueous electrolytes near model electrodes. Particular focus is on the effect of electrode patterning on the structure of interfacial electrolytes, and on the potential drop between the solid electrodes and the bulk electrolytes. The latter is estimated by numerically integrating the Poisson equation using the charge densities due to water and ions accumulated near the interface as input. We considered uniform and patterned electrodes, both positively and negatively charged. The uniformly charged electrodes are modeled as graphite. The patterned ones are obtained by removing carbon atoms from the top-most graphene layer, yielding nanoscopic squares and stripes patterns. For simplicity, the patterned electrodes are effectively simulated as insulators (the charge remains localized on the top-most layer of carbon atoms). Our simulations show that the patterns alter the structure of water and the accumulation of ions at the liquid-solid interfaces. Using aqueous NaCl solutions, we found that while the capacitance calculated for three positively charged electrodes did not change much, that calculated for the negatively charged electrodes significantly increased upon patterning. We find that both water structure and orientation, as well as ion accumulation affect the capacitance. As electrode patterning affects differently water structure and ion accumulation, it might be possible to observe ion-specific effects. These results could be useful for advancing our understanding of electric double layer capacitors, capacitive desalination processes, as well as of fundamental interfacial electrolytes properties.

36 MATERIALS SCIENCE↗

Molecular hydrogen solvated in water – A computational study

The aqueous hydrogen molecule is studied with molecular dynamics simulations at ambient temperature and pressure conditions, using a newly developed flexible and polarizable H{sub 2} molecule model. The design and implementation of this model, compatible with an existing flexible and polarizable force field for water, is presented in detail. The structure of the hydration layer suggests that first-shell water molecules accommodate the H{sub 2} molecule without major structural distortions and two-dimensional, radial-angular distribution functions indicate that as opposed to strictly tangential, the orientation of these water molecules is such that the solute is solvated with one of the free electron pairs of H{sub 2}O. The calculated self-diffusion coefficient of H{sub 2}(aq) agrees very well with experimental results and the time dependence of mean square displacement suggests the presence of caging on a time scale corresponding to hydrogen bond network vibrations in liquid water. Orientational correlation function of H{sub 2} experiences an extremely short-scale decay, making the H{sub 2}–H{sub 2}O interaction potential essentially isotropic by virtue of rotational averaging. The inclusion of explicit polarizability in the model allows for the calculation of Raman spectra that agree very well with available experimental data on H{sub 2}(aq) under differing pressure conditions, including accurate reproduction of the experimentally noted trends with solute pressure or concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Mobile Interfaces: Liquids as a Perfect Structural Material for Multifunctional, Antifouling Surfaces

Life creates some of its most robust, extreme surface materials not from solids but from liquids: a purely liquid interface, stabilized by underlying nanotexture, makes carnivorous plant leaves ultraslippery, the eye optically perfect and dirt-resistant, our knees lubricated and pressure-tolerant, and insect feet reversibly adhesive and shape-adaptive. Novel liquid surfaces based on this idea have recently been shown to display unprecedented omniphobic, self-healing, anti-ice, antifouling, optical, and adaptive properties. In this Perspective, we present a framework and a path forward for developing and designing such liquid surfaces into sophisticated, versatile multifunctional materials. Drawing on concepts from solid materials design and fluid dynamics, we outline how the continuous dynamics, responsiveness, and multiscale patternability of a liquid surface layer can be harnessed to create a wide range of unique, active interfacial functions able to operate in harsh, changing environments not achievable with static solids. We discuss how, in partnership with the underlying substrate, the liquid surface can be programmed to adaptively and reversibly reconfigure from a defect-free, molecularly smooth, transparent interface through a range of finely tuned liquid topographies in response to environmental stimuli. With nearly unlimited design possibilities and unmatched interfacial properties, liquid materials as long-term stable interfaces yet in their fully liquid state may potentially transform surface design everywhere from medicine to architecture to energy infrastructure.

Grinthal, A↗

Characterization of Slosh Damping for Ortho-Grid and Iso-Grid Internal Tank Structures

Grid stiffened tank structures such as Ortho-Grid and Iso-Grid are widely used in cryogenic tanks for providing stiffening to the tank while reducing mass, compared to tank walls of constant cross-section. If the structure is internal to the tank, it will positively affect the fluid dynamic behavior of the liquid propellant, in regard to fluid slosh damping. As NASA and commercial companies endeavor to explore the solar system, vehicles will by necessity become more mass efficient, and design margin will be reduced where possible. Therefore, if the damping characteristics of the Ortho-Grid and Iso-Grid structure is understood, their positive damping effect can be taken into account in the systems design process. Historically, damping by internal structures has been characterized by rules of thumb and for Ortho-Grid, empirical design tools intended for slosh baffles of much larger cross-section have been used. There is little or no information available to characterize the slosh behavior of Iso-Grid internal structure. Therefore, to take advantage of these structures for their positive damping effects, there is much need for obtaining additional data and tools to characterize them. Recently, the NASA Marshall Space Flight Center conducted both sub-scale testing and computational fluid dynamics (CFD) simulations of slosh damping for Ortho-Grid and Iso-Grid tanks for cylindrical tanks containing water. Enhanced grid meshing techniques were applied to the geometrically detailed and complex Ortho-Grid and Iso-Grid structures. The Loci-STREAM CFD program with the Volume of Fluid Method module for tracking and locating the water-air fluid interface was used to conduct the simulations. The CFD simulations were validated with the test data and new empirical models for predicting damping and frequency of Ortho-Grid and Iso-Grid structures were generated.

Douglas G. Westra↗