Modeling of Upward Subcooled Flow Boiling of Refrigerant-113 in a Vertical Annulus

Document Type : Research Paper


1 Mechanical Engineering Department, University of Tabriz , Iran

2 Mechanical Engineering Department, University of Tabriz, Iran


In this paper, a modified two-fluid model has been adopted to simulate the process of upward vertical subcooled flow boiling of refrigerant R-113 in a vertical annular channel at low pressure. The modified model considers temperature dependent properties and saturation temperature variation and was validated against a number of published low-pressure subcooled boiling experiments. The results show good agreement with experimental data including radial profiles of void fraction, phase velocities and liquid temperature. A sinusoidal axial distribution of wall heat flux was applied as well as constant wall heat flux. Results show that by increasing the wall heat flux, the bubble boundary layer will become thicker and the profiles of axial liquid velocity will gradually depart from those of single-phase flow.


Main Subjects

       Kurul, N., and Podowski, M. Z., “On the Modeling of Multi-dimensional Effects in Boiling Channels,” in: 27th National Heat Transfer Conference, Minneapolis, (1991).
[2]   Kurul, N., and Podowski, M.Z., “Multidimensional Effects in Forced Convection Subcooled Boiling,” in: Proceedings of the Ninth International Heat Transfer Conference, Jerusalem, Israel, pp. 21–26, (1990).
[3]   Mimouni, S., Boucker, M., Laviéville, J., Guelfi, A., and Bestion, D., “Modelling and Computation of Cavitation and Boiling Bubbly Flows with the NEPTUNE CFD Code,” Nuclear Engineering and Design, Vol. 238, pp. 680–692, (2008).
[4]   Krepper, E., Koncar, B., and Egorov, Y., “CFD Modeling of Subcooled Boiling-concept, Validation and Application to Fuel Assembly Design,” Nuclear Engineering and Design, Vol. 237, pp. 716–731, (2007).
[5]   Lai, J.C., and Farouk, B., “Numerical Simulation of Subcooled Boiling and Heat Transfer in Vertical Ducts,” International Journal of Heat and Mass Transfer, Vol. 36, pp. 1541–1551, (1993).
[6]   Anglart, H., and Nylund, O., “CFD Application to Prediction of Void Distribution in Two-phase Bubbly Flows in Rod Bundles,” Nuclear Engineering and Design, Vol. 163, pp. 81–98, (1996).
[7]   Hari, S., and Hassan, Y.A., “Refinement of the RELAP5/MOD3.2 Subcooled Boiling Model for Low-pressure Conditions,” In: Proceedings of the 8th International Conference ICONE-8, Baltimore, MD, USA, (2000).
[8]   Koncar, B., Kljenak, I., and Mavko, B., “Modelling of Local Two-phase Flow Parameters in Upward Subcooled Flow Boiling at Low Pressure” International Journal of Heat and Mass Transfer, Vol. 47, pp. 1499–1513, (2004).
[9]   Hainoun, A., Hicken, E., and Wolters, J., “Modeling of Void Formation in the Subcooled Boiling Regime in the ATHLET Code to Simulate Flow Instability for Research Reactors,” Nuclear Engineering and Design, Vol. 167, pp. 175–191, (1996). 
[10] Janssens-Maenhout, G., Knebel, J.U., and Mueller, U., “Subcooled Nucleate Boiling at Low Pressure and Low Heat Flux,” in: Proceedings of 3rd International Conference on Multiphase Flow, Lyon, France, (1998).
[11] Lee, T.H., Park, G.C., and Lee, D.J., “Local Flow Characteristics of Subcooled Boiling Flow of Water in a Vertical Concentric Annulus,” International Journal of Multiphase Flow, Vol. 28, pp. 1351–1368, (2002).
[12] Legendre, D., Boree, J., and Magnaudet, J., “Thermal and Dynamic Evolution of a Spherical Bubble Moving Steadily in a Superheated or Subcooled Liquid,” Physics of Fluids, Vol. 10, pp. 1256–1272, (1998).
[13] Kljenak, I., Park, G.C., Mavko, B., and Lee, T., “Bubble-tracking Modeling of Subcooled Nucleate Boiling in a Vertical Annulus,” in: Proceedings of the 12th International Heat Transfer Conference, Grenoble, France, pp. 767–772, (2002).
[14] Stpierre, C.C., and Bankoff, S.G., “Vapor Volume Profiles in Developing Two-phase Flow,” International Journal of Heat and Mass Transfer, Vol. 10, pp. 237–249, (1967).
[15] Sekoguchi, K., Tanaka, O., Esaki, S., Katsuki, N., and Nakasatomi, M., “Prediction Method of Flow Patterns in Subcooled and Low Quality Boiling Regions,” Bulletin of the Japan Society of Mechanical Engineers. Eng, Vol. 24, pp. 834–841, (1981).
[16] Bartel, M.D., “Experimental Investigation of Subcooled Boiling,” M.Ss. Thesis in Nuclear Engineering, Purdue University, West Lafayette, IN, USA, (1999).
[17] Chen, E., Li, Y., and Cheng, X., “CFD Simulation of Upward Subcooled Boiling Flow of Refrigerant-113 using the Two-fluid Model,” Applied Thermal Engineering, Vol. 29, pp.  2508–2517, (2009).
[18] Ishii, M., and Mishima, K., ‘‘Two-fluid Model and Hydrodynamic Constitutive Relations,’’ Nuclear Engineering and Design, Vol. 82, pp. 107–126, (1984).
[19] ANSYS-CFX Development Team, “Solver Theory, Multiphase Flow Theory,” CFX-12 Documentation, ANSYS, (2009).
[20] Lemmert, M., and Chawla, J.M., “Influence of Flow Velocity on Surface Boiling Heat Transfer Coefficient,” in: E. Hahne, U. Grigull (Eds.), Heat Transfer in Boiling, Academic Press and Hemisphere, (1977).
[21] Ivey, H.J., “Relationships Between Bubble Frequency, Departure Diameter and Rise Velocity in Nucleate Boiling,” International Journal of Heat and Mass Transfer, Vol. 10, pp. 1023–1040, (1967).
[22] Tolubinski, V. I., and Kostanchuk, D. M., “Vapour Bubbles Growth Rate and Heat Transfer Intensity at Subcooled Water Boiling”, 4th International Heat Transfer Conference, Paris, France, (1970).
[23] Ranz, W. E., and Marshall, W. R., ‘‘Evaporation from Drops,’’ Chemical Engineering Progress, Vol. 48, pp. 142–180, (1952).       
[24]      Ishii, M., and Zuber, N., “Drag Coefficient and Relative Velocity in Bubbly, Droplet or Particulate Flows”, AIChE ( The American Institute of Chemical Engineering) Journal, Vol. 25, pp. 843-855, (1979).
[25] Tomiyama, A., “Struggle with Computational Bubble Dynamics”, 3rd International Conference on Multiphase Flow, Lyon, France, pp. 1-18, (1998).
[26] Antal, S. P., Lahey, R. T., and Flaherty, J. E., “Analysis of Phase Distribution in Fully Developed Laminar Bubbly Two-phase Flow”, International Journal of Multiphase Flow, Vol. 7, pp. 635-652, (1991).
[27] Rhie, C.M. and Chow, W.L., “A Numerical Study of the Turbulent Flow Past an Isolated Airfoil with Trailing Edge Separation”, AIAA J.,  pp. 82-0998, (1982).
[28] Kang, S., and Roy, R.P., “Vapor Phase Measurements in Subcooled Boiling Flow,” Journal of Heat Transfer, Vol. 124, pp. 1207–1209, (2002).
[29] Koncar, B., and Krepper, E., “CFD Simulation of Convective Flow Boiling of Refrigerant in a Vertical Annulus,” Nuclear Engineering and Design, Vol. 238, pp. 693–706, (2008).
[30] Sato, Y., Sadatomi, M., and Sekoguchi, K., “Momentum and Heat Transfer in Two-phase Bubble Flow,” International Journal of Multiphase Flow, Vol. 7, pp. 167–177, (1981).
[31] Lucas, D., Krepper, E., and Prasser, H.M., “Prediction of Radial Gas Profiles in Vertical Pipe Flow on Basis of the Bubble Size Distribution,” International Journal of Thermal Science, Vol. 40, pp. 217–225, (2001).
[32] Bartel, M.D., Ishii, M., Masukawa, T., Mi, Y., and Situ, R., “Interfacial Area Measurements in Subcooled Flow Boiling,” Nuclear Engineering and Design, Vol. 210, pp. 135–155, (2001).