Abstract

Thermosiphon is a kind of heat transfer equipment with high thermal conductivity. It can transfer heat for long distance by boiling and condensing. Thus, it is important to study the heat transfer character of the boiling process in the evaporator. In this paper, we have compared mainly seven different heat transfer formulas by other researches and our experiment results. Such formulas include Rohsenow formula (1951), Cooper formula (1984), Labuntsov formula (1973), Kutateladeze formula (1990), Kruzhilin formula (1947), Imura formula (1979) and Hasna formula (2017). Rohsenow formula is derived in theoretical for nucleate boiling process in pool boiling. Cooper formula is mainly for refrigerating coolants. Imura formula is an empirical formula based on the experiment result of the two phase closed thermosiphon. Hasna formula has the same pattern as the Cooper formula but it is used for separated heat pipe system. Labuntsov formula, Kutateladeze formula and Kruzhilin formula are all for the nucleate boiling. Four coolants are used for comparison including R744, R134a, R22 and R410A. The results show that, Imura formula’s calculating results are lower than the experiment results. Cooper formula accords very well with the experiment results for Freon coolant, but for R744, the relative error is as large as 62%. Rosenow formula’s average error is ~29%, and the variation trend is close to the experiment results.

1 INTRODUCTION

The separated heat pipe system which can also be called loop thermosiphon is a kind of heat transfer equipment. There are mainly two kinds of thermosiphon, as can be seen in Figure 1. The first kind is the separated heat pipe system and the second is the two phase closed thermosiphon. Thermosiphon mainly contains four parts: evaporator, condenser, vapor pipe and liquid pipe. Inside the evaporator, coolant absorbs heat and changes its phase from liquid to vapor, then vapor comes to the condenser through vapor pipe, in the condenser, vapor releases heat to heat sink and becomes liquid. Finally, liquid coolant flows back to the evaporator through liquid pipe. The condenser is placed higher than the evaporator, such system can be driven by gravity instead of pump.

The scheme of the separated heat pipe system and two phase closed thermosiphon, adapted from reference [1].
Figure 1.

The scheme of the separated heat pipe system and two phase closed thermosiphon, adapted from reference [1].

Although these two kinds of thermosiphon are all driven by gravity, there is a little difference in the flow patterns. For two phase closed thermosiphon, the vapor and liquid flows in the opposite direction while for the separated heat pipe system, the flowing direction is the same. In separated heat pipe system, which is also called loop heat pipe system, the coolant flows as a circle. As a result, from evaporator to condenser and from condenser to the evaporator is the same direction, also it can be regarded as the counterclockwise circulation.

Thus, the boiling pattern in the two phase closed thermosiphon is more likely to be the pool boiling, while the separated heat pipe system is the combination of pool boiling and convection boiling. As a result, it is important to study the boiling heat transfer coefficient for the separated heat pipe system.

There are four kinds of coolants used as the working fluid for the heat pipe system in our paper, including: R744, R134a, R410A and R22. Among these four kinds of coolants, R744 (carbon dioxide) is an environmental friendly coolant, however as the working pressure is much higher than Freon, it is not used generally, R22 is a kind of traditional coolant which will be obsoleted in the future because it will threaten the ozone sphere. R134a and R410A are all environmental friendly coolant. Detailed information of different coolants can be seen in Table 1.

Table 1.

Detail information of different coolants.

CoolantODPP(MPa)@20°CCritical pressureRelative molecular weight
R74405.737.3844.0
R220.0450.914.97486.5
R134a00.574.067102.0
R410A01.454.9586.0
CoolantODPP(MPa)@20°CCritical pressureRelative molecular weight
R74405.737.3844.0
R220.0450.914.97486.5
R134a00.574.067102.0
R410A01.454.9586.0

Note 1: R410A is a kind of non-azeotropic refrigerants, thus the critical pressure is estimated in the T-S diagram.

Note 2: GWP means global warming potential and ODP means Ozone Depression Potential. ODP larger than 0 means such coolant can destroy the ozone layer. The lower the GWP, the more the environmental friendly.

