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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Probing the Self-Assembly dynamics of cellulose nanocrystals by X-ray photon correlation spectroscopy

Charge-stabilized colloidal cellulose nanocrystals (CNCs) can self-assemble into higher-ordered chiral nematic structures by varying the volume fraction. The assembly process exhibits distinct dynamics during the isotropic to liquid crystal phase transition, which can be elucidated using X-ray photon correlation spectroscopy (XPCS). Anionic CNCs were dispersed in propylene glycol (PG) and water spanning a range of volume fractions, encompassing several phase transitions. Coupled with traditional characterization techniques, XPCS was conducted to monitor the dynamic evolution of the different phases. Additionally, simulated XPCS results were obtained using colloidal rods and compared with the experimental data, offering additional insights into the dynamic behavior of the system. The results indicate that the particle dynamics of CNCs undergo a stepped decay in three stages during the self-assembly process in PG, coinciding with the observed phases. The phase transitions are associated with a total drop of Brownian diffusion rates by four orders of magnitude, a decrease of more than a thousand times slower than expected in an ideal system of repulsive Brownian rods. Given the similarity in the phase behaviors in CNCs dispersed in PG and in water, we hypothesize that these dynamic behaviors can be extrapolated to other polar solvent environments. Importantly, these findings represent the direct measurement of CNC dynamics using XPCS, underscoring the feasibility of directly assessing the dynamic behavior of other rod-like colloidal suspensions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Realistic fuel spray modeling for gasoline direct injection engine applications

Fuel spray modeling plays a critical role during modern gasoline direct injection (GDI) engine development due to fuel injection’s dominant impact on engine performance and emissions as well as the complex physical processes involved. In engineering three-dimensional (3D) computational fluid dynamics (CFD) simulations, the liquid-phase fuel atomization, evaporation, and mixing are usually modeled with the discrete droplet model (DDM) adopting a Lagrangian approach for multiphase CFD simulations. To this end, general practices heavily depend on the reduced order characterization of the injector nozzle flow. However, such simplified injector modeling may lead to insufficient representations of the complex spray dynamics. To tackle this problem, this study proposes a novel workflow to numerically evaluate GDI sub-cooled and flash-boiling sprays under engine-relevant conditions using a side-mounted GDI injector together with real gasoline fuel properties. The workflow introduces a one-way coupling (OWC) method leveraging high-fidelity nozzle flow simulations to provide realistic boundary conditions to the Lagrangian injector model. The proposed workflow was first verified in a constant volume chamber (CVC) environment and then implemented in a practical GDI engine setup to study spray morphology, fuel-air mixing, and wall-wetting propensity. In addition, detailed comparison was performed between the OWC method and the conventional rate of injection (ROI) routine. In conclusion, quantitative analysis of spray characteristics was conducted to highlight possible source of discrepancies of the conventional ROI method.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2025

Data-driven picosecond X-ray imaging for quantitative plasma-induced shock characterization

Imaging dynamic events, especially shockwave behavior, is key to advancing high-energy-density (HED) research. Recent advances in fourth- and fifth-generation X-ray light sources allow for high-resolution imaging of fast phenomena, but limited beam time necessitates maximizing data acquisition. We present a benchtop-scale pulsed plasma device submerged in liquid heptane, capable of generating dynamic events at rates exceeding 10 Hz, supporting the field’s data-driven goals by producing large, high-quality imaging datasets. Using X-ray phase contrast imaging (XPCI) at the Advanced Photon Source, we imaged weak shockwaves (Mach ~ 1.2) in heptane interacting with plasma-induced cavitation bubbles, causing deviation from Rankine-Hugoniot behavior; to our knowledge, this represents the first direct imaging of such interaction. Our quantitative analysis offers insight into weak shock phenomena and energy-focusing applications in pulsed plasmas. These results highlight the potential for large datasets to advance dynamic HED research at current light source facilities, and have implications for fields such as inertial confinement fusion, plasma-enhanced chemical processing, and biomedical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Gelation and Re-entrance in Mixtures of Soft Colloids and Linear Polymers of Equal Size

