-
PDF
- Split View
-
Views
-
Cite
Cite
Tao Ding, Jian min Wang, Chun Ou yang, Han wen Cao, Zhi guang He, Zhen Li, Comparison work about different empirical formulas for the boiling heat transfer coefficient in separated heat pipe system, International Journal of Low-Carbon Technologies, Volume 14, Issue 2, June 2019, Pages 103–107, https://doi.org/10.1093/ijlct/ctz007
- Share Icon Share
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].](https://oup.silverchair-cdn.com/oup/backfile/Content_public/Journal/ijlct/14/2/10.1093_ijlct_ctz007/1/m_ctz007f01.jpeg?Expires=1747919422&Signature=pTqe2rC82qaNS-ldAxNK3IkNygRF2nG~DFUKDqUPQpSxzzat6Ns9dXZITUvS9kNsuShQkWCUkqWo0KJM67i7dGuUIDFFpk0hY-fKvrY3l13Gr0z1iCuT-bCF6oCUloxsB~EP6J7B7sOdRDdwomHJ635a2Y7XC57Zdi8S5GdQnDPZlYI7uAPOC0fHVcSdOuzIYD80tZiH4ZQxwYlxS-b54LGRcoMEI-PIk4hG12Z5~~pq8Dnnwvx5u1e7KTMFo0vZoutbgjAPhBpFLdmUV36GIwBNJhy4iyA5EvOrKlo--hkOOoMiazc1TPk0d9VAla6knRLa-8~yoNtMi19JH6aGVA__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
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.
Coolant . | ODP . | P(MPa)@20°C . | Critical pressure . | Relative molecular weight . |
---|---|---|---|---|
R744 | 0 | 5.73 | 7.38 | 44.0 |
R22 | 0.045 | 0.91 | 4.974 | 86.5 |
R134a | 0 | 0.57 | 4.067 | 102.0 |
R410A | 0 | 1.45 | 4.95 | 86.0 |
Coolant . | ODP . | P(MPa)@20°C . | Critical pressure . | Relative molecular weight . |
---|---|---|---|---|
R744 | 0 | 5.73 | 7.38 | 44.0 |
R22 | 0.045 | 0.91 | 4.974 | 86.5 |
R134a | 0 | 0.57 | 4.067 | 102.0 |
R410A | 0 | 1.45 | 4.95 | 86.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.
Coolant . | ODP . | P(MPa)@20°C . | Critical pressure . | Relative molecular weight . |
---|---|---|---|---|
R744 | 0 | 5.73 | 7.38 | 44.0 |
R22 | 0.045 | 0.91 | 4.974 | 86.5 |
R134a | 0 | 0.57 | 4.067 | 102.0 |
R410A | 0 | 1.45 | 4.95 | 86.0 |
Coolant . | ODP . | P(MPa)@20°C . | Critical pressure . | Relative molecular weight . |
---|---|---|---|---|
R744 | 0 | 5.73 | 7.38 | 44.0 |
R22 | 0.045 | 0.91 | 4.974 | 86.5 |
R134a | 0 | 0.57 | 4.067 | 102.0 |
R410A | 0 | 1.45 | 4.95 | 86.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.
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:
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).
Coolant . | Imura . | Kutateladze . | Labuntsov . | Rohsenow . | Kruzhilin . | Hasna . | Cooper . |
---|---|---|---|---|---|---|---|
R134a | 0.34 | 0.28 | 0.48 | 0.09 | 0.59 | 0.24 | 0.28 |
R22 | 0.36 | 0.31 | 0.40 | 0.43 | 0.55 | 0.28 | 0.17 |
R410A | 0.45 | 0.31 | 0.37 | 0.32 | 0.60 | 0.25 | 0.25 |
R744 | 0.39 | 0.62 | 0.56 | 0.22 | 0.58 | 0.99 | 0.62 |
Average | 0.38 | 0.41 | 0.46 | 0.29 | 0.58 | 0.54 | 0.37 |
Coolant . | Imura . | Kutateladze . | Labuntsov . | Rohsenow . | Kruzhilin . | Hasna . | Cooper . |
---|---|---|---|---|---|---|---|
R134a | 0.34 | 0.28 | 0.48 | 0.09 | 0.59 | 0.24 | 0.28 |
R22 | 0.36 | 0.31 | 0.40 | 0.43 | 0.55 | 0.28 | 0.17 |
R410A | 0.45 | 0.31 | 0.37 | 0.32 | 0.60 | 0.25 | 0.25 |
R744 | 0.39 | 0.62 | 0.56 | 0.22 | 0.58 | 0.99 | 0.62 |
Average | 0.38 | 0.41 | 0.46 | 0.29 | 0.58 | 0.54 | 0.37 |
Coolant . | Imura . | Kutateladze . | Labuntsov . | Rohsenow . | Kruzhilin . | Hasna . | Cooper . |
---|---|---|---|---|---|---|---|
R134a | 0.34 | 0.28 | 0.48 | 0.09 | 0.59 | 0.24 | 0.28 |
R22 | 0.36 | 0.31 | 0.40 | 0.43 | 0.55 | 0.28 | 0.17 |
R410A | 0.45 | 0.31 | 0.37 | 0.32 | 0.60 | 0.25 | 0.25 |
R744 | 0.39 | 0.62 | 0.56 | 0.22 | 0.58 | 0.99 | 0.62 |
Average | 0.38 | 0.41 | 0.46 | 0.29 | 0.58 | 0.54 | 0.37 |
Coolant . | Imura . | Kutateladze . | Labuntsov . | Rohsenow . | Kruzhilin . | Hasna . | Cooper . |
---|---|---|---|---|---|---|---|
R134a | 0.34 | 0.28 | 0.48 | 0.09 | 0.59 | 0.24 | 0.28 |
R22 | 0.36 | 0.31 | 0.40 | 0.43 | 0.55 | 0.28 | 0.17 |
R410A | 0.45 | 0.31 | 0.37 | 0.32 | 0.60 | 0.25 | 0.25 |
R744 | 0.39 | 0.62 | 0.56 | 0.22 | 0.58 | 0.99 | 0.62 |
Average | 0.38 | 0.41 | 0.46 | 0.29 | 0.58 | 0.54 | 0.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.
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
APPENDIX
Nomenclature . | |
---|---|
Nu | Nusselt number |
Re | Renold number |
Pr | Prandtl number |
Density | |
The bubble’s diameter | |
Surface tension | |
Latent heat | |
h | Heat transfer coefficient |
cp | specific heat at constant pressure |
q,qh | Heat flux |
Viscosity | |
Heat conductivity | |
Working pressure/critical pressure | |
M | Relative molecular weight |
Pa | Atmospheric pressure |
l | Liquid |
v | Vapor |
Nomenclature . | |
---|---|
Nu | Nusselt number |
Re | Renold number |
Pr | Prandtl number |
Density | |
The bubble’s diameter | |
Surface tension | |
Latent heat | |
h | Heat transfer coefficient |
cp | specific heat at constant pressure |
q,qh | Heat flux |
Viscosity | |
Heat conductivity | |
Working pressure/critical pressure | |
M | Relative molecular weight |
Pa | Atmospheric pressure |
l | Liquid |
v | Vapor |
Nomenclature . | |
---|---|
Nu | Nusselt number |
Re | Renold number |
Pr | Prandtl number |
Density | |
The bubble’s diameter | |
Surface tension | |
Latent heat | |
h | Heat transfer coefficient |
cp | specific heat at constant pressure |
q,qh | Heat flux |
Viscosity | |
Heat conductivity | |
Working pressure/critical pressure | |
M | Relative molecular weight |
Pa | Atmospheric pressure |
l | Liquid |
v | Vapor |
Nomenclature . | |
---|---|
Nu | Nusselt number |
Re | Renold number |
Pr | Prandtl number |
Density | |
The bubble’s diameter | |
Surface tension | |
Latent heat | |
h | Heat transfer coefficient |
cp | specific heat at constant pressure |
q,qh | Heat flux |
Viscosity | |
Heat conductivity | |
Working pressure/critical pressure | |
M | Relative molecular weight |
Pa | Atmospheric pressure |
l | Liquid |
v | Vapor |