Introduction To Heat Transfer
Introduction To Heat Transfer
1 Introduction to
Heat Transfer
Questions-Answers
Que 1.1. What are the mechanisms of heat transfer ? How are
Answer
A. Mechanisms of Heat Transfer : Following are the three modes of
heat transfer :
a. Mechanism of Heat Transfer through Conduction :
1. Thermal conduction is a mechanism of heat propagation from region of
higher temperature to a region of lower temperature within a medium
or between different mediums in direct physical contact.
Heat & Mass Transfer 1–3 B (ME-Sem-5)
Answer
A. Thermal Conductivity :
1. The amount of heat conducted through a body of unit area, and unit
thickness in unit time when the difference in temperature between the
faces causing heat flow is unity.
2. The unit of thermal conductivity is W/mK or W/m °C.
Heat & Mass Transfer 1–5 B (ME-Sem-5)
Answer
1. In actual conditions, more than one mode of heat transfer is involved.
2. For example, consider a plate of emissivity maintained at temperature
Ts. The heat transfer from this plate will take place through both
convection and radiation.
3. The heat loss from this plate is
q = qconv + qrad
q = hA(Ts – T) + A (Ts4 – T4 )
T
qconv qrad
Ts
Fig. 1.3.1.
Introduction to Heat Transfer 1–6 B (ME-Sem-5)
Que 1.4. A brick (k = 1.2 W/m-K) wall 0.15 m thick separates hot
combustion gases of a furnace from the outside ambient air which
is at 25 °C. The outer surface temperature of the brick wall is found
to be 100 °C. If the natural convection heat transfer coefficient on
the outside of the brick wall is 20 W/m2-K and its emissivity is 0.8,
calculate the inner surface temperature of the brick wall.
AKTU 2013-14, Marks 10
Answer
Answer
(ts ti )
– 45 × 0.54 × = 22 × 0.54 × (310 – 15) + 250
L
(310 ti )
– 45 × 0.54 × = 22 × 0.54 × 295 + 250
0.025
972 (ti – 310) = 3754.6
3754.6
ti = 310 = 313.86 °C
972
Qrad.
tf = 15 °C
Air
(fluid)
Qconv. = hA(ts - tf) ts = 310 °C
L = 25 mm
Answer
1. Consider an infinitesimal rectangular parallelopiped (volume element)
of sides dx, dy, and dz parallel, respectively, to the three axes (X, Y, Z) in
a medium in which temperature is varying with location and time as
shown in Fig. 1.6.1.
2. Let, t = Temperature at the left face ABCD; this
temperature may be assumed uniform over the
entire surface, since the area of this face can be
made arbitrarily small, and
dt
= Rate of change of temperature along X-direction.
dx
t
Then, dx = Change of temperature through distance dx,
x
and
t
t dx = Temperature on the right face EFGH (at a
x
distance dx from the left face ABCD).
Further, let, kx, ky, kz = Thermal conductivities along X, Y and Z axes.
Introduction to Heat Transfer 1–8 B (ME-Sem-5)
Q(y + dy ) Qz
A (x , y , z )
D H
y
O X
z
x C G dy
Z Q(x + dx )
Elemental Qx A (x , y , z )
volume E
(rectangular
parallelopiped)
dz
Q(z + dz) B F
dx
Qg = qgdxdydx
Qy
Fig. 1.6.1. Elemental volume for three-dimensional heat
conduction analysis - Cartesian coordinates.
3. Quantity of heat flowing into element from the left face ABCD during
the time interval d in X-direction is given by
t
Heat influx, Qx = – kx (dydz)
d ...(1.6.1)
x
4. During the same time interval d the heat flowing out of the right face
of control volume (EFGH) will be
Heat efflux, Q(x + dx) = Qx
(Qx )dx ...(1.6.2)
x
5. Heat accumulation in the element due to heat flow in X-direction,
subtracting eq. (1.6.2) from eq. (1.6.1), we get
dQx = Qx Qx (Qx )dx
x
= (Qx ) dx
x
t
= kx ( dydz) d dx
x x
t
= kx dxdydzd
x x
6. Similarly the heat accumulated due to heat flow by conduction along Y
and Z-directions in time d will be :
dQy = ky t dxdydzd
y y
Heat & Mass Transfer 1–9 B (ME-Sem-5)
t
dQz = kz dxdydzd
z z
7. So, net heat accumulated in the element due to conduction of heat from
all the coordinate directions is given as,
t t
= kx dxdydzd + ky dxdydzd
x x y y
t
+ kz dxdydzd
z z
t t t
= kx ky kz dxdydzd
x x y y z z
...(1.6.3)
8. The total heat generated in the element is given by
Qg = qg(dxdydz)d ...(1.6.4)
9. Energy stored in element,
t
= ( dxdydz)c
d ...(1.6.5)
10. Now using energy balance for the element, we have
t t t
x kx x y ky y z kz z dxdydzd + qg(dxdydz)d
t
= ( dxdydz)c d
11. Dividing both side by dxdydzd, we have
t t t t
kx ky kz + qg = c
x x y y z z
This is known as the general heat conduction equation for ‘non-
homogeneous material’, ‘self heat generating’ and ‘unsteady three-
dimensional heat flow’.
12. In case of homogeneous and isotropic material (kx = ky = kz= k),
2 t 2 t 2 t qg c t
=
x 2 z2 z2 k k
2 t 2 t 2 t qg 1 t
or =
x 2 y 2 z2 k
k
where, = = Thermal diffusivity
.c
Introduction to Heat Transfer 1–10 B (ME-Sem-5)
Answer
1. Consider an elemental volume having the coordinate (r, , z) for three
dimensional heat conduction analysis as shown in Fig. 1.7.1.
2. The volume of the element,
dv = r d dr dz.
Let qg = Heat generation (uniform) per unit volume
per unit time.
k = Thermal conductivity.
= Density.
c = Specific heat capacity.
Q(z + dz)
d z
rd
Z
Q ( + d )
(r , , z ) Qr (r – )
A plane
dr
dz (z – ) Elemental
plane
Q volume
z Y (r – z)
A (r , , z )
r plane dz
Q(r + dr )
dr r
X
Qz
Fig. 1.7.1. Elemental volume for three-dimensional heat
conduction analysis-cylindrical coordinates.
2 t 1 t
dQr = k(dr rd dz) 2 d
r r r
4. Heat flow in tangential direction i.e., (r-z) plane
t
Heat influx, Q = – k(dr dz) d
r
Heat efflux, Q( + d) = Q + (Q) rd
r
Heat stored in the element,
dQ = Q – Q( + d)
t
=– (Q) rd = – k(dr dz) r d rd
r r
1 2t
dQ = k(dr rd dz) d
r 2 2
5. Heat flow in axial direction i.e., (r-) plane
t
Heat influx, Qz = – k(rd dr) d
z
Heat efflux, Q(z + dz) = Qz + (Qz) dz
z
Heat stored in the element,
dQz = Qz – Q(z + dz)
t
=– k( rd dr ) z d dz
z
2 t
dQz = k(dr rddz) d
z2
6. Net heat stored in the element
2 t 1 t 1 2 t 2 t
= k. dr rddz 2 · 2 2 2 d
r r r r z
2 t 1 t 1 2 t 2t q g c t 1 t
or = = ...(1.7.1)
r 2 r r r 2 2 z2 k k
This is the general heat conduction equation in cylindrical coordinates.
1. For steady-state, unidirectional heat flow in the radial direction, and
with no internal heat generation, eq. (1.7.1) reduces to
2 t 1 t
2 =0
r r r
1 d dt
or r =0
r dr dr
1
Since 0
r
d dt dt
r = 0 or r = constant
dr dr dr
Que 1.8. Derive the general heat conduction equation in
spherical coordinates.
Answer
1. Consider an elemental volume having the coordinates (r, , ), for three
dimensional heat conduction analysis, as shown in Fig. 1.8.1.
2. The volume of the element = dr rdrsin d
3. Let, qg = Heat generation (uniform) per unit volume per unit time.
We assume that thermal conductivity (k), density () and specific heat
(c) are uniform.
4. Heat flow through r- plane (-direction) :
t
Heat influx, Q = – k (dr rd) d ...(1.8.1)
r sin
Heat efflux, Q( + d) = Q + (Q ) r sin d ...(1.8.2)
r sin
Heat accumulated in the element due to heat flow in the -direction,
dQ = Q – Q( + d)
1
=– (Q ) r sin d
r sin
1 1 t
=– k( dr rd) d r sin d
r sin r sin
1 2t
= k (dr rd r sin d) d
r sin 2
2 2
Heat efflux, Q( + d) = Q +
(Q ) rd ...(1.8.4)
r
Heat accumulated in the element due to heat flow in the
-direction,
dQ = Q – Q(+d) [Subtracting eq. (1.8.4) and (1.8.3)]
=– (Q) rd
r
t
=– k( dr r sin d) d r d
r r
k dr rd rd t
= sin d
r r
1 t
= k (dr rd r sin d) sin d
r 2 sin
6. Heat flow in - plane (r-direction) :
t
Heat influx, Qr = – k (rd r sin d) ...(1.8.5)
r
Heat efflux, Q(r + dr) = Qr +
(Qr)dr ...(1.8.6)
r
Heat accumulation in the element due to heat flow in the
r-direction,
dQr = Qr – Q(r + dr) = – (Qr ) dr
r
t
=– k (rd r sin d) d dr
r r
2 t
= k d sin d dr r d
r r
1 2 t
= k (dr rd r sin d) 2 r d
r r r
Z Q( + d)
A( r, , ) r sin d
Qr
r Q( + d)
A (r , , ) Elemental
rd
Y volume
Q Q(r + dr )
X dr
Q
Fig. 1.8.1. Elemental volume for three-dimensional heat
conduction analysis-spherical coordinates.
7. Net heat accumulated in the element
Introduction to Heat Transfer 1–14 B (ME-Sem-5)
1 2t 1 t 1 2 t
= k dr rd r sin d 2 sin 2 r d
r 2
sin 2
2
r sin r r r
8. The total heat generated within the element is given by,
Qg = qg (dr rd r sin d) d
9. The increase in thermal energy in the element
t
d
= (dr rd r sin d) c
10. From the principle of conservation of energy,
1 2t 1 t 1 2 t
k dr rd r sin d 2 2 sin 2 r d
r sin 2
r sin r r r
t
+ qg (dr rd r sin d) d = (dr rd r sin d) c d
11. Dividing both sides by k(dr rd r sin d) d, we get
1 2t 1 t 1 2 t qg c t
2 2 sin 2 r =
2
r sin
2
r sin r r r k k
2
1 t 1 t 1 2 t q g 1 t
or 2 2 2
2 sin 2 r =
r sin r sin r r r k d
...(1.8.7)
This equation is known as general heat conduction equation in spherical
coordinates.
12. In case there are no heat sources present and the heat flow is steady
and one-dimensional, then eq. (1.8.7) reduces to
1 d 2 dt
r =0
r 2 dr dr
PART-2
Initial and Boundary Conditions, Simple and Composite
Systems in Rectangular, Cylindrical and Spherical Coordinates
With and Without Energy Generation.
Questions-Answers
Answer
A. Initial Conditions :
1. The initial conditions describe the temperature distribution in a medium
at the initial moment of time, and these are needed only for the transient
(time-dependent) problems.
2. The initial conditions can be expressed as
At = 0; t = t (x, y, z) ...(1.9.1)
3. For a uniform initial temperature distribution, a simple but typical form
of the above identity can be recast as
At = 0; t = t0 (x, y, z) ...(1.9.2)
B. Boundary Conditions :
1. The boundary conditions refer to physical conditions existing at the
boundaries of the medium, and specify the temperature or the heat flow
at the surface of the body.
C. Types of Boundary Conditions :
The different types of boundary conditions applied to heat conduction
problems are as follows :
a. Boundary Condition of First K ind (Prescribed Surface
Temperature) :
1. The temperature distribution, ts is given at a boundary surface for each
moment of time.
ts = t(x, y, z, )
2. A typical example of boundary condition of the first kind for a slab is
shown in Fig. 1.9.1.
t(x, y, ) = 0, at x = 0
t(x, y, ) = 0, at y = 0
t(x, y, ) = f1(y), at x = a
t(x, y, ) = f2(x), at y = b
Y
ty = b = f2(x)
b
t x=0 = 0 t(x, y, ) tx =a = f1(y)
O a X
ty = 0 = 0
Fig. 1.9.1. Boundary condition of the first kind.
Introduction to Heat Transfer 1–16 B (ME-Sem-5)
t
= F0
x x=0 t(x, )
or X
– k t
= q0
x x=0
3. Similarly at x = L,
t
– k = h2 (tx = L – t2)
x x L
t
or k h2 t = h2t2 = F2
x xL
t1
tx =L
t
qx = h1(t 1 – tx = 0) = – k tx t
x x = 0 =0 qx = h2 (tx = L – t2) = – k
x x = L
t2
Answer
1. Consider a plane wall of homogeneous and isotropic material through
which heat is flowing only in X-direction. The general heat conduction
equation is,
2t 2t 2t qg 1 t
2
+ 2
+ 2
+ = ...(1.10.1)
x y z k
t
2. If the heat conduction takes place under steady state 0 , one-
2 t 2t
dimensional 2 2 0 and with no internal heat generation
y z
qg
0 . Then, the eq. (1.10.1) reduced to,
k
2t
=0 ...(1.10.2)
x 2
T
t1 Plane wall
k
Q dt Q
x Q t1 t2 Q
dx
t2 L
(Rth)cond. =
1 kA
L X
2
(a) (b )
Fig. 1.10.1. Heat conduction through a plane wall.
3. By integrating the eq. (1.10.2) twice, we get
dt
= C1 and t = C1 x + C2 ...(1.10.3)
dx
Where C1 and C2 are arbitrary constants.
4. By applying boundary conditions in eq. (1.10.3),
At, x = 0, t = t1
At, x = L, t = t2
t t
We get, t1 = C2 and C1 = 2 1
L
Introduction to Heat Transfer 1–18 B (ME-Sem-5)
Answer
1. Consider the transmission of heat through a composite wall consisting
of a number of slabs as shown in Fig. 1.11.1.
2. Let LA, LB, LC = Thicknesses of slabs A, B and C respectively
(also called path length),
kA, kB, kC = Thermal conductivities of the slabs A, B, and C
respectively,
t1, t4 (t1 > t4) = Temperatures at the wall surfaces 1 and 4
respectively, and
t2, t3 = Temperatures at the interfaces 2 and 3
respectively.
3. Since the quantity of heat transmitted per unit time through each slab/
layer is same, so we have
kA A(t1 t2 ) kB A(t2 t3 ) kC A(t3 t4 )
Q=
LA LB LC
(Assuming that there is a perfect contact between the layers and no
temperature drop occurs across the interface between the materials).
4. Rearranging the above expression, we get
QLA
t1 – t2 = ...(1.11.1)
kA A
QLB
t2 – t3 = ...(1.11.2)
kB A
QLC
t3 – t4 = ...(1.11.3)
kC A
Heat & Mass Transfer 1–19 B (ME-Sem-5)
t1 Interfaces
t2
Q t3 Q
Temperature
profile
A B C t4
kA kB kC
1 2 3 4
LA LB LC
(a)
Q t1 t2 t3 t4 Q
LA LB LC
Rth – A = , Rth – B = , Rth – C =
k AA kBA kC A
( b)
Fig. 1.11.1. Steady state conduction through a composite wall.
A(t1 t4 ) ( t1 t4 ) (t1 t4 )
or, Q = =
LA LB LC LA LB LC [ Rth A Rth B Rth C ]
k k k k A k A k A
A B c A B C
6. If the composite wall consists of n slab/layers, then
[t t ]
Q = 1 n ( n1)
L
1 kA
Que 1.12. Derive the expression for steady state heat conduction
through a hollow cylinder with uniform thermal conductivity,
whose inner surface is exposed to hot fluid and outer surface is
exposed to a cold fluid.
Introduction to Heat Transfer 1–20 B (ME-Sem-5)
Answer
1. Consider a hollow cylinder made of material having constant thermal
conductivity and insulated at both ends.
2. The general heat conduction equation in cylindrical co-ordinate is,
2 t 1 t 1 2t 2t qg 1 t
2 + + 2 2
+ + = ...(1.12.1)
dr r r r z2 k
t
3. For steady state 0 , unidirectional (t f [, x]) heat flow with no
internal heat generation (qg = 0), then the eq. (1.12.1) becomes
d 2 t 1 dt
=0
dr 2 r dr
1 d dt
or r =0
r dr dr
1
Since 0,
r
d dt
r =0 ...(1.12.2)
dr dr
(Cold fluid (t2)
Hollow cylinder
(length = l ) Element
No heat flows in the axial
Hot
r2 direction
fluid (t1)
r
r1
t1 dt
t2
Q t1 t2 dr Q
r
Rth = 1 ln 2
2kL r1
Fig. 1.12.1.
At, r= r1, t = t1
At, r= r2, t = t2
We get, t1 = C ln r1 + C1 ...(1.12.4)
and t2 = C ln r2 + C1 ...(1.12.5)
6. From eq. (1.12.4) and eq. (1.12.5), we get
t1 t2
C=–
r
ln 2
r
1
t1 t2
and C1 = t1 + ln r1
ln (r2 / r1 )
7. Substituting the value of C and C1 in eq. (1.12.3) we get
(t1 t2 ) (t t )
t = t1 + ln (r1) – 1 2 ln (r)
r2 r
ln ln 2
r1 r1
r
ln
t t1 r1
or = ...(1.12.6)
t2 t1 r
ln 2
r1
This gives the temperature distribution in a hollow cylinder which is
varying logarithmically and is a function of radius only.
8. Now, the heat transfer rate is given as,
r
( t2 t1 ) ln
dt d
t1 r1
Q = – kA = – kA
dr dr r
ln 2
r1
[From eq. (1.12.6)]
(t1 t2 ) (t1 t2 ) t
= 2 kL = =
r2 ln (r2 / r1 ) Rth
ln
r1 2 kL
ln (r2 / r1 )
Rth 2kL
t1 t2
So, Q=
ln (r2 / r1 )
2kL
This is the expression for rate of heat transfer through hollow cylinder
which is insulated at both ends.
Answer
1. Consider heat flows radially through a composite cylinder as shown in
Fig. 1.13.1.
2. Let, th and tc = The temperature of hot fluid and cold fluid
f f
respectively,
hhf and hcf = Convective heat transfer coefficient of hot fluid
and cold fluid respectively,
kA and kB = Thermal conductivity of inside layer A and of
outside layer B, and
t1, t2, t3 = Temperatures at point 1, 2 and 3.
3. The rate of heat transfer,
kA (2 L)( t1 t2 )
Q = hh 2r1 L (th – t1) =
f f r
ln 2
r1
kB 2 L (t2 t3 )
= = hc 2r3 L (t3– tc )
r f f
ln 3
r 2
4. By rearranging, we get
Q
th – t1 = ...(1.13.1)
f hhf r1 2 L
Cold fluid
(air) tcf
Q
hcf
B
Hot thf A
fluid
hhf t1 t2 t3
r
r2 1
r3
Fig. 1.13.1.
Q
kA 2L
t1 – t2 = ...(1.13.2)
r
ln 2
r1
Heat & Mass Transfer 1–23 B (ME-Sem-5)
Q
kB 2L
t2 – t3 = ...(1.13.3)
r
ln 3
r2
Q
t3 – tcf = ...(1.13.4)
hcf r3 2L
5. Adding eq. (1.13.1), eq. (1.13.2), eq. (1.13.3) and eq. (1.13.4), we get
1 1 1 1
h r k
k
h
Q hf 1 A B c f r3
= th – tc
2L r r f f
ln 2 ln 3
r
1 r
2
2L (thf tcf )
or Q=
1 1 1 1
h r k
k
h
A B c f r3
hf 1
r r
ln 2 ln 3
r r
1 2
2L (thf tcf )
or Q=
r r
ln 2 ln 3
1 r1 r2 1
hh r1 kA kB hcf r3
f
6. If there are n concentric cylinders, then
2L (thf tcf )
Q= ...(1.13.5)
1 n
1 rn 1 1
ln
hhf r1 n1 kn rn hcf rn 1
7. If inside and outside heat transfer coefficient are not considered, then
eq. (1.13.5) can be written as,
2 L(t t( n1) )
Q= n
1 r
ln n 1
k
n 1 n rn
Que 1.14. Derive expressions under one dimensional steady state
heat conduction for temperature distribution for the sphere.
OR
Starting with an energy balance on a spherical shell volume
element, derive the one dimensional transient heat conduction
equation for a sphere with constant thermal conductivity and no
heat generation [Fig. (1.14.1)].
Introduction to Heat Transfer 1–24 B (ME-Sem-5)
r + r
r
O R r
Fig. 1.14.1.
Answer
1. Let thermal conductivity of a hollow sphere is uniform.
2. Let r1 and r2 are inner and outer radius and t1 and t2 are temperatures
of inner surface and outer surface and k be thermal conductivity of the
material with the given temperature range.
Q, Heat flows radially
outwards ( t1 > t2)
Hollow sphere
dr
Element
r2
r
r1 t1 t2
Q t1 t2 Q
(r2 – r1)
Rth =
4kr1r2
Fig. 1.14.2. Steady state conduction through a hollow sphere.
3. Applying the general heat conduction equation in spherical co-ordinates,
1 2 t 1 2t 1 t q g 1 t
r 2 sin =
r r r r sin 2 r 2 sin
2 2
k
...(1.14.1)
t
4. For steady state
0 , unidirectional heat flow in the radial
direction and with no heat generation, eq. (1.14.1) reduces to
1 2 dt
r =0
r 2 r dr
Heat & Mass Transfer 1–25 B (ME-Sem-5)
d 2 dt 1
r =0 2 0
dr dr r
dt
or r2 = Constant = C ...(1.14.2)
dr
5. On integrating eq. (1.14.2) twice, we have
C
t= + C1 ...(1.14.3)
r
6. Applying boundary condition,
At r = r1, t = t1
At r = r2, t = t2
C
We get, t1 = + C1 ...(1.14.4)
r1
C
and t2 = + C1 ...(1.14.5)
r2
7. From eq. (1.14.4) and (1.14.5), we get
(t1 t2 ) r1 r2
C=
r1 r2
(t1 t2 ) r1 r2
and C1 = t1 +
r1 (r1 r2 )
8. Substituting the values of C and C1 in eq. (1.14.3), we have
(t1 t2 ) r1 r2 (t t ) r r
t=– + t1 + 1 2 1 2
r (r1 r2 ) r1 (r1 r2 )
(t1 t2 ) (t1 t2 )
t=– + t1 +
1 1 1 1
r r1
r2 r1 r2 r1
t1 t2 1 1
t = t1 + ...(1.14.6)
1 1 r1 r
r r
2 1
1 1
t t1 r r1
or =
t2 t1 1 1
r2 r1
t t1 r r r1
or = 2 ...(1.14.7)
t2 t1 r r2 r1
Eq. (1.14.7) shows that temperature distribution associated with radial
conduction through a sphere is represented by a hyperbola.
9. Now, from Fourier’s equation,
Introduction to Heat Transfer 1–26 B (ME-Sem-5)
Q = kA dt k4 r 2 d t1 t1 t2 1 1
dr dr 1 1 r1
r
r2 r1
t t 1 t t 1 ( t t )rr
= k4 r2 1 2 2 = k4 r 2 1 2 2 = 4 k 1 2 1 2
1 1 r r1 r2 r (r1 r2 )
r2 r1 r1r2
t1 t2 t t
Q= = 1 2
(r2 r1 ) Rth
4 kr1r2
r2 r1
Where, Rth =
4 kr1 r2
Answer
ho t4
to
Air film
t3
k3
t2
k2
Physical system
k1 t1 St
r2 fil eam
m
r1
0.5
8
cm
cm
3
cm
5
cm
Thermal circuit
ti t1 t2 t3 t4 to
q q
Ri R1 R2 R3 Ro
Fig. 1.15.1.
ti t4 300 50
2. Heat loss per unit length of the pipe, q =
R 3.014
q = 82.95 W/m
Que 1.16. Consider a 0.8 m high and 1.5 m wide double pane window
consisting of two 4 mm thick layers of glass (k = 0.78 W/m-°C) separated
by a 10 mm wide stagnant air space (k = 0.026 W/m-°C). Determine
the steady rate of heat transfer through this double pane window
and the temperature of its inner surface for a day during which the
room is maintained at 20 °C while the temperature of the outdoors
is – 10 °C. Take the convection heat transfer coefficients on the
inner and outer surface of the window to be h1 = 10 W/m2-°C and
h2 = 40 W/m2-°C, which includes the effects of radiation.
AKTU 2014-15, Marks 10
Answer
ta t
Q=
1 x1 x x 1
2 3
ha A k1 A k2 A k3 A h A
20 (10)
=
1 0.004 0.01 0.004 1
10 1.2 0.78 1.2 0.026 1.2 0.78 1.2 40 1.2
= 69.24 W
2. Also, Q = hiAs(ti – t1)
69.24 = 10 × 1.2 × (20 – t1)
t1 = 14.23 °C
Answer
40 °C
k2 k1 k2
20 °C
12 cm
Fig. 1.17.1.
Heat & Mass Transfer 1–29 B (ME-Sem-5)
1.2 cm
40 °C 1.2 cm
k2 k2
k1
20 °C
1.2 cm
12 cm (a) (b)
Fig. 1.17.3.
1. Thermal circuit for the given situation is shown below :
t1= 40 °C t4= 20 °C
R2 R1 R2
Rb
Rb
Rb
Rb
Fig. 1.17.4. Thermal circuit for the system.
Introduction to Heat Transfer 1–30 B (ME-Sem-5)
Lb 0.18
2. Now, Rb = 2
= 39.78
kb Ab 1.2
40
4 100
1 2 1 2
= {Rs = R1 + R2 + R2}
Reffective Rb Rs Rb
4 1 4 Rs Rb 4 0.945 39.78
= = = 1.16
Rb Rs Rb Rs 0.945 39.78
So, Reff = 1/1.16 = 0.86
t 40 – 20
3. Heat flux, q2 = = 23.2 W/m2
Reff 0.86
4. Change in flux = q2 – q1 = 23.2 – 21.164 = 2.036 W/m2
Answer
1 1 ln (r2 / r1 ) ln (r3 / r2 ) 1
Rth =
2L
hi ri k 1 k 2 h0 r0
2.75 cm
ti
to
2.5 cm
3.25 cm
Fig. 1.18.1.
t 200 25
3. We know that, Q = =
Rth 0.3216
L
Q
= 544.154 W/m ...(1.18.1)
L
4. Let temperature at interface is t,
t
Q=
Rth
200 t
Q=
1 ln (r2 / r1 )
(hiri )(2L) 2k1 L
200 t
=
1.369 10 3 0.9279 10 3
L L
Q 200 t
= ...(1.18.2)
L 2.2969 10 3
5. Equating eq. (1.18.1) and eq. (1.18.2), we get
200 t
544.154 =
2.2969 10 3
t = 198.75 °C.
6. Overall heat transfer coefficient,
Q = UAt
Q Q 544.154 L
U=
At 2r3 Lt 6
2 L 175
100
= 8.248 W/m2 -K
Introduction to Heat Transfer 1–32 B (ME-Sem-5)
Answer
1500 °C
(t1) 90 °C
(t5)
220 150 3 mm
mm mm 50
mm
Fig. 1.19.1.
Answer
1. Consider a plane wall of thickness L, of uniform thermal conductivity
k and in which heat sources are uniformly distributed in the whole
volume.
2. Let the wall surface are maintained at temperature t1 and t2.
3. Assume that heat flow is one-dimensional, under steady state conditions
and there is a uniform volumetric heat generation within the wall.
4. Consider an element of thickness dx at a distance x from the left face
of the wall.
5. Heat conducted in at a distance x,
dt
Qx = – kA
dx
Heat generated in the element,
Qg = qg Adx
Heat conducted out at a distance (x + dx),
d
Q(x + dx) = Qx +(Qx) dx
dx
6. From energy balance, we get,
Qx + Qg = Qx + dx
d
Qx + Qg = Qx + (Q ) dx
dx x
d
Qg = (Q ) dx
dx x
Introduction to Heat Transfer 1–34 B (ME-Sem-5)
2
d dt dt
qg A dx = kA dx = – kA dx
dx dx dx2
d 2t q
or + g =0 ...(1.20.1)
dx 2 k
7. Integrating eq. (1.20.1) we get
dt q
= – g x + C1 ...(1.20.2)
dx k
8. On integrating eq. (1.20.2), we get.
qg 2
t=– x + C1 x + C2 ...(1.20.3)
2k
9. Now, if both the surfaces are at same temperature then,
at, x = 0, t = t1 = tw
and at, x = L, t = t2 = tw
T
Element
Plane
tmax
wall
tw tw
Qx + dx
Qx
Qg
x=0 x=L
x dx x = L /2
L
Fig. 1.21.1.
10. By applying above boundary conditions we in eq. (1.20.3), get,
C2 = tw
q
and C1 = g L
2k
11. Substituting the value of C1 and C2 in eq. (1.20.3) :
qg qg
t=– x2 +
L x + tw
2k 2k
12. For the location of maximum temperature :
dt dt q
= 0 or = g (L – 2x) = 0
dx dx 2k
Heat & Mass Transfer 1–35 B (ME-Sem-5)
qg L
Since, 0, therefore, (L – 2x) = 0 or x =
2k 2
L
13. The maximum temperature occurs at x = i.e., mid plane of plane
2
wall and its value equals
qg
tmax = ( L x) x + tw
2 k x L
2
q L L
tmax = g L + tw
2 k 2 2
qg 2
tmax = L + tw ...(1.20.4)
8k
14. Heat transfer takes place towards both the surfaces and for each surface
it is given by :
dt qg
Q = – kA = – kA ( L 2 x)
dx x 0 or x L 2k x 0 or x L
AL
i.e., Q= qg
2
15. When both the surfaces are considered
AL
Q= 2× qg = AL qg
2
16. Also heat conducted to each wall surface is further dissipated to the
surroundings at temperature ta,
AL
Thus, qg = hA (tw – ta)
2
q
or tw = ta + g L ...(1.20.5)
2h
17. Substituting value of eq. (1.20.5) in eq (1.20.4), we get
L L2
tmax = ta + qg
2h 8 k
Answer
1. Consider a cylindrical rod in which one-dimensional radial conduction
is taking place under steady state condition.
2. Heat conducted in at radius r,
dt
Qr = – k 2rL
dr
Introduction to Heat Transfer 1–36 B (ME-Sem-5)
r dr
qg
L
ta
= k
4 = 2
=1–
tmax tw qg 2 R R
·R
4k
2
t tw r
i.e., = 1–
tmax tw R
13. Also,
Introduction to Heat Transfer 1–38 B (ME-Sem-5)
Answer
1. Consider one dimensional radial conduction of heat, under steady state
conditions through a sphere having uniform heat generation.
2. Let, R = Outside radius of sphere,
k = Thermal conductivity (uniform),
qg = Uniform heat generation per unit volume, per
unit time within the solid,
tw = Temperature of the outside surface (wall) of
the sphere, and
ta = Ambient temperature.
3. Consider an element at radius r and thickness, dr as shown in
Fig. 1.22.1.
Element
qg tw
ta
r dr
R
Fig. 1.22.1. Sphere with uniform heat generation.
4. Heat conducted in at radius r,
Heat & Mass Transfer 1–39 B (ME-Sem-5)
dt dt
k 4r 2
Qr = kA
dr dr
Heat generated in the element,
Qg = qg × A × dr = qg × 4r2 × dr
Heat conducted out at radius (r + dr)
d
Q(r + dr) = Qr (Qr ) dr
dr
5. Under steady state conditions, we have
Qr + Qg = Q(r + dr)
d
= Qr (Qr ) dr
dr
d
Qg = (Qr ) dr
dr
d dt
qg × 4r2 × dr = 4 k r 2 dr
dr dr
d 2 dt
qg × 4r2 × dr = – 4 k r dr
dr dr
1 d 2 dt qg
r =0
r 2 dr dr k
1 2 d 2t dt qg
r 2 2r =0
r2 dr dr k
d 2 t 2 dt qg
=0
dr 2 r dr k
d2t dt qg r
r 2 2 =0
dr dr k
d 2 t dt dt qg r
or r 2 =0
dr dr dr k
dt dt dt qg r
or r =0 ...(1.22.1)
dr dr dr k
6. Integrating the eq. (1.22.1), we get
dt q r2
r t g = C1
dr k 2
d q r2
or ( rt) g = C1 ...(1.22.2)
dr k 2
7. Integrating the eq. (1.22.2), we get
r3
qg
rt = C1r + C2
...(1.22.3)
k 6
8. At the centre of sphere, r = 0
C2 = 0
Introduction to Heat Transfer 1–40 B (ME-Sem-5)
qg R
or tw = ta ...(1.22.6)
3h
15. Putting the value of eq. (1.22.6) in eq. (1.22.5), we get
The maximum temperature,
qg qg
tmax = ta R R2
3h 6k
PART-3
Concept of Thermal Resistance, Analogy between Heat and
Electricity Flow, Thermal Contact Resistance, Overall Heat
Transfer Coefficient, Critical Radius of Insulation.
Questions-Answers
Answer
A. Analogous System :
1. When two physical systems are described by similar equations and
have similar boundary conditions, then the system are said to be
analogous.
Introduction to Heat Transfer 1–42 B (ME-Sem-5)
2. The heat transfer process may be compared by analogy with the flow
of electricity in an electrical system.
B. Thermal Resistance :
i. For Conduction :
As per Ohm’s law :
Potential difference dV
Current, I= =
Electrical resistance R
By analogy, for conduction
Temperature difference (dt)
Q=
dx / kA
dx
Rth =
kA
Rth = Thermal conduction resistance.
i R i
Q Rth Q
dx
Rth =
kA
Fig. 1.23.1.
ii. For Convection : Convection heat transfer is given by Newton’s law
of cooling i.e.,
Qconv. = h A(ts – tf)
(t t )
Qconv. = s f
1 / hA
1
Rth =
conv. hA
1
Rth =
ts hA tf
Fig. 1.23.2.
iii. For Radiation : Radiation heat transfer is governed by Stefan
Boltzmann law i.e.,
Q = AT4
or Q = A (T14 – T24) (For real surface)
(T1 T2 )
or Q=
1
A (T1 T2 )(T12 T22 )
1
So, (Rth)Rad. =
A (T1 T2 )(T12 T22 )
Heat & Mass Transfer 1–43 B (ME-Sem-5)
Answer
Answer
1. In real systems, due to surface roughness and void spaces, the contact
surfaces touch only at discrete locations. Thus, the area available for
flow of heat at the interface will be small as compared to geometric
face area. Due to this reduced area and presence of air voids, a large
resistance to heat flow at the interface occurs. This resistance is known
as thermal contact resistance.
Mathematically,
t t
Rc =
Q qA
2. It causes temperature drop between two materials at the interface.
Que 1.26. Derive an expression for overall heat transfer coefficient
through a plane wall.
Answer
1. Heat transfer between hot fluid and surface-1 of the wall take place
due to convection. The amount of heat transfer is given as,
Introduction to Heat Transfer 1–44 B (ME-Sem-5)
t1
Q
h cf
Hot hhf t2 Cold
fluid fluid film
film
tcf
1 2
L
thf t1 t2 tcf
1 L 1
h hf . A k.A hcf . A
Que 1.27. Derive the expression for steady state overall heat
transfer coefficient for a composite sphere whose inner surface is
at higher temperature and outer surface is at lower temperature.
Answer
1. Considering Fig. 1.27.1 as cross-section of a composite sphere, the heat
flow equation can written as follows
4 kA r1r2 (t1 t2 )
Q = hh 4r12 (th – t1) =
f f (r2 r1 )
4 kB r2 r3 (t2 t3 )
= = hc 4 r32 (t3 – tc )
(r3 r2 ) f f
2. By rearranging, we get
Q
(th – t1) = ...(1.27.1)
f 4 r12 hhf
Q (r2 r1 )
(t1 – t2) = ...(1.27.2)
4 kA r1 r2
Q (r3 r2 )
(t2 – t3) = ...(1.27.3)
4 kB r2 r3
Q
(t3 – tc ) = ...(1.27.4)
f 4 r32 hcf
3. Adding eq. (1.27.1), (1.27.2), (1.27.3) and (1.27.4), we get
Q 1 (r r ) (r r ) 1
(th – tc ) = 2
2 1 3 2 2
f f 4 h r k r r k r r h
hf 1 A 1 2 B 2 3 cf r3
(thf tcf )
So, Q= ...(1.27.5)
1 1 (r2 r1 ) (r3 r2 ) 1
4 hhf r12 kA r1 r2 kB r2 r3 hcf r32
6. If U is the overall heat transfer coefficient. Then, heat transfer is given
as,
Q = U A(th – tc ) ...(1.27.6)
f f
7. Comparing the eq. (1.27.5) and (1.27.6), we get
Introduction to Heat Transfer 1–46 B (ME-Sem-5)
1
4
U= 1 (r2 r1 ) r3 r2 1
A
hh r12 kA r1 r2
kBr2 r3 hcf r32
f
Cold fluid (air) t cf
Q
hcf
B tcf
Hot t A
fluid hf
hhf t1 t2 t3
r1
r2
r3
Fig. 1.27.1. Steady state conduction through a composite sphere.
Answer
1. The thickness upto which heat flow rate increases and after which
heat rate flow decreases is termed as critical thickness of insulation.
2. In case of cylinders and spheres it is called critical radius (rc).
i. For cylinder, rc = k / h
ii. For sphere, rc = 2k / h
3. The addition of insulation always increases the conductive thermal
resistance.
4. But when the total thermal resistance is made of conductive thermal
resistance [(Rth)cond.] and convective thermal resistance [(Rth)conv.], the
addition of insulation in some cases may reduce the convective thermal
resistance due to increase in surface area, as in the case of a cylinder
and sphere, and the total thermal resistance may actually decrease
resulting in increased heat flow.
Que 1.29. Prove that the critical radius of insulation for cylinder
k 2k
is, rc = and for sphere is, rc = ·
ho ho
Heat & Mass Transfer 1–47 B (ME-Sem-5)
Answer
A. Critical Radius of Insulation for Cylinder :
1. Consider a solid cylinder of radius r 1 insulated with an insulation of
thickness (r 2 – r 1).
2. Then the rate of heat transfer from the surface of the solid cylinder
to the surroundings is,
2 L(t1 tair )
Q= ...(1.29.1)
r
ln 2
r1 1
k ho r2
Solid cylinder
k h0
t1
tair
r1 r2 – r1
r2
ln (r2 / r1 )
3. From eq. (1.29.1) it is evident that as r2 increases, the factor
k
1
increases but the factor decreases. Then, the required condition
ho r2
is,
r2
ln
d r1 1
=0
dr2 k ho r2
1 1 1 1
+ = 0
k r2 ho r22
Introduction to Heat Transfer 1–48 B (ME-Sem-5)
1 1
– =0
k ho r2
hor2 = k
k
r2 = (rc) =
ho
B. Critical Radius of Insulation for Sphere :
1. Consider a solid sphere which is insulated by a thickness of insulation
(r2 – r1).
2. Then the rate of heat transfer from sphere to surrounding is given by,
(t1 tair )
Q=
(r2 r1 ) 1
4 kr r
1 2 4 r22ho
3. For critical radius :
d (r2 r1 ) 1
=0
dr2 4 kr1 r2 4 r22 ho
1 2
=0
kr22 r23 ho
2kr2 = r23ho
2
2k
r2 = rc =
ho
Insulation
Solid sphere
k h0
t1 tair
r1
r2 r2 – r1
Answer
6.5
Given : r1 = = 3.25 mm = 0.00325 m
2
k = 0.174 W/m-°C, ho = 8.722 W/m2-°C, tair = 20 °C, t1 = 60 °C.
To Find : i. Minimum thickness of insulation.
ii. Heat loss per meter length.
iii. Percentage increase in the heat dissipation.
1. We know that, the critical radius of insulation of cylinder
k 0.174
rc = = = 0.01995 m = 19.95 mm
ho 8.722
2. Minimum thickness of insulation = rc – r1 = 19.95 – 3.25 = 16.7 mm
3. Heat loss without insulation,
2L (t1 tair ) 2 1 (60 20)
Q1 = = = 7.124 W/m
1 1
ho r1 8.722 0.00325
4. Heat loss with insulation when critical thickness is used,
2 L (t1 tair ) 2 1 (60 20)
Q2 = =
r 0.01995
ln c ln
0.00325 1
r1 1
k ho rc 0.174 8.722 0.01995
= 15.537 W/m
5. Percentage increase in heat dissipation,
Q2 Q1 (15.537 7.124)
= × 100 = × 100 = 118.09 %
Q1 7.124
Heat & Mass Transfer 2–1 B (ME-Sem-5)
2 Fins &
Transient Conduction
Questions-Answers
Que 2.1. What do you understand by fins and also explain types
of fins ?
Answer
A. Fins :
1. These are extended surfaces used for increasing the heat transfer rate
from the surfaces whenever the available surface is found inadequate
to transfer the required quantity of heat with the available temperature
drop and convective heat transfer coefficient.
2. Heat transfer through fin is because of conduction and convection
between its boundary and surrounding.
B. Types of Fins : Fins are of two types :
1. Fins of uniform cross-sectional area :
a. Fins of rectangular profile, and
b. Fins of circular profile.
2. Fins of non-uniform cross-sectional area :
a. Longitudinal fin of trapezoidal profile,
Heat & Mass Transfer 2–3 B (ME-Sem-5)
Que 2.2. What is the reason for the widespread use of fins on
Answer
1. Generally fins are used to provide extra surface area to increase rate of
heat transfer from existing surface.
2. In most of practical situation the heat transfer fluids involves liquid and
gas.
3. In such a situation of heat transfer between liquid and gas, as gas is
having less value of heat transfer coefficient it offers greater resistance
for the heat transfer.
4. So to increase the rate of heat transfer on the gas side, extra surface
area is provided by using fins.
5. Use of fins on surfaces is found in case of air-cooled cylinder, condenser
tubes in refrigerator.
Answer
A. Temperature Distribution and Heat Transfer Rate for
Rectangular Fins :
1. For the proper design of fins, the knowledge of temperature distribution
along the fin is necessary.
2. For analysis of heat flow through fin, following assumptions are made :
i. Steady state heat conduction.
Fins & Transient Conduction 2–4 B (ME-Sem-5)
P
Qconv.
t0 l
k y
x
dx
Qx Q(x + dx ) b
ACS
h, ta
Wall
Element
x
Fig. 2.3.1. Rectangular fin of uniform cross-section.
dt
Qx + dx = – kACS
dx x dx
9. Heat convected out of the element between the plane x and (x + dx),
Qconv. = h (P dx) (t – ta) ...(2.3.1)
10. Applying the energy balance on the element, we have
Heat & Mass Transfer 2–5 B (ME-Sem-5)
Qx = Qx + dx + Qconv.
dt dt
– kA CS = – kA CS + h (P dx) (t – ta) ...(2.3.2)
dx x dx x dx
11. By Taylor’s expansion,
dt dt d dt d 2 dt ( dx)2
dx =
dx 2 ... ...(2.3.3)
x dx dx x dx dx x dx dx x 2
Substituting this in eq. (2.3.2), we have
dt dt d dt
– kA CS = – kA CS – kA CS dx
dx x dx x dx dx x
d 2 dt (dx)2
– kA CS ... h( P dx)(t – ta )
dx2 dx x 2
Neglecting higher term as dx 0,
dt dt d2t
– kA CS =– kA CS – kA CS 2 dx + h (P dx) (t – ta)
dx x dx x dx
d2t
kA CS 2 dx – h (P dx) (t – ta) = 0
dx
d2t hP
or – (t – ta ) 0 ...(2.3.4)
dx2 kACS
12. If temperature excess is,
= t – ta ...(2.3.5)
Since, ambient temperature ta is constant
d dt
=
dx dx
d 2 d 2t
or = ...(2.3.6)
dx 2 dx 2
13. Putting the value of eq. (2.3.5) and eq. (2.3.6) in eq. (2.3.4), we get
d 2
– m2 = 0 ...(2.3.7)
dx 2
hP
Where, m=
kACS
14. Eq. (2.3.4) and eq. (2.3.7) shows the general differential heat conduction
equation for fins of uniform cross-sectional area.
15. The general solution of this linear and homogeneous second order
differential equation is,
= C1 emx + C2 e–mx ...(2.3.8)
or t – ta = C1 emx + C2 e–mx
Fins & Transient Conduction 2–6 B (ME-Sem-5)
dt
Qfin = – kACS
dx x 0
dt
Putting the value of from eq. (2.3.12) in above equation, we
dx x 0
get
Qfin = Ph kACS (t0 – ta)
Case III : The Fin is of Finite Length and Losses Heat by Convection :
1. The boundary conditions are :
i. At x = 0, = 0
ii. Heat conducted to the fin at (x = l) = Heat convected from the end
to surrounding.
dt
i.e., – kA CS = h ASA(t – ta)
dx x l
Where, ASA is surface area from which the convective heat transfer
takes place and at the tip of fin, A CS = ASA.
dt h
Thus, =– at x = l
dx k
2. Applying these boundary condition in eq. (2.3.8), we get
h
t ta cosh[ m(l x)] [sinh{m (l x)}]
= = km
o to ta h
cosh(ml) [sinh(ml)]
km
...(2.3.13)
3. Differentiating eq. (2.3.13) w.r.t. x, we get
h
dt sinh {ml} km cosh {ml}
dx = – m (to ta ) ...(2.3.14)
x 0 cosh {ml} h sinh {ml}
km
4. Rate of heat flow through the fin is given by :
dt
Qfin = – kACS
dx x 0
dt
Putting the value of from eq. (2.3.14) in above equation, we
dx x 0
get
h
tanh( ml) km
Qfin = Ph kACS (to – ta)
1 h tanh (ml)
km
Fins & Transient Conduction 2–8 B (ME-Sem-5)
B. Advantages :
1. By using the fin, heat transfer rate can be increased without any
preventive maintenance.
2. It is the cheapest way for increasing the heat transfer rate from the hot
bodies.
C. Applications :
1. Economisers for steam power plants,
2. Convectors for steam and hot-water heating systems,
3. Radiators for automobiles, and
4. Air cooled engine cylinder heads.
Answer
A. Efficiency of Fin :
1. It is defined as the ratio of the actual heat transferred by the fin to the
maximum heat transferable by fin, if entire fin area were at base
temperature.
Actual heat transfer by the fin(Qfin )
fin =
Maximum heat that would be transferred if whole surface
of the fin is maintained at the base temperature (Qmax )
2. For a fin which is infinitely long,
PhkACS (to ta ) kACS 1
fin =
Phl(to ta ) hPl 2 ml
3. For a fin which is insulated at the tip,
Que 2.5. Two pin fins are identical, except that the diameter of
one of them is twice the diameter of the other. For which fin will the
(a) effectiveness, and (b) fin efficiency be higher ? Explain.
AKTU 2014-15, Marks 05
Answer
Let the diameter of fin I be d1 and diameter of fin II be d2.
Given : d1 = 2d2
a. Effectiveness :
1. The effectiveness of a fin is given as follows,
h cosh mL mk sinh mL
kAm [ ts t ]
mk cosh mL h sinh mL
=
hA (ts t )
kA1 m(ts t ) h cosh mL mk sinh mL
fin hA(ts t ) mk cosh mL h sinh mL
2. Now, I =
fin II kA2 m(ts t ) h cosh mL mk sinh mL
hA(ts t ) mk cosh mL h sinh mL
2
A d1
fin I = 1 4
fin II
A2 d 2
4 2
fin I (2 d2 )2
= 4
fin II 2
d
4 2
fin I = 4 fin II
3. So, effectiveness of fin I is higher than effectiveness of fin II.
b. Efficiency :
1. The efficiency of a fin is given as follows,
1 A
fin = (t t ) dA
A [ ts t ] 0
1 A
2. Now,
fin I
=
A1 [ ts t ] 0 [t t ] dA A2
fin II 1 A
A2 [ ts t ]
0 [t t ] dA A1
2 2
fin I d2 d2
= 4 4
fin II 2
d (2 d2 )2
4 1 4
Fins & Transient Conduction 2–10 B (ME-Sem-5)
fin I 1 1
= fin I fin II [fin I = 25 % of fin II]
fin II 4 4
3. So, efficiency of fin II is higher than efficiency of fin I.
Answer
ta
to
0.1 m
0.02 m
Fig. 2.6.1.
hP 10 d 10 4
1. We know that, m=
kACS 2 20 0.02
20 d
4
m = 10
2. Rate of heat transfer,
Qfin = PhkACS (to ta ) tanh (ml)
t ta 1
=
to ta cosh (ml)
t 50 1
=
300 50 cosh (10 0.1)
Temperature at fin tip, t = 212 ºC.
Que 2.7. Aluminium fins 1.5 cm wide and 1.0 mm thick are placed
on 2.5 cm diameter tube to dissipate the heat. The tube surface
temperature is 170 °C and the ambient fluid temperature is 25 °C.
Calculate the heat loss per fin for h = 130 W/m2-°C for aluminium.
AKTU 2014-15, Marks 10
Answer
Answer
A. Thermometer Well : It is defined as a small tube welded radially into
a pipeline through which a fluid whose temperature is to be measured
is flowing.
B. Estimation of the Error in Temperature Measurement in a
Thermometer Well :
1. Let, l = Length of the well/tube,
d = Internal diameter of the well/tube,
= Thickness of well/tube,
Fins & Transient Conduction 2–12 B (ME-Sem-5)
Pipeline
t0 (at x = 0)
oil
Fluid flow
tf
tl (at x = l )
d
l
Well
temperature
x
Fig. 2.8.2. Temperature variation in well.
Heat & Mass Transfer 2–13 B (ME-Sem-5)
6. At x = l, we have
tl tf cosh [ m( l l)] 1
= [Thermometric error]
to tf cosh (ml) cosh (ml)
Where, tl = Temperature recorded by the thermometer at the
bottom of the well.
7. Now, perimeter of the well,
P = (d + 2) d
and cross-sectional area,
Acs = d
P d 1
=
Acs d
hP h
8. Then, m=
kAcs k
Thus, the temperature measured by the thermometer is not affected
by the diameter of the well.
Answer
h 30
1. We know that, m= = 24.5
k 50 0.001
2. Now, ml = 24.5 × 0.14 = 3.43
3. The true temperature of air is found from the relation,
tl tf 1
=
to t f cosh (ml)
100 tf 1 1
= = 0.0647
50 tf cosh (3.43) 15.4545
(100 – tf) = 0.0647 (50 – tf)
100 – tf = 3.235 – 0.0647 tf
tf = 103.46 °C.
Heat & Mass Transfer 3–1 B (ME-Sem-5)
3 Forced &
Natural Convection
Questions-Answers
Answer
1. The layer adjacent to the boundary is known as boundary layer.
Boundary layer is formed whenever there is relative motion between
the boundary and the fluid.
2. When a real fluid (viscous fluid) flows past a stationary solid boundary,
a layer of fluid which comes in contact with the boundary surface,
adheres to it and condition of no slip occurs.
3. Thus the layer of fluid which cannot slip away from the boundary surface
undergoes retardation, this retarded layer further causes retardation
for the adjacent layers of the fluid, thereby developing a small region in
the immediate vicinity of the boundary surface in which the velocity of
the flowing fluid increases rapidly from zero at the boundary surface
and approaches the velocity of main stream.
4. According to boundary layer theory, the extensive fluid medium around
bodies moving in fluids can be divided into following two regions :
Heat & Mass Transfer 3–3 B (ME-Sem-5)
i. A thin layer adjoining the boundary called the boundary layer where
the viscous shear takes place.
ii. A region outside the boundary layer where the flow behaviour is quite
like that of an ideal fluid and the potential flow theory is applicable.
Answer
The equation for boundary layer thickness can be obtained in terms of
displacement, momentum and energy that are discussed as below :
a. Displacement Thickness :
1. It can be defined as the distance, measured perpendicular to the boundary
by which the main/free stream is displaced on account of formation of
boundary layer. It is denoted by *.
Boundary layer U
dy u = 0.99U Velocity
distribution
U
y
Stationary plate
Fig. 3.2.1. Displacement thickness.
2. Let fluid of density flow past a stationary plate with velocity U as
shown in Fig. 3.2.1.
3. Consider an elementary strip of thickness dy at a distance y from the
plate.
4. Mass flow per second through the elementary strip = udy
5. Mass flow per second through elementary strip, if the plate was not
there = Udy
6. Reduction of mass flow rate through elementary strip
= (U – u) dy
7. Total reduction of mass flow rate due to introduction of plate
= (U u) dy ...(3.2.1)
0
8. Let the plate is displaced by a distance * and velocity of flow for the
distance * is equal to the main/free stream velocity. Then, loss of mass
of fluid/sec flowing through the distance *
= U* ...(3.2.2)
Forced & Natural Convection 3–4 B (ME-Sem-5)
Answer
1. Whenever a flow of fluid takes place past a heated or cold surface, a
temperature field is set up in the field next to the surface. If the
surface of the plate is hotter than fluid, the temperature distribution
will be as shown in the Fig. 3.3.1. The zone or the layer wherein the
temperature field exists is called the thermal boundary layer.
Limit of thermal
boundary layer
Region
affected by Temperature
U, t heating
Y profile, t(y)
Free th
stream ( t s – t )
t (y ) Plate
surface (ts)
X
t t
Laminar Turbulent
Fig. 3.3.1. Thermal boundary layer formed during
flow of cool fluid over a warm plate.
2. Due to the exchange of heat between the plate and the fluid,
temperature gradient occurs.
Forced & Natural Convection 3–6 B (ME-Sem-5)
Que 3.4. Derive the differential energy equation for flow past a
flat plate. Also write down assumptions of it.
OR
Derive an expression for energy equation of thermal boundary layer
over flat plate. AKTU 2016-17, Marks 7.5
Answer
A. Assumptions :
1. The flow is steady and incompressible.
t ts
2. Fluid’s properties evaluated at the film temperature, tf = are
2
constant.
3. The body forces, viscous heating and conduction in the flow direction
are negligible.
B. Derivation :
1. For an element of dimensions (dx × dy × Unit depth) in the boundary
layer the quantities of energy entering and leaving are shown in
Fig. 3.4.1.
–kdx t + t dy
y y y
Viscous heat
Ey2 = v+ v dy dx C t + t dy
generation y p
y
u 2
dxdy
y
dx
7. Now, for steady state condition, the algebraic sum of heat due to
convection, conduction and viscous effect equals to zero. Thus
2
t u t v 2 t u
C p u t v t dxdy + k dxdy + µ dxdy = 0
x x y y y2 y
2
t t u v 2t u
or – Cp u v t k 2 = 0 ...(3.4.1)
x y x y y y
u v
8. For a two-dimensional boundary layer flow, = 0 and
x y
therefore eq. (3.4.1) can be written as,
2
t t k 2t u
u v =
x y C p y2 C p y
which represents the differential energy equation for flow past a flat
plate.
9. If the heat generation due to viscous effects is negligible, then the
energy equation can be re-written as
t t k 2t 2t
u v = =
x y C p y2 y2
k
Where, = = Thermal diffusivity
C p
Answer
1. Consider a control volume as shown in Fig. 3.5.1.
Hydrodynamic
boundary layer
Qh Thermal
U , t boundary layer
B C
H
Q x + dx
Qx
Y th
Control
volume
X A D
xo Q cond.
x dx Thin flat plate
H H
ii. Mass of fluid leaving through face CD = u dy dx
u dy
x 0 0
iii. Mass of fluid entering the control volume through face BC
H
= 0 u dy
x
0
H
0
H
u dy dx – u dy
H
= u dy dx
x 0
4. Heat influx through the face AB,
Qx = Mass × Specific heat × Temperature
or, Qx = 0
H
u dy C pt C p ut dy
0
H
H H
C p ut dy dx
= C p 0 ut dy
x 0
9. After simplification and rearrangement, we have
d H k t t
(t – t)u dy =
dx 0
...(3.5.1)
C p y y 0 y y 0
Eq. (3.5.1) is the integral equation for the boundary layer for constant
properties and constant free stream temperature t.
Forced & Natural Convection 3–10 B (ME-Sem-5)
Answer
1. As we know that, the inter-relationship between fluid friction and
Newton’s law of viscosity is given as,
du
0 = ...(3.6.1)
dy
2. Heat flow along Y-direction follows the Fourier equation
dt
Q = kA ...(3.6.2)
dy
3. When Pr is unity, temperature and velocity profile are identical (for
most of the gases, 0.6 < Pr < 1.0)
C p k
i.e., = 1 or = Cp ...(3.6.3)
k
4. By combining eq. (3.6.1), eq. (3.6.2) and eq. (3.6.3), we get
dt
Q = – C p A 0
du
5. Separating the variables and integrating within the limits :
At the plate surface : u = 0 and t = ts
At the outer edge of boundary layer : u = U and t = t
Q U t
du = – t dt
C p A 0 0 s
Q
U = (ts – t)
C p A 0
Q 0 Cp
= ...(3.6.4)
A (ts t ) U
Q
6. But, = hx ...(3.6.5)
A (ts t )
1
and, U 2
0 = C fx ...(3.6.6)
2
7. Substituting eq. (3.6.5) and eq. (3.6.6) in eq. (3.6.4), we get
Heat & Mass Transfer 3–11 B (ME-Sem-5)
1 C p Cfx
hx = Cfx U 2 ( C p U )
2 U 2
hx C
= fx
Cp U 2
hx
is called the Stanton number Stx.
Cp U
8. It represents the Nusselt number divided by the product of the Reynolds
and Prandtl number, i.e.,
Nux C fx
= Stx ...(3.6.7)
Rex Pr 2
Eq. (3.6.7) is called the Reynolds analogy.
9. Further in case of laminar boundary layer on a flat plate, we have
hx x
Nux = 0.332 ( Rex )1/ 2 ( Pr)1/ 3 ...(3.6.8)
k
10. Dividing both sides of the eq. (3.6.8) by Rex (Pr)1/3, we get
Nux 0.332 Cfx
= ...(3.6.9)
Rex ( Pr)1/ 3 ( Rex )1/ 2 2
11. The LHS of the eq. (3.6.9) can be rewritten as
Nux Nux
= ( Pr)2/ 3 Stx (Pr)2/ 3
Rex ( Pr)1/ 3 Rex Pr
C fx
Stx (Pr)2/3 = ...(3.6.10)
2
Eq. (3.6.10) shows the inter-relationship between heat and momentum
transfer and designated as Colburn analogy.
12. For Pr = 1, the Reynolds and Colburn analogies are same.
Answer
1. Fig. 3.7.1 shows flow over a plate.
Turbulent boundary layer
Mixed G
Laminar boundary boundary F
layer layer
U
E Laminar
sub-layer
B C D
A
Laminar Transition Turbulent
zone zone zone
Leading edge
Fig. 3.7.1.
Forced & Natural Convection 3–12 B (ME-Sem-5)
2. If the length of the plate is more than the distance x, calculated from
equation, Ux/ = 5 × 105, the thickness of boundary layer will go on
increasing in the downstream direction.
3. Then the laminar boundary layer becomes unstable and motion of fluid
within it is disturbed and irregular which leads to a transition from
laminar to turbulent boundary layer.
4. This short length over which the boundary layer flow changes from
laminar to turbulent is called transition zone or mixed zone, which is
shown by distance BC.
5. Further downstream the transition zone, the boundary layer is turbulent
and continues to grow in thickness.
6. This layer of boundary is called turbulent and continues to grow in
thickness which is shown by the portion FG.
PART-2
Flow over a Flat Plate, Flow across a Single Cylinder and a Sphere,
Flow Inside Ducts, Thermal Entrance Region, Empirical Heat
Transfer Relations, Relation between Fluid Friction and Heat
Transfer, Liquid Metal Heat Transfer.
Questions-Answers
Answer
a. Flow over Flat Surface : The general equation giving to the local heat
transfer coefficient for turbulent flow over flat plate is given by,
Heat & Mass Transfer 3–13 B (ME-Sem-5)
Answer
1. The following dimensionless numbers are used in forced convection :
hL
i. Nusselt number, Nu =
k
L v
ii. Reynold’s number, Re =
C p
iii. Prandtl number, Pr =
k
h
iv. Stanton number, St =
C p v
2. In order to determine the value of convection coefficient h, the following
equations are used :
Nu = f1 (Re, Pr) = C1 (Re)m (Pr)n
St = f2 (Re, Pr) = C2 (Re)a (Pr)b
Que 3.10. Describe the relation between fluid friction and heat
transfer. How is the average friction and heat transfer coefficients
determined in flow over a flat plate ? AKTU 2014-15, Marks 10
Answer
A. Relation between Fluid Friction and Heat Transfer : Refer Q. 3.6,
Page 3–10B, Unit-3.
B. Expression for Average Friction :
1. The local skin friction coefficient is given as
0 0.646
Cfx =
1 2 Rex
U
2
2. Average value of skin friction coefficient,
1 L 1 L 0.646 dx
Cf = C fx dx =
L 0 L 0 U / x
1.292
or, Cf = 1.292
LU ReL
Forced & Natural Convection 3–14 B (ME-Sem-5)
LU
Where, ReL = is Reynolds number based on total length L of the
flat plate.
C. Expression for Heat Transfer Coefficient :
Q dt
1. We know that, = hx(ts – t) = – k
A dy y 0
k(dt / dy) y0
or, hx =
ts t
dt 3 t t
2. But, dy = s
y 0 2 th
3 t t
k s
2 th 3k 3k 1
hx = = ...(3.10.1)
(ts t ) 2 th 2 r
3/ 4 1/ 3
4.64 x
3. Substituting r = 0.975 1 x0 and = in eq. (3.10.1), we
( Pr ) x
1/ 3
Rex
get
3k ( Pr )1/ 3 Rex
hx =
2 3 / 4 1/ 3 4.64 x
x
0.975 1 0
x
k 1
hx = 0.332( Pr)1/ 3 ( Re)1/ 2 1/ 3
x x0 3/ 4
1
x
4. When the plate is heated over the whole length i.e., x0 = 0, we have
k
hx = 0.332 ( Pr )1/ 3 ( Rex )1/ 2
x
Que 3.11. Explain drag force and drag coefficient.
Answer
A. Drag Force (FD) :
1. The drag force on elemental area
= Force due to pressure in the direction of fluid motion
+ Force due to shear stress in the direction
of fluid motion
= pdA cos + 0dA cos (90° – ) = pdA cos + 0dA sin
2. The term p cos dA is called the pressure drag or form drag while the
term 0 sin dA is called the friction drag or skin drag or shear drag.
0dA
FL FR
pdA
U X 90°– FD
X
Stationary body
Fig. 3.11.1.
B. Local Coefficient of Drag : It is defined as the ratio of the shear stress
1
0 to the quantity U2. It is denoted by C*D.
2
0
Hence, C*D =
1
U 2
2
C. Average Coefficient of Drag : It is defined as the ratio of the total drag
force to the quantity 1/2 AU2. It is also called coefficient of drag and it
is denoted by CD.
FD
Hence, CD =
1 / 2 AU 2
Where, A = Area of the surface (or plate),
U = Free stream velocity, and
= Mass density of fluid.
Que 3.12. Air at 27 °C and 1 atm flows over a flat plate at a velocity
3 m/s. The plate is heated over its entire length to a temperature of
70 °C. Calculate the heat transferred if the plate length is 45 cm
and width is 1 m. Properties of air, = 17.36 × 10 –6 m 2/s ,
k = 0.02749 W/m-K , C p = 1.006 kJ/kg-K , Pr = 0.7. Take,
NuL = 0.664 ReL0.5 Pr 1/ 3 . AKTU 2013-14, Marks 10
Answer
UL 3 0.45
1. Reynolds number, Re = = 77764.98
17 . 36 10 6
Re < 5 × 105, hence flow is laminar.
2. Average Nusselt number,
NuL = 0.664 Re0.5 Pr1/3
hL
= 0.664 (77764.98)0.5 (0.7)1/3
k
0.02749 0.664 (77764.98)0.5 (0.7)1/3
h =
0.45
h = 10.044 W/m2-ºC
3. Heat transfer, Q = h As(ts – t)
= 10.044 × 0.45 × 1 × (70 – 27) = 194.351 W
Answer
h = 72.57 W/m2-K
3. Heat transferred per unit length of the tube,
Q = hAs(tw – t) = 72.57 × d × 1 × (220 – 200)
= 72.57 × × 0.03 × 20 = 136.79 W/m
4. Bulk temperature increase of air over a 4 m length of the tube
Q = mCp (t) = (AU) Cp (t)
4 × 136.79 = 1.493 × × (0.03)2 × 12 × 1.025 × 103 × (t)
4
t = 42.15 ºC
Que 3.14. Air is flowing over a flat plate 5 m long and 2.5 m
wide with a velocity of 4 m/s at 15 °C. If = 1.208 kg/m 3 and
= 1.47 × 10–5 m 2/s, calculate the length of plate over which the
boundary layer is laminar and thickness of the boundary layer
(laminar), shear stress at the location where boundary layer ceases
to be laminar and the total drag force on the both sides of that
portion of the plate where boundary layer is laminar.
AKTU 2014-15, Marks 10
Answer
Given : L = 5 m, w = 2.5 m, U = 4 m/s, = 1.208 kg/m3,
= 1.47 × 10–5 m2/s
To Find : i. Length of plate over which the boundary layer is
laminar, and thickness of the boundary layer
(laminar).
ii. Shear stress at the location where boundary layer
ceases to be laminar.
iii. Total drag force on the both sides of that portion of
plate where boundary layer is laminar.
UL 45
1. Reynold number, Re = = 1.361 × 106.
1.47 10 5
Hence on the front portion, boundary layer is laminar and on the rear,
it is turbulent.
Ux
Re = 5 10 5
4x
= 5 × 105
1.47 10 –5
5 105 1.47 10 5
x= = 1.837 m
4
Hence the boundary layer is laminar on 1.837 m length of the plate.
5x
2. Thickness of the boundary layer (laminar), =
Re
Forced & Natural Convection 3–18 B (ME-Sem-5)
5 1.837
= = 0.01299 m or 12.99 mm
5 105
0.664 0.664
3. Local coefficient of drag, Cfx = = 0.000939
Re 5 105
1
Shear stress, 0 = Cfx U 2
2
1
= 0.000939 × × 1.208 × 42 = 0.00907 N/m2
2
4. Total drag force on both sides of plate,
1
FD = 2C f AU 2
2
Where, Cf = Average coefficient of drag (or skin friction)
1.328 1.328
= = 1.878 × 10–3
Re 5 105
and area of the plate, A = xL = 1.837 × 2.5 = 4.59 m2
1
FD = 2 × 1.878 × 10–3 × × 1.208 × 4.59 × 42
2
= 0.167 N
Que 3.15. A flat plate is 2 m long, 0.8 m wide and 3 mm thick. Its
density and specific heat is 3000 kg/m3 and 700 J/kg-K respectively.
The plate is having initial temperature of 80 °C. A stream of air at
20 °C blown over both surfaces of the plate along its width, at a
velocity of 2 m/s. Calculate rate of heat dissipation from the plate
and initial rate of cooling of the plate. The properties of air are :
= 1.09 kg/m3, k = 0.028 W/m-°C, C = 1007 J/kg-K,
= 2.03 × 10 – 5 kg/m-s, Pr = 0.698
Nu a = 0.664 (Rex)0.5 (Pr)0.33 for Laminar Flow
Nu a = 0.0336 (Rex)0.8 – 836 (Pr)0.33 for Turbulent Flow
AKTU 2015-16, Marks 7.5
Answer
Que 3.16. Castor oil at 25 °C flows at a velocity of 0.1 m/s past a flat
plate in a certain process. If the plate is 4.5 m long and is maintained
at a uniform temperature of 95 °C. Calculate the following using
exact solution :
i. The hydrodynamic and thermal boundary layer thickness on
one side of the plate,
ii. The total drag force per unit width on one side of the plate,
iii. The local heat transfer coefficient at the trailing edge, and the
heat transfer rate. AKTU 2017-18, Marks 10
Forced & Natural Convection 3–20 B (ME-Sem-5)
Answer
Answer
1. The example of liquid metal are sodium, potassium, lead, etc. and in the
liquid metal, Prandtl number is very less than 1 i.e., Pr << 1.
2. For constant temperature,
Nu = 4.8 + 0.0156 Pe0.85 Pr0.08
Where, Pe = Re Pr
3. For constant heat flux condition,
Nu = 4.82 + 0.0185 Pe0.827
PART-3
Physical Mechanism of Natural Convection, Buoyant Force
Empirical Heat Transfer Relations for Natural Convection over
Vertical Plates and Cylinders, Horizontal Plates and Cylinders
and Sphere, Combined Free and Forced Convection.
Questions-Answers
Answer
A. Buoyant Force : Whenever a body is immersed in a fluid, an upward
force is exerted by the fluid on the body. This upward force is equal to
the weight of the fluid displaced by the body and is called the buoyant
force or simply buoyancy.
G
B
FB
Fig. 3.18.1.
B. Centre of Buoyancy :
1. It is defined as a point about which the force of buoyancy is supposed to
act.
2. As the force of buoyancy is a vertical force and is equal to the weight of
the fluid displaced by the body, the centre of buoyancy will be the centre
of gravity of the fluid displaced.
3. In Fig. 3.18.1, FB is the force of buoyancy and point B is the centre of
buoyancy.
Answer
1. In free convection, the additional force is the body force g and hence
the momentum equation gets modified to
u u 2u p
u v = y 2 – g – x ...(3.19.1)
x y
Heat & Mass Transfer 3–23 B (ME-Sem-5)
Answer
A. Different Dimensionless Numbers used in Free Convection :
The dimensionless numbers used for empirical relations for heat transfer
in natural convection are :
hL
1. Nusselt number, (Nu) =
k
L3 g t
2. Grashoff number, (Gr) =
2
Cp
3. Prandtl number, (Pr) =
k
B. Empirical Heat Transfer Relations for Free Convection :
a. Vertical Plates : The commonly used relation are :
i. For Laminar Flow :
NuL = 0.59 (Gr Pr)1/4 for (104 < Gr Pr < 109)
ii. For Turbulent Flow :
NuL = 0.10 (Gr Pr)1/3 for (109 < Gr Pr < 1012)
Forced & Natural Convection 3–24 B (ME-Sem-5)
b. Horizontal Plates :
Case I : Upper surface heated or the lower surface cooled :
i. For Laminar Flow :
NuL = 0.54 (Gr Pr)1/4 for (105 < Gr Pr 2 × 107)
ii. For Turbulent Flow :
NuL = 0.14 (Gr Pr)1/3 for (2 × 107 < Gr Pr 3 × 1010)
Case II. Lower surface heated or upper surface cooled :
i. For Laminar Flow :
NuL = 0.27 (Gr Pr)1/4 for (3 × 105 < Gr Pr 3 × 1010)
ii. For Turbulent Flow :
NuL = 0.107 (Gr Pr)1/3 for (7 × 106 < Gr Pr 11 × 1010)
c. Cylinders : The relation for cylinders is given by,
i. For Laminar Flow :
NuL = 0.53 (Gr Pr)1/4 for (104 < Gr Pr < 109)
ii. For Turbulent Flow :
NuL = 0.13 (Gr Pr)1/3 for (109 < Gr Pr < 1012)
d. Spheres : The relation for sphere is given by,
NuL = 2 + 0.43 (Gr Pr)1/4 for (1 < Gr Pr < 105)
Answer
1. The heat transfer coefficient in case of natural or free convection, like
forced convection depends upon the variables v,, k, , Cp and L or D.
2. Since the fluid circulation in free convection is owing to difference in
density between the various fluid layers due to temperature gradient
and not by external agency, therefore, velocity v is no longer an
independent variable but depends upon the following factors :
i. t i.e., the difference of temperatures between the heated surface
and the undisturbed fluid.
ii. i.e., coefficient of volume expansion of the fluid.
iii. g i.e., acceleration due to gravity.
(gt is considered as one physical factor.)
3. Thus heat transfer coefficient ‘h’ may be expressed as follows :
h = f (, L, , Cp, k, gt)
or f1 (, L, , k, h, Cp, gt) = 0 ...(3.21.1)
The parameter (gt) represents the buoyant force and has the
dimensions of [LT–2].
Heat & Mass Transfer 3–25 B (ME-Sem-5)
( gt) 2 L3 ( gt) L3
or 3 = = Gr
2 2
7. Putting the values of 1, 2 and 3 in eq. (3.21.2), we have
Nu = f(Gr, Pr)
Que 3.22. A horizontal pipe 1 ft (0.3048 m) in diameter is maintained
at a temperature of 250 °C in a room where the ambient air is at
15 °C. Calculate the free convection heat loss per meter of length.
AKTU 2014-15, Marks 10
Answer
Que 3.23. A 350 mm long glass plate is hung vertically in the air at
24 °C while its temperature is maintained at 80 °C. Calculate the
boundary layer thickness at the trailing edge of the plate. If a similar
plate is placed in a wind tunnel and air is blown over it at a velocity
of 5 m/s, find the boundary layer thickness at its trailing edge. Also
determine the average heat transfer coefficient, for natural and
forced convection for the above mentioned data.
AKTU 2017-18, Marks 10
Answer
hL
NuL = = 0.664 (Re)1/2 (Pr)1/3
k
= 0.664 (9.505 × 104)1/2 (0.7)1/3 = 181.78
ii. Heat transfer coefficient,
k 28.15 103
h= 181.78 × 181.78
L 0.35
2
= 14.62 W/m -°C
4 Thermal Radiation
Questions-Answers
Answer
A. Radiation Heat Transfer :
1. It is defined as the transfer of energy from a system in the form of
electromagnetic mechanism which is caused by the temperature
difference.
2. Factors on which rate of emission of radiation by a body depends are as
follows :
i. Temperature of the surface,
ii. The nature of the surface, and
iii. Wavelength of radiation.
B. Surface Emission Properties :
a. Total Emissive Power (E) :
1. It is defined as the total amount of radiation emitted by a body per unit
area and time.
2. According to Stefan Boltzmann “the total emissive power of a black
body is proportional to the fourth power of its absolute temperature.”
Eb = t4 W/m2
Heat & Mass Transfer 4–3 B (ME-Sem-5)
Eb = At4 W
Where, = Stefan-Boltzmann constant.
= 5.67 × 10–8 W/m2-K4
b. Monochromatic Emissive Power (E ) : It is defined as the rate of
energy radiated per unit area of the surface per unit wavelength.
Total emissive power is :
E= E d W/ m 2
0
c. Emission from real surface-emissivity : The emissive power from a
real surface is given by
E = At4 W
Where, E = Emissivity of the material.
d. Intensity of radiation.
e. Radiation density and pressure.
f. Radiosity (J) : It refers to all of the radiant energy leaving a surface.
Que 4.2. What do you mean by black body ? Also write the
properties of black body.
Answer
A. Black Body :
1. The body that absorbs all the radiant energy reaching its surface is
known as black body.
2. For a black body : = 1, = 0, = 0
B. Properties of Black Body :
1. It absorbs all the incident radiation falling on it and does not transmit or
reflect regardless of wavelength and direction.
2. It emits maximum amount of thermal radiations at all wavelengths at
any specified temperature.
3. It is a diffuse emitter (i.e., the radiation emitted by a black body is
independent of direction).
Answer
a. White Body :
1. A body which reflects all incident radiation falling on it is known as
white body.
2. For a white body : = 1, = 0 and = 0.
b. Gray Body :
1. If the radiative properties , , of a body are assumed to be uniform
over the entire wavelength spectrum, then such a body is called gray
body.
Thermal Radiation 4–4 B (ME-Sem-5)
Answer
i. Emissivity () : It is defined as the ability of a body to radiate heat. It is
also defined as the ratio of the emissive power of a body to the emissive
power of a black body of equal temperature.
ii. Absorptivity () : It is the ratio of amount of radiation absorbed to the
total incident radiation.
= Ga / G
iii. Reflectivity () : It is the ratio of amount of radiation reflected to the
total incident radiation.
= Gr / G
iv. Transmissivity () : It is the ratio of amount of radiation transmitted
to the total incident radiation.
= Gt / G
Answer
Given : = 2 or =, Q = 660 W/m2, = 3 or = , A = 5 m2
2 3
To Find : Energy absorbed, transmitted and reflected.
1. We know that, + + = 1,
Then, + + =1
3 2
6 3 2
=1
6
11 = 6
Heat & Mass Transfer 4–5 B (ME-Sem-5)
6
=
11
6 3
So, = or
22 11
6 2
and = or
33 11
2. Energy absorbed per unit area,
6
= 660 × = 660 × = 360 W/m2
11
Total energy absorbed = 360 × 5 = 1800 W
3. Energy reflected per unit area,
2
= 660 × = 660 × = 120 W/m2
11
Total energy reflected = 120 × 5 = 600 W
4. Energy transmitted per unit area,
3
= 660 × = 660 × = 180 W/m2
11
Total energy transmitted = 180 × 5 = 900 W
Answer
a. Kirchhoff’s Law :
1. Kirchhoff’s law states that the ratio of total emissive power (E) to
absorptivity () is constant for all bodies which are in thermal equilibrium
with their environment.
Walls having uniform temperature
Large body
(2)
Eb A1 Hollow space
Small
body (1)
Fig. 4.6.1.
2. At equilibrium, the energy absorbed by body is equal to energy emitted
by body.
For body 1 : A1E1 = 1A1Eb ...(4.6.1)
Thermal Radiation 4–6 B (ME-Sem-5)
d ( E ) b
=0
d
d ( E ) b d C1 5
or
d
= 0
d C
exp 2 1
t
C2 6 5 C2 C2 1
exp t 1 (5 C1 ) C1 exp t t 2
or =0
2
C2
exp t 1
C 1 C
or – 5 C1–6 exp 2 + 5 C1–6 + C1C2–5 2 exp 2 = 0
t t t
4. Dividing both sides by 5 C1–6, we get
C 1 1 C
exp 2 1 C2 exp 2 = 0 ...(4.6.5)
t 5 t t
5. Solving eq. (4.6.5) by trial and error method, we get
C2 C2
= = 4.965
t max t
C2 1.439 104
max t = 2898 mK
4.965 4.965
max t = 2898 mK
Answer
A. Intensity of Radiation : It is defined as the rate of energy leaving a
surface in a given direction per unit solid angle per unit area of emitting
surface normal to the mean direction in space.
B. Proof :
1. Let us consider radiation from the elementary area dA 1 at the centre
of a sphere. Suppose this radiation is absorbed by a second elemental
area dA 2, a portion of the hemispherical surface.
2. The projected area of dA 1 on a plane perpendicular to the line joining
dA 1 and dA 2 = dA 1 cos .
dA
3. The solid angle subtended by dA 2 = 2 2
r
dQ12
The intensity of radiation, I = ...(4.7.1)
dA
dA1 cos 2 2
r
Where, dQ12 is the rate of radiation heat transfer from dA 1 to dA 2.
Thermal Radiation 4–8 B (ME-Sem-5)
2
dA2 = r sin d d
r si
n
rd
r
d
r sin d
dA1
d
Fig. 4.7.1.
4. Since, dA 2 = rd (r sin d)
or dA 2 = r2 sin d d ...(4.7.2)
5. From eq. (4.7.1) and eq. (4.7.2), we get
dQ12 = I dA1 sin cos d d
6. The total radiation through the hemisphere is given by
2
Q = IdA1 2
sin cos d d
0 0
2
= 2IdA1 sin cos d
0
2
= IdA1 sin 2 d
0
or, Q = I dA1
7. Since, Q = E dA 1
EdA1 = I dA 1
or, E = I
E
I=
The above equation shows that intensity of radiation is 1/ times of the
total emissive power.
PART-2
Radiation Exchange between Diffuse Non-Black Bodies in an
Enclosure, Radiation Shield, Radiation Combined with
Conduction and Convection, Absorption and Emission in
Gaseous Medium, Solar Radiation, Green House Effect.
Heat & Mass Transfer 4–9 B (ME-Sem-5)
Questions-Answers
Answer
A. Different Cases of Radiation Heat Exchange between
Non-Black Bodies :
1. Infinite parallel plates.
2. Infinite long concentric cylinders.
3. Small gray bodies.
4. Small body in a large enclosure.
B. Small Body in a Large Enclosure :
1. Consider a small body placed in a large enclosure.
2. The large gray enclosure acts like a black body, absorbing practically
all the radiation incident upon it and reflecting negligibly small energy
back to the small gray body.
3. In this case F1 – 2 = 1, since all the radiations emitted by the small body
would be intercepted by the outer large enclosure.
4. Thus, energy emitted by the enclosure 1 = A11 t14
Energy emitted by the enclosure 2 = A22 t24
Energy incident upon the small body 1 = F2–1 A22 t24
Thermal Radiation 4–10 B (ME-Sem-5)
Que 4.9. What do you mean by shape factor ? Write its salient
Answer
A. Shape Factor : It is defined as the fraction of radiative energy that is
diffused from one surface and strikes the other surface directly with
no intervening reflections.
B. Salient Features of Shape Factor :
1. When heat transfer takes place between two bodies, the shape factor
relation is given by,
A1F1–2 = A2F2–1
This is called reciprocity theorem.
2. For flat and convex surfaces, the shape factor with respect to itself is
zero because the energy leaving the surface never returns on the
same surface
i.e., F1–1 = F2–2 = F3–3 = ....Fii = 0.
3. For concave surfaces, the shape factor with respect to itself is never
zero.
4. The shape factor is purely a function of geometrical parameters.
C. Proof of Reciprocity Theorem :
1. Considering heat exchange between elementary areas dA 1 and dA 2 of
two black radiating bodies, separated by a non absorbing medium and
having areas A1 and A2 and temperatures t1 and t2 respectively.
2. Let d1 be angle subtended at dA 1 by dA 2 and d2 angle subtended at
dA 2 by dA 1. Then
dA2 cos 2
d1 =
r2
dA1 cos 1
and d2 =
r2
Heat & Mass Transfer 4–11 B (ME-Sem-5)
A 2,t2
dA 2
Normal to dA 2
r
2
1 Normal to dA 1
A1, t1
dA 1
Answer
A. Diffuse Emitter :
1. A surface whose emittance is same for all directions is called a diffuse
emitter.
Heat & Mass Transfer 4–13 B (ME-Sem-5)
x
dA
Fig. 4.10.1. Diffuse emitter.
B. Radiation Shape Factor : Refer Q. 4.9, Page 4–10B, Unit-4.
Que 4.11. Consider two concentric cylinders having diameters
10 cm and 20 cm and a length of 20 cm. Designating the open ends of
the cylinders as surfaces 3 and 4, estimate the shape factor, F3–4.
AKTU 2014-15, Marks 10
Answer
A1F1-3 = A3F3-1
(10)(20)
F3-1 = 0.0874 = 0.233
(20 2 102 ) / 4
9. Then from eq. (4.11.1),
F3-4 = 1 – 0.233 – 0.6901 = 0.0769
Answer
1. Since all the radiation being emitted by the small sphere is incident upon
and absorbed by the inner surface of the large sphere, therefore,
configuration factor between 1 and 2 is F1–2 = 1.
2. Now, A1 F1–2 = A2 F2–1 (Reciprocity theorem)
4r12 × 1 = 4 r22 × F2–1
2 2
r 0.03
F2–1 = 1 = 0.0278
r2 0.18
Thus 2.78 % of the heat emitted from the inner surface of the large
sphere is incident upon the small sphere and absorbed by it.
3. Also from energy balance for the large sphere,
F2–1 + F2–2 = 1
F2–2 = 1– F2–1 = 1 – 0.0278 = 0.9722
Thus, 97.22 % of emission from the large sphere is absorbed by the
inner surface of the sphere itself.
4. Net interchange of heat between the two spheres is,
Qnet = F1–2 A1 (t14 – t24)
= 1 × 4 × (0.03)2 × 5.67 × 10–8 × [(573)4 – (288)4]
Qnet = 64.72 W
Heat & Mass Transfer 4–15 B (ME-Sem-5)
Answer
A. Radiation Shield : The shield that reduces the radiation heat transfer
by effectively increasing the surface resistances without actually
removing any heat from the system are known as radiation shield.
B. Derivation :
1. Let us consider two parallel plates, 1 and 2, each of area A(A1 = A2 = A)
at temperature t1 and t2 respectively with a radiation shield placed between
them as shown in Fig. 4.13.1.
2. With no radiation shield, the net heat exchange between the parallel
plates is given by,
A (t14 t24 )
(Q12)net = ...(4.13.1)
1 1
1
1 2
Radiation
1 shield 2
1 3 3 2
(Q13)net (Q32)net
1 – 1 1 1 – 3 1 – 3 1 1 – 2
A11 A1F1–3 A13 A33 A3F3–2 A22
5. Since the radiation shield does not deliver or remove heat from the
system,
(Q13)net = (Q32)net
A (t14 t34 ) A (t34 t24 )
= ...(4.13.2)
1 1 1 1
1 1
1 3 3 2
6. On simplifying eq. (4.13.2), we get
1 1 1 1
t14 1 t24 1
t34 = 3 2 1 3
...(4.13.3)
1 1 1 1
1 1
3 2 1 3
7. On putting value of t3 in the left side of eq. (4.13.2), we have
A (t14 t24 )
(Q12)net = ...(4.13.4)
1 1 1 1
1 1
1 3 3 2
8. Dividing eq. (4.13.4) by eq. (4.13.1), we get
1 1
1
(Q12 )net with shield = 1 2
...(4.13.5)
(Q12 )net without shield 1 1 1 1 1 1
3 2
1 3
9. If 1 = 2 = 3, then the eq. (4.13.5) reduces to
1
(Q32)net = (Q12)net
2
Que 4.14. What do you understand by “Radiation Shield” ? If n
numbers of radiation shields are used between two parallel plates,
determine the relation for heat transfer using n shields in terms of
original heat transfer (without shields). Take emissivities of all
surfaces to be . AKTU 2013-14, Marks 10
Answer
A. Radiation Shield : Refer Q. 4.13, Page 4–15B, Unit-4.
B. Derivation :
1. In the general case where there are n shields, all the surface resistances
would be the same, since the emissivities are equal.
2. There will be two surface resistances for each shield and one for each
heat transfer surface.
3. There will also be (n + 1) ‘space resistances’ but the configuration factor
is unity for each infinite parallel plane.
1 –
Total resistance (R)n–shields = (2n 2) (n 1) A
Heat & Mass Transfer 4–17 B (ME-Sem-5)
2
= (n 1) 1 A ...(4.14.1)
4. The radiant heat transfer rate between two infinitely large parallel
plates separated by n-shields is,
1
(Q) n-shields = A (t14 – t24) ...(4.14.2)
2
(n 1) – 1
5. For no shields (n = 0), the resistance is given by,
2
(R)without shield = – 1 A ...(4.14.3)
A (t14 – t24 )
and thus, (Q)without shields = ...(4.14.4)
(2 / – 1)
6. From eq. (4.14.1) to eq. (4.14.4), we get
(Q)nshields ( R) without shields 1
= =
(Q) without shields ( R)n-shields n1
Answer
25
Given : r1 = 12.5 mm = 0.0125 m,
2
300
r3 = = 150 mm = 0.15 m, t1 = 660 + 273 = 933 K,
2
t2 = 27 + 273 = 300 K, 1 = 0.8, 3 (shield) = 0.2.
To Find : i. Reduction in heat transfer.
ii. Temperature of shield, t3.
1. Considering L is the length of the heater, the heat lost by the heater to
the room is given by,
Q = A1 1 [t14 – t24] ...(4.15.1)
Q = 2 r1 L × 0.8 × 5.67 × 10 – 8 × [(933)4 – (300)4]
= 2 × 0.0125 × L × 0.8 × 5.67 × (9334 – 3004) × 10 – 8
= 0.356 × L × (7577.5 – 81)
Q
or, q= = 0.356 × (7577.5 – 81) = 2668.7 W 2.67 kW/m
L
Thermal Radiation 4–18 B (ME-Sem-5)
2. When the cylinder is enclosed in a radiation shield then the heat flow is
given by
A1 (t14 t34 )
Q = A3 3 ( t34 t24 )
1 1 r
1 1
1 3 r3
3. As heat lost by heater to shield is further lost by shield to the room,
where suffix ‘3’ belongs to shield.
Q 2r1 (t14 t34 ) 2r3 (t34 t24 )
q = ...(4.15.2)
L 1 1 r 1
1 1
1 3 r3 3
r1 (t14 t34 ) r (t 4 t24 )
or = 3 3 ...(4.15.3)
1 1 r 1
1 1
1 3 r3 3
4. Substituting the given values in eq. (4.15.3), we get
0.0125 (9334 t34 ) 0.15 (t34 3004 )
=
1 1 0.0125 1
1
0.8 0.2 0.15 0.2
0.0125 (9334 t34 )
= 0.03 (t34 – 3004)
1.58
9334 – t34 = 3.792 (t34 – 3004) = 3.792 t34 – 3.792 ×3004
4.792 t34 = 9334 + 3.792 × 3004 = 7.885 × 1011
t34 = 1.64 × 1011 or t3 = 636.4 K or 363.4 °C
5. Substituting the value of t3 in eq. (4.15.2), we get
2 0.15 5.67 10 8 [(636.4)4 (300)4 ]
q =
1
0.2
= 1.0688 × (1640.3 – 81) = 1666.6 W/m = 1.67 kW/m
6. Percentage reduction in heat flow
q q 2.67 1.67
q = 100 × 100 = 37.45 %
q 2.67
Answer
Given : p = 0.8, = 0.05
To Find : Number of shields.
Heat & Mass Transfer 4–19 B (ME-Sem-5)
p p
Fig. 4.16.1.
1. Let, there are n shields between two parallel planes.
2. Then, increase in surface resistance = 2n and increase in space
resistance = n + 1.
3. Total resistance with n shields,
1 P 1 1 1 P
(Rth)n = 2n (n 1)
P F P
1 P 1 n 1
(Rth)n = 2 2n
P F
4. On putting the given values,
1 0.8 1 0.05 n 1
(Rth)n = 2 2n
0.8 0.05 1
( For large parallel planes; F = 1)
= 0.5 + 38n + n + 1
(Rth)n = 1.5 + 39n
5. Heat exchange per unit area with shields between two surfaces,
(t14 t24 )
(q)n =
( Rth ) n
(t14 t24 )
(q)n = ...(4.16.1)
1.5 39n
6. When, there is no shield between surfaces, heat transfer per unit
area,
( t14 t24 )
q=
1 P 1
2
P F
( t14 t24 ) (t14 t24 )
q=
1 0.8 1.5
2 1
0.8
7. Since, it is given that
1
(q)n = ×q
75
(t14 t24 ) 1 (t14 t24 )
=
1.5 3.9 n 75 1.5
Thermal Radiation 4–20 B (ME-Sem-5)
Que 4.17. Consider two large parallel plates one at t1 = 727 °C with
emis sivity 1 = 0.8 and other at 227 °C with emiss ivity
2 = 0.4. An aluminium radiation shield with an emissivity,
s = 0.05 on both sides is placed between the plates. Calculate the
percentage reduction in heat transfer rate between the two plates
as a result of shield. Use = 5.67 × 10–8 W/m2K4.
AKTU 2017-18, Marks 10
Answer
Given : t1 = 727 °C = 727 + 273 = 1000 K, 1 = 0.8,
t2 = 227 °C = 227 + 273 = 500 K, 2 = 0.4,
s = 3 = 0.05, = 5.67 × 10–8 W/m2 K4
To Find : Percentage reduction in heat transfer.
1. Without shield,
Q12 (per unit area)
1000 4 500 4
5.67
(t14 t24 ) 100 100 53156 = 19329 W
= =
1 1 1 1 2.75
1 1
1 2 0.8 0.4
2. With shield,
(Q13)net = (Q32 ) net
A ( t14 t34 ) A ( t34 t24 )
=
1 1 1 1
1 1
1 3 3 2
4 4 4 4
1000 t3 t3 500
100 100 = 100 100
1 1 1 1
1 1
0.8 0.05 0.05 0.4
10000 x 4 x 4 625 t3
= Assuming, x
1.25 20 1 20 2.5 1 100
20.25 4
10000 – x4 = ( x 625)
21.5
10000 – x4 = 0.942 (x4 – 625)
= 0.942 x4 – 588.75
1.942 x4 = 10588.75 or x = 8.59
T3 = 100 × 8.59 = 859 K
Heat & Mass Transfer 4–21 B (ME-Sem-5)
Radiation shield
Plate Plate
1 3 3 2
1 3 3 2
1000 K 500 K
Fig. 4.17.1.
1000 4 859 4
5.67
3. (Q13)net (per unit area) = 100 100
1 1
1
0.8 0.05
4555.3
= × 5.67 = 1273.2 W
20.25
4. Reduction in heat flow due to shield,
(Q12)net – (Q13)net = 19329 – 1273.2 = 18055.8 W
18055.8
5. Percentage reduction in heat transfer = = 93.41 %
19329
Que 4.18. Derive the expression for radiation combined with
convection.
Answer
1. Total heat transfer by both radiation and convection is given by,
q = qconv. + qrad.
= hconv. (tg – tw) + hrad. (tg – tw)
= (hconv. + hrad.) (tg – tw)
Where, tg = Gas temperature,
tw = Wall temperature,
hrad. = Radiation heat transfer coefficient, and
hconv. = Convection heat transfer coefficient.
2. The radiation heat transfer coefficient (h rad) is defined as,
(tg – tw) hrad. = [tg4 – tw4]
(t g4 tw4 )
or hrad. = = (tg2 + tw2) (tg + tw)
(t g – tw )
3. Now, the hrad. value for the case of two large parallel plates is given by
Q ( t14 t24 )
= = hrad. (t1 – t2)
A 1 1
–1
1 2
Thermal Radiation 4–22 B (ME-Sem-5)
( t12 t2 2 ) (t1 t2 )
hrad. =
1 1
–1
1 2
Answer
1. Consider a beam of monochromatic radiation with intensity I0 entering
a gas layer of thickness L as shown in Fig. 4.19.1.
Gas layer
I0 IL
x
dx
L
Fig. 4.19.1. Monochromatic radiation
passing through an absorbing gas.
2. As the beam passes through the gas layer, its intensity gets reduced and
the decrease is given by
dIx = – k Ixdx ...(4.19.1)
Where, Ix = Monochromatic intensity at a distance x, and
k = Monochromatic absorption coefficient,
3. Integrating eq. (4.19.1) between the limits x = 0 and x = L, we get
x L dI x L
x
x0 Ix = x 0 – k dx
I L
ln = – kL
Io
IL = Io e–kL ...(4.19.2)
Where, IL = Radiation intensity at x = L.
4. The ratio IL / Io is the monochromatic transmittivity of the gas.
5. In general, the gases do not reflect radiant energy i.e., their reflectivity
is zero, therefore,
+ = 1 or, = 1 – = 1 – e–k L
The quantity (1 – e–k L) represents the monochromatic absorptivity of
the gas.
Que 4.20. Write a short note on solar radiation.
Heat & Mass Transfer 4–23 B (ME-Sem-5)
Answer
1. Solar radiation consists of energy that is emitted by various layers of the
sun.
2. It is a function of geometrical relation between the surface and the sun
which is continuously changing on a daily and annual basis.
3. Solar radiation reaching the earth surface is partly reflected, some of
the energy is absorbed and the remainder being shared thermally to the
land masses and oceans, chemically in vegetation through the process
of photosynthesis and mechanically in the form of wind.
4. The approximate distribution of the flow of sun’s energy to the earth’s
surface is :
i. 9 % is scattered.
ii. 15 % is absorbed in the atmosphere and out of it 4 % reaches the earth’s
surface by convection.
iii. 43 % is transmitted to the earth directly and by diffuse radiation.
iv. 33 % is reflected back to space.
Answer
A. Green House Effect :
1. The green house effect is a process by which radiative energy leaving a
planetary surface is absorbed by atmospheric gases like carbon dioxide,
methane, water vapours, etc., called green house gases.
2. They transfer this energy to other components of the atmosphere and it
is re-radiated in all directions.
3. This transfer of energy to the surface increases the temperature of
environment and leads to global warming.
B. Reason to Consider it as a Matter of Great Concern : It is matter
of great concern among environmental scientists because of the following
reasons :
1. If the emission of green house gases continues, the global temperature
will increase.
2. It results in floods.
3. Global warming will lead to changes in the rainfall pattern.
Heat & Mass Transfer 5–1 B (ME-Sem-5)
5 Heat Exchanger
Questions-Answers
Answer
A. Heat Exchanger :
1. It is a device in which heat transfer takes place between two fluids at
different temperature i.e., hot fluid and cold fluid.
2. Example : Boilers, condensers, intercoolers, preheaters, regenerators,
automobile radiators, etc.
Heat & Mass Transfer 5–3 B (ME-Sem-5)
Non-condensable
gas
Hot water
Container
Steam
Cold water
Fig. 5.1.1. Heat exchanger.
B. Classification of Heat Exchanger : The classification of heat
exchanger is as follows :
a. Classification according to Nature of Heat Exchange :
i. Direct Contact Heat Exchangers :
1. In direct contact heat exchanger, the heat exchange takes place by
direct mixing of hot and cold fluids.
2. Example : Cooling tower, jet condenser, etc.
ii. Indirect Contact Heat Exchangers :
1. In indirect contact heat exchanger, the heat exchange takes place
by indirect contact of hot and cold fluids i.e., there is a wall between
two fluids.
2. Example : Gas turbine, IC engines, etc.
b. Classification according to Relative Direction of Fluid Motion :
i. Parallel Flow Heat Exchangers :
1. In a parallel flow heat exchanger, the hot and cold fluid flow in the
same directions.
2. Example : Water heaters, oil coolers, etc.
t
tc1 th1
(temperature)
Hot fluid
Cold
th1 t h2
Hot th2
Cold tc2
Cold fluid
tc1 L (length)
tc2
Cold
th1 Hot t h2
Cold
t
t c1
th1
Hot fluid
tc2
th2
L
Fig. 5.1.3. Counter flow heat exchanger.
Hot fluid
(in) Hot fluid
(Unmixed (out)
stream)
Hot
fluid
(in)
Baffles
Tube fluid
ii. Condenser :
1. In condenser, the hot fluid remains at constant temperature while
the temperature of cold fluid increases from inlet to outlet.
t
t
t h1 = t h2
th1 th2 th1
th = constant = t h– t c
tc2
tc = constant th2
= th – tc
tc1 tc1 tc2
tc1 = tc2
L L
1 2 1 2
(a) Temperature distribution (b) Temperature distribution
in a condenser. in an evaporator.
Fig. 5.1.6.
Answer
A. Fouling :
1. The phenomenon of rust formation and deposition of fluid impurities
is called fouling.
2. Due to these surface deposits the thermal resistance is increased and
eventually the performance of the heat exchanger lowers.
B. Fouling Factor :
1. Since it is difficult to ascertain the thickness and thermal conductivity
of the scale deposits, the effect of scale on heat flow is considered by
specifying an equivalent scale heat transfer coefficient hs.
2. The reciprocal of scale heat transfer coefficient, hs is called the fouling
factor, Rf .
1
Thus Rf =
hs
3. Fouling factors are determined experimentally by testing the heat
exchanger in both the clean and dirty conditions. The fouling factor, Rf
is thus defined as :
1 1 1
Rf =
h
s U dirty U clean
4. If hsi and hso be the heat transfer coefficients for the scale deposited on
the inside and outside surfaces respectively, then the thermal resistance
Heat & Mass Transfer 5–7 B (ME-Sem-5)
to scale formation on the inside surface (Rsi) and outside surface (Rso)
are given by
1
Rsi =
Ai hsi
1
Rso =
A0 hso
Answer
A. Overall Heat Transfer Coefficient : It is defined as a quantity which
gives the heat transmitted per unit area per unit time per degree
temperature difference between the bulk fluids on each side of the
substance.
Q = UAT
Where, U = Overall heat transfer coefficient.
B. Expression for Plane Wall and Cylindrical Surface :
1. The overall heat transfer coefficient for two fluids separated by a plane
wall is given by :
1
U=
1 l 1
hi k h0
2. The overall heat transfer coefficient for fluids which are separated by a
tube wall is given by :
i. For inner surface :
1
Ui =
1 ri r0 ri 1
ln
hi k ri r0 h0
ii. For outer surface :
1
U0 =
r0 1 r0 r0 1
r h k ln r h
i c i 0
Answer
A. Assumptions for LMTD Expression :
1. The overall heat transfer coefficient (U) is constant.
2. Steady flow condition.
3. Constant specific heat and mass flow rates of both fluids.
4. Heat exchanger should be perfectly insulated.
5. Change in potential and kinetic energies are negligible.
6. Axial conduction along the tubes of the heat exchanger is negligible.
B. Expression for LMTD for Parallel Flow Heat Exchanger :
Annulus surrounding
the pipe
Cold Pipe
Hot
Cold
(a) Flow arrangement.
th1 Ho
t fl
uid dth
Temperature
th2
1 = t h1 – t c1 dQ 2 = th2 – tc2
tc2
uid
d fl
C ol
dtc
dA
tc1
1 Area 2
(b) Temperature distribution
Fig. 5.4.1. Calculation of LMTD for a parallel flow heat exchanger.
1 1
dth – dtc = – dQ
C
h C c
1 1
d = – dQ ...(5.4.2)
C
h C c
1 1
d = – UdA (th – tc)
Ch Cc
1 1
or d = – UdA
C
h C c
d 1 1
or = – UdA ...(5.4.3)
C
h C c
4. Integrating the eq. (5.4.3) between inlet and outlet conditions, we get
2
d 1 1 A A
= UdA
Ch Cc
A 0
1
1 1
or ln 2 = UA ...(5.4.4)
1 C
h C c
5. Now, the total heat transfer rate between the two fluids is given by,
Q = Ch (th1 – th2) = Cc (tc2 – tc1)
1 t t
or = h1 h 2 ...(5.4.5)
Ch Q
1 t t
= c2 c1 ...(5.4.6)
Cc Q
1 1
6. Substituting the value of and from eq. (5.4.5) and eq. (5.4.6) in
Ch Cc
eq. (5.4.4), we get
t t t t
ln 2 = – UA h1 h 2 c2 c1
1 Q Q
UA UA
= [( th 2 tc 2 ) (th1 tc1 )] = (2 1 )
Q Q
Heat Exchanger 5–10 B (ME-Sem-5)
UA (2 1 )
Q=
ln 2
1
or Q = UA m
Where, m = Logarithmic mean temperature difference
(LMTD).
1 2
m = 2 1
2
ln ln 1
1 2
Que 5.5. Derive expression for LMTD for counter flow heat
exchanger.
Answer
Annulus surrounding
the pipe
Cold Pipe
Hot
Cold
th1 Ho
tf
lu i
d
1
th
tc2 C dth
o ld
flu
Temperature
id th2
dQ
tc 2
tc1
dA
dtc
Area
(b) Temperature distribution.
Fig. 5.5.1. Calculation of LMTD for a counter flow heat exchanger.
1. Let us consider an elementary area dA of the heat exchanger. The rate
of flow of heat through this elementary area dA is given by,
dQ = UdA (th – tc) = UdAt ...(5.5.1)
2. The energy balance over dA is,
dQ = – mh C ph dth = – mc C pc dtc ...(5.5.2)
Heat & Mass Transfer 5–11 B (ME-Sem-5)
dQ dQ
and, dtc =
mc C pc Cc
1 1
dth – dtc = – dQ
C
h Cc
1 1
or d = – dQ ...(5.5.3)
C
h C c
4. Putting value of dQ from eq. (5.5.1) in eq. (5.5.3),
1 1 1 1
d = – UdA (th – tc) = – UdA
C
h Cc C h Cc
d 1 1
or, = – UdA ...(5.5.4)
C
h Cc
5. On integrating eq. (5.5.4) from A = 0 to A = A, we have
1 1
ln 2 = – UA ...(5.5.5)
1 C
h C c
6. Now, the total heat transfer rate between the two fluids is given by,
Q = Ch (th1 – th2) = Cc (tc2 – tc1)
1 t t
or = h1 h 2
Ch Q
1 t t
or = c2 c1
Cc Q
1 1
7. Substituting the values of and in eq. (5.5.5), we get
Ch Cc
t t t t
ln 2 = – UA h1 h 2 c2 c1
1 Q Q
UA
= [(th1 tc2 ) (th 2 tc1 )]
Q
UA UA
= (1 2 ) (2 1 )
Q Q
UA (2 1 )
or Q=
ln 2
1
or Q = UA m
Where, m = Logarithmic mean temperature difference
(LMTD).
Heat Exchanger 5–12 B (ME-Sem-5)
2 1 2
m = 1
2
ln ln 1
1 2
Answer
1. In the thermal analysis of various types of heat exchangers LMTD
method has been used.
2. This method is pretty simple and can be used in the design of heat
exchangers when all the terminal temperatures are known or easily
determined.
3. The difficulty arises if the temperatures of the fluids entering or leaving
the heat exchanger are not known.
4. This type of situation is encountered in the selection of a heat exchanger
or when the exchanger is to be run at off design conditions.
5. In such cases, it is preferable to utilise an altogether different method
known as the effectiveness-NTU method.
Answer
A. Significance of NTU Method :
1. A heat exchanger can be designed by the LMTD, when inlet and outlet
conditions are specified.
2. When the problem is to determine the inlet or exit temperature for a
particular heat exchanger, the analysis is performed easily by using a
method based on effectiveness of the heat exchanger (concept first
proposed by Nusselt) and number of transfer units (NTU).
B. Effectiveness for Counter Flow Heat Exchanger :
1. The heat exchange dQ through an area dA of the heat exchanger is
given by
dQ = UdA (th – tc) ...(5.7.1)
= – mC ph dth = – mC pc dtc
= – Ch dth = – Cc dtc ...(5.7.2)
Heat & Mass Transfer 5–13 B (ME-Sem-5)
t t 1 1
ln h2 c1 = UA
t t
h1 c2 C
c Ch
t t UA C
ln h2 c1 = 1 c
th1 tc2 Cc Ch
th2 tc1
= exp [(UA/Cc) {1 – (Cc/Ch)}] ...(5.7.5)
th1 tc2
5. The expression for effectiveness is :
Ch (th1 th2 ) Cc (tc2 tc1 )
= =
Cmin (th1 tc1 ) Cmin (th1 tc1 )
Cmin (th1 tc1 )
Hence, th2 = th1 – ...(5.7.6)
Ch
Cmin (th1 tc1 )
tc2 = tc1 + ...(5.7.7)
Cc
6. Substituting values from eq. (5.7.6) and eq. (5.7.7) in eq. (5.7.5) we get,
Answer
1. The heat exchange dQ through an area dA is given by,
dQ = UdA (th – tc) ...(5.8.1)
Heat & Mass Transfer 5–15 B (ME-Sem-5)
= – m Cph dth = m Cpc dtc
= – Ch dth = Ccdtc ...(5.8.2)
2. From eq. (5.8.2), we get
dQ dQ
dth = and dtc = ...(5.8.3)
Ch Cc
1 1
d(th – tc) = – dQ ...(5.8.4)
Ch Cc
3. Putting the value of dQ from eq. (5.8.1) in eq. (5.8.4), we have
d(th tc ) 1 1
= – UdA ...(5.8.5)
(th tc ) Ch Cc
4. On integrating the eq. (5.8.5), we get
(th2 tc2 ) 1 1
ln = – UA
(th1 tc1 ) Ch Cc
(th2 tc2 ) UA Ch
ln
( t t )
= 1
h1 c1 Ch Cc
(th2 tc2 ) UA C
h
( t t )
= exp 1 ...(5.8.6)
h1 c1 Ch Cc
5. The expression of effectiveness is :
Ch (th1 th2 ) Cc (tc2 tc1 )
= =
Cmin ( th1 tc1 ) Cmin ( th1 tc1 )
Cmin (th1 tc1 )
So, th2 = th1 – ...(5.8.7)
Ch
Cmin (th1 tc1 )
tc2 = tc1 + ...(5.8.8)
Cc
6. On putting value of th2 and tc2 from eq. (5.8.7) and (5.8.8) in eq. (5.8.6),
we get
1 1 1
(th1 tc1 ) Cmin (th1 tc1 )
(th1 tc1 ) C
h Cc
UA C
= exp 1 h
C
h Cc
1 1 UA Ch
or 1 – Cmin = exp 1
C
h Cc C
h Cc
UA Ch
1 exp 1 C
Ch c
or = ...(5.8.9)
1 1
Cmin
Ch Cc
Heat Exchanger 5–16 B (ME-Sem-5)
7. If Cc > Ch, then Cmin = Ch and Cmax = Cc, hence eq. (5.8.9) will be,
UA Cmin
1 exp 1
Cmin Cmax
= ...(5.8.10)
C
1 min
Cmax
8. If Cc < Ch, then Cmin = Cc and Cmax = Ch, hence eq. (5.8.9) will be,
UA Cmax
1 exp 1 C
= Cmax min
...(5.8.11)
C
1 min
Cmax
9. On rearranging eq. (5.8.10) and eq. (5.8.11), we have
UA Cmin
1 exp 1 C
C
min max
= ...(5.8.12)
C
1 min
Cmax
Cmin UA
Since, R= and = NTU
Cmax Cmin
1 exp (NTU)(1 R)
=
1 R
Que 5.9. After a long time in service, a counter flow oil cooler is
checked to ascertain if its performance has deteriorated due to
fouling. In the test a standard oil flowing at 2.0 kg/s is cooled from
420 K to 380 K by a water supply of 1.0 kg/s at 300 K at inlet. If the
heat transfer surface is 3.33 m2 and the design value of the overall
heat transfer coefficient is 930 W/m 2-K, how much has it been
reduced by fouling ? Take Cp of oil as 2330 J/kg-K, Cp of water
4 174 J/kg-K. AKTU 2013-14, Marks 10
Answer
.
Given : mh = 2 kg/s, th1 = 420 K, th2 = 380 K, A = 3.33 m2,
.
mc = 1 kg/s, tc2 = 300 K, U = 930 W/m2-K, (Cp)h = 2330 J/kg-K,
(Cp)c = 4174 J/kg-K
To Find : Reduction in overall heat transfer coefficient due to fouling.
th
380 tc1
t(K)
tc
300
L
0
Fig. 5.9.1.
tc1 = 344.7 K
2. We know that,
U U 930 721
4. Reduction in U due to fouling = 100 = 22.5 %
U 930
Answer
650 W/m2-°C. Determine the rate of heat transfer and the length of
the tube required to achieve it. AKTU 2015-16, Marks 7.5
Answer
Given : CPo = 2100 J/kg-ºC, toi = 20 ºC, toe = 60 ºC, mo = 0.3 kg/s
d = 2 cm = 0.02 m, twi = 130 ºC, U = 650 W/m2-ºC, hfg = 2174 kJ/kg
To Find : i. Rate of heat transfer.
ii. Length of the tube.
1. We know that, heat transfer rate
Q = mo CPo (toe – toi)
= 0.3 × 2100 × (60 – 20) = 25200 W
Q = 25.2 kW
1 2
2. LMTD is given as, m =
ln 1
2
Where, 1 = twi – toi = 130 – 20 = 110 ºC
and 2 = twi – toe = 130 – 60 = 70 ºC
110 70
m = = 88.5 ºC
110
ln
70
3. We know that, Q= UAm
Q 25200
A=
U m 650 88.5
A= 0.438 m2
4. Also, A= dl
A 0.438
l=
d 2.0 10 2
l = 6.97 m
l 7.0 m
Que 5.12. A steam condenser is transferring 250 kW of thermal
energy at a condensing temperature of 65 °C. The cooling water
enters the condenser at 20 °C with a flow rate of 7500 kg/hr. Calculate
the LMTD. If overall heat transfer coefficient for the condenser
surface is 1250 W/m2 -°C. What surface area is required to handle
this load ? What error would be introduced if the arithmetic mean
temperature difference is used rather than the log mean temperature
difference ? AKTU 2016-17, Marks 10
Heat Exchanger 5–20 B (ME-Sem-5)
Answer
.
Given : Q = 250 kW, th1 = th2 = 65 °C, tc1 = 20 °C, mc = 7500 kg/hr,
U = 1250 W/m2-°C.
To Find : i. LMTD.
ii. Surface area.
iii. Error.
1. For the condensing steam, the temperature remains constant
throughout the flow passage, i.e.,
th1 = th2 = 65 °C
2. We know that energy balance on the cooling water gives,
.
Q = mc Cc (tc2 – tc1)
7500
250 = 4.186 × (tc2 – 20)
3600
Exit temperature of cooling water,
250 3600
tc2 = 20 + = 48.67 °C
7500 4.186
3. Log mean temperature difference,
1 2
m =
log e (1 / 2 )
Where, 1 = th1 – tc1 = 65 – 20 = 45 °C
2 = th2 – tc2 = 65 – 48.67 = 16.33 °C
45 16.33
m = = 28.28 °C
log e (45 / 16.33)
4. Now, we know that Q = UAm
Area required to handle the load is,
Q 250 103
A= = 7.07 m2
U m 1250 28.28
5. The arithmetic mean temperature difference is
1 2 45 16.33
= = 30.66 °C
2 2
Q 250 103
A = = 6.52 m2
U 1250 30.66
A A 7.07 6.52
6. Error = 100 × 100 = 7.77 %
A 7.07
Heat & Mass Transfer 5–21 B (ME-Sem-5)
Answer
20000
h =
Given : Cph = 3.3 kJ/kg-K, m = 5.56 kg/s, th1 = 120 °C,
3600
50000
c =
m = 13.89 kg/s, Cpc = 4.186 kJ/kg-K, tc1 = 20 °C, A = 10 m2,
3600
U = 1050 W/m2-C°
To Find : a. The effectiveness of the heat exchanger.
b. The output temperature of water and chemical.
Water
Chemical
Water
Fig. 5.13.1.
Heat Exchanger 5–22 B (ME-Sem-5)
1 0.471
= = 0.402
1.316
Ch (th1 th2 )
5. Also, =
Cmin ( th1 tc1 )
(120 th 2 )
or, 0.402 = ( Ch = Cmin)
(120 20)
or, th2 = 120 – 0.402 (120 – 20) = 79.8 °C
6. Substituting the value of th2 in eq. (5.13.1), we get
(120 – 79.8) = 3.17 (tc2 – 20)
(120 79.8)
tc2 = + 20 = 32.7 °C
3.17
Que 5.14. In a counter flow double pipe heat exchanger, water is
heated from 25 °C to 65 °C by an oil with a specific heat of
1.45 kJ/kg-K and mass flow rate of 0.9 kg/s. The oil is cooled from
230 °C to 160 °C. If the overall heat transfer coefficient is
420 W/m2-°C. Calculate the following :
i. the rate of heat transfer,
ii. the mass flow rate of water, and
iii. the surface area of the heat exchanger.
AKTU 2017-18, Marks 10
Answer
.
Given : tc1 = 25 °C, tc2 = 65 °C, Cph = 1.45 kJ/kg-K, mh = 0.9 kg/s,
th1 = 230 °C, th2 = 160 °C, U = 420 W/m2-°C
To Find : i. The rate of heat transfer.
ii. The mass flow rate of water.
iii. The surface area of the heat exchanger.
1. The rate of heat transfer,
Q = mh C ph (th1 th2 )
Q = 0.9 × 1.45 × (230 – 160) = 91.35 kJ/s
2. Heat lost by oil (hot fluid) = Heat gained by water (cold fluid)
Heat & Mass Transfer 5–23 B (ME-Sem-5)
mh C ph (th1 th2 ) = mh C pc (tc 2 tc1 )
91.35 = mc × 4.187 × (65 – 25)
91.35
Mass flow rate of water, mc = = 0.545 kg/s
4.187 (65 25)
tc2 = 65 °C Water 25 °C(tc1)
th1 = 230 °C Oil 160 °C( th2)
tc2 = 65 °C Water 25 °C(tc1)
(a) Flow arrangement
th1 = 230 °C
1 Oil
(h o
t fl
tc2 = 65 °C uid
)
W ate
r ( c ol
d fluid th2 = 160 °C
) 2
tc1 = 25 °C
Area/Length
(b) Temperature distribution.
Fig. 5.14.1. Counter-flow heat exchanger.
3. Logarithmic mean temperature difference (LMTD) is given by,
1 2 (th1 tc2 ) ( th2 tc1 )
m = =
ln (1 / 2 ) ln [(th1 tc2 ) / (th2 tc1 )]
(230 65) (160 25)
=
ln [(230 65) / (160 25)]
165 135
m = =149.5 °C
ln [(165 / 135)]
4. We know that, Q = UAm
Q 91.35 103
A= = 1.45 m2
U m 420 149.5
PART-2
Introduction to Condensation Phenomena, Heat Transfer Relations
for Laminar Film Condensation on Vertical Surface and on Outside
& Inside of a Horizontal Tube, Effect of Non-Condensable Gases,
Dropwise Condensation, Boiling Modes, Pool Boiling, Hystersis
in Boiling Curve, Forced Convection Boiling, Heat Pipes.
Heat & Mass Transfer SQ–1 B (ME-Sem-5)
1 Introduction to
Heat Transfer
(2 Marks Questions)
1.11. Describe the modes of heat transfer for heat loss from a
thermos flask.
Ans.
1. Convection heat loss from fluid to flask inner surface.
2. Conduction heat loss from inner surface of flask to outermost
surface of flask.
3. Radiation heat loss from outer surface of flask to atmosphere.
2 Marks Questions SQ–4 B (ME-Sem-5)
2 Fins and
Transient Conduction
(2 Marks Questions)
hAs
that larger the quantity , faster the exponential term will
VC
approach zero or more rapid will be the response of the
temperature measuring device.
2.10. Define sensitivity.
Ans. Sensitivity is defined as the time required by a thermocouple to
reach its 63.2 % of the value of initial temperature difference is
called its sensitivity.
2.11. Why are thin and closely spaced fins always preferred ?
Ans. Since,
Pk
fin =
hA
If the ratio of P (perimeter) and A (cross-sectional area) is increased
the effectiveness of the fin is improved. Due to this reason, thin
and closely spaced fins are preferred.
2 Marks Questions SQ–6 B (ME-Sem-5)
gtL3
Gr =
2
Heat & Mass Transfer SQ–9 B (ME-Sem-5)
4 Thermal Radiation
(2 Marks Questions)
4.6. Define all radiative properties of the surface and also explain
relation between them. AKTU 2016-17, Marks 02
Ans.
1. Following are the radiation properties of surfaces :
i. Total emissive power (E) :
Eb = T4 W/m2
Eb = AT4 W
ii. Monochromatic emissive power.
E= E d W/ m 2
0
iii. Intensity of radiation.
iv. Radiosity.
v. Irradiation.
2. Relation between total emissive power and monochromatic emissive
power is given as :
E
=
Eb
4.16. Why does a concave surface have a shape factor with itself ?
Ans. A concave surface has a shape factor with itself because the radiant
energy coming out from the one part of the surface is intercepted
by the other part of the same surface. The shape factor of a surface
with respect to itself is F1–1.
4.17. Why a flat or convex surface has a zero shape factor with
respect to itself ?
Ans. For a flat or convex surface, the shape factor with respect to itself
is zero (i.e., F1–1 = 0). This is due to the fact that for any part of flat
or convex surface, one cannot see/view any other part of the same
surface.
Heat & Mass Transfer SQ–13 B (ME-Sem-5)
5 Heat Exchanger
(2 Marks Questions)
5.11. What is the difference between pool boiling and flow boiling ?
AKTU 2016-17, Marks 02
Ans.
S. No. Pool Boiling Flow Boiling
1. The value of critical heat flux The value of critical heat flux of
of water is 1.3 MW/m2 at 1 water reaches upto 35 MW/m2
atm. at 1 atm.
2. It is easy. It is more complicated.
5.12. What are the uses of heat pipe ? AKTU 2016-17, Marks 02
Ans.
1. Spacecraft,
2. Cooking,
3. Permafrost cooling, and
4. Computer system, etc.