Abstract
In this paper a numerical simulation of upward two-phase flow evaporation in a vertical tube has been studied by considering water as working fluid. To this end, the computational fluid dynamic simulations of this system are performed with heat and mass transfer mechanisms due to energy transfer during the phase change interaction near the heat transfer surface. The volume of fluid model in an available Eulerian-Eulerian approach based on finite volume method is utilized and the mass source term in conservation of mass equation is implemented using a user defined function. The characteristics of water flow boiling such as void fraction and heat transfer coefficient distribution are investigated. The main cause of fluctuations on heat transfer coefficient and volume fraction is velocity increment in the vapor phase rather than the liquid phase. The case study of this research including convective heat transfer coefficient and tube diameter are considered as a parametric study. The operating conditions are considered at high pressure in saturation temperature and the physical properties of water are determined by considering system’s inlet temperature and pressure in saturation conditions. Good agreement is achieved between the numerical and the experimental values of heat transfer coefficients.









Similar content being viewed by others
Abbreviations
- cl, cg :
-
Empirical coefficients used in Eq. 12, 13 (1/s).
- D:
-
Diameter (m).
- E:
-
Enthalpy (J/kg).
- \( {\overrightarrow{\mathrm{F}}}_{\upsigma} \) :
-
Volumetric surface tension force (N/m^3).
- \( \overrightarrow{\mathrm{g}} \) :
-
Gravitational acceleration (m/s^2).
- h:
-
Convective heat transfer coefficient (W/m^2-K).
- hlg :
-
Latent heat (J/kg).
- k:
-
Curvatures of phases.
- ke :
-
Effective thermal conductivity (W/m-K).
- L:
-
Length (m).
- P:
-
Pressure (Pa).
- \( \dot{\mathrm{Q}} \) :
-
Heat transfer rate source term (W/m^3).
- S:
-
Mass source term (kg/m^3-s).
- T:
-
Temperature (K).
- t:
-
Time (s).
- \( \overrightarrow{\mathrm{u}} \) :
-
Velocity field (m/s).
- z:
-
Axial location (m)
- σ:
-
Interfacial surface tension (N/m).
- α:
-
Volume fraction.
- ρ:
-
Density (kg/m^3).
- μ:
-
Viscosity (kg/m-s)
- l:
-
Related to the liquid phase.
- g:
-
Related to the vapor phase.
- sat:
-
Saturation.
References
Yang Z, Peng X, Ye P (2008) Numerical and experimental investigation of two phase flow during boiling in a coiled tube. Int J Heat Mass Transf 51:1003–1016
Zambrana J, Leo T, Perez-del-Notario P (2008) Vertical tube length calculation based on available heat transfer coefficient expressions for the subcooled flow boiling region. Appl Therm Eng 28:499–513
Fang C, David M, Rogacs A, Goodson K (2010) Volume of fluid simulation of boiling two-phase flow in a vapor-venting microchannel. Front Heat Mass Transfer 1
Wei J, Pan L, Chen D, Zhang H, Xu J, Huang Y (2011) Numerical simulation of bubble behaviors in subcooled flow boiling under swing motion. Nuclear Eng Des 241:2898–2908
Kouhikamali R, Hassanpour B, Javaherdeh K (2012) Numerical simulation of forced convective evaporation in thermal desalination units with vertical tubes. Desalin Water Treat
Park, Seouk I (2010) Numerical analysis for flow, heat and mass transfer in film flow along a vertical fluted tube. Int J Heat Mass Transf 53:309–319
Chen Q, Xu J, Sun D, Cao Z, Xie J, Xing F (2013) Numerical simulation of the modulated flow pattern for vertical upflows by the phase separation concept. Int J Multiphase Flow 56:105–118
Ma C, Bothe D (2013) Numerical modeling of thermocapillary two-phase flows with evaporation using a two-scalar approach for heat transfer. J Comput Phys 233:552–573
Ratkovich N, Majumder S, Bentzen T (2012) Empirical correlations and CFD simulations of vertical two-phase gas–liquid (Newtonian and non-Newtonian) slug flow compared against experimental data of void fraction. Chem Eng Res Des
Meng M, Yang Z, Duan Y, Chen Y (2013) Boiling flow of R141b in vertical and inclined Serpentine Tubes. Int J Heat Mass Transf 57:312–320
Liu Y, Cui J, Li W (2011) A two-phase, two-component model for vertical upward gas–liquid annular flow. Int J Heat Fluid Flow 32:796–804
Liu Y, Li W, Quan S (2011) A self-standing two-fluid CFD model for vertical upward two-phase annular flow. Nucl Eng Des 241:1636–1642
Vazquez-Ramirez E, Riesco-Avila J, Polley G (2013) Two-phase flow and heat transfer in horizontal tube bundles fitted with baffles of vertical cut. Appl Therm Eng 50:1274–1279
Maqbool M, Palm B, Khodabandeh R (2013) Investigation of two phase heat transfer and pressure drop of propane in a vertical circular minichannel. Exp Thermal Fluid Sci 46:120–130
Chiapero E, Fernandino M, Dorao C (2014) Experimental results on boiling heat transfer coefficient, frictional pressure drop and flow patterns for R134a at a saturation temperature of 34 °C. Int J Refrig 40:317–327
Tsui Y, Lin S (2016) Three-dimensional modeling of fluid dynamics and heat transfer for two-fluid or phase change flows. Int J Heat Mass Transf 93:337–348
Hardt S, Wondra F (2008) Evaporation model for interfacial flows based on a continuum-field representation of the source terms. J Comput Phys 227:5871–5895
Pan Z, Weibel J, Garimella S (2016) A saturated-interface-volume phase change model for simulating flow boiling. Int J Heat Mass Transf 93:945–956
Kafeel K, Turan A (2013) Axi-symmetric simulation of a two phase vertical thermosyphon using Eulerian two-fluid methodology. Heat Mass Transf 49:1089–1099
Lorstad D, Fuchs L (2004) High-order surface tension VOF-model for 3D bubble flows with high density ratio. J Comput Phys 200:153–176
Das A, Das P (2015) Modeling of liquid–vapor phase change using smoothed particle hydrodynamics. J Comput Phys 303:125–145
Collier J, Thome J (1994) Convective Boiling and Condensation. Clarendon Press, Oxford
Lakehal D, Meier M, Fulgosi M (2002) Interface tracking towards the direct simulation of heat and mass transfer in multiphase flows. Int J Heat Fluid Flow 23:242–257
Lorstad D, Francois M, Shyy W, Fuchs L (2004) Assessment of volume of fluid and immersed boundary methods for droplet computations. Int J Numer Methods Fluids 46:109–125
Benson, David J (2002) Volume of fluid interface reconstruction methods for multi-material problems. Appl Mech Rev 55(2):151–165
Brackbill J, Kothe D, Zemach C (1992) A continuum method for modeling surface tension. J Comput Phys 100:335–354
Just T, Kelm S (1986) Die Industry, pp 38–76
Issa R (1985) Solution of the implicitly discretized fluid flow equations by operator-splitting. J Comput Phys 62:40–65
Versteeg H, Malalasekera W (2007) An Introduction to Computational Fluid Dynamics, Second edition ed., Pearson Education Limited
Menter F, Kuntz M, Langtry R (2003) Ten Years of industrial Experience with the SST Turbulence Model. Turbul Heat Mass Transf 4:625–632
Mosaad M, Johannsen K (1989) Experimental Study of Steady-State Film Boiling Heat Transfer of Subcooled Water Flowing Upwards in a Vertical Tube. Exp Thermal Fluid Sci 2:477–493
Kawaguchi K, Okui K, Asai T, Hasegawa Y (2006) The heat transfer and pressure drop characteristics of the finned tube banks in forced convection (Effect of fin height on heat transfer characteristics). Heat Transfer-Asian Res 35(3):194–208
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Hasanpour, B., Irandoost, M.S., Hassani, M. et al. Numerical investigation of saturated upward flow boiling of water in a vertical tube using VOF model: effect of different boundary conditions. Heat Mass Transfer 54, 1925–1936 (2018). https://doi.org/10.1007/s00231-018-2289-3
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00231-018-2289-3