Condensation Inside Horizontal Tubes: State of The Problem and Analysis of Research Results
Condensation Inside Horizontal Tubes: State of The Problem and Analysis of Research Results
net/publication/331431457
CITATIONS READS
3 2,417
4 authors, including:
Volodymyr Sereda
National University of Water and Environmental Engineering
19 PUBLICATIONS 134 CITATIONS
SEE PROFILE
All content following this page was uploaded by Volodymyr Sereda on 01 March 2019.
Key words: Heat transfer, Condensation, Smooth tube, Horizontal tube, Inside, Re-
view, Refrigerant
Abstract
This paper presents analysis of the study on in-tube condensation performed until 2010.
It is found that in all published studies has been measured only tubes with average perimeter
and heat transfer coefficients at certain length (minimum 0,5 m) during condensation inside
smooth tubes that did not provide sufficient data concerning the features of heat transfer dur-
ing condensation, namely impact of vapour quality, vapour velocity, heat flax, condensate
accumulation on heat transfer. The disadvantages of the existing methods of heat transfer re-
search, of experimental data presentation, and of design correlations are pointed out.
Introduction
155
15th Workshop on Transport Phenomena in Two-Phase flow, 2011
Investigation methods
In all works beginning from [9] up to now ([8] is the last one known to the authors) the
method of investigation of heat transfer under condensation inside horizontal tubes is the
same and insignificantly differs only by diameters and lengths of the experimental sections,
number of the thermocouples to measure wall temperature (in some works heat-transfer coef-
ficients are determined by Wilson’s method, i.e., without measurement of wall temperature),
and the methods of surface cooling. The minimal length l of tested section took place in [10]
and equaled to 0.3 m. In other works l = 0.5 to 1.5 m. Wall temperature of tested tubes was
measured in three or four points around tube perimeter and in several cross-sections along the
length. Thus, in effect in these works perimeter-average heat-transfer coefficient α ϕ was
measured in specified points along the tube.
Schematic drawing of the experimental apparatus (Fig. 1) included evaporator 1, vapor
compressor 2, preliminary section of condensation 3 to create a certain vapor content x at the inlet
to tested section 4, condenser 5 to condense vapor downstream tested section 4, inspection sec-
tions 6, pressure and temperature gauges at the inlet and outlet of all condensation sections, and
condensate flow meters 7. In some works tested section had several subsections (minimal length
of each was 0.3 m) cooled by means of a separate water circuit.
1 – Evaporator
(water-cooled);
2 – Compressor;
3 – Pre-condenser (water-
cooled);
4 – Test section
(water-cooled);
5 – Post-condenser (water-
cooled);
6 – Sight glass;
7 – Mass flow meter;
p – Pressure tap;
T – Thermocouple.
Fig. 1. Schematic drawing of the experimental apparatus.
By setting small (0.1 to 0.2) changes in vapor content in tested section the authors sup-
posed that they measure heat-transfer coefficients for local specified values of G and x. Such
approach to the measurements and assessment of locality of α has the following disadvan-
tages.
The first disadvantage. To provide properly substantiated small variations in ∆х and G in
the experiments it is necessary to maintain heat flux density in the tested section decreasing as G
is decreasing and length of the tested section is growing under fixed ∆х. It follows from trivial
Gdhlv
correlation qmax = Δxmax . For example, with reference to R410A (hlv = 1.6·105 kJ/kg) and d
4l
= 8 mm at l = 1 m and ∆хmax = 0.1 the following takes place:
at G = 100 kg/m2/s qmax 2.5·103 W/m2;
and at G = 200 kg/m2/s qmax 5·103 W/m2.
156
15th Workshop on Transport Phenomena in Two-Phase flow, 2011
Actual value of q in industrial condensers by several times exceeds 2500 W/m2. It means
that effect of q on α has not been virtually investigated.
The second disadvantage of the used method is connected with the first. In real heat ex-
changers with long tubes (4 to 9 m) at vapor inlet into the tube, where х > 0.7 to 0.5, heat-
transfer coefficients could be several times higher than at the end of the tube, where x → 0 .
Hence, heat flux q at the tube inlet is also significantly higher than at its end. As a result, be-
havior and rate of heat transfer either along the tube or depending upon variation in х in real
heat exchangers significantly differ from those obtained under condensation, when low values
of q are fixed.
The third uncertainty of existing methods of investigation lies in the following. Due to
impact of gravity on flow of condensate film virtually at х < 1 under all values of G conden-
sate drains along tube perimeter. When tested sections are longer than 0.3 m, ϕ-average heat-
transfer coefficients α ϕ vary even under negligible variations in G and x along tube length.
Therefore, depending upon the values of G, q, x, l, and physical properties of vapor the
values of heat-transfer coefficients estimated by the data based on the applied methods of
measurement can be substantially different.
157
15th Workshop on Transport Phenomena in Two-Phase flow, 2011
much higher than in the lower one. It is reflected in works [12] and [19], where just local α ϕ
and qϕ were measured.
Under low G and х, neither impact of х nor unambiguous impact of G on α were ob-
served that is seen from the analysis of work [7] and demonstrated in Fig. 2. It is also clear, if to
keep in mind that in the region of condensation by Nusselt the heat transfer depends only upon
q (∆T) similarly to condensation in outer surface of horizontal tube.
1 – G=300 kg/(m2·s);
2 – G=250 kg/(m2·s);
4 – G=175 kg/(m2·s);
5 – G=150 kg/(m2·s);
6 – G=125 kg/(m2·s);
7 – G=100 kg/(m2·s);
8 – G = 75 kg/(m2·s).
Fig. 2. α plotted against
х for condensation of R134a at
different G.
Theoretical solution of the problem of condensation inside horizontal tubes was carried
out for two modes of two-phase flow: for the mode with prevailing impact of gravity force (in
the works this mode is called stratified) and for annular mode of phase flow, when gravity is
neglected.
In all works, for the stratified mode Nusselt’s correlation [15] derived for condensation
on outer surface of horizontal tubу is used:
⎛ λ 3ρ ( ρ − ρv ) ghlv ⎞
0.25
α st = 0.725 ⎜ l l l ⎟ , (1)
⎝ μl d ΔT ⎠
with various empirical corrections taking into account availability of stream in the lower part
of tube, where heat transfer is determined by formula for convective heat transfer, for exam-
ple, in [20].
Theoretical solution for annular mode of phase flow is analogous to those carried out in
different works for condensation inside vertical tube, e.g., in [16]. In this case, in the majority
of works beginning with [21 to 23] turbulent flow of condensate film and vapor is considered
and the structure of defining parameters, as well as degree of their effect on heat transfer in all
correlations are very similar (see Table 1).
158
15th Workshop on Transport Phenomena in Two-Phase flow, 2011
Nu = 0.023Rel0.8 Prl0.4 ×
5. Tang et al. ⎡ ⎛ x ⎞ ⎤
0.836
Frso > 7 [25]
× ⎢1 + 4.863 ⎜ − ln( pred ) ⎟ ⎥
⎣⎢ ⎝ 1 − x ⎠ ⎦⎥
Cavallini & λ
8. α = l 0.05Prl1 3 Ree0.8q [28]
Zecchin d
⎛J ⎞
α an − st = ( α an , JG = 2.5 − α st ) ⎜ G ⎟ + α st
⎝ 2.5 ⎠
−1
⎡ ⎛1− x ⎞ ⎤
0.268
J G < 2,5
9. Cavallini α st = 0.725 ⎢1 + 0.82 ⎜ ⎟ ⎥ × [20]
⎢⎣ ⎝ x ⎠ ⎥⎦ X tt < 1.6
⎡ λ 3ρ ( ρ − ρv ) ghlv ⎤
0.25
⎛1− θ ⎞
×⎢ l l l ⎥ + αl ⎜ ⎟
⎣ μ l d ΔT ⎦ ⎝ π ⎠
α f θ + ( 2π − θ ) α с
α lp = ,
2π
⎡ ρ ( ρ − ρv ) ghlv λ 3l ⎤
10. Thome α f = 0.728 ⎢ l l ⎥ [5]
⎣ μv dq ⎦
λ
αc = 0.003Rel0.74 Prl0.5 l fi ;
δ
159
15th Workshop on Transport Phenomena in Two-Phase flow, 2011
In [31], five correlations from Table 1 (3, 6, and 8 to 10) for condensation of refriger-
ant R290 at р = 2 MPa and din = 11.9 mm, three values of mass velocity (50, 125, and 200
kg/m2/s) and three values of q (2000, 8000, and 19000 W/m2) are compared. Even in such
example, limited by physical and technical parameters, variety in the values of α calculated by
different correlations reaches 60%, as seen in Fig. 3.
3 – Dobson-Chato [26];
4 – Cavallini and
Zecchin [28];
5 – Shah [23].
Fig. 3. α plotted against х
for condensing R290
at q = 2000 W/m2 and G =
100 kg/(m2·s).
Poor compliance of various correlations from Table 1 with the experimental results
obtained by different authors should be noted. For example, Kaushik and Azer [27] propose
correlation (7) that generalizes their experimental data on R113 to within ±30%. While in the
work by Fujii [24] the data on condensation of R22, R134a, and R123 are compared with de-
sign correlations (3), (7), and (8) and, as is shown, these correlations give heat-transfer coeffi-
cients lower by 100% and more and only own correlation (4) generalizes the experimental
data to within below 20%.
In [3], the experimental data for various refrigerants are compared as a function α = f ( x ) at
the values of G from 74.6 to 644 kg/(m2·s). It is found that both good compliance of all design
correlations with the experimental data and discrepancy of some results up to 100%, despite
each of the authors of the design correlations demonstrates convergence with own experimen-
tal data within 30%. Fig. 4 from [3], where for the same G (600 kg/(m2·s)) the design and ex-
perimental data for refrigerants R32 and R410A are compared, proves that the proposed de-
sign correlations are not of generalizing type. It is seen that the same correlations satisfacto-
rily comply both with experimental data and with each other for R410A and substantially dis-
agree (up to 100%) for R32.
In [4] and [5], a “new” model for calculation of heat transfer is proposed. It is based
on a simplified flow pattern map with two selected zones: the zone, where heat transfer de-
pends on ∆Т (this zone is closer to the stratified mode of phase flow, i.e., J G ≤ J GT ) and the
zone, where heat transfer does not depend on ∆Т ( J G > J GT ). Such approach is not correct at
least because in an annular mode of phase flow, particularly under laminar film flow of con-
densate, heat transfer decreases with growing ∆Т.
160
15th Workshop on Transport Phenomena in Two-Phase flow, 2011
α,
1 – experimental
data;
2 – Dobson-Chato
[26];
3 – Shah [23];
4 – Fujii [24];
а)
α, 5 – Tang et al.
[25];
6 – Сavallini et al.
[20];
7 – Bivens and
Yokoseki;
8 – Tandon et al.
b)
Fig. 4. α plotted against х for condensation of R410A (a) and R32 (b)
at G = 600 kg/(m2·s).
161
15th Workshop on Transport Phenomena in Two-Phase flow, 2011
The experimental data of many authors are compared in [4] with the proposed calcula-
tion model as the plots of Nucalc Nuexp = f ( x ) . Analysis of them shows that deviation of the
experimental data from those designed both for the refrigerants and for water steam can exceed
60%.
3 – Shah [23];
4 – Dobson-Chato [26];
In works [7] and [8], the attempts to improve convergence of the experimental data with
calculation models [4] and [5] due to refinement of the transition boundaries of the zones with
heat transfer dependent and independent on ∆Т. Partially improved convergence of the calcula-
tions with the experiments on R134a within the range of G from 75 to 300 kg/m2/s and very low
q (1200 to 3600 W/m2) was obtained.
In Fig. 5, the experimental data on water steam obtained in [11] are compared with the
design correlations (3, 6, 9, and 10) from Table 1. It is distinguishing feature of these data that
in [11] all of the data on x, q, ∆T, G, and α are given in tabular form. Since length of the tubes
l = 2.5 to 12 m, then we calculated the data both at х = 0 to 1.0 and at ∆х not exceeding 0.2 to
0.3. As seen from Fig. 5, there is divergence of all correlations except for [23] that is close to
formula by Boyko.
Conclusion
162
15th Workshop on Transport Phenomena in Two-Phase flow, 2011
from the point of view of heat transfer under flow of condensate film on the greater part of the
surface, because such analysis does not explain impact of either heat flux or temperature drop
∆Т.
It is shown that the existing correlations for calculation of local and mean coefficients
of heat transfer inside smooth tubes do not correspond with the data obtained in the experi-
ments on different fluids. Discrepancy of calculated and experimental data reaches 30 to 70%
and more.
Nomenclature
163
15th Workshop on Transport Phenomena in Two-Phase flow, 2011
Subscripts:
eq equivalent
in – tube inlet;
out – tube outlet;
l – liquid phase;
lp local perimeter
v – vapor phase;
f - film;
φ – tube perimeter;
st - stratified;
an - annular;
tr - transient;
с – convective;
m – mean;
calc – calculated;
exp – experimental
Superscripts
− indicates average
References
1. Webb R.L. Principles of Enhanced Heat Transfer. – N.Y.: John Wiley & Sons,
1994. – 556 p.
2. Rifert V.G., Smirnov H.G. Condensation Heat Transfer Enhancement. – South-
hampton, Boston.: WIT-Press, 2004. – 372 p.
3. Garcia-Valladares O. Review of In-Tube Condensation Heat Transfer Correla-
tions for Smooth and Microfin Tubes // Heat Transfer Engineering. – 2003. – vol.
24(4). – P. 6-24.
4. Cavallini A., Censi G, Del Col D., et al. Condensation of refrigerants in smooth
tubes: a new heat transfer model for heat exchanger design // 3rd International
Conference on Heat Transfer, Fluid Mechanics and Thermodynamics. – 21-24 of
June 2004, Cape Town, South Africa.
5. Ei Hajal, Thome J.R., Cavallini A. Condensation in horizontal tubes. Part 1: Two-
phase flow pattern map // International Journal of Heat and Mass Transfer. –
2003. – vol. 46. – P. 3349-3363.
6. Olivier J.A., Liebenberg L., Thome J.R., Meyer J.P. Heat transfer, pressure drop,
and flow pattern recognition during condensation inside smooth, helical micro-fin,
and herringbone tubes // International Journal of Heat and Mass Transfer. – 2007.
– vol. 30. – P. 609-623.
7. Suliman R., Liebenberg L., Meyer J.P. Improved flow pattern map for accurate
prediction of the heat transfer coefficients during condensation of R-134a in
smooth horizontal tubes and within the low-mass flux range // International Jour-
nal of Heat and Mass Transfer. – 2009. – vol. 52. – P. 5701-5711.
8. van Rooyen E., Christians M., Liebenberg L., Meyer J.P. Probabilistic flow pat-
tern-based heat transfer correlation for condensing intermittent flow of refriger-
ants in smooth horizontal tubes // International Journal of Heat and Mass Transfer.
– 2010. – vol. 53. – P. 1446-1460.
164
15th Workshop on Transport Phenomena in Two-Phase flow, 2011
9. Azer N.Z., Ahis L.V., Soliman H.M. Local Heat Transfer Coefficients during An-
nular-Flow Condensation // ASHRAE Transaction. – 1979. – vol. 78, Part 2. – P.
1446-1460.
10. Islam M.A., Miyara A. Liquid film and droplet flow behavior and heat transfer
characteristics of herringbone microfin tubes // International Journal of Refrigera-
tion. – 2007. – vol. 30. – P. 5701-5711.
11. Boyko L.D. Heat Transfer in Condensing Vapor inside Tubes (in Russian)// Heat
Transfer in the Elements of Power Plants (Tyeploobmyen v elementakh ener-
geticheskikh ustanovok). – 1966. – P. 197-212.
12. Volkov D.I. Heat Transfer in Condensing Vapor inside Horizontal Tubes (in Rus-
sian) //TsKTI Transactions (Trudy TSKTI). – 1965. – Issue 57.
13. Konsetov V.V. Heat Transfer in Condensing Vapor inside Horizontal Tubes (in
Russian) // Journal of Engineering Physics (Inzhenyerno-fizicheskiy zhurnal). –
1960. – #6.
14. Gorodinskaya S.L. Heat Transfer in Condensing Vapor inside Horizontal Tubes
(in Russian) // Proceedings of KPI (Izvyestiya KPI). – 1954. – vol. 13. – P. 158-
170.
15. Nusselt W. Die Oberflachen Kondensation des Wasserdampfs // Zs. des Vereines
deutscher Ingenieure. – 1946. – Bd. 60. – S. 568-575.
16. Duckler A.E. Fluid mechanics and heat transfer in vertical falling film systems //
Chem. Engng. Prog. Syrup. Ser. – 1960. – vol. 56. – P. 1-10.
17. Sparrow К. Turbulent Flow, Heat, and Mass Exchange in Tube with Surface Ex-
haust // Heat Transfer. – 1970. – vol. 92. – P.121-131.
18. Rifert V.G., Tobilyevich A.N., Sardak A.I. Flow and Heat Transfer Modes in Con-
densing inside Horizontal Tubes (in Russian)// Proceedings of the Academy of
Sciences of the USSR. Power Engineering and Transport (Energetika i transport
). – 1985. – vol. 4. – P. 101-109.
19. Rifert V.G. Heat transfer and flow model of phases in laminar vapor condensation
inside a horizontal tube // International Journal of Heat and Mass Transfer. –
1988. – vol. 31. – P. 547-558.
20. Cavallini A., Censi G., Del Col D., et al. Condensation inside and outside smooth
and enhanced tubes – a review of recent research // International Journal of Re-
frigeration. – 2003. – vol. 26. – P. 373-392.
21. Kutateladze S.S. Heat Transfer in Film Condensing of Vapor inside Horizontal
Tube (in Russian)// In articles: Problems of Heat Transfer and Hydraulics of Two-
Phase Media (Voprosy tyeplootdachi i gidravliki dvukhfaznykh sred). Moscow:
Gosenergoizdat. – 1961.
22. Anan’yev E.P., Boyko L.D., Kruzhilin G.H. Heat transfer in the presence of steam
condensation in a horizontal tube // International Heat transfer conference. – Part
II. – 1961.
23. Shah M.M. A general correlation for heat transfer during film condensation inside
pipes // Int. J. Heat Mass Transfer. – 1979. – vol. 22. – P. 547-556.
24. Fujii T. Enhancement to Condensing Heat Transfer – New Developments // Jour-
nal of Enhanced Heat Transfer. – 1995. – vol. 2. – P. 127-137.
25. Tang L., Ohadi M. M., Johnson A. T. Flow Condensation in Smooth and Microfin
Tubes with HCFC-22, HFC-134a, and HFC-410 Refrigerants. Part II: Design
Equations // Journal of Enhanced Heat Transfer. – 2000. - vol. 7. – P. 311-325.
165
15th Workshop on Transport Phenomena in Two-Phase flow, 2011
166