Zhu 1
Zhu 1
Zhu 1
Abstract
This paper presents a two-dimensional analytical model for low-temperature cylindrical heat pipes. A closed-form
solution which incorporates liquidvapor interfacial hydrodynamic coupling and non-Darcian transport through the
porous wick for the rst time, is obtained for predicting the vapor and liquid velocity and pressure distributions. In
addition, the steady-state vapor and wall temperatures for a given input heat load in the evaporator region and a
convective boundary condition in the condenser region, are obtained. The eects of liquidvapor interfacial
hydrodynamic coupling and non-Darcian transport through the porous wick on the vapor and liquid velocity and
pressure distributions as well as the heat pipe capillary limit are discussed and assessed. The analytical solutions of
the axial vapor and wall temperature distributions, the vapor and liquid pressure distributions, and the centerline
vapor velocities compare very well with both experimental and numerical results. This work constitutes for the rst
time a comprehensive analytical solution which provides closed form solutions for the vapor and liquid ow as well
as the operating temperature and the maximum heat removal capability of the heat pipe. # 1999 Elsevier Science
Ltd. All rights reserved.
1. Introduction
Heat pipes are currently used in a wide variety of
heat transfer related applications. Analyses of heat
pipe operations, both analytical and numerical, have
been performed extensively by many investigators.
Almost all of the analytical studies have been concentrated on the dynamics of vapor ow. Liquid ow and
the liquidvapor coupling were mostly neglected in
analytical studies because of their complexity. A comprehensive analytical model for the overall simulation
of steady-state heat pipe operation is not available in
open literature. Cao and Faghri [1] investigated the
eects of heat pipe wall and the porous wick on the
heat pipe operation. They concluded that it is important to include the porous wick and the wall in heat
pipe analysis and to treat the entire heat pipe as a
single system rather than to analyze the vapor ow
alone. Rosenfeld [2] also reported the importance of
heat transfer within the wall and the porous wick in
the case of an asymmetric heat input.
Due to the diculty of obtaining an analytical solution for overall heat pipe operation, more and more
0017-9310/99/$ - see front matter # 1999 Elsevier Science Ltd. All rights reserved.
PII: S 0 0 1 7 - 9 3 1 0 ( 9 9 ) 0 0 0 1 7 - 4
3406
Nomenclature
A1
A2
B
C
D
G1
G2
h
hfg
K
ke
kwall
L
La
Lc
Le
M1
M2
p
pc
Q
r
rc
Ro
Rv
Rw
Re
T
Tb
u
ui
Uv
Ul
v
v1
v2
x
Greek symbols
E
porosity of the wick
p
g
porous wick shape parameter ( E=K)
m
dynamic viscosity [N s/m2]
r
density [kg/m3]
sl
surface tension of the working liquid [N/m].
Subscripts
i
liquidvapor interface
l
liquid phase
v
vapor phase.
Superscript
+
dimensionless quantity.
numerical models have been developed. Some comprehensive numerical models [35] cover both the vapor
ow and the liquid ow. The boundary and inertial
eects were included in these numerical models by
applying the generalized momentum equation in porous medium to describe the liquid ow in heat pipes.
The coupling of the liquid and vapor momentum
equations was also incorporated in these models by
applying either the LaplaceYoung equation [3,4] or
the momentum jump condition [5] at the liquidvapor
interface. In these models, the matching conditions of
velocity and shear stress at the liquidvapor interface
were neglected by assuming a non-slip condition and
neglecting the interfacial drag. The eects of these
matching conditions as well as the boundary and inertial eects on heat pipe operation have not been investigated in any of the previous studies.
Vafai et al. [612] have developed comprehensive
pseudo-three-dimensional analytical models for asymmetrical at-shaped, including both disk-shaped and
at-plate, heat pipes. They incorporated liquid ow,
secondary vapor ow and the eects of liquidvapor
hydrodynamic coupling and non-Darcian transport in
their models [811]. Their results show that, for the
at-shaped heat pipes, while the eects of liquidvapor
interfacial hydrodynamic coupling are negligible,
neglecting the boundary and inertial eects can lead to
signicant error in predicting liquid ow and the maximum heat transfer capability of the heat pipe.
3407
Fig. 1. Schematic of the heat pipe and the coordinate system used in the analysis.
3408
0
@x
@r
r
@ uv
@ uv
@ pv
@ 2 uv 1 @ uv
rv uv
vv
mv
r @r
@x
@r
@x
@ r2
@ pv
0
@r
rv vv x, r Rv rl vl x, r Rv
8
0RxRLe
< rv v1 ,
0,
Le RxRLe La
:
rv v2 ,
Le La RxRL
ml @ 2 ul 1 @ ul
m
r FE
@ pl
l ul l 1=2 j ul j ul
r @r
K
E @ r2
K
@x
and
mv
uv x, r Rv ul x, r Rv ui x
@ ul 1 @
rul 0
r @r
@x
@ uv
@ ul
m
l
@ r rRv
@ r rRv
Q
2rv pRv Le hfg
uv vv ul vl 0
x L:
uv vv ul vl 0,
r 0:
vv 0,
@ uv
@ pv
0
@r
@r
10
pv pl
11
12
In the above governing equations, rv is the vapor density and rl the liquid density, mv the vapor dynamic
viscosity, ml the liquid dynamic viscosity, K and E are
the permeability and porosity of the wick and F is a
geometric function based on the porous wick structure
and is calculated using the expression outlined in Vafai
[16], i.e., F=1.75/Z150E 3/2.
r Rw :
ul vl 0
13
3. Analytical solution
An in-depth integral analysis along with the method
of matched asymptotic expansions is employed to
obtain the closed-form analytical solution for the
vapor and liquid velocity and pressure distributions.
The temperature values are then obtained through the
use of a simple conduction model.
3.1. Vapor velocity prole
The following velocity prole is utilized for the
vapor ow within the heat pipe:
uv x, r Uv xa0 a1 r a2 r2
14
where
1
Uv x
pR2v
Rv
0
uv x, r2pr dr
15
8
r Rv
>
>
>
Ul x ui x Ul x exp p
>
>
K=E
>
<
ul x,r
Ul x
>
>
>
>
r Rw
>
>
Ul x exp p
:
K=E
3409
r Rv
0R p R1
K=E
outer solution
r Rw
1R p R0
K=E
17
r Rv
0R p R1
K=E
outer solution
18
r Rw
1R p R0
K=E
Eq.
(18) can be written in a more compact form, by utilizing the fact that the thickness of the interface regions,
p
K=E, is much smaller than (RwRv)/2. This results in
8
>
r Rv
>
>
p
u
x
U
x
exp
U
>
l
i
l
<
K=E
ul x,r
>
r
R
w
>
>
Ul x exp p
>
:
K=E
Rv RrR
uv x, r ui x 2Uv x ui x 1
r
Rv
2 #
19
Rw Rv
RrRRw
2
Rw Rv
2
16
Uv x
8
>
>
>
>
>
>
>
>
<
2v1
x,
Rv
2v1
Le ,
Rv
>
>
>
>
>
>
2v2
>
>
L x,
:
Rv
0RxRLe
Le RxRLe La
Le La RxRL
20
3410
Bv1 x,
Ul x
Bv1 Le ,
:
Bv2 L x,
0RxRLe
Le RxRLe La
Le La RxRL
21
where
4
Uv x
24
m gRv
p
The fact that 1/g= K=EW(RwRv)/2 is used in the
derivation of Eq. (24). Substituting Eqs. (20) and (21)
into Eq. (24) results in the following expression for
ui(x):
ui x Ul x
ui x
Cv1 x,
Cv1 Le ,
:
Cv2 L x,
0RxRLe
Le RxRLe La
Le La RxRL
25
where
!
2Rv
8r
22
1 2 2
2
Rv
m g Rv
p
where g= E=K is the shape parameter of the porous
wick, and r +=rl/rv, m +=ml/mv. If the hydrodynamic
coupling is neglected at the liquidvapor interface, the
following expression is obtained for B based on a noslip condition at the liquidvapor interface:
B
8
<
r R2w
2Rv
B 2
r Rw R2v
23
G1 v21 M1 v1 x 2 ,
pv x pv 0
M1 v1 L2e 2M1 v1 Le x,
:
G1 v22 M1 v2 x L2 M1 v2 L La Lc ,
G1 v21
8
B
m gR2v
26
27
where
G1
rv
D2 7D 16
3R2v
28
4mv
D 2
R3v
29
M1
D CRv
30
If the hydrodynamic coupling at the liquidvapor interface is neglected, Eq. (27) reduces to
8
>
>
>
>
>
>
>
>
>
<
16rv 2 8mv
2
v
v
1 x ,
3R2v 1
R3v
16rv 2 8mv
16m
v
L2e 3 v v1 Le x,
pv x pv 0
1
2 1
3
3R
R
Rv
>
v
v
>
>
>
>
>
16rv 2 8mv
8mv
>
2
>
>
: 3R2 v2 R3 v2 x L R3 v2 L La Lc ,
v
v
v
0RxRLe
Le RxRLe La
Le La RxRL
31
2La L2c
33
and
ml B
1
2
1
2K
gRv 1 Rw =Rv
"
M2
rl FE
1
CC 2B
B2
3K 1=2
gRv 1 Rw =Rv 2
ml
1
2C
B
gRv 1 Rw =Rv 2
2K
M2 0
0RxRLe
Le RxRLe La
Le La RxRL
32
34
Q 2pRo Lc hTwall,c Tb
35
G2 v1 x 2 M2 v21 x 3 ,
G2 v1 Le 2x Le M2 v21 L2e 3x 2Le ,
:
G2 v2 x L2 L La Lc M2 v22 x L3 L 2La L2c ,
pl 0 pv L G2 v2 L La Lc M2 v22 L
G2
8
<
where
3411
!
36
37
Q
Tv Tb
2pLc
1
hRo
38
lnRo =Rw
lnRw =Rv
kwall
keff
39
and
8
"
#
>
Q
lnRo =Rw
lnRw =Rv
Lc
1
>
>
>
1
0RxRLe
Tb
>
>
2pLc
kwall
keff
hRo
Le
>
>
>
<
Q
lnRo =Rw
lnRw =Rv
1
Twall x
Le RxRLe La
Tb
>
>
2pLc
kwall
keff
hRo
>
>
>
>
>
Q
>
>
:
Tb
Le La RxRL
2phRo Lc
40
3412
2sl
rc
41
2sl
rc
42
Substituting Eqs. (9), (31) and (33) into Eq. (42) and
solving for the maximum heat transport capillary limit
yields
Qmax
q
1
A22 8A1 sl =rc A2
2A1
43
M2 L 2La
12p2 R2v r2v h2fg
44
A2
G2 2M1 L La
4pRv rv hfg
45
2sl
A 2 rc
46
3413
Fig. 2. Comparison of the calculated vapor and wall temperature distributions with the numerical results given by Tournier and ElGenk [5] and the experimental results given by Huang et al. [13].
3414
Fig. 3. Comparison of the calculated vapor and liquid pressure distributions with the numerical results given by Tournier and ElGenk [5].
Fig. 4. Variations of the mean vapor velocity and the centerline vapor velocity along the heat pipe.
Fig. 5. Variations of maximum liquid velocity and interfacial velocity along the heat pipe.
Fig. 6. The eects of liquidvapor interfacial hydrodynamic coupling on the maximum liquid velocities.
3415
3416
Fig. 7. The eects of liquidvapor interfacial hydrodynamic coupling on vapor and liquid pressures.
Fig. 8. Boundary and inertial eects on the vapor and liquid pressure distributions.
6 show that the interfacial eects lead to a larger maximum liquid velocity. The larger the injection Reynolds
number, the larger the dierence caused by the interfacial eects. Fig. 7 shows the eects of liquidvapor
interfacial hydrodynamic coupling on the vapor and
liquid pressures. The interfacial eect on the vapor
pressure is negligible. However, neglecting liquid
vapor interfacial hydrodynamic coupling will lead to a
smaller liquid pressure drop along the heat pipe.
The boundary and inertial eects on vapor and
liquid pressure proles are shown in Fig. 8. As indicated in Eq. (27), the pressure distribution in the vapor
phase is not aected by the boundary and inertial
eects. However, neglecting the boundary and inertial
eects leads to an underestimation of the liquid pressure. The larger the Reynolds number, the larger the
error involved in the calculation of liquid pressures by
using a Darcian model. This is in agreement with the
results by Vafai and Tien [15].
Fig. 9 shows the interfacial and the boundary and
inertial eects on the heat transport capillary limit of
the heat pipe. It can be seen in Fig. 9 that either
neglecting the interfacial eects or neglecting the
boundary and inertial eects can lead to an overestimation of the capillary limit. This is because a smaller
overall liquid pressure drop along the heat pipe is
3417
Fig. 9. The eects of liquidvapor interfacial hydrodynamic coupling and the boundary and inertial eects on the heat pipe capillary limit.
3418
vapor interfacial hydrodynamic coupling and nonDarcian transport through the porous wick were also
analyzed. The results show that the interfacial eects
are small and can be neglected. However, substantial
errors can occur when using Darcy's law in calculating
the liquid pressure distributions as well as the maximum heat removal capability of the heat pipe.
Acknowledgement
The
grant
(DE-F602-93ER61612)
by
the
Department of Energy is acknowledged and greatly appreciated.
References
[1] Y. Cao, A. Faghri, Transient two-dimensional compressible analysis for high temperature heat pipes with
pulsed heat input, Numerical Heat Transfer 114 (1990)
10281035.
[2] J.H. Rosenfeld, Modeling of heat transfer into a heat
pipe for a localized heat input zone, in: AIChE
Symposium Series, no. 257, vol. 83, 1987, pp. 7176.
[3] V.H. Ransome, H. Chow, ATHENA heat pipe transient
model, in: Trans. 4th Symposium on Space Nuclear
Power Systems, CONF-870102-Summs., Albuquerque,
NM, 1987, pp. 389392.
[4] J.M. Doster, M.L. Hall Numerical modeling of hightemperature liquid-metal heat pipes, 1989 Joint ASME/
A.I.Ch.E National Heat Transfer Conference,
Philadelphia, PA, vol. 89-HT-13, 1989, pp. 19.
[5] J.M. Tournier, M.S. El-Genk, A heat pipe transient
analysis model, Int. J. Heat Mass Transfer 37 (1994)
753762.
[6] K. Vafai, W. Wang, Analysis of ow and heat transfer
characteristics of an asymmetrical at plate heat pipe,
Int. J. Heat Mass Transfer 35 (1992) 20872099.
[7] K. Vafai, N. Zhu, W. Wang, Analysis of assymetrical
disk-shaped and at plate heat pipes, ASME J. Heat
Transfer 117 (1995) 209218.
[8] N. Zhu, K. Vafai, The eects of liquidvapor coupling
and non-Darcian transport on asymmetrical disk-shaped
heat pipes, Int. J. Heat Mass Transfer 39 (1996) 2095
2113.
[9] N. Zhu, K. Vafai, Optimization of asymmetrical diskshaped heat pipes, AIAA J. Thermophys. Heat Transfer
10 (1996) 179182.
[10] N. Zhu, K. Vafai, Numerical and analytical investigation of vapor ow in a disk-shaped heat pipe incorporating secondary ow, Int. J. Heat Mass Transfer 40
(1997) 28872900.
[11] N. Zhu, K. Vafai, Analytical modeling of the startup
characteristics of asymmetrical at-plate and diskshaped heat pipes, Int. J. Heat Mass Transfer 41 (1998)
26192637.
[12] N. Zhu, K. Vafai, Vapor and liquid ow in an asymmetrical at plate heat pipe: a three-dimensional analytical and numerical investigation, Int. J. Heat Mass
Transfer 41 (1998) 159174.
[13] L. Huang, M.S. El-Genk, J-M. Tournier, Transient performance of an inclined water heat pipe with a screen
wick, ASME National Heat Transfer Conference,
Atlanta GA, Heat Pipes and Capillary Pumped Loops,
HTD-vol. 236, 1993, pp. 8792.
[14] C.A. Busse, F.C. Prenger, Numerical analysis of the
vapor ow in cylindrical heat pipes, in: Proc. 5th Int.
Heat Pipe Conf., Part 1, JIATEC, Tokyo, Japan, 1984,
pp. 214219.
[15] K. Vafai, C.L. Tien, Boundary and inertia eects on
ow and heat transfer in porous media, Int. J. Heat
Mass Transfer 24 (1981) 195203.
[16] K. Vafai, Convective ow and heat transfer in variable
porosity media, J. Fluid Mech 147 (1984) 233259.
[17] K. Vafai, R. Thiyagaraja, Analysis of ow and heat
transfer at the interface region of a porous medium, Int.
J. Mass Transfer 30 (1987) 13911405.
[18] A. Faghri, Performance characteristics of a concentric
annular heat pipe: Part IIvapor ow analysis, ASME
J. Heat Transfer 111 (1989) 851857.
[19] C.A. Busse, Pressure drop in the vapor phase of long
heat pipes, in: Proc. Thermionic Conversion Specialist
Conference, Palo Alto, CA, 1967, pp. 391398.
[20] L. Huang, M.S. El-Genk, Experimental investigation of
transient operation of a water heat pipe, in: Proceedings
of 10th Symposium on Space Nuclear Power and
Propulsion, Albuquerque, NM, 1993, pp. 365374.
[21] S.W. Chi, Heat Pipe Theory and Practice, Hemisphere,
Washington, DC, 1976.
[22] J-M. Tournier, M.S. El-Genk, A.J. Juhasz, Heat-pipe
transient model for space applications, in: Proceedings
of 8th Symposium on Space Nuclear Power Systems,
Albuquerque, NM, 1991, pp. 857868.