Note 3: Since the GWP of each coolant in different references have different values, we have not listed the GWP in the table. However, the GWP of R744 is much lower than other Freon coolants.

Table 1.

Detail information of different coolants.

CoolantODPP(MPa)@20°CCritical pressureRelative molecular weight
R74405.737.3844.0
R220.0450.914.97486.5
R134a00.574.067102.0
R410A01.454.9586.0
CoolantODPP(MPa)@20°CCritical pressureRelative molecular weight
R74405.737.3844.0
R220.0450.914.97486.5
R134a00.574.067102.0
R410A01.454.9586.0

Note 1: R410A is a kind of non-azeotropic refrigerants, thus the critical pressure is estimated in the T-S diagram.

Note 2: GWP means global warming potential and ODP means Ozone Depression Potential. ODP larger than 0 means such coolant can destroy the ozone layer. The lower the GWP, the more the environmental friendly.

Note 3: Since the GWP of each coolant in different references have different values, we have not listed the GWP in the table. However, the GWP of R744 is much lower than other Freon coolants.

In order to compare different empirical formulas, we have done an experimental work. In the experiment, the boiling heat transfer coefficient in the separated heat pipe system has been studied. The heating power is ~100–700 W, the tube inside diameter is 8 mm, the length of the evaporator is 1.515 m. For these four coolants, the heat transfer coefficient of R134a has the lowest uncertainty (3.48–8.55%), while R744 has the highest (10.04–25.18%). The experiment details can be seen in reference [1].

2 INTRODUCTION TO DIFFERENT EMPIRICAL FORMULAS AND CALCULATION ERROR ANALYSIS

2.1 Introduction to different empirical formulas

The paper has compared mainly seven different heat transfer formulas and our experiment results. Such formulas include Rohsenow [2] formula (1951), Cooper [3] formula (1984), Labuntsov [4, 5] formula (1973), Kutateladeze [5, 6] formula (1990), Kruzhilin [5, 7] formula (1947), Imura [8] formula (1979) and Hasna [9] formula (2017). Rohsenow formula is derived in theoretical for nucleate boiling process in pool boiling.

Rohsenow’s formula is derived in theoretical [2, 10]. In his theory, the boiling process can be regarded as a special kind of liquid heat transfer process. As the heat is mainly transferred directly from heating surface to liquid and bubbles may enhanced the heat transfer process. Thus, the Nusselt number (Nu) can be written as equation (1).
(1)
Here, the diameter and velocity of the fluid needs to know before calculating the Renold number (Re). According to Rohsenow’s [2] theory, the bubble will leave the heating surface when the surface tension (Ft) is equal to buoyancy force (Ff).
(2)
Thus, the leaving diameter is Db in equation (3), and V is the characteristic velocity in Re number.
(3)
Now, substitute Db and V in equation (1), we can get equation (4), which is the Rohsenow’s formula:
(4)
Cooper [3] has studied the pool boiling heat transfer coefficient for refrigerating coolants and has given his formula, see equation (5). pr means working pressure/critical pressure. M is the Relative molecular weight. Detailed information can be seen in Table 1. ε is the surface roughness, for common pipe, ε is ~0.3–0.4, in our calculation, as the tube surface roughness is unknown, thus the lgε is omitted. Cooper’s formula is written as equation (6). Such ε can be omitted because the heat transfer coefficient results calculated by equations (5 and 6) are similar. For example, for R744, when heating power is 109 W, the relative difference is ~5%, when heating power is 696 W, the relative difference is ~2%. But for R134a, the relative difference is higher which is ~16–19%.
(5)
(6)
Hasna [9] has obtained his formula based on the cooper’s formula. Thus, Hasna’s formula has the same pattern as the Cooper formula but it is used for separated heat pipe system. In Hasna’s experiment, the working fluid is water and the heat flux changes from ~2 W/cm2 to 20 W/cm2.
(7)
Imura [8] formula is an empirical formula based on the experiment results of the two phased closed thermosiphon. For the two phase closed thermosiphon, as the vapor flows upward and liquid flows downward, the boiling pattern inside the evaporator is more likely to be the pool boiling. In his formula, the relationship between the h and thermophysical property is considered. In equation (8), subscript l means liquid, g means vapor, cp is the specific heat at constant pressure, q is the heat flux, r is the latent heat, μ is viscosity, g is the gravity accelerate, ρ is density, k is the thermal conductivity.
(8)
Labuntsov formula [4, 5] (equation (9)), Kutateladeze formula [5, 6] (equation 10a–c) and Kruzhilin formula [5, 7] (equation (11)) are all for the nucleate boiling.
(9)
(10)
(11)

