A sharp-interface phase change model for a mass-conservative interface tracking method

A new phase-change model has been developed for a mass-conservative interface tracking method. The mass transfer rate is directly calculated from the heat flux at the liquid-vapor interface, and the phase change takes place only in the cells which include this interface. As a consequence of the sharpness of the mass transfer rate distribution, the velocity jump across the interface can be captured, and high accuracy can be maintained. The method has been implemented in an incompressible Navier-Stokes equations solver employing a projection method based on a staggered finite-volume algorithm on Cartesian grids. The model has been verified for one-dimensional phase-change problems and a three-dimensional simulation of a growing vapor bubble in a superheated liquid under zero gravity condition. The computed results agree with theoretical solutions, and the accuracy of the model is confirmed to be of second-order in space using a grid refinement study. A three-dimensional simulation of a rising vapor bubble in a superheated liquid under gravity has been performed as a validation case, and good agreement with experimental data is obtained for the bubble growth rate. As a demonstration of the applicability of the method to engineering problems, a nucleate boiling simulation is presented with a comparison to experimental data. Good agreement is obtained for the bubble shapes and the bubble departure period. In all the simulation cases, strict mass conservation is satisfied.

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

[2]  Samuel W. J. Welch,et al.  Local simulation of two-phase flows including interface tracking with mass transfer , 1995 .

[3]  Arnoldo Badillo,et al.  Quantitative phase-field modeling for boiling phenomena. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[4]  Jules Thibault,et al.  Direct numerical simulation of interphase heat and mass transfer in multicomponent vapour–liquid flows , 2010 .

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

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

[7]  Dipankar Sanyal,et al.  A Sphericosymmetric VOF Approach for Investigating Immiscible Two-Phase Systems with One Liquid Phase , 2006 .

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

[9]  G. Kreiss,et al.  A conservative level set method for two phase flow II , 2005, Journal of Computational Physics.

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

[11]  Gretar Tryggvason,et al.  Computations of Boiling Flows , 2004 .

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

[13]  G. Son.  A Numerical Method for Bubble Motion with Phase Change , 2001 .

[14]  Timothy G. Trucano,et al.  Verification and validation. , 2005 .

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

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

[17]  Dipankar Sanyal,et al.  An integral approach for simulation of vapour film dynamics around a spherical surface , 2009 .

[18]  C. A. Busse,et al.  Analysis of the heat transfer coefficient of grooved heat pipe evaporator walls , 1992 .

[19]  R. Schrage.  A theoretical study of interphase mass transfer , 1953 .

[20]  E. G. Keshock,et al.  Effects of reduced gravity on nucleate boiling bubble dynamics in saturated water , 1964 .

[21]  P. Roe.  CHARACTERISTIC-BASED SCHEMES FOR THE EULER EQUATIONS , 1986 .

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

[23]  G. Son,et al.  Numerical simulation of nucleate boiling on a horizontal surface at high heat fluxes , 2008 .

[24]  Matthew W. Williams,et al.  A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework , 2006, J. Comput. Phys..

[25]  Gretar Tryggvason,et al.  Numerical simulation of dendritic solidification with convection: two-dimensional geometry , 2002 .

[26]  Gretar Tryggvason,et al.  Direct numerical simulations of flows with phase change , 2005 .

[27]  Peter Stephan,et al.  CFD Simulation of Boiling Flows Using the Volume-of-Fluid Method within OpenFOAM , 2009 .

[28]  Gaurav Tomar,et al.  Numerical simulation of bubble growth in film boiling using a coupled level-set and volume-of-fluid method , 2005 .

[29]  T. Yabe,et al.  Exactly conservative semi-Lagrangian scheme for multi-dimensional hyperbolic equations with directional splitting technique , 2001 .

[30]  V. S. Vaidhyanathan,et al.  Transport phenomena , 2005, Experientia.

[31]  G. Son,et al.  Numerical Simulation of Bubble Merger Process on a Single Nucleation Site During Pool Nucleate Boiling , 2002 .

[32]  L. W. Florschuetz,et al.  Growth rates of free vapor bubbles in liquids at uniform superheats under normal and zero gravity conditions , 1969 .

[33]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

[34]  L. E. Scriven.  On the dynamics of phase growth , 1995 .

[35]  Bojan Niceno,et al.  A conservative local interface sharpening scheme for the constrained interpolation profile method , 2012 .

[36]  Abhijit Mukherjee,et al.  Numerical Study of Lateral Merger of Vapor Bubbles During Nucleate Pool Boiling , 2003 .

[37]  O. Lebaigue,et al.  The second gradient method for the direct numerical simulation of liquid—vapor flows with phase change , 2001 .

[38]  Richard Edwin Sonntag,et al.  Fundamentals of Thermodynamics , 1998 .

[39]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

[40]  Samuel W. J. Welch,et al.  NUMERICAL COMPUTATION OF FILM BOILING INCLUDING CONJUGATE HEAT TRANSFER , 2002 .

[41]  Edip Can,et al.  A level set method for vapor bubble dynamics , 2012, J. Comput. Phys..

[42]  Vijay K. Dhir,et al.  Numerical Simulations of the Dynamics and Heat Transfer Associated With a Single Bubble in Subcooled Pool Boiling , 2010 .

[43]  Vijay K. Dhir,et al.  Three-dimensional simulation of saturated film boiling on a horizontal cylinder , 2008 .

[44]  Stanley J. Kleis,et al.  A Volume of Fluid Phase Change Model on Adaptive Octree Grids , 2011 .

[45]  V. Dhir.  Mechanistic Prediction of Nucleate Boiling Heat Transfer-Achievable or a Hopeless Task? , 2006 .

[46]  Yanhua Yang,et al.  Numerical simulation of film boiling on a sphere with a volume of fluid interface tracking method , 2008 .

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

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

[49]  Samuel W. J. Welch,et al.  Planar Simulation of Bubble Growth in Film Boiling in Near-Critical Water Using a Variant of the VOF Method , 2004 .

[50]  Ronald Fedkiw,et al.  Extrapolation in the Normal Direction , 2003 .

[51]  Nobuhide Kasagi,et al.  Modeling and Design of Micro-Grooved Flat Plate Evaporator , 2005 .

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

[53]  Arnoldo Badillo,et al.  Quantitative phase-field modeling for wetting phenomena. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[54]  A. Chorin.  Numerical Solution of the Navier-Stokes Equations* , 1989 .

[55]  Vijay K. Dhir,et al.  A Level Set Method for Analysis of Film Boiling on an Immersed Solid Surface , 2007 .

[56]  Patrick J. Roache,et al.  Verification and Validation in Computational Science and Engineering , 1998 .

[57]  Yohei Sato,et al.  A new contact line treatment for a conservative level set method , 2012, J. Comput. Phys..