Liquid-vapor Interface Research Articles

Soft wetting: an analytical model for pillar topography- and softness-dependent droplet depinning force.

The extent to which a droplet pins on a textured substrate is determined by the dynamics of the contact line and the liquid-vapor interface. However, the synergistic contribution of contact line sliding and interface distortion to the droplet depinning force remains unknown. More strikingly, current models fail to predict the depinning force per unit length of droplets on soft pillar arrays. Therefore, we fabricate soft pillar arrays with varying geometrical dimensions and mechanical properties and measure the depinning forces per unit length by allowing droplets to evaporate on such substrates. We then analyze the decrease in excess Gibbs free energy of the apparent droplet caused by the detachment of the droplet boundary from the previously pinned pillars. In contrast to prior notions, based on the measured decreases in excess Gibbs free energy, we find that the coefficient, that governs the ratio of interface distortion's contribution to the depinning force to that of the sliding contact line, increases with a decrease in pillar packing density. By considering the combined contribution from contact line sliding, liquid-vapor interface distortion, and pillar deflection, we introduce an analytical model to predict the droplet depinning force per unit length and corroborate the model using experimental data reported in this and prior studies.

Surface accumulation and acid-base equilibrium of phenol at the liquid-vapor interface.

We have investigated the surfactant properties of phenol in aqueous solution as a function of pH and bulk concentration using liquid-jet photoelectron spectroscopy (LJ-PES) and surface tension measurements. The emphasis of this work is on the determination of the Gibbs free energy of adsorption and surface excess of phenol and its conjugate base phenolate at the bulk pKa (9.99), which can be determined for each species using photoelectron spectroscopy. These values are compared to those obtained in measurements well below and well above the pKa, where pure phenol or phenolate, respectively, are the dominant species, and where the Gibbs free energy of adsorption determined from surface tension and LJ-PES data are in excellent agreement. At the bulk pKa the surface-sensitive LJ-PES measurements show a deviation of the expected phenol/phenolate ratio in favor of phenol, i.e., an apparent upward shift of the at the surface. In addition, the Gibbs free energies of adsorption determined by LJ-PES at the bulk pKa for phenol and phenolate deviate from those observed for the pure solutions. We discuss these observations in view of the different surface propensity of phenol and phenolate as well as potential cooperative interactions between them in the near-surface region.

Stability and structure of the aqueous LiTFSI-LiCl interface.

It has recently been demonstrated that aqueous lithium bis(trifluoromethanesulfonyl)imide (LiTFSI) and lithium chloride (LiCl) solutions can form stable liquid-liquid biphasic systems when both electrolyte phases have sufficiently high concentrations. In this work, we combine molecular dynamics simulations and experimental analysis to investigate what drives the formation of the interface and how the interfacial molecular structure correlates with its thermodynamic stability. We observe that at the liquid-vapour interface, TFSI- anions exhibit surfactant-like properties, leading to a reduction in surface tension and an increase in interfacial thickness. In contrast, the interfacial stability of the LiTFSI-LiCl biphasic systems increases with the concentration of both salts, as evidenced by the increasing surface tension and decreasing interfacial thickness. The opposing effects that the ionic concentration has on the thermodynamic stability of the different interfaces are linked to the anions' interfacial adsorption/desorption, which in turn affects the number and strength of water-water hydrogen bonds, the interfacial molecular structure and the diffusion of cations across the interface. Finally, calculations and experiments indicate that the liquid-liquid separation is driven primarily by the concentration of LiCl, and is the result of a 'salting out' effect.

Research on the Unsteady Flow and Vortex Characteristics of Cavitation at the Tongue in Centrifugal Pump