2.2 Error between experiment and formula analysis

In order to compare the error between our experiment results and the other researcher’s formulas, we have defined the following errors:

The first error is the calculation error, as can be seen in equation (12). Here, hcal means the calculation results by other researcher’s formulas and hexp is our experiment results.
(12)
The second error is the root-mean-square error. As for a given coolant, there are seven heating fluxes. For each condition, it has a calculation error, in order to know the average error of the seven conditions, the root-mean-square error can be used. In equation (13), N is the total calculation number.
(13)
In fact, the above mentioned two errors cannot represent the stability of the empirical formula. For example, for one formula, the calculation error is higher for some condition while it is lower at the other condition. But for another formula, at all condition, the error may roughly be the same. The root-mean-square error may be similar for these two formulas but in fact, the stability of the second empirical formula may be better than the first one. As a result, we use the standard deviation to evaluate the stability, as can be seen in equation (14). When calculation, we can use the function STDVP in excel to calculate the standard deviation. There are four coolants and each coolant has seven conditions, thus the standard deviation of one formula is for 28 conditions.
(14)

3 RESULTS AND ANALYSIS

According to the calculation results:

For R134a coolant, Imura’s formula is ~34% lower than the experiment results. Besides, with the increasing of the heating power, the error also comes higher. Rohsenow’s equation has the best results, the average error is just 9%. For R22, Cooper’s results are more close to the experiment, Kruzhilin’s variation trends are more similar but the results are lower than the experiment results.

For R410A, Imura’s formula is ~45% lower than the experiment results. The variation trend is just the similar with R134a. The variation trends of Kutateladze formula (average error: 31%), Labuntsov formula (average error: 37%), Kruzhilin formula (average error: 60%), and Cooper formula (average error: 25%) are similar, for higher heating power, the error comes lower. For Hasna’s formula, his experiment’s heating flux range is from ∼2 to 20 W/cm2, and our experiment’s heating flux range is from 0.26 to 1.84 W/cm2, thus, when the heating power is ~700 W, the Hasna’s formula may be more accurate to calculate our experiment condition, the results also shows that, the higher heating power, the lower the error between Hasna’s formula and our experiment results (Table 2).

Table 2.

The average relative error between formula and experiment results.

CoolantImuraKutateladzeLabuntsovRohsenowKruzhilinHasnaCooper
R134a0.340.280.480.090.590.240.28
R220.360.310.400.430.550.280.17
R410A0.450.310.370.320.600.250.25
R7440.390.620.560.220.580.990.62
Average0.380.410.460.290.580.540.37
CoolantImuraKutateladzeLabuntsovRohsenowKruzhilinHasnaCooper
R134a0.340.280.480.090.590.240.28
R220.360.310.400.430.550.280.17
R410A0.450.310.370.320.600.250.25
R7440.390.620.560.220.580.990.62
Average0.380.410.460.290.580.540.37
Table 2.

The average relative error between formula and experiment results.

CoolantImuraKutateladzeLabuntsovRohsenowKruzhilinHasnaCooper
R134a0.340.280.480.090.590.240.28
R220.360.310.400.430.550.280.17
R410A0.450.310.370.320.600.250.25
R7440.390.620.560.220.580.990.62
Average0.380.410.460.290.580.540.37
CoolantImuraKutateladzeLabuntsovRohsenowKruzhilinHasnaCooper
R134a0.340.280.480.090.590.240.28
R220.360.310.400.430.550.280.17
R410A0.450.310.370.320.600.250.25
R7440.390.620.560.220.580.990.62
Average0.380.410.460.290.580.540.37

