Numerical predictions of flow boiling characteristics: Current status, model setup and CFD modeling for different non-uniform heating profiles

Abstract A detailed analysis of two-phase flow boiling characteristics inside high pressure systems is presented focusing on non-uniform axial heating profiles. For this purpose, a detailed numerical model has been developed after presenting the current status of the use of CFD techniques in flow boiling predictions. User defined functions written in C++ were compiled and hooked to the software in order to account for mass interaction between phases using the Eulerian multiphase flow model. The modeled domain is a 2 m long stainless steel pipe with inside and outside diameters of 15.4 mm, 25.4 mm, respectively. The base imposed uniform heat flux was 345.6 kW/m 2 , for mass flow rate of water of 0.161 kg/s and at a temperature of 200 °C. The model was validated against range of experimental data and the results are very promising for the use of CFD in flow boiling characterization. Effects of increasing uniform heat flux were considered for different increments of 30, 50 and 75% as reference to the basic applied heat flux. The influences of heat flux profile in the axial direction were investigated while maintaining the same total power. Different heat flux profiles of linearly increasing, linearly decreasing, sine, and cosine shapes were considered.

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

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

[3]  Anandaroop Bhattacharya,et al.  Numerical modelling of boiling heat transfer in microchannels , 2009 .

[4]  Kazuo Ikeda,et al.  Development of Mitsubishi High Thermal Performance Grid , 2002 .

[5]  V. I. Tolubinsky,et al.  VAPOUR BUBBLES GROWTH RATE AND HEAT TRANSFER INTENSITY AT SUBCOOLED WATER BOILING , 1970 .

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

[7]  J. Flaherty,et al.  Analysis of phase distribution in fully developed laminar bubbly two-phase flow , 1991 .

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

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

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

[11]  E. Lorenzini,et al.  A simplified method proposal for practical determination of aperiodic two-phase flow instability , 1981 .

[12]  Soon Heung Chang,et al.  CHF experiment and CFD analysis in a 2 × 3 rod bundle with mixing vane , 2009 .

[13]  Rongguo Huang,et al.  Numerical investigation of boiling flow of nitrogen in a vertical tube using the two-fluid model , 2006 .

[14]  Yanzhong Li,et al.  CFD simulation of upward subcooled boiling flow of refrigerant-113 using the two-fluid model , 2009 .

[15]  K. Ikeda,et al.  Single-phase CFD applicability for estimating fluid hot-spot locations in a 5 × 5 fuel rod bundle , 2006 .