Liquid mixtures composed of colloidal particles and much smaller non-adsorbing linear homopolymers can undergo a gelation transition due to polymer-mediated depletion forces. We now show that the addition of linear polymers to suspensions of soft colloids having the same hydrodynamic size yields a liquid-togel-to-re-entrant liquid transition. In particular, the dynamic state diagram of 1,4-polybutadiene star-linear polymer mixtures was determined with the help of linear viscoelastic and small-angle Xray scattering experiments. While keeping the star polymers below their nominal overlap concentration, a gel was formed upon increasing the linear polymer content. Further addition of linear chains yielded a re-entrant liquid. This unexpected behavior was rationalized by the interplay of three possible phenomena: (i) depletion interactions, driven by the size disparity between the stars and the polymer length scale which is the mesh size of its entanglement network; (ii) colloidal deswelling due to the increased osmotic pressure exerted onto the stars; and (iii) a concomitant progressive suppression of the depletion efficiency on increasing the polymer concentration due to reduced mesh size, hence a smaller range of attraction. Our results unveil an exciting new way to tailor the flow of soft colloids and highlight a largely unexplored path to engineer soft colloidal mixtures.

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

Pressure Drop with Moving Contact Lines and Dynamic Contact Angles in a Hydrophobic Round Minichannel: Visualization via Synchrotron X-ray Imaging and Verification of Experimental Correlations

In hydrophobic mini- and microchannels, slug flow with moving contact lines is typically generated under various two-phase flow conditions. There is a significant pressure drop in this flow pattern with moving contact lines, which is closely related to the dynamic contact angles. Researchers have investigated dynamic contact angles experimentally for decades, but due to the limitations of visualization techniques, these experiments have typically been conducted in low Weber number regions (We < 10 –3 ). In this study, we clearly visualized the dynamic contact angles of a liquid slug in high Weber number regions (10 –3 < We <1) via synchrotron X-ray imaging with high temporal (~1000 fps) and spatial (~ 2 μm/pixel) resolutions. We precisely measured the pressure drop with moving contact lines in a hydrophobic minichannel (inner diameter = 1.018 mm). On the basis of our experimental data, we verified previous correlations for dynamic contact angles and explored the relationship between pressure drop with moving contact lines and dynamic contact angles.

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

Ion Correlations and Partial Ionicities in the Lamellar Phases of Block Copolymeric Ionic Liquids

Recently, significant interest has arisen on the impact of dynamical ion correlations on the conductivity and transport properties of polymeric electrolyte materials. It has been hypothesized that confining ion motion to narrow channels may reduce such ion correlations and enhance the resulting ionic conductivity. Motivated by such considerations, in this study we used a multiscale simulation framework to study the dynamical ion correlations in the microphase-separated lamella phase of block copolymeric ionic liquids and compare with the corresponding results for homopolymeric systems. We probed the influence of ion correlations through the partial ionicity, Δ, which quantifies the ratio of true conductivity to the ideal, Nernst–Einstein conductivity for the anion-related contributions. Consistent with our original hypothesis, our results demonstrate that the partial ionicity relating to the mobile anions is much larger in the lamella phases of block copolymers compared to that in homopolymers. Analysis of the distinct conductivity contributions demonstrates that such results arise as a result of an intricate compensation among the nonideal dynamical correlations relating to anions in lamella phases. Overall, our results suggest that self-assembled phases of block copolymers may provide an avenue to tune the dynamical ion correlations in polymer electrolyte systems.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Pressure fields produced by single-bubble collapse near a corner

Damage produced by repeated bubble collapse to neighboring rigid objects in hydraulic systems is an important consequence of cavitation. Although bubble collapse near a single wall has received significant attention in the past, few studies exist on the dynamics of bubbles collapsing near a corner, i.e., two flat rigid surfaces intersecting at a right angle. In this work we quantify the pressure fields produced by a single bubble collapsing near two perpendicular rigid walls. Using a high-order accurate shock- and interface-capturing method to solve the three-dimensional compressible Navier-Stokes equations for gas and liquid flows, we simulate the dynamics of a single bubble collapsing at different initial stand-off distances from the two walls. In contrast to a bubble collapsing near a single wall, the collapse of bubbles within a critical stand-off distance is not symmetric about the bisecting plane due to the interaction between the bubble and the second wall. The second wall affects the pressure produced during the collapse in the following ways: (i) For bubbles initially located sufficiently close to both walls, the reentrant jet produced during collapse no longer points in the direction normal to the closest wall but at an angle toward the corner, (ii) the part of the emitted shock with the highest amplitude propagates in line with the jet, and (iii) the bubble migrates in that same direction during its collapse with a dependence on the stand-off distance, consistent with predictions made using Kelvin impulse. The location of maximum pressure along the walls is measured for the different initial stand-off distances. Using acoustic arguments, we find a semiempirical relationship to predict the initial stand-off distances for which the maximum pressure occurs in the corner. Here we find that when the bubble is sufficiently close to equidistant from each boundary, the maximum pressure is observed in the corner due to the water-hammer and implosion shocks reflecting off the boundaries and intersecting in the corner. We also show that when bubbles are initially attached to either wall the wall pressure produced can be significantly increased compared to bubbles detached from either wall.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Phase diagram of the TIP4P/Ice water model by enhanced sampling simulations

