Two- and three-dimensional numerical simulations of sessile droplet evaporating on heated substrate (original) (raw)
Related papers
International Journal of Heat and Mass Transfer, 2007
A numerical investigation of the evaporation process of n-heptane and water liquid droplets impinging onto a hot substrate is presented. Three different temperatures are investigated, covering flow regimes below and above Leidenfrost temperature. The Navier–Stokes equations expressing the flow distribution of the liquid and gas phases, coupled with the Volume of Fluid Method (VOF) for tracking the liquid–gas interface, are solved numerically using the finite volume methodology. Both two-dimensional axisymmetric and fully three-dimensional domains are utilized. An evaporation model coupled with the VOF methodology predicts the vapor blanket height between the evaporating droplet and the substrate, for cases with substrate temperature above the Leidenfrost point, and the formation of vapor bubbles in the region of nucleate boiling regime. The results are compared with available experimental data indicating the outcome of the impingement and the droplet shape during the impingement process, while additional information for the droplet evaporation rate and the temperature and vapor concentration fields is provided by the computational model.
Modeling of Evaporation From a Sessile Constant Shape Droplet
ASME 2017 15th International Conference on Nanochannels, Microchannels, and Minichannels, 2017
In this study, a computational model for the evaporation from a sessile liquid droplet fed from the center to keep the diameter of the droplet constant is presented. The continuity, momentum and energy equations are solved with temperature dependent thermo-physical properties using COMSOL Multi-physics. At the surface of the droplet, convective heat and evaporative mass fluxes are assigned. Since the flow field is affected by evaporative flux, an iterative scheme is built and the computation is automated using COMSOL-MATLAB interface. Correlations are implemented to predict the convective heat transfer coefficients and evaporative flux. Three different wall temperatures are used in simulations. The results show that the flow inside the droplet is dominated by buoyancy when the effect of the thermo-capillarity is neglected. The resulting flow generates a circulation pattern emerging from the entrance to the apex, along the surface of the droplet to the bottom heated wall and back to ...