A Ghost Fluid/Level Set Method for boiling flows and liquid evaporation: Application to the Leidenfrost effect

The development of numerical methods for the direct numerical simulation of two-phase flows with phase change, in the framework of interface capturing or interface tracking methods, is the main topic of this study. We propose a novel numerical method, which allows dealing with both evaporation and boiling at the interface between a liquid and a gas. Indeed, in some specific situations involving very heterogeneous thermodynamic conditions at the interface, the distinction between boiling and evaporation is not always possible. For instance, it can occur for a Leidenfrost droplet; a water drop levitating above a hot plate whose temperature is much higher than the boiling temperature. In this case, boiling occurs in the film of saturated vapor which is entrapped between the bottom of the drop and the plate, whereas the top of the water droplet evaporates in contact of ambient air. The situation can also be ambiguous for a superheated droplet or at the contact line between a liquid and a hot wall whose temperature is higher than the saturation temperature of the liquid. In these situations, the interface temperature can locally reach the saturation temperature (boiling point), for instance near a contact line, and be cooler in other places. Thus, boiling and evaporation can occur simultaneously on different regions of the same liquid interface or occur successively at different times of the history of an evaporating droplet. Standard numerical methods are not able to perform computations in these transient regimes, therefore, we propose in this paper a novel numerical method to achieve this challenging task. Finally, we present several accuracy validations against theoretical solutions and experimental results to strengthen the relevance of this new method.

[1]  R. Fedkiw,et al.  A boundary condition capturing method for incompressible flame discontinuities , 2001 .

[2]  Ronald Fedkiw,et al.  A Boundary Condition Capturing Method for Multiphase Incompressible Flow , 2000, J. Sci. Comput..

[3]  Yohei Sato,et al.  A sharp-interface phase change model for a mass-conservative interface tracking method , 2013, J. Comput. Phys..

[4]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[5]  S. Zaleski,et al.  DIRECT NUMERICAL SIMULATION OF FREE-SURFACE AND INTERFACIAL FLOW , 1999 .

[6]  Frédéric Gibou,et al.  Efficient symmetric discretization for the Poisson, heat and Stefan-type problems with Robin boundary conditions , 2010, J. Comput. Phys..

[7]  M. Ishii,et al.  Thermo-Fluid Dynamics of Two-Phase Flow , 2007 .

[8]  S. Osher,et al.  A level set approach for computing solutions to incompressible two-phase flow , 1994 .

[9]  R. Fedkiw,et al.  A Boundary Condition Capturing Method for Poisson's Equation on Irregular Domains , 2000 .

[10]  Gihun Son,et al.  A Level-Set Method for Analysis of Particle Motion in an Evaporating Microdroplet , 2015 .

[11]  S. Welch,et al.  A Volume of Fluid Based Method for Fluid Flows with Phase Change , 2000 .

[12]  Liang-Shih Fan,et al.  Three-dimensional direct numerical simulation for film-boiling contact of moving particle and liquid droplet , 2006 .

[13]  G. Son,et al.  Numerical Simulation of Film Boiling Near Critical Pressures With a Level Set Method , 1998 .

[14]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[15]  S. Hardt,et al.  Evaporation model for interfacial flows based on a continuum-field representation of the source terms , 2008, J. Comput. Phys..

[16]  Dieter Bothe,et al.  Numerical modeling of thermocapillary two-phase flows with evaporation using a two-scalar approach for heat transfer , 2013, J. Comput. Phys..

[17]  Frédéric Risso,et al.  Effect of rising motion on the damped shape oscillations of drops and bubbles , 2013 .

[18]  L. Fan,et al.  Three-dimensional simulation of impingement of a liquid droplet on a flat surface in the Leidenfrost regime , 2005 .

[19]  V. Carey Liquid-Vapor Phase-Change Phenomena: An Introduction to the Thermophysics of Vaporization and Condensation Processes in Heat Transfer Equipment, Third Edition , 2020 .

[20]  Brian Spalding,et al.  Combustion and mass transfer , 1979 .

[21]  Frédéric Risso,et al.  On the computation of viscous terms for incompressible two-phase flows with Level Set/Ghost Fluid Method , 2015, J. Comput. Phys..

[22]  Michel Gradeck,et al.  Energy balance of droplets impinging onto a wall heated above the Leidenfrost temperature , 2013 .

[23]  Peter Stephan,et al.  Numerical simulation of the transient heat transfer during nucleate boiling of refrigerant HFE-7100 , 2010 .

[24]  Fabrice Lemoine,et al.  High-speed shadow imagery to characterize the size and velocity of the secondary droplets produced by drop impacts onto a heated surface , 2013 .

[25]  Sébastien Tanguy,et al.  Benchmarks and numerical methods for the simulation of boiling flows , 2014, J. Comput. Phys..

[26]  G. Son,et al.  Dynamics and Heat Transfer Associated With a Single Bubble During Nucleate Boiling on a Horizontal Surface , 1999 .

[27]  David F. Fletcher,et al.  A hydrodynamic and thermodynamic simulation of droplet impacts on hot surfaces, Part I: theoretical model , 2001 .

[28]  Li-Tien Cheng,et al.  A second-order-accurate symmetric discretization of the Poisson equation on irregular domains , 2002 .

[29]  J. Dendy Black box multigrid , 1982 .

[30]  T. Aslam A partial differential equation approach to multidimensional extrapolation , 2004 .

[31]  Frédéric Gibou,et al.  A level set based sharp interface method for the multiphase incompressible Navier-Stokes equations with phase change , 2007, J. Comput. Phys..

[32]  Mark Sussman,et al.  A sharp interface method for incompressible two-phase flows , 2007, J. Comput. Phys..

[33]  D. Juric,et al.  Computations of Boiling Flows , 1998 .

[34]  Cornelis Vuik,et al.  Fast and robust solvers for pressure-correction in bubbly flow problems , 2008, J. Comput. Phys..

[35]  Kenneth K. Kuo,et al.  A ghost fluid method for compressible reacting flows with phase change , 2013, J. Comput. Phys..

[36]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[37]  Ronald Fedkiw,et al.  High Resolution Sharp Computational Methods for Elliptic and Parabolic Problems in Complex Geometries , 2013, J. Sci. Comput..

[38]  Sébastien Tanguy,et al.  A Level Set Method for vaporizing two-phase flows , 2007, J. Comput. Phys..

[39]  Xianguo Li,et al.  A mass transfer correlation for droplet evaporation in high-temperature flows , 1991 .

[40]  Bernhard Weigand,et al.  Direct numerical simulation of evaporating droplets , 2008, J. Comput. Phys..

[41]  Nikos Nikolopoulos,et al.  A numerical investigation of the evaporation process of a liquid droplet impinging onto a hot substrate , 2007 .

[42]  Gretar Tryggvason,et al.  A front tracking method for computations of boiling in complex geometries , 2004 .

[43]  R. Fedkiw,et al.  A fourth order accurate discretization for the Laplace and heat equations on arbitrary domains, with applications to the Stefan problem , 2005 .

[44]  B. Duret,et al.  DNS analysis of turbulent mixing in two-phase flows , 2012 .

[45]  Frédéric Gibou,et al.  Geometric integration over irregular domains with application to level-set methods , 2007, J. Comput. Phys..

[46]  David F. Fletcher,et al.  A hydrodynamic and thermodynamic simulation of droplet impacts on hot surfaces, Part II: validation and applications , 2001 .

[47]  Frédéric Gibou,et al.  Robust second-order accurate discretizations of the multi-dimensional Heaviside and Dirac delta functions , 2008, J. Comput. Phys..