J. R. Simonson (Auth.) - Engineering Heat Transfer-Palgrave Macmillan UK (1988)
J. R. Simonson (Auth.) - Engineering Heat Transfer-Palgrave Macmillan UK (1988)
J. R. Simonson (Auth.) - Engineering Heat Transfer-Palgrave Macmillan UK (1988)
J. R. Simonson
Senior Lecturer
The City University, London
SECOND EDITION
M
MACMILLAN
EDUCATION
© J. R. Simonson 1975,1988
Published by
MACMILLAN EDUCATION LTD
Houndmills, Basingstoke, Hampshire RG21 2XS
and London
Companies and representatives
throughout the world
Nomenclature xi
1 Introduction 1
5 Transient Conduction 58
5.1 The Uniform Temperature, or Lumped Capacity, System 58
5.2 The Solution of Transient Conduction Problems in
One Dimension 61
5.3 Two-dimensional Transient Conduction 67
5.4 Periodic Temperature Changes at a Surface 68
vi CONTENTS
6 Forced Convection: Boundary Layer Principles 78
6.1 Introduction 78
6.2 Equations of the Laminar Boundary Layer on a
Flat Plate 80
6.3 Laminar Forced Convection on a Flat Plate 87
6.4 Laminar Forced Convection in a Tube 92
Index 261
Preface to First Edition
a distance increment
A area
b, 1, t, w linear dimension
c capacity ratio of heat exchanger
C,K constants of integration
Cd average friction factor
Cf skin friction coefficient
Cp specific heat at constant pressure
Cp volumetric specific heat at constant pressure
d diameter
E effectiveness of heat exchanger
f friction factor
F geometric configuration factor
~ geometric emissivity factor
/o drag factor
g gravitational acceleration
G irradiation, mass velocity
Gz Graetz number, Re Pr(djx)
hR convection coefficient
H product hA
hrg latent enthalpy of evaporation
hr radiation coefficient
i current density
I current
I intensity of radiation
J radiosity
k thermal conductivity
L,D,T, W linear dimension
L,M,T,O dimensions of length, mass, time, temperature
m mass flow, or mass in transient conduction
m mass flow, where a non-flow m also occurs
n coordinate direction
n frequency of temperature variation
NTU number of transfer units
p,P,flp pressure, difference of pressure
p perimeter
PN plate number
ql heat transfer per unit area and time
q heat generation per unit volume and time
Q heat transfer per unit time, or a physical variable in dimension-
less analysis
xi
XII NOMENCLATURE
Dimensionless groups
F Fourier number, Ata.ja 2
Gr Grashofnumber,~g8p2 z3tp.2
J Colburn /-factor, St. Pr 2 13
Nu 1 Nusselt number, hl/k
NOMENCLATURE Xlll
Suffices
a at axis of tube
b black body
b limit of laminar sub-boundary layer
c cold fluid
c convection
d, /,X length terms used in dimensionless groups
e equivalent
f fluid
h hot fluid, heated length
i, 0 inlet, outlet (in heat exchangers)
l liquid
m mean value
M metal, in heat-exchanger wall
n direction of component
0 datum length
p constant pressure
r radial direction, or radial position
R radiation
s surroundings, of free stream
sat. saturated temperature
t temperature, turbulent
v constant volume, vapour
w wall
x,y,z direction of component
(J angular component
A monochromatic
Superscript
average value
1
Introduction
The rate of heat flow (in only the x-direction, see Fig. 1.1) is propor-
tional to the product of the area of flow and the temperature gradient,
the constant of proportionality being the thermal conductivity k
which is a property of the material. The negative sign results from
the convention of defining a positive heat flow in the direction of a
negative temperature gradient. The property k may be a function of
temperature and direction of heat flow. Materials with directional
dependence of thermal conductivity are said to be anisotropic.
The units involved depend on the system chosen. In the SI system,
the unit of heat or internal energy is the joule, hence rate of heat
transfer is measured in Jjs or W. However, the kilojoule, (kJ), and
kilowatt, (kW~ are accepted multiples of the SI unit, and to be consis-
tent with general usage in thermodynamics, the kJ and kW are the
preferred units in this book. With the area in m 2 and the temperature
gradient in K/m, the units of k are kW /(m K). This follows the British
Standards recommendation 4 for the presentation of complex units.
4 ENGINEERING HEAT TRANSFER
uniform temperatures over each face
I 41
Q = -kA dt
" dx
Fig. 1.1. Follrier's hlw for ou-tlinuuio•fll stellll:y st11te co1Uiuctio11 ;, 11
pl~~~~e sld of 11111terill/.
fluid, t 2
Q = hA(t 1 - t 2)
Fig. 1.2. Dillgrtllllllllltic representation of convection from a flat plate,
i/butrating the use of Newton's equation.
Table 1.1
Convection system Range of h,* kW/(m 2 K)
REFERENCES
1. Fourier, J. B. Theorie analytique de Ia chaleur, Paris, 1822. Translated by
A. Freeman, Dover Publications, New York, 1955.
2. Biot, J. B. Bibliotheque Britannique, Vol. 27, 310 (1804).
3. Biot, J. B. Traite de physique, Vol. 4, 669 (1816).
4. British Standards Institution, The Use of SJ Units, PD5686: 1972.
5. Newton, I. Phil. Trans., Roy. Soc., London, Vol. 22,824 (1701).
6. Boltzmann, L. Wiedemanns Annalen, Vol. 22, 291 (1884).
7. Stefan, J. Sitzungsber. Akad. Wiss. Wien. Math.-naturw. Kl., Vol. 79, 391
(1879).
2
The equations of heat conduction
q = Q, = _kat (2.2)
II A on
where q, is the heat ftux in heat units per unit time and per unit
area in the n-direction. This is a vector quantity since it has magni-
tude and direction. The greatest heat flux at an isothermal surface
will always occur along the normal to that surface.
Heat conduction within a solid may be visualized as a heat flux
which varies with direction and position throughout the material.
This follows from the fact that temperature within the solid is a
function of position coordinates of the system (e.g., x, y, z). In addi-
tion, temperature may be a function of time, (t), so in general
t = f(x, y, z, t).
The problem of determining the magnitude of heat conduction
8
THE EQUATIONS OF HEAT CONDUCTION 9
resolves itself to finding first the isotherms within the system and
the way in which their positions vary with time. In steady state
conduction the isotherms remain stationary with time, and one may
visualize a large number of isothermal surfaces throughout the
system, differing incrementally in temperature. The heat flux
normal to any one surface will vary with position depending on the
distance between surfaces. It is then necessary to sum the heat flow
through the boundary surfaces if internal heat sources are present, or,
if not, through any isothermal surface. In unsteady conduction the
problem is complicated by the fact that isothermal surfaces are no
longer fixed, and the rate at which heat is being stored must be taken
into account.
Before taking the first step, which is to develop the equation for
temperature as a function of position and time, it is opportune to
introduce some facts about different conducting materials.
Solid materials may be divided into two groups, metallic and
non-metallic, for which there is a marked contrast in the values of
conductivity. The Appendix lists properties for some of the more
useful materials. The high values of conductivity of metals are
attributable to the well ordered crystalline structure of the material.
The close arrangement of molecules permits a rapid transfer of
energy and, in addition, free electrons play a considerable part.
Metals such as copper which are good electrical conductors also
conduct heat well. There is also a marked similarity between
conduction heat transfer and the flow of electricity, and the electrical
analogy is often used in the solution of conduction problems.
In contrast, non-metals do not have a well ordered crystalline
structure and, in addition, are often porous in nature. Thus energy
transfer between molecules is seriously impeded, and the values of
conductivity are much lower. The small pores within the material,
being full of air, further restrict the flow of heat since gases are poor
conductors. This is because molecules of a gas are relatively widely
spaced and the transfer of energy depends on collisions between
these molecules.
The thermal conductivities of most substances V< cry with tempera-
ture, and for accuracy such variation should be allowed for in
conduction problems. However, this is a complication which may
be ignored in an introductory study of the subject because the
variation with temperature is not great. Over a reasonable tem-
perature range the relationship between conductivity k, and tern-
10 ENGINEERING HEAT TRANSFER
perature t may be assumed linear:
k, = k 0 (1 + ext) (2.3)
where k0 is the conductivity at temperature t 0 , and ex is a constant.
In most practical applications it is sufficient to assume a mean
uniform value for conductivity.
A complication more serious than temperature variation of con-
ductivity occurs in certain engineering materials, viz., that conduc-
tivity may vary with the direction of heat flow. This arises commonly
in laminated materials used in electrical engineering. Thus the
conductivity parallel to the laminates is different to the value
perpendicular to the laminates. Most types of wood also exhibit
this property, the conductivity parallel to the grain being different
to that across the grain. Conducting materials exhibiting this
property are said to be anisotropic. In the absence of this property
the material is said to be isotropic. For anisotropic materials the
analysis of conduction is more difficult and is not included in this
introductory text The basic ideas are given by Eckert and Drake, 1
and the general treatment may be found in the work of Carslaw and
Jaeger. 2
Differential equations of the temperature field will now be
developed in two coordinate systems, Cartesian and cylindrical.
ar
dQ,. = - kdydz ox
dQ(z+dz)
dQ,
ot
- kdydz ox- kdxdydz 0 x 2
o2 t
Therefore the net rate of heat flow into the element in the x-direction
is
(2.4)
In a similar manner, the net rates of heat flow into the element in the
y- and z-directions are given by
ot 2
dQY - dQ(y+dy) = k dX dy dz iJyl (2.5)
iJlt
dQZ - dQ(z+dz) = k dx dy dz OZ2 (2.6)
The total rate of heat flow into the element is the sum of the right-
hand sides of equations (2.4), (2.5), and (2.6), which is
o2 t o2 t o2
k dx dy dz (iJx2 + oy2 + oz2
t)
In addition to heat flowing into and out ofthe element, the possi-
bilities of heat being generated within the element (e.g., due to the
flow of electricity) and of heat being stored within the element (as
in the case of unsteady conduction) have to be considered.
12 ENGINEERING HEAT TRANSFER
If q' is the rate at which heat is generated per unit volume, the
rate of heat generation within the element is
q'dxdydz
The rate at which heat is being stored within the element is
governed by the rate of temperature change otjot. If pis the density
of the material and eP the specific heat at constant pressure, the
rate of heat storage will be
dx dy dz peP
at
at
This assumes the element may expand or contract freely at constant
pressure. The rate at which heat is being stored within the element
is equal to the sum of the net rate of heat flow into the element and
the rate of heat generation within it, hence:
ar (o2 t o2 t o2 t)
pep ot = k ox 2 + oy 2 + oz 2 + q'
(2.8)
and
(2.10)
and, consequently,
(2.11)
and
d2t
dx2 = 0 (2.12)
Hence
ot o2 t
dQ, - dQ 1,+dr) = k dz dr dO or + k dz rdO or 2 dr (2.13)
r dQ(r+dr)
In a similar manner, the net heat flowing into the element in the
circumferential direction is found to be
iJ2t
dQ 11 - dQ111 +dlll = k dr dz , 2062 rd(J (2.14)
Then an energy balance for the element leads to the general dif-
ferential equation for heat flow in three dimensions in a cylindrical
coordinate system. i.e.,
at (o 2t 1 at 1 o2t o2t)
Ot = IX or 2 + Orr + r2 iJ(J 2 + oz 2 + pep
q'
(2.16)
d 2 t +! dt = 0 (2.18)
dr 2 r dr
PROBLEM
Show that the general equation of heat conduction in spherical coordinates
is given by
REFERENCES
l. Eckert, E. R. G., and Drake, R. M. Introduction to the Transfer of Heat
and Mass, 2nd ed., McGraw-Hill, New York (1959).
2. Carslaw, H. S., ·and Jaeger, J.C. Conduction of Heat in Solids, Oxford
University Press (1947).
3
One-dimensional steady state
conduction
qx
• X
tl
R=
tl
-, dt
~ j_
- ~
t
tz
q,. dx
XI
X
Xz
16
ONE-DIMENSIONAL STEADY STATE CONDUCTION 17
and hence
- kA(t 2 - t 1)
Q =---- (3.1)
X (x2 - x.)
where (x 2 - x.) is the thickness of the slab and A is the area of the
slab. Using consistent units, the heat transfer calculated will be in
heat units per unit time.
The same result will be obtained if the appropriate differential
equation is integrated. Integration twice of equation (2.12),
gives
(3.2)
where C 1 and C 2 are constants of integration to be determined from
the boundary conditions, i.e., the temperatures at x 1 and x 2 •
Equation (3.2) indicates that the temperature variation through the
slab is linear. The temperature gradient from equation (3.2) used in
Fourier's law gives equation (3.1). Equation (3.1) may be re-written
as
Qx k(t1 - t2)
-A = qx = x - x ' (3.3)
2 1
\
ftu;~
--, t
tfl I
Ifluid 2
·~ q. I qc + q, = q"
l x~'>
' I
x2 \:
tf2
kl k2 X31
Fig. 3.2. A nudtipk pltute slllb separating two jblids, one-dimemioiUII stetuly
state colfliMctio11.
ONE-DIMENSIONAL STEADY STATE CONDUCTION 19
As a convenience, the radiation exchange between the wall and
fluid or some other surface beyond the fluid may be expressed by an
analogous equation
(3.5)
Since radiation exchanges are a function of the fourth power of the
absolute temperatures involved, the radiation coefficient ~ is
heavily temperature-dependent.
The total heat flow or conduction flux q" from the wall by convec-
tion and radiation is found by adding (3.4) and (3.5):
q" = q. + qr = - hc(tf - tw) - hR(tr - tw)
(3.6)
A multiple slab of two layers of conductivities k 1 and k2 which
separates two fluids f1 and f2 at temperatures tfl and tr 2 is now
considered. For exchange between fluid f 1 and wall surface at t 1 :
qx = -(he+ hR)fl(tl - tfl) (3.7)
For conduction through the two layers of material:
- kl(t2 - tt} - k2(t3- t2)
q" = x2 = (3.8)
- x1 x3 - x2
For exchange between the wall surface at t3 and the fluid f2 :
q" = - (h. + hR)r 2 (tr 2 - t 3) (3.9)
q" [
1
+ hR)fl
+
x2 - x1
+
x3 - x2
+ -c:----:---
1 J
(h. kl k2 (h 0 + hR)f2
= - (tf2- ted
and hence
(3.10)
where
1 1 x2 - x1 x3 - x2 1
- = + + + (3.11)
U (he + hR)fl k1 k2 (he + hR)C2
1/U is the overall thermal resistance per unit area between fluids
and U is the overall heat transfer coefficient. The resistances to heat
20 ENGINEERING HEAT TRANSFER
flow due to convection and radiation act in parallel and the
resistances due to the conducting layers act in series. The heat flow is
calculated from (3.10), once the overall coefficient U is found from
(3.11), and interface temperatures follow from (3.7), (3.8), and (3.9).
Sometimes in composite structures slabs of differing thermal
conductivity are present as shown in Fig. 3.3. This situation may be
treated one-dimensionally provided it is assumed that the y-z faces
ofthe intermediate slabs have uniform temperatures. The total resis-
tance may be deduced by adding the intermediate resistances in
parallel before adding the others in series.
-,...-
A2 k2
-Q
AI A3
kl
k3 ks
A4 k4
I I I I
XI x2 XJ x4
II 12 13 14
II R2 14
Rl
R4
-Q
Fig. 3.3. Olle-dimeiiSiolllll at~tuly attlte colflblctioll ;, tl a~riea-ptll'tllkl
ayat~m of altlba; tz tUUl t3 tu~ •lliform temp~rat•rea ;, the y--z pl1111e.
EXAMPLE 3.1
The reduction of heat loss from buildings is of very great practical
and economic importance. The Chartered Institution of Building
Services Engineers in the CIBS Guide Book A, give values of overall
heat transfer coefficient (known as U values and expressed as
W/(m 2 K) and not kW/(m 2 K)) for various types ofwal~ window, and
roof. Some typical values are given in Appendix 3. Use the U values
given below to calculate the heat transfer rate through a house
ONE-DIMENSIONAL STEADY STATE CONDUCTION 21
structure in cases (i) and (ii). In both cases the wall area is 110m2 ,
wood-frame window total area 14m2 , upstairs ceiling area 36m 2 ,
environment temperature difference 21 oc.
Case (i) 335 mm solid brick wall, U = 1·5 W/(m 2 K); pitched
roof with felt, foil-backed board ceiling, U = 1·5 W/(m 2 K);
single-glazed windows, U = 4·3 W/(m 2 K).
Case (ii) 335 mm solid wall plus 30 mm foam board lining,
k = 0·026 W /(m K); pitched roof as before plus 50 mm glass-fibre
insulation U = 0·5; double-glazed windows, U = 2·5 W/(m 2 K).
Solution. For parallel heat flow through walls, windows and roof, in
case (i)
Q = 21(110 X 1·5 + 36 X 1·5 + 14 X 4·3) = 5860 W
= 5·86 kW
In case (ii), the thermal resistance of the insulated wall is the original
resistance plus the insulation resistance which equals 1/1·5 +
0·03/0·026 = 1·82.
New U value = 1/1·82 = 0·55
Q = 21(110 X 0·55 + 36 X 0·5 + 14 X 2·5) = 2380 W
= 2·38 kW
A saving of 3·48 kW is achieved. Actual heating requirements will be
greater than the figures calculated on account of air changes, and
some losses through the ground floor.
EXAMPLE 3.2
3 157·32
7 107·92
11 99·56
15 103·96
20 114·85
25 126·87
heA - -- ----heR
hRA ___ _-
----hRR
-heR
-hRR
Troom
conv. conv.
rad. rad.
10 PRitU,"PLEASF. TYPE CONVECT! Oil COEFF!C!EIITo AIR SIDE [II W/112 K:"
20 IUPUT HCA
30 PR [lllo"PLEAH TYPE RAD!AT!OII COrFF !CI EUT, AIR SIDE Ill W/H2 K:"
40 Ill PUT liP A
50 rRif~T,"PLEASE TYPE CONVECTION COHF!C!ENTo ROOM SIDE [II W/H2 K:"
60 IUPIJT IICR
70 PRIIITo"PLEASE TYPE RADIAT!O~ COF.FFICIEIITo ~001-1 SIDE Ill W/"2 K :••
81) Ill PUT IIRR
90 PAftiT,"PLEASE TYPE WINDOW Fr.AIH AREA, [~ 112: ..
1r.0 triPUT WF.\
110 PRINT,"rLEASE TYPE ~[liD OW FRAIIF. TIIICY.NESS, [II H:u
120 !IIPUT IHT
130 rR!IIIo"PLEASE TYPE Wir~DOU GLASS AP.EA, IU 112: ..
140 [II PUT \IGA
150 PRIUT,"PLEA$[ TYPE WI~DOW GLASS TH!CKUESS, I~ MM:"
160 !UPHT W~T
165 UGT•WGT/10no.n
170 PRI~To"PLEASE TYP~ AIR GAP TH!CK"E5$o IH MM:"
180 IHPUT AGT
1r5 AGT•AGT/1000.0
19(\ TKF•0.166
1?5 TKA•O.O?~
200 TKG•0.76?
205 PRINTo"PLEASE TYPE ROOM TEMPERATURF IN OCG C:"
210 IUPUT TROQn
215 rni~T,"PL~A5E TYPE E~TER~AL AIR T~MPEQATURE IN DEG C:"
220 IUPUT TAIR
23(' r1=1.0/(IIH*(HCA+HRA))
240 P2=WFT/cwrA•Tr.r>
250 R3•1.0/(WFA•(HCR+HRR))
260 ~4=1.0/CWGA•CHCA+~RA))
270 R5•2.0•UGT/CWGA•TKGl+1.0/(T~A•WGA/AGT+(AGT-0.007l•350.0l
280 ~6•1.0/CWGA•(HCR+HRR))
290 RF•P1+R2+R3
300 r.W•R4+P5+R6
31C HOV•RF•RH/(RU+Rf)
320 QT•(TROOM-TAIR)/ROV
330 Of•(TROOM-TAIRl/RF
340 QW•(TROOM-TAIRl/RW
350 T1•TROO"-(R2+R3l•(TROOM-TAIRl/r.r
360 T2•TROOM-R3•CTROOM-TA!R)/~F
370 T3•TROO"-(R5oR6)•(TROOM-TA!P)/RW
380 T4•TROO"-R6•(TROOH-TAINI/R~
390 PRINT," RESULTS FOR AlP ~AP T~!CKNESS:"oAGT
400 PA!NTo"T,TAL HEAT THAIISFEP RATC 1:1 W:"oQT
410 rRUIT,''IIEAT TRI\tiSFER RAT£, FRAr"!F., IU ·.~:",QF
420 rAtfiT,''II~AT TR4fiSFER RATE, WifiDOW, I~ W:",OW
43~ rRJ"T,"OUT5tP( GLASS TEHPERATIJR~, I~ DEG C:",J3
440 rRIHT,"I:ISIDE GLASS TEMP(nATUR(, Ill DEG C:'',T4
450 PRINTo"QUTS!DE FRA~E TEMPERATURCo !U DEG C:"oT1
460 PRINT,"l~S!D[ fRAME TE~PENATURCo I~ PEG C:"oT2
470 PRI~T," DO YOU WISH TO PO A RfPEAT CALCULATIO~? 1-YESo2-NO:"
480 INPUT IX:
490 IFC!X~<2>GO TO 10
500 STOP
24 ENGINEERING HEAT TRANSFER
and
(3.12)
ONE-DIMENSIONAL STEADY STATE CONDUCTION 25
EXAMPLE 3.3
The heat flux through a plane slab 0·1 m thick is 146kW/m2 for
surface temperatures of 120° and 30°C. Find the value and sign of
IX in the thermal conductivity function given that k0 = 0·16 kW/(mK)
Here A is the surface area at the radius r, and obviously A will vary
between the inner and outer radii. It is therefore convenient to
consider a unit length of cylinder, when the radial heat transfer is
26 ENGINEERING HEAT TRANSFER
dt
Q, = - k(2nr) dr (3.13)
2nr is the area per unit length. Since the same quantity Q, is flowing
through a steadily increasing cylindrical area, the temperature
gradient must decrease with increasing radius. Integrating:
Q = - 2nk (t 2 - td (3.14)
r In r2/r I
By analogy with Ohm's law, the thermal resistance per unit length
of cylinder in this case is [In (r 2/r 1)]/2nk.
-Q,
• c=:J •
ft In r21 'I t2
R=--
2nk
d2t + ! dt = 0 ((2.18))
dr2 r dr
ONE-DIMENSIONAL STEADY STATE CONDUCTION 27
1 In r2!r1 In r3!r2 1
2nrlhcl 2nkl 2nk2 2n r3 hc 2
Fig. 3.6. Steady state radial conduction in concentric cylinders separating
two fluids.
28 ENGINEERING HEAT TRANSFER
ing two fluids f1 and f2 . It is assumed that heat transfer at the
surfaces is in each case predominantly by convection. Considering
unit length, at the inside surface:
Q, = - 2nr 1 hct(t 1 - tfl)
The same quantity is conducted through the two layers, hence
Q, = - 2nk (t2 - _!.!2 = - 2nk2 (t3 - t2)
1 In r /r In r 3 /r 2
2 1
or
Q, = - U(tr2 - tn) (3.17)
where
1 1 In r 2/r 1 In r 3 /r 2
U = 2nr 1 hct + 2nk 1 + 2nk + 2nr hc
2 3 2
U is the overall heat transfer coefficient per unit length between the
two fluids and 1/U is the thermal resistance of unit length and consists
of the sum of individual thermal resistances to conducted and
convected heat flow.
or
k = CritlCa1 rad"IUS
ro = h 0 0
(3.18)
This is the value of outer radius for which heat transfer through the
system will be a maximum. It follows that if for a given k and h the
external radius is less than k/h, then increasing the thickness of
insulation up to a radius of k/h will increase, rather than decrease
the heat loss from the system. The situation is likely to arise if k has
a relatively high value, and h a relatively low value.
(3.19)
and the overall heat transfer coefficient for a double spherical layer
separating two fluids f 1 and f2 is
1 1 r2 - r1 1 r3 - r2 1 )
1/( 2
4nr1hc1
+ 4nk1
-- 0
r1r2
+ 4n
-k-
2
0
r2r3
+ 4nr32h c2
EXAMPLE 3.4
25 0·125 0·607
40 0·140 0·620
55 0·155 0·626
70 0·170 0·628
90 0·190 0·625
200 0·300 0·577
10 R1=0.1
20 T1=200.0
30 TA=20.~
40 PRINT,"PLEASE TYPE INSULATIOH THICKUESSo MM:"
50 l:tPUT T•f
60 R2=TH/1~00.0+R1
70 PRI~To"PLEASE TYPE IHSULATIOH THERMAL CONDUCTIVITY, W/M DEG C:"
80 llfPU T TK
90 PRIHTo"PLEASE TYPE OUTER RADIUS COHVECTIOU COEFFICIENT, W/M2 DEG C:"
100 llfPUT H
11~ RCRIT•TK/H
120 RESIS=1.0/(6.2832•R2•Hl+ALOG<R2/ R1l/(6.2832•TK)
13U QfLOW•(T1-TAl/RESIS
1~0 PRII~T,''RESULTS FOR Atl OUTER RADIUS OF, H:'',R2
150 PRI~T,"CRITICAL RADIUS fOR THE COEFFICIEIHS USED:"oRCRIT
160 PRINT,~tiEAT FLO~, U/11 PIPE LEHGTII:",OFLOU
170 PRIUT,"OD YOU UISH TO DO A REPEAT QUH! 1-YES; 2-HO"
180 IriPUT IX::
190 lf(IX~<2l&O TO 40
200 STOP
(3.20)
h(t
a x1
-t)=-(-k(dt))
a dx x,
(3.21)
32 ENGINEERI NG HEAT TRANSFER
and
hb (t -t
x2 a
)=- k(dt)
dx ~
(3.22)
q = -
X
k(~)
dx X
(3.23)
r
Fig. 3.8. Radial conduction in a cylindrical layer with internal heat generation.
34 ENGINEERING HEAT TRANSFER
of rat which dt/dr = 0 gives the position of the maximum tempera-
ture, and this substituted in (3.24) gives the value of the maximum
temperature.
The following example illustrates the way in which ohmic heating
problems may be solved.
EXAMPLE 3.5
An internally cooled copper conductor of 4 em outer diameter and
1·5 em inner diameter carries a current density of 5000 amp/cm 2 •
The temperature of the inner surface is maintained at 70°C, and it
may be assumed that no heat transfer takes place through insulation
surrounding the copper. Determine the equation for temperature
distribution through the copper, hence find the maximum tempera-
ture of the copper, the radius at which it occurs, and the heat
transfer rate internally. Check that this is equal to the total energy
generation in the conductor. For copper, take k = 0·38 kW/(mK)
and the resistivity p = 2 x 10- 11 ohm metre.
Solution. If i is the current density,
q' = pi 2 = 2 X 10- 11 X (5000 X 104 ) 2 kWjm 3
= 5 x 104 kW/m 3
The boundary conditions are that at r = 0·75 em, t = 70°C and
that at r = 2 em, dt/dr = 0. This is because the heat transfer is zero
at r = 2 em. It follows that the maximum temperature also occurs
at r = 2 em. The constants of integration in equation (3.24) may
now be found.
dt q'r C7
-= --+-=0 at
dr 2k r
5 X 104 X 0·()2 + C7 = O
2 X 0·38 0·02
- 1318 + C 7 j0{)2 = 0
c, = 26·3
C 8 is given by
70 = - 45 xX ()-10
4
(0·75) 2 + 26·3ln (Q-0075) + C 8
38 x 100
ONE-DIMENSIONAL STEADY STATE CONDUCTION 35
-1·85- 128·7 + C8
C 8 = 200·6
The equation for temperature is therefore:
t = - 32,900r2 + 26· 3 In r + 200·6
with r in metres.
The maximum temperature occurs at the outer radius. Substituting
r = 0·02 m in the above equation gives
- 5 X J0 4 X 0·0075
------~--
26·3
+- -
2 X 0·38 0·0075
- 494 + 3510
+ 3016
The heat transfer internally is in the direction of negative radius.
hence
Q(-r) = - (-kA:~)
+ 0·38 X (2n X 0·0075) X 3016
53·9 kW/m length
This result may be checked since all the heat generated m the
conductor must be dissipated internally.
Q(-rl = (volume/m length) x q'
= n(0·02 2 - 0·0075 2 ) X 5 X \0 4
= 53·9 kW/m length.
PROBLEMS
I. The walls of a refrigerator for a shop consist of slag wool 0·1522 m thick
sandwiched between sheet iron. 0·0794 em thick, on one side and asbestos
36 ENGINEERING HEAT TRANSFER
board, 0·953 em thick, on the other. The total surface effective for heat transfer
is 37·2 m 2 The atmospheric temperature is 18·3°C and the temperature in
the cold room is - 3·9oC.
The thermal conductivity of iron, slag wooL and asbestos board may be
taken as 69·1, 0·346, 1·21 x 10- 3 respectively and the surface heat transfer
coefficient as 1·705 x 10- 3 ; in kW, m, K, units.
Compute the heat leakage into the refrigerator. (Ans. 0·51 kW.) (King's
College, London}.
2. A spherical container 1·22 m internal diameter is made of sheet metal of
negligible thermal resistance and covered by cork insulation 0·457 m thick.
The interior contains a liquefied gas at - 62·2°C for which a surface heat
transfer coefficient of 1·06 k W/(m 2 K) may be considered to apply. The
atmospheric temperature is 18·3oC. Moisture vapour permeates the cork and
freezes at a suitable position to form an ice barrier. The mean surface co-
efficient for the outside may be regarded as 0·021 kW/(m 2 K). Calculate the
thickness of the ice assuming that the conduction characteristics of the cork
remain constant throughout.
Assume the thermal conductivity for cork is 43·2 X 10- 6 kW/(m K). (Ans.
Q-305 m.) (Queen Mary College, London).
3. A 30 mm diameter pipe at 100oc is losing heat by natural convection to
the atmosphere at 20°C at the rate of 0·1 k W per m length. It is required to
cut down this loss to 0·05 k Wjm. Two insulating materials A and Bare avail-
able. There is sufficient of A to use it at the rate of 3·14 x w- 3 m 3 jm length,
and ofB to use it at the rate of 4·0 x 10- 3 m 3 frn. The thermal conductivities
of A and B are 0·005 and Q-001 kW/(m K) respectively. Js it possible to
achieve the required degree of insulation? Assume the convection coefficient
applicableto the bare pipe is also applicable to the outer surface of insulation.
(Ans. B inside 0·0437 kWjm; A inside 0·0742 kWjm.) (The City University)
6. The inner surface of a 0·23 m furnace wall is at 800°C. The outer surface
convects to the atmosphere at 21 oc, with a coefficient of 0·012 kW/(m 2 K).
The conductivity of the furnace wall is 870 x 10- 6 k W/(m K). To cut down
heat loss, an additional wall 0·23 m thick of insulating brick is added on the
outside. having a conductivity of 260 x 10- 6 kW/(m K). For the same outer
surface coefficient, calculate the percentage reduction of heat loss, the brick
interface temperature, and the brick outer surface temperature. (Ans. 7H %.
633°C, 73°c.)
7. A wide copper strip 10 mm thick carries an electric current of density
50 amp/mm 2 . The heat generated is dissipated from the two wide faces
of the strip by convection. On the left-hand face the convection co-
efficient is 5 kW/m 2 K and on the right it is 10 kW/m 2 K. The surround-
ings are at 25°C.
You may assume for one-dimensional conduction with heat genera-
tion, that
8. The walls of a house consist of two skins of brick each 115 mm thick,
separated by an air gap of 50 mm. The heat transfer coefficients on the
outside and inside walls of the house, respectively, are:
outside inside
by convection 7 W/m 2 K; by convection 3 W/m 2 K
and
by radiation 5 W/m 2 K; by radiation 4 W/m 2 K
38 ENGINEERING HEAT TRANSFER
Across the air gap, the overall heat transfer coefficient is 4 Wfm 2 K,
taking into account conduction, radiation and convection effects.
To reduce heat loss the air gap is filled with insulation, so that heat
transfer across this region is now by conduction only.
(a) For a wall area of 70m 2 calculate, for the uninsulated wall, the
annual heat loss in kWh, assuming an annual mean ambient temper-
ature of l2°C and a mean internal temperature of 22° C maintained
24 hours a day.
(b) Calculate the percentage reduction in annual heat loss across the
walls by the insertion of insulation, of thermal conductivity 0·03 2
W/mK.
The thermal conductivity of brick is 0·81 W/m K.
(Ans. (a) 8066·3 kWh. (b) 63·32 per cent.)
(4.2)
TWO-DIMENSIONAL STEADY STATE CONDUCTION 41
y
t specified at boundary
m'
I t specified at
boundary
or
(4.3)
ot) + (ooy2t)
(ox2
2
0
2
0
= tl + tz + tJ + t4 -
a2
4to = 0 (4.5)
or simply
(4.6)
Q(4-0l = - k(a x 1)
(to - tJ = k(t 4 - t 0 )kW/m (4.8)
a
(4.9)
EXAMPLE4.1
Establish the temperatures a-k in the duct shown, by relaxation, and
calculate the conduction heat transfer through the duct. k for the duct
is 0·1 kW/(mK).
Mesh point: a b c d e f g h j k
Initial residual: -4 -1 +9 +22 +7 +4 +6 +5 -2 0
Operation:
+6atd=31° +15 -2 +13
+4 ate= 45° +3 -1 +2
+3ate=44° +5 +1 +7
+2 at f = 48° +3 -1 +8
+2 at g =49° +1 0 +7
+2 at h = 50° +2 -1 0
+1 at d = 32° 0 +1 +4
+1 ate= 45° +2 0 +2
-1 at a =49° 0 +2
+1 atf= 49° +1 -2 +3
+1 at g = 50° 0 +2 0 +2 +1 -1 -1 0 0 0
d e f g h k
-i-
c
goo .. 1
L- . . ;.;o:~ l __ ----
Fig. Extlmple 4.1. Unique ptlrt ofrecltlllflllW tiMet.
EXAMPLE4.2
The figure shows a hollow square section, with the centre hollow rotated
through 45°. A square mesh is superimposed so that temperatures in the
section may be determined. Then by considering heat flow along mesh
lines terminating at the boundaries, the heat flow into and out of the
section is found, per unit length. A simple program in BASIC is given to
carry out this task. The principles involved may be deduced by studying
the listing. All points are solved to form one complete iteration, and
iterations continue until all points have changed in temperature by less
than 0·1 °C since the previous iteration. In this example this is achieved
in 9 iterations. The average heat flow for a temperature drop of 80°C is
1510·54W/m, so that the actual surface values vary by 6·67 per cent.
This may be attributed to the coarse mesh with the inside surface set at
45° to the mesh lines.
2 3 4 5
""" "6 7 8 9 10
""" """ 11 12 13 14
"""
POl'~T TE11PERATURE
n.2n0r.ooone 02
0.20000DOOE 02
0.2oonnorJoE 02
0.znonoonoe oz
0.200fl0f!'l0E 02
6 0.26?1449~E 02
O.BSS2788E U2
r. 0.4090~473£ 02
~ n.,30~52H4E nz
,,
,n n. 5403'-!'\·J~e G?
0.477:.5155E 02
1~ 0.r>17~5517E 02
1~ C. 771,~7700E G2
,.,
14 0.10000UOOE 03
O.:HJSL.U663E 02
H fl.10000000E 03
tlUMUER OF ~ELAXATIOrJS
ll<AT TRAUSHW !II 0.14097312E 04 IIATTS/~1 LEt/GTII
IIEAT TRANSFER OUT 0.16113510E C4 WATTS/~ LENGTH
TWO-DIMENSIONAL STEADY STATE CONDUCTION 47
The simplest possible case has been considered so far, i.e. that of fixed
boundary temperatures. If these are uniform then the boundary is
said to be isothermal. What happens beyond the boundary to create
the isothermal condition is outside the scope of the problem, and
in this sense the exercise is rather unrealistic. A boundary which is
convecting or radiating, or perhaps is insulated, is the more practical
situation. Mesh points occurring on such boundaries will have
temperature relationships other than (4.6) or (4.7) for points in the
field Also, internal boundaries between different conducting
materials may exist With these and other complexities such as the
boundary shape, the work soon becomes too complicated to be
treated by relaxation methods. However, whatever these complexi-
ties may be, the problem always reduces to solving a set of simul-
taneous equations, and two computer-based methods are available.
The first is that of a direct solution using the Gaussian elimination
method 3 , and the second is an iterative solution. The essentials of the
second method will be described since one basic program with minor
changes may be applied to a wide range of problems.
4
electrical equivalent
fluid
2
electrical equivalent
tl + tJ
:. - 2- + t 2 + (ha/k)tc- (2 + ha/k)t0 = 0 (4.11)
fluid
.....---.----.. 0
~~ tr
I 1 I -
L----l--
I· a .I T h
l_2 _r 2
electrical equivalent
-\ a/2 f-
Fig. 4.4. External boundary corner witlr con11ection.
r----~-
1a
r-
l
I 2
fluid
1 a/2 I
I
I 4 I 1
L ____ L_
I
I
3
~ 3 I
L ____ j
J=1, 10 TF1,HCONI
2
3
'"' I'
!'-,~
"'""~ -~1- l
f- 5
lines of symmetry
Fig. 4.6. Field for an elementary BASIC program-part of a hollow duct.
TFI, TF2, fluid temperatures; HCON I, HCON 2, convection coefficients;
TK, thermal conductivity.
(a)
vd
(b)+~
~w~l I
1 =-
wt
pi
v = s.-pt v
I
(4.15)
where Si is the 'shape factor'. The shape factor for this rectangular
element is wjl. Any other geometrical shape which passed the same
current for the same voltage drop would have the same shape
factor.
A similar equation may now be written for a geometrically similar
element, length L, width W and thickness T, along which heat is
conducting. (}is the temperature difference and k the conductivity,
hence from Fourier's law,
Q = k WT (} = k S TO (4.16)
L q
where Sq = WjL and is the shape factor for the element. As the
elements are geometrically similar, the shape factors are equal.
Dividing (4.16) by (4.15):
(4.17)
Q will be in heat units per unit length and time, with T as unit
thickness. Thus from measurements of I and V and from a knowledge
of(} and k for the conducting problem, the heat ftow may be calcu-
lated. Equation (4.17) is valid for any geometrical shape, provided
the prototype and· model are geometrically similar, when the shape
factors are equal. pjt is a property of the conducting paper and is
supplied by the manufacturers. It has a value of about 2000 ohms
per square. (R = plfwt = p/t for a square, regardless of its size.)
Conducting paper is rolled during manufacture, and this can result
in an anisotropic effect. To account for this let the element in Fig.
4.7b have an equal resistance in both principal directions. Then
p 1 1jwt = p 2 wjlt, where p1 is now the resistivity in the direction of
current flow shown. It follows that ljw = J(p 2 / p 1). The model in
Fig. 4.7a would then be made rectangular to the extent indicated by
measurements of p 1 and p2 •
Structures of composite material having differing thermal
conductivities may be simulated by punching regular holes in the
paper, to produce an area having an increased resistance. Convecting
boundaries may be included either by cutting a band of paper outside
TWO-DIMENSIONAL STEADY STATE CONDUCTION 55
the boundary into strips normal to the boundary, or by adding carbon
resistors Rc at the edge 4 , so that Rc/R = (1/h)/(L/k WT) where his the
convection coefficient.
Three-dimensional shapes of rectangular form may be modelled
using a number of sheets of paper to represent layers in the third
coordinate direction, with additional resistances joining the centres
of corresponding elements. Complicated three-dimensional shapes
may be simulated in electrolytic tanks 5 .
PROBLEMS
1. The diagram shows a plan view of the vertical insulation round the walls
of a liquefied natural-gas storage tank. The inside and outside surface temp-
eratures of the insulation are - 161 oand + 1oc. Calculate the heat transfer
rate into the tank per metre height. Treat points a and b as fixed boundary
temperatures of- 53° and -107°C and assume one-dimensional conduction
through the tank sides. Take k for the insulation as 50 x 10- 6 kW/(m K).
The mesh size is 0·2 m. (Ans. 0·1764 kW/m.) (The City University).
-161°C
I. 2·4m
.I
2. The diagram represents a square section with 81 mesh points.
Boundary temperatures and convection coefficients are subscripted to
allow for varying input information. Write a BASIC program to calculate
temperatures at the 81 points. Note that there are only three separate
physical situations: points 2, 3, 4 and 5 are similar, also points 6, 7, 8
and 9 are similar. Point 1 applies throughout the interior field. To check
the program, the same convection coefficient on all faces will produce a
symmetrical temperature field. Subsequent runs may be made with
alternative boundary conditions.
56 ENGINEERING HEAT TRANSFER
Tl(l)= 10.0
HI(!)
--
T!(J) = 10.0 T2(J) = 10.0
HI(J) + H2(J)
--I
4
T2(1) = 100.0
H2(1)
1 2 3
15
4 9 10 12
to-1 W/(m:tK
11 16 /1
2 0°C /1
/ 14 8
/
/
13
5 6 7
5. The figure shows a unique 1/8th part of a hollow square duct having
inside and outside surface temperatures of 150°C and 15°C. Using the
initial values of temperature given, relax until residuals have values not
greater than± 3. Calculate the heat transfer rate through the duct perm
length b¥ taking the average of the conduction rates at the 15°C, 70° C
and 150 C isothermals. The thermal conductivity of the material of the
duct is 4 Wjm K. (Ans. 4240 W/m.)
15°
T
" "'L4° 35° I 46o I 510 I 520
!"-.
"" 60°
I'.
\.85° 195°
I
J!9
"'."'
I
1
II£
""
150°
I
'
REFERENCES
58
TRANSIENT CONDUCTION 59
(}
---+-----===-----==---
1
I
/
l
j et heating curve
time
t
8
cooling curve
dt
1-
Fig. 5.1. Heating t111d cooling curtJesfor lumped systems.
The heating curve in Fig. 5.1 is the result of an internal heat source,
and the maximum temperature rise (}max is achieved when the rate of
cooling is equal to the rate of heat release. For a volume V and heat
release rate q'
q'V = hA (}max (5.2)
At some temperature(} the energy balance with a density p is given by
pVcP d(} = (q'V - hA(}) dt
:. - t 1 (q'V
= [ - -ln - hA(}) ]6t (5.3)
pVcP hA 0
60 ENGINEERING HEAT TRANSFER
where 01 is the temperature rise at time t.
e-lrr = 1 - hA()t
q'V
Introducing ()max from (5.2) this then gives
01 = ()max (1 - e- 1rr) (5.4)
When()= 0, (5.3) becomes pVcPd() = q'Vdt
. (d(J) q' V hA ()max () _,.
· · dt ll=o = pVcP = pVcP = T
Hence the time constant will appear on the heating curve as shown in
Fig. 5.1. This provides an experimental method for measuring T and
deducing, for example, the convection coefficient, h.
EXAMPLE 5.1
An electric motor is designed to operate on a repetitive load duty
between temperature limits of 30° and 55°C, assuming an ambient of
15°C Heat dissipation on load is 0·38 kW, and the off-load period is
294 seconds. The motor has an effective mass of 3·5 kg, specific heat
of 0·45 kJ/(kg K~ and a surface convection coefficient of 0·15 kW/
(m 2 K). Determine: (i) the cooling area to be provided; (ii) The maxi-
mum temperature that would be achieved if the motor ran indefini-
tely; (iii) The duration of the allowable load period. (The City
University)
Solution. The motor cools from 55° to 30°C in 294 seconds, hence the
time constant may be found: (3(}-15)/(55-15) = e- 1rr wheret = 294
therefore t/T = 0·982 and T = 300 sees. But
T = mcJhA = (3·5 x 0·45)/(0·15 x A)
:. A = (3·5 x 0·45)/(0·15 x 300) = 0-()35 m 2
At the maximum temperature rise, all energy dissipation is convected
away, hence ()max hA = 0·38
.". ()max = 0·38/(0·15 X 0·035) = 72·5
TRANSIENT CONDUCTION 61
The maximum temperature is therefore 72·5 + 15·0 = 87·5°C.
To establish the load period, use equation (5.4) to determine the
times to reach temperature rises of 40° and 15° from ambient Thus:
40 = 72·5 (1 - e-•·' 300) :. t 1 = 241 sees
15 = 72·5(1- e-•z/ 300) t = 69·6secs
:. 2
:: = « ( ;:~) ((2.8))
For the general problem, numerical procedures will be described.
It is necessary to replace equation (2.8) by a finite difference
relationship. Figure 5.2 shows a plane slab uniformly divided into
sub-slabs of thickness a, with a temperature contour at some time t0 •
Recalling the argument of section 4.1 it will be seen that the tempera-
tures t 3 • 0 , t 4 • 0 , and t 5 • 0 are related
(fPt)
ox2 4,0
(5.5)
With a forward time step, the fmite difference relationship for (otfo t)
(at)
IS
= t4,1 - t4,o
Ot 4,0 L1t
where _t 4 • 1 is the temperature at point 4 at time tl' which is At after t 0 •
Equatton (2.8) can now be replaced by
t4,1 - t4,o = oc(t3,o + ts,o- 2t4,o)
(5.6)
L1t a2
I
I
I
I
I
II I lt4o'1
I I I
1-Q-1 sub-slab
I I
EXAMPLE 5.2
The surfaces of a brick wall, 300 mm thick initially 20°C throughout,
rise in temperature at a constant rate of 10°C every 2500 seconds.
Dividing the wall into six equal slabs find the temperature distribu-
tion in the wall after 1<f seconds. Use F = !, a = Q-05 x 10- s m 2 fsec.
TRANSIENT CONDUCTION 63
. t
. - 1,1
= (2k
pcPa2
at) t 2,0
+ (2h
pcPa
at) t t,o
+t
1,0
( 1 _ 2k
pcPa2
At _2hpcPaAt)
:. t 1 , 1 = 2Ft 2 , 0 + 2FBtr,o + t 1 , 0 (1- 2F(l +B)) (5.8)
where B = hajk, which is the Biot number of a sub-slab, and F is the
Fourier number as before. A new stability criterion applies to the
t.
boundary equation, i.e., F(l + B) ~ for the coefficient of t1. 0 to
remain positive. This means that F < t and the numerical procedure
using equations (5.7) and (5.8) is more complicated throughout the
whole field.
EXAMPLE 5.3
120W/m 2 K 15 W/m 2 K
1--
~..--~...--~...-L--L---L--L---L--
25°-500°C
in 200 sec f . - - - ,...-f--,...-f--,...-
f.4-----50 mm - - - - - . J
Fig. Example 5.3.
10 DIME"SIO~ TI11J,TNEWI11)
20 A=0.005
30 TK=0,19
40 DEN=577,0
50 CP=816.0
60 111=120,0
70 112=15.0
80 T1=25.0
90 T2=25.0
100 TTOP=500,0
110 TRT=200,0
120 TLIM=1800.0
130 TPRINT=200.0
140 DTPRI"T=200,0
150 D1•111•A/TK
160 B2=H2•A/TK
170 f•TK/IDEn•CP•A•AJ
180 TIIH1=1,0/I?..O•F•I1.0+01 ))
190 TIME2=1.0/12,0•F•I1,0+B2ll
200 IFITIHE1>TIME2JGO TO 230
210 TIME=TinE1
220 GO TO 240
230 TIME=TIHE2
Z40 PRINT,"TIHE STEP IS:'',TIME
250 PRI~T,"~O YOU WISH TO ROUND OFF AT A SHALLER VALUE OF TIME STEP?"
260 PRINT,"1 - YES; 2 - NO''
270 [~PUT IX%
280 IF!IX%>1JGO TO 300
290 !~PUT TI"E
300 F=F•TitlE
310 F2•2.0•F
320 FB12=2.0•F•01
330 F022=2.0•F•02
340 FT1=1.0-2.0•F•I1,0+B1l
350 FT2=1.0-2.0•F•(1,0+92J
360 FT=1.0-2.0•F
370 DTR=I500,0-25.0l•TIME/TRT
380 TR UII=O, 0
390 T1=T1+DTR
395 IFIT1>TTOPJTHE~ T1=TTOP
400 TNEW!1J=f2•TI2J+F012•T1+111)•FT1
410 FOR 1=2,10
420 THEW!I:J=F•!TIIZ-1J+T!IX+1Jl+TIIZJ•Fr
430 ~IEXT IZ
440 TNEWI11l=F201(10l+FB220J2+T!11)•FT2
450 TRU"=TRUN+TIME
460 FOR 1~=1,11
470 TII%)•TNEWIIZJ
480 NEXT U
490 IF!TRUN<TPRINTJGO TO 390
500 II EAT Ill= ( T 1- T ( 1 ) ) • H 1
510 HEATOUT•!TI11J-T2J•H2
66 ENGINEERING HEAT TRANSFER
520 HSTORE=O.O
530 FOR 1%=1,11
510 HSTORE=HSTORE+(T(I%l-25.0l•A•DEN •CP
550 IIEXT 17.
560 PAifiT,"IIEAT TO DOOR ='',IIEATIH," WATTS"
561 PRI~T,"IIEAT FRO~ DOOR ='',ti£ATOUT,'' WATTS"
562 PRINT,''IIEAT STORED :",HSTORE,'' JOULES''
565 PRINT,"TEMPERARURES -"
570 PRINT T<l:l FOR 1%=1,11
580 TPRINT=TPRINT+DTPRINT
590 IF!TPRINT>TLI11lGO TO 620
600 GO TO 3?0
620 STOP
It may be seen from Fig. 5.2 that the averaging of temperatures can be
carried out by drawing. This is the basis of the Binder-Schmidt
method 1 • 2 • Graphical constructions are also possible for convecting
boundaries. These rely on the fact that the conduction rate at the
boundary must equate to the convection rate. Various cases are con-
sidered by Hsu 3 •
Ot
at
= IX (0ox
2
t
2
+ iJ 2 t)
oy 2
(5.9)
and referring back to the nomenclature of Fig. 4.1, the finite difference
relationship for a field point with a forward time step can be seen to
be
to, 1 - to o = IX (tl,O + tz,o + t3,o + t4,o - 4to,o) ( 5 .lO)
Llt a2
This is re-arranged to give
to,t = F(tl,o + tz,o + t3,o + t4,o) + to,o(l - 4F) (5.11)
with the stability requirement that F ~ ±-
In transient work, alternate finite difference relationships having
backward time steps may be used. The equivalent of equation (5.10)
would be:
to,l- to,o = IX(tl,l + tz,t + t3,1 + t41- 4to,t) ( 5 . 12 )
.D.t a2
Then t 0 ,0 is the only known temperature, and equations for all points
must be solved simultaneously to obtain the temperatures after the
next time step. However, there is no stability restriction in this case.
68 ENGINEERING HEAT TRANSFER
ao = oc(a2o)2 (5.14)
at ax
70 ENGINEERING HEAT TRANSFER
oe = p
~
2
2 Ce-pxsin(2nnt- qx) + pqCe-pxcos(2n:nt- qx)
t•X
+ pqCe-P" cos(2n:nt- qx)- q 2 Ce-P" sin(2nnt- qx)
Hence, equation (5.14) becomes, noting that ce-P" may be cancelled
from all terms:
pq = nn/IX
(} = C sin 2nnt
and the boundary condition gives
(} = (lJm)o sin 2nnt
Thus comparing these two equations shows that C = (lJm)o. The
final solution is therefore
(5.17)
This equation shows that the maximum variation of 8 decreases
exponentially with x, the distance into the solid, according to the
equation
(5.18)
The general form of the result given by equation (5 .17) is shown in
Figs. 5.3 and 5.4. In Fig. 5.3 the temperature variation with distance
distance x
x=O
representation of a semi-infinite solid
00
J
X=O -+--------..-00
00
Fig. 5.3. Established temperature variation vs. distance into solid, at t =J/2n.
t------period 1/n----~
~ Jc:lJ
Fig. 5.4. Temperature variation with time at x = 0, and at depth x.
At=~-
( 1 )0·5 (5.19)
2 nnrx.
At for a complete temperature wave of length X is 1/n hence, for
a complete wave, the wave-length is given by
2n = X(nn/rx.)0 "5
X = 2(nrx./n)0 "5 (5.20)
Also, the velocity of propagation of the temperature wave into the
solid is
U = -X = 2(nnrx.) 0 .5 (5.21)
1/n
The ratio of maximum temperature variations about the mean may
be compared using equation (5.18). Thus:
(0 )
= exp[- x(nn/rx.) 0 "5 ]
(0m)o
m "
TRANSIENT CONDUCTION 73
Jf it is required to determine the distance x at which (Om)" has
decreased to a certain percentage of (Om)0 , this result may be re-
arranged to give
(5.22)
oO)
(OX = - (Om)o{1tn/oc)0 "5 (sin 21rnt + cos 2nnt)
x=O
it follows that
( ~0 )
uX x=O
= - (Om)o{27tn/oc)0 "5 sin(27tnt + 7t/4)
and hence,
q = k(Om)o{27tn/oc) 0 "5 sin(27tnt + 1t/4) (5.23)
From equation (5.23) it will be seen that the surface heat transfer
rate varies sinusoidally and with the same frequency as the surface
temperature, but leading by a period of 1/Sn. The total heat transfer
at the wall is given by
Jq dt = Jk(Om) {27tn/oc)
0 0 "5 sin(27tnt + 7t/4) dt
= -k(Om)o{1/27tnoc)0 "5 cos(27tnt + 7t/4) heat/unit area (5.24)
Thus the energy stored, as represented by an integral of heat transfer
rate at the surface, also varies with the same frequency, but it will
be found that it lags behind the surface temperature variation by a
period of 1/Sn. Further, it will be apparent that the surface heat
transfer is both to and from the solid and that the energy stored is
74 ENGINEERING HEAT TRANSFER
in sequence both positive and negative relative to the mean
temperature.
EXAMPLE5.4
In a cyclic heating process the inside of a furnace wall is subjected to
a sinusoidal temperature variation. The temperature rises from a
minimum of 100°C to a maximum of 750°C in 3 hours. With k =
692 x 10- 6 kW/(m K) and ex = 0·0516 x 10- 5 m 2 /sec, determine:
(i) the velocity of the temperature wave penetrating the wall;
(ii) the time lag of the wave function at a depth of0·2 m compared with
the surface; (iii) the maximum and minimum temperatures at a depth
of 0·2 m; (iv) the temperature at a depth of 0·2 m when the surface
temperature is a maximum; and the surface temperature when the
temperature at a depth of0·2 m is a maximum (The City University).
PROBLEMS
1. Steel strip of thickness 1·27 em emerges from a rolling mill at a temperature
of 538°C and with a velocity of 2·44 mjsec. The strip is cooled in such a way
that its surface temperature falls linearly with distance from the mill at a rate
of ll0°Cjm.
Derive a finite difference method for dealing with this case oftransient heat
conduction assuming that heat flows only in the direction normal to the strip
faces.
Subdividing the strip into six increments of thickness, determine the tem-
perature distribution in the strip and the heat flux from the surface at a
position 2·74 m from the mill. (For steel take thermal conductivity 43·3
x 10- 3 kW/(m K), thermal diffusivity 0·98 x 10- 5 m 2 js.) (Ans. 3·41 x 103
kWj(m 2 )(University of Manchester).
4. Given the differential equation (otfor:) = rx.(o 2 tfox 2 ), for unsteady conduc-
tion in a 'one-dimensional' wall, show that the temperature t•.p+ 1 at some
section n and time instant (p + 1) can be calculated approximately from
76 ENGINEERING HEAT TRANSFER
600 x w- 6
-~---~----
A 1500 0·84
B 600 X 10- 6 1200 1·60
c 280 X 10- 6 750 1·10
The screen is 5·3 em thick and is divided into 5 increments. It is initially l5°C
throughout; the temperature of one face rises linearly by 20oc per minute.
Regard the other face as insulated. Determine, for the chosen material, and
by a numerical technique, the time before the insulated face temperature
starts to rise, and the temperature of this face after 18 minutes. (Ans. B,
15 mins, 24·4°C) (The City University).
10. An internal combustion engine runs at 2500 r.p.m. The thermal diffusivity
of the carbon steel ofthe cylinder walls is H6 x 10- s m 2 /s. The temperature
of the cylinder wall varies sinusoidally between 5000 and l00°C. Assuming
that the cylinder wall behaves like a semi-infinite solid, determine the depth
into the wall in em at which the temperature amplitude has decreased to I per
cent of the surface value, and plot the heat transfer rate at the wall surface
over a complete cycle. k = 40 x 10- 3 kW/(mk). (Ans. 0·194cm, limits
±2·68 x 104 kW/mz.)
REFERENCES
1. Binder, L. Dissertation, Mtinchen (1910).
2. Schmidt, E. Festschrift zum siebzigsten Geburtstag August Foppl, Springer,
Berlin (1924).
3. Hsu, S. T. Engineering Heat Transfer, D. Van Nostrand Company, Inc.,
New.York, 103, (1963).
4. Bayley, F. J., Owen, J. M. and Turner, A. B., Heat Transfer, Nelson (1972).
6
Forced convection: boundary
layer principles
6.1 Introduction
In Chapter 1 it has been shown that to evaluate convection heat
transfer, the magnitude of the coefficient h in Newton's equation
has to be found. The study of convection centres round the behaviour
ofthe fluid flowing past a surface, and the subject matter divides itself
under various headings concerned with the type of flow situation
or the 'method of analysis. This chapter shows how the convection
free stream, v5
laminar transition • : • turbulent
:---------'-====---;--= I
-----x
pv.x
--~
5 X lOs (6.1)
Jl.
where Jl. is the coefficient of molecular viscosity. The boundary layer
exists as a result of the action of viscous shear within the fluid, the
shear stress being proportional to the velocity gradient
dvx
T =Jl. dy (6.2)
-Et-~~tl~~:>f)-8-
f--- starting length Ideveloped
fully fully
developed
laminar turbulent
flow flow
Fig. 6.2. BoiUUimy layer growth in a tube.
+,
/.}/"
I ' .........
I ·--..... ·--......
-"-~ ...........
~t:
(b)
heanransfer from
fluid to waU
(a)
heat transfer from
wall to fluid
Fig. 6.3. Exlllllpk!l oftemperlltlll'e bollllllary layer!l.
pvxdy
and
FORCED CONVECTION: BOUNDARY LAYER PRINCIPLES 81
respectively and hence the net flow into the element in the x-direction
is
- povxd
- X dy
OX
Similarly the net flow into the volume in the y-direction is
OV
- p.::...:.l.dydx
oy
The total net flow in must be zero, hence
- p(OVx
ox + ~)dxdy
oy =0
~~+~=0 (6.5)
ox oy
The equation of momentum arises from the application of New-
ton's second law of motion to the element, assuming the fluid is
p dy
Fig. 6.4. Elemellt of bo11111lary lllyer for co11tiludty flllll mo~Uittllm IHUfulee.
82 ENGINEERING HEAT TRANSFER
-kd ~ at a2t )
-kdy ( -+-dx
!Yax ax ax2
-kdx.1..!..
ely
p.dx[av, + ~(av,)dy]
ay ay ay
so that the net viscous shear in the x-direction is p.dx(o 2 v)oy 2 )dy.
The pressure force on the left face is pdy, and on the right - [p +
FORCED CONVECTION: BOUNDARY LAYER PRINCIPLES 83
(opjox)dx] dy giving a net pressure force in the direction of motion
of -(opjox)dxdy. Equating the sum ofthe net forces to the momen-
tum flow out of the control volume in the x-direction gives, after
neglecting second-order differentials and using the continuity
equation: ( avx OVX) a2 vx op 6 6)
p vx ax + vy oy = ll iJy2 - ox ( .
- [pcP(vY:; + ~t + ~; :; dy)]dxdy = 0
Using the continuity equation and neglecting the second-order terms:
at at ( iJ 2t iJ 2t)
Vx iJx + Vy iJy = IX OX2 + iJy2 (6.7)
at at (o 2t)
vx iJx + vY iJy = IX iJy2
v is the kinematic viscosity or momentum diffusivity, p.jp, and
vjiX = (p./p)/(k/peP) = p.c~k, which is called the Prandtl number, Pr.
If v = IX, then Pr = 1, and the pair of equations will lead to identical
non-dimensionalised solutions of vx and t as functions of y. The
Prandtl number is the ratio of fluid properties controlling the
84 ENGINEERING HEAT TRANSFER
velocity and temperature distributions, and it can vary between
around 4 x 10- 3 for a liquid metal to the order of 4 x Hf for a
viscous oil.
6.2.2 The Integral Momentum and Energy Equations of the Laminar
Boundary Layer
To consider the motion in the boundary layer, an elemental control
volume is chosen that extends from the wall to just beyond the limit of
the boundary layer in they-direction, is dx thick in the x-direction, and
has unit depth in the z-direction. This is shown in Fig. 6.6. An equation
is sought which relates the net momentum outflow in the x-direction to
the net force acting in the x-direction.
fo pv~ dy
6
J"
d opvxdydx
dx
This is the difference between the fluid leaving across face CD and
entering across face AB. The fluid entering across face BD has a
velocity vs in the x-direction, hence the flow of momentum into
the control volume in the x-direction is
Vsdx
d J" PVx dy dx
0
J"
d o pv; dy dx - vsdx
dx d o PVx dy dx J"
Pressure forces will act on faces AB and CD, and a shear force
will act on face AC. There will be no shear force on face BD since
this is at the limit of the boundary layer and dvx/dy = 0. The net
force acting on the control volume in the x-direction will be
Px(> - ( px + -dpx )
dx dx (> - T
w
dx = - dpx dx - r dx
~
dx w
(6.8)
dx
d J" PVx(Vs -
0 vx) dy = Tw (6.9)
Fluid will also enter the control volume across face BD at the rate
d fy.
dx o PVxdydx
Again this is the difference between the flow rate out at face CD
and in at face AB. The enthalpy flow will be
-k dx( 01 )
oy y=o
FORCED CONVECTION: BOUNDARY LAYER PRINCIPLES 87
For conservation of energy:
d
cpts -d
X
J"" PV.x dy dx - -dd JY· pcPtv.x dy dx -
0 X 0
k dx (ot)
,--
uy y=O
= 0 (6.10)
Beyond the limit of the temperature boundary layer, the temperature
is constant at t., and hence the integration need only be taken up
to y :r <5 1• Equation (6.10) therefore gives
d
-d-
J,) pv.x(v. - v.x) dy = <w
X 0
88 ENGINEERING HEAT TRANSFER
140 v
(Jd(J = -·-dx
13 v.
On integration
(J 2 l40vx
-=--+C
2 13v.
C = 0, since (J = 0 at x = 0
(J 2 = 280vx
l3v.
or
4·64
(6.13)
X RetX
This result, due to Pohlhausen,4 is required later on in considering
the integral energy equation.
The temperature distribution in the boundary layer is assumed to
follow a similar law to the velocity distribution. Thus:
(} = (t - tw) = dy + ey2 + fy3
where, again, d, e, and f are constants. The boundary conditions
are that at y = b1 (the thickness of the temperature boundary
layer), (} = 0, and also (il(}joy)t}, = 0. Also, from equation (6.7) it
follows that (o 2 0joy 2 )y=o = 0 because vx and vY are both zero at
y = 0. From these conditions it follows that
e
o.
= 0, != - 2(J3
I
and hence
(6.14)
FORCED CONVECTION: BOUNDARY LAYER PRINCIPLES 89
Turning to the integral energy equation, the substitutions
e= (t - tw) and e. = (t. - twl are made to give
j_fb, (e. -
dx o
e)vx dy - a (:e)y y:O
= 0 (6.15)
ae)
(oy 3e.
a y;O = a2b 1
This result is substituted in equation (6.15) together with the ex-
e
pressions for and vx to give:
3e.
= a2).J
This then leads to
--=a--
d 3.A.J.
dx 20
3
2.A.Jv.
d 2 lOa
dx (). b) = v•.A.J
22 e+
dx
22 dJ
dx
= 1oa
v.A.J
(6.16)
( ~~1 ) 3 = ()-93
Pr
+M
xt
noting that oc/v = Pr, and M is a constant of integration. The thick-
ness of the thermal boundary layer will be 0 at the beginning of the
heated section, at x = xh, say, and hence
M = _ 0·934
Pr
and finally :
hence ~~ = (<>-93)t
Pr
}
i = p~t
(6.17)
or, approximately,
r~u X
Fig. 6.B. Ltuniluu forced cOIIW!Ctioll 011 a flat p/tlte. Helltillg COIIIIIIeiiCel .lia
from tl¥ letuling ··dge.
FORCED CONVECTION: BOUNDARY LAYER PRINCIPLES 91
Using the equations for !50 !5, and the value of(o8joy)y=o it is possible
to determine the heat transfer at the wall, since
~--t, ~
lt,.-.i
Fig. 6.9. To illustrate the significance of the Nusselt number.
how ~5; is defined. The heat transfer at the wall is qw = h8, and may
be expw;sed as qw = (k/!5;)8,. It follows that h = k/15; and hence
hx x
Nu=k=-;5;
1i 1
=-
fx h dx
X 0
92 ENGINEERING HEAT TRANSFER
where
h = 0·332k ( ~~r Pr-i, from equation (6.18)
fi k 0· 332
= - - - --
(v·x) t Prt
X f V
ExAMPLE 6.1
Air flows at 5 m/s along a flat plate maintained at 77°C. The bulk
air temperature is 27°C. Determine at 0·1 m ,0·5 m and 1·0 m from
the leading edge, the velocity and temperature boundary layer
thicknesses, and the local and average convection coefficients. Use
mean properties of air from Table A6.
{Jmm
I I
· dx1
;,
Fig. 6.10. Element of laminar flow in a tube.
The sum of this and the right-hand side of (6.20) is zero, hence on
cancelling terms,
~~(/')=peP~ (6.2.1)
vr or or k ox
This is the energy equation for laminar flow in a tube. With a con-
stant wall heat flux qw, and constant fluid properties, the temperature
of the fluid (at any radius) must increase linearly in the direction of
flow. so that
ot
ox = constant
~=
v.
1- (!_)2
rw
where r w is the wall radius.
This result must be substituted into equation (6.21) before inte-
grating. Hence, after re-arrangement,
~ (r01 ) = ~a ~v.
or or ox
[1 - (!_)
rw
]r
2
rot = ~ ~v. ( r2 -
or a OX 2 4r;
_c_) + c 1
FORCED CONVECTION: BOUNDARY LAYER PRINCIPLES 95
and after a second integration, gives
ot ('- 2 - -r 4- )
t = -1 -va + C 1 In r + C 2 (6.22)
ex ox 4 16r!
c2 = t
w
-!CX ~v
OX a
3';
16
t = ~v r [!(!._)
!cxoxaw4rw
2
2
- _!_(!._)
16rw
4
- 2_]
16
+ t
w
(6.23)
(6.25)
The equation for heat transfer at the wall may be obtained by con-
sidering the temperature gradient at r = rw. Thus, from equation
(6.25)
96 ENGINEERING HEAT TRANSFER
and
(6.26)
(}
r- 2nr dr pvc/}
=-o_ _ _ __
q = _ k(d(}) = 4k(}a = 4k . 72 (}
w dr rw 3r w 3r w 44 m
This is now equivalent to
FORCED CONVECTION: BOUNDARY LAYER PRINCIPLES 97
Hence
h=~.n
3rw 44
and
hd 4k 72 2rw
Nu 11 =-=-·-·-
k 3rw 44 k
8 72
= 3.44
= 4·36 (6.27)
These results are independent of Reynolds number because, for
fully developed flow, the boundary layer thickness is equal to the
tube radius.
ExAMPLE 6.2
Water at a mean temperature of 40°C flows at a mean velocity of
o-1 m/s in a 3 mm bore tube having a constant wall flux of 1-Q k W jm 2 •
Determine the temperature of the water as a function of radius,
using equation (6.23). Use fluid properties from Table AS.
at =
.·. ox qw1tdj(mass flow x cP)
1·0 X 1t X 3 4 X 106
= 4·178 X 1000 X -1t-X-32,.--x_0_·_1_x_9_9_4-·6
21tr 2 v
= 1trw2 va - ~
4 = 1trw2 va ;2
:. V8 = 2vm
With IX= 15·1 x 10- 8, equation (6.23) becomes:
15·1 10 4 rw
2
t = 108 X 3·22 X 0·1 X 2 X 1.526 [!(!_) - _!_(!_)
16 'w
4
- ~]
16
+t w
2
= 9·6[!(!_) - _!_(!_)
4 r.., 16 rw
4
- ~]
16
+ tw
At r = rw> t = tw; at r = 0, t = tw - 1·8°C.
PROBLEMS
1. Derive the heat flow equation of the boundary layer
d f'
dx o (0. - O)U dy = rx (d(J)
dy o
and apply this equation to 'slug' flow of a liquid metal along a plate of uniform
temperature to find the thickness of the temperature boundary layer. (J would
be the liquid metal temperature relative to the plate temperature. Assume
that the temperature ptofile in the boundary layer can be described by an
equation of the form
0
~=
j[ 2x
(4- x)R"
2
J
FORCED CONVECTION: BOUNDARY LAYER PRINCIPLES 99
Show that for a liquid metal of P = 0·01 the temperature boundary layer
thickness is approximately equal toM. (University of Bristol).
Show that the temperature profile can be put into dimensionless form as
~ = .! = 1 - ~(!..)2 + !(!..)4
to - tw Bo 3 R 3 R
where t, 10 , and tw are the local, axial, and wall temperatures respectively,
and R is the radius of the tube. Also show that the Nusselt number
qwd 8
80 k 3
Explain, by writing down the initial equations, how you would derive the
Nusselt number qwdf(Jmk, where (Jm is the bulk temperature of the fluid relative
to the wall. (University of Bristol).
3. Show that if a flat plate has a heated section commencing at xh from the
leading edge, the local Nusselt number at distance x from the leading edge,
(x > xh), is given by:
Nu" = 0·332 Re! Prt(I - (xh/x)*)-t
Determine the velocity and th~rmal boundary layer thicknesses and the local
heat transfer rate at 1 m from the leading edge of a plate heated 0·5 m from
the leading edge, for air at 27aC flowing over the plate at 0·5 mjs, if the
temperature of the heated section is 127°C. (Ans. ~ = 0·0298 m, ~. =
0·0243 m, 0·184 kWjm 2 .)
4. The velocity in the boundary layer of a stream of air flowing over a flat
plate can be represented by
~ = ~(~)- ~(~r
where U is the main stream velocity, u the velocity at a distance y from the
100 ENGINEERING HEAT TRANSFER
flat plate within the boundary layer of thickness b. The variation of boundary
layer thickness along the plate may be taken as
b/x = 4·64(Re,rt
If the plate is heated to maintain its surface at constant temperature show
that the average Nusselt number over a distance x from the leading edge of
the hot plate is
Nu = 0·66(Pr)t(Re,.,>t
(University of Leeds).
REFERENCES
1. Bayley, F. J., Owen, J. M. and Turner, A. B. Heat Transfer, Nelson (1972).
2. Klinruin, T. von, Z. angew. Math. u. Mech., Vol. 1, 233 (1921).
3. Eckert, E. R. G. and Drake, R. M. Analysis of Heat and Mass Transfer,
McGraw-Hill, New York (1972).
4. Pohlhausen, K. Z. angew. Math. u. Mech., Vol. 1, 252 (1921).
7
Forced convection: Reynolds
analogy and dimensional analysis
dv
-r, = pe dy (7.2)
101
102 ENGINEERING HEAT TRANSFER
When turbulent flow exists, the viscous shear stress is also present
which may be added to <1 • The total shear stress in turbulent flow
is thus
dv
• = p(v + e) dy (7.3)
Cd = 0·455 (7.6)
(log Rexf·SS
Equation (7.6) is an empirical relationship, 5 which takes into account
the laminar and turbulent portions of the boundary layer.
The ratio of the velocity at the limit of the laminar sublayer to
the free stream velocity is also of importance, as will be seen later;
this is a function of the Reynolds number at x from the leading edge:
vb 2·12
v, (Rex)O·l
(7.7)
FORCED CONVECTION: REYNOLDS ANALOGY 103
and vb 2·44
(7.10)
Vm = (Red)t
q k dt
(7.13)
T ,u dv
This gives the ratio of qj-r at some arbitrary plane in the flow.
Noting that qj-r has the same value anywhere in the y-direction,
it is possible to express qw/Tw at the wall in terms of free stream and
wall temperatures and velocities.
Thus
(7.14)
v.
Fig. 7.1. Velocity and temperature distributions for the simple Reynolds
tuUllogy.
Equations (7.14) and (7.15) for laminar and turbulent flow are
clearly identical if Pr = 1, i.e., if cP = k/ Jl., or Jl.Cp/k = 1. Re-arranging
equation (7.15) gives
where e. = (t. - tw), and where h is the surface heat transfer co-
efficient.
Substituting the skin friction coefficient Cf gives
or
h Cf
(7.16)
pv,cP = 2
This is one form of the result obtained from the simple Reynolds
analogy; it gives the convection coefficient h in terms of the skin
friction coefficient Cf. h/pv,cP is the Stanton number St. It is the
N usselt number divided by the product of the Reynolds and Prandtl
numbers. Further re-arrangement is possible; for example, con-
106 ENGINEERING HEAT TRANSFER
sidering laminar flow at dist~nce x from the leading edge of a flat
plate, both sides of(7.16) are multiplied by xfk to give
hx Cf pv.xcP
T=T-k-
But cpJ.L/k = 1, or cP/k = 1/J.L, hence
hx Cf pv.x
---
k 2 J.l
or
Cf
Nu X =-Re
2 X (7.17)
for laminar flow on a flat plate. This result may be compared with
equation (6.18) obtained by consideration of the integral boundary
layer equations. If Pr = 1 in this equation then the result is
Nux = 0·332(Rex)t
Reynolds analogy may also be applied to flow in tubes, and for this
o.
purpose and v. in the above analysis may be replaced by the mean
values ()m and vm, since the velocity and temperature distributions
are identical. The linear dimension is now the diameter of the tube, d.
The relationship will be
hd Cf pvmd
---
k 2 J.l
or
(7.19)
(7.21)
v.
and hence
(7.22)
and
(), vb
1 + -(Pr- 1)
v,
and eliminating (tb - tw)/vb between this result and equation (7.22)
gives
Pr
Jl.Vs 1 + vb (Pr - 1)
v,
(7.23)
v.
For turbulent flow on flat plates, both sides are multiplied by x/k
and 11 is introduced to the right-hand side to give
qwx pv.x Cp/1 Cf 1
O.k = 7
k 2 . 1 + vb (Pr _ 1)
v.
Cf RexPr
Nux = 2 ----'-'---
1 + vb (Pr - 1)
v.
Also, for turbulent flow on flat plates, equations (7.5) and (7.7) are
introduced to give t
0·0292RexPr
NUx = 1 + 2·12Re;io(Pr - 1) (7·25 )
Nu 11 = Cf. Re11Pr
2 vb
1 + -(Pr- 1)
Vm
EXAMPLE 7.1
Compare the heat transfer coefficients for water flowing at an
average fluid temperature of 100°C, and at a velocity of 0·232 m/s
in a 2·54 em bore pipe, using the simple Reynolds analogy, equation
(7.20), and the Prandtl-Taylor modification, equation (7.28). At
100°C, Pr = 1·74, k = 0·68 X w-J kWj(m K), and v = 0·0294 X
10-s m 2 js.
The first answer is thus 13·5 per cent lower than the second, which
may be assumed more correct. This solution is for flow in smooth
pipes.
then
</J 2( 1T: 1 • 1T: 2 • · · · ' 7t(n- k)) = I
Each n term will be composed of the Q variables, in the general
form
1T: = Q~Q~Q3 ... Q~,
and will be dimensionless. The set of n terms will include all inde-
pendent dimensionless groupings of the variables. No n term can
be formed by combining other n terms. A set of equations for
a, b, c, ... , x is obtained by equating the sum of the exponents of
each independent dimension to zero. This will yield k equations
for n unknowns. One method of solution is to choose values for
(n - k) of the exponents in each term. The selected exponents must
112 ENGINEERING HEAT TRANSFER
be independent, which can be shown to be true if the determinant
formed from the coefficients of the others does not vanish.
An alternative procedure is to select k of the Q variables and to
combine them in turn with each of the other (n - k) Q variables.
The selection k of the Q quantities must together involve all the
independent dimensions, but they must not form a dimensionless
group by themselves. Further, each of the (n - k) Q variables in
each n term is given the exponent 1. This facilitates the algebra,
as will be seen, and is allowable since it only amounts to reducing
the n term by some unknown root. Thus, if there were six Q variables
and four independent dimensions, the two 1t terms would be:
7tt = Q~·Q~•Q)•Qi•Qs
7t2 = Q~2Q~2Q32Qi2Q6
(L)a'
T
( L)"'(_!!__).'( M)d'__!!_
LTO LT UTO
_= r
Nud 1·86 LRed Pr :;lJ (/J.)0-14
(~t
IJ.w (ref. 12) (7.30)
As the entry length x in the Graetz number increases, the value of Nud
approaches the constant value of 3·66. This compares with the value of
4·36 (equation 6.27), deduced from an assumed parabolic velocity
distribution in fully developed flow.
In recent years, fresh analyses of laminar flow have been undertaken
by means of numerical solutions of the momentum, energy and conti-
nuity equations. The work of Collins 15 has shown a close agreement
with the results of Test 13 , equation 7.31, and a considerable divergence
from Sieder and Tate, equation 7 .30.
116 ENGINEERING HEAT TRANSFER
EXAMPLE 7.2
A heat exchanger thermal wheel consists of a porous disc of available
disc area A, having a porosity P. The mass flow to pass through the disc
is m, and the flow area A x P consists of N holes of diameter d. The
flow of gas is laminar through the holes and it is required to find the
group hAh/(mCp) for the wheel,for a range of hole sizes, where his the
laminar heat transfer coefficient, and Ah is the area for heat transfer.
(This group is the NTU, Number of Transfer Units; see also Chapter 12
and Example 12.3.)
andm =pPA V
V =m/(pPA)
m
andRe= d/(P A ll)
Ah = 4PA L/d
In the following program listing, the heat transfer coefficient is calcula·
1
ted from Nud = 1·86(Re Pr d/LP, and WD = the wheel thickness L,
W = the mass flow rate, CP = the fluid specific heat, DEN= the fluid
density p, PR = the fluid Prandtl number, TK =the fluid thermal con-
ductivity, VIS = the fluid viscosity, and D is the hole diameter, d.
180 GZ=RE•P't•D/WD
190 MOL=P•A•4.0/(3.1417•D•D)
200 LPRINT " Matrix Hole Diameter=";
210 LPRINT USING "#I.#I##";D;
220 LPRINT " m, ~umber of Holes =";
230 LPRINT USING "#I##II":MOL
240 LPRI~T" Reynol1s Number ".,
250 LPRINT USING "#####I.#";RE;
260 LPRI~T " Graetz Nu~ber =";
270 LPRINT USING "II##U.U";Gz
280 LPRINT" Heat Transfer Coefficient :";
290 LPRINT USING "I###I.#";H;LPRINT " W/m2 K"
300 LPRINT " Heat Transfer Area =";
310 LP~INT USING "#####.#";AREA;
323 LORINT " mz, NTU =";
330 LPRINT USING "#II.#";TU
340 LPRINT
350 D=D+n.rnn5
360 IF(D<0.3035)GO TO 14n
370 STOP
where n = 0·4 for heating, and 0·3 for cooling, and Re11 > 10,000.
This is for fully developed flow, i.e., (xjd) > 60, and all fluid pro-
perties are at the arithmetic mean bulk temperature.
For both larger temperature differences and a wider range of
Prandtl number:
(7.34)
In this equation (}7 < Pr < 16,700, and all other details are as
before, with J.lw taken at the wall temperature.
In many situations where tube lengths are relatively short, fully
developed flow is not achieved, and the following relationship may be
used:
Nud = 0·036 Re/" 8 p,t (d/L)o-oss (7 .35)
velocity
profile
Fig. 7.3. Normalized temperature and velocity profiles for flow in a tube at
very low values of Pr.
EXAMPLE 7.3
Freon at a mean bulk temperature of -1oac flows at 0·2 m/s in a
20 mm bore pipe. The freon is heated by a constant wall heat flux
from the pipe, and the surface temperature is 15°C above the mean
fluid temperature. Calculate the length of pipe for a heat transfer
rate of 1·5 kW. Use fluid properties from table AS.
Pr = 4·0, J.t = 31·6 X 10- 5 Pas. At +5°C, J.t = 28·8 X 10- 5 Pas.
A comparison of results using equations (7 .33) and (7 .34) may be
obtained. Re = 20 x 0·2 x 105 /1000 x 0·0221 = 18,100. Therefore
Re 0 "8 = 2547. Pr = 4·0, hence Pr 0 "4 = 1·74 and p,t = 1·588.
(J.t/J.tw} 0 " 14 = (31·6/28·8} 0 "14 = 1·013
From equation (7 .33),Nud = 0·023 x 2547 x 1·74 = 102·0
From equation (7 .34),Nud = 0·027 x 2547 x 1·588 x 1·013 = 110·7
Using the second result, which is 8·5 per cent larger than the first,
h = 110·7 X 72·7 X w- 6 X 103 ;20 = 0·402 kW;(m 2 K)
The pipe length required is calculated from Q = ndLh(tw - t1)
where tw and t 1 are the wall and fluid temperatures, hence
120 ENGINEERING HEAT TRANSFER
L = 1·5j(1t X 20 X 10- 3 X 0·402 X 15) = 3·96 m
PROBLEMS
I. The expression, Stanton number = ! x friction factor, may be derived
from the simple Reynolds analogy. Briefly explain this analogy, discussing
any assumptions made and stating limitations to the application of the above
expression.
Air at a mean pressure of 6·9 bar and a mean temperature of 65·5oc flows
through a pipe of 0·051 m internal diameter at a mean velocity of 6·1 mjs.
The inner surface ofthe pipe is maintained at a constant temperature and the
pressure drop along a 9·14 m length of pipe is 0·545 bar. Determine: (a) the
Stanton number, and (b) the mean surface heat transfer coefficient. (Ans.
0·288, 12·56kW/(m2 K.) (University of London).
~ ~[I + a(~r -
= 1)]
which gives the heat transfer per unit area and time, H, in terms of the drag
force per unit area, F, and in which Pr denotes the Prandtl number CJJ./k;
the other symbols having their usual meaning. (a= vJv•. )
Use the Taylor-Prandtl equation to show mathematically the following
deductions, and explain them in simple terms:
(a) For gases the Taylor-Prandtl equation approximates closely to the
8. Air at mean conditions of 510°C, 1·013 bar, and 6·09 mjs flows through
a thin 2·54 em diameter copper tube in surroundings at 272°C.
(a) At what rate, per metre length, will the tube lose heat?
(b) What would be the reduction of heat loss if 2·54 em of lagging with
k = 173 x 10- 6 kW/(m K)were applied to thetube?Take N 4 = 0·023(R4) 0 . 8
P 0 . 33 with all the properties taken at the bulk air temperature. Assume the
surface heat transfer coefficient from the outside of the unlagged and lagged
122 ENGINEERING HEAT TRANSFER
tube to be 17-Q and 11-3 x 10- 3 kW/(m 2 K) respectively. (Ans. 0·174 kW/m,
32 per cent.) (University of Bristol).
10. Explain and derive the simple Reynolds analogy between heat transfer
and fluid friction. Outline the Prandtl-Taylor modification to the simple
theory.
2·49 kg/s of air is to be heated from 15 to 75oc using a shell and tube heat
exchanger. The tubes which are 3·17 em in diameter have condensing steam
on the outside and the tube wall temperature may be taken as 100°C. Specify
the number of tubes in parallel and their length if the maximum allowable
pressure drop is 12·7 em of water.
Assume that/= 0·079Re-t and that the air has the following properties:
density 1-123 kgjm 3, kinematic viscosity 1·725 x 10- 5 m 2 js. (To solve this
problem, see also chapter 13.) (Ans. 94 tubes, 3·75 m long.) (University of
Leeds).
REFERENCES
385-393 (1968).
14. Hausen, H. VDIZ, No.4, 91 (1943).
15. Collins, M. W. 'Finite difference analysis for developing laminar
flow in circular tubes applied to forced and combined convection',
Int. J. Num. Meth. Eng., Vol. 15,381-404 (1980).
16. McAdams, W. H. Trans. AIChE, Vol. 36, 1 (1940).
17. Chapman, A. J. Heat Transfer, 3rd ed., The Macmillan Company,
New York (1974).
18. Lubarsky, B. and Kaufman, S. J. NACA Tech. Note 3336 (1955).
19. Seban, R. A. and Shimazaki, T. T. Trans. ASME, Vol. 73, 803
(1951).
8
Natural convection
I
direction of
induced motion
Fig. 8.1. Natural convection boundary layer on a vertical flat plate.
124
NATURAL CONVECTION 125
from a vertical plate, it is possible to obtain a solution ofthe bound-
ary layer equations of motion and energy, if a body force term is
included. This approach is limited in general application and, conse-
quently, the method of dimensional analysis will be used.
(p. - p)g
and Ps is related to p by
Ps = p(l +PO)
where P is the coefficient of cubical expansion of the fluid. Thus
the buoyancy force is
[p(l + po) - P]g = pgpO (8.1)
JI 1 ,lI
Horizontal cylinders
J
Vertical surfaces
(8.4)
128 ENGINEERING HEAT TRANSFER
when 109 < GrdPr < 10 12 (turbulent flow). Below GrdPr = 104 ,
it is not possible to express results by a simple relationship, and
ultimately the Nusselt number decreases to a value of 0·4. At these
low values of GrdPr the boundary layer thickness becomes appre-
ciable in comparison with the cylinder diameter, and in the case
of very fine wires heat transfer occurs in the limit by conduction
through a stagnant film. Fluid properties are evaluated at the average
of the surface and bulk fluid temperatures, which is the mean film
temperature. If the surface temperature is unknown, a trial and
error solution is necessary to find h from a known heat transfer rate.
Nu 1 = 0·129(Gr,Pr) 0 ·3 3 (8.6)
when Gr1Pr > 108 (turbulent flow). Thus turbulent flow is possible
in this geometrical arrangement.
The converse arrangement is the hot surface above a cold fluid,
or hot surface facing downwards, and a hot fluid above a cold
surface, or cold surface facing upwards. In either case, it is obvious
that convective motion is severely restricted since the surface itself
prevents vertical movement. Laminar motion only has been
observed, and the recommendation is
Nu 1 = 0·25(Gr1Pr) 0 "25 (8.11)
when Gr1Pr > 10 5 . Fluid properties are again taken at the mean
film temperature.
8.3.4 Additional Geometries
For two vertical parallel plates at the same temperature, convecting to
air:
130 ENGINEERING HEAT TRANSFER
where lis the height of the cylinder and dis the diameter. In the case of
small rectangular solid objects there are generally insufficient vertical
and horizontal surfaces to be treated separately, and the natural cooling
of such objects in closed or open environments is a fairly common
occurrence. The interaction of separate boundary layer growth means
that an overall characteristic length is less than either of the actual
principal dimensions, so that
where the characteristic length is x, and /h and lv are the principal hori-
zontal and vertical dimensions, respectively. Then:
-
Nux = 0.6 ( Grx Pr)0·25 (ref. 8) (8.16)
for 1if < GrPr < 109
8.3.5 Approximate Formulae for use with Air
A great deal of natural convection work involves air as the fluid medium
and the fluid properties of air do not vary greatly over limited tempera-
NATURAL CONVECTION 131
ture ranges. Thus it is possible to derive simple formulae from equations
(8.3) to (8.6) and (8.9) to (8.11) in which the physical properties in the
Nusselt, Grashof, and Prandtl numbers are grouped together and
assumed constant. From equation (8.2a)
h = constan{kl-b(pg:
2CprJ!'f/3b-t
= constant x 0"1 3 b- 1 (8.17)
It will have been noted that b = 0·25 in laminar flow and 0·33 in
turbulent flow, so that the index of I is -0·25 in laminar flow and 0
in turbulent flow. The simplified expressions become
and
h = C(J0 ' 33 in turbulent flow (8.19)
(oy·ls
h = 0·00124(}0 "33 turbulent flow
Vertical surfaces { h ~ 000141 I laminar flow
I= height
h = 0·00131(}0 "33 turbulent flow
Horizontal surfaces
(oy·ls
I = length of side
ExAMPLE 8.1
An oil filled electric heating panel has the form of a thin vertical
rectangle, 2m long by 0·8 m high. It convects freely from both
surfaces. The surface temperature is 85°C and the surrounding air
temperature 20°C. Calculate the rate of heat transfer by natural
convection, and compare the result with that obtained from the
simplified formula for air.
Fluid properties at the average of surface and bulk air tempera-
tures, 53°C, are f3 = 1/326, Pr = 0·702, p = 1·087kg/m 3 , Jl. =
1·965 x 10- 5 Pas and k = 28·1 x 10- 6 kW/(m K).
Solution. The characteristic linear dimension is the panel height,
Q-8m. The product (Gr1Pr) must first be found.(}= 85- 20 = 65.
Gr1Pr = [(1/326) x 9·81 x 65 x 1·087 2 x 0·8 3 /(1·965 x 10- 5 ) 2 ]
X ()-702 = 2·15 X 109
Hence the flow is turbulent, and Nu 1 = 0·129(Gr1Pr)0 ' 33 may be used.
(Gr1Pr) 0 ' 33 = 1·29 x 103 Nu 1 = 0·166 x 103
and
0·166 X 103 X 28·1
1i = 0.8 x 106 = 5·85 x 10- 3 kW/(m 2 K)
tion which is from the top and bottom edges of the bars are negligible.
(Ans. 0·35 m.) (Queen Mary College, London).
REFERENCES
136
SEPARATED FLOW CONVECTION 137
(9.1)
Such a form is useful since it has been found generally that / 0 can
be related to the Reynolds number of flow. In the determination of
/ 0 from tip, the values of L, d, and Vm have to be defined in relation
to the particular geometry. The symbol / 0 is used to indicate that
it represents essentially a drag loss rather than a loss due to viscous
shear.
0·01 /
...
~0·005
.
::,... ""•
~
...:, ... 1" .2'!:
.
,_.,. :
0·002 q1GGER-FULL_
eX ·:11-PART
" -NONE
0·001 I I II II
0·005 0·01 0·02 0·05 0·1 0·2
fo
Re > 5,000
Fig. 9.1. Relations/rip between the /-factor and fo for flow in a compkx
system. By courtesy of H. Schenck, Jr., and The American Society of NaMI
Engineers, Inc.
ExAMPLE 9.1
Solution. From Fig. 9.1 the J-factor is o-01. The J-factor is given
by equation (7.27): J = StPrt, for mean conditions where St =
'Ft./pv.cP. At a mean temperature of 71°C, p for air is 0·998 kgfm 3, cP =
1·009 kJ/(kg K), Pr = 0·697. Hence the heat transfer coefficient is
given by
'ft. X (0·697)i = 0·0 1
0·998 X 30 X 1·009
:. 'Ft.= 0·386kW/(m2 K)
SEPARATED FLOW CONVECTION 139
LL \ J I ~ ,{ ........ ~ /'-
' '
I I \ _1_ _1
,
J I I
L \ I ! l '- '-./ '/
disc-and-doughnut baffles segmental baffles
Fig. 9.2. Types of sheU baffle; see also Fig. 13.2 which shows doughnut and
900 segmental baffles.
dis the tube diameter, Gmax = mass velocity = p x v where vis the
velocity through the smallest free-flow area between tubes, fluid
properties are evaluated at mean of wall and bulk fluid temperatures,
and 10 <Red < 40,000. For a more exact dependence on Reynolds
number, the numerical constant 0.33 may be replaced by 0·193 for
4000 < (dGmaxfJ.J.) < 40,000 and 0·0266 for 40,000 < (dGmaxfJ.J.) <
400,000. The corresponding indices of the Reynolds number are 0·618
and 0·805 (ref. 8). A much more detailed analysis for staggered and
in-line tube arrangements of different spacings was carried out by
Grimison. 9
On the shell side of shell and tube heat exchangers, two relation-
ships proposed by Donohue 10 may be used. The baffle arrangements
quoted are illustrated in Fig. 9.2.
For disc-and-doughnut baffles:
Nu4 =
dG
0·033d~' 6 ( ---:
)0·6 (Pr)t (:w)o·t4 (9.5)
For segmental baffles, the (0·033 d~' 6 ) in (9·5) is replaced by 0·25. Note
that G. = .J(GwG), where Gw = mass velocity through the baffle
window area, and G. = mass velocity based on flow area at the
diameter of the shell. Fluid properties are evaluated at the fluid
bulk temperature, with the exception of f.l.w which is at the tube wall
temperature. It is important to note that in using equation (9.5) all
terms are dimensionless groups except for (0·033 d~' 6 ). Here d. is an
equivalent diameter = 4(SrSv - nd 2/4)/(nd) where ST =tube trans-
SEPARATED FLOW CONVECTION 141
verse spacing, Sv = tube vertical spacing, d = tube diameter, and
d. is in mm.
The above equations give only very basic correlations of cross-
flow convection, for further information the reader is referred to
Kays and London. 11
EXAMPLE 9.2
In a shell and tube heat exchanger, the tubes are 25·4 mm diameter
and are spaced at 5(}8 mm centres both horizontally and vertically.
Water flows at 24·6 kg/sin the shell, and the baffie window area is
0·0125 m 2 and the net shell area is 0·05 m 2 • The water bulk tempera-
ture is 60°C and the tube wall temperature is 90°C. Calculate the
shell side heat transfer coefficient.
Solution. Property values of water are taken from Table AS. Thus,
= p x v, and at 60°C Jl. = 47·0 x 10- 5 , at 90°C Jl. = 31·9 x
Jl.
10- 5 Pas. Pr = 3·02, k = 651 x 10- 6 kW/(m K).
Equation 9.5 will be used. First calculate d., the equivalent
diameter.
(-
Jl. )0·14 = (-47 )0·14 = 1·056
Jl.w 31·9
3. Hydrogen passes through a staggered bank of 200 tubes, 1·8 m long, and
25·4 mm diameter. The mass velocity is 1·5 kgj(m 2 s). Calculate the rate of
heat transfer for a mean gas temperature of 373 K and a tube surface to gas
temperature difference of 50 K. Calculate also the heat transfer rate if air
at twice the mass velocity is substituted for hydrogen. (Ans. 499 kW, 70·2 kW.)
4. Carbon dioxide flows in the shell side of a shell and tube heat exchanger.
There are 36 tubes 15 mm diameter by 2m long. The shell area for flow is
0·025 m2 and the baffie window area is 0·0125 m 2 • The vertical and hori-
zontal spacing of the tubes is 22·5 mm between centres. The mass flow of
carbon dioxide is 0·6 kgfs at a mean temperature of 400 K. The mean tube
surface temperature is 300 K. Calculate the convective heat transfer coef-
ficient on the shell side of the tubes and the heat transfer rate. (Ans. 0·168
kW/(m 2 K), 56·8 kW.)
SEPARATED FLOW CONVECTION 143
REFERENCES
1. Schmidt, E. and Wenner, K. Forschung, Gebiete lngenieurw., Vol. 12,65
(1933).
2. Snyder, N. W. Chem. Eng. Progr., Symposium Series, Vol. 49, No.5, 11
(1953).
3. Schenck, H. Jnr. Heat Transfer Engineering, Longmans, Green and Co.
Ltd. (1960).
4. Schenck, H. Jnr. J. Amer. Soc. Naval Eng., Vol. 69, 767 (1957).
5. Douglas, M. J. M. and Churchill, S. W. Chem. Eng. Propr., Symposium
Series, Vol. 52, No. 18, 23 (1956).
6. Hsu, S. T. Engineering Heat Transfer, D. Van Nostrand Company, Inc.,
Princeton (1963).
7. Colburn, A. P. Trans. AIChE, Vol. 29, 174 (1933).
8. Holman, J. P. Heat Transfer, McGraw-Hill Book Company, Inc.,
New York (1981).
9. Grimison, E. D. Trans. ASME, Vol. 59, 583 (1937).
10. Donohue, D. A. Ind. Eng. Chem., Vol. 41, 2499 (1949).
11. Kays, W. M. and London, A. L. Compact Heat Exchangers, McGraw-
Hill Book Company, Inc., New York (1964).
10
Convection with phase change
144
CONVECTION WITH PHASE CHANGE 145
OVx
r =p.-
., oy
These shear stresses act over an area 1 x dx. Balancing the forces
146 ENGINEERING HEAT TRANSFER
1------y
'
b --j (dv,) = ..()
dy y=l
gives
d2vx pg
dy2 = ---;;
and on integration,
pgy2
v, = - 2il +ely+ C2
The boundary conditions are that v, = 0 at y = 0 and dvJdy = 0
at y = a, the thickness of the film. Hence
pga
c2 = 0, and - - + cl = 0
Jl.
The equation for v, is thus
v
X
=- pg(y2
Jl 2
- ya) (10.1)
CONVECTION WITH PHASE CHANGE 147
J: J:- p~g(~ -
film thickness b. Thus
2
m= pvxdy = yb) dy
or
dm kOw
(10.4)
dx hrgb
Equations (10.2) and (10.4) may be combined to give
p 2 gb 2 db = k(Jw
J1 dx hrgb
This result may be integrated between the top of the surface down
to x to give
or
(10.5)
148 ENGINEERING HEAT TRANSFER
= 8·08 kW/(m 2 K)
CONVECTION WITH PHASE CHANGE 149
hd = 0·770(~)*
hJC d
If (x/d) is 75, say, it follows that h4 = 2·26 h". Thus over twice the
fluid is condensed with the tubes arranged horizontally, h" being
the coefficient for the vertical tube. Condensation inside a tube is a
process of some interest since it occurs in refrigeration and heat pump
condensers, but it is a process of considerable complexity and is
conveniently described by empirical relationships. For low vapour
]*
velocities Chato4 gives:
~d p!
3 (p 1 - Pv) g (hrg + 0·68 Cpl(tsat - tw ))
Nud = 0·555 (10.9)
k1J.L1 (tsat- tw)
8, (wire- fluid)
Fig. 1 0.2. The boiling curve, after Farber and Scorah. 1
(10.10)
where
cP1 = specific heat of saturated liquid
hr1 = enthalpy of vapourisation
Pr1 = Prandtl number of saturated liquid
Jl1 = viscosity of saturated liquid
p1 = density of saturated liquid
Pv = density of saturated vapour
a = surface tension of liquid vapour interface
9 = heated surface saturation temperature difference
Q/A = heat flux per unit area
g = gravitational acceleration
C.r =experimental constant
Order ofmagi
nitude of co-
Quality Mechanisms efficients I
kW/(m2 Kl
Convection to
. . uperheated I· 7
I Gas
- vapour
I X= 10 Convection 17
Dry wall
only
1
I
! Fog or
I dispersed
I liquid
I
Annular Comectivc I 15 to
or Boiling
U1
260
climhing 1 depending on
I film velocity and
!
heat flux.
I Churn
or
0 Bulk boiling (but reducing
to 0· 3-6
0
I
slug
I· roth)
C>
..
oo.•o 1 with film
Saturated boiling)
I ~
or I •.•
. ·. nucleate
huhbl) boiling
l
I Suh~coolcd II 7=0 Suh-conled
nuclcatt: hodin'
licjllld I ConH~ction to
6
"-illl'r
Fig. 1 0.3. Flow and boiling regimes in a vertical heated tube. From data of
Firman, Gardener, and Clapp. 10 By courtesy of the Institution of
Mechanical Engineers
ExAMPLE 10.2
Using the Rohsenow equation, calculate the heat transfer coeffi-
cient for boiling when water boils at atmospheric pressure in
a copper pan with the copper surface at l20°C, and compare with
the convection coefficient for water flowing in a 40 mm diameter
tube at 1 m/s under the same conditions, using equation (7 .34).
Use Csr = 0·013, cP 1 = 4·216 kJ/(kg K), h~ = 2256·7 kJfkg, Pr1 =
1·74, J1.1 = 279 x 10- 6 Pas, p1 = 957 kgfm , Pv = 0·598 kg/m 3 , a=
0·0587 N/m. At a mean temperature of ll0°C, p1 = 950 kg/m3 ,
J1. 1 = 252 x 106 Pas, Pr = 1·56; k = 684 x 10- 6 kW/(mK) and at
120oC J1.1 = 230 X 10- 6 Pas.
154 ENGINEERING HEAT TRANSFER
Solution. The Rohsenow equation will give Q/ A from which h may
J(
be found. Thus:
4·216 X 20 = ()-Ol3 [ Q/A X 106 X 0·0587 \] 0 · 3 3
2256·7 X (1·74) 1"7 279 X 2256·7 9·81(957-0·598)1
Q/A = 358·0 kW/m 2
PROBLEMS
To solve Question I see also chapter 12.
1. An air heater consists of horizontal tubes 30 mm diameter and 23 mm
bore arranged in vertical banks of twenty. Air passes inside the tubes and
is heated from 32°C to 143°C by saturated steam at 180°C which passes
over the tubes. The mean air velocity is 23 m/s and the air flow 3·82 kg/s.
Calculate the number and length of tubes required. The heat transfer coef-
ficient for saturated steam to tube surface (h,.) can be found from
h = 0·725(k~p;ghr,)t kW/(m 2 K)
•• NdJl/it
where the suffix c denotes condensate properties evaluated at the saturation
temperature, g is the gravitational acceleration in m/s 2 , N is the number of
horizontal tubes in a vertical bank, d is the outside diameter in m, t is the
temperature difference between the saturated vapour and the tube surface
and may be assumed to be 11 a c. The other symbols have their usual
meaning. (Ans. 400 tubes, 2·34 m.) (Unirersity of Glasgow).
5. Outline the Nusselt theory of film condensation, indicating the steps which
lead to the following formula for the average surface heat transfer coefficient
hm during the condensation of a saturated vapour on a plane vertical surface:
N = hmL = 0·943(p1gL3hcg)*
Urn K J1Kf1T
Lis the height ofthe surface, g the acceleration due to gravity, hc 1 the enthalpy
of evaporation, 11 T the difference between the temperatures of the vapour
and the surface and p, Jl, and K are respectively the density, absolute viscosity,
and thermal conductivity of the condensate at the saturation temperature.
Saturated steam at 149°C is to be condensed in a cylinder of diameter
1·217 m and length 0·305 m, having its axis vertical. The curved wall is main-
tained at 10°C by external coolant and no condensation takes place on
the two horizontal surfaces. The steam is fed in through a pipe in the top
surface of the cylinder.
Determine the initial average surface heat transfer coefficient, and estimate
the time taken to fill the container with water which may be assumed to
remain at 149°C. (Ans. 4·85 kW/(m 2 K), 0·976 h.) (University of Cambridge).
REFERENCES
1. Nusselt, W. Z. d. Ver. deutsch. lng., Vol. 60, 541 (1916).
2. Nusselt, W. Z. d. Ver. deutsch. lng., Vol. 60, 569 (1916).
3. Hsu, S. T. Engineering Heat Transfer, D. Van Nostrand Company, Inc.,
Princeton (1963).
4. Chato, J. C. J. Am. Soc. Refrig. Air Cond. Eng., Feb., 52 (1962).
5. Akers, W. W., Deans, H. A. and Crosser, 0. K. Chem. Eng. Progr.,
Symposium Series, Vol. 55, No. 29, 171 (1959).
6. Carpenter, F. G. and Colburn, A.P. 'General Discussion on Heat
Transfer', I. Mech. E. London (19 51).
7. Farber, E. A. and Scorah, R. L. Trans. ASME., Vol. 70,369 (1948).
8. Rohsenow, W. M., Trans. ASME, Vol. 74,969 (1952).
9. Rohsenow, W. M. and Griffith, P. AIChE-ASME Heat Transfer
Symposium, Louisville, Ky (1955).
156 ENGINEERING HEAT TRANSFER
10. Firman, E. C., Gardner, G. C., and Clapp, R. M. /. Mech. E.
Symposium on Boiling Heat Transfer, Manchester, Review Paper 1
(1965).
11. Bagley, R. I. Mech. E, Symposium on Boiling Heat Transfer,
Manchester, Paper 13 (1965).
12. Jakob, M. Heat Transfer, Vol. 2, John Wiley, New York (1957).
11
Extended surfaces
Convection from a solid surface to a surrounding fluid is limited by
the area of that surface. It would seem reasonable, therefore, that if
the surface area could be extended, then a gain in total heat transfer
would be achieved. This is done by adding fins to the surface. Heat
transfer is then by conduction along the fin, and by convection from
the surface of the fin. It is likely that the convection coefficient of
the basic surface will be altered by the addition of fins, due to the
new flow pattern involved and the fact that the temperature of the
fin surface will not be uniform. Though the average surface tempera-
ture is reduced by the addition of fins, the total heat transfer is
increased. In the treatment that follows it is assumed that the
convection coefficient is known. The Nusselt numbers of finned
surfaces may be determined experimentally.
There are various types of fin, the most common being the
straight fin, the spine, and annular fin. The straight fin is rectangular
in shape and generally of uniform cross-section, and the spine is
simply a short thin rod protruding from the surface. Annular fins
are often found ifthe primary or basic surface is cylindrical. Examples
are to be found in heat exchangers and air-cooled petrol engines.
Extended surface nuclear fuel cans are shown in Fig. 11.1. These
are both straight and spiral in form.
Only the straight fin and spine will be considered here in detail.
Fins of non-uniform cross-se_ction and annular fins are more
complex mathematically, and the reader is referred elsewhere for
details. 1 •2 •3
157
158 ENGINEERING HEAT TRANSFER
transfer in a fin it is assumed that the thickness, or diameter of the
spine, is small compared with the length. Conduction along the fin
may then be assumed to be one-dimensional. The conduction and
convection heat transfers involved are shown in Fig. 11.3. Two
important dimensions of fins are their area of cross-section A, and
their perimeter P. In the straight fin it is convenient to assume that
a is small compared with L. Thus:
Fig. 11.1 Magnox fuel element, as used in the U.K. gas cooled reactor pro-
gramme. Note the external extended surface configuration. This illustration is
reproduced by courtesy of Nuclear Fuels, pic.
EXTENDED SURFACES 159
d 2t hP
dx2 - kA (t - t.) = 0
Surroundings
t,
Qo - - Qx+dx
.-----------+-~--------~
-- - X
x=O x + dx
Fig. 11.3. Heat transfer from an extended surface. This diagram implies heat
transfer from the fin to the surroundings. For the reverse case an inverted tem-
perature profile would be obtained.
d 2 (} hP (} _ O
dx 2 - kA -
(11.1)
where
m = (:~r (11.2)
(11.4)
0 cosh m(l - x)
-= (11.5)
00 coshml
Even though it was assumed that (dO/dx)(x=n = 0, the temperature
at the end of the fin is still above t 1 , and is given by
0 _ Oo
(11.6)
1 - coshml
This is obtained by putting x = lin equation ( 11.5).
The total heat transfer from the fin is obtained by considering
the conduction into the fin at the root. Thus:
Qo = - kA(dO)
dx x=O
- kA(d(}) = + hA01
dx x=l
- k(mC 1 eml - mC 2 e- ml) = + h01 ( 11.8)
The constants C 1 and C 2 may now be obtained by solving equations
(11.3) and (11.8). Substituting for C2 in (11.8) and eliminating 01 by
using (11.1):
-k[mC 1 eml- m(0 0 - C 1 )e-m1] = +h[C 1 em 1 + (0 0 - C 1)e-m1]
This then gives:
Oo[e-ml- (h/km)e-m 1]
CI - -~-'---;------o--=----,-
- (eml + e-ml) + (h/km)(eml _ e-ml)
and
C - Oo[eml + (h/km)eml]
2 - (eml + e ml) + (h/km)(eml _ e ml)
and on substituting back into equation (11.1) gives
(} em(l-x) + e-m(l-x) + (h/km)[em(l-x)- e-m(l-x)]
Oo = (eml + e-ml) + (h/km)(eml- e-ml)
(} - Oo
(11.10)
~-
cosh ml + (h/km) sinh ml
The heat transfer from the fin is obtained as before by considering
(dOjdx)x= 0 • Thus
Q0 = - kA(d(})
dx x=O
EXAMPLE 11.1
~'--
r-
40mm
-ll-
lmm
Clearly, the requirement for increasing heat transfer from the surface is
that h/km < 1, so that Q0 is greater than hA 8 0 which is the heat
transfer from the surface having no fins. The possibility of using fins for
insulation arises when h/km > I, and Q0 < h A 8 0 • It will be found that
to achieve a value of h/km > 1, the resulting fins are so thick that the
surface is virtually completely covered in insulation.
Q0 mkA60 tanh ml
'lr = Q~ = Plh6 0
This reduces to
EXTENDED SURFACES 165
Q0 tanhml
(11.13)
Q~ = ml
If the fin which has a significant end heat transfer is compared with
the ideal fin as defined by ( 11.12) then
tanh ml + h/km
'1r = (11.14)
ml + (hl/k) tanh ml
The fin effectiveness is a useful idea in relation to the next topic
to be considered, the overall heat transfer coefficients of surfaces
which have fins. In Chapter 3, overall coefficient were derived for
plane and cylindrical surfaces. Similar coefficients can be written
for surfaces, both plane and cylindrical, on which fins have been
added.
In the derivation of ( 11.13) it is seen that
('1rPl)h0 0 = mkA00 tanh ml
so '1r may be interpreted as the fraction of fin area which may be
regarded as being at 00 all over for purposes of calculating heat
transfer. A function 'Ire is now introduced which is the fraction of
area of a finned surface at 00 • If A 8 and All are the total fm surface
area and fin root area per unit area of primary or basic surface,
respectively, then the total area of surface at 00 is 1 - All + 'frA 5 .
As the total area is 1 - All + A 5 , the ratio of actual to ideal heat
transfer from a finned surface is
Fins are often added to only one surface to reduce the thermal
resistance on that side. However, Fig. 11.4 shows a plane surface
166 ENGINEERING HEAT TRANSFER
fluid a fluid b
t. tb
Heij!ht of
~I ,.. Height of
...,
a-fins b-fins
(b)
Fig. 11.4. (a) Plane finned surfaces &eparating two fluids. (b) Temperature pro-
files for (a). Heat transfer is to the a-fins from the fluid at t0 and from the b-fim
to the fluid at tb-
with fins both sides. Heat transfer to the a-fins per unit plane area is
given by:
EXTENDED SURFACES 167
(11.17)
In this analysis, the group (1 - AR + 'TeAs). means that AR, '1c and As
all refer to the a-fins, and similarly for the b-fins.
A similar result may be obtained for a tube finned internally and
externally, as shown in Fig. 11.5. Unit length of tube may be con-
sidered and the following three equations for heat transfer may be
written:
U _ 1/{ 1 +lnr 2 /r 1
L - 27tr 1(1 - Aa + '1rAJ.h. 27tk
UA - 1/{r1(1 - Aa
r2
+ 'lrAJ.h.
+ r21n rJrt
k
+ (1 - Aa ! '7rAJbhb}
(11.19)
This has units of kW/(m 2 K) where the area refers to the external tube
surface (primary) area. The corresponding result for the unfmned
surface is
(11.20)
EXAMPLE 11.2
are 6-Q and 1-Q kW/(m 2 K) and the thermal conductivity is 0·04
kW/(m K). A similar tube has 20 axial fins 14 mm high by 2 mm thick.
Find the overall coefficient in both cases, and in each case state which
thermal resistance is controlling.
Solution. In the first case, equation (11.18) simplified for no fins
becomes equation (3.17) with only one conduction term.
. U 1/{ 1 In 1·25 1 }
• . L = 27t X ()-01 X 6-() + 21t X ()-04 + 27t X ()-0125 X 1·0
= 1/(2·65 + 0·887 + 12·75)
= 1/16·3 = Q-0614 kW/(m K)
Although the inside and outside resistances are now similar, the
outside one is just still controlling.
170 ENGINEERING HEAT TRANSFER
tr
1--r--.---r----,----.,t,......ho--,---,.-,1_
I • I • I • I • b
~~~~~~~~T
2h 0
tl + t2 + ( bk tc -
x) (2 + bk
2h x)
to = 0
0
(11.22)
X=(T(I-1)+0·5•BTF(N)•TF)/(1·0+0·5•BTF(N)) (11.23)
EXTENDED SURFACES 171
Hence it is seen that the root node is at the surface temperature T(l)
and BTF(N) refers to the end node at I = N.
Fig. 11.7 shows the layout of a flat rectangular fin having a circular
or elliptical root, which is approximated to the rectangular grid.
Such a fin is usually symmetrical, so only one quarter need be con-
sidered. The general BASIC program in Chapter 4 is suitable for a
outline of root
,...----,
: I
1 12 I
,----:-- -~-l--.
I I I :
I 13 t1 1
I I to I Y
L ___ _j _ _ _ L ___ .J
1
I I f
___I I_
I
X
element thickness b, : 14
convecting to tr on both L ___ ..J
faces
Fig. 11. 7. F1at rectangular fin with circular root.
kyb ( - to
t1 -- +tJ--- to) + k b (t2---to+t4-- -to)
X
X X y y
+ 2hxy(tc - t0 ) =0 (11.24)
Other equations are needed for side and corner boundary nodes,
obtained from similar energy balances.
172 ENGINEERING HEAT TRANSFER
PROBLEMS
1. The diagram shows the cross-section of a nuclear reactor fuel element
consisting of a uranium fuel rod 28 mm diameter contained in a magnox
can which has longitudinal finning on its external surface. The fuel rod/can
interface temperature is 430°C and the heat release rate is 65·6 kW per m
length. Calculate the maximum temperature within the fuel rod, the tern-
fuel rod
urrounding
graphite
longitudinal fins
magnox can
perature drop across the 2·5 mm wall of magnox, and the total surface
area of longitudinal tinning per unit length, given that the fin root is 40 per
cent of the external can surface area, fin efficiency is 90 per cent (based on
fins + splitters which also act as fins), mean coolant gas temperature 280°C,
surface convection coefficient 0·8 kW/(m 2 K), k for magnox 0·15 kW/(m K~
k for uranium fuel 0·0325kW/(mK). (Ans.: 591°C, 11·5°C, 0·589m 2/m.)
(The City University).
174 ENGINEERING HEAT TRANSFER
3. A bar simulating a gas turbine blade, cooled at the root, is 10·2 em long
and has a cross-sectional area (A) 1·93 cm 2 , and a perimeter (p) of 7·6 em.
Gas at 815°C streams across it, and one end is cooled to 483oC. The mean
heat transfer coefficient for the gas flow conditions can be assumed constant
over the surface at ()-284 kW/(m 2 K), and the thermal conductivity of the
material of the bar (k) is 26 x 10- 3 kW/(m K~ Show that
9" cosh mL(1 - x/L)
9, = coshmL
where
9" = t1 - t", 9, = t1 - t, m = / :~
L is the length of the bar and x the distance from the cooled end, t 1 being the
gas temperature and t_, and t, the appropriate blade temperatures. Hence
find the heat passing to the cooled end ofthe bar in kW. Neglect the heat lost
from the uncooled end of the bar. (Ans.: 0·11 kW.) (Queen Mary College,
London).
4. Heat flows from a body A along a wire of diameter d and length I, the other
end of the wire being connected to a body B. The conductivity of the wire is
k and the surface coefficient of heat transfer ia 2 kd.
The temperature of the body A is maintained at 8A above the temperature
of the environment, and the temperature of the body B is 88 above the tem-
perature of the environment. Derive an expression for the temperature of
the wire at x from A, and deduce the particular values of 88 for which
(a) heat flow into B is one-half of the heat flow from A,
(b) heat flow into B is zero. (University of Oxford).
REFERENCES
1. Jakob, M. Heat Transfer, Vol. 1, John Wiley and Sons, Inc., New York
(1949).
2. Eckert, E. R. G., and Drake, R. M. Analysis of Heat and Mass Transfer,
McGraw-Hill Book Company, Inc., New York (1972).
3. Chapman, A. J. Heat Transfer, 3rd ed., The Macmillan Company, New
York ( 1974).
12
Heat exchangers
Much of the basic conduction and convection theory finds its
greatest application in the heat exchanger. Whenever it is necessary
to transfer energy from one fluid to another in large quantities,
some form of heat exchanger is used. The most common form of
heat exchanger is that in which two fluid streams pass through in
steady flow, and heat transfer takes place through a separating wall.
Mechanisms involved are therefore convection to or from the solid
surface and conduction through the wall. The wall may be corrugated
or finned to increase turbulence and the heat transfer area.
The thermal capacity of a heat exchanger is usually kept small,
and is of significance only in transient conditions. However, a
regenerative type of heat exchanger does have a large thermal
capacity matrix through which the hot and cold fluids _pass alternately.
By this means energy is transferred indirectly between the fluids.
This chapter is concerned only with non-regenerative heat ex-
changers in which the fluids are separated. Other types of heat
exchanger exist in which the fluids mix. These include cooling
towers and jet condensers, for example. The basic principles will
be considered in relation to the simplest types only.
~~~;1 I .1 !J -~~~~j
(c) Tube and shell exchanger in parallel flow.
Tube fluid may be reversed to give
oount" How. I~
IIU
(I f--
tl ~
(d) Two-pass tube and shell exchanger.
Flow is partly counter, partly parallel.
Fig. 12.1. Basic in-line heat exchangers
of the cold fluid. Suffixes 1 and 2 are used for inlet and outlet of
individual streams, and ei is the temperature difference between
fluids at the inlet end and eo the difference at the outlet end of the
exchanger. An important term in heat exchanger theory is the
capacity ratio C. It is a ratio of the products of mass flow rate and
'I-0 0
(b) Converging at inlet end (c) Diverging at inlet end
Counter flow
thl th I th2 _l
r ·rI
I Oo
0; I
0; t c2
.l j_
tel t, I
(d) Boiling (e) Condensing
(12.1)
180 ENGINEERING HEAT TRANSFER
(12.2)
(12.3)
and
E = thl - thz when ( 12.4)
thl -tel
These definitions may be used in either counter or parallel flow,
but the value of E will be lower in parallel flow.
Temperature distributions with a change of phase are also
shown in Fig. 12.3. These will occur in boiling, Fig. 12.3d, and
condensing, Fig. 12.3e. Only the phase change takes place in the
exchanger, so the temperature of the boiling or condensing fluid
does not change. The temperature distributions are the same for
both parallel and counter flow. The capacity ratio C becomes 0
HEAT EXCHANGERS 181
for both boiling and condensing since the larger thermal capacity
is in each case infinite. This follows, since by definition, cP = dh/dt
= oo when dt = 0. Equations (12.3) and (12.4) may be used in
condensing and boiling, respectively.
The other limit of capacity ratio is C = 1 and occurs when the
thermal capacities of the two streams are equal. This is not illustrated,
but it results in the temperature distributions being parallel straight
lines in the case of counter flow, () being a constant over the whole
heat exchange area.
(12.6)
( 12.7)
()m is a mean temperature difference between the fluids, and U A
and U L are mean coefficients, in kW/(m 2 K) and kW/(mK) or equiv-
alent units, applicable over the entire area A or length L of the ex-
changer. It is general practice to work in terms of the external
surface area of the tubes in heat exchanger design, and the overall
coefficient UA in terms of this area is given by equations (11.19)
for finned surfaces and ( 11.20) for plain surfaces.
12.2.1 Counter and Parallel Flow
If the mass flow rates and inlet and outlet temperatures are known, the
heat transfer Q will be known, but further details of the exchanger
cannot be specified until Om is known. Bm can be derived as follows:
182 ENGINEERING HEAT TRANSFER
(12.9)
(12.10)
Parallel flow
l A or L
0 dA or dL
Counter flow
Fig. 12.4. For the determination of logarithmic mean temperature difference
dO = d(th - tc)
= dth- dtc (12.11)
HEAT EXCHANGERS 183
If d8 is expressed in terms of dQ from (12.9) and (12.10),
(12.12)
and
(12.13)
since dth is negative in both cases, and dtc is positive for parallel
flow and negative for counter flow; It follows that dO/dQ has
constant but different values for parallel and counter flow, and
therefore
(12.14)
with + for parallel flow and - for counter flow. Further, dQ from
equation (12.8) may be substituted in (12.12) and (12.13) to give
(12.15)
EXAMPLE 12.1
0·2 kg/s of an alcohol is to be cooled from 75 to 35°C in a counter
flow heat exchanger. Cooling water enters the exchanger at l2°C
and at the rate of 0·16 kg/s. The convection coefficient between the
alcohol and the tube wall is 0·34 kW/(m 2K), and between the
tube wall and the water, 0·225 kW/(m 2K). The tubes may be assumed
thin. cP for the alcohol is 2· 52 kJ/(kg K) and for water is 4·187
kJ/(kgK).
Calculate the capacity ratio, the effectiveness, and the area of
the heat exchanger surface.
The heat exchange area may be found from equation (12.5). To find
uA = 0·1355 kW/(t;n 2 K)
Equation (13.5) gives 20·15 = UAAOm = 0·1355 x A x 28
A = 5·31 m 2
'"- 1·0
~ 0·9
0 "' 1'\
.l! 0·8
.B=o-7
1\
~
\81 "!i i\~ ~
~ 0·6 II
~
0
u 0·5 6 . \ ll
0 0·1 U· ij·j ·4 ur:> U·6 0·7 0·8 0·9 1·0
t,2- t,l
t.l - ttl
Z = tsl - ts2 = (mc.,)tube
t,2 - ttl (mcp)shell
Fig. 12.5. Logarithmic temperature difference correction factor for cross flow,
one fluid mixed, one fluid unmixed. From R. A. Bowman, A. E. Mueller, and
W. M. Nagle, Trans. ASME, Vol. 62, p. 283 (1940). By permission of the American
Society of Mechanical Engineers.
and,
(12.19)
= -dQ(l- C)
dO
0
Integrating:
= -NTU(l- C)
188 ENGINEERING HEAT TRANSFER
= 1-£ =e-NTU(l-C)
1- CE
from the right-hand side, above. This final result is now rearranged
to give
1 _ e-NTU(l-C)
E = 1 - C e NTU(l c> (12.21)
If mecpe had been assumed the smaller quantity, the same equation
would have been obtained, where E, NTU, and C would have
then been defined by the alternative expressions.
A relationship exists, then, between E, NTU, and C given by
equation (12.21). Using this result it is possible to determine outlet
temperatures tel and th 1 , and Q, the overall heat transfer for a given
design, without using a trial and error solution.
EXAMPLE 12.2
Determine the effectiveness and fluid outlet temperature of an oil
cooler handling 0·5 kg/s of oil at an inlet temperature of 130°C. The
mean specific heat is 2·22 k1/(kg K). ()-3 kg/s of water entering at 15°C
passes in counter flow at a rate of 0· 3 kg/s. The heat transfer surface
area is 2·4 m1 and the overall heat transfer coefficient is known to be
1·53 kW/(m1 K)
Solution. The thermal capacities are: oil, 0·5 x 2·22 = 1-11 kJ/(s K),
water, 0·3 x 4·182 = 1·255 kJ/(s K)
c = 1-11/1·255 = Q-885
HEAT EXCHANGERS 189
and,
1·53 X 2·4
NTU = = 3·31
1·11
Then,
1_ e-3·31(1-o 885)
1 - e- 0 . 38 0·316
1 - 0·885 e 0 ' 38 = 0·395 = 0· 8
1300 - th2
130° - 15°
(oil outlet)
By enthalpy balance
1·11 X (130- 38)
(tc2 - tel) = = 81·5 K
1· 255
(water outlet)
Thus, the effectiveness is the same for both counter and parallel
flow.
The other limiting value is C = 1 for equal thermal capacities
and, in this case, for parallel flow equation (12.22) gives
1 _ e-2NTU
E=--2-- (12.24)
and also
Also
E may be written as
100
unmae _1
l 1.-- 1--
/ ".... ~ - --
- - - ~~ -
0·25
80
~ ;.::: ~ --
- - --
/ 0·75
~60 ~ ~ ,.,r- 1-33
I
.:: 40 ~
-
~
~ f..- I'Cmc.J,I,..d
(mc,),...laod
=1
u L
~
l.l.l 20 I
I
2 3 4 5
Number of transfer units, NTU = U AA/(mcp)mln
=f=:;~-
Fig. 12.6. Effectiveness vs. NTU for a cross-flow exchanger, one fluid mixed,
one fluid unmixed. From Compact Heat Exchangers, by W. M. Kays and A. L.
London, McGraw-Hill Book Company, Inc., New York (1958). Used by
permission of McGraw-Hill Book Company.
(12.26)
Transient: c)
(m• p 1 (t 3,0 - t 1,0 ) -- H 1 [tM,o - t1,0 + t3,0J (12.27)
2
- (mCv)t(t3,1 - t3,o)
(12.28)
Mass within
--
f-....
I
I
I
tl J_...-
---- lc·"-.
(mcph
.....
t4
tI
~
1
t3
M
(2
t _ [ 2NTUt
3' 0 - 2 + NTU1
Jt M,o
+ [22 +~NTU
NTU
1]
ft,o
(12.29)
1
J 2 J
t4,o =
[22+NTU
2
NTUz tM,o +
[22 ~NTU
+ NTUz t2,o
(12.30)
In order to calculate t 3 and t 4 after the next time step, an equation for
a new tM must be obtained. Thus:
(12.31)
194 ENGINEERING HEAT TRANSFER
and hence
(12.32)
where (mCp)M is the product of mass and specific heat of the metal
element. For steady state, an enthalpy balance holds:
(12.33)
From this equation, and equation (12.26), the following relations for
t 3 and t 4 are obtained:
t =[ 2+NTU+CxNTU
3
2NTU ]t 2
+ [1 - 2NTU
2+NTU+CxNTU
]t 1
(12.34)
[ 2xCxNTU
t 4 = 2+NTU+CxNTU t 1 +
J c-
f, 2xCxNTU
2+NTU+CxNTU tz
J
(12.35)
where C = (m Cp) 1 /(m Cp)z and NTU = H/(m Cp) 1
The above equations may be used in the analysis of transient and
steady state behaviour of cross-flow heat exchangers. As shown in Fig.
12.7 outlet temperatures from one node become the inlet temperatures
of the next node, so that the field of nodes may be solved in simple
computer routines. Examples are given in Ref. 4. The equations to be
used in any analysis are (12.29), (12.30), (12.32), (12.34), and (12.35).
In each case the last term of the equation must remain positive for the
solution to be stable. This sets a limit on the values of the variables
occurring.
t2 =[ 2NTU1
2 + NTU1 + C1 x NTU1
JrA+[l- 2NTU1
2 +NTU1 + C1 x NTU1
Jr 1
(12.36)
tB = [
2xC 1 xNTU1
2 + NTU1 + C1 x NTU1
~
t1 + 1 -
2xC1 xNTU1J t
2 +NTU 1 + C1 xNTU1 A
J
(12.37)
!z
./
/ /
/.
,::· /
/
j', /
I ·, -y· / /
/
/
/
i i /
/
/
I /
(mcp)l 1.::" /
gas flow '- , '- . I / / ...- / t~
'- !..- (mcp)M
matrix flow
For the element shown, the gas stream entering at t 1 is in heat exchange
with a moving solid matrix stream entering at t A• and having a thermal
capacity rate of (m Cp)M. In the two equations above:
196 ENGINEERING HEAT TRANSFER
Later in the cycle, the matrix is in heat exchange with a second gas
stream, and similar equations may be written, including the terms:
Usually
NTU2 = NTU 1 and C2 = C 1
The elements of the form considered in Fig. 12.8 are pieced together
as shown diagrammatically in Fig. 12.9. Imagine Fig. 12.9 as a fixed
grid or series of control volumes with the fluids and matrix flowing
through in the direction of the arrows. The actual thermal wheel is
shown schematically in Fig. 12.1 0. The matrix passes through sealing
zones between the two gas streams, and it is assumed that there is no
change in the matrix temperature across this zone. However, as the
matrix passes out of the second gas stream, the matrix temperature at
all nodes must equal those values at inlet to the first gas stream, for
steady state conditions. This means that for a correct solution an itera-
tive procedure is necessary, since the matrix inlet temperature tA, tK,
tL and tM must be assumed or guessed at the start of the analysis.
Values obtained at outlet in the first iteration are inserted back in the
inlet and the process repeated until a close convergence is obtained.
Heat transfer between the gas streams is obtained by summing up the
enthalpy loss or gain of each fluid stream, and average gas stream outlet
temperatures may be found to give the regenerator effectiveness. For
steady state conditions the cooling of the matrix by the cold stream
must always equal the heating by the hot stream, and hence the capacity
ratio of the two gas streams is usually (though not necessarily) unity. In
this case it follows that C 1 = C2 and NTU 1 = NTU2 •
197
!
HEAT EXCHANGERS
,..,, ), t (mcph
I 14 I 15 I 16 I 17 I 18 I 18 118 I 18
110 Ill I 12 I 13
I6 17 I8 t9
12 t3 14 t5
IB tc lo IF IG IH lj lA
tc
EXAMPLE 12.3
400 TIJX,J%l~f5•AIJ4,JX-1J+F6•TIJX-1,1Xl
404 UEX T IX
405 IHXT J%
410 AA(1,1l~F3•TZ+F4•AI10,10l
420 TT(1,1l~F7•AI10,10l+F8•T2
430 FOR n ~z TO 10
440 AAIJ%,1l~F3•TTIJ%-1,1)+f4•AI10-J4+1,10l
4SO TTIJ~,1l~F7•AI1u-J~+1,10l+f3•TTIJ~-1,1l
4S5 llEXT J4
460 fOR 14•2 TO 10
470 AA(1,J:l•F3•T2+F4•AA(1,1~-1)
480 TT(1,J~l=F7•AA(1,J%-1)+f~•TZ
485 NEXT I~
490 FOR J%~7. TO 10
500 FOR IX•Z TO 10
510 AAIJ%,Jil~FJ•TTIJ%-1,JXl+F4•AAIJX,J~-1)
520 TTIJ%,1~)•f7•AAIJ%,1X-1l+FB•TTIJ%-1,1%l
524 IIEXT n
525 :lEXT J:':
530 X•O
540 FOR J:•1,TO 10
5SO X•X+AA(JX,10l;NEXT J%
560 Z•AOS(X-Y)
570 lf(l<0.001lGO TO 610
580 Y•X
590 ITERX•ITER%+1
600 GO TO 290
610 TOUTZ•O.O
620 FOR 1%•1 TO 10
630 TUUTZ•TOUTl+TT(10,1%)/10.0
635 NEXT IX
640 Eff•(TOUT2-T2)/(T1-T2l
650 LPRINT " Gas Strea• Outlet TeMperatures:"
660 LPRINT USING •nNNN.MH";T(1Q,J%) FOR IX•1 TO 10
670 LP~INT " Air Strea• Outlet Te•perdtorrs:"
680 LPRINT U~ING •nnnn.gn•;Tr<10,JXl FOR 1%•1 TO 10
690 LPRINT '' Effectiveness ~ ";LPRINT USIUG"I.#II";EFF;
700 LPRI~T '' for '':LP~IUT U~lNG"nU.N#I":TII:LPRlUT'' Trdnsfer Units"
710 LPRINT ''lluYber of Iterations to Conver~e • ";
720 LPRIUT USIJ4ti ''U~M":tTEit
730 GO TO 7?0
740 LPRINT " Check Data - C•lculMtions Unstable:"
750 LPI«INT " f1 a ";LPNIIIT IJSlriG "##.d#";Fl;
760 LPRIUT •• FJ • ":LPHINT US!UG "nR.MH":F3:
770 LPR!UT " F5 • ":LPRINT USliJG "##,#N";rs;
780 LP:liur •• F7 =- ":LPRtNT ustuG "nn.M"":F7
790 STOP
'i-;lil
1.0 C=
0.1
;:;s 0.5
6-t····~-l·~
0.8 1.0
.L.,~~~~
E
0.6
:~·.·.·~~~ f1T
0.4 !), '
'
L
t
1-
j
)..
+-
+-
,='·
,_,_
tr +·
0.2
i --1=
0 5 10 15 20
NTU
Fig. Example 12.3. Thermal Wheel characteristics.
200 ENGINEERING HEAT TRANSFER
12.5 Plate Heat Exchangers
The plate type of heat exchanger is basically of the in-line type, but the
construction is very different from the conventional shell and tube
concept. A plate heat exchanger consists of a frame in which a
number of heat-transfer plates are supported and clamped between
a header and a follower. Each plate has four ports and the edges of
the plates and ports are sealed by gaskets so that hot and cold
fluids flow in alternate passages formed between the plates. This
means the fluids flow in very thin streams having a high heat-transfer
Fig. 12.11. A typical flow diagram of a plate heat exchanger showing a two·
pass arrangement (diagram by courtesy of the A.P. V. Company Ltd).
Fig. 12.12. A Paraflow type Rl45 plate heat exchanger, capable of accepting
up to 955 m 3 per hour at 1 (). 7 bar, and up to 130°C; plate size is 2122 X 849 mm
(photograph courtesy of the A.P. V. Company Ltd).
where PN is the plate number, U..1 a. For mhcph being the minimum
capacity rate, or for equal rates as defined previously, it is seen from
equations (12.19), (12.5) , and (12.7) that (th 1 - th 2 )/8m =NTU, the
number of transfer units, and hence
me
n=~ x NTU (12.38)
PN
The performance of a particular plate design can be expressed
graphically in terms of the plate number, the NTU value, and the
pressure drop plotted against the plate rate, or the mass flow rate
across a plate, see Fig. 12.13. Separate curves would exist for different
202 ENGINEERING HEAT TRANSFER
capacity ratios, and from such information for various plate designs,
the required unit for a particular duty can be selected. Certain
correction factors have to be introduced, on account of concurrency
and other effects which depend upon the particular plate arrange-
ment, and on account of uneven distribution along the plate pack
due to pressure losses along the ports. For exactness liquid proper-
ties have also to be considered, and separate relationships would
apply to laminar and transitional flow.
pressure
loss
plate rate
Fig. 12.13. Plate number, pressure loss and NTU characteristics of a plate
heat exchanger for turbulent flow.
temperature at any given time, and the theory of Section 5.1 applies.
Figure 12.14 shows the temperature time curve for the heat exchange
process, and from equation (5.1), the temperature-time curve is given
by:
8 2 /8 1 = exp (- t/T) (12.39)
where t is the heating time between temperature differences of 8 1 and
8 2 and T is the time constant mCp/ hA, where m is the mass of heated
and stirred fluid, CP is the specific heat, his the convection coefficient
between the pipe and fluid, and A is the area of the heating surface
(e.g., A= rrdL for a pipe) in contact with the heated fluid.
Batch temperature
Time ~
PROBLEMS
1. A tubular heater of the counter flow type is used to heat 1·26 kg/s of
fuel oil of specific heat 3·14 kJ/(kg K) from 10° to 26·7°C. Heat is supplied by
means of 1·51 kgfs of water which enters the heater at 82°C.
(a) Derive an equation relating the temperatures of oil and water at any
section of the heater.
(b) Determine the necessary surface if the rate of heat transfer is 1-135
kW/(m 2 K). (Ans.: 1·013 m 2) (University College, London).
204 ENGINEERING HEAT TRANSFER
2. In a test on a steam condenser the rate of flow of cooling water was varied
whilst the condensation temperature was maintained constant. The following
results were obtained:
Overall heat transfer coefficient K, kW/(m 2 K) 2·7 2·98 3·39 3·59
Water velocity V, m/s 0·986 1·27 1·83 2-16
Assuming the surface coefficient on the water side to be proportional to
yo·&, determine from an appropriate graph, the mean value of the steam side
surface coefficient. The thickness of the metal wall is 0·122 em and thermal
conductivity of tube material 0·111 kW/(m K). (Ans.: 6·04 kW/(m 2 K.)
(University of Manchester).
water 30"C
water 30°C
-----<---tll!llll!lllll!llllll!llllll water 80"C
fluid
90oC
(b)
chamber of a gas turbine cycle. The pressure ratio of the cycle is 5: 1 and the
heating fluid is the exhaust from the turbine which expands the gas from
65o•c with an isentropic efficiency of 82 per cent. If the air conditions initially
are 1·013 bar and 21 •c and the isentropic efficiency of the compressor
is 80 per cent, calculate the area of heat exchanger for a total fluid mass flow
of 22·7 kgjs.
Assume a logarithmic mean temperature difference and constant specific
heat of 1·0 for the air and 1·09 kJ/(kg K) for the products. y = 1·4 for air
and products. (Ans.: 424m 2 .) (University of Manchester).
10. Define the terms 'effectiveness' and 'number of transfer units' as applied
to heat exchangers stating any assumptions involved. Obtain a relationship
between effectiveness and number of transfer units for a counter-current heat
exchanger and plot this relationship when the ratio of the stream heat
capacities is 0·5.
20·15 kgjs of an oil fraction at a temperature of 121 oc is to be cooled in a
simple counter-current heat exchanger using 5·04 kg/s of water initially at
1o•c. The exchanger contains 200 tubes each 4·87 m long and 1·97 em outside
diameter; the resulting heat transfer coefficient referred to the outside tube
area is 0·34 kW/(m 2 K). If the specific heat of the oil is 2·094 kJ/(kg K)
calculate the exit temperature of the oil. (Ans.: 90·8°C.) (University of Leeds).
REFERENCES
I. Smith, D. M. Engineering, Vol. 138, 479, 606 (1934).
2. Bowman, R. A., Mueller, A. C., and Nagle, W. M. Trans. ASME, Vol. 62,
283 (1940).
3. Kays. W. M., and London, A. L. Compact Heat Exchangers, McGraw-Hill
Book Company, Inc .. New York (1964).
HEAT EXCHANGERS 207
4. Simonson, J. R. 'Transient and steady state analysis of cross flow
heat exchangers by programs in Fortran', Trans. Inst. Chern. Engrs.,
Vol. 55,53-58 (1977).
13
The laws of black- and grey-body
radiation
208
THE LAWS OF BLACK- AND GREY-BODY RADIATION 209
qb = s= qb .. d.A. = (jT
4 (13.5)
The variation of qb;. with wavelength was established by Planck 2 in
his quantum theory of electromagnetic radiation, thus
q -
c 1A. -5
b.l. - exp (C 2 jA.T)- I
where ). = wavelength, J.lm, T = absolute temperature, C 1 = 3·743
x 105 kWJ.1 4 fm 2 , C 2 = 1·439 x 104 J.1K. The form of the varia-
tion of qbx is shown in Fig. 13.2, and it is seen that there is a peak
value of qb.l. which occurs at a wave length which is related to the
absolute temperature by Wien's displacement law:
).max T = 2897·6!-1 K
Real materials that are not black will have monochromatic
emittances that are different from qbA.• and hence it is useful to
define a monochromatic emissivity e;. by the equation
q;, = e;.qb'-
or
(13.6)
212 ENGINEERING HEAT TRANSFER
The black and non-black emittances which give e;. are measured at
the same temperature. In general, e;. is a function of wave-length,
temperature and direction. Real surfaces often exhibit directional
variation in emissive power, thus non-electrically conducting
materials emit more in the normal direction whereas for conducting
materials often the reverse is true. For practical calculations, quoted
emissivities are total hemispherical values. Most real materials
exhibit some variation in e;. with wave-length. These are known as
selective emitters. However, there is a second type of ideal surface,
known as a grey surface, where the emissivity is constant with
wave-length. Some real materials approximate closely to this ideal,
but the concept reduces calculations to the extent that it is worth-
while to accept the error introduced in exchange for the simplifica-
THE LAWS OF BLACK- AND GREY-BODY RADIATION 213
Fig. 13.3. Comparison of the emission of black, grey, and selective emitting
surfaces; ex= q-,.,/qbx.
But
from (13.8)
Therefore
e=oc {13.9)
13 .S Intensity of Radiation
The radiation from a unit area of black body is qb = uT4 • For
diffuse radiation from a small flat area of black surface dA, the
entire emittance Qb must pass through a hemispherical surface
surrounding the emitting area. It is necessary to consider the
distribution of radiant energy per unit area over the spherical
surface, before calculations can be made of radiation exchanges.
The intensity of black-body radiation, /, is the radiation emitted
per unit time and unit solid angle subtended at the source, and per
unit area of emitting surface normal to the mean direction in space,
and may be expressed as
(13.10)
(13.11)
and from Fig. 13.6 it is seen that dA 2 = rdcp (r sin cfJ d8) =
r 2 sin cjJ dc/J dO. Hence
then
f"'="'2s6=2"
Qb = I dA 1 tJ>=O B=o sin cos dcjJ dO cjJ c/>
normal/
apart, and the line joining their centres makes angles </l 1 and </l 2
to their normals. Each element of area subtends a solid angle at the
centre of the other; these are dw 1 subtended at dA 1 by dA 2 , and
dw 2 subtended at dA 2 by dA 1 . The solid angles are given by:
d w 2 = dA 1 cos </l1
an d
r2
From equation ~13.11) the radiant energy emitted by dA 1 that
impinges on dA 2 is given by:
(13.13)
Since both surfaces are black this energy is absorbed by dA 2 • A
similar quantity of energy is also radiated by dA 2 and absorbed by
dA 1 expressed as
(13.14)
218 ENGINEERING HEAT TRANSFER
fi
total energy radiated by A 1 that falls on A 2 is given by
f
Hence the fraction of energy radiated by A 1 that falls on A 2 is
ff
that falls on A 1 is given by
l
and the total energy radiated by A 2 is A 2 aTi, so that
= Fz-1 (13.18)
and
(13.19)
From equations (13.16) and (13.18) it is seen that F 1 _ 2 andF2 _ 1 are
simply related:
(13.20)
The net radiation exchange from equations (13 .17) and ( 13 .19) can
be expressed in terms of either configuration factor, thus
Qb(12) = F1-2A1a(Ti- Ti)
(13.21)
It is necessary to know or to be able to calculate configuration
factors before black-body radiation exchanges can be determined.
Only a few results will be considered here, and the reader is referred
elsewhere for further information on this subject. 1 •5 · 6
13.6.1 Examples of the Black-Body Geometric Configuration Factor
(i) Cases where F 1 _ 2 = l. The simplest case is when surface A 1 is
entirely convex and is completely enclosed by A 2 . Then F 1 _ 2 must
be 1, since all the energy radiated by A 1 must fall on A 2 . It follows
also that F2 _ 1 is Ad A 2 . In this case. the net black-body radiation
exchange is
(13.22)
Another simple example is when surfaces A 1 and A 2 are parallel
and large, and radiation occurs across the gap between them, so
that in this case A 1 = A 2 and all radiation emitted by one falls on
the other if edge effects are neglected. Hence,
Fl-2 = F2-t = 1
Concentric surfaces may be included if the gap between them is
small so that little error is introduced by the small difference
between the area of A 1 and A 2 . The net radiation exchange is again
given by equation (13.22).
220 ENGINEERING HEAT TRANSFER
y;f
·----~---
dAl . 2
i
f---1- r-dl
L--L I L ___j
Fig. 13.8. The thermocouple configuration factor.
f
Applying equation (13 .16) gives
Qb<J- 2 > = dA 1 cos </1 2 ~nR dl
Qb(ll A1 A, nr
=
1[ U
4 (R2 + F)t
J+L
-L = (R2
L
+ U)t (13.23)
THE LAWS OF BLACK- AND GREY-BODY RADIATION 221
ExAMPLE 13.1
. 0·3
I
0·4
0·6
~ 0·2 /' /' 1·0
/'
I·~
~
0· )/' :;;;-- ...... 4·0
6·0 \0·
0 20·0
0·1 0·2 0·3 0·4 0·6 0·81·0 2·0 3·0 4·0 6·0" 8·0 10
R 2 = WjD
Fig. 13.10. Configuration factor for perpendicular rectangles with a common
edge. (From A. J. Chapman, Heat Transfer, The MacMillan Company, New
York (1974). By permission of the publisher.)
+ A1F1-3 a (Ti- TD
(13.24)
For any number of surfaces, the configuration factor from one to the
remainder is the sum of the component factors. This result is useful
since unknown factors may be deduced from known standard cases.
I
I
I
I
I
I I
/ ) - - - - - - - - - - - _L __ _
/'
/
/
/
//
............
EXAMPLE 13.2
Qb(12)=1;
THE LAWS OF BLACK-AND GREY-BODY RADIATION 225
J = eQb + pG
1 _ G = 1 _ J - eQb
p
and since p + e = 1 for opaque surfaces, this reduces to
e
J- G = -(Qb- J)
p
If two surfaces only are involved, and these form an enclosure, this
is also the net energy exchange between them, Q( 1 2 J, and the equa-
tion may be compared with Ohm's Jaw so that QbiA, which is uT4 ,
and J IA are potentials and PI Ae is the resistance. The correspond-
ing circuit element for either surface is shown in Fig. 13.13.
R = _!!__
Ae
o----c=J--o
V= QJA V= JfA
Fig. 13.13.
Jl J2)
Q<l21 = ( A! - A2 A!Fl-2
R=-l-
A•F•-2
o---c=:::J--0
V=J 1/A 1 V=J 2/A 2
Fig. 13.14.
~+ 1 +~
A1e1 A1F 1_ 2 A 2e2
1
R=--
A•Ft-2
1
= A2F2-t
Fig. 13.15. Complete circuit for radiation exchange between two grey surfaces
forming an enclosure.
1 (1-e 1 ) 1 A 1 (1-e 2 )
~-2 = -e-1- + F1-2 + A2 _e_2_
1 1_ 2 2
This result may be used for any two surfaces of area A 1 and A 2 ,
provided they form an enclosure, for which the configuration factor
is F 1 -2- Equation (13.26) simplifies for the special case of infinite
parallel or concentric grey planes for which F 1 _ 2 = 1 and A 1 = A 2 •
Then,
~-2 = - - - - - (13.27)
1 1
-+- -1
~>1 ~>2
EXAMPLE 13.3
A small oven measures 0·4 m by 0·5 m by 0·3 m high. The floor of
the oven receives radiation from all the walls and roof which are at
300°C and have an emissivity of 0·8. The floor is maintained at
150°C and has an emissivity of0·6. Calculate the radiation exchange.
Fig. 13.16. Electrical analogy of gas radiation involving two grey surfaces.
EXAMPLE 13.4
A flat plate solar collector has a selective surface with an absorptivity
of 0·92 and an emissivity of 0·15. The coefficient for heat loss by
convection is 0·003 kW/(m 2 K). If the area is 20m2 calculate the
rate of energy collection, and the collection efficiency at a time when
the solar insolation is 800 W /m 2, for a collector temperature of
60°C and an ambient temperature of 18°C.
PROBLEMS
~-2 = ---------
1/Et +(At /A2)(1/€2-1)
5. Two large parallel plates, 1 and 2, having emissivities on their inner faces
of 0·5 and 0·8 are maintained at 300° and 100oc respectively. A third plate
having unknown emissivities on its faces A and B is placed between the
other two plates. When face A is pointing towards plate 1, the third plate
reaches an equilibrium of 278°C. When the third plate is turned round so
that face B is pointing towards plate 1, its equilibrium temperature drops to
140°C. Determine the emissivities of the two faces A and B. (Ans. Face A,
e = 0·916, Face B, e = 0·102.) (The City University).
234 ENGINEERING HEAT TRANSFER
6. A plane surface of emissivity € 1 and at temperature T 1 radiates to a
parallel shield at temperature T 2 • The shield has an emissivity e 2 facing
the plane surface, and an emissivity e3 facing the surroundings at T3 ,
which may be assumed large. Show that the radiation from the plane
surface to the surroundings, via the shield, per unit area, is given by:
REFERENCES
I. Jakob, M. Heat Transfer, Vol. I, John Wiley and Sons, Inc., New York
(1949).
2. Planck, M. The Theory of Heat Radiation (Translation) Dover (1959).
3. Kirchhoff, G. Ostwalds Klassiker d. exakten Wissens .. 100, Leipzig (1898).
4. Lambert. J. H. Photometria (1860).
5. Sparrow, E. M. and Cess, R. D. Radiation Heat Transfer. Brooks/Cole
Publishing Company (1970).
6. Chapman, A. J. Heat Transfer, 3rd ed. The Macmillan Company, New
York (1974).
7. Oppenheim. A. K. Amer. Soc. Mech. Engs. Paper 54-A75 (1954).
8. Hottel, H. C. Notes on Radiant Heat Transmission, Chern. Eng. Dept..
M.I.T. (1951 ).
9. McAdams, W. H. Heat Transmission (Chapter 4, Hottel, H. C.), 3rd ed.,
McGraw-Hill Book Company, Inc., New York (1954).
Appendix 1: Heat Transfer
Literature
The following is a list of journals, proceedings, and bibliography
which may be consulted in order to keep abreast of the most recently
published work in heat transfer.
The International Journal of Heat and Mass Transfer, Pergamon
Press, monthly
The Journal of Mechanical Engineering Science, The Institution of
Mechnical Engineers, bi-monthly.
Journal of Heat Transfer, Transactions of the American Society of
Mechanical Engineers, Series C, quarterly.
Proceedings of the International Heat Transfer Conferences, e.g.,
4th 1970 (Paris~ 5th 1974(Tokyo), Elsevier Publishing Company,
Amsterdam.
Progress in Heat and Mass Transfer, Monograph Series of the
International Journal of Heat and Mass Transfer, Pergamon
Press.
Advances in Heat Transfer, Academic Press, New York.
Proceedings of the Heat Transfer and Fluid Mechanics Institute,
Stanford University Press, California.
Heat Bibliography, HMSO London, annual.
Reports of the National Engineering Laboratory, East Kilbride,
(available on request).
The Engineering Index, Engineering Index, Inc., New York.
Applied Science and Technology Index, The H. W. Wilson Company,
New York.
The British Technology Index, The Library Association, London.
ISMEC Bulletin, Information Service in Mechanical Engineering.
The Institution of Mechanical Engineers.
Science Abstracts A, Physics Abstracts, The Institution of Electrical
Engineers.
Science Abstracts B, Electrical and Electronic Abstracts, The
Institution of Electrical Engineers and The Institute of Electrical
and Electronic Engineers, Inc.
236
Appendix 2: Units and
Conversion Factors
Mass 1 kg = 2·2046lb
Length 1m = 3·2808 ft
Time 1 s = 2·778 x 1o- 4 h
Temperature 1 K = 1·8 aRankine
Derived SI units are:
Force 1 N = 0·2248lbf (1 newton = 1 kg m/s 2 )
Pressure 1 Pa = 14·5 x 10- 5 lbf/in 2 (1 pascal= 1 N/m 2 )
1 bar = 105 Pa = 14·5lbf/in2
Density 1 kgfm 3 = 0·06243 lb/ft 3
Specific
volume 1 m 3 fkg = 16·0179 ft 3 flb
Energy 1 J = 1 Nm; 1 kJ = 103 Nm = 737·6ft lbf
Power 1 W = 1 Nmjs; 1 kW = 737·6 ft lbf/s = 1·341 h.p.
237
N
w
00
p CP X103 /{ X103 IX
Aluminium, pure 2707 896 204 8-42 X 10- 5 206 215 229 249
Duralumin, 94-96 AI, 3-5 Cu 2787 883 164 6·68 182 194
Lead 11.370 130 34·6 2·34 33-4 31·5 29·8
N Iron, pure 7897 452 72-7 2·03 67·5 62·3 55·4 48·5 39·8
w
'CJ Iron, wrought, C < 0·5% 7849 460 58·9 1·63 57·1 51·9 48·5 45·0 36·4
Iron, cast, C ::::: 4% 7272 419 51·9 1·70
Carbon steel, C ::::: 0·5% 7833 465 53-7 1-47 51·9 48·5 45·0 41·5 34·6
Carbon steel. C = 1·5% 7753 486 36-4 0·97 36·3 36·3 34·6 32-9 31·2
Nickel steel, 10% 7945 460 26·0 0·72
Nickel steel, 30% 8073 460 12·1 0·33
Nickel steel, 50% 8266 460 13-8 0·36
Nickel steel, 70% 8506 460 26·0 0·67
Nickel steel, 90% 8762 460 46·7 1·16
Chrome steel, I % 7865 460 60·6 1·67 55-4 51·9 46·7 41·5 36·4
Chrome steel, 5 % 7833 460 39·8 1·11 38·1 36·4 36·4 32·9 29·4
Chrome steel, 10% 7785 460 31·2 0·87 31·2 31·2 29·4 29·4 31·2
Cr-Ni steel, 18% Cr, 8% Ni 7817 460 16·3 0-44 17·3 17·3 19·0 19·0 22·5
Ni-Cr steel, 20% Ni, 15% Cr 7865 460 14·0 0·39 15·1 15·1 16·3 17·3 19·0
Manganese steel, 2% 7865 460 38·1 1·05 36·4 36-4 36·4 34·6 32·9
tv
-1>-
0
Table A. I. Continued
p cP x 103 k X 103 IX
tT'l
100 200 300 400 600
zQ
(~~) (~2) (•q z
(k~K) (!:) tT'l
tT'l
~
Tungsten steel, 2% 7961 444 62-3 1·76 x w-~ 58·9 53·7 48·5 45·0 36·4
Silicon steel, 2 % 7673 460 31·2 0·89
z
Q
Copper, pure 8954 383 386 ll·2 379 374 369 364 353 ::z::
Bronze, 75 Cu, 25 Sn 8660 343 26·0 0·86 tT'l
Brass, 70 Cu, 30 Zn 8522 385 111 3-41 128 144 147 147 >
....,
German silver, 62 Cu 15 Ni, 22 Zn 8618 394 24·9 0·73 31·2 39·8 45·0 48·5 ....,
~
Constantan, 60 Cu, 40 Ni 8922 410 22·7 0·61 22·2 26·0 >
Magnesium, pure 1746 1013 171 9·71 168 163 158 z
(/)
Molybdenum 10.220 251 123 4·79 ll8 ll4 111 109 106 "l"l
Nickel, 99·9% pure 8906 446 90·0 2·27 83·1 72-7 64·0 58·9 tT'l
~
Silver, 99·9% pure 10.520 234 407 16·6 415 374 362 360
Tungsten 19,350 134 163 6·27 151 142 133 126 113
Zinc, pure 7144 384 112 4·ll 109 106 100 93-5
Tin, pure 7304 227 64·0 3-88 58·9 57·1
Adapted from Table A-1, E. R. G. Eckert and R. M. Drake, Jr., Heat and Mass Transfer, McGraw-Hill Book Company, New York (1959).
APPENDIX 3 241
cP x 103 p t k X loJ ac
CP X 1()3 p t k X 103 ~
Adapted from A. J. Chapman, Heat Transfer, The Macmillan Company, New York
(1960); L. S. Marks, Mechanical Engineers' Handbook, 5th ed., McGraw-Hill Book
Company, Inc., New York (1951); W. H. McAdams, Heat Transmission. 3rd ed.,
McGraw-Hill Book Company, Inc., New York (1954); and E. R. G. Eckert and R. M.
Drake, Jr., Heat and Mass Transfer, McGraw-Hill Book Company, Inc., New York
(1959).
APPENDIX 3 243
k(Wj(mK))
From the CIBS Guide Book A, The Chartered Institution of Building Services
Engineers, London. The above U values and thermal conductivities are a brief
extract only (used by permission of the Institution).
Table A.4. Physical Properties of some Common Low Melting Point Metals
N
Adapted from Table' 1&-1, J. G. Knudsen and D. L. Katz, Fluid Dynamics and Heat Transfer, McGraw-Hill Book Company, Inc., New York t;
(1958).
N
-1>-
0\
t p cP x t()l v (% Pr {J
I kxW
(OC) (kgfm 3 1 kJ/(kgK) (m 2/s) I kW/(mK) (m 2/s) (1/K) m
z
C'l
Water (H 2 0) z
tT1
0 1002 4218 0·179 X 10-S ()-552 13-1 X 10- 8 13-6 m
:;e
20 1001 4182 ()-101 ()-597 14·3 7-()2 I ()-18 X 10- 3
40 994·6 4178 ()-0658 ()-628 15·1 4·34
zC'l
60 985-4 4184 00477 ()-651 15·5 3-()2 =
4196 m
80 974·1 00364 0·668 16·4 2·22 >
96()-6 00294 ()-680
....,
100 4216 16·8 1·74 ....,
120 945·3 4250 00247 ()-685 1H 1·446 :;e
140 928·3 4283 00214 0·684 17·2 1·241 >
160 909·7 4342 00189 ()-680 17-3 I -o99 z
V>
()-675 1-()04 'Tj
180 889·0 4417 00173 17·2
200 866·7 4505 00160 ()-665 17·1 ()-937
m
:;e
220 842-4 4610 00149 ()-653 16·8 ()-891
240 815·7 4756 00143 ()-635 16·4 ()-871
260 785-9 4949 00137 ()-611 15·6 ()-874
280 752·5 5208 00135 ()-580 14·8 o-910
300 714·3 5728 00135 ()-540 13-2 1-()19
Methyl Chloride (CH 3 Cl)
-50 1053 1476 0·032o x to- s 0·215 13·9 x 10- 8 2·31
-40 1033 1483 0·0318 0·209 13-7 2·32
-30 1017 1492 0·0314 0·202 13-4 2·35
-20 999·4 1504 Q-0309 0·196 no 2·38
-10 981·4 1519 0·0306 0·187 12·6 2·43
0 962-4 1538 0·0302 0·178 12·1 2·49
10 942·4 1560 D-0297 0·171 11·7 2·55
20 923·3 1586 Q-0292 0·163 Il-l 2·63
30 903-1 1616 0·0287 Q-154 10·6 2·72
40 883·1 1650 0·0281 D-144 9·96 2·83 >
'"t:l
50 861·2 1689 Q-0274 0·133 9·21 2·97 '"t:l
ttl
Freon (CCI 2 F 2 )
z
t:i
-50 1547 875·0 o-o310 x 10- s Q-0675 5·01 x 10- 8 6·2 2·63 x 10- 3 ><
w
-40 1519 884·7 D-0279 0-0692 5·13 5-4
-30 1490 895·6 0·0253 0·0692 5·26 4·8
-20 1461 907·3 0·0235 0.0710 5·39 4·4
-10 1430 920·3 D-0221 0.0727 5·50 4·0
0 1397 934·5 Q-0214 0.0727 5·57 3·8
10 1364 949·6 Q-0203 0·0727 5·60 3·6
20 1330 965·9 Q-0198 0·0727 5·60 3·5
30 1295 983·5 0·0194 0·0710 5·60 3·5
40 1257 1002 Q-0191 0·0692 5·55 3·5
50 1216 1022 0·0189 0·0675 5·44 3·5
w
~
-.J
t0
00
"""
Table A.5. Continued
t p c, X J()l v k X 103 IX Pr p
("C) (kl/m 3 ) kJ/(kg K) (m 1/s) kW/(mK) (m 2 /s) (1/K)
tT1
z
C1
Glycerin (C 3 H 5 (0Hh) z
0 1276 2261 8·31 x w- 3 0·282 9·83 x to- s 84·7 X 103 tT1
tT1
10 1270 2320 3-()() Q-284 9·65 31·0 ::e
20 1264 2387 H7 Q-286 9·47 12·5 I 0·504 x w- 3 z
30 1258 2445 ()-50 Q-286 9·29 5·38 C1
40 1252 2512 Q-22 0·286 9·13 2-45 ::I:
tT1
50 1245 2583 Q-15 0·287 8·93 1·63 >
....,
....,
Ethylene glycol (C 2 H 4 (0Hh) ::e
0 1130 2294 5-75 x w- 5 Q-242 9·34 X JO-S 615 >
z
~
20 1117 2382 1·92 Q-249 9·39 204 1 o-648 x w- 3 'Tl
40 1101 2474 Q-869 Q-256 9·39 93 tT1
60 1088 2562 Q-475 Q-260 9·31 51
::e
80 1078 2650 Q-298 Q-261 9·21 32-4
100 1059 2742 Q-203 Q-263 9-()8 22-4
Engine oil (unused)
0 899 1796 4·28 X JO-l ()-147 9·11 x lO-a 47,100
20 888 1880 ()-9() ()-145 8·72 10,400 I o-102 x w- 3
40 876 1964 ()-24 ()-144 8·33 2870
60 864 2047 00839 ()-140 8-()() lOSO
80 852 2131 00375 ()-138 7·69 490
100 840 2219 00203 ()-137 7-38 276
120 829 2307 00123 ()-135 HO 175
140 817 2395 00080 0·133 6·86 116
160 806 2483 00056 ()-132 6·63 84 >
:g
Mercwy(Hg) tr1
zt1
0 13.630 140-3 00124 x 10-s 8·21 430 X lO-a 00288
20 13.580 139-4 00114 8-69 461 00249 )·82 X 10- 4 xw
so 13,510 138-6 00104 9·40 S02 00207
100 13,390 137-3 000928 1o-5 571 00162
ISO 13,260 136·5 0-()()853 11·5 635 ()-0134
200 13,150 136·1 000802 12·3 691 ()-0116
2SO 13,030 135-7 000764 13-l 740 ()-0103
316 12.8SO 134.() 000673 14.() 815 0-()()83
Adapted from Table A-3, E. R. G. Eckert and R. M. Drake, Jr., Heat and Mass Transfer, McGraw-Hill Book Company, Inc., New York (1959).
N
~
\0
N
Table A.6. Thermal Properties of Gases at Atmospheric Pressure Vl
0
T p C0 X 1()3 v k X J()l Q( Pr
Jl
(oK) (kg/m 3 ) kJ/(kg K) (m 2/s) kW/(mK) (m 2/s) Pas
Air tr1
250 1·413 1005 ()-949 X 10-S ()-0223 1·32 X 10-S 1·60 X w-s z
0·722
.300 H77 1006 1·57 ()-0262 2·22 1·85 ()-708 z0
350 ()-998 1009 2-()8 ()-0300 tr1
2·98 2·08 0·697 tr1
400 Q-883 1014 2·59 0·0337 3-76 2·29 0·689 :;tl
450 ()-783 1021 2·89 ()-0371 4·22 2·48 ()-683 -z
500 ()-705 1030 3-79 ()-0404 5·57 2·67 0
0·680
550 ()-642 1039 4·43 ()-0436 6·53 ::I:
2·85 0·680 tr1
600 ()-588 1055 5·13 0·0466 7·51 3·02 0·680 >
650 ()-543 1063 o-l
5·85 ()-0495 8·58 3·18 0·682
()-503 o-l
700 1075 6-63 ()-0523 9·67 3-33 0·684 :;tl
750 ()-471 1086 7-39 ()-0551 10·8 3-48 0·686
800 o-441 1098 8·23 ()-0578 12-() 3·63 ()-689
zen>
850 o-415 1ll0 9-()7 ()-0603 'TI
IH 3·77 ()-692 tr1
900 Q-392 ll21 9·93 ()-0628 14·3 3·90 ()-696 :;tl
950 ()-372 ll32 1()-8 ()-0653 15·5 4·02 0·699
1000 ()-352 1142 11·8 ()-0675 16·8 4·15 ()-702
1100 Q-320 ll61 13·7 ()-0723 19·5 4·40 0·706
1200 ()-295 ll79 15·7 ()-0763 22-() 4·63 0·714
1.300 ()-271 1197 17·9 ()-0803 24·8 4·85 0·722
Hydrogen
250 Q-0981 14,060 8·06 X lO-S 0·156 11·3 x w- s 7·92 x w- 6 G-713
300 0·0819 14,320 10·9 0·182 15·5 8·96 Q-706
350 0·0702 14.440 14·2 Q-206 20·3 9·95 0·697
400 O-o614 14.490 17-7 0·229 25·7 1Q-9 0·690
450 0·0546 14,500 21·6 0·251 31·6 11·8 0·682
500 0·0492 14.510 25·7 0·272 38·2 12·6 0·675
550 0·0447 14.330 30·2 0·293 45·2 13·5 0·668
600 0·0408 14,540 35·0 0·315 53·1 14·3 0·664
650 0·0349 14,570 45·5 0·351 69·0 15·9 0·659
700 0·0306 14,680 56·9 0·384 85·6 17-4 0·664 >
"'C
750 0·0272 14,820 69·0 0·412 102 18·8 0·676 "'C
800 0·0245 14.970 82·2 0·440 120 20·2 0·686 z0tr.l
850 0-()223 15.170 96·5 0·464 137 21·5 0·703
Oxygen ><
-
w
200 1·956 913·1 Q-795 X lO-S 0-()182 1·0? X JO-S 14·9 x w- 6 0·745
250 1·562 915·6 1·144 Q-0226 1·58 17·9 0·725
300 1·301 920·3 1·586 Q-0267 2·24 2Q-6 0·709
350 1-113 929·0 2·080 Q-0307 2·97 23·2 0·702
400 Q-976 942-Q 2·618 Q-0346 3-77 25·5 0·695
450 0·868 956·7 3·199 Q-0383 4·61 27·8 Q-694
500 Q-780 972-2 3·834 Q-0417 5·50 29·9 0·697
550 0·710 988·1 4·505 Q-0452 6"·44 32·0 0·700
600 Q-650 1004 5·214 Q-0483 7-40 33·9 0·704
N
01
IV
...,.
Table A.6. Continued IV
T p c, X 1()3 v k X 103 IX IJ Pr
("K) (kl/m 3 ) k.J/(kg K) (m 2 /s) kW/(m K) (m 1 /s) Pas
--
Nitrogen ttl
200 1·711 1043 o-757 x w- 5 00182 1-o2 x w- 5 12·9 x w- 6 (}747 0
z
300 1-142 1041 1·563 00262 2·21 17·8 (}713 zttl
-
400 0·854 1046 2·574 00333 3-74 22·0 0·691 ttl
500 0·682 1056 3-766 00398 5·53 25·7 0·684 ~
600 0·569 1076 5·119 (}0458 7-49 29·1 (}686 -z
700 0·493 1097 6·512 00512 9·47 32·1 0·691 0
800 (}428 1123 8·145 00561 11·7 34·8 (}700 :::
ttl
900 (}380 1146 9·106 Q-0607 13-9 37·5 (}711 >
o-,l
1000 0·341 1168 11·72 00648 16·3 4(}0 (}724 o-,l
1100 (}311 1186 13·60 (}0685 18·6 42-3 (}736 ~
1200 (}285 1204 15·61 00719 20·9 44·5 (}748 >
ztil
Carbon dioxide "!1
250 2·166 803-9 Q-581 x w- 5 00129 0·740 x w- 5 12·6 x w-" (}793 ttl
~
300 1·797 87(}9 (}832 00166 1-()6 15·0 0·770
350 1·536 9002 1-119 00205 1·48 17·2 (}755
400 1·342 942-() 1·439 00246 1-95 19·3 (}738
450 H92 979·7 1·790 00290 2-48 21·3 (}721
soo 1-()73 1013 2-167 (}0335 3-()8 23-3 (}702
sso (}974 1047 2·574 00382 3-75 25-1 (}685
600 (}894 1076 3002 00431 4·48 26·8 (}668
Carbon monoxide
250 ()-841 1043 1-128 X 10- 5 (){)214 1·51 X 10- 5 15·4 X 10- 6 Q-750
300 1-139 1042 1·567 (){)253 2-13 17·8 ()-737
350 ()-974 1043 2-()62 (){)288 2·84 2()-1 ()-728
400 ()-854 1048 2·599 ()-0323 3·61 22·2 ()-722
450 ()-762 1055 3-188 ()-()436 4·44 24·2 ()-718
500 ()-682 1063 3·819 (){)386 5-33 26·1 ()-718
550 0·620 1076 4·496 ()-0416 6·24 27·9 ()-721
600 ()-568 1088 5·206 ()-0445 7-19 29·6 ()-724
Water vapour
>
380 ()-586 2060 ()-216 X 10- 4 (){)246 2·04 X 10- 5 12·7 X 10- 6 1-()60 :g
()-554 2014 ()-242 (){)261 2·24 13-4 1·040 ttl
400 zt)
450 ()-490 1980 ()-311 (){)299 3·07 15-3 1·010
500 ()-441 1985 ()-386 (){)339 3·87 17·0 ()-996 ><
550 0·400 1997 ()-470 (){)379 4·75 18·8 ()-991 w
600 0·365 2026 ()-566 ()-()422 5·73 2()-7 ()-986
650 ()-338 2056 ()-664 (H)464 6·66 22·5 ()-995
700 0·314 2085 0·772 (){)505 7-12 24·3 HXlO
750 ()-293 2119 ()-888 (){)549 8·83 26·0 1·005
800 ()-274 2152 1-()20 (){)592 1(){) 27·9 1-()10
850 Q-258 2186 H52 ()-0637 11·3 29·7 1-()19
Adapted from Table A-4, E. R. G. Eckert and R. M. Drake, Jr., Heat and Mass Transfer, McGraw-Hill Book Company, Inc., New York (1959).
(Note: At pressures other than atmospheric, the density can be determined from the ideal gas equation, p == p/RT. Hence at any given tem-
perature p = p 0 (pjp 0 ) where Po is atmospheric pressure and Po is given in the table. k, p, and c• may be assumed independent of pressure. tv
v and 11 are inversely proportional to the density; hence at a given temperature are inversely proportional to the pressure.) ~
254 ENGINEERING HEAT TRANSFER
Ref. t Emissivity
(OC)
Aluminium:
Highly polished plate, 98·3% pure 11 237-576 Q-039--0057
Rough polish 1 100 Q-18
Commercial sheet 1 100 {)-()9
Heavily oxidized 2 93-505 Q-20-0-31
At-surfaced roofing 5 38 Q-216
Brass:
Highly polished, 73-2 Cu, 26·7 Zn 11 247-357 0028-0031
Polished 1 100 006
Rollcd plate, natural surface 10 22 006
Chromium. polished 1 100 0075
Copper:
Carefully polished electrolytic copper 6 80 Q-018
Polished 1 100 Q-052
Molten 3 1076-1278 Q-16-o-13
Iron and steel :
Steel, polished 1 100 0066
Iron, polished 12 427-1028 Q-14-0-38
Cast iron, polished 9 200 Q-21
Cast iron, newly turned lO 22 ()-44
Wrought iron, highly polished 16 38-249 Q-28
Iron plate, completely rusted lO 19 ()-69
Sheet steel, shiny oxide layer 10 24 Q-82
Steel plate, rough 5 38-372 Q-94-0-97
Cast iron, molten 15 1300-1400 Q-29
Steel, molten 7 1522-1650 Q-43-o-40
Stainless steel, polished 1 100 Q-074
Lead, grey oxidized 10 24 Q-28
Magnesium oxide 8 278-827 ()-55-Q-20
Nichrome wire, bright 14 49-1000 ()-65-Q-79
Nickel-silver, polished 1 100 0·135
Platinum filament 4 27-1230 0036-o-192
Silver, polished, pure 11 227--fJ27 ()-02....()-032
Tin, bright tinned iron lO 23 ()-043, ()-064
Tungsten filament 18 3320 Q-39
Zinc, galvanized sheet iron, fairly
bright lO 28 Q-23
APPENDIX 3 255
Ref. t Emissivity
c·q
Asbestos board 10 23 0·96
Brick:
Red, rough 10 21 (}93
Building 14 1000 0·45
Fireclay 14 1000 0·75
Magnesite, refractory 14 1000 0·38
Candle soot 17 97-272 (}952
Lampblack, other blacks 14 5{}-1000 0·96
Graphite, pressed, filed surface 8 249-516 (}98
Concrete tiles 14 1000 (}63
Enamel, white fused, on iron 10 19 0·90
Glass, smooth 10 22 0·94
Oak, planed 10 21 0·90
Flat black lacquer 5 38-94 0·96--0·98
Oil paints. 16 different, all colours 13 100 0·92-0·96
Aluminium paints, various 13 100 0·27-0·67
Radiator paint, bronze 1 100 0·51
Paper, thin, pasted on blackened plate 10 19 (}92, 0·94
Plaster, rough lime 16 1{}-87 (}91
Roofing paper 10 21 (}.91
Water (calculated from spectral data) {}-100 0·95-0·963
REFERENCES
I. Barnes, B. T., Forsythe, W. E., and Adams, E. Q. J. Opt. Soc. Amer.,
Vol. 37, 804 (1947).
2. Binkley, E. R., private communication (1933).
3. Burgess, G. K. Nat/. Bur. Stand., Bull. 6, Sci. paper 121, Ill (1909).
4. Davisson, C., and Weeks, J. R. Jr. J. OpL Soc. Amer., Vol. 8, 581 (1924).
5. Heilman, R. H. Trans. ASME, FSP 51,287 (1929).
6. Hoffman, K. Z. Physik, Vol. 14, 310 (1923).
7. Knowles, D., and Sarjant, R. J. J. Iron and Steel Inst. (London), Vol. 155,
577 (1947).
8. Pirani, M. J. Sci. Instrum., Vol. 16, 12 (1939).
9. Randolf, C. F., and Overhaltzer, M. J. Phys. Rev., Vol. 2, 144 (1913).
10. Schmidt, E. Gesundh-Ing., Beiheft 20, Reihe 1, 1-23 (1927).
256 ENGINEERING HEAT TRANSFER
()o02
()oOl
()o008
()o006
2500K
Fig. Al. Emissivity of carbon dioxide; adapted from W. H. McAdams
Heat Transfer, McGraw·Hill Book Company, 3rd ed., New York (1954);
by permission of the publishers.
257
ENGINEERING HEAT TRANSFER
()-8
()-6
()-5
()-4
o-3
()o()l
<roo!
~~K~~~K~~~~--~~~--~~-
0
PH20 PH.p PHil
Pco2 + PH 20 Pco2+ PHp Pco2 + 1 JH20
Fig. AS. Adllptetlfrom W. H. McAtlams, Heat Tra11smissioll, McGraw-Hill
Book ColllfHIIIy, New York (1954); by permusio11 of the p•blisurs. For
Iiiia ofcouttutt P002 L + P820 L, ilf m btu, 1-1·5 m btu, :Z-HJ m btu, 3-0·6
m btu, 4-0·5 m btu, 6-0•:Z m hr, 7-0·1 m btu.
260 ENGINEERING HEAT TRANSFER
Then the correction factors are applied as in the case of emissivity
determination, and finally the mutual absorption correction is
similarly made.
ExAMPLE
absorptivity laminar 78
definition of 209 separation of 136
of black body 210 sub-layer 78, 107
of grey body 213 thermal 80
Akers, W. W. 149 thickness of 90
algebra, configuration factor, in thickness of 8 7-8
radiation 222-4 turbulent 78-9
analogy, Reynolds 101-11, see velocity distribution in 79
also Reynolds analogy velocity distribution in 79-87
analogy in complex flow 137 boundary mesh points 47-9
analogy of conduction 52-5 British Nuclear Fuels, plc 158
analogy of radiation 224-8, Buckingham's pi theorem 111
230-1 building materials, thermal con-
anisotropic materials 10 ductivities of 243
Farber, E. A. 150
Deans, H. A. 149 Fenner, R. T. 52
INDEX 263
film, equivalent conducting 91 Grashof number, definition of
fin analysis, program list 172 126
finite difference relationships grey body 212
in steady state conduction 42, emission 213
47-50 Griffith, P. 151
in transient conduction 61 , Grimison, E. D. 140
63,67
finned surface heat, definition of 2
equivalent effectiveness of heat exchangers
164-5 basic types of 176-7
overall coefficient of 165-8 batch 202-3
fins cross flow 177, 191-4
conduction in 158-62,170-2 determination of performance
effectiveness of 164-5 of 181-203
limit of usefulness of 164 in counter and parallel flow
numerical relationships in 181-91
170-3 in cross flow 191-4
temperature distribution in effectiveness of
160-3 at limiting value of capacity
fire-resistant door analysis 64-7 ratio 190
program list 65 in counter flow 187-8
Firman, E. C. 152 in cross flow 190-1
forced convection in parallel flow 189
definition of 4 general discussion of 176
dimensional analysis of 111-15 in-line 176-7
empirical results of 115, 118, thermal wheel 116, 194-9
139-40 transfer units 18 5-91
in laminar flow 78-98 heat flux 8
flat plates 87-92 heat sink, transistor 163
in tubes 92-8 heat transfer across boundary layer
in tubulent flow 117-19 in laminar flow 103
forces, buoyancy 4, 124-5 in turbulent flow 103
Fourier number, definition of 62 heat transfer coefficient 19, 28-9,
Fourier's law 3, 8 165-8, see also convection
friction coefficient coefficient and Nusselt number
for flat plates 102 in complex flow system 137-9
for tubes 103 in fins 158-64
in liquid metals 118-19
Gardner, G. C. 152 in uniform temperature system
gas emissivities 257-60 58-61
gas radiation, non-luminous 228- heat transfer in building structures
31 20-1
gases, thermal properties of horizontal surfaces, natural con-
250-3 vection 12 7-8
Gaussian elimination method 4 7 Hottel, H. C. 224, 229
Gauss-Siedel iterative method 52 Hsu,S.T. 67,127,139,153
Graetz number 115
graphical solution of transient insulation, critical thickness of
conduction 67 28-30
264 INDEX