In order to study the unsteady flow and vortex characteristics of tongue cavitation, numerical calculation is carried out for the whole flow channel given different conditions. Then, the calculation results are verified experimentally. The results show that after the occurrence of tongue cavitation in the centrifugal pump, it deteriorates with the decrease of NPSHa. However, when NPSHa is reduced to 3.78 m, it does not change significantly anymore with the decrease of NPSHa. The extrusion of fluid by the vapor at the tongue promotes the formation of the separation vortex, and the re-jet flow caused by the separation vortex leads to vapor shedding. The frequency of cavitation shedding is consistent with the frequency of vortex shedding. In the vorticity transport equation, the relative vortex stretching term and the relative vortex dilatation term dominate the vapor shedding by controlling the change in vorticity. The baroclinic torque term mainly affects the change of vorticity at the vapor-liquid interface, but to a much lesser extent than the first two terms.

Energetic description of the liquid-vapor interface of water with organic coating molecules

In situ measurements revealed the presence of a myriad of organic chemical compounds in cloud droplets. Among them, a significant fraction is composed of polar or amphiphilic compounds and these...

Quantitative characterization of the pseudo-boiling contribution to supercritical heat transfer

This paper explores the supercritical heat transfer mechanism by characterizing the boiling contribution ratio qb/q, where qb is the boiling heat flux and q is the applied heat flux. Experiments are performed using nickel–chromium wire in 15 °C liquid carbon dioxide at 5.2, 7.6, 9.0, and 11.0 MPa. The evaporation heat flux qe is the amount of heat used for vapor generation, while qb is the heat transfer in the bulk liquid due to the disturbance of the flow/temperature field by vapor–liquid interface motion. A data processing procedure is developed to measure qb/q from the captured images. Similar trends appear for both supercritical pseudo-boiling and subcritical boiling. The evaporation-like regime at supercritical pressures reaches qb/q = 0.21–0.43, while the film boiling (evaporation) regime achieves qb/q = 0.08. In the supercritical-boiling-like regime, qb/q increases sharply from 0.19 to 0.65, whereas in the subcritical-nucleate-boiling regime, qb/q maintains a value of 0.30 followed by a rapid rise to 0.68 under a vigorous bubble merging and departing mechanism. At both subcritical and supercritical pressures, the heat transfer deteriorates in the evaporation regime, but is significantly enhanced by phase-change-induced flow/temperature field perturbations. The boiling curves differ in the two pressure domains. At supercritical pressures, natural convection transitions smoothly to the evaporation-like regime, then to the boiling-like regime. At subcritical pressures, a steep transition from natural convection to nucleate boiling occurs, and then, film boiling is induced through the action of surface tension. The above findings complete the inverse boiling curves in the two pressure domains.

Subcooled flow boiling in ultrahigh-aspect-ratio microchannels for high heat flux cooling

In present work, the numerical simulation and experimental investigations were performed on the subcooled flow boiling behavior in ultrahigh-aspect-ratio microchannel heat sink (AR = 25). As the heat flux is low, the unique discotic bubble expands the condensation area while ensuring an adequate evaporation area in the channel. When the heat flux increases, the violently fluctuating liquid-vapor interface forms a throat shape that effectively reduces the cross-section of vapor flow. Furthermore, the flow pattern of annular flow at the channel bottom and nucleate bubble at the channel top realizes a combined heat transfer mechanism (nucleate boiling and liquid film evaporation occurs simultaneously). Moreover, the phenomenon of “droplet inside the bubble” caused by bubble coalescence also significantly enhances the thermal process. The fluctuation of local pressure drop in UARM acts gently, which effectively inhibits the boiling instability, and its average value can be reduced by 68.9%. The maximum CHF reaches 1258.1 W/cm2 under a mass flux of 410.2 kg/(m2·s) based on the experimental data, indicating a tremendous potential for the efficient heat dissipation.

Stability of the Liquid-Vapor Interface under the Combined Influence of Normal Vibrations and an Electric Field

Stability of the Liquid-Vapor Interface under the Combined Influence of Normal Vibrations and an Electric Field

Constrained Vapor Bubble Experiment (CVB) in the Light Microscopy Module (LMM)