In this work, we studied the phase diagram for the TIP4P/Ice water model using enhanced sampling molecular dynamics simulations. Our approach is based on the calculation of ice–liquid free energy differences from biased coexistence simulations that reversibly sample the melting and growth of layers of ice. We computed a total of 19 melting points for five different ice polymorphs, which are in excellent agreement with the melting lines obtained from the integration of the Clausius–Clapeyron equation. For proton-ordered and fully proton-disordered ice phases, the results are in very good agreement with previous calculations based on thermodynamic integration. For the partially proton-disordered ice III, we find a large increase in stability that is in line with previous observations using direct coexistence simulations for the TIP4P/2005 model. This issue highlights the robustness of the approach employed here for ice polymorphs with diverse degrees of proton disorder. Our approach is general and can be applied to the calculation of other complex phase diagrams.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Dynamic condensation model of rolling droplets for high-performance heat transfer

Condensation is a critical process during vapor-liquid phase change in relation to heat transfer. To achieve a high heat transfer coefficient, the classical model for dropwise condensation requires a low contact angle and low contact angle hysteresis, failing to align with experimental observations on a hydrophobic and slippery quasi-liquid surface (QLS). Here, we report a dynamic condensation model that incorporates high-frequency condensate removal by emphasizing the role of timescale during droplet growth and shedding. Our model agrees well with the experimental result that a surface with high contact angle and low contact angle hysteresis promotes condensation, particularly during rolling-propelled condensate removal. Particle image velocimetry reveals that rolling droplets on a hydrophobic QLS exhibit 4-fold higher shedding speeds than the sliding droplets on a hydrophilic QLS, leading to significant heat transfer enhancement. This work deepens our theoretical understanding of condensation heat transfer and provides advanced physics-informed design rationales for water and energy systems.

42 ENGINEERING↗
DOE OSTI2021

Structures of glass-forming liquids by x-ray scattering: Glycerol, xylitol, and D-sorbitol

In this work, synchrotron X-ray scattering has been used to investigate three liquid polyalcohols of different sizes (glycerol, xylitol and D-sorbitol) from above the glass transition temperatures $T_g$ to below. We focus on two structural orders: the association of the polar OH groups by hydrogen bonds (HBs) and the packing of the non-polar hydrocarbon groups. We find that the two structural orders evolve very differently reflecting the different natures of bonding. Upon cooling from 400 K, the O$\cdots$O correlation at 2.8 Å increases significantly in all three systems, indicating more HBs, until kinetic arrest at $T_g$; the increase is well described by an equilibrium between bonded and non-bonded OH with ΔH = 9.1 kJ/mol and ΔS = 13.4 J/mol/K. When heated above $T_g$, glycerol loses the fewest HBs per OH, followed by xylitol and by D-sorbitol, in the same order the number of OH groups per molecule increases (3, 5, and 6). The pair correlation functions of all three liquids show exponentially damped density modulations of wavelength 4.5 Å, which are associated with the main scattering peak and with the intermolecular C$\cdots$C correlation. In this respect, glycerol is the most ordered with the most persistent density ripples, followed by D-sorbitol and by xylitol. Heating above $T_g$ causes faster damping of the density ripples with the rate of change being the slowest in xylitol, followed by glycerol and by D-sorbitol. Given the different dynamic fragility of the three liquids (glycerol being the strongest and D-sorbitol the most fragile), we relate our results to the current theories of the structural origin for the difference. We find that the fragility difference is better understood on the basis of the thermal stability of HB clusters than that of the structure associated with the main scattering peak.

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

