1D SS Heat Conduction
1D SS Heat Conduction
1D SS Heat Conduction
1
Steady Heat Conduction in Plane Walls
Considerable Temperature Difference between the inner and the outer surfaces of the wall (significant temperature gradient
in the x-direction)
Assuming heat transfer is the only in 1-dimension (1D):
Steady state conduction
𝑑 𝑑𝑇 𝑑𝑇 𝑇 𝑥 = 𝐶1 𝑥 + 𝐶2
𝑘 =0 = 𝐶1
𝑑𝑥 𝑑𝑥 𝑑𝑥 2
Steady Heat Conduction in Plane Walls
𝑇 𝑥 = 𝐶1 𝑥 + 𝐶2
𝐴𝑡 𝑥 = 0, 𝑇 = 𝑇1 𝐴𝑡 𝑥 = 𝐿, 𝑇 = 𝑇2
𝑞𝑥
𝑇2 − 𝑇1 𝑇2 − 𝑇1
𝐶2 = 𝑇1 𝐶1 = 𝑇 𝑥 = 𝑥 + 𝑇1
𝐿 𝐿
Temperature distribution across a wall is linear and is independent of thermal
conductivity
Heat Conduction Rate through the plane wall using Fourier’s Eq.
𝑑𝑇 𝑑𝑇 𝑇2 − 𝑇1 𝑑𝑇 𝑇2 − 𝑇1
𝑞𝑥 = −𝑘𝐴 = 𝐶1 = 𝑞𝑥 = −𝑘𝐴 = −𝑘𝐴
𝑑𝑥 𝑑𝑥 𝐿 𝑑𝑥 𝐿
Once the rate of heat conduction is available, temperature T(x) at any location x can be determined by replacing T2 in above Eq.
by T, and L by x
Steady Heat Conduction in Plane Walls
𝑑𝑇 𝑇2 − 𝑇1
𝑞𝑥 = −𝑘𝐴 = −𝑘𝐴
𝑑𝑥 𝐿
Heat Conduction rate through a plane wall can be rearranged as
𝑇1 − 𝑇2
𝑞𝑥 = (𝑊)
𝑅𝑤𝑎𝑙𝑙
Where Rwall is the conduction resistance
𝐿
𝑅𝑤𝑎𝑙𝑙 = (℃Τ𝑊)
𝑘𝐴
Thermal Resistance of a medium depends on the geometry and the thermal properties of the
medium
Steady Heat Conduction in Plane Walls
Analogy to Electrical Current Flow
𝐿
𝑅𝑤𝑎𝑙𝑙 = (℃Τ𝑊)
𝑘𝐴
𝑉1 − 𝑉2
Above Eq. is analogous to the relation for electric current flow I, expressed as: 𝐼 =
𝑅𝑒
𝑇1 − 𝑇2
𝑞𝑥 =
𝑅𝑤𝑎𝑙𝑙
𝑇𝑠 − 𝑇𝑠𝑢𝑟𝑟
𝑞𝑟𝑎𝑑 = ℎ𝑟𝑎𝑑 𝐴𝑠 𝑇𝑠 − 𝑇𝑠𝑢𝑟𝑟 =
𝑅𝑟𝑎𝑑
1
𝑅𝑟𝑎𝑑 = (℃Τ𝑊)
ℎ𝑟𝑎𝑑 𝐴𝑠
Radiation and Convection Resistance
A surface exposed to the surrounding might involves convection and radiation
simultaneously
𝑞𝑐𝑜𝑛𝑣
Convection and Radiation Resistances are
parallel to each other
𝑞
𝑇1 − 𝑇∞ A1 1 k1
𝑞𝑥 = 3
𝑅𝑡𝑜𝑡𝑎𝑙
T1 A3 h, T∞
A2 2 k2
k3
𝑅𝑡𝑜𝑡𝑎𝑙 = 𝑅12 + 𝑅3 + 𝑅𝑐𝑜𝑛𝑣
𝑅1 𝑅2 L1=L2 L3
= + 𝑅3 + 𝑅𝑐𝑜𝑛𝑣
𝑅1 + 𝑅2
𝑞1
𝑞𝑥 𝑞𝑥
𝐿1 𝐿2 𝐿3 1
𝑅1 = 𝑅2 = 𝑅3 = 𝑅𝑐𝑜𝑛𝑣 = 𝑞2
𝑘1 𝐴1 𝑘2 𝐴2 𝑘3 𝐴3 ℎ𝐴3
The Overall Heat-transfer Coefficient
While dealing with the problems of fluid to fluid heat transfer across a metal
boundary, it is usual to adopt an overall heat transfer coefficient U Thf Cold
fluid
film
𝑅𝑎𝑡𝑒 𝑜𝑓 𝑅𝑎𝑡𝑒 𝑜𝑓 𝑅𝑎𝑡𝑒 𝑜𝑓 T1
ℎ𝑒𝑎𝑡 𝑐𝑜𝑛𝑣𝑒𝑐𝑡𝑖𝑜𝑛 = ℎ𝑒𝑎𝑡 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑖𝑜𝑛 = ℎ𝑒𝑎𝑡 𝑐𝑜𝑛𝑣𝑒𝑐𝑡𝑖𝑜𝑛 hcf
𝑖𝑛𝑡𝑜 𝑡ℎ𝑒 𝑤𝑎𝑙𝑙 𝑡ℎ𝑟𝑜𝑢𝑔ℎ 𝑡ℎ𝑒 𝑤𝑎𝑙𝑙 𝑓𝑟𝑜𝑚 𝑡ℎ𝑒 𝑤𝑎𝑙𝑙 q q
1 1 𝐿 0.008
𝑅𝑖 = 𝑅𝑐𝑜𝑛𝑣.1 = = = 0.08333 ℃Τ𝑊 𝑅𝑔𝑙𝑎𝑠𝑠 = 𝑅𝑐𝑜𝑛𝑑 = =
ℎ1 𝐴 10 × 1.2 𝑘𝐴 0.78 × 1.2
= 0.00855 ℃Τ𝑊
1 1
𝑅0 = 𝑅𝑐𝑜𝑛𝑣.2 = = = 0.02083 ℃Τ𝑊 𝑅𝑡𝑜𝑡𝑎𝑙 = 𝑅𝑐𝑜𝑛𝑣.1 + 𝑅𝑐𝑜𝑛𝑑 + 𝑅𝑐𝑜𝑛𝑣.2
ℎ2 𝐴 40 × 1.2
= 0.08333 + 0.00855 + 0.02083
= 0.1127 ℃Τ𝑊
𝑇∞1 − 𝑇∞2 20 − −10
𝑞= = = 266𝑊
𝑅𝑡𝑜𝑡𝑎𝑙 0.1127
𝑇∞1 − 𝑇𝑖 20 − 𝑇𝑖
𝑞= = = 266 𝑊 𝑇𝑖 = 20 − 266 × 0.08333 = −2.2℃
𝑅𝑐𝑜𝑛𝑣,1 0.08333
Steady State Heat Conduction
Heat conduction in a wall with heat generation
Consider plane wall with steady state 1D uniform heat generation per unit volume, 𝒒.ሶ
𝑇2 − 𝑇1 𝑞ሶ 2 𝑇1 + 𝑇2
𝐶1 = 𝐶2 = 𝐿 +
2𝐿 2𝑘 2
𝑇 −𝐿 = 𝑇1 𝑇 𝐿 = 𝑇2
ሶ 2
𝑞𝐿 𝑥2 𝑇2 − 𝑇1 𝑥 𝑇1 + 𝑇2
𝑇 𝑥 = 1− 2 + +
2𝑘 𝐿 2 𝐿 2
ሶ 2
𝑞𝐿 𝑥2 𝑇2 − 𝑇1 𝑥 𝑇1 + 𝑇2
Steady State Heat Conduction 𝑇 𝑥 =
2𝑘
1− 2 +
𝐿 2 𝐿
+
2
No heat transfer across this plane, and it may be represented by the adiabatic surface
ሶ 2
𝑞𝐿 𝑥2
𝑇 𝑥 = 1 − 2 + 𝑇𝑠
2𝑘 𝐿
Also applies to plane walls that are perfectly insulated on one side (x =0) and maintained at a
fixed temperature Ts on the other side (x =L).
The temperature equation in terms of T0 can be written as:
𝑇 𝑥 − 𝑇𝑜 𝑥 2
Normally, it is the temperature of an adjoining fluid, T∞, and not Ts,
= which is known. It then becomes necessary to relate Ts to T∞
𝑇𝑠 − 𝑇𝑜 𝐿
Steady State Heat Conduction ሶ 2
𝑞𝐿 𝑥2
𝑇 𝑥 = 1 − 2 + 𝑇𝑠
2𝑘 𝐿
Heat conduction in a wall with heat generation
𝑇 𝑥 − 𝑇𝑜 𝑥 2
= Normally, it is the temperature of an adjoining fluid, T∞, and not Ts,
𝑇𝑠 − 𝑇∞ 𝐿 which is known. It then becomes necessary to relate Ts to T∞
𝑑𝑇
−𝑘 = ℎ 𝑇𝑠 − 𝑇∞
𝑑𝑥 𝑥=𝐿
𝑑𝑇
−𝑘 = 𝑞𝐿
ሶ
𝑑𝑥 𝑥=𝐿
𝑞𝐿ሶ
𝑇𝑠 = 𝑇∞ +
ℎ
EXAMPLE
A plane wall is a composite of two materials, A and B. The wall of material A has uniform heat generation = 1.5x106
W/m3, kA =75 W/mK, and thickness LA = 50 mm. The wall material B has no generation with kB= 150 W/mK and
thickness LB = 20 mm. The inner surface of material A is well insulated, while the outer surface of material B is cooled
by a water stream with T ∞= 30°C and h= 1000 W/m2 K.
1. Sketch the temperature distribution that exists in the composite under steady state conditions.
2. Determine the temperature T0 of the insulated surface and the temperature T2 of the cooled surface.
Assumptions:
• Steady State conditions
• 1D conduction in x-direction
T1 dT
T2
dr
Steady State Heat Conduction
Heat conduction through a Hollow Cylinder
Area through which heat is transmitted A = 2πrL Heat flows
radially
Path length = dr (over which temp fall is dT)
outwards
𝑑𝑇 𝑑𝑇 𝑑𝑟
𝑞 = −𝑘𝐴 = −𝑘 2𝜋𝑟𝐿 𝑞 = −𝑘 2𝜋𝐿 𝑑𝑇
𝑑𝑟 𝑑𝑟 𝑟 r2
𝑟2 𝑇2
𝑑𝑟 dr r1 No heat flow
𝑞න = −𝑘2𝜋𝐿 න 𝑑𝑇 in axial
𝑟 direction
𝑟1 𝑇1
r
𝑟2 𝑇2
𝑞 ln 𝑟 𝑟1 = −𝑘2𝜋𝐿 𝑇 𝑇1 T1 dT
𝑟2
𝑞 ln = −𝑘2𝜋𝐿 𝑇2 − 𝑇1 = 𝑘2𝜋𝐿 𝑇1 − 𝑇2 T2
𝑟1 dr
Steady State Heat Conduction
Heat conduction through a Hollow Cylinder
𝑟2 Heat flows
𝑞 ln = −𝑘2𝜋𝐿 𝑇2 − 𝑇1 = 𝑘2𝜋𝐿 𝑇1 − 𝑇2 radially
𝑟1 outwards
𝑘2𝜋𝐿 𝑇1 − 𝑇2 𝑇1 − 𝑇2 r2
𝑞= =
ln 𝑟2 Τ𝑟1 ln 𝑟2 Τ𝑟1
2𝜋𝑘𝐿 dr r1 No heat flow
in axial
direction
ln 𝑟2 Τ𝑟1 r
𝑅𝑐𝑦𝑙 =
2𝜋𝑘𝐿
T1 dT
ln 𝑂𝑢𝑡𝑒𝑟 𝑟𝑎𝑑𝑖𝑢𝑠Τ𝐼𝑛𝑛𝑒𝑟 𝑟𝑎𝑑𝑖𝑢𝑠
=
2𝜋 𝑡ℎ𝑒𝑟𝑚𝑎𝑙 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑖𝑣𝑖𝑡𝑦 (𝐿𝑒𝑛𝑔𝑡ℎ)
T2
dr
Steady State Heat Conduction
Temperature distribution in a Hollow Cylinder
𝑘2𝜋𝐿 𝑇1 − 𝑇2 𝑘2𝜋𝐿 𝑇 − 𝑇2
𝑞= =
ln 𝑟2 Τ𝑟1 ln 𝑟2 Τ𝑟
𝑇1 − 𝑇2 𝑇 − 𝑇2 𝑇1 − 𝑇2
= = 𝑇= × ln 𝑟2 Τ𝑟 + 𝑇2
ln 𝑟2 Τ𝑟1 ln 𝑟2 Τ𝑟 ln 𝑟2 Τ𝑟1
(ln 𝑟2 − ln 𝑟) −(ln 𝑟2 − ln 𝑟)
𝑇= 𝑇1 − 𝑇2 + 𝑇2 𝑇= 𝑇1 − 𝑇2 + 𝑇2
(ln 𝑟2 − ln 𝑟1 ) −(ln 𝑟2 − ln 𝑟1 )
(ln 𝑟 − ln 𝑟2 ) 𝑇1 − 𝑇2
𝑇= 𝑇1 − 𝑇2 + 𝑇2 𝑇= × ln 𝑟Τ𝑟2 + 𝑇2
(ln 𝑟1 − ln 𝑟2 ) ln 𝑟1 Τ𝑟2
Steady State Heat Conduction
Heat conduction through a Composite Cylinder Cold Fluid (air)
Tcf
𝑞
𝑞 = ℎℎ𝑓 . 2𝜋𝑟1 𝐿. (𝑇ℎ𝑓 − 𝑇1 ) 𝑇ℎ𝑓 − 𝑇1 = q
ℎℎ𝑓 . 2𝜋𝑟1 𝐿
B
hcf
A
𝑞
𝑘𝐴 . 2𝜋𝐿 𝑇1 − 𝑇2 𝑇1 − 𝑇2 =
= 𝑘𝐴 . 2𝜋𝐿
ln(𝑟2 Τ𝑟1 )
ln(𝑟2 Τ𝑟1 ) hhf
Thf T1 T2 T3 Tcf
𝑘𝐵 . 2𝜋𝐿 𝑇2 − 𝑇3 𝑞 Hot Fluid
= 𝑇2 − 𝑇3 =
𝑘𝐵 . 2𝜋𝐿
ln(𝑟3 Τ𝑟2 )
ln(𝑟3 Τ𝑟2 )
𝑞
= ℎ𝑐𝑓 . 2𝜋𝑟3 𝐿. (𝑇3 − 𝑇𝑐𝑓 ) 𝑇3 − 𝑇𝑐𝑓 =
ℎ𝑐𝑓 . 2𝜋𝑟3 𝐿 r1
r2
r3
Steady State Heat Conduction
Heat conduction through a Composite Cylinder
𝑞 𝑞 𝑞 𝑞
𝑇ℎ𝑓 − 𝑇1 = 𝑇1 − 𝑇2 = 𝑇2 − 𝑇3 =
ℎℎ𝑓 . 2𝜋𝑟1 𝐿 𝑘𝐴 . 2𝜋𝐿 𝑘𝐵 . 2𝜋𝐿 𝑇3 − 𝑇𝑐𝑓 = ℎ𝑐𝑓 . 2𝜋𝑟3 𝐿
ln(𝑟2 Τ𝑟1 ) ln(𝑟3 Τ𝑟2 )
𝑞 1 1 1 1 2𝜋𝐿 (𝑇ℎ𝑓 − 𝑇𝑐𝑓 )
𝑇ℎ𝑓 − 𝑇𝑐𝑓 = + + + 𝑞=
2𝜋𝐿 ℎℎ𝑓 . 𝑟1 𝑘𝐴 𝑘𝐵 ℎ𝑐𝑓 . 𝑟3 1 1 1 1
+ + +
ln(𝑟2 Τ𝑟1 ) ln(𝑟3 Τ𝑟2 ) ℎℎ𝑓 . 𝑟1 𝑘𝐴 𝑘𝐵 ℎ𝑐𝑓 . 𝑟3
ln(𝑟2 Τ𝑟1 ) ln(𝑟3 Τ𝑟2 )
2𝜋𝐿 (𝑇ℎ𝑓 − 𝑇𝑐𝑓 )
𝑞=
1 ln(𝑟2 Τ𝑟1 ) ln(𝑟3 Τ𝑟2 ) 1
+ + +
ℎℎ𝑓 . 𝑟1 𝑘𝐴 𝑘𝐵 ℎ𝑐𝑓 . 𝑟3
Consider a long cylinder with steady state 1D uniform heat generation per unit volume, with constant thermal conductivity
1𝑑 𝑑𝑇 𝑞ሶ 𝑑 𝑑𝑇 𝑞ሶ 𝑑𝑇 𝑞ሶ 𝑑𝑇 𝑞ሶ 2
𝑟 + =0 𝑟 =− 𝑟 𝑑 𝑟 = − 𝑟𝑑𝑟 𝑟 = − 𝑟 + 𝐶1
𝑟 𝑑𝑟 𝑑𝑟 𝑘 𝑑𝑟 𝑑𝑟 𝑘 𝑑𝑟 𝑘 𝑑𝑟 2𝑘
𝑑𝑇 𝑞ሶ 𝐶1 𝑞ሶ 𝐶1 𝑞ሶ 2
=− 𝑟+ 𝑑𝑇 = − 𝑟 + 𝑑𝑟 𝑇(𝑟) = − 𝑟 + 𝐶1 ln 𝑟 + 𝐶2
𝑑𝑟 2𝑘 𝑟 2𝑘 𝑟 4𝑘
Steady State Heat Conduction
Heat conduction through a Cylinder with Energy Generation
𝑞ሶ 2 𝑑𝑇
𝑇(𝑟) = − 𝑟 + 𝐶1 ln 𝑟 + 𝐶2 𝐵𝑜𝑢𝑛𝑑𝑎𝑟𝑦 𝐶𝑜𝑛𝑑𝑖𝑡𝑖𝑜𝑛𝑠: ቤ = 0 𝑎𝑛𝑑 𝑇 𝑟𝑜 = 𝑇𝑠
4𝑘 𝑑𝑟 𝑟=0
𝑑𝑇 𝑞ሶ 𝐶1
=− 𝑟+ Placing f𝑖𝑟𝑠𝑡 𝑏𝑜𝑢𝑛𝑑𝑎𝑟𝑦 𝑐𝑜𝑛𝑑𝑖𝑡𝑖𝑜𝑛 𝐶1 = 0
𝑑𝑟 2𝑘 𝑟
𝑞ሶ 2
Placing 𝑠𝑒𝑐𝑜𝑛𝑑 𝑏𝑜𝑢𝑛𝑑𝑎𝑟𝑦 𝑐𝑜𝑛𝑑𝑖𝑡𝑖𝑜𝑛 𝐶2 = 𝑟 + 𝑇𝑠
4𝑘 𝑜
ሶ 𝑜2
𝑞𝑟 𝑟2
𝑇(𝑟) = 1 − 2 + 𝑇𝑠
4𝑘 𝑟𝑜
ሶ 𝑜2
𝑞𝑟 ሶ 𝑜2
𝑞𝑟
𝐴𝑡 𝑚𝑖𝑑𝑝𝑙𝑎𝑛𝑒 𝑟 = 0; 𝑇 0 = 𝑇𝑜 = + 𝑇𝑠 𝑇𝑜 − 𝑇𝑠 =
4𝑘 4𝑘 𝑂𝑣𝑒𝑟𝑎𝑙𝑙 𝐸𝑛𝑒𝑟𝑔𝑦 𝐵𝑎𝑙𝑎𝑛𝑐𝑒:
𝑞ሶ 𝜋𝑟𝑜2 𝐿 = ℎ 2𝜋𝑟𝑜 𝐿 𝑇𝑠 − 𝑇∞
𝑟2 𝑇 𝑟 − 𝑇𝑠 𝑟2
𝑇 𝑟 − 𝑇𝑠 = 𝑇𝑜 − 𝑇𝑠 1− 2 = 1− 2
𝑟𝑜 𝑇𝑜 − 𝑇𝑠 𝑟𝑜 𝑞𝑟
ሶ 𝑜
𝑇𝑠 = + 𝑇∞
2ℎ
Spheres (Heat Transfer)
Steady State Heat Conduction
Heat conduction through a hollow sphere
Area through which the heat is transmitted, A = 4πr2
𝑟2 𝑇2
2
𝑑𝑇 𝑑𝑟 q (Heat flows Hollow Sphere
𝑞 = −𝑘. 4𝜋𝑟 . 𝑞 න 2 = −4𝜋𝑘 න 𝑑𝑇 radially
𝑑𝑟 𝑟1 𝑟 𝑇1 outwards, T1>T2
Element
𝑟2 r2
𝑟 −2+1 1 1
𝑞 = −4𝜋𝑘 𝑇
𝑇2 −𝑞 − = −4𝜋𝑘 𝑇2 − 𝑇1 r
−2 + 1 𝑇1 𝑟2 𝑟1 dr
𝑟1
r1
𝑞 𝑟2 − 𝑟1 T1 T2
= 4𝜋𝑘 𝑇1 − 𝑇2
𝑟1 𝑟2
4𝜋𝑘𝑟1 𝑟2 𝑇1 − 𝑇2
𝑞=
𝑟2 − 𝑟1
𝑇1 − 𝑇2
𝑞=
𝑟2 − 𝑟1
4𝜋𝑘𝑟1 𝑟2 𝑟2 − 𝑟1 (𝑜𝑢𝑡𝑒𝑟 𝑟𝑎𝑑𝑖𝑢𝑠 − 𝑖𝑛𝑛𝑒𝑟 𝑟𝑎𝑑𝑖𝑢𝑠)
𝑅𝑆𝑝ℎ𝑒𝑟𝑒 = =
4𝜋𝑘𝑟1 𝑟2 4𝜋 𝑜𝑢𝑡𝑒𝑟 𝑟𝑎𝑑𝑖𝑢𝑠 𝑖𝑛𝑛𝑒𝑟 𝑟𝑎𝑑𝑖𝑢𝑠 𝑡ℎ𝑒𝑟𝑚𝑎𝑙 𝑐𝑜𝑛𝑑𝑢𝑐𝑡𝑖𝑣𝑖𝑡𝑦
Steady State Heat Conduction
Temperature distribution in a hollow sphere
1 1 1 1
−𝑞 − = −4𝜋𝑘 𝑇2 − 𝑇1 𝑞 − = 4𝜋𝑘 𝑇1 − 𝑇2
𝑟2 𝑟1 𝑟1 𝑟2
𝑇1 − 𝑇2 𝑇1 − 𝑇 𝑇1 − 𝑇2 𝑇1 − 𝑇 𝑟1 𝑟2 𝑟 − 𝑟1
𝑞 = 4𝜋𝑘 = 4𝜋𝑘 = 𝑇1 − 𝑇2 × × = 𝑇1 − 𝑇
1 1 1 1 𝑟2 − 𝑟1 𝑟 − 𝑟1 𝑟2 − 𝑟1 𝑟1 𝑟
𝑟1 − 𝑟2 𝑟1 − 𝑟 𝑟1 𝑟2 𝑟1 𝑟
𝑟1 𝑟2 𝑟 − 𝑟1 𝑟2 𝑟 − 𝑟1 𝑟
𝑇 = 𝑇1 − 𝑇1 − 𝑇2 × × 𝑇 = 𝑇1 − 𝑇1 − 𝑇2 × Dividing by 𝑟2
𝑟2 − 𝑟1 𝑟1 𝑟 𝑟(𝑟2 − 𝑟1 ) 2
𝑟 − 𝑟1
𝑇 = 𝑇1 − 𝑇1 − 𝑇2 × 𝑟 𝑟
Dividing by 𝑟
𝑟2 (𝑟2 − 𝑟1 )
1 − (𝑟1 Τ𝑟)
𝑇 = 𝑇1 − 𝑇1 − 𝑇2 ×
1 − 𝑟1 Τ𝑟2
Steady State Heat Conduction
Heat conduction through a Composite Sphere
𝑞 𝑞 𝑟2 − 𝑟1 𝑞 𝑟3 − 𝑟2 𝑞
𝑇ℎ𝑓 − 𝑇1 = 𝑇1 − 𝑇2 = 𝑇2 − 𝑇3 = 𝑇3 − 𝑇𝑐𝑓 =
ℎℎ𝑓 . 4𝜋𝑟1 2 4𝜋𝑘𝐴 𝑟1 𝑟2 4𝜋𝑘𝐵 𝑟2 𝑟3 ℎ𝑐𝑓 . 4𝜋𝑟3 2
𝑞 1 𝑟2 − 𝑟1 𝑟3 − 𝑟2 1
𝑇ℎ𝑓 − 𝑇𝑐𝑓 = 2
+ + +
4𝜋 ℎℎ𝑓 . 𝑟1 𝑘𝐴 𝑟1 𝑟2 𝑘𝐵 𝑟2 𝑟3 ℎ𝑐𝑓 . 𝑟3 2
4𝜋 𝑇ℎ𝑓 − 𝑇𝑐𝑓
𝑞=
1 𝑟2 − 𝑟1 𝑟3 − 𝑟2 1
+ + +
ℎℎ𝑓 . 𝑟1 2 𝑘𝐴 𝑟1 𝑟2 𝑘𝐵 𝑟2 𝑟3 ℎ𝑐𝑓 . 𝑟3 2
T2 is not known therefore we start with an assumption for T2 and later on validate the accuracy of
our assumption.
T2 must be between 0°C and 22°C and it should be close to 0°C, since heat transfer coefficient within
the tank is considerably higher. Therefor assuming the T2 value as 5 °C
ℎ𝑟𝑎𝑑 = 1 5.67 × 10−8 278 + 295 2782 + 2952 = 5.34 𝑊 Τ𝑚2 . 𝐾
1 1 𝑟2 − 𝑟1 1.52 − 1.50
𝑅𝑐𝑜𝑛𝑣,𝑖 = = = 0.000442 ℃Τ𝑊 𝑅𝑆𝑝ℎ𝑒𝑟𝑒 = =
ℎ1 𝐴1 80 × 28.3 4𝜋𝑘𝑟1 𝑟2 4𝜋1 × 15 × 1.52 × 1.50
= 0.000047 ℃Τ𝑊
1 1
𝑅𝑐𝑜𝑛𝑣,𝑜 = = = 0.00345 ℃Τ𝑊 1 1
ℎ2 𝐴2 10 × 29 𝑅𝑟𝑎𝑑 = = = 0.00646 ℃Τ𝑊
ℎ𝑟𝑎𝑑 𝐴2 5.34 × 29
EXAMPLE Heat Transfer to a Spherical Container
𝑇∞2 − 𝑇2
𝑞= = 22 − 8029 × 0.00225 = 4℃
Verification of our assumption 𝑅𝑒𝑞𝑢𝑣
𝑤ℎ𝑖𝑐ℎ 𝑖𝑠 𝑠𝑢𝑓𝑓𝑖𝑐𝑒𝑛𝑡𝑙𝑦 𝑐𝑙𝑜𝑠𝑒 𝑠𝑜 𝑛𝑜 𝑛𝑒𝑒𝑑 𝑡𝑜 𝑟𝑒𝑝𝑒𝑎𝑡 𝑡ℎ𝑒 𝑐𝑎𝑙𝑐𝑢𝑙𝑎𝑡𝑖𝑜𝑛𝑠
673,300
𝐴𝑚𝑜𝑢𝑛𝑡 𝑜𝑓 𝑚𝑒𝑙𝑡𝑒𝑑 𝑖𝑐𝑒 = = 2079𝑘𝑔
333.7
Critical Radius for
Insulation
Steady State Heat Conduction
Critical Radius of Insulation
Adding more insulation to a wall always decreases heat transfer
Adding insulation to a cylindrical pipe or a spherical shell, however, is a different matter- WHY???
Adding insulation increases Conduction Resistance of the insulation layer but decreases the Convection
Resistance of the surface because of the increase in the outer surface area for convection
k
Heat Transfer from the pipe may increase or decrease, depending on which
r1 Rins Rconv
effect dominates
T1 T∞
A cylindrical pipe of outer radius r1 whose outer surface temperature T1 h
is maintained constant Insulation r2
o Pipe is covered with an insulator (k and r2)
o Convection heat transfer at T∞ and h
𝑇1 − 𝑇∞ 𝑇1 − 𝑇∞
Rate of heat transfer from insulated pipe to the surrounding air: 𝑞 = =
𝑅𝑖𝑛𝑠 + 𝑅𝑐𝑜𝑛𝑣 ln(𝑟2 Τ𝑟1 ) + 1
2𝜋𝐿𝑘 ℎ(2𝜋𝑟2 𝐿)
Steady State Heat Conduction
Critical Radius of Insulation
𝑇1 − 𝑇∞ 𝑇1 − 𝑇∞
Rate of heat transfer from insulated pipe to the surrounding air: 𝑞 = =
𝑅𝑖𝑛𝑠 + 𝑅𝑐𝑜𝑛𝑣 ln(𝑟2 Τ𝑟1 ) 1
+
2𝜋𝐿𝑘 ℎ(2𝜋𝑟2 𝐿)
Variation of the heat transfer rate with the outer radius of the insulation r2 is shown
it is evident that as r2 increases, the factor ln(r2/r1)/k increases but the factor 1/hr2
decreases
q
value of r2 at which q reaches a maximum is determined by:
1 1 1 1 𝑘
. = 𝑟2 = 𝑟𝑐 =
𝑘 𝑟2 ℎ 𝑟22 ℎ 0 r1 rcr=k/h r2
Steady State Heat Conduction
Critical Radius of Insulation
𝑘
𝑟𝑐𝑟,𝑐𝑦𝑙𝑖𝑛𝑑𝑒𝑟 = (𝑚) represents condition for minimum resistance and consequently maximum heat flow rate
ℎ
Insulation radius at which resistance to heat flow is minimum is called the ‘Critical Radius’
Critical Radius rc is dependent on the thermal quantities k and h and is independent of r1 (i.e., cylinder radius)
Lowest value of h encountered in practice is about 5 W/m2·°C for the case of natural convection of gases, and k
of common insulating materials is about 0.05 W/m2·°C
Critical Radius would be much less in forced convection, often less than 1 mm, because of much larger h
values associated with forced convection
EXAMPLE Heat Loss from an Insulated Electric Wire
A 3-mm-diameter and 5-m-long electric wire is tightly wrapped with a 2-mm thick plastic cover whose thermal conductivity is k
0.15 W/m · °C. Electrical measurements indicate that a current of 10 A passes through the wire and there is a voltage drop of 8 V
along the wire. If the insulated wire is exposed to a medium at T 30°C with a heat transfer coefficient of h 12 W/m2 · °C,
determine the temperature at the interface of the wire and the plastic cover in steady operation. Also determine whether
doubling the thickness of the plastic cover will increase or decrease this interface temperature.
1 1
𝑞 = 𝑉𝐼 = 8 × 10 = 80𝑊 𝑅𝑐𝑜𝑛𝑣 = = = 0.76 ℃Τ𝑊
ℎ𝐴2 12 × 0.110
q
ln(𝑟2 Τ𝑟1 ) ln(3.5Τ1.5)
𝑅𝑝𝑙𝑎𝑠𝑡𝑖𝑐 = = = 0.18 ℃Τ𝑊 𝑅𝑡𝑜𝑡𝑎𝑙 = 0.76 + 0.18 = 0.94 ℃Τ𝑊
2𝜋𝐿𝑘 2𝜋(5)(0.15)
𝑇1 − 𝑇∞
𝑞= 𝑇1 = 𝑞𝑅𝑡𝑜𝑡𝑎𝑙 + 𝑇∞ = 80 0.94 + 30 = 105℃
q 𝑅𝑡𝑜𝑡𝑎𝑙
𝑘 0.15
𝑟𝑐𝑟 = = = 12.5𝑚𝑚 ≫ 5.5𝑚𝑚 (2 × 𝑝𝑙𝑎𝑠𝑡𝑖𝑐 𝑐𝑜𝑣𝑒𝑟)
ℎ 12
𝑟𝑐𝑟 = 12.5 𝑚𝑚 represents condition for minimum resistance, which means increasing the plastic
cover thickness will decrease the thermal resistance until the 𝑟𝑐𝑟 goes beyond 12.5mm
Heat Transfer from a surface at a temperature Ts to the surrounding medium at T∞ is given by Newton’s Law of
cooling:
𝑞𝑐𝑜𝑛𝑣 = ℎ𝐴𝑠 𝑇𝑠 − 𝑇∞ (𝑊)
When temperatures Ts and T∞ are fixed by design considerations, there are two ways to increase the rate of heat
transfer:
• to increase the convection heat transfer coefficient h or
• to increase the surface area As
Increasing h may require installation of a pump or fan, or replacing the existing one with a larger one, which may
or may not be practical or adequate
Fin applications
• Cooling engine heads on motorcycles
• Cooling of electric power transformer
• To promote heat exchange between air
and working fluid in air conditioner
Fin configurations
• Straight fin
• Annular fin
• Pin fin
Steady State Heat Conduction
Heat Transfer from Extended Surfaces
Assumptions
𝑑 𝑑𝑇 𝑑 𝑑𝑇 ℎ 𝑑𝐴𝑠
0 = −𝑘 𝐴𝑐 𝑑𝑥 + ℎ𝑑𝐴𝑠 𝑇𝑠 − 𝑇∞ 𝐴𝑐 − 𝑇 − 𝑇∞ = 0
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑘 𝑑𝑥 𝑠
𝑑2𝑇 1 ℎ 𝑑𝐴𝑠
2
− 𝑇𝑠 − 𝑇∞ = 0
𝑑𝑥 𝐴𝑐 𝑘 𝑑𝑥
𝐴𝑠 = 𝑃x 𝑑𝐴𝑠
=𝑃
𝑑𝑥
𝑑2𝑇 ℎ𝑃
2
− 𝑇𝑠 − 𝑇∞ = 0
𝑑𝑥 𝑘𝐴𝑐
Steady State Heat Conduction
Heat Transfer from Extended Surfaces
Fins of Uniform Cross-Sectional Area
To simplify the form of this equation, we transform the dependent variable by defining an excess temperature as:
𝑑2𝑇 ℎ𝑃 𝑑2𝜃 2𝜃 = 0
ℎ𝑃
− 𝑇𝑠 − 𝑇∞ = 0 − 𝑚 𝑤ℎ𝑒𝑟𝑒 𝑚2 ≡
𝑑𝑥 2 𝑘𝐴𝑐 𝑑𝑥 2 𝑘𝐴𝑐
Eq. is a linear, homogeneous, second-order differential equation with constant coefficients
𝑑𝑇 𝑑𝜃
ℎ𝐴𝑐 𝑇 𝐿 − 𝑇∞ = −𝑘𝐴𝑐 ቤ ℎ𝜃(𝐿) = −𝑘 ቤ
𝑑𝑥 𝑥=𝐿 𝑑𝑥 𝑥=𝐿
Substituting the fin base
𝜃𝑏 = 𝐶1 + 𝐶2 ①
and fin tip boundary conditions
𝑑𝑇 𝑑𝜃
𝑞𝑓 = 𝑞𝑏 = −𝑘𝐴𝑐 ቤ = −𝑘𝐴𝑐 ቤ
𝑑𝑥 𝑥=0 𝑑𝑥 𝑥=0
Hence the expression for fin heat transfer rate from temperature
Distribution is as follows:
𝑞𝑓 = ℎ𝑃𝑘𝐴𝑐 𝜃𝑏 tanh 𝑚𝐿
𝜃 𝑥 = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
Steady State Heat Conduction
Heat Transfer from Extended Surfaces
The constant temperature boundary condition is concerned with prescribed temperature at the fin tip.
The boundary condition at fin tip 𝜃 𝐿 = 𝜃𝐿
Substituting the fin base Solving ①and ②for C1 and C2 and plugging in the θ(x)
equation (temperature distribution)
𝜃𝑏 = 𝐶1 + 𝐶2 ①
and fin tip boundary conditions 𝜃(𝑥) 𝜃𝐿 Τ𝜃𝑏 sinh 𝑚𝐿 + sinh 𝑚(𝐿 − 𝑥)
𝐶1 𝑒 𝑚𝐿 + 𝐶2 𝑒 −𝑚𝐿 = 𝜃𝐿 ② =
𝜃𝑏 si𝑛ℎ 𝑚𝐿
cosh 𝑚𝐿 − 𝜃𝐿 Τ𝜃𝑏
𝑞𝑓 = ℎ𝑃𝑘𝐴𝑐 𝜃𝑏
sinh 𝑚𝐿 𝜃 𝑥 = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
Steady State Heat Conduction
Heat Transfer from Extended Surfaces
The Infinite fin length boundary condition is concerned with very long fin (L= ∞).
The boundary condition at Fin tip: q ∞ = 0
Substituting the fin base Solving ①and ②for C1 and C2 and plugging in the θ(x)
equation (temperature distribution)
𝜃𝑏 = 𝐶1 + 𝐶2 ①
and fin tip boundary conditions 𝜃(𝑥)
𝐶2 𝑒 −𝑚∞ − 𝐶1 𝑒 𝑚∞ =0 ② = 𝑒 −𝑚𝑥
𝜃𝑏
𝑞𝑓 = ℎ𝑃𝑘𝐴𝑐 𝜃𝑏
𝜃 𝑥 = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
Steady State Heat Conduction
Heat Transfer from Extended Surfaces
Summarized Table
Steady State Heat Conduction
Fin Performance Parameters: Fin Effectiveness
Fin effectiveness (ɛf): The ratio of the fin heat transfer rate to the heat transfer rate that
would exist without the fin 𝑇
𝑞𝑓
𝑞𝑓 𝑞𝑓 𝑏
𝜀𝑓 = = 𝐴𝑏
ℎ𝐴𝑏 𝑇𝑏 − 𝑇∞ ℎ𝐴𝑏 𝜃𝑏
Where Ab is the cross-sectional area of the base.
𝑞𝑓 𝑞𝑓 𝑞𝑓
𝜼𝑓 = = =
𝑞𝑚𝑎𝑥 ℎ𝐴𝑓 𝑇𝑏 − 𝑇∞ ℎ𝐴𝑓 𝜃𝑏
ሶ 2
𝑞𝐿
𝑇 0 = 𝑇𝑜 = + 𝑇𝑏
2𝑘
𝑇𝑜
𝑇𝑏 𝑇∞ 𝑞ሶ 𝜋𝑟 2 𝐿
𝑞ሶ 𝜋𝑟 2 𝐿 = ℎ𝑃𝑘𝐴𝑐 𝑇𝑏 − 𝑇∞ 𝑇𝑏 = + 𝑇∞
𝑞ሶ Fin of infinite length ℎ𝑃𝑘𝐴𝑐
15mm
x
(7.5 × 106 )(𝜋 × 0.00252 × 0.015)
𝑞𝑓 = 𝑞𝑥 (𝐿) Case 4 for Fins (Infinite length) 𝑇𝑏 = + 20 = 272℃
𝑞𝑓 = ℎ𝑃𝑘𝐴𝑐 𝜃𝑏 = ℎ𝑃𝑘𝐴𝑐 𝑇𝑏 − 𝑇∞ (10 × 𝜋 × 0.005 × 25 × 𝜋 × 0.00252 )
ℎ𝑃 250 × 0.11
𝑚= = = 47.87 𝑚−1
𝑘𝐴𝑐 20 × 0.0006
𝜃𝑏
𝜃(𝐿) =
cosh 𝑚𝐿
𝑇𝑏 − 𝑇∞ 𝑇𝑏 − 𝑇∞ 300 − 1200
𝑇𝐿 − 𝑇∞ = 𝑇𝐿 = 𝑇∞ + = 1200 + = 1037℃ < maximum allowable blade
cosh 𝑚𝐿 cosh 𝑚𝐿 cosh(47.87 × 0.05) temperature of 1050°C
Turbine blades mounted to a rotating disc in a gas turbine engine are exposed to a gas stream that is at T=1200°C and maintains a
convection coefficient of h=250 W/m2 K over the blade. The blades, which are fabricated from Inconel, k=20 W/m K, have a
length of L=50 mm. The blade profile has a uniform cross-sectional area of Ac=6x10-4 m2 and a perimeter of P =110 mm. A
proposed blade-cooling scheme, which involves routing air through the supporting disc, is able to maintain the base of each blade
at a temperature of Tb=300°C.
(a) If the maximum allowable blade temperature is 1050°C and the blade tip may be assumed to be adiabatic, is the proposed
cooling scheme satisfactory?
(b) For the proposed cooling scheme, what is the rate at which heat is transferred from each blade to the coolant?
Case 2- Adiabatic boundary condition
ℎ𝑃 250 × 0.11
𝑚= = = 47.87 𝑚−1
𝑘𝐴𝑐 20 × 0.0006
𝑞𝑓 = ℎ𝑃𝑘𝐴𝑐 𝜃𝑏 tanh 𝑚𝐿 = (250 × 0.11 × 20 × 6 × 10−4 ) 300 − 1200 tanh 47.87 × 0.05 = − 508𝑊
𝑞𝑏 = 𝑞𝑓 = 508𝑊
PRACTICE PROBLEMS
(Incropera, 6th edition)
Derive correlation of Critical thickness of insulation for sphere
Problem: 3.56