Abstract This short article describes the major findings from the CVB experiment performed in the LMM on the International Space Station from 2010–2012. CVB was the first experiment to run in the new facility and focused on understanding the heat transfer and fluid mechanics occurring inside a wickless miniature heat pipe. The LMM was used to map the location of the vapor-liquid interface inside the device and to measure the film thickness profile on the walls of the device. Several interesting and unexpected phenomena were observed in microgravity including flooding of the heater end with liquid as the heat input increased, explosive nucleation of vapor bubbles at the heater end in the shortest version of the heat pipe tested, condensation on highly superheated surfaces, and the spontaneous formation of rip currents as the device tried to enhance the contact line area available for evaporation of the liquid.

A multiscale analytical-numerical method for the coupled heat and mass transfer in the extended meniscus region considering thin-film evaporation in microchannels

A Computational Fluid Dynamics (CFD) framework is established to investigate the microscale transport phenomena in the evaporating extended meniscus region formed in rectangular microchannels. A wide range of microchannel widths (W = 1 to 250 µm) and wall superheats (ΔT = 0.01 to 5.0 K) are considered. Prior to performing CFD simulations, it is first necessary to employ thin-film evaporation modeling based on an augmented Young–Laplace model and kinetic-theory based model for the evaporating mass flux across curved surfaces to find the exact shape of the liquid-vapor interface along with disjoining and capillary pressure distributions. Two widely-used methods for the determination of boundary conditions (BCs) at the beginning of the thin-film region (i.e., setting δ0″ = 0 and finding δ0′ by the far-field BC or setting δ0′ = 0 and finding δ0″ by the far-field BC), where δ is film thickness, are thoroughly examined and the first approach is recommended because of its better convergence and ease of implementation in the thin-film evaporation model. Then, the two-dimensional domain of the extended meniscus is generated and discretized to simulate the evaporating liquid flow through the UDF programming in ANSYS Fluent. The developed framework is shown to be a simple yet powerful and practical method for accurately predicting the evaporation mass and heat fluxes from extended menisci in microchannels. The CFD simulation results indicate that the previous one-dimensional models which represent the state of the art in the literature in the field are not able to predict the rate of evaporative heat transfer from the extended meniscus in microchannels with an acceptable degree of accuracy. Based on the CFD simulations results, a multiple regression analysis is employed to establish several simple thermal resistance correlations. The correlations can be straightforwardly integrated into macroscale mathematical models of micro and miniature heat pipes for analyzing their thermal characteristics.

Liquid-Vapor Interface of Aqueous Ethylene Glycol Solutions: A Molecular Dynamics Study.

While the organic constituent in an aqueous binary solution enriches its liquid-vapor (l-v) interface, the extent of enrichment can depend nonlinearly on its mole fraction. A microscopic quantification and rationalization of this behavior are crucial to understand the dependence of properties such as surface tension and evaporation rate of the solution on its composition. Extensive all-atom molecular dynamics simulations of aqueous ethylene glycol (EG) solutions show that the composition of the solution at the l-v interface deviates the most from that in the bulk solution at an EG mole fraction of 0.3. The population of EG molecules with their central C-C dihedral in the gauche conformation was found to be higher at the l-v interface than that in the bulk solution to facilitate the orientation of its hydrophobic methyl groups toward the vapor phase. Free energy calculations reveal that in dilute EG solutions, an EG molecule is most stable at the l-v interface. The behavior of vapor pressure in aqueous EG solutions is ideal and follows Raoult's law, while in contrast, the aqueous solution of dimethyl sulfoxide does not. A rationale for the same is provided through the orientational distribution of interfacial water molecules in the respective solutions.

Enhanced interfacial boiling of impacting droplets upon vibratory surfaces