Free-Surface Liquid Lithium Flow Modeling and Stability Analysis for Fusion Applications

Liquid metal plasma facing components are considered an attractive design choice for fusion devices including pilot plants. Virtual prototyping of such devices includes modeling of free-surface flow of the electrically conductive liquid, which requires computational fluid dynamics (CFD) and magnetohydrodynamics (MHD) simulations. Numerical tools capable of simulating flows and heat transfer in the free-surface MHD flow were developed at PPPL based on the customized ANSYS CFX. Here, MHD is introduced using a magnetic vector potential approach. Free-surface flow capabilities are available in the code and were tested. Special stabilization procedures were derived and applied to improve convergence of the momentum equations with the source terms due to the Lorentz force and surface tension. Important characteristics of the fusion-relevant liquid metal flow is free surface smoothness and stability. Heat flux from the plasma impacts the liquid surface at a very acute angle, so any change of the free surface from axisymmetry can dramatically increase the local heat flux density and thus create excessive evaporation of liquid lithium into the plasma, which is detrimental to operations. Stability analysis of the liquid metal film flow was performed to determine applicable flow regimes. Thin film flow along horizontal wall is considered. Effects of gravity, magnetic field, and surface tension are included in the analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electrochemical, thermodynamic, and physical properties of tetradecyltrihexylphosphonium ([P 6,6, 6,14 ] + ) and methyl-propyl piperidinium containing ionic liquids and their propylene carbonate solutions

Here, the viscosity, conductivity, electrochemical window and related thermodynamic properties such as excess volume and dynamic viscosity deviation of two phosphonium ionic liquids were measured and calculated for “as-supplied” and dried neat liquids and also v/v mixtures of 99/1, 95/5, 90/10, 75/25, and 50/50 of the ionic liquids with propylene carbonate (PC). Tetradecyltrihexylphosphonium dicyanamide ([P 6,6,6,14 ] + dicyanamide), Cyphos 105, and tetradecyltrihexylphosphonium bis(trifluororomethanesulfonyl)imide (called bistriflimide) ([P 6,6,6,14 ] + bistriflimide), Cyphos 109, ionic liquids were studied. Generally speaking, there were slight differences in the measured properties of the wet and dried IL solutions which were reflected in differences if the calculated properties. The measured and calculated properties were compared with those collected and derived from a piperidinium-based ionic liquid, methyl-propyl piperidinium bistriflimide. Arrhenius plots for both the viscosity and conductivity were linear, with the linearity of the Litovitz plots being slightly higher. The viscosity and conductivity of the phosphonium solutions yielded Walden Plots that differed from Walden plots of piperidinium solutions and other ionic liquid solutions. For the phosphonium based ionic liquids, the ionicity was found to increase with increasing dilution in propylene carbonate. The electrochemical windows of the ionic liquids were determined as a function of the concentration of the ionic liquid present in the solution. The size of the window generally initially decreased and then increased with the addition of PC for the [P 6,6,6,14 ] + dicyanamide ionic liquid increasing from 3.4 V in the 95/5 v/v (volume Cyphos 105) to volume of PC solvent) to 3.7 in the 50/50 dried IL case. The electrochemical windows of the dried liquids were slightly larger 3.5 V for the 95/5 v/v and 4.2 V for the 50/50 v/v Cyphos 105/PC solutions. The electrochemical window was larger and the conductivity was higher for the piperidinium IL solutions, but since the ionicity increases with dilution, the phosphonium based IL solutions may prove favorable for processes such as electrochemical reduction.

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

Rotational Dynamics and Transition Mechanisms of Surface-Adsorbed Proteins

