Research on the effect of Reynolds correlation in natural-convection film condensation

Research on the effect of Reynolds correlation in natural-convection film condensation 1 first-author 1 1 corresp Corresponding author: cnst_xizhen@163.com Corresponding author: cnst_xizhen@163.com cnst_xizhen@163.com 1 1 2 Institute of Nuclear and New Energy Technology, Collaborative Innovation Center of Advanced Nuclear Energy Technology, Key Laboratory of Advanced Reactor Engineering and Safety of Ministry of Education, Tsinghua University, Beijing 100084, China Institute of Nuclear and New Energy Technology, Collaborative Innovation Center of Advanced Nuclear Energy Technology, Key Laboratory of Advanced Reactor Engineering and Safety of Ministry of Education, Tsinghua University, Beijing 100084, China Engineering Physics Institute, University of Wisconsin Madison, Madison 53706, USA Engineering Physics Institute, University of Wisconsin Madison, Madison 53706, USA Film condensation is a vital phenomenon in the nuclear engineering applications, such as the gas-steam pressurizer design, and heat removing on containment in the case of postulated accident. Reynolds number in film condensation can be calculated from either the mass relation or the energy relation, but few researches have distinguished the difference between them at present. This paper studies the effect of Reynolds correlation in the natural convection film condensation on the outer tube. The general forms of the heat transfer coefficient correlation of film condensation are developed in different flow regimes. By simultaneously solving a set of the heat transfer coefficient correlations with Re mass and Re energy , the general expressions for Re mass and Re energy and the relation between the corresponding heat transfer coefficients are obtained. In the laminar and wave-free flow regime, Re mass and Re energy are equivalent, while in the laminar and wavy flow regime, Re mass is much smaller than Re energy , and the deviation of the corresponding average heat transfer coefficients is about 30% at the maximum. In the turbulent flow regime, the relation of Re mass and Re energy is greatly influenced by Prandtl number. The relative deviation of their average heat transfer coefficients are the non-linear function of Reynolds number and Prandtl number. Compared with experimental results, the heat transfer coefficient calculated from Re energy is more accurate. Film condensation Reynolds correlation Heat transfer coefficient Natural convection 1 INTRODUCTION 1 Film condensation is a classical and complicated physical phenomenon, which widely exists in nuclear engineering, chemical engineering, aerospace engineering and many other fields. Many researchers have made efforts to get the accurate heat transfer coefficient correlations. The complexity in film condensation involves the following aspects: the non-equilibrium characteristics, the different modes of condensation, vapor and liquid in different flow regimes, various wall geometries and different kinds of vapors and liquids[ 1 1 ]. Originally, Nusselt developed a heat transfer coefficient correlation for laminar film condensation [ 2 2 ]. Then, Nusselt’s theory was extensively improved and extended by considering the film subcooling, the inertia and drag, and the shear stress [ 3 3 - 6 6 3-6 ]. As an important factor in the film condensation, suction effect was found to influence the heat transfer coefficient, as high as 20% at the maximum [ 7 7 ]. Recently, the film condensation heat transfer in the presence of different species of non-condensable gas has been extensively studied, where the boundary layer method and the diffusion layer model were adopted [ 8 8 - 11 11 8-11 ]. Furthermore, Wu et al. [ 12 12 ] improved the diffusion layer model based on real gas state equation, and expanded its application to the gas-steam pressurizer. In terms of the analytical model, Sin [ 13 13 ] developed the analytical model of film condensation in Passive Containment Cooling System(PCCS), and Mosayebidorcheh et al. [ 14 14 ] used the new hybrid method to seek the analytical solution of the steady state condensation film on the inclined rotating disk. Many researchers verified Nusselt’s theory of film condensation with the experimental results. However, large deviation occurred in the laminar and wavy flow regime and the turbulent flow regime. For this reason, many modified heat transfer coefficient correlations in film condensation were proposed [ 15 15 - 18 18 15-18 ]. Some semi-theoretical correlations for film condensation based on the experimental method were developed [ 19 19 , 20 20 ]. One of the classic cases was developed by Butterworth [ 21 21 , 22 22 ]. In the Butterworth’s correlation, the whole flow regime is divided into three according to the value of Reynolds number: the laminar and wave-free flow regime, the laminar and wavy flow regime, and the turbulent flow regime. The corresponding correlations of the heat transfer coefficient in different flow regimes are given based on massive experimental databases. Many experiments have also been conducted to investigate the film condensation on the horizontal and vertical plane, inside and outside the tube, in different flow regimes, on declined condensation plane with different inclination angles, in presence of non-condensable gas, and in different applications [ 23 23 - 27 27 23-27 ]. Most heat transfer coefficient correlations for the film condensation given by previous researchers are the functions of Reynolds number, or adding the modified term itself as a function of Reynolds number. Thus, the accurate calculation of Reynolds number is considered vital to obtain an accurate heat transfer coefficient. According to the basic governing equations of the film condensation, the expressions of the mass flow rate of the condensate based on the mass relation and the energy relation can be obtained respectively. They provided two approaches to calculate the heat transfer coefficient of the film condensation. However, few researches have distinguished the differences between them, and improper applications are usually found in public researches. In this paper, the effect of Reynolds correlation in the natural convection film condensation on outer tube is studied. Based on the general assumption of the modified term of film condensation heat transfer coefficient, the general forms of heat transfer coefficient are derived in different flow regimes. We also discover the difference between the two approaches to calculate the heat transfer coefficient of the film condensation corresponding to the different Reynolds correlations by using the analytical method and the iterative solution method. The general expressions for Re mass and Re energy and the relation between the corresponding heat transfer coefficients are obtained. Re mass and Re energy are graphically shown to see the difference. They are also used to correlate the experimental data from the open literature. It is found that the approach based on Re energy is more accurate. 2 FILM CONDENSATION THEORY 1 2.1 Fundamentals of Condensation Heat Transfer 2 Generally, for the film condensation outside the tube, Reynolds number of the liquid film outside tube in the experiment is originally defined as <math display="block" id="M1"><mrow><mtext>Re</mtext><mo>=</mo><mfrac><mrow><mn>4</mn><mtext>Γ</mtext></mrow><mrow><msub><mi>μ</mi><mtext>l</mtext></msub></mrow></mfrac><mo>.</mo></mrow></math> The mass rate of condensation Г is measurable in the experiment while calculated in the following theoretical analysis. Based on the Nusselt’s theory, the expression of the condensate mass flow rate can be [ 28 28 ] <math display="block" id="M2"><mrow><mtext>Γ</mtext><mo>=</mo><mfrac><mrow><msub><mi>ρ</mi><mtext>l</mtext></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>l</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>g</mtext></msub><mo stretchy="false">)</mo><mtext>g</mtext><msup><mi>δ</mi><mn>3</mn></msup></mrow><mrow><mn>3</mn><msub><mi>μ</mi><mtext>l</mtext></msub></mrow></mfrac><mo>.</mo></mrow></math> Substituting Eq.(2) into Eq.(1), and it gives <math display="block" id="M3"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>mass</mtext></mrow></msub><mo>=</mo><mfrac><mrow><mn>4</mn><msub><mi>ρ</mi><mtext>l</mtext></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>l</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>g</mtext></msub><mo stretchy="false">)</mo><mtext>g</mtext><msup><mi>δ</mi><mn>3</mn></msup></mrow><mrow><mn>3</mn><msubsup><mi>μ</mi><mtext>l</mtext><mtext>2</mtext></msubsup></mrow></mfrac><mo>.</mo></mrow></math> The mass flow rate of the film condensate has another expression based on the energy relation [ 29 29 ], namely <math display="block" id="M4"><mrow><mtext>Γ</mtext><mo>=</mo><mfrac><mrow><msub><mover accent="true"><mi>h</mi><mo>¯</mo></mover><mrow><mtext>film</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>T</mi><mrow><mtext>sat</mtext></mrow></msub><mo>−</mo><msub><mi>T</mi><mtext>w</mtext></msub><mo stretchy="false">)</mo><mi>L</mi></mrow><mrow><msubsup><mi>h</mi><mrow><msub><mrow/><mrow><mtext>fg</mtext></mrow></msub></mrow><mtext>"</mtext></msubsup></mrow></mfrac><mo>.</mo></mrow></math> By substituting Eq.(4) into Eq.(1), we obtain Reynolds correlation based on energy relation as <math display="block" id="M5"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>energy</mtext></mrow></msub><mo>=</mo><mfrac><mrow><mn>4</mn><msub><mover accent="true"><mi>h</mi><mo>¯</mo></mover><mrow><mtext>film</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>T</mi><mrow><mtext>sat</mtext></mrow></msub><mo>−</mo><msub><mi>T</mi><mtext>w</mtext></msub><mo stretchy="false">)</mo><mi>L</mi></mrow><mrow><msub><mi>μ</mi><mtext>f</mtext></msub><msubsup><mi>h</mi><mrow><mtext>fg</mtext></mrow><mtext>"</mtext></msubsup></mrow></mfrac><mo>.</mo></mrow></math> According to Nusselt’s model, the correlations of heat transfer coefficient in the laminar and wavy flow regime and the turbulent flow regime are developed by adding the modified term. In this paper, we take full account of the difference of variable coefficients in the modified term and the uncertainty of the demarcation points of the different flow regimes, and obtain the general expressions of the heat transfer coefficient in three flow-regimes film condensation model. The condensation film is divided into three regimes: the laminar and wave-free flow regime (0 < Re ≤ Re 1 ), the laminar and wavy flow regime ( Re 1 < Re ≤ Re 2 ), and the turbulent flow regime ( Re > Re 2 ), as shown in Fig.1 Fig.1 . The demarcation points of the film flow regimes ( Re 1 , Re 2 ) are variable. In the laminar and wavy flow regime the modified term ( T m ) is a function of Reynolds number ( Re ). In the turbulent flow regime, T m is a function of Reynolds number and Prandtl number ( Re , Pr ). The coefficients in the modified terms are semi-empirical parameters, and usually can be treated as variable. In the laminar and wavy-free flow regime (0 < Re ≤ Re 1 ), Nusselt’s equation is <math display="block" id="M6"><mrow><msub><mi>h</mi><mrow><mtext>film</mtext></mrow></msub><mo>=</mo><msup><mrow><mrow><mo>[</mo><mrow><mfrac><mrow><msub><mi>ρ</mi><mtext>l</mtext></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>l</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>g</mtext></msub><mo stretchy="false">)</mo><mtext>g</mtext><msubsup><mi>h</mi><mrow><msub><mrow/><mrow><mtext>fg</mtext></mrow></msub></mrow><mtext>"</mtext></msubsup><msubsup><mi>k</mi><mrow><msub><mrow/><mtext>l</mtext></msub></mrow><mtext>3</mtext></msubsup></mrow><mrow><mn>4</mn><msub><mi>μ</mi><mtext>l</mtext></msub><mi>L</mi><mo stretchy="false">(</mo><msub><mi>T</mi><mrow><mtext>sat</mtext></mrow></msub><mo>−</mo><msub><mi>T</mi><mtext>w</mtext></msub><mo stretchy="false">)</mo></mrow></mfrac></mrow><mo>]</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mn>4</mn></mrow></msup><mo>.</mo></mrow></math> In terms of the original Reynolds number, it can also be expressed as <math display="block" id="M7"><mrow><mfrac><mrow><msub><mi>h</mi><mrow><mtext>film</mtext></mrow></msub></mrow><mrow><msub><mi>k</mi><mtext>l</mtext></msub></mrow></mfrac><msup><mrow><mrow><mo>[</mo><mrow><mfrac><mrow><msubsup><mi>μ</mi><mtext>l</mtext><mtext>2</mtext></msubsup></mrow><mrow><msub><mi>ρ</mi><mtext>l</mtext></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>l</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>g</mtext></msub><mo stretchy="false">)</mo><mtext>g</mtext></mrow></mfrac></mrow><mo>]</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>=</mo><mn>1.10</mn><msup><mrow><mtext>Re</mtext></mrow><mrow><mo>−</mo><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>.</mo></mrow></math> In the laminar and wavy flow regime ( Re 1 < Re ≤ Re 2 ) and in the turbulent flow regime ( Re > Re 2 ), the modified terms are added to the basic equation of local coefficient of heat transfer in the laminar and wavy-free flow regime. The local heat transfer coefficients in all flow regimes are <math display="block" id="M8"><mrow><mfrac><mrow><msub><mi>h</mi><mrow><mtext>film</mtext></mrow></msub></mrow><mrow><msub><mi>k</mi><mtext>l</mtext></msub></mrow></mfrac><msup><mrow><mrow><mo>[</mo><mrow><mfrac><mrow><msubsup><mi>μ</mi><mtext>l</mtext><mtext>2</mtext></msubsup></mrow><mrow><msub><mi>ρ</mi><mtext>l</mtext></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>l</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>g</mtext></msub><mo stretchy="false">)</mo><mtext>g</mtext></mrow></mfrac></mrow><mo>]</mo></mrow></mrow><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>=</mo><mn>1.10</mn><msup><mrow><mtext>Re</mtext></mrow><mrow><mo>−</mo><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>⋅</mo><msub><mi>T</mi><mtext>m</mtext></msub><mo>,</mo></mrow></math> where the general expressions of T m are <math display="block" id="M9"><mrow><msub><mi>T</mi><mtext>m</mtext></msub><mo>=</mo><mrow><mo>{</mo><mrow><mtable><mtr><mtd><mtable columnalign="left"><mtr><mtd><mn>1</mn><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mtext> </mtext><mn>0</mn><mo><</mo><mtext>Re</mtext><mo>≤</mo><msub><mtext>Re</mtext><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>m</mi><mn>1</mn></msub><msup><mtext>Re</mtext><mrow><msub><mi>n</mi><mtext>1</mtext></msub></mrow></msup><mtable><mtr><mtd><mrow/></mtd><mtd><mrow><msub><mrow><mtext>Re</mtext></mrow><mn>1</mn></msub><mo><</mo><mtext>Re</mtext><mo>≤</mo><msub><mrow><mtext>Re</mtext></mrow><mn>2</mn></msub></mrow></mtd></mtr></mtable></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mrow><msub><mi>m</mi><mn>2</mn></msub><msup><mrow><mtext>Re</mtext></mrow><mrow><msub><mi>n</mi><mtext>2</mtext></msub></mrow></msup><msup><mrow><mi>Pr</mi></mrow><mi>s</mi></msup><mtable><mtr><mtd><mrow/></mtd><mtd><mrow><mtext>Re</mtext><mo>></mo><msub><mrow><mtext>Re</mtext></mrow><mn>2</mn></msub></mrow></mtd><mtd><mrow/></mtd></mtr></mtable></mrow></mtd></mtr></mtable></mrow></mrow><mo>.</mo></mrow></math> The average film condensation heat transfer coefficient is obtained by the integral average method. <math display="block" id="M10"><mrow><mfrac><mn>1</mn><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mtext>film</mtext></mrow></msub></mrow></mfrac><mo>=</mo><mfrac><mn>1</mn><mrow><mtext>Re</mtext></mrow></mfrac><mstyle displaystyle="true"><mrow><munderover><mo>∫</mo><mn>0</mn><mrow><mtext>Re</mtext></mrow></munderover><mrow><mfrac><mn>1</mn><mrow><msub><mi>h</mi><mrow><mtext>film</mtext></mrow></msub></mrow></mfrac></mrow></mrow></mstyle><mi>d</mi><mtext>Re</mtext><mo>.</mo></mrow></math> The average heat transfer coefficients in different flow regimes are http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10335737&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10335734&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10335737&type= where <math display="block" id="M11"><mrow><msub><mrow><mtext>Re</mtext></mrow><mtext>m</mtext></msub><mo>=</mo><mtext> </mtext><mfrac><mn>1</mn><mrow><mn>1.1</mn><msub><mi>m</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>1</mn></msub><mo stretchy="false">)</mo></mrow></mfrac><mo stretchy="false">(</mo><msubsup><mrow><mtext>Re</mtext></mrow><mn>2</mn><mrow><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>1</mn></msub><mo stretchy="false">)</mo></mrow></msubsup><mo>−</mo><msubsup><mrow><mtext>Re</mtext></mrow><mn>1</mn><mrow><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>1</mn></msub><mo stretchy="false">)</mo></mrow></msubsup><mo stretchy="false">)</mo><mo>+</mo><mn>0.68</mn><msubsup><mrow><mtext>Re</mtext></mrow><mn>1</mn><mrow><mfrac><mn>4</mn><mn>3</mn></mfrac></mrow></msubsup><mo>.</mo></mrow></math> The subcooling of the condensate and temperature jump across the film [ 28 28 ] should be taken into account when the latent heat, h fg ’ , is calculated, as given by <math display="block" id="M12"><mrow><msubsup><mi>h</mi><mrow><mtext>fg</mtext></mrow><mtext>"</mtext></msubsup><mo>=</mo><msub><mi>h</mi><mrow><mtext>fg</mtext></mrow></msub><mrow><mo>[</mo><mrow><mn>1</mn><mo>+</mo><mn>0.68</mn><mfrac><mrow><msub><mi>C</mi><mrow><mtext>pl</mtext></mrow></msub><mo stretchy="false">(</mo><msub><mi>T</mi><mrow><mtext>sat</mtext></mrow></msub><mo>−</mo><msub><mi>T</mi><mtext>w</mtext></msub><mo stretchy="false">)</mo></mrow><mrow><msub><mi>h</mi><mrow><mtext>fg</mtext></mrow></msub></mrow></mfrac></mrow><mo>]</mo></mrow><mo>.</mo></mrow></math> The temperature jump across the condensate is considered when calculating the temperature of liquid film [ 28 28 ]. <math display="block" id="M13"><mrow><msub><mi>T</mi><mrow><mtext>film</mtext></mrow></msub><mo>=</mo><msub><mi>T</mi><mtext>w</mtext></msub><mo>+</mo><mn>0.25</mn><mo stretchy="false">(</mo><msub><mi>T</mi><mrow><mtext>sat</mtext></mrow></msub><mo>−</mo><msub><mi>T</mi><mtext>w</mtext></msub><mo stretchy="false">)</mo><mo>.</mo></mrow></math> As seen from Eqs.(3) and (5), Re energy contains the average heat transfer coefficient <math id="M14"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mtext>film</mtext></mrow></msub></mrow></math> , whereas, Re mass is a function of the thickness of the film condensate, which corresponds to the local heat transfer coefficient h film . In consideration of Eqs.(7) and (11), it is indicated that it provides two approaches to calculate the heat transfer coefficient of the film condensation and Reynolds number, such as pair parameters of ( <math id="M15"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mtext>film</mtext></mrow></msub></mrow></math> , Re energy ) and ( h film , Re mass ). However，the fundamental difference between these two approaches still remains unclear when applied to calculate the film condensation heat transfer coefficient. 2.2 Theoretical analysis of Re mass and Re energy 2 Under the same thermal hydraulic parameter condition，these two approaches of calculating the film condensation heat transfer coefficient are used to obtain the general expressions for Re mass and Re energy in different flow regimes. In the laminar and wave-free flow regime, combining Eqs.(3), (5), (6) and (11), Re mass and Re energy are equivalent and given as <math display="block" id="M16"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>mass</mtext></mrow></msub><mo>=</mo><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>energy</mtext></mrow></msub><mo>=</mo><mn>3.776</mn><msup><mrow><mo stretchy="false">[</mo><mfrac><mrow><msub><mi>ρ</mi><mtext>l</mtext></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>l</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>g</mtext></msub><mo stretchy="false">)</mo><mtext>g</mtext><msubsup><mi>k</mi><mtext>l</mtext><mtext>3</mtext></msubsup><msup><mi>L</mi><mn>3</mn></msup><msup><mrow><mo stretchy="false">(</mo><msub><mi>T</mi><mrow><mtext>sat</mtext></mrow></msub><mo>−</mo><msub><mi>T</mi><mtext>w</mtext></msub><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow><mrow><msubsup><mi>μ</mi><mtext>l</mtext><mtext>5</mtext></msubsup><msubsup><mi>h</mi><mrow><mtext>fg</mtext></mrow><mtext>"</mtext></msubsup><msup><mrow/><mtext>3</mtext></msup></mrow></mfrac><mo stretchy="false">]</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>4</mn></mrow></msup><mo>.</mo></mrow></math> In laminar and wavy flow regime, the similar approach in the laminar and wavy-free flow regime is also presented. Eqs.(3), (7) and (9) may be combined to provide <math display="block" id="M17"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>mass</mtext></mrow></msub><mo>=</mo><msub><mi>C</mi><mrow><mn>11</mn></mrow></msub><msup><mrow><mo stretchy="false">(</mo><mfrac><mrow><msub><mi>ρ</mi><mtext>l</mtext></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>l</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mi>g</mi></msub><mo stretchy="false">)</mo><mtext>g</mtext><msubsup><mi>k</mi><mtext>l</mtext><mtext>3</mtext></msubsup><msup><mi>L</mi><mn>3</mn></msup><msup><mrow><mo stretchy="false">(</mo><msub><mi>T</mi><mrow><mtext>sat</mtext></mrow></msub><mo>−</mo><msub><mi>T</mi><mtext>w</mtext></msub><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow><mrow><msubsup><mi>μ</mi><mtext>l</mtext><mtext>5</mtext></msubsup><msubsup><mi>h</mi><mrow><mtext>fg</mtext></mrow><mtext>"</mtext></msubsup><msup><mrow/><mtext>3</mtext></msup></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><mfrac><mn>1</mn><mrow><mn>4</mn><mo stretchy="false">(</mo><mn>3</mn><msub><mi>n</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn><mi>）</mi></mrow></mfrac></mrow></msup><mo>.</mo></mrow></math> Inserting Eqs.(5) and (11), it provides an expression for Re energy . <math display="block" id="M18"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>energy</mtext></mrow></msub><mo>=</mo><msup><mrow><mo stretchy="false">[</mo><mfrac><mn>4</mn><mrow><msub><mi>C</mi><mrow><mn>13</mn></mrow></msub></mrow></mfrac><msup><mrow><mo stretchy="false">(</mo><mfrac><mrow><msub><mi>ρ</mi><mtext>l</mtext></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>l</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>g</mtext></msub><mo stretchy="false">)</mo><mtext>g</mtext><msup><mi>L</mi><mn>3</mn></msup><msubsup><mi>k</mi><mrow><msub><mrow/><mtext>l</mtext></msub></mrow><mtext>3</mtext></msubsup><msup><mrow><mo stretchy="false">(</mo><msub><mi>T</mi><mrow><mtext>sat</mtext></mrow></msub><mo>−</mo><msub><mi>T</mi><mtext>w</mtext></msub><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow><mrow><msubsup><mi>μ</mi><mtext>l</mtext><mtext>5</mtext></msubsup><msubsup><mi>h</mi><mrow><mtext>fg</mtext></mrow><mtext>"</mtext></msubsup><msup><mrow/><mtext>3</mtext></msup></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>+</mo><mfrac><mrow><msub><mi>C</mi><mrow><mn>12</mn></mrow></msub></mrow><mrow><msub><mi>C</mi><mrow><mn>13</mn></mrow></msub></mrow></mfrac><mo stretchy="false">]</mo></mrow><mrow><mfrac><mn>1</mn><mrow><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></mfrac></mrow></msup><mo>,</mo></mrow></math> where <math id="M19"><mrow><msub><mi>C</mi><mrow><mn>11</mn></mrow></msub><mtext>=(</mtext><mfrac><mn>4</mn><mrow><mn>3</mn><msubsup><mi>m</mi><mn>1</mn><mn>3</mn></msubsup></mrow></mfrac><msup><mo stretchy="false">)</mo><mrow><mfrac><mn>1</mn><mrow><mn>3</mn><msub><mi>n</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn></mrow></mfrac></mrow></msup><msup><mn>4</mn><mrow><mfrac><mn>3</mn><mrow><mn>4</mn><mo stretchy="false">(</mo><mn>3</mn><msub><mi>n</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mfrac></mrow></msup><mo>,</mo><msub><mi>C</mi><mrow><mn>12</mn></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1.1</mn><msub><mi>m</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>1</mn></msub><mo stretchy="false">)</mo></mrow></mfrac><mo>,</mo><msub><mi>C</mi><mrow><mn>13</mn></mrow></msub><mo>=</mo><mn>0.68</mn><msubsup><mrow><mtext>Re</mtext></mrow><mn>1</mn><mrow><mfrac><mn>4</mn><mn>3</mn></mfrac></mrow></msubsup><mo>−</mo><mfrac><mrow><msubsup><mrow><mtext>Re</mtext></mrow><mn>1</mn><mrow><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>1</mn></msub><mo stretchy="false">)</mo></mrow></msubsup></mrow><mrow><mn>1.1</mn><msub><mi>m</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>1</mn></msub><mo stretchy="false">)</mo></mrow></mfrac><mo>.</mo></mrow></math> Substituting Eq.(17) into Eq.(16), the relation between Re mass and Re energy in the laminar and wavy flow regime would be given as <math display="block" id="M20"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>mass</mtext></mrow></msub><mo>=</mo><msub><mi>C</mi><mrow><mn>11</mn></mrow></msub><msup><mrow><mo stretchy="false">(</mo><mfrac><mrow><msub><mi>C</mi><mrow><mn>12</mn></mrow></msub></mrow><mn>4</mn></mfrac><msubsup><mrow><mtext>Re</mtext></mrow><mrow><mi mathvariant="normal">energy</mi></mrow><mrow><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>1</mn></msub><mo stretchy="false">)</mo></mrow></msubsup><mo>−</mo><mfrac><mrow><msub><mi>C</mi><mrow><mn>13</mn></mrow></msub></mrow><mn>4</mn></mfrac><mo stretchy="false">)</mo></mrow><mrow><mfrac><mn>3</mn><mrow><mn>4</mn><mo stretchy="false">(</mo><mn>3</mn><msub><mi>n</mi><mn>1</mn></msub><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mfrac></mrow></msup><mo>.</mo></mrow></math> In the turbulent flow regime, Re mass and Re energy are given as <math display="block" id="M21"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>mass</mtext></mrow></msub><mo>=</mo><msub><mi>C</mi><mrow><mn>21</mn></mrow></msub><msup><mrow><mo stretchy="false">[</mo><msup><mrow><mtext>Pr</mtext></mrow><mrow><mtext>-3s</mtext></mrow></msup><msup><mrow><mo stretchy="false">(</mo><mfrac><mrow><msub><mi>ρ</mi><mtext>l</mtext></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>l</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>g</mtext></msub><mo stretchy="false">)</mo><mtext>g</mtext><msubsup><mi>k</mi><mtext>l</mtext><mtext>3</mtext></msubsup><msup><mi>L</mi><mn>3</mn></msup><msup><mrow><mo stretchy="false">(</mo><msub><mi>T</mi><mrow><mtext>sat</mtext></mrow></msub><mo>−</mo><msub><mi>T</mi><mtext>w</mtext></msub><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow><mrow><msubsup><mi>μ</mi><mtext>l</mtext><mtext>5</mtext></msubsup><msubsup><mi>h</mi><mrow><mtext>fg</mtext></mrow><mtext>"</mtext></msubsup><msup><mrow/><mn>3</mn></msup></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><mfrac><mn>1</mn><mn>4</mn></mfrac></mrow></msup><mo stretchy="false">]</mo></mrow><mrow><mfrac><mn>1</mn><mrow><mn>3</mn><msub><mi>n</mi><mn>2</mn></msub><mo>+</mo><mn>1</mn></mrow></mfrac></mrow></msup><mo>,</mo></mrow></math> <math display="block" id="M22"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>energy</mtext></mrow></msub><mo>=</mo><msup><mrow><mo stretchy="false">[</mo><mfrac><mrow><mn>4</mn><msup><mrow><mtext>Pr</mtext></mrow><mi>s</mi></msup></mrow><mrow><msub><mi>C</mi><mrow><mn>21</mn></mrow></msub></mrow></mfrac><msup><mrow><mo stretchy="false">(</mo><mfrac><mrow><msub><mi>ρ</mi><mtext>l</mtext></msub><mo stretchy="false">(</mo><msub><mi>ρ</mi><mtext>l</mtext></msub><mo>−</mo><msub><mi>ρ</mi><mtext>g</mtext></msub><mo stretchy="false">)</mo><mtext>g</mtext><msup><mi>L</mi><mn>3</mn></msup><msubsup><mi>k</mi><mtext>l</mtext><mtext>3</mtext></msubsup><msup><mrow><mo stretchy="false">(</mo><msub><mi>T</mi><mrow><mtext>sat</mtext></mrow></msub><mo>−</mo><msub><mi>T</mi><mtext>w</mtext></msub><mo stretchy="false">)</mo></mrow><mn>3</mn></msup></mrow><mrow><msubsup><mi>μ</mi><mtext>l</mtext><mtext>5</mtext></msubsup><msubsup><mi>h</mi><mrow><mtext>fg</mtext></mrow><mtext>"</mtext></msubsup><msup><mrow/><mtext>3</mtext></msup></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><mn>1</mn><mo>/</mo><mn>3</mn></mrow></msup><mo>−</mo><mfrac><mrow><msub><mi>C</mi><mrow><mn>23</mn></mrow></msub></mrow><mrow><msub><mi>C</mi><mrow><mn>22</mn></mrow></msub></mrow></mfrac><msup><mrow><mtext>Pr</mtext></mrow><mtext>s</mtext></msup><mo>+</mo><msubsup><mrow><mtext>Re</mtext></mrow><mn>2</mn><mrow><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>2</mn></msub><mo stretchy="false">)</mo></mrow></msubsup><mo stretchy="false">]</mo></mrow><mrow><mfrac><mn>1</mn><mrow><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>2</mn></msub><mo stretchy="false">)</mo></mrow></mfrac></mrow></msup><mo>,</mo></mrow></math> Where <math id="M23"><mrow><msub><mi>C</mi><mrow><mn>21</mn></mrow></msub><mo>=</mo><msup><mrow><mtext>(</mtext><mfrac><mn>4</mn><mrow><mn>3</mn><msubsup><mi>m</mi><mn>2</mn><mn>3</mn></msubsup></mrow></mfrac><mo stretchy="false">)</mo></mrow><mrow><mfrac><mn>1</mn><mrow><mn>3</mn><msub><mi>n</mi><mn>2</mn></msub><mo>+</mo><mn>1</mn></mrow></mfrac></mrow></msup><msup><mn>4</mn><mrow><mfrac><mn>3</mn><mrow><mn>4</mn><mo stretchy="false">(</mo><mn>3</mn><msub><mi>n</mi><mn>2</mn></msub><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mfrac></mrow></msup><mo>,</mo><msub><mi>C</mi><mrow><mn>22</mn></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1.1</mn><msub><mi>m</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>2</mn></msub><mo stretchy="false">)</mo></mrow></mfrac><mo>,</mo><msub><mi>C</mi><mrow><mn>23</mn></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1.1</mn><msub><mi>m</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>1</mn></msub><mo stretchy="false">)</mo></mrow></mfrac><mo stretchy="false">(</mo><msubsup><mrow><mtext>Re</mtext></mrow><mn>2</mn><mrow><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></msubsup><mo>−</mo><msubsup><mrow><mtext>Re</mtext></mrow><mn>1</mn><mrow><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>1</mn></msub></mrow></msubsup><mo stretchy="false">)</mo><mo>+</mo><mn>0.68</mn><msubsup><mrow><mtext>Re</mtext></mrow><mn>1</mn><mrow><mfrac><mn>4</mn><mn>3</mn></mfrac></mrow></msubsup></mrow></math> By substituting Eq.(20) into Eq.(19), we can get the relation between Re mass and Re energy in the turbulent flow regime. <math display="block" id="M24"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>mass</mtext></mrow></msub><mo>=</mo><msub><mi>C</mi><mrow><mn>21</mn></mrow></msub><msup><mrow><mo stretchy="false">[</mo><mfrac><mrow><msub><mi>C</mi><mrow><mn>22</mn></mrow></msub></mrow><mn>4</mn></mfrac><msup><mrow><mtext>Pr</mtext></mrow><mrow><mtext>-5</mtext><mi>s</mi></mrow></msup><mo stretchy="false">(</mo><msubsup><mrow><mtext>Re</mtext></mrow><mrow><mtext>energy</mtext></mrow><mrow><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>2</mn></msub><mo stretchy="false">)</mo></mrow></msubsup><mo>−</mo><msubsup><mrow><mtext>Re</mtext></mrow><mn>2</mn><mrow><mo stretchy="false">(</mo><mfrac><mn>4</mn><mn>3</mn></mfrac><mo>−</mo><msub><mi>n</mi><mn>2</mn></msub><mo stretchy="false">)</mo></mrow></msubsup><mo stretchy="false">)</mo><mo>+</mo><mfrac><mrow><msub><mi>C</mi><mrow><mn>23</mn></mrow></msub></mrow><mn>4</mn></mfrac><msup><mrow><mtext>Pr</mtext></mrow><mrow><mtext>-4</mtext><mi>s</mi></mrow></msup><mo stretchy="false">]</mo></mrow><mrow><mfrac><mn>3</mn><mrow><mn>4</mn><mo stretchy="false">(</mo><mn>3</mn><msub><mi>n</mi><mn>2</mn></msub><mo>+</mo><mn>1</mn><mo stretchy="false">)</mo></mrow></mfrac></mrow></msup><mo>.</mo></mrow></math> 2.3 Relation of Re mass and Re energy in Butterworth’s case 2 For further analysis, the coefficients of the Butterworth’s correlation are adopted as an example for Eqs.(18) and (21) ( Table 1 Table 1 ). 1 1 Parameter 1 1 m 1 1 1 n 1 1 1 m 2 1 1 n 2 1 1 s 1 1 Re 1 1 1 Re 2 1 1 Value 1 1 0.6868 1 1 0.11 1 1 0.2091 1 1 7/12 1 1 1/2 1 1 30 1 1 1600 By substituting all these parameters into Eqs.(18) and (21), the relations between Re mass and Re energy in different flow regimes can be obtained. In the laminar and wave-free flow regime <math display="block" id="M25"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>mass</mtext></mrow></msub><mo>=</mo><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>energy</mtext></mrow></msub><mo>.</mo></mrow></math> In the laminar and wavy flow regime <math display="block" id="M26"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>mass</mtext></mrow></msub><mo>=</mo><msup><mrow><mo stretchy="false">(</mo><mn>7.12</mn><msubsup><mrow><mtext>Re</mtext></mrow><mrow><mtext>energy</mtext></mrow><mrow><mn>1.22</mn></mrow></msubsup><mo>−</mo><mn>34.29</mn><mo stretchy="false">)</mo></mrow><mrow><mn>0.564</mn></mrow></msup><mo>.</mo></mrow></math> In the turbulent flow regime <math display="block" id="M27"><mrow><msub><mrow><mtext>Re</mtext></mrow><mrow><mtext>mass</mtext></mrow></msub><mo>=</mo><mn>228.40</mn><msup><mrow><mo stretchy="false">(</mo><msubsup><mrow><mtext>Re</mtext></mrow><mrow><mtext>energy</mtext></mrow><mrow><mtext>0</mtext><mtext>.75</mtext></mrow></msubsup><msup><mrow><mtext>Pr</mtext></mrow><mrow><mo>−</mo><mn>2.5</mn></mrow></msup><mo>+</mo><mn>150.86</mn><msup><mrow><mtext>Pr</mtext></mrow><mrow><mo>−</mo><mn>2.0</mn></mrow></msup><mo>−</mo><mn>253</mn><msup><mrow><mtext>Pr</mtext></mrow><mrow><mo>−</mo><mn>2.5</mn></mrow></msup><mo stretchy="false">)</mo></mrow><mrow><mn>0.2727</mn></mrow></msup><mo>.</mo></mrow></math> Under the same thermal hydraulic parameter condition, the relation of Re mass and Re energy is significantly different in each flow regime. Obviously, R e mass and Re energy are equivalent in the laminar and wavy-free regime. However, Fig.2 Fig.2 shows the relation between Re mass and Re energy in the laminar and wavy flow regime and the turbulent flow regime, described by Eqs.(23) and (24). From Fig.2 (a) Fig.2 (a) , we can see there is a significant difference between Re mass and Re energy in the laminar and wavy regime. Under the same thermal hydraulic parameter condition, Re energy keeps higher than Re mass (maximally 4 times). Prandtl number does not play a role in this flow regime. In the turbulent flow regime, Prandtl number significantly influences the relation between Re mass and Re energy (see the line of Re mass ~ Re energy in Fig.2 (b) Fig.2 (b) ). Considering the corosspoint of Re mass = Re energy and Re mass ~ Re energy , there is an obvious demarcation point at Prandtl number, 0.35. When Prandtl number is less than 0.35, the lines of Re mass = Re energy and Re mass ~ Re energy intersect at one point, and Re mass < Re energy and Re mass > Re energy are presented on either side of the coresspoint. When Prandtl number is larger than 0.35, the line of Re mass ~ Re energy deviates heavenly from that of Re mass = Re energy as Prandtl number increasing, and Re mass < Re energy . Especially, when Prandtl number increases to about 0.8 or more, 30 < Re mass < 1600, corresponding to the laminar and wavy regime; whereas, Re energy > 1600, the corresponding flow regime is the turbulent flow regime. In conclusion, Prandtl number greatly influences the relation between Re mass and Re energy , and their large deviation occurs in the turbulent flow regime. We can also obtain the heat transfer coefficient and Reynolds number, such as ( h film , Re mass ) and ( <math id="M28"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mtext>film</mtext></mrow></msub></mrow></math> , Re energy ), in different flow regimes by using the iterative solution method. The iterative calculation flow chart is shown in Fig.3 Fig.3 . According to the result of h film and Re mass , we can obtain <math id="M29"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mtext>film</mtext></mrow></msub></mrow></math> by using the average method of Eq.(10). Define the relative deviation of the average heat transfer coefficient as: <math display="block" id="M30"><mrow><mi>ε</mi><mo>=</mo><mfrac><mrow><msub><mover accent="true"><mi>h</mi><mo>¯</mo></mover><mrow><mtext>mass</mtext></mrow></msub><mtext> -</mtext><msub><mover accent="true"><mi>h</mi><mo>¯</mo></mover><mrow><mtext>energy</mtext></mrow></msub></mrow><mrow><msub><mover accent="true"><mi>h</mi><mo>¯</mo></mover><mrow><mtext>energy</mtext></mrow></msub></mrow></mfrac><mo>×</mo><mtext>100%</mtext><mo>.</mo></mrow></math> This definition of relative deviation not only presents its relative magnitude, but also directly shows whether it is positive or negative. According to Eqs.(22) and (25), ε = 0 in the laminar and wave-free flow regime, and it means that <math id="M31"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mtext>film</mtext></mrow></msub></mrow></math> and <math id="M32"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mi>e</mi><mi>n</mi><mi>e</mi><mi>r</mi><mi>g</mi><mi>y</mi></mrow></msub></mrow></math> are equivalent. In the laminar and wavy flow regime, the value of ε is positive, as shown in Fig.4 (a) Fig.4 (a) . The maximum deviation is as high as 30%, corresponding to Re energy of 1600. It means that <math id="M33"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mi>m</mi><mi>a</mi><mi>s</mi><mi>s</mi></mrow></msub></mrow></math> > <math id="M34"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mi>e</mi><mi>n</mi><mi>e</mi><mi>r</mi><mi>g</mi><mi>y</mi></mrow></msub></mrow></math> exists in the whole laminar and wavy flow regime. Fig.4 Fig.4 (b) shows a complicated situation in the turbulent flow regime, and ε is a highly non-linear function of Reynolds number and Prandtl number. When Prandtl number is less than 0.35, there are two intersections of the curve of ε and the horizontal coordinate. The value of ε is alternately negative, positive and negative, with the increase of Reynolds number; <math id="M35"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mi>m</mi><mi>a</mi><mi>s</mi><mi>s</mi></mrow></msub></mrow></math> < <math id="M36"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mi>e</mi><mi>n</mi><mi>e</mi><mi>r</mi><mi>g</mi><mi>y</mi></mrow></msub></mrow></math> and <math id="M37"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mi>m</mi><mi>a</mi><mi>s</mi><mi>s</mi></mrow></msub></mrow></math> > <math id="M38"><mrow><msub><mrow><mover accent="true"><mi>h</mi><mo stretchy="true">¯</mo></mover></mrow><mrow><mi>e</mi><mi>n</mi><mi>e</mi><mi>r</mi><mi>g</mi><mi>y</mi></mrow></msub></mrow></math> alternately present in low, moderate and high Reynolds number regime. The minimum value of ε is larger than -60%, corresponding to Pr =0.1 and Re energy =1600. When Prandtl number is larger than 0.35, ε monotonously decreases from positive value to negative value as Reynolds number increasing, and there is only one crosspoint of ε =0. When Prandtl number increases from 0.35 to 4.0, the value of ε increases from 0 to 35% in the low Reynolds number region (1600 < Re < 4000), the maximum value of ε is less than 60%, and the Re energy value corresponding to the crosspoint of ε =0 decreases from 12000 to 4000. In conclusion, when 0.1 < Pr < 4.0, ε is generally within the range of -60% ~ +60%, and it is a highly non-linear function of Reynolds number and Prandtl number in the turbulent flow regime. 3 EXPERIMENT VALIDATION 1 Experimental results from Bum-Jin Chung [ 15 15 ] and G.M.Hebbard and W.L.Badger [ 18 18 ] are used to verify the accuracy of the two Reynolds correlations. The experiment set by B.J. Chung [ 15 15 ] focuses on the film wise and drop wise condensation of steam on short inclined plates. Table 2 Table 2 shows the experiment conditions. 1 1 Parameters in the experiment 1 1 Value 1 1 Species of gas 1 1 Vapor 1 1 Total pressure 1 1 Atmospheric pressure 1 1 Temperature difference between bulk and wall 1 1 10-35K 1 1 Inclined angle 1 1 90 o 1 1 Effective length 1 1 0.09m 1 1 Convection in bulk 1 1 Natural The result is shown in Figs.5 Figs.5 and 6 6 . The heat flux and Reynolds number vary according to the temperature difference between wall and bulk. The result of the Nusselt’s theory case is obtained from the basic Nusselt's theory in the laminar and wave-free flow regime (Eqs.(6) or (7)), and the modified term ( T m ) is not taken into account in all flow regimes. It is just a reference case. It shows that the heat flux calculated with Re energy has a good agreement with the experimental data, especially in high temperature difference region. The result obtained with Re mass over-predicted the test results at an average deviation of 11.0%. The heat flux calculated from Nusselt’s theory is lower than the experimental results, and the deviation increases as ∆ T wb increasing. However, when ∆ T wb is small (10K < ∆ T wb < 16K) in this experiment, Nusselt’s theory can effectively predict the heat flux, while the method using Re energy slightly over-predicted the results. Major reason is that Reynolds number corresponding to this condition is less than 80 ( Fig.6 Fig.6 ), and Nusselt’s theory can be effective in the wavy-weak flow regime. Re mass and Re energy at the same ∆ T wb are quite different, as shown in Fig.5 Fig.5 . The average deviation Reynolds number between Re energy and Nusselt’s theory is within 10%. It is consistent with the theoretical analysis in the part II-C. In conclusion, the calculating approach based on Re energy promises a better performance in the low Reynolds number regime. Another pure steam condensation experiment performed by G.M. Hebbard and W.L. Badger [ 15 15 ] is selected to assess the Reynolds correlation in the laminar and wavy flow regime and the turbulent flow regime. The experimental conditions are listed in Table 3 Table 3 . 1 1 Parameters in the experiment 1 1 Value 1 1 species of gas 1 1 Vapor 1 1 Steam temperature 1 1 333.15-393.15K 1 1 Temperature difference between bulk and wall 1 1 5-25K 1 1 Inclined angle 1 1 90 o 1 1 Effective Length of tube 1 1 3.657m 1 1 Diameter of tube 1 1 0.0254m When the bulk temperature varies from 353.15 K to 393.15 K, Reynolds number is generally located in the laminar and wavy flow regime and the turbulent flow regime. According to the value of Reynolds number, different correlations are used for various flow regimes to determine the film condensation heat transfer coefficient. The result is shown in Fig. 7 Fig. 7 . Form Fig.7 Fig.7 , in the case of different steam temperature, the value curve of the heat transfer coefficient with Re energy is always located between the value curve of the heat transfer coefficient with Re mass and that of the heat transfer coefficient with Nusselt’s theory. Furthermore, the heat transfer coefficient with Re mass is maximum. The experimental data is also sandwiched between the value curves of the heat transfer coefficient with Re mass and Nusselt’s theory, and accords well with the heat transfer coefficient with Re energy . It is indicated that when the steam temperature changes, the calculations with Re energy always show a great accuracy. However, the heat coefficients are over-predicted by Re mass and underestimated by Nusselt’s theory. It is necessary to point out that when ∆ T wb is small (3K < ∆ T wb < 7K), the predictions based on Re mass fit the experimental data better, whereas when ∆ T wb is higher (7K < ∆ T wb < 25K), namely in the turbulent flow regime, the results based on Re energy are more accurate. In this experiment, the pure steam in the bulk is not solely under natural convection, but slowly driven by external force. This is the main reason why the predictions based on Re mass fit the experimental data better when ∆ T wb is small. When ∆ T wb is small, the condensation film is thin, and the weak forced convection is still influential and enhances the heat transfer coefficient. When ∆ T wb is larger and the condensation film becomes thicker, this effect becomes insignificant as the forced convection is no longer dominant. Thus, when ∆ T wb is small, the experiment results are slightly higher than the heat transfer coefficient of natural convection, which is similar to the prediction of Re energy . But they are closer to the lager-predicted result calculated from Re mass . In the natural convection of film condensation, the prediction of Re energy is more accurate. In addition, there may be a gradual transition between Re mass and Re energy cases when ∆ T wb is small, and the experimental value locates in the transition regime. The validation of this gradual transition will be our further study in the future. According to the previous analysis, Reynolds correlation influences not only the heat transfer coefficient, but also the calculation of the value itself, as they are defined based on different principles, namely the mass relation and the energy relation. The approach based on Re energy is a better form of formula when calculating the condensation heat transfer coefficient. 4 CONCLUSION 1 This paper mainly discusses the effect of the Reynolds correlation on the prediction of heat transfer coefficients in the natural-convection film condensation in different flow regimes. For the natural convection film condensation on the outer wall of tube, the general forms of the heat transfer coefficient correlation are developed in different flow regimes, which can be applied in a wide range of conditions. The difference between the two approaches to calculating the heat transfer coefficient of the film condensation corresponding to Re mass and Re energy is distinguished by using the analytical method and the iterative solution method. The general expressions for Re mass and Re energy and the relation of the corresponding average heat transfer coefficients are also obtained for three different flow regimes. In the laminar and wave-free flow regime, Re mass and Re energy are equivalent. In the laminar and wavy flow regime, Re mass is much smaller than Re energy , and the maximal deviation of the average heat transfer coefficients from the two Reynolds correlations is as high as about 30%. In the turbulent flow regime, Prandtl number greatly influences the relation between Re mass and Re energy . Compared with the experimental results, the approach based on Re energy turns out to be a good choice to calculate the heat transfer coefficient. Reference J.W. Rose , Condensation heat transfer fundamentals . Chem Eng Rese Design , 1998 , 76 ( 2 ): 143 - 152 . DOI: 10.1205/026387698524712 http://doi.org/10.1205/026387698524712 Condensation heat transfer fundamentals W. Nusselt , The condensation of steam on cooled surfaces . Z Ver Dtsch Ing , 1916 , 60 : 541 - 546 . The condensation of steam on cooled surfaces W.M. Rohsenow W , Heat transfer and temperature distribution in laminar film condensation . Trans. Asme , 1956 , 78 : 1645 - 1648 . Heat transfer and temperature distribution in laminar film condensation E.M Sparrow , J.L. Gregg . A boundary layer treatment of laminar film condensation . J Heat Transfer , 1959 , 81 ( 1 ): 13 - 18 . A boundary layer treatment of laminar film condensation W.J. Minkowycz , E.M. Sparrow . Condensation heat transfer in the presence of noncondensables, interfacial resistance, superheating, variable properties, and diffusion . Intern J Heat Mass Trans , 1966 , 9 ( 10 ): 1125 - 1144 . DOI: 10.1016/0017-9310(66)90035-4 http://doi.org/10.1016/0017-9310(66)90035-4 Condensation heat transfer in the presence of noncondensables, interfacial resistance, superheating, variable properties, and diffusion E.M. Sparrow , S.H. Lin , Condensation heat transfer in the presence of a non-condensable gas . J Heat Trans , 1964 , 86 ( 3 ): 430 - 436 . DOI: 10.1115/1.3688714 http://doi.org/10.1115/1.3688714 Condensation heat transfer in the presence of a non-condensable gas L.E. Herranz , A. Campo , Adequacy of the heat-mass transfer analogy to simulate containment atmospheric cooling in the new generation of advanced nuclear reactors: experimental confirmation . Nucl technol , 2002 , 139 ( 3 ): 221 - 232 . http://www.ans.org/store/article-3315/ http://www.ans.org/store/article-3315/ Adequacy of the heat-mass transfer analogy to simulate containment atmospheric cooling in the new generation of advanced nuclear reactors: experimental confirmation J.C.Y. Koh , An integral treatment of two-phase boundary layer in film condensation . J Heat Trans , 1961 , 83 ( 3 ): 359 - 362 . DOI: 10.1115/1.3682285 http://doi.org/10.1115/1.3682285 An integral treatment of two-phase boundary layer in film condensation V.E. Denny , A.F. Mills , V.J. Jusionis , Laminar film condensation from a steam-air mixture undergoing forced flow down a vertical surface . J Heat Trans , 1971 , 93 ( 3 ): 297 - 304 . DOI: 10.1115/1.3449814 http://doi.org/10.1115/1.3449814 Laminar film condensation from a steam-air mixture undergoing forced flow down a vertical surface M.H. Kim , M.L. Corradini , Modeling of condensation heat transfer in a reactor containment . Nucl Eng Design , 1990 , 118 ( 2 ): 193 - 212 . DOI: 10.1016/0029-5493(90)90057-5 http://doi.org/10.1016/0029-5493(90)90057-5 Modeling of condensation heat transfer in a reactor containment P.F. Peterson , Theoretical basis for the Uchida correlation for condensation in reactor containments . Nucl Eng Design , 1996 , 162 ( 2 ): 301 - 306 . DOI: 10.1016/0029-5493(95)01125-0 http://doi.org/10.1016/0029-5493(95)01125-0 Theoretical basis for the Uchida correlation for condensation in reactor containments L. Wu , Y. Liu , H.J. Jia , Improved diffusion layer model based on real gas state equation for high pressure condensation . Annals of Nuclear Energy , 2015 , 85 , 444 - 452 . DOI: 10.1016/j.anucene.2015.05.006 http://doi.org/10.1016/j.anucene.2015.05.006 Improved diffusion layer model based on real gas state equation for high pressure condensation S. Kima , Y.G. Lee , D.W. Jerng , Laminar film condensation of saturated vapor on an isothermal vertical cylinder . International J Heat Mass Trans , 2015 , 83 : 545 - 551 . DOI: 10.1016/j.ijheatmasstransfer.2014.12.009 http://doi.org/10.1016/j.ijheatmasstransfer.2014.12.009 Laminar film condensation of saturated vapor on an isothermal vertical cylinder S. Mosayebidorcheh , T. Mosayebidorcheh , MM. Rashidi , Analytical solution of the steady state condensation film on the inclined rotating disk by a new hybrid method . Scientific Research and Essays , 2014 , 9 ( 12 ), 557 - 565 . DOI: 10.13140/2.1.4464.2241 http://doi.org/10.13140/2.1.4464.2241 Analytical solution of the steady state condensation film on the inclined rotating disk by a new hybrid method B.J. Chung , M.C. Kim , M. Ahmadinejad , Film-wise and drop-wise condensation of steam on short inclined plates . J Mech Sci Technol , 2008 , 22 ( 1 ): 127 - 133 . http://link.springer.com/article/10.1007/s12206-007-1015-8 http://link.springer.com/article/10.1007/s12206-007-1015-8 Film-wise and drop-wise condensation of steam on short inclined plates L. Slegers , R.A. Seban , Laminar film condensation of steam containing small concentrations of air . Intern Journal of Heat and Mass Transfer , 1970 , 13 ( 12 ): 1941 - 1947 . DOI: 10.1016/0017-9310(70)90094-3 http://doi.org/10.1016/0017-9310(70)90094-3 Laminar film condensation of steam containing small concentrations of air S.S. Kutateladze , I.I. Gogonin , Heat transfer in film condensation of slowly moving vapor . Intern J Heat Mass Trans , 1979 , 22 ( 12 ): 1593 - 1599 . DOI: 10.1016/0017-9310(79)90075-9 http://doi.org/10.1016/0017-9310(79)90075-9 Heat transfer in film condensation of slowly moving vapor G.M. Hebbard , W.L. Badger , Steam-Film Heat Transfer Coefficients for Vertical Tubes . Industrial & Engineering Chemistry , 1934 , 26 ( 4 ): 420 - 424 . DOI: 10.1021/ie50292a013 http://doi.org/10.1021/ie50292a013 Steam-Film Heat Transfer Coefficients for Vertical Tubes J.W. Rose , Fundamentals of condensation heat transfer-Laminar film condensation . JSME Intern J , 1988 , 31 : 357 - 375 . http://ci.nii.ac.jp/naid/110002494180/ http://ci.nii.ac.jp/naid/110002494180/ Fundamentals of condensation heat transfer-Laminar film condensation X.Y. Wei , X.D. Fang , R.R. Shi , A comparative study of heat transfer coefficients for film condensation . Energy Sci and Technol , 2012 , 3 ( 1 ): 1 - 9 . http://www.cscanada.net/index.php/est/article/view/2344 http://www.cscanada.net/index.php/est/article/view/2344 A comparative study of heat transfer coefficients for film condensation G. Neiils , S. Kliein , Heat transfer . Cambridge University Press (USA) , 2009 : 799 - 806 . Heat transfer H. Uchida , A. Oyama , Y. Togo , Evaluation of post-incident cooling systems of light water power reactors . Tokyo University ., 1964 . http://www.osti.gov/scitech/biblio/4023463 http://www.osti.gov/scitech/biblio/4023463 Evaluation of post-incident cooling systems of light water power reactors A.A. Dehbi , The effects of noncondensable gases on steam condensation under turbulent natural convection conditions . Massachusetts Institute of Technology , 1991 . http://dspace.mit.edu/handle/1721.1/13905 http://dspace.mit.edu/handle/1721.1/13905 The effects of noncondensable gases on steam condensation under turbulent natural convection conditions M.H. Anderson , L.E. Herranz , M.L. Corradini . Experimental analysis of heat transfer within the AP600 containment under postulated accident conditions . Nucl Eng Design , 1998 , 185 ( 2 ): 153 - 172 . DOI: 10.1016/S0029-5493(98)00232-5 http://doi.org/10.1016/S0029-5493(98)00232-5 Experimental analysis of heat transfer within the AP600 containment under postulated accident conditions J.W. Kim , Y.G. Lee , H.K. Ahn , et al ., Condensation heat transfer characteristic in the presence of noncondensable gas on natural convection at high pressure . Nucl Eng Design , 2009 , 239 ( 4 ): 688 - 698 . DOI: 10.1016/j.nucengdes.2008.12.011 http://doi.org/10.1016/j.nucengdes.2008.12.011 Condensation heat transfer characteristic in the presence of noncondensable gas on natural convection at high pressure R.I. Pashkevich , P.V. Muratov , Film condensation in a large diameter tube with upward steam flow . Intern J Heat Mass Trans , 2015 , 81 : 804 - 810 . DOI: 10.1016/j.ijheatmasstransfer.2014.11.001 http://doi.org/10.1016/j.ijheatmasstransfer.2014.11.001 Film condensation in a large diameter tube with upward steam flow M. Kondo ,  H. Nakamura ,  Y. Kukita , et al ., Primary-side two-phase flow and heat transfer characteristics of a horizontal-tube PCCS condenser . In: 14th International Conference on Nuclear Engineering , Miami . 2006 , 685 - 693 . DOI: 10.1115/ICONE14-89652 http://doi.org/10.1115/ICONE14-89652 Primary-side two-phase flow and heat transfer characteristics of a horizontal-tube PCCS condenser J.G. Collier , J.R. Thome , Convective boiling and condensation . Oxford University Press(UK) , 1994 : 445 - 450 . Convective boiling and condensation T.L. Bergman , A.S. Lavine , F.P. Incropera , et al ., Fundamentals of Heat and Mass Transfer , 7 th Edition, John Wiley & Sons, Inc. , New Jersey(USA) , 2011 : 675 - 679 . Fundamentals of Heat and Mass Transfer 10.1007/s41365-017-0240-9 This work was supported by Doctoral Fund of Ministry of Education of China (No. 20130002120016), the Science Found for Creative Groups of National Natural Science Foundation of China (No. 51321002), and the project of "Research and validation of key technology and device for NHR200-II nuclear heating reactor" 2016-11-09 2017-04-20 2017-04-21 2017-05-06 2017-06-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2017 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-017-0240-9 http://dx.doi.org/10.1007/s41365-017-0240-9 . 2017 NOMENCLATURE C variable coefficient C p constant pressure specific heat, J/(K·kg) G acceleration of gravity, m/s 2 h condensation heat transfer coefficient, W/(m 2 ·K) h fg latent heat of vaporization, J/kg h fg ’ modified latent heat of vaporization, J/kg <math id="M39"><mrow><mi>h</mi><mo stretchy="false">(</mo><mo>_</mo><mo stretchy="false">)</mo></mrow></math> average heat transfer coefficient, W/(m 2 ·K) k thermal conductivity, W/m·K L length, m m variable coefficient n variable coefficient Pr Prandtl number Re Reynolds number s variable coefficient T temperature, K Symbols μ viscosity, Pa·s Г mass flow rate of condensate, kg/s ε coefficient ρ density, kg/m 3 δ thickness of condensation film, m △ difference Subscripts 1 laminar and wavy flow regime 2 Turbulent flow regime b bulk energy based on energy relation f liquid phase g gas phase l liquid phase mass based on mass relation sat saturation w wall Nuclear Science and Techniques 06 28