HypothesisDespite the flourishing studies of droplet interfacial boiling, the boiling upon vibratory surfaces, which may cause vigorous liquid–vapor-solid interactions, has rarely been investigated. Enhanced boiling normally can be gained from rapid removal of vapor and disturbance of liquid–vapor interface. We hypothesize that the vibratory surfaces enhance both effects with new intriguing phenomena and thus, attain an enhanced boiling heat transfer. ExperimentsWe experimentally investigated the impacting fluid dynamics and coupled heat transfer patterns of multiple droplets and a single droplet impinging on still and vibratory surfaces of various materials and different wettability. FindingsThe boiling under vibratory surfaces with increased vibration velocity amplitude and enhanced wettability can be enhanced by 80% in heat transfer coefficient and Nusselt number, which is attributed to several reasons: shortened bubble lifespan, thinner and smaller bubbles, and enhanced disturbances in liquid–vapor interfaces. The vibration also delays the Leidenfrost point when the droplet impacts a descending surface, which shows that the droplet impact moment (vibration phase angle) is particularly crucial. The descending surface releases the generated vapor actively and facilitates liquid–solid contact, thereby delaying the Leidenfrost. From fundamentals to application, this article strengthens our understanding of vibrated interfacial boiling in scenarios closer to multiple natural processes and practical industries.

White Light Reflectance Spectroscopy as a novel and noninvasive method for in situ wetting monitoring in membrane distillation

Wetting, in membrane distillation (MD), is a major drawback appearing as a dynamic mechanism consisting of multiple stages. Normally, the liquid-air interface is placed on top of the membrane surface and pores are only occupied with vapor. If liquid starts to penetrate, surface wetting is followed by partial wetting that develops until the fully wetted state occurs, compromising the operation. Usually, only the fully wetted state is assumed, indirectly, as liquid from the feedwater stream, containing soluble salts, passes through the pores and contaminates the permeate, thus increasing its conductivity. In this work, we present White Light Reflectance Spectroscopy (WLRS) as a novel method for in-situ wetting monitoring in MD. In WLRS, light incidents vertically onto the membrane and reflectance is continuously recorded. Thus, it is possible to observe any movement of the liquid-vapor interface in real-time, via changes in the reflectance spectrum, without compromising the structural integrity of the membrane nor the module. Notably, the proposed method is applicable to all standard MD modes, without operating limitations and restrictions. With this noninvasive technique we can predict and observe all wetting stages, not only in the presence of surfactants in the feed stream, but also due to hydraulic pressure increase.

Enhanced electrolytic immersion cooling for thermal crisis mitigation in high-energy–density systems

Motivated by the increasing need for effective cooling solutions in high-energy–density systems, this experimental study presents the two-phase cooling of a superheated SS 316L sample via immersion quenching in saturated deionized (DI) water at atmospheric pressure conditions with and without a DC electric field. We investigate the effect of the applied electric field, electrode polarity, and in-situ oxidation on quenching characteristics such as the cooling profile, vapour layer behaviour, minimum film boiling temperature (Tmin), and heat transfer rate. The cooling curves of samples quenched with the application of an electric field shift towards the left compared to the sample quenched without an electric field. The cathode sample at 200 V exhibits 33 % faster cooling than the bare sample (0 V). Overall, the in-situ oxidised SS 316L cathode sample at 200 V exhibits a 55 % reduction in film boiling duration, and Tmin increased from 268 °C to 322 °C compared to the bare sample (0 V). The visualisation studies highlight that the liquid–vapour interface experiences a series of oscillations followed by temporal collapse due to electrostatic attraction and electrolytic activity. The obtained results show that hydrogen-rich vapour bubbles increase heat transfer performance. The evolution of hydrogen and its adsorption at the sample surface reduces the activation energy for bubble nucleation and improves the bubble density via liquid pumping. These insights open the pathway for employing hydrogen bubbles for handling ultra-high thermal loads in high-energy density systems, and the specific case of a revised design of concentrating solar receiver is considered based on the present findings.

CFD modelling of the non-isobaric evaporation of cryogenic liquids in storage tanks