Assembly of biomolecules at solid-water interfaces requires molecules to traverse complex orientation-dependent energy landscapes through processes that are poorly understood, largely due to the dearth of in-situ single molecule measurements and statistical analyses of the rotational dynamics that define directional selection. Emerging capabilities in high-speed atomic force microscopy and machine learning have allowed us to directly determine the orientational energy landscape and observe and quantify the rotational dynamics for protein nanorods on the surface of muscovite mica under a variety of conditions. Comparisons with kinetic Monte Carlo simulations show that the transition rates between adjacent orientation-specific energetic minima can largely be understood through traditional models of in-plane Brownian rotation across a biased energy landscape, with resulting transition rates that are exponential in the energy-barriers between states. However, transitions between more distant angular states are decoupled from barrier height, with jump-size distributions showing a power-law decay that is characteristic of a non-classical Levy-flight random walk, indicating that large jumps are enabled by alternative modes of motion via activated states. The findings provide new insights into the dynamics of biomolecules at solid-liquid interfaces that lead to self-assembly, epitaxial matching and other orientationally anisotropic outcomes and define a general procedure for exploring such dynamics with implications for hybrid biomolecular-inorganic materials design.

Orientational energy landscapes, Rotational dynami↗
DOE OSTI2022

Plasma dynamics, instabilities and OH generation in a pulsed atmospheric pressure plasma with liquid cathode: a diagnostic study

Abstract While plasma–liquid interactions have been an important focus in the plasma research community, the impact of the strong coupling between plasma and liquid on plasma properties and processes remains not fully understood. In this work, we report on the impact of the applied voltage, pulse width and liquid conductivity on the plasma morphology and the OH generation for a positive pulsed DC atmospheric pressure plasma jet with He–0.1% H 2 O mixture interacting with a liquid cathode. We adopted diagnostic techniques of fast imaging, 2D laser induced fluorescence of OH and Thomson scattering spectroscopy. We show that plasma instabilities and enhanced evaporation occur and have a significant impact on the OH generation. At elevated plasma energies, it is found that the plasma contracts due to a thermal instability through Ohmic heating and the contraction coincides with a depletion in the OH density in the core due to electron impact dissociation. For lower plasma energies, the instability is suppressed/delayed by the equivalent series resistor of the liquid electrode. An estimation of the energy flux from the plasma to the liquid shows that the energy flux of the ions released into the liquid by positive ion hydration is dominant, and significantly larger than the energy needed to evaporate sufficient amount of water to account for the measured H 2 O concentration increase near the plasma–liquid interface.

Physics↗
DOE OSTI2024

Molecular dynamics simulations of uranyl and plutonyl cations in a task-specific ionic liquid

Ionic liquids (ILs) are a unique class of solvents with potential applications in advanced separation technologies relevant to the nuclear industry. ILs are salts with low melting points and a wide range of tunable physical properties, such as viscosity, hydrophobiciy, conductivity, and liquidus range. ILs have negligible vapor pressure, are often non-flammable, and can have high thermal stability and a wide electrochemical window, making them attractive for use in separations processes relevant to the nuclear industry. Metal salts generally have a low solubility in ILs; however, by incorporating new functional groups into the IL cation or anion that promote complexation with the metal, the solubility can be greatly increased. One such task-specific ionic liquid (TSIL) is 1-carboxy-N, N, N-trimethylglycine bis(trifluoromethylsulfonyl)imide ([Hbet][Tf 2 N]). Water, which is detrimental for electrochemical separations, is a common impurity in ILs and can coordinate with actinyl cations, particularly in ILs containing only weakly coordinating components. Understanding the behavior of actinides in TSIL/water mixtures on a molecular level is vital for designing improved separations processes. Classical molecular dynamics simulations of uranyl(VI) and plutonyl(VI) in 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf 2 N]) with deprotonated Hbet (betaine) and water have been performed to understand the coordination and dynamics of the actinyl cations. We find that betaine is a much stronger ligand than water and prefers to coordinate the metal in a bidentate manner. Potential of mean force simulations yield a relative free energy for betaine coordination of approximately -120 to -90 kJ/mol in mixtures with water. As the amount of betaine coordinated to the actinide increases, the diffusion coefficient of the actinyl cation decreases. Moreover, the betaine ligand is able to bridge between two metal centers, resulting in dimeric complexes with actinide–actinide distances of ~5 Å. Potential of mean force simulations show that these structures are stable, with relative free energies of up to -40 kJ/mol. The crystal structure for [(UO 2 ) 2 (bet) 6 (H 2 O) 2 ][Tf 2 N] 4 shows that the betaine bridges between two uranium atoms to form dimeric complexes similar to those found in our simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