PART 3
INTRODUCTION TO ENGINEERING HEAT TRANSFER
HT-1
Introduction to Engineering Heat Transfer
These notes provide an introduction to engineering heat transfer. Heat transfer processes set limits
to the performance of aerospace components and systems and the subject is one of an enormous
range of application. The notes are intended to describe the three types of heat transfer and provide
basic tools to enable the readers to estimate the magnitude of heat transfer rates in realistic aerospace
applications. There are also a number of excellent texts on the subject; some accessible references
which expand the discussion in the notes are listen in the bibliography.
HT-2
Table of Tables
Table 2.1: Thermal conductivity at room temperature for some metals and non-metals ............. HT-7
Table 2.2: Utility of plane slab approximation..........................................................................HT-17
Table 9.1: Total emittances for different surfaces [from: A Heat Transfer Textbook, J. Lienhard ]HT-63
HT-3
Table of Figures
Figure 1.1: Conduction heat transfer ......................................................................................... HT-5
Figure 2.1: Heat transfer along a bar ......................................................................................... HT-6
Figure 2.2: One-dimensional heat conduction ........................................................................... HT-8
Figure 2.3: Temperature boundary conditions for a slab............................................................ HT-9
Figure 2.4: Temperature distribution through a slab .................................................................HT-10
Figure 2.5: Heat transfer across a composite slab (series thermal resistance) ............................HT-11
Figure 2.6: Heat transfer for a wall with dissimilar materials (Parallel thermal resistance)........HT-12
Figure 2.7: Heat transfer through an insulated wall ..................................................................HT-11
Figure 2.8: Temperature distribution through an insulated wall................................................HT-13
Figure 2.9: Cylindrical shell geometry notation........................................................................HT-14
Figure 2.10: Spherical shell......................................................................................................HT-17
Figure 3.1: Turbine blade heat transfer configuration ...............................................................HT-18
Figure 3.2: Temperature and velocity distributions near a surface. ...........................................HT-19
Figure 3.3: Velocity profile near a surface................................................................................HT-20
Figure 3.4: Momentum and energy exchange in turbulent flow. ...............................................HT-20
Figure 3.5: Heat exchanger configurations ...............................................................................HT-23
Figure 3.6: Wall with convective heat transfer .........................................................................HT-25
Figure 3.7: Cylinder in a flowing fluid .....................................................................................HT-26
Figure 3.8: Critical radius of insulation ....................................................................................HT-29
Figure 3.9: Effect of the Biot Number [hL / k
body
] on the temperature distributions in the solid and in
the fluid for convective cooling of a body. Note that k
body
is the thermal conductivity of the
body, not of the fluid.........................................................................................................HT-31
Figure 3.10: Temperature distribution in a convectively cooled cylinder for different values of Biot
number, Bi; r
2
/ r
1
= 2 [from: A Heat Transfer Textbook, John H. Lienhard] .....................HT-32
Figure 4.1: Slab with heat sources (a) overall configuration, (b) elementary slice.....................HT-32
Figure 4.2: Temperature distribution for slab with distributed heat sources ..............................HT-34
Figure 5.1: Geometry of heat transfer fin .................................................................................HT-35
Figure 5.2: Element of fin showing heat transfer......................................................................HT-36
Figure 5.3: The temperature distribution, tip temperature, and heat flux in a straight one-
dimensional fin with the tip insulated. [From: Lienhard, A Heat Transfer Textbook, Prentice-
Hall publishers].................................................................................................................HT-40
Figure 6.1: Temperature variation in an object cooled by a flowing fluid .................................HT-41
Figure 6.2: Voltage change in an R-C circuit............................................................................HT-42
Figure 8.1: Concentric tube heat exchangers. (a) Parallel flow. (b) Counterflow.......................HT-44
Figure 8.2: Cross-flow heat exchangers. (a) Finned with both fluids unmixed. (b) Unfinned with one
fluid mixed and the other unmixed ....................................................................................HT-45
Figure 8.3: Geometry for heat transfer between two fluids .......................................................HT-45
Figure 8.4: Counterflow heat exchanger...................................................................................HT-46
Figure 8.5: Fluid temperature distribution along the tube with uniform wall temperature .........HT-46
Figure 9.1: Radiation Surface Properties ..................................................................................HT-52
Figure 9.2: Emissive power of a black body at several temperatures - predicted and observed..HT-53
Figure 9.3: A cavity with a small hole (approximates a black body).........................................HT-54
Figure 9.4: A small black body inside a cavity .........................................................................HT-54
Figure 9.5: Path of a photon between two gray surfaces ...........................................................HT-55
HT-4
Figure 9.6: Thermocouple used to measure temperature...........................................................HT-59
Figure 9.7: Effect of radiation heat transfer on measured temperature. .....................................HT-59
Figure 9.8: Shielding a thermocouple to reduce radiation heat transfer error ............................HT-60
Figure 9.9: Radiation between two bodies................................................................................HT-60
Figure 9.10: Radiation between two arbitrary surfaces .............................................................HT-61
Figure 9.11: Radiation heat transfer for concentric cylinders or spheres ...................................HT-62
Figure 9.12: View Factors for Three - Dimensional Geometries [from: Fundamentals of Heat
Transfer, F.P. Incropera and D.P. DeWitt, John Wiley and Sons]......................................HT-64
Figure 9.13: Fig. 13.4--View factor for aligned parallel rectangles [from: Fundamentals of Heat
Transfer, F.P. Incropera and D.P. DeWitt, John Wiley and Sons]......................................HT-65
Figure 9.14: Fig 13.5--View factor for coaxial parallel disk [from: Fundamentals of Heat Transfer,
F.P. Incropera and D.P. DeWitt, John Wiley and Sons] .....................................................HT-65
Figure 9.15: Fig 13.6--View factor for perpendicular rectangles with a common edge .............HT-66
HT-5
1.0 Heat Transfer Modes
Heat transfer processes are classified into three types. The first is conduction, which is defined
as transfer of heat occurring through intervening matter without bulk motion of the matter. Figure
1.1 shows the process pictorially. A solid (a block of metal, say) has one surface at a high
temperature and one at a lower temperature. This type of heat conduction can occur, for example,
through a turbine blade in a jet engine. The outside surface, which is exposed to gases from the
combustor, is at a higher temperature than the inside surface, which has cooling air next to it. The
level of the wall temperature is critical for a turbine blade.
T
high
T
low
Solid
Heat “flows” to right ( q
&
)
Figure 1.1: Conduction heat transfer
The second heat transfer process is convection, or heat transfer due to a flowing fluid. The
fluid can be a gas or a liquid; both have applications in aerospace technology. In convection heat
transfer, the heat is moved through bulk transfer of a non-uniform temperature fluid.
The third process is radiation or transmission of energy through space without the necessary
presence of matter. Radiation is the only method for heat transfer in space. Radiation can be
important even in situations in which there is an intervening medium; a familiar example is the heat
transfer from a glowing piece of metal or from a fire.
Muddy points
How do we quantify the contribution of each mode of heat transfer in a given situation?
(MP HT.1)
2.0 Conduction Heat Transfer
We will start by examining conduction heat transfer. We must first determine how to relate the
heat transfer to other properties (either mechanical, thermal, or geometrical). The answer to this is
rooted in experiment, but it can be motivated by considering heat flow along a "bar" between two
heat reservoirs at T
A
, T
B
as shown in Figure 2.1. It is plausible that the heat transfer rate Q
&
, is a
HT-6
function of the temperature of the two reservoirs, the bar geometry and the bar properties. (Are there
other factors that should be considered? If so, what?). This can be expressed as
Q
&
= f
1
(T
A
, T
B
, bar geometry, bar properties) (2.1)
It also seems reasonable to postulate that Q
&
should depend on the temperature difference T
A
- T
B
. If
T
A
– T
B
is zero, then the heat transfer should also be zero. The temperature dependence can therefore
be expressed as
Q
&
= f
2
[ (T
A
- T
B
), T
A
, bar geometry, bar properties] (2.2)
L
T
B
T
A
Q
&
Figure 2.1: Heat transfer along a bar
An argument for the general form of f
2
can be made from physical considerations. One
requirement, as said, is f
2
= 0 if T
A
= T
B
. Using a MacLaurin series expansion, as follows:
f( T) f(0)
f
(T)
T
0
?
?
?=+
?
?
+L (2.3)
If we define ?T = T
A
– T
B
and f = f
2
, we find that (for small T
A
– T
B
),
f(T T ) Q f(0)
f
(T T )
TT .
2A B 2
2
AB
T
A
T
B
0
AB
?== +
?
??
?
()
+
?
?=
L (2.4)
We know that f
2
(0) = 0 . The derivative evaluated at T
A
= T
B
(thermal equilibrium) is a measurable
property of the bar. In addition, we know that QTT
f
TT
AB
2
AB
?
>>
?
??
()
>00 if or . It also seems
reasonable that if we had two bars of the same area, we would have twice the heat transfer, so that
we can postulate that Q
&
is proportional to the area. Finally, although the argument is by no means
rigorous, experience leads us to believe that as L increases Q
&
should get smaller. All of these lead
to the generalization (made by Fourier in 1807) that, for the bar, the derivative in equation (2.4) has
the form
HT-7
?
??
()
=
?=
f
TT
kA
L
2
AB
T
A
T
B
0
. (2.5)
In equation (2.5), k is a proportionality factor that is a function of the material and the
temperature, A is the cross-sectional area and L is the length of the bar. In the limit for any
temperature difference ?T across a length ?x as both L, T
A
- T
B
→ 0, we can say
()()
dx
dT
kA
L
TT
kA
L
TT
kAQ
ABBA
?=
?
?=
?
=
&
. (2.6)
A more useful quantity to work with is the heat transfer per unit area, defined as
q
A
Q
&
&
= . (2.7)
The quantity
q&
is called the heat flux and its units are Watts/m
2
. The expression in (2.6) can
be written in terms of heat flux as
dx
dT
kq ?=& . (2.8)
Equation 2.8 is the one-dimensional form of Fourier's law of heat conduction. The
proportionality constant k is called the thermal conductivity. Its units are W / m-K. Thermal
conductivity is a well-tabulated property for a large number of materials. Some values for familiar
materials are given in Table 1; others can be found in the references. The thermal conductivity is a
function of temperature and the values shown in Table 1 are for room temperature.
Table 2.1: Thermal conductivity at room temperature for some metals and non-metals
Metals Ag Cu Al Fe Steel
k [W/m-K] 420 390 200 70 50
Non-metals H
2
0 Air Engine oil H
2
Brick Wood Cork
k [W/m-K] 0.6 0.026 0.15 0.18 0.4 -0 .5 0.2 0.04
HT-8
2.1 Steady-State One-Dimensional Conduction
()xQ
&
()dxxQ +
&
dx
x
Insulated
(no heat transfer)
Figure 2.2: One-dimensional heat conduction
For one-dimensional heat conduction (temperature depending on one variable only), we can
devise a basic description of the process. The first law in control volume form (steady flow energy
equation) with no shaft work and no mass flow reduces to the statement that Q
&
Σ for all surfaces = 0
(no heat transfer on top or bottom of figure 2.2). From equation (2.8), the heat transfer rate in at the
left (at x) is
˙
Qx kA
dT
dx
x
()=?
?
?
?
?
. (2.9)
The heat transfer rate on the right is
˙˙
˙
Qx dx Qx
dQ
dx
dx
x
+()= ()++L. (2.10)
Using the conditions on the overall heat flow and the expressions in (2.9) and (2.10)
˙˙
˙
QQ
Q
xx
d
dx
xdx 0()? ()+ () +
?
?
?
?
?
?
=L . (2.11)
Taking the limit as dx approaches zero we obtain
dQ x
dx
˙
()
=0 , (2.12a)
or
HT-9
d
dx
kA
dT
dx
0
?
?
?
?
?
?
= . (2.12b)
If k is constant (i.e. if the properties of the bar are independent of temperature), this reduces to
d
dx
A
dT
dx
0
?
?
?
?
?
?
= (2.13a)
or (using the chain rule)
dT
dx
1
A
dA
dx
dT
dx
0
2
2
+
?
?
?
?
?
?
= . (2.13b)
Equations (2.13a) or (2.13b) describe the temperature field for quasi-one-dimensional steady state
(no time dependence) heat transfer. We now apply this to some examples.
Example 2.1: Heat transfer through a plane slab
x
T = T
1
T = T
2
Slab
x = 0 x = L
Figure 2.3: Temperature boundary conditions for a slab
For this configuration, the area is not a function of x, i.e. A = constant. Equation (2.13) thus became
0
2
2
=
dx
Td
. (2.14)
Equation (2.14) can be integrated immediately to yield
a
dx
dT
= (2.15)
HT-10
and baxT += . (2.16)
Equation (2.16) is an expression for the temperature field where a and b are constants of integration.
For a second order equation, such as (2.14), we need two boundary conditions to determine a and b.
One such set of boundary conditions can be the specification of the temperatures at both sides of the
slab as shown in Figure 2.3, say T (0) = T
1
; T (L) = T
2
.
The condition T (0) = T
1
implies that b = T
1
. The condition T
2
= T (L) implies that T
2
= aL + T
1
, or
L
TT
a
12
?
= .
With these expressions for a and b the temperature distribution can be written as
Tx T
TT
L
x
1
21
()
=+
??
?
?
?
?
?
. (2.17)
This linear variation in temperature is shown in Figure 2.4 for a situation in which T
1
> T
2
.
T
2
T
1
T
x
Figure 2.4: Temperature distribution through a slab
The heat flux
q&
is also of interest. This is given by
()
constant
12
=
?
?=?=
L
TT
k
dx
dT
kq&
. (2.18)
Muddy points
How specific do we need to be about when the one-dimensional assumption is valid? Is it
enough to say that dA/dx is small? (MP HT.2)
Why is the thermal conductivity of light gases such as helium (monoatomic) or hydrogen
(diatomic) much higher than heavier gases such as argon (monoatomic) or nitrogen
(diatomic)? (MP HT.3)
HT-11
2.2 Thermal Resistance Circuits
There is an electrical analogy with conduction heat transfer that can be exploited in problem
solving. The analog of Q
&
is current, and the analog of the temperature difference, T
1
- T
2
, is voltage
difference. From this perspective the slab is a pure resistance to heat transfer and we can define
R
TT
Q
21
?
=
&
(2.19)
where R = L/kA, the thermal resistance. The thermal resistance R increases as L increases, as A
decreases, and as k decreases.
The concept of a thermal resistance circuit allows ready analysis of problems such as a composite
slab (composite planar heat transfer surface). In the composite slab shown in Figure 2.5, the heat
flux is constant with x. The resistances are in series and sum to R = R
1
+ R
2.
If T
L
is the temperature
at the left, and T
R
is the temperature at the right, the heat transfer rate is given by
21
RR
TT
R
TT
Q
RLRL
+
?
=
?
=
&
. (2.20)
1 2
T
L
Q
&
T
R
R
1
R
2
x
Figure 2.5: Heat transfer across a composite slab (series thermal resistance)
Another example is a wall with a dissimilar material such as a bolt in an insulating layer. In
this case, the heat transfer resistances are in parallel. Figure 2.6 shows the physical configuration,
the heat transfer paths and the thermal resistance circuit.
HT-12
R
2
Q
&
k
2
k
1
k
1
R
1
model
Figure 2.6: Heat transfer for a wall with dissimilar materials (Parallel thermal resistance)
For this situation, the total heat flux Q
&
is made up of the heat flux in the two parallel paths:
21
QQQ
&&&
+= with the total resistance given by:
21
111
RRR
+= . (2.21)
More complex configurations can also be examined; for example, a brick wall with insulation
on both sides.
T
2
Brick
0.1 m
T
1
= 150 °C
T
4
= 10 °C
T
2
T
3
Insulation
0.03 m
T
3
T
1
T
4
R
1
R
2
R
3
Figure 2.7: Heat transfer through an insulated wall
The overall thermal resistance is given by
33
3
22
2
11
1
321
Ak
L
Ak
L
Ak
L
RRRR ++=++=
.
(2.22)
Some representative values for the brick and insulation thermal conductivity are:
HT-13
k
brick
= k
2
= 0.7 W/m-K
k
insulation
= k
1
= k
3
= 0.07 W/m-K
Using these values, and noting that A
1
= A
2
= A
3
= A, we obtain:
K/W m 0.42
K W/m0.07
m 0.03
2
1
1
31
====
k
L
ARAR
K/W m 0.14
K W/m0.7
m 0.1
2
2
2
2
===
k
L
AR .
This is a series circuit so
2
2
41
W/m142
K/W m 0.98
K 140
hroughout constant t ==
?
===
RA
TT
A
Q
q
&
&
x
41
4
TT
TT
?
?
0
1.0
1
2
3
4
Figure 2.8: Temperature distribution through an insulated wall
The temperature is continuous in the wall and the intermediate temperatures can be found
from applying the resistance equation across each slab, since Q
&
is constant across the slab. For
example, to find T
2
:
2
1
21
W/m142=
?
=
AR
TT
q&
This yields T
1
– T
2
= 60 K or T
2
= 90 °C.
The same procedure gives T
3
= 70 °C. As sketched in Figure 2.8, the larger drop is across the
insulating layer even though the brick layer is much thicker.
Muddy points
What do you mean by continuous? (MP HT.4)
Why is temperature continuous in the composite wall problem? Why is it continuous at the
interface between two materials? (MP HT.5)
HT-14
Why is the temperature gradient dT/dx not continuous? (MP HT.6)
Why is ?T the same for the two elements in a parallel thermal circuit? Doesn't the relative
area of the bolt to the wood matter? (MP HT.7)
2.3 Steady Quasi-One-Dimensional Heat Flow in Non-Planar Geometry
The quasi one-dimensional equation that has been developed can also be applied to non-planar
geometries. An important case is a cylindrical shell, a geometry often encountered in situations
where fluids are pumped and heat is transferred. The configuration is shown in Figure 2.9.
r
1
control volume
r
1
r
2
r
2
Figure 2.9: Cylindrical shell geometry notation
For a steady axisymmetric configuration, the temperature depends only on a single coordinate (r)
and Equation (2.12b) can be written as
k
d
dr
Ar
dT
dr
0
()
?
?
?
?
?
?
= (2.23)
or, since A = 2π r,
d
dr
r
dT
dr
0
?
?
?
?
?
?
= . (2.24)
The steady-flow energy equation (no flow, no work) tells us that
outin
QQ
&&
=
or
0=
dr
Qd
&
(2.25)
The heat transfer rate per unit length is given by
Qk2 r
dT
dr
?
=? ? π .
HT-15
Equation (2.24) is a second order differential equation for T. Integrating this equation once gives
a
dr
dT
r = . (2.26)
where a is a constant of integration. Equation (2.26) can be written as
r
dr
adT = (2.27)
where both sides of equation (2.27) are exact differentials. It is useful to cast this equation in terms
of a dimensionless normalized spatial variable so we can deal with quantities of order unity. To do
this, divide through by the inner radius, r
1
()
()
1
1
/
/
rr
rrd
adT = (2.28)
Integrating (2.28) yields
Ta
r
r
b
1
=
?
?
?
?
?
?
+ln . (2.29)
To find the constants of integration a and b, boundary conditions are needed. These will be taken to
be known temperatures T
1
and T
2
at r
1
and r
2
respectively. Applying T = T
1
at r = r
1
gives T
1
= b.
Applying T = T
2
at r = r
2
yields
1
1
2
2
ln T
r
r
aT += ,
or
()
12
12
/ln rr
TT
a
?
= .
The temperature distribution is thus
()
()
()
1
12
1
12
/ln
/ln
T
rr
rr
TTT +?=
. (2.30)
As said, it is generally useful to put expressions such as (2.30) into non-dimensional and
normalized form so that we can deal with numbers of order unity (this also helps in checking
whether results are consistent). If convenient, having an answer that goes to zero at one limit is also
useful from the perspective of ensuring the answer makes sense. Equation (2.30) can be put in non-
dimensional form as
HT-16
()
()
12
1
12
1
/ln
/ln
rr
rr
TT
TT
=
?
?
. (2.31)
The heat transfer rate, Q
&
, is given by
()
()
()
()
12
21
112
12
1
/ln
2
1
/ln
2
rr
TTk
rrr
TT
kr
dr
dT
kAQ
?
=
?
?=?= ππ
&
per unit length. The thermal resistance R is given by
()
k
rr
R
π2
/ln
12
= (2.32)
R
TT
Q
21
?
=
&
.
The cylindrical geometry can be viewed as a limiting case of the planar slab problem. To
make the connection, consider the case when
1
1
12
<<
?
r
rr
. From the series expansion for ln (1 + x)
we recall that
ln 1+ x x -
x
2
+
x
3
+
23
()
≈ K (2.33)
(Look it up, try it numerically, or use the binomial theorem on the series below and integrate term by
term. K++?=
+
2
1
1
1
xx
x
)
The logarithms in Equation (2.31) can thus be written as
ln ln1
rr
r
rr
r
r
r
rr
r
1
1
1
1
2
1
21
1
+
??
?
?
?
?
?
?
?
?
?
and (2.34)
in the limit of (r
2
– r
1
) << r
1
. Using these expressions in equation (2.30) gives
()
()
()
1
12
1
12
T
rr
rr
TTT +
?
?
?=
. (2.35)
With the substitution of r – r
1
= x, and r
2
– r
1
= L we obtain
()
L
x
TTTT
121
?+= (2.36)
HT-17
which is the same as equation (2.17). The plane slab is thus the limiting case of the cylinder if (r -
r
1
) / r << 1, where the heat transfer can be regarded as taking place in (approximately) a planar slab.
To see when this is appropriate, consider the expansion
()
x
x+1ln
, which is the ratio of heat flux for a
cylinder and a plane slab.
Table 2.2: Utility of plane slab approximation
x .1 .2 .3 .4 .5
()
x
x+1ln .95 .91 .87 .84 .81
For < 10% error, the ratio of thickness to inner radius should be less than 0.2, and for 20% error, the
thickness to inner radius should be less than 0.5.
A second example is the spherical shell with specified temperatures T (r
1
) = T
1
and T (r
2
) =
T
2
, as sketched in Figure 2.10.
r
1
T
2
T
1
r
2
Figure 2.10: Spherical shell
The area is now
2
4)( rrA π= , so the equation for the temperature field is
d
dr
r
dT
dr
0
2
?
?
?
?
?
?
= . (2.37)
Integrating equation (2.37) once yields
2
/ ra
dr
dT
= . (2.38)
Integrating again gives
HT-18
b
r
a
T +?=
or, normalizing the spatial variable
()
b
rr
a
T +
′
=
1
/
(2.39)
where a′ and b are constants of integration. As before, we specify the temperatures at r = r
1
and r =
r
2
. Use of the first boundary condition gives () baTrT +
′==
11
. Applying the second boundary
condition gives
()
()
b
rr
a
TrT +
′
==
12
22
/
Solving for a′ and b,
.
/1
/1
21
21
1
21
21
rr
TT
Tb
rr
TT
a
?
?
?=
?
?
=′
(2.40)
In non-dimensional form the temperature distribution is thus:
()
()
21
1
21
1
/1
/1
rr
rr
TT
TT
?
?
=
?
?
(2.41)
HT-19
3.0 Convective Heat Transfer
The second type of heat transfer to be examined is convection, where a key problem is
determining the boundary conditions at a surface exposed to a flowing fluid. An example is the wall
temperature in a turbine blade because turbine temperatures are critical as far as creep (and thus
blade) life. A view of the problem is given in Figure 3.1, which shows a cross-sectional view of a
turbine blade. There are three different types of cooling indicated, all meant to ensure that the metal
is kept at a temperature much lower than that of the combustor exit flow in which the turbine blade
operates. In this case, the turbine wall temperature is not known and must be found as part of the
solution to the problem.
Figure 3.1: Turbine blade heat transfer configuration
To find the turbine wall temperature, we need to analyze convective heat transfer, which means we
need to examine some features of the fluid motion near a surface. The conditions near a surface are
illustrated schematically in Figure 3.2.
T
c
∞
Ve locity
distribution;
c = 0 at surface
c (velocity)
δ′
T
w
T
∞
y
y
Figure 3.2: Temperature and velocity distributions near a surface.
HT-20
In a region of thickness δ′, there is a thin "film" of slowly moving fluid through which most of the
temperature difference occurs. Outside this layer, T is roughly uniform (this defines δ′). The heat
flux can thus be expressed as
q
Q
A
kT T
d
w
?
?
∞
==
?
()
′
(3.1)
It cannot be emphasized enough that this is a very crude picture. The general concept, however, is
correct, in that close to the wall, there is a thin layer in which heat is transferred basically by
conduction. Outside of this region is high mixing. The difficulty is that the thickness of the layer is
not a fluid property. It depends on velocity (Reynolds number), structure of the wall surface,
pressure gradient and Mach number. Generally δ′ is not known and needs to be found and it is
customary to calculate the heat transfer using [k
fluid
/ δ′]. This quantity has the symbol h and is
known as the convective heat transfer coefficient. The units of h are W/m
2
K. The convective heat
transfer coefficient is defined by
q
Q
A
hT T
w
?
?
∞
== ?
()
(3.2)
Equation 3.2 is often called Newton’s Law of Cooling. For many situations of practical interest, the
quantity h is still known mainly through experiments.
Muddy points
How do we know that δ' is not a fluid property? (MP HT.8)
3.1 The Reynolds Analogy
We describe the physical mechanism for the heat transfer coefficient in a turbulent boundary layer
because most aerospace vehicle applications have turbulent boundary layers. The treatment closely
follows that in Eckert and Drake (1959). Very near the wall, the fluid motion is smooth and laminar,
and molecular conduction and shear are important. The shear stress, τ, at a plane is given by
τμ =
dy
dc
(where μ is the dynamic viscosity), and the heat flux by
dy
dT
kq ?=&
. The latter is the same
expression that was used for a solid. The boundary layer is a region in which the velocity is lower
than the free stream as shown in Figures 3.2 and 3.3. In a turbulent boundary layer, the dominant
mechanisms of shear stress and heat transfer change in nature as one moves away from the wall.
HT-21
plane
Figure 3.3: Velocity profile near a surface
As one moves away from the wall (but still in the boundary layer), the flow is turbulent. The fluid
particles move in random directions and the transfer of momentum and energy is mainly through
interchange of fluid particles, shown schematically in Figure 3.4.
2
2
a
a
1
1
m′ c
p
T′
m′ c
p
T
Figure 3.4: Momentum and energy exchanges in turbulent flow.
With reference to Figure 3.4, because of the turbulent velocity field, a fluid mass m′ penetrates the
plane aa per unit time and unit area. In steady flow, the same amount crosses aa from the other side.
Fluid moving up transports heat m′ c
p
T. Fluid moving down transports m′ c
p
T′ downwards. If T >
T′, there is a turbulent downwards heat flow, ˙q
t
, given by ˙qmcTT
tp
= ′′?() that results.
Fluid moving up also has momentum m′ c and fluid moving down has momentum m′ c′. The net
flux of momentum down per unit area and time is therefore m′ (c′ - c). This net flux of momentum
per unit area and time is a force per unit area or stress, given by
tm'cc
t
= ′?() (3.3)
Based on these considerations, the relation between heat flux and shear stress at plane aa is
HT-22
˙q
t
=
′?
′?
?
?
?
?
τ
tp
c
TT
cc
(3.4)
or (again approximately)
˙q
ttp
c
dT
dc
=τ (3.5)
since the locations of planes 1-1 and 2-2 are arbitrary. For the laminar region, the heat flux towards
the wall is
˙q=τ
μ
kdT
dc
The same relationship is applicable in laminar or turbulent flow if
p
c
k
=
μ
or, expressed slightly
differently,
1
/
/
===
α
υ
ρ
ρμ
p
p
ckk
c
where υ is the kinematic viscosity, and α is the thermal diffusivity.
The quantity μc
p
/k is known as the Prandtl number (Pr), after the man who first presented the idea of
the boundary layer and was one of the pioneers of modern fluid mechanics. For gases, Prandtl
numbers are in fact close to unity and for air Pr = 0.71 at room temperature. The Prandtl number
varies little over a wide range of temperatures; approximately 3% from 300-2000 K.
We want a relation between the values at the wall (at which T = T
w
and c = 0) and those in the free
stream. To get this, we integrate the expression for dT from the wall to the free stream
dT
1
c
dc
p
=
˙q
τ
(3.6)
where the relation between heat transfer and shear stress has been taken as the same for both the
laminar and the turbulent portions of the boundary layer. The assumption being made is that the
mechanisms of heat and momentum transfer are similar. Equation (3.6) can be integrated from the
wall to the freestream (conditions "at ∞"):
dT
c
dc
p
ww
1
∞∞
∫∫
=
?
?
?
?
˙q
τ
(3.7)
where
˙q
τ
and c
p
are assumed constant.
HT-23
Carrying out the integration yields
TT
c
c
w
w
wp
∞
∞
?=
˙q
τ
(3.8)
where c
∞
is the velocity and c
p
is the specific heat.
In equation (3.8),
w
q& is the heat flux to the wall
and τ
w
is the shear stress at the wall. The relation between skin friction (shear stress) at the wall and
heat transfer is thus
˙q
w
pw
w
2
cT Tc
c
ρ
τ
ρ
∞∞ ∞
∞∞
?()
= . (3.9)
The quantity
τ
ρ
w
2
1/2 c
∞∞
is known as the skin friction coefficient and is denoted by C
f
. The skin
friction coefficient has been tabulated (or computed) for a large number of situations. If we define a
non-dimensional quantity
˙q
w
pw
w
pw p
cT Tc
hT T
cT Tc
h
ccρρρ
∞∞ ∞
∞
∞∞ ∞∞∞
?()
=
?()
?()
= = St,
known as the Stanton Number, we can write an expression for the heat transfer coefficient, h as
hcc
C
2
p
f
≈
∞∞
ρ . (3.10)
Equation (3.10) provides a useful estimate of h, or
w
q& , based on knowing the skin friction, or drag.
The direct relationship between the Stanton Number and the skin friction coefficient is
St =
C
2
f
The relation between the heat transfer and the skin friction coefficient
˙q
w
wp w
cT T
c
≈
?()
∞
∞
τ
is known as the Reynolds analogy between shear stress and heat transfer. The Reynolds analogy is
extremely useful in obtaining a first approximation for heat transfer in situations in which the shear
stress is "known".
HT-24
An example of the use of the Reynolds analogy is in analysis of a heat exchanger. One type of heat
exchanger has an array of tubes with one fluid flowing inside and another fluid flowing outside, with
the objective of transferring heat between them. To begin, we need to examine the flow resistance
of a tube. For fully developed flow in a tube, it is more appropriate to use an average velocity c and
a bulk temperature T
B
. Thus, an approximate relation for the heat transfer is
˙q
wwp
Bw
c
TT
c
≈
?
τ . (3.11)
The fluid resistance (drag) is all due to shear forces and is given by τ
w
A
w
= D, where A
w
is the tube
“wetted” area (perimeter x length). The total heat transfer, Q
&
, is Aq
w
&
w
, so that
c
TT
cDQ
WB
p
?
=
&
(3.12)
The power, P, to drive the flow through a resistance is given by the product of the drag and the
velocity, cD , so that
˙
Q
P
cT T
c
pB w
=
?
()
2
(3.13)
The mass flow rate is given by
Acm ρ=&
where A is the cross sectional area. For given mass flow
rate and overall heat transfer rate, the power scales as c
2
or as 1/A
2
, i.e.
P
cT T A
2
pB w
2
∝
?()
˙
˙Qm
2
1
ρ
(3.14)
Equations (3.13) and (3.14) show that to decrease the power dissipated, we need to decrease c ,
which can be accomplished by increasing the cross-sectional area. Two possible heat exchanger
configurations are sketched in Figure 3.5; the one on the right will have a lower loss.
heat exchangerheat exchanger
high loss
lower loss
m
?
diffuser
Figure 3.5: Heat exchanger configurations
HT-25
To recap, there is an approximate relation between skin friction (momentum flux to the wall) and
heat transfer called the Reynolds analogy that provides a useful way to estimate heat transfer rates in
situations in which the skin friction is known. The relation is expressed by
St =
C
2
f
(3.15a)
or
flux momentum convected
wallflux to momentum
fluxheat convected
wallflux toheat
= (3.15b)
or
˙q
w
pw
w
2
cc T T
c
ρ
τ
ρ
∞∞ ∞
∞∞
?()
≈ (3.15c)
The Reynolds analogy can be used to give information about scaling of various effects as well as
initial estimates for heat transfer. It is emphasized that it is a useful tool based on a hypothesis about
the mechanism of heat transfer and shear stress and not a physical law.
Muddy points
What is the "analogy" that we are discussing? Is it that the equations are similar? (MP HT.9)
In what situations does the Reynolds analogy "not work"? (MP HT.10)
3.2 Combined Conduction and Convection
We can now analyze problems in which both conduction and convection occur, starting with a wall
cooled by flowing fluid on each side. As discussed, a description of the convective heat transfer can
be given explicitly as
˙
˙
Q
q
A
hT T
w
== ?()
∞
(3.16)
This could represent a model of a turbine blade with internal cooling. Figure 3.6 shows the
configuration.
HT-26
T
1
T
2
T
w1
T
w2
L
1
δ′
2
δ′
T
Figure 3.6: Wall with convective heat transfer
The heat transfer in fluid 1 is given by
˙
Q
A
hT T
1w 1
1
=?
()
,
which is the heat transfer per unit area to the fluid. The heat transfer in fluid (2) is similarly given by
˙
Q
A
hT T
22 w2
=?().
Across the wall, we have
˙
Q
A
k
L
TT
w2 w1
=?().
The quantity Q
&
/A is the same in all of these expressions. Putting them all together to write the
known overall temperature drop yields a relation between heat transfer and overall temperature drop,
T
2
– T
1
:
TT TT T T T T
A
1
h
L
k
1
h
21 2w2 w2w1 w11
12
?= ?()+?()+?()=++
?
?
?
?
?
?
˙
Q
. (3.17)
We can define a thermal resistance, R, as before, such that
˙
Q=
?()TT
R
21
,
HT-27
where R is given by
R
1
hA
L
Ak
1
hA
12
=++. (3.18)
Equation (3.18) is the thermal resistance for a solid wall with convection heat transfer on each side.
For a turbine blade in a gas turbine engine, cooling is a critical consideration. In terms of Figure 3.6,
T
2
is the combustor exit (turbine inlet) temperature and T
1
is the temperature at the compressor exit.
We wish to find T
w
2
because this is the highest metal temperature. From (3.17), the wall
temperature can be written as
T
w2
= T
2
-
˙
Q
Ah
2
=?
?
T
TT
R
1
Ah
2
21
2
(3.19)
Using the expression for the thermal resistance, the wall temperatures can be expressed in terms of
heat transfer coefficients and wall properties as
1
2
1
2
12
22
++
?
?=
k
Lh
h
h
TT
TT
w
(3.20)
Equation (3.20) provides some basic design guidelines. The goal is to have a low value of T
w
2
.
This means h
1
/h
2
should be large, k should be large (but we may not have much flexibility in choice
of material) and L should be small. One way to achieve the first of these is to have h
2
low (for
example, to flow cooling air out as in Figure 3.1 to shield the surface).
A second example of combined conduction and convection is given by a cylinder exposed to a
flowing fluid. The geometry is shown in Figure 3.7.
c
∞
T
∞
r
1
r
2
r
1
Figure 3.7: Cylinder in a flowing fluid
HT-28
For the cylinder the heat flux at the outer surface is given by
˙
˙
q
Q
== ?()=
∞
A
hT T r r
w 2
at
The boundary condition at the inner surface could either be a heat flux condition or a temperature
specification; we use the latter to simplify the algebra. Thus,
11
at rrTT == . This is a model for the
heat transfer in a pipe of radius r
1
surrounded by insulation of thickness r
2
- r
1
. The solution for a
cylindrical region was given in Section 2.3 as
Tr a
r
r
b() ln=+
1
Use of the boundary condition T (r
1
) = T
1
yields b = T
1
.
At the interface between the cylinder and the fluid, r = r
2
, the temperature and the heat flow are
continuous. (Question: Why is this? How would you argue the point?)
˙q=? =? = +
?
?
?
?
?
?
?
?
?
?
?
?
?
∞
k
dT
dr
k
a
r
h aln
r
r
TT
2
2
1
1
(3.21)
Plugging the form of the temperature distribution in the cylinder into Equation (3.21) yields
?+
?
?
?
?
?
?
=?()
∞
a
k
r
hln
r
r
hT T
2
2
1
1
.
The constant of integration, a, is
a
hT T
k
r
h
r
r
TT
k
hr
r
r
1
2
2
1
1
2
2
1
=
??
()
+
=?
?
()
+
∞∞
ln ln
,
and the expression for the temperature is, in normalized non-dimensional form
TT
TT
r/r
k
hr
r/r
1
1
1
2
21
?
?
=
()
+
()
∞
ln
ln
. (3.22)
heat flux just inside cylinder surface heat flux to fluid
HT-29
The heat flow per unit length,
˙
Q, is given by
˙
Q=
?()
+ ()
∞
2T Tk
k
hr
r/r
1
2
21
π
ln
(3.23)
The units in Equation (3.23) are W / m-s.
A problem of interest is choosing the thickness of insulation to minimize the heat loss for a fixed
temperature difference T
1
- T
∞
between the inside of the pipe and the flowing fluid far away from
the pipe. (T
1
- T
∞
is the driving temperature distribution for the pipe). To understand the behavior of
the heat transfer we examine the denominator in Equation (3.23) as r
2
varies. The thickness of
insulation that gives maximum heat transfer is given by
d
dr
k
hr
r
r
0
22
2
1
+
?
?
?
?
?
?
=ln (3.24)
(Question: How do we know this is a maximum?)
From Equation (3.24), the value of r
2
for maximum Q
&
is thus
(r
2
)
maximum heat transfer
= k/h. (3.25)
If r
2
is less than this, we can add insulation and increase heat loss. To understand why this occurs,
consider Figure 3.8, which shows a schematic of the thermal resistance and the heat transfer. As r
2
increases from a value less than r
2
= k/h, two effects take place. First, the thickness of the insulation
increases, tending to drop the heat transfer because the temperature gradient decreases. Secondly,
the area of the outside surface of the insulation increases, tending to increase the heat transfer. The
second of these is (loosely) associated with the k/hr
2
term, the first with the ln(r
2
/r
1
) term. There are
thus two competing effects which combine to give a maximum
˙
Q at r
2
= k/h.
HT-30
R
Q
?
r
2
ln
r
r
2
1
k
hr
2
k
h
Figure 3.8: Critical radius of insulation
Muddy points
In the expression
1
h.A
, what is A? (MP HT.11)
It seems that we have simplified convection a lot. Is finding the heat transfer coefficient, h,
really difficult? (MP HT.12)
What does the "K" in the contact resistance formula stand for? (MP HT.13)
In the equation for the temperature in a cylinder (3.22), what is "r"? (MP HT.14)
3.3 Dimensionless Numbers and Analysis of Results
Phenomena in fluid flow and heat transfer depend on dimensionless parameters. The Mach number
and the Reynolds number are two you have already seen. These parameters give information as to
the relevant flow regimes of a given solution. Casting equations in dimensionless form helps show
the generality of application to a broad class of situations (rather than just one set of dimensional
parameters). It is generally good practice to use non-dimensional numbers, forms of equations, and
results presentation whenever possible. The results for heat transfer from the cylinder are already in
dimensionless form but we can carry the idea even further. For the cylinder:
T-T
TT
r/r
k/hr r /r
1
1
1
221
'
<
=
()
+
()
ln
ln
(3.26)
The parameter
hr
k
2
or
hL
k
, where L is a relevant length for the particular problem of interest, is
called the Biot number denoted by Bi. In terms of this parameter,
HT-31
T-T
TT
r/r
1
Bi
r/r
1
1
1
21
∞
?
=
()
+
()
ln
ln
(3.27)
The size of the Biot number gives a key to the regimes in which different features are dominant. For
Bi >> 1 the convection heat transfer process offers little resistance to heat transfer. There is thus
only a small ?T outside (i.e. T (r
2
) close to T
∞
) compared to the ?T through the solid with a limiting
behavior of
TT
TT
r/r
r/r
11
21
?
?
=
∞
ln
ln
as Bi goes to infinity. This is much like the situation with an external temperature specified.
For Bi << 1 the conduction heat transfer process offers little resistance to heat transfer. The
temperature difference in the body (i.e. from r
1
to r
2
) is small compared to the external temperature
difference, T
1
- T
∞
. In this situation, the limiting case is
TT
TT
Bi r/r
1
1
1
?
?
=
()
∞
ln
In this regime there is approximately uniform temperature in the cylinder. The size of the Biot
number thus indicates the regimes where the different effects become important.
Figure 3.9 shows the general effect of Biot number on temperature distribution. Figure 3.10 is a plot
of the temperature distribution in the cylinder for values of Bi = 0.1, 1.0 and 10.0.
HT-32
Figure 3.9: Effect of the Biot Number [hL / k
body
] on the temperature distributions in the solid
and in the fluid for convective cooling of a body. Note that k
body
is the thermal conductivity of
the body, not of the fluid.
[adapted from: A Heat Transfer Textbook, John H. Lienhard, Prentice-Hall Publishers, 1980]
HT-33
Figure 3.10: Temperature distribution in a convectively cooled cylinder for different values of
Biot number, Bi; r
2
/ r
1
= 2 [from: A Heat Transfer Textbook, John H. Lienhard]
4.0 Temperature Distributions in the Presence of Heat Sources
There are a number of situations in which there are sources of heat in the domain of interest.
Examples are:
1) Electrical heaters where electrical energy is converted resistively into heat
2) Nuclear power supplies
3) Propellants where chemical energy is the source
These situations can be analyzed by looking at a model problem of a slab with heat sources α (W/m
3
)
distributed throughout. We take the outside walls to be at temperature T
w
. and we will determine the
maximum internal temperature.
T
w
T
w
(a)
heat
sources
3
m
W
α
x
x x + dx
q&
dx
dx
qd
q
&
&+
Slice at x, x+dx
(b)
Figure 4.1: Slab with heat sources (a) overall configuration, (b) elementary slice
With reference to Figure 4.1(b), a steady-state energy balance yields an equation for the heat flux,
q&
.
HT-34
q adx q
dq
dx
dx 0
??
?
+?+
?
?
?
?
?
?
?
?
= (4.1)
or
dq
dx
?
=α. (4.2)
There is a change in heat flux with x due to the presence of the heat sources. The equation for the
temperature is
0/
2
2
=+ k
dx
Td
α (4.3)
Equation (4.3) can be integrated once,
ax
kdx
dT
+?=
α
(4.4)
and again to give
baxx
k
T ++?=
2
2
α
(4.5)
where a and b are constants of integration. The boundary conditions imposed are () ()
w
TLTT ==0 .
Substituting these into Equation (4.5) gives
w
Tb = and
k
L
a
2
α
= . The temperature distribution is thus
W
TLx
kk
x
T ++?=
22
2
αα
. (4.6)
Writing (4.6) in a normalized, non-dimensional fashion gives a form that exhibits in a more useful
manner the way in which the different parameters enter the problem:
TT
L/k
1
2
x
L
x
L
W
2
2
2
?
=?
?
?
?
?
?
?
α
(4.7)
This distribution is sketched in Figure 4.2. It is symmetric about the mid-plane at
2
L
x = , with half
the energy due to the sources exiting the slab on each side.
HT-35
TT
aL /k
W
2
??
?
?
?
?
?
0.125
X
L
0 0.5 1.0
Figure 4.2:Temperature distribution for slab with distributed heat sources
The heat flux at the side of the slab, x = 0, can be found by differentiating the temperature
distribution and evaluating at x = 0 : < = <
£
¤
2
¥
|
′
= <
=
k
dT
dx
k
L
2K
1
L
L/2
x 0
2
_
_ .
This is half of the total heat generated within the slab. The magnitude of the heat flux is the same at
x = L, although the direction is opposite.
A related problem would be one in which there were heat flux (rather than temperature) boundary
conditions at x = 0 and x = L, so that T
w
is not known. We again determine the maximum
temperature. At x = 0 and L, the heat flux and temperature are continuous so
< = <
()
=
'
k
dT
dx
hT T x 0,Lat . (4.8)
Referring to the temperature distribution of Equation (4.6) gives for the two terms in Equation (4.8),
k
dT
dx
k-
x
k
axka=+
£
¤
¥
|
= < +()
_
_ (4.9)
hT T h
x
k
xbT
2
<
()
= < ++ <
£
¤
2
¥
|
′ '