A new CFD model relevant to non-isobaric cryogen evaporation has been developed. It considers the heat influx, from the surroundings to the liquid and vapour phases, and the heat transfer between the phases. The phase change is modelled considering two contributions: an interfacial energy balance far from the tank wall, and by direct wall-to-liquid conduction near the wall. The model provides well resolved liquid and vapour temperature and velocity profiles, as well as pressure build-up and evaporation rates. It shows that below the vapour-liquid interface, liquid natural convection is dampened by thermal stratification, which is driven by the increase of the interfacial temperature due to increase in pressure. The pressure build-up is mainly governed by the heating of the vapour, particularly at the beginning of the evaporation. The model was used to analyse liquid velocity profiles for two different sets of experimental data for evaporation of liquid nitrogen by optimizing empirical parameters, namely the wall heat partitioning fraction and the overall heat transfer coefficients. The results show that a large fraction of the vapour heat ingress is not transferred to the vapour phase bulk. Instead, it is transported through the walls and transferred to the liquid at the liquid-vapour-wall contact point, driving phase change. Heat conduction from the walls to the interface has a significant effect, even for low heat fluxes, and cannot be neglected.

Comparative thermokinetic study of vegetable oils and deoxygenated derivatives

This study shows the usefulness of thermokinetic analysis in the characterization of vegetable oils and their deoxygenated derivatives (bioparaffin mixtures). The approach allows obtaining information about the complexity and nature of the process rate-controlling steps (inter- or intra-molecular), to establish quality criteria and to estimate the cohesive forces in liquids and fuels in general. The study also reveals the influence of the gaseous products, Ps, on the limit temperature, Tl, at which the liquid–vapor interface coexists or vanishes. The presence of multiple steps and stages, and the energy barrier variability between ∼ 120 and > 250 kJ mol−1 in the starting vegetable oils, demonstrates the importance of the origin and composition of triglyceride and the effect of unsaturation degree (C = C)n in the fatty acid degradation pathway. Besides, the isoconversional analysis of the deoxygenated products (bioparaffins) describes a kinetics controlled by a single-step, with almost constant activation energy barriers ∼ 70 kJ mol−1, related to the inter-molecular cohesive forces, usually E(α) < 100 kJ mol−1, equivalent to the vaporization enthalpy, ΔvH. The results are validated by a complementary thermokinetic analysis of pure solvents (n-decane and n-hexadecane). Finally, the study describes several experimental observations of great relevance and interest to understand and control the nature of thermally activated processes in vegetable oils and biofuels.

Smoothed particle hydrodynamics simulations of the evaporation of suspended liquid droplets

The ordinary evaporation and explosive vaporization of equilibrium, van der Waals, liquid drops subjected to uniform heating at supercritical temperatures are investigated by means of numerical simulations with the aid of a modified version of the DualSPHysics code. The models include the effects of surface tension, thermocapillary forces, mass transfer across the interface, and liquid–vapor interface dynamics by means of a diffuse-interface description. In contrast to previous simulations in this line, a new non-classical source term has been added to the internal energy equation to deal with the vaporization rate through the diffuse interface. This term is related to the diffusion of the latent heat in the interface zone and is, therefore, necessary for a correct physical description of the liquid–vapor interface structure. As the heating temperature increases the drops undergo surface evaporation, nucleation of an interior vapor bubble, nucleation followed by fragmentation of the liquid, and explosive vaporization. Heating at supercritical temperatures brings the drop out of equilibrium and forces it to rapid quenching into either the miscibility gap, where it undergoes surface evaporation by spinodal decomposition, or the metastable region bounded by the binodal and spinodal curves, where it nucleates a vapor bubble. The results also indicate that at comparable heating, drops of lower density experience faster evaporation rates than drops of higher density.

On the mechanics of thermal collapse of a vapor bubble translating in an isothermal subcooled liquid

