CFD for subcooled flow boiling: Coupling wall boiling and population balance models

Abstract In this work we investigate the present capabilities of CFD for wall boiling. The computational model used combines the Euler/Euler two-phase flow description with heat flux partitioning. Very similar modelling was previously applied to boiling water under high pressure conditions relevant to nuclear power systems. Similar conditions in terms of the relevant non-dimensional numbers have been realized in the DEBORA tests using dichlorodifluoromethane (R12) as the working fluid. This facilitated measurements of radial profiles for gas volume fraction, gas velocity, liquid temperature and bubble size. Essential for the momentum, mass and energy exchange between the phases is an adequate description of the interfacial area or respectively the bubble size. In a previous study ( Krepper and Rzehak, 2011 ) it was shown that a monodisperse bubble size representation is not sufficient to this end. Therefore, in the present work a population balance approach is used, where bubbles are generated at the wall with a certain size that subsequently evolves due to both condensation/evaporation and coalescence/breakup processes. The results show the potential of this approach which is able to describe the observed bubble size increase after leaving the wall as well as the change of gas volume fraction profile from wall to core peaking with increasing inlet temperature.

[1]  H. Blanch,et al.  Bubble coalescence and break‐up in air‐sparged bubble columns , 1990 .

[2]  D. Kenning,et al.  Subcooled flow boiling at high heat flux , 1985 .

[3]  N. Zuber,et al.  Drag coefficient and relative velocity in bubbly, droplet or particulate flows , 1979 .

[4]  S. Pope Turbulent Flows: FUNDAMENTALS , 2000 .

[5]  Nikolay Ivanov Kolev UNIQUENESS OF THE ELEMENTARY PHYSICS DRIVING HETEROGENEOUS NUCLEATE BOILING AND FLASHING , 2006 .

[6]  A. K. Agrawal,et al.  Shape of liquid drops moving in liquid media , 1966 .

[7]  F. Menter Two-equation eddy-viscosity turbulence models for engineering applications , 1994 .

[8]  W. E. Ranz,et al.  Evaporation from drops , 1952 .

[9]  Eckhard Krepper,et al.  The inhomogeneous MUSIG model for the simulation of polydispersed flows , 2008 .

[10]  Eckhard Krepper,et al.  CFD for subcooled flow boiling: Simulation of DEBORA experiments , 2011 .

[11]  H. Svendsen,et al.  Theoretical model for drop and bubble breakup in turbulent dispersions , 1996 .

[12]  Etienne Manon,et al.  Contribution à l'analyse et à la modélisation locale des écoulements bouillants sous-saturés dans les conditions des réacteurs à eau sous pression , 2000 .

[13]  Eckhard Krepper,et al.  A population balance approach considering heat and mass transfer—Experiments and CFD simulations , 2011 .

[14]  Michael Z. Podowski,et al.  MULTIDIMENSIONAL EFFECTS IN FORCED CONVECTION SUBCOOLED BOILING , 1990 .

[15]  P. Griffith,et al.  The mechanism of heat transfer in nucleate pool boiling—Part I: Bubble initiaton, growth and departure , 1965 .

[16]  Peter Griffith,et al.  THE MECHANISM OF HEAT TRANSFER IN NUCLEATE POOL BOILING. Technical Report No. 19 , 1962 .

[17]  F. White Viscous Fluid Flow , 1974 .

[18]  Yassin A. Hassan,et al.  A two-equation turbulence model of turbulent bubbly flows , 2001 .

[19]  G. Cubizolles,et al.  LOCAL MEASUREMENTS ON FLOW BOILING OF REFRIGERANT 12 IN A VERTICAL TUBE , 2001 .

[20]  J. Converti,et al.  A model for turbulent polydisperse two-phase flow in vertical channels , 2003 .

[21]  Eckhard Krepper,et al.  Use of models for lift, wall and turbulent dispersion forces acting on bubbles for poly-disperse flows , 2007 .

[22]  Michio Sadatomi,et al.  Momentum and heat transfer in two-phase bubble flow—I. Theory , 1981 .

[23]  Günter Brenn,et al.  Modelling of the Near-Wall Liquid Velocity Field in Subcooled Boiling Flow , 2005 .

[24]  O. C. Jones,et al.  The Prediction of Two-Phase Turbulence and Phase Distribution Phenomena Using a K-κ Model , 1989 .

[25]  C. Lifante,et al.  Coupling of wall boiling with discrete population balance model , 2012 .

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

[27]  Th. Frank,et al.  Extension of the inhomogeneous MUSIG model for bubble condensation , 2011 .

[28]  B. A. Kader Temperature and concentration profiles in fully turbulent boundary layers , 1981 .

[29]  Wei Yao,et al.  Prediction of Parameters Distribution of Upward Boiling Two-Phase Flow With Two-Fluid Models , 2002 .

[30]  H. J. Ivey Relationships between bubble frequency, departure diameter and rise velocity in nucleate boiling , 1967 .

[31]  Janusz A. Szymczyk,et al.  Flow field around growing and rising vapour bubble by PIV measurement , 2005, J. Vis..

[32]  Th. Frank,et al.  VALIDATION OF CFD MODELS FOR MONO- AND POLYDISPERSE AIR-WATER TWO-PHASE FLOWS IN PIPES , 2008 .

[33]  D. Drew,et al.  Theory of Multicomponent Fluids , 1998 .

[34]  Horst-Michael Prasser,et al.  Steam bubble condensation in sub-cooled water in case of co-current vertical pipe flow , 2007 .

[35]  S. Becker,et al.  Modelling and simulation of the dynamic flow behaviour in a bubble column , 2001 .

[36]  I. Zun,et al.  The transverse migration of bubbles influenced by walls in vertical bubbly flow , 1980 .

[37]  Eckhard Krepper,et al.  CFD modelling of subcooled boiling—Concept, validation and application to fuel assembly design , 2007 .

[38]  B. Mikic,et al.  A New Correlation of Pool-Boiling Data Including the Effect of Heating Surface Characteristics , 1969 .

[39]  Boštjan Končar,et al.  Wall function approach for boiling two-phase flows , 2010 .

[40]  Robert Cole,et al.  A photographic study of pool boiling in the region of the critical heat flux , 1960 .

[41]  Hidesada Tamai,et al.  Transverse migration of single bubbles in simple shear flows , 2002 .