A numerical-experimental investigation of pinned, water microdroplet evaporation on an isothermally heated flat surface was conducted by solving the axisymmetric and threedimensional time-dependent governing equations of continuity, momentum, energy, and species. The explicit volume of fluid (VOF) model with dynamic meshing and variable-time stepping on parallel processors was used to capture the time-dependent liquid-gas interface motion throughout the computational domain. The continuum surface force (CSF), the gravitational body force and Schrage's molecular kinetic-based evaporation model were included in the governing equations. The numerical model included temperature-and species-dependent thermodynamic and transport properties. A neutrally hydrophobic droplet was heated on an isothermal surface. Flow is generated in both liquid and gaseous phases due to the combined effects of surface tension gradients, evaporation, and heat conduction and convection. Throughout evaporation, the droplet exhibited higher temperature and curvature at the contact line and lower temperature and curvature at the droplet apex. These scalar gradients along the liquid-gas interface lead to thermalMarangoni driven re-circulating flow. This re-circulating flow causes the wall heat flux (q wall ) to peak both at the centerline and contact line, while the wall shear stress ( wall ) peaks at the contact line only. With decreasing accommodation coefficient (), the re-circulating flow is weakened. Consequently, q wall is reduced while wall remains almost unaffected, and the droplet lifetime () is extended. Thermal resistance is controlled by the evaporating liquid-gas interface and the time-averaged overall heat transfer coefficient ( ) decreases with decreasing The evaporation rate for water is sensitive to small changes in . Both 2 predictions and measurements indicate that the net effect of decreasing curvature ( 0 ) by increasing the droplet size is to increase . There are two reasons for this to happen: 1) larger droplets require more total energy than smaller droplets for full evaporation; 2) larger droplets, which consequently are exposed to lower surface tension forces than smaller droplets, develop weaker thermal-Marangoni flows yielding reduced flow re-circulation, which in turn diminishes q wall and . Moreover, the numerical results under-and overpredict for smaller and larger droplets, respectively, due to lack of de-pinning mechanism and the uncertainty in . Nonetheless, good qualitative comparisons are observed between simulations and experiments. The larger droplets (i.e., smaller 0 ) exhibit vapor bubble formation in their centerline. The experimental results indicate the three-dimensional flows at late stages of evaporation are coupled with the contact line receding motion. Lastly, the three-dimensional model was validated for early and intermediate stages of evaporation.
NomenclatureCa = Capillary number = V 0 / CHF = critical heat flux [W/m 2 ] c p = specific heat capacity at constant pressure [J/...