For R744, Imura’s formula is also lower than our experiment results. And the errors are within 32–42%. Which means the errors of the seven condition are very stable. For Kutateladze formula, Labuntsov formula and Cooper formula, when the heating power is lower, the error is also not very large. However, when the heating power comes larger, the error also comes really large. Take the cooper’s formula for example, if the heating power is 109 W, the error between our experiment result and Cooper’s formula is ~9%, however, when the heating power is 696 W, the error is 114%. Here is the reason why at higher heating power, the formula’s error comes larger. As for these formula, they are mainly for pool boiling. However, for the boiling in the heat pipe system, it is more likely to be the combination of pool boiling and convection boiling. For heat pipe system, higher heating power also means higher flowing velocity of the liquid inside the evaporator, higher flow rate may enhance the heat transfer coefficient at first, however when the flow rate comes much higher, it may bother the bubble’s formation, thus, the heat transfer coefficient in the evaporator comes lower than pool boiling (Figure 2).

The results of calculation for different formulas and experiment.
Figure 2.

The results of calculation for different formulas and experiment.

Here, we can make a comparison between there seven formulas:

  • Imura’s formula:

    The calculation results for every coolants are all lower than experiment results. The average error is 38%, besides, the standard deviation is just as low as 0.08. According to the Imura’s equation, it considered the relation between heat transfer coefficient and many thermophysical properties of the coolants, which is very comprehensive. As for the reason why all the calculation results are all lower than the experiment result, it maybe the flow pattern. As for the combination of the flow and pool boiling, the heat transfer maybe enhanced by the coolant flow, however, when the flow rate comes much higher, the heat transfer maybe weaken. In fact, maybe the Imura’s equation does not consider the enhancing effect by the convection boiling. All in all, it is for two phase closed thermosiphon, the convection or flow factor is not considered.

  • Kutateladze’s formula, Labuntsov’s formula, Kruzhilin’s formula:

    All of these three formulas are for the pool boiling. The total average error of Kruzhilin is 58%. For Freon coolants, the errors of Kutateladze’s formula and Labuntsov’s formula are not very large, but for R744, for higher heating power, the error can even as high as 128% (Kutateladze, 696 W). In a word, these formulars are for pool boiling, which is not very suitable for our experiment thus, the error is also very big.

    The standard deviations of Kutateladze’s formula, Labuntsov’s formula, Kruzhilin’s formula are 0.39, 0.41 and 0.05 which means Kruzhilin’s formula has the best stability.

  • Cooper’s formula and Hasna’s formula:

    Cooper’s formula is a classical and widely used equation. Hasna’s formula is obtained based on the form of Cooper. Cooper’s formula is more suitable for Freon coolants. The standard deviation of Cooper for Freon is 0.1, but for all coolants it is 0.37. For Hasna’s equation, the higher the heating power, the lower the error. Besides, it is mainly for water in sub-atmospheric. Thus it is not suitable for R744.

  • Rohsenow’s formula:

    Rohsenow’s formula is derived in theoretical. For different coolant-heating surface, it has different Cwl and n. In our calculation, we let Cwl = 0.0049 and n = 1.7 (n-pentane-Lapped cooper). In fact, if there is the experiment results for cooper and each different coolant, the results can be more accurate. All in all, this formula is suitable for all the four coolants and the average error of all coolants is 29%. The standard deviations is 0.18.

4 CONCLUSION

At different working condition, different formulas show different results. These formulas include Rohsenow formula (1951), Cooper formula (1984), Labuntsov formula (1973), Kutateladeze formula (1990), Kruzhilin formula (1947), Imura formula (1979) and Hasna formula (2017). All in all, Imura, Rohsenow, Cooper and Hasna’s formula shows the better results.

Imura formula’s calculating results are all less than the experiment results. For example, it is 34% less than the experiment results for R134a while, for R744, it is 39% less than the experiment results, but the variation trend is very close to the experiment results. Cooper formula accords very well with the experiment results for Freon coolant, the relative average error is ~17–28%, but for R744, the relative error is as large as 62%. Rohsenow’s formula average error is ~29%, and the variation trend is close to the experiment results. For Hasna’s equation, when calculating the Freon coolants, the higher the heating power, the lower the error. However, it is not suitable for R744.

