2 Lectures Me315 071 Chapter 3
2 Lectures Me315 071 Chapter 3
Objectives:
1. To determine expressions for the temperature distribution and heat transfer rate in common
(planar, cylindrical, and spherical) geometries.
2. To introduce the concept of thermal resistance and the use of thermal circuits to model heat
flow.
3. To consider the use of extended surfaces (fins) for heat transfer enhancement.
T,1
Ts,1
T,2
x x=L
One-dimensional (1-D) heat conduction occurs in the plane wall. Therefore, the temperature
in the wall is a function of x only and heat transfer is exclusively in the x direction.
T 0 T s,1
(2)
T L T s , 2
Integrate equation (1) twice (assuming k is constant) and apply the boundary conditions to
obtain the temperature distribution as
T x T s, 2 T s,1 Lx T s,1 (3)
Note: Equation (3) shows that for 1-D, steady-state conduction in a plane wall with no heat
generation and constant k, the temperature varies linearly with x.
The conduction heat transfer rate is obtained by using the Fourier’s law:
qx k A
dT k A
dx
L
T s,1 T s, 2 (4)
qx
qx
A
k
L
T s,1 T s, 2 (5)
The thermal resistance for convection is obtained from Newton’s law of cooling
( q h ATs T ) as
T s T 1
R t , conv (8)
q hA
The thermal circuit for heat transfer through the plane wall of section 3.1.1 is shown below.
In terms of the overall temperature difference, the heat transfer rate may also be
expressed as
T,1 T, 2
qx (10)
R tot
1 L 1
where R tot (resistances in series) (11)
h1 A k A h2 A
Surface radiation and convection resistances act in parallel, and if T = Tsur, they may be
combined to obtain a single, effective surface resistance.
T,1
T2
T3
Hot fluid LA LB LC
T, 1, h1 T,4
kA kB kC
Cold fluid
A T, 4, h4
B C
1 LA LC 1
LB
h1 A KA A KC A h4 A
KB A
The 1-D heat transfer rate for the system may be expressed as
T,1 T, 4
qx (13)
Rt
where
LA LB LC
Rt
1 1
R tot (14)
h1 A kAA kB A kC A h4 A
With composite systems it is often convenient to work with an overall heat transfer
coefficient, U, which is analogous to Newton’s law of cooling, i.e.
qx U A T (16)
where ΔT is the overall temperature difference.
L L
R2
K2 A 2
2
T1 T2 T1 T2
1
L
R1
K1 A1
1
where R tot
1 R1 1 R 2
Alternatively, the heat transfer rate can be calculated as the sum of heat transfer rates in
the individual materials, i.e.
T1 T 2 T1 T 2
q x q 1x q 2 x (19)
R1 R2
Composite walls, such as the one shown below, may be characterized by series-parallel
configurations.
Area A
LE LF = LG LH
kF F
T1 kH T2
kE kG
E H
G
Although the heat flow in the composite wall above is multidimensional, approximate
solution can be obtained by assuming one-dimensional heat transfer. Two different thermal
circuits (shown below) may be used:
a)
b)
The actual value of q lies between the values obtained with circuits (a) and (b).
Thermal contact resistance is dependent upon the solid materials, surface roughness, contact
pressure, temperature, and interfacial fluid.
The values of thermal contact resistance for some surface pairs are given in Tables 3.1 and
3.2.
The standard approach for conduction analysis involves solving the heat equation to obtain
the temperature distribution and then applying Fourier’s law to obtain the heat transfer rate.
An alternative conduction analysis can be used to obtain the heat transfer rate if the following
conditions are satisfied: steady-state, 1-D, and no heat generation.
Consider conduction in the system shown below.
Under steady-state conditions with no heat generation and no heat loss from the sides, the
conservation of energy requirements on a differential element gives qx = qx+dx. This
shows that qx is independent of x.
For the general case in which A = A(x) and k = k(T), the integral form of Fourier’s law
may be expressed as
x dx T
qx k T dT
x o Ax
(21)
To
If, at point x = xo, To = T(xo) is known, integration of the above equation provides the
functional form of T(x).
If, at another point x = x1, T1 = T(x) is known, integration of the equation between xo and
x1 provides an expression from which qx may be computed.
For the specific case in which area A is uniform and k is independent of T, the above
equation reduces to
x
qx k T (22)
A
where Δx = x1 - xo and ΔT = T1 – To.
Cylindrical and spherical systems often experience temperature gradients in the radial
direction only and may therefore be treated as one dimensional.
For steady-state conditions with no heat generation, the heat equation for the system is
1 d dT
k r 0
r d r d r
(23)
Integrating the heat equation (assuming constant k) and using the boundary conditions
yield
T s,1 T s, 2 r
T r ln T s, 2
ln r1 r 2 r2
(25)
The heat transfer rate is obtained by using the temperature distribution with Fourier’s
law:
k 2 r L
dT dT
qr k A
dr dr
2 L k T s ,1 T s , 2 (26)
ln r 2 r1
Equation (26) shows that the heat transfer rate qr is a constant in the radial direction.
From equation (26), the thermal resistance for radial conduction in a cylindrical wall is
R t , cond
T ln r 2 r1
(27)
qr 2 L k
Neglecting interfacial contact resistances, the heat transfer rate may be expressed as
(28)
T,1 T, 4
qr
1
ln r 2 r1 ln r 3 r 2 ln r 4 r 3
1
2 r1 L h1 2 k A L 2 k B L 2 k C L 2 r 4 L h 4
If U is defined in terms of the inside area A1, equations (28) and (29) may be equated to
yield
1
U1 (30)
1 r1 r 2 r1 r1 r3
r 4 r1 1
ln ln ln
h1 kA r1 k B r 2 k C r3 r4 h 4
Similar equations could be written for U2, U3, etc. Note that
U 1 A1 U 2 A 2 U 3 A 3 U 4 A 4 R t 1 (31)
Critical radius
For a plane wall exposed to a fluid, an increase in the thickness of the wall results in an
increase in the conduction resistance Rk = L/(kA) but does not change the convection
resistance Rc. Hence, the heat transfer rate will reduce as the wall thickness increases.
However, for geometries with non-constant cross-sectional area (e.g. a cylinder), increase in
the wall thickness does not always bring about a decrease in the heat transfer rate.
The critical radius of insulation for a cylinder exposed to convection is
k
rcr
h
where k is the of thermal conductivity of the insulation material and h is the convection heat
transfer coefficient on the insulation. This is illustrated in example 3.5.
k 4 r 2
dT dT
q r k A (32)
dr dr
Conservation of energy requires that qr = qr+dr for steady-state, 1-D conduction with no
heat generation. Therefore, qr is a constant independent of r. Equation (32) may therefore
be expressed in integral form as
qr r2 d r Ts , 2
k T dT
4 r
(33)
1 r
2 Ts , 1
Spherical composites may be treated the same way as composite walls and cylinders.
The critical radius of insulation for a sphere is
2k
rcr
h
See Table 3.3 for a summary of the steady-state solutions to the heat equation with no
generation for plane, cylindrical and spherical walls.
Thermal energy generation within a medium is a result of conversion from some other energy
forms.
A common source of energy generation is electrical energy due to a current I passing through
a medium of electrical resistance Re.
The rate of energy generation is
E g I 2 R e (36)
Other sources of energy conversion are nuclear, chemical, and electromagnetic energy.
For constant thermal conductivity k, the heat equation for the system is
d 2T q
0 (38)
dx2 k
Integrating the heat equation twice and using the boundary conditions yield the
temperature distribution
q L 2 x 2 T s, 2 T s,1 x T s,1 T s, 2
T x 1 (40)
2k L2
2 L 2
When both surfaces of the plane wall are maintained at the same temperature Ts,1 = Ts,2 = Ts,
the temperature distribution is symmetrical about the midplane.
q L 2 x 2
T x 1 Ts (41)
2k L2
The maximum temperature exists at the midplane
q L 2
T 0 T o Ts (42)
2k
T x To
2
x
(43)
T s To L
k
dT
d x xL
h T s T (44)
At steady-state, the rate of heat generation within the cylinder must equal the rate at
which heat is convected from the surface of the cylinder to a moving fluid. This condition
allows the surface temperature to be maintained at a fixed value of Ts.
For constant thermal conductivity k, the heat equation for the system is
1 d dT q
r 0
r d r d r k
(45)
T r o T s
dT
0 and (46)
d r r 0
Integrating the heat equation twice and using the boundary conditions yield the
temperature distribution
q r o r 2
2
T r 1 Ts (47)
4k r 2
o
2
q r o L h 2 r o L T s T (50)
The term extended surface is commonly used to depict an important special case involving
heat transfer by conduction within a solid and heat transfer by convection (and/or radiation)
from the boundaries of the solid.
A strut that connects two walls at different temperatures as shown in Figure 3.11 is an
example of extended surface.
The most frequent application of conduction-convection effects is one in which an extended
surface is used specifically to enhance heat transfer between a solid and an adjoining fluid.
Such an extended surface is termed a fin. [See Figure 3.12.]
Examples of fin applications can be found in the engine heads on motorcycles and lawn
mowers, electric power transformers and the condenser of an airconditioner.
There are different fin configurations (see figure 3.14).
d 2T 1 d Ac d T 1 h d As
T T 0
2 (52)
A d x d x A k d x
dx c c
where Ac is the cross-sectional area and dAs is the surface area of the differential element.
For those fins, Ac is a constant and As = Px, where As is the surface area measured from
the base and P is the fin perimeter. Therefore,
dAc/dx = 0 and dAs/dx = P (53)
Substitute (53) into (52) to obtain
d 2T
hP
T T 0 (54)
dx 2 k A c
To simplify the form of the equation, let us define a new variable called excess
temperature θ as
θ (x) = Tb - T (55)
Substitute (55) into (54) to obtain
(56)
d 2
m 2 0
dx 2
where
hP
m2 (57)
k Ac
Equation (56) is the differential equation for the variation of excess temperature from the
base to the tip of the fin. It is a linear, homogenous, 2nd- order differential equation with
constant coefficients. Its general solution is of the form
x C 1 e m x C 2 e m x (58)
Two boundary conditions are needed to evaluate the constants C1 and C2 in equation (58).
The first boundary condition is that of the temperature at the base of the fin (x = 0), i.e.
θ(0) = Tb - T = θb
The second boundary condition is specified at the fin tip (x = L). Four different possible
cases at the tip are considered.
Case A: convection heat transfer from the fin tip
Case B: adiabatic fin tip
Case C: prescribed temperature at the fin tip
Case D: very long fin.
The fin tip boundary condition corresponding to each case is presented in Table 3.4. The
table also contains the resulting temperature distribution and heat transfer rate.
Note: The heat transfer rate is obtained by applying either the Fourier's law at the fin base
or Newton's law of cooling over the fin surface.
The fin effectiveness f is the ratio of the fin heat transfer rate to the heat transfer rate that
would exist without the fin.
qf
f (59)
h A c, b b
The thermal resistance due to convection at the exposed base (if there were no fin) is
1
R t, b (62)
h A c, b
Note: The usefulness of fin can also be assessed with the use of a parameter called the Biot number, Bi,
hl
Bi
k
where l is the ratio of the fin surface area to the fin perimeter, l = A/P.
As an approximation, the Biot number can be taken as the ratio between the conduction resistance
and the convection resistance
R t , cond
Bi
R t , conv
For a fin to be considered effective, the following condition is to be satisfied
Bi 0.1
Efficiency
Another measure of fin thermal performance is the fin efficiency f, defined as
qf qf
f (64)
q max h A f b
where Af is the surface area of the fin and qmax is the heat transfer rate if the entire fin surface
were at the base temperature.
For example, for a straight fin of uniform cross-section and an adiabatic tip, the fin efficiency
is
M tanh mL tanh mL
f (65)
h P L b mL
Corrected length
In lieu of the cumbersome expression for heat transfer from a straight fin with an active tip
(case A of Table 3.4), approximate prediction may be obtained by using the adiabatic tip
result (case B of Table 3.4) with a corrected fin length of the form
Lc = L + (t/2) for a rectangular fin
Lc = L + (D/4) for a pin fin
Hence, for fin with convection at the tip, the heat transfer rate and efficiency may be
approximated as
q f M tanh m Lc (66)
tanh mL c
f (67)
m Lc
qf
f (68)
h A f b
where At is the sum of the surface areas of the fins and the exposed portion of the base
(often termed the prime surface). qt is the total heat transfer rate from the surface area At.
If there are N fins in the array, each of surface area Af and the area of the prime surface is
designated as Ab, the total surface area is
At N A f Ab (70)
The total rate of heat transfer by convection is
q t N f h Af b h Ab b (71)
o 1
NAf
At
1 f (73)
Once o is known, the total heat transfer rate for the fin array may be calculated using
equation (69).
From (69), an expression may be inferred for the thermal resistance of a fin array
b 1
R t, o (74)
qt o h At
where Rt,o is an effective resistance that accounts for parallel heat flow paths by
conduction/convection in the fins and by convection from the prime surface.
The figures below illustrate the thermal circuits corresponding to the parallel paths and their
representation in terms of an effective resistance. Two possible cases are considered:
a) Fins that are integral with the base
b) Fins that are attached to the base (in which case the thermal contact resistance between fin
and base has to be accounted for)