Abstract: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.

Criterion for Reversal of Thermal Marangoni Flow in Drying Drops

Temperature Distribution Along The Surface Of Evaporating Droplets

Radial-Velocity Profile along the Surface of Evaporating Liquid Droplets

“Phase Diagram” of Surface Temperature Distribution of Sessile Droplets and the Effects of Evaporative Cooling

Evaporating sessile droplets: solutal Marangoni effects overwhelm thermal Marangoni flow

Numerical Investigation of Wafer Drying Induced by the Thermal Marangoni Effect

Marangoni Flow In An Evaporating Water Droplet

Role of surfactant-induced Marangoni effects in droplet dynamics on a solid surface in shear flow

Instability Patterns of Marangoni Flow in Evaporating Droplets on Lyophobic Surface

Study of Sessile Droplet Evaporation on Wall with Radial Temperature Gradient

Theoretical Analysis of a Sessile Evaporating Droplet on a Curved Substrate with an Interfacial Cooling Effect

Sessile Droplet Evaporation on Wall with Radial Temperature Gradient

Influence of surface wettability on transport mechanisms governing water droplet evaporation

Influence of thermal properties on hydrothermal waves in evaporating sessile droplets

Marangoni Convection Instability Inside an Evaporating Droplet on an Inclined Substrate

Autothermotaxis of volatile drops

Self-sustained Movement of Condensate Drop on Uniform Temperature Surface During Marangoni Dropwise Condensation

Three-Dimensional Numerical Simulation on Marangoni Convection in a Sessile Water Droplet Evaporating in Its Vapor at Low Pressure

Evaporation Flux Distribution of Drops on a Hydrophilic or Hydrophobic Flat Surface by Molecular Simulations.

Numerical Investigation of Terminal Residual Thickness of the Water Film and Surface Flow in the Entrainment Region in Marangoni Drying