FUNDING

This paper is supported by National Key R&D Program of China (2016YFB0601600).

REFERENCES

1

Ding
T
,
Cao
HW
,
He
ZG
, et al.
Boiling heat transfer characteristics of the R744 coolant in the evaporator of the separated heat pipe system
.
Int. J. Heat Mass Transfer
2017
;
113
:
1254
64
.

2

Rohsenow
WM
. A method of correlating heat transfer data for surface boiling of liquids. Cambridge, MA: MIT Division of Industrial Cooporation,
1951
.

3

Cooper
MG
.
Saturation nucleate pool boiling—a simple correlation
.
Inst. Chem. Eng. Symp. Ser
.
1984
;
86
:
785
93
.

4

Labuntsov
DA
.
Heat transfer problems with nucleate boiling of liquids
.
Therm. Eng.(USSR)(Engl. Transl.)
1973
;
19
:
21
8
.

5

Pioro
IL
,
Rohsenow
W
,
Doerffer
SS
.
Nucleate pool-boiling heat transfer. II: assessment of prediction methods[J]
.
Int. J. Heat Mass Transfer
2004
;
47
:
5045
57
.

6

Kutateladze
SS
.
Heat Transfer and Hydrodynamic Resistance: Handbook
.
Energoatomizdat Publishing House
,
1990
. Chapter 12.7 (in Russian).

7

Kruzhilin
GN
.
Free-convection transfer of heat from a horizontal plate and boiling liquid
.
Doklady AN SSSR (Reports of the USSR Academy of Sciences)
1947
;
58
:
1657
60
(in Russian).

8

Imura
H
,
Kusuda
H
,
Ogata
JI
, et al.
Heat transfer in two-phase closed-type thermosyphons
.
JSME Trans
1979
;
45
:
712
22
.

9

Louahlia-Gualous
H
,
Le Masson
S
,
Chahed
A
.
An experimental study of evaporation and condensation heat transfer coefficients for looped thermosyphon
.
Appl Therm Eng
2017
;
110
:
931
40
.

10

Bergman
TL
,
Incropera
FP
,
DeWitt
DP
, et al.
Fundamentals of Heat and Mass Transfer
.
John Wiley & Sons
,
2011
.

APPENDIX

Nomenclature
NuNusselt number
ReRenold number
PrPrandtl number
ρDensity
DbThe bubble’s diameter
σSurface tension
hfg,rLatent heat
hHeat transfer coefficient
cpspecific heat at constant pressure
q,qhHeat flux
μViscosity
k,λHeat conductivity
prWorking pressure/critical pressure
MRelative molecular weight
PaAtmospheric pressure
lLiquid
vVapor
Nomenclature
NuNusselt number
ReRenold number
PrPrandtl number
ρDensity
DbThe bubble’s diameter
σSurface tension
hfg,rLatent heat
hHeat transfer coefficient
cpspecific heat at constant pressure
q,qhHeat flux
μViscosity
k,λHeat conductivity
prWorking pressure/critical pressure
MRelative molecular weight
PaAtmospheric pressure
lLiquid
vVapor
Nomenclature
NuNusselt number
ReRenold number
PrPrandtl number
ρDensity
DbThe bubble’s diameter
σSurface tension
hfg,rLatent heat
hHeat transfer coefficient
cpspecific heat at constant pressure
q,qhHeat flux
μViscosity
k,λHeat conductivity
prWorking pressure/critical pressure
MRelative molecular weight
PaAtmospheric pressure
lLiquid
vVapor
Nomenclature
NuNusselt number
ReRenold number
PrPrandtl number
ρDensity
DbThe bubble’s diameter
σSurface tension
hfg,rLatent heat
hHeat transfer coefficient
cpspecific heat at constant pressure
q,qhHeat flux
μViscosity
k,λHeat conductivity
prWorking pressure/critical pressure
MRelative molecular weight
PaAtmospheric pressure
lLiquid
vVapor

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial re-use, please contact [email protected]