Numerical modeling of thermocapillary two-phase flows with evaporation using a two-scalar approach for heat transfer

A one-field model is derived from the sharp interface continuum mechanical balances for two-phase evaporative and thermocapillary flows. Emphasis is put on a clear distinction of the different velocities at the interface which appear due to phase transfer. The one-field model is solved numerically within a Finite Volume scheme and the interface is captured using an extended Volume of Fluid method, where the interface is reconstructed linearly with the PLIC technique. The numerical heat transfer is based on a two-scalar approach where two separate temperature fields are used for the temperature inside the two phases. This results in an accurate treatment of the interfacial heat transfer, specifically the interface temperature which is crucial numerically, both for evaporation and thermocapillarity. The method is validated for two-phase heat conduction, with analytical solution in case of no evaporation and with experimental measurement in case of incorporated evaporation effect. The method is applied to realistic cases dealing with non-uniformly heated thin liquid films, i.e. liquid films on (i) structured heated substrates and (ii) locally heated substrates. The numerical predictions in terms of flow pattern, surface deformation, temperature and velocity are compared with experiments conducted at the Universite Libre de Bruxelles for (i) and at the Technische Universitat Darmstadt for (ii). Qualitative agreement is achieved and shows the potential of this approach to simulate thermocapillary flows with dynamically deformable interfaces combined with evaporation.

[1]  Suresh V. Garimella,et al.  MOTION OF INTERACTING GAS BUBBLES IN A VISCOUS LIQUID INCLUDING WALL EFFECTS AND EVAPORATION , 1997 .

[2]  Murray Rudman,et al.  VOLUME-OF-FLUID CALCULATION OF HEAT OR MASS TRANSFER ACROSS DEFORMING INTERFACES IN TWO-FLUID FLOW , 2002 .

[3]  John C. Slattery,et al.  Advanced transport phenomena , 1999 .

[4]  J. Brackbill,et al.  A continuum method for modeling surface tension , 1992 .

[5]  Jorge J. Moré,et al.  User Guide for Minpack-1 , 1980 .

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

[7]  Martin Wörner,et al.  Volume-averaged conservation equations for volume-of-fluid interface tracking , 2001 .

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

[9]  W. Sabisch,et al.  Dreidimensionale numerische Simulation der Dynamik von aufsteigenden Einzelblasen und Blasenschwärmen mit einer Volume-of-Fluid-Methode [online] , 2000 .

[10]  Peter Stephan,et al.  Modification and extension of a standard volume-of-fluid solver for simulating boiling heat transfer , 2010 .

[11]  Oleg Kabov,et al.  Study of thermocapillary film rupture using a fiber optical thickness probe , 2007 .

[12]  Dieter Bothe,et al.  Direct numerical simulation of thermocapillary flow based on the Volume of Fluid method , 2011 .

[13]  Peter Stephan,et al.  A numerical model for the thermocapillary flow and heat transfer in a thin liquid film on a microstructured wall , 2007 .

[14]  S. Zaleski,et al.  Modelling Merging and Fragmentation in Multiphase Flows with SURFER , 1994 .

[15]  Oleg Kabov,et al.  Shear-driven flows of locally heated liquid films , 2008 .

[16]  Dieter Bothe,et al.  Volume of Fluid Direct Numerical Simulation of Heat and Mass Transfer using Sharp Temperature and Concentration Fields , 2011 .

[17]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[18]  P. Sagaut.  BOOK REVIEW: Large Eddy Simulation for Incompressible Flows. An Introduction , 2001 .

[19]  Ilya Prigogine,et al.  Surface tension and adsorption , 1966 .

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

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

[22]  Manolis Gavaises,et al.  Numerical investigation of the evaporation of two-component droplets , 2011 .

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

[24]  Dieter Bothe,et al.  A VOF-Based Conservative Method for the Simulation of Reactive Mass Transfer from Rising Bubbles , 2011 .

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

[26]  Dieter Bothe,et al.  3D Numerical Modeling of Soluble Surfactant at Fluidic Interfaces Based on the Volume-of-Fluid Method , 2009 .

[27]  Suman Chakraborty,et al.  Mechanics Over Micro and Nano Scales , 2011 .

[28]  Rho-Taek Jung,et al.  Numerical simulation of high Schmidt number flow over a droplet by using moving unstructured mesh , 2005 .

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

[30]  R. Clift,et al.  Bubbles, Drops, and Particles , 1978 .

[31]  W. Rider,et al.  Reconstructing Volume Tracking , 1998 .

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

[33]  Hugo A. Jakobsen,et al.  Chemical Reactor Modeling: Multiphase Reactive Flows , 2008 .

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

[35]  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..

[36]  Dmitry V. Zaitsev,et al.  An experimental modeling of gravity effect on rupture of a locally heated liquid film , 2007 .