In this work, we numerically simulate bubble radius and velocity in a thermal collapse of a vapor bubble translating in an isothermal subcooled liquid based on solving the equations of energy and motion simultaneously. The mechanics of the collapse in the thermal regime plays a significant role in the thermal convection from the vapor-liquid interface of a translating bubble. In our mathematical model, the equation of motion considers the added mass force due to phase change and shape change, while the equation of energy includes thermal conduction and convection. The simulations include both spherical and non-spherical bubbles, the effect of initial bubble diameter and degree of subcooling. The simulations have been performed in three modes of heat transfer (conduction, convection and combined conduction-convection) to observe their effect on the bubble dynamics. The numerical results have been compared with the two classes of thermal collapse experiments (Type A and B) and CFD analysis reported in literature. They agree well with the experiments and CFD-analysis. We have seen a remarkable effect of the initial shape of the bubble in terms of aspect ratio on the time scale and modes of the heat transfer and the collapse time. An interesting concern is that the mechanism of the collapse and the modes of heat transfer depends on how the collapse is initiated in the experiments. Type A and B thermal collapses are mainly classified based on the bubble motion and the modes of thermal transport from the interface.

Enhanced boiling heat transfer on three-dimensional hybrid micropillar array surfaces

Boiling heat transfer is widely used in numerous industrial systems because of its excellent heat dissipation performance. The critical heat flux (CHF) represents the upper limit of efficient heat removal. The most widely applied CHF model is Zuber's hydrodynamic instability model, which posits that the interference between vapor and liquid counterflows results in the CHF. In this paper, a three-dimensional (3D) hybrid micropillar array surface for hindering the liquid–vapor interface is proposed to enhance the CHF and heat transfer coefficient (HTC). The 3D hybrid surface comprises a tall hydrophobic (HB) pillar array surrounded by a short hydrophilic (HP) pillar array. The tall HB pillar array acts as a potential nucleation site for increasing the HTC, and the short HP pillar array acts as a wicking structure for facilitating liquid rewetting. In addition, the height difference between the HB and HP pillar arrays can potentially reduce the interference between liquid and vapor counterflows and increase the CHF. Relative to a plain silica surface, the 3D hybrid surface increases the CHF and HTC by 107% and 360%, respectively. Theoretical modeling suggests that the CHF on the 3D hybrid surface was still limited by pool hydrodynamics.

Molecular dynamics simulation of interfacial heat transfer behavior during the boiling of low-boiling-point organic fluid

Boiling heat transfer of low-boiling-point working fluid is one of the promising heat dissipation technologies commonly used in electronic equipment cooling. The study on interfacial heat transfer mechanism is of great significance to improve both heat removal capacity and heat dissipation efficiency. In this work, 1,1,1,2-tetrafluoroethane (R134a) was used as the working fluid, and the interfacial boiling behavior of R134a under different initial thicknesses of liquid film (δ0), solid-liquid interaction force and initial temperature (T0) was analyzed by molecular dynamics (MD) method. In the simulation, δ0 changes from 2.5 to 7.5 nm, the energy coefficient (α) characterizing solid-liquid interaction force changes from 0.25 to 4, T0 ranges from 180 to 200 K, and the temperature difference between the solid substrate and the liquid film is 220 K (i.e., the superheat for boiling). The results show that only δ0 could significantly alter the boiling mode of R134a.When δ0 = 2.5 nm, it shows thin film boiling, otherwise as δ0 gets thicker, it shifts into explosive boiling mode. Under the explosive boiling mode, the total thermal resistance (Rtot) consists of three parts, i.e., the solid interface thermal resistance (RW), thermal resistance of vapor layer (RV), and vapor-liquid interface thermal resistance above layer (RV-L). RV contributes the most to Rtot. However, during the thin film boiling, only RW contributes to Rtot. With the increase of δ0 and α, Rtot presents the similar trend, i.e., increasing first and decreasing afterwards. When δ0 = 7.5 nm and α is 1, Rtot reaches to the highest value of 3.49×10−8 K‧m2/W. For the influence of T0, it oppositely affects Rtot. Overall, the increase of α and T0 and the decrease of δ0 are all beneficial to accelerate the boiling of R134a.