Zhenhai Pan

Assessment of Water Droplet Evaporation Mechanisms on Hydrophobic and Superhydrophobic Substrates

Evaporation rates are predicted and important transport mechanisms identified for evaporation of water droplets on hydrophobic (contact angle ~110°) and superhydrophobic (contact angle ~160°) substrates. Analytical models for droplet evaporation in the literature are usually simplified to include only vapor diffusion in the gas domain, and the system is assumed to be isothermal. In the comprehensive model developed in this study, evaporative cooling of the interface is accounted for, and vapor concentration is coupled to local temperature at the interface. Conjugate heat and mass transfer are solved in the solid substrate, liquid droplet, and surrounding gas. Buoyancy-driven convective flows in the droplet and vapor domains are also simulated. The influences of evaporative cooling and convection on the evaporation characteristics are determined quantitatively. The liquid-vapor interface temperature drop induced by evaporative cooling suppresses evaporation, while gas-phase natural convection acts to enhance evaporation. While the effects of these competing transport mechanisms are observed to counterbalance for evaporation on a hydrophobic surface, the stronger influence of evaporative cooling on a superhydrophobic surface accounts for an overprediction of experimental evaporation rates by ~20% with vapor diffusion-based models. The local evaporation fluxes along the liquid-vapor interface for both hydrophobic and superhydrophobic substrates are investigated. The highest local evaporation flux occurs at the three-phase contact line region due to proximity to the higher temperature substrate, rather than at the relatively colder droplet top; vapor diffusion-based models predict the opposite. The numerically calculated evaporation rates agree with experimental results to within 2% for superhydrophobic substrates and 3% for hydrophobic substrates. The large deviations between past analytical models and the experimental data are therefore reconciled with the comprehensive model developed here.

Influence of Surface Wettability on Transport Mechanisms Governing Water Droplet Evaporation

Prediction and manipulation of the evaporation of small droplets is a fundamental problem with importance in a variety of microfluidic, microfabrication, and biomedical applications. A vapor-diffusion-based model has been widely employed to predict the interfacial evaporation rate; however, its scope of applicability is limited due to incorporation of a number of simplifying assumptions of the physical behavior. Two key transport mechanisms besides vapor diffusion-evaporative cooling and natural convection in the surrounding gas-are investigated here as a function of the substrate wettability using an augmented droplet evaporation model. Three regimes are distinguished by the instantaneous contact angle (CA). In Regime I (CA ≲ 60°), the flat droplet shape results in a small thermal resistance between the liquid-vapor interface and substrate, which mitigates the effect of evaporative cooling; upward gas-phase natural convection enhances evaporation. In Regime II (60 ≲ CA ≲ 90°), evaporative cooling at the interface suppresses evaporation with increasing contact angle and counterbalances the gas-phase convection enhancement. Because effects of the evaporative cooling and gas-phase convection mechanisms largely neutralize each other, the vapor-diffusion-based model can predict the overall evaporation rates in this regime. In Regime III (CA ≳ 90°), evaporative cooling suppresses the evaporation rate significantly and reverses entirely the direction of natural convection induced by vapor concentration gradients in the gas phase. Delineation of these counteracting mechanisms reconciles previous debate (founded on single-surface experiments or models that consider only a subset of the governing transport mechanisms) regarding the applicability of the classic vapor-diffusion model. The vapor diffusion-based model cannot predict the local evaporation flux along the interface for high contact angle (CA ≥ 90°) when evaporative cooling is strong and the temperature gradient along the interface determines the peak local evaporation flux.

Spurious Current Suppression in VOF-CSF Simulation of Slug Flow through Small Channels

A numerical treatment is proposed to minimize the creation of unphysical, spurious currents in modeling liquid–gas slug flow using the volume of fluid-continuum surface force (VOF-CSF) method. An elongated gas slug drawn into a small circular channel initially filled with liquid is considered. To suppress spurious currents formed by numerical errors in calculation of the surface tension force at small capillary numbers (Ca < 0.01), an artificial relative reference frame is specified with motion in a direction opposite to the flow. An increase in the local relative velocity magnitude near the interface is demonstrated to be the key mechanism for spurious current suppression. A comparison of simulations performed with and without this treatment shows that spurious currents are eliminated at Ca = 0.0029; liquid film thickness, gas slug velocity, and liquid-phase circulation near the leading slug interface are preserved and the computed values agree with the literature. This demonstrates that the proposed moving reference frame method does not influence the computed physical phenomena of interest while suppressing unphysical spurious velocities.

A Cost-Effective Modeling Approach for Simulating Phase Change and Flow Boiling in Microchannels

High-fidelity simulation of flow boiling in microchannels remains a challenging problem, but the increasing interest in applications of microscale two-phase transport highlight its importance. In this paper, a volume of fluid (VOF)-based flow boiling model is proposed with computational expense-saving features that enable cost-effective simulation of two-phase flow and heat transfer in realistic geometries. The vapor and liquid phases are distinguished using a color function which represents the local volume fraction of the tracked phase. Mass conservation is satisfied by solving the transport equations for both phases with a finite-volume approach. In order to predict phase change at the liquid-vapor interface, evaporative heat and mass source terms are calculated using a novel, saturated-interface-volume phase change model that fixes the interface at the saturation temperature at each time step to achieve stability. Numerical oscillation of the evaporation source terms is thus eliminated and a non-iterative time advancement scheme can be adopted to reduce computational cost. The reference frame is set to move with the vapor slug to artificially increase the local velocity magnitude in the thin liquid film region in the relative frame, which reduces the influence of numerical errors resulting from calculation of the surface tension force, and thus suppresses the development of spurious currents. This allows use of non-uniform meshes that can efficiently resolve high-aspect-ratio geometries and flow features and significantly reduces the overall numerical expense. The proposed model is used to simulate the growth of a vapor bubble in a heated 2D axisymmetric microchannel. The bubble motion, bubble growth rate, liquid film thickness, and local heat transfer coefficient along the wall are compared against previous numerical studies.Copyright

Numerical Study of Water Droplet Evaporation on a Superhydrophobic Surface

A comprehensive numerical model is developed to predict evaporation of a water droplet from an unheated superhydrophobic substrate. Analytical models that only consider vapor diffusion in the gas domain, and assume the system to be isothermal, over-predict the evaporation rates by ∼25% compared to experiments conducted on such surfaces. The current model solves for conjugate heat and mass transfer in the solid substrate, liquid droplet, and surrounding gas. Evaporative cooling of the interface is accounted for, and vapor concentration is coupled to local temperature at the interface. Buoyancy-driven convective flows in the droplet and vapor domains are also simulated. A droplet evaporating in a constant-contact-angle mode with an initial volume of 3 μl and contact angle of 160 deg is considered at an ambient temperature of 21°C and 29% relative humidity, to match conditions of related experiments. The interface cooling effect suppresses the evaporation rate significantly; however, natural convection in the gas and liquid domains has a negligible impact on the evaporation rate. The local evaporation flux along the droplet interface predicted by the model is compared to that predicted by an analytical diffusion-based model. The numerically calculated total evaporation rate agrees with experimental results to within 2%. The large deviations between past analytical models and the experimental data on superhydrophobic surfaces are reconciled.Copyright