Extended Surfaces (Fin) : Heat Transfer Lectures Chemical Engineering Department University of Technology, Iraq
Extended Surfaces (Fin) : Heat Transfer Lectures Chemical Engineering Department University of Technology, Iraq
h, T∞
𝒅𝒒𝒄𝒐𝒏𝒗. = 𝒉𝑨ሗሺ𝑻 − 𝑻∞ ሻ
To qx T qx+dx TL
Z
dx
L
x
Base
Figure 1: Schematic diagram illustrating one-dimensional conduction and convection through a rectangular fin
We approach the problem by making an energy balance at steady state operating conditions
on an element of the fin of thickness dx as shown in the figure. Thus
Assumptions:
• Steady-state conditions.
• One-dimensional conduction along the fin
• Constant thermal conductivity.
• Negligible radiation exchange with surroundings.
38 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝑞𝑥 = 𝑞𝑥+𝑑𝑥 + ℎ𝐴ሗሺ𝑇 − 𝑇∞ ሻ
𝜕
𝑞𝑥 = 𝑞𝑥 + ሺ𝑞 ሻ𝑑𝑥 + ℎ𝐴ሗሺ𝑇 − 𝑇∞ ሻ
𝜕𝑥 𝑥
𝜕
𝑞𝑥 − 𝑞𝑥 − ሺ𝑞 ሻ𝑑𝑥 = ℎ𝐴ሗሺ𝑇 − 𝑇∞ ሻ
𝜕𝑥 𝑥
𝜕
− ሺ𝑞 ሻ𝑑𝑥 = ℎ𝐴ሗሺ𝑇 − 𝑇∞ ሻ
𝜕𝑥 𝑥
𝜕𝑇
𝑆𝑖𝑛𝑐𝑒, 𝑞𝑥 = −𝑘𝐴 , 𝐹𝑜𝑢𝑟𝑖𝑒𝑟 ′ 𝑠 𝑙𝑎𝑤, 𝐵𝑦 𝑠𝑢𝑏𝑠𝑡𝑖𝑡𝑢𝑡𝑖𝑛𝑔 𝑡ℎ𝑒 𝑎𝑏𝑜𝑣𝑒 𝑒𝑞𝑢𝑎𝑡𝑖𝑜𝑛, 𝑤𝑒 𝑜𝑏𝑡𝑎𝑖𝑛
𝜕𝑥
𝜕 𝜕𝑇
− (−𝑘𝐴 ) 𝑑𝑥 = ℎ𝐴ሗሺ𝑇 − 𝑇∞ ሻ
𝜕𝑥 𝜕𝑥
𝜕2𝑇
+𝑘𝐴 𝑑𝑥 = ℎ𝐴ሗሺ𝑇 − 𝑇∞ ሻ (1)
𝜕𝑥 2
ሗ + 2𝑏ሻ𝑑𝑥
𝑆𝑖𝑛𝑐𝑒, 𝐴ሗ = 𝑃𝑑𝑥 = ሺ2𝑍
𝜕2𝑇
𝑘𝐴 𝑑𝑥 = ℎ 𝑃𝑑𝑥ሺ𝑇 − 𝑇∞ ሻ
𝜕𝑥 2
𝜕2𝑇
𝑘𝐴 2 = ℎ𝑃 ሺ𝑇 − 𝑇∞ ሻ
𝜕𝑥
𝜕2𝑇 ℎ𝑃
2
= ሺ𝑇 − 𝑇∞ ሻ (2)
𝜕𝑥 𝑘𝐴
𝐿𝑒𝑡 𝜃 = 𝑇 − 𝑇∞
𝜕2𝜃 ℎ𝑃
2
= 𝜃 (3)
𝜕𝑥 𝑘𝐴
ℎ𝑃
𝐿𝑒𝑡 𝑚2 =
𝑘𝐴
39 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
Where
h: heat transfer coefficient, W/m2. oC
P: perimeter of rectangular fin (2Z+2b), m
k: thermal conductivity of the fin, W/m. oC
A: cross-sectional area of the rectangular fin (bZ), m2
𝜕2𝜃
= 𝑚2 𝜃
𝜕𝑥 2
𝜕2𝜃
− 𝑚2 𝜃 = 0 (4)
𝜕𝑥 2
The above equation (i.e., Eq.4) is a linear, homogeneous, second order differential equation
with constant coefficients. Its general solution is of the form
𝐴+𝐵 𝐴−𝐵
𝑤ℎ𝑒𝑟𝑒 𝐶1 = 𝑎𝑛𝑑 𝐶2 =
2 2
𝑒 𝑚𝑥 − 𝑒 −𝑚𝑥
sinh 𝑚𝑥 =
2
𝑒 𝑚𝑥 + 𝑒 −𝑚𝑥
cosh 𝑚𝑥 =
2
Case 1: The fin is of finite length and loses heat by convection from its end. Thus
𝐵. 𝐶. 2, 𝑎𝑡 𝑥 = 𝐿, 𝑇 = 𝑇𝐿 ⇒ 𝜃 = 𝑇 − 𝑇∞ = 𝑇𝐿 − 𝑇∞ = 𝜃𝐿
40 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
Case 2: The fin is very long, and the temperature at the end of the fin is essentially that of the
surrounding fluid. Thus
𝐵. 𝐶. 3, 𝑎𝑡 𝑥 = ∞, 𝑇 = 𝑇∞ ⇒ 𝜃 = 𝑇 − 𝑇∞ = 𝑇∞ − 𝑇∞ = 0
𝜕𝑇 𝜕𝜃
𝐵. 𝐶. 4, 𝑎𝑡 𝑥 = 𝐿, = =0
𝜕𝑥 𝜕𝑥
𝜕𝜃 𝜕𝜃
𝐵. 𝐶. 5, 𝑎𝑡 𝑥 = 𝐿, −𝑘𝐴 | = ℎ𝐴𝜃𝐿 ⇒ −𝑘 | = ℎ𝜃𝐿
𝜕𝑥 𝑥=𝐿 𝜕𝑥 𝑥=𝐿
To find the temperature distribution for the rectangular fin, we need to find the C1 and
C2. Therefore, we will use the above boundary conditions to find that. So, let us start from the
general solution as follows:
𝐵. 𝐶. 1, 𝑎𝑡 𝑥 = 0, 𝑇 = 𝑇𝑜 ⇒ 𝜃 = 𝑇 − 𝑇∞ = 𝑇𝑜 − 𝑇∞ = 𝜃𝑜
𝜃𝑜 = 𝐶1 𝑒ሺ0ሻ + 𝐶2 𝑒−ሺ0ሻ = 𝐶1 + 𝐶2 ⇒ 𝐶1 = 𝜃𝑜 − 𝐶2
𝐶1 = 𝜃𝑜 − 𝐶2 (6)
By substituting B.C.3 into Eq. 5, we obtain
𝐵. 𝐶. 3, 𝑎𝑡 𝑥 = ∞, 𝑇 = 𝑇∞ ⇒ 𝜃 = 𝑇 − 𝑇∞ = 𝑇∞ − 𝑇∞ = 0
𝜃ሺ𝑥ሻ = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
𝐶2 0
0 = 𝐶1 𝑒 ∞ + 𝐶2 𝑒 −∞ ⇒ 0 = 𝐶1 𝑒 ∞ + +∞
⇒ 0 = 𝐶1 𝑒 ∞ + 0 ⇒ 𝐶1 = ∞ = 0
𝑒 𝑒
𝐶1 = 𝜃𝑜 − 𝐶2 ⇒ 0 = 𝜃𝑜 − 𝐶2 ⇒ 𝐶2 = 𝜃𝑜
Now by substituting the values of 𝐶1 and 𝐶2 into the general solution equation (Eq.5), we obtain
41 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝜃ሺ𝑥ሻ = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
𝜃ሺ𝑥ሻ = ሺ0ሻ𝑒 𝑚𝑥 + 𝜃𝑜 𝑒 −𝑚𝑥
Therefore, we obtain
𝜃 = 𝜃𝑜 𝑒 −𝑚𝑥
𝜃 𝑇 − 𝑇∞ ℎ𝑝
−√ 𝑥
= 𝑒 −𝑚𝑥 = = 𝑒 𝑘𝐴 (7)
𝜃𝑜 𝑇𝑜 − 𝑇∞
Note: Eq.7 represents the temperature distribution in the fin. This temperature distribution was
obtained by using B.C.1 and B.C.3. Therefore, the temperature distribution will be different if we
use different boundary conditions. We will see that in the next section.
Now let us see what the temperature distribution will be if we use the 1st and 2 nd boundary
conditions (Case 1).
So let us start with the general solution
𝜃ሺ𝑥ሻ = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥 (5)
By substituting B.C1 into Eq. 5, we obtain
𝐵. 𝐶. 1, 𝑎𝑡 𝑥 = 0, 𝑇 = 𝑇𝑜 ⇒ 𝜃 = 𝑇 − 𝑇∞ = 𝑇𝑜 − 𝑇∞ = 𝜃𝑜
𝜃𝑜 = 𝐶1 𝑒ሺ0ሻ + 𝐶2 𝑒−ሺ0ሻ = 𝐶1 + 𝐶2 ⇒ 𝐶1 = 𝜃𝑜 − 𝐶2
𝐶1 = 𝜃𝑜 − 𝐶2 (6)
By substituting B.C.2 into Eq. 5, we obtain
𝐵. 𝐶. 2, 𝑎𝑡 𝑥 = 𝐿, 𝑇 = 𝑇𝐿 ⇒ 𝜃 = 𝑇 − 𝑇∞ = 𝑇𝐿 − 𝑇∞ = 𝜃𝐿
𝜃ሺ𝑥ሻ = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
𝜃𝐿 = 𝐶1 𝑒𝑚𝐿 + 𝐶2 𝑒−𝑚𝐿
𝑠𝑖𝑛𝑐𝑒, 𝐶1 = 𝜃𝑜 − 𝐶2
𝜃𝐿 = ሺ𝜃𝑜 − 𝐶2 ሻ𝑒𝑚𝐿 + 𝐶2 𝑒−𝑚𝐿
𝐶2 𝑒 𝑚𝐿 − 𝐶2 𝑒 −𝑚𝐿 = 𝜃𝑜 𝑒 𝑚𝐿 − 𝜃𝐿
ሺ𝑒 𝑚𝐿 − 𝑒 −𝑚𝐿 ሻ𝐶2 = 𝜃𝑜 𝑒 𝑚𝐿 − 𝜃𝐿
42 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝜃𝑜 𝑒 𝑚𝐿 − 𝜃𝐿
𝐶2 =
𝑒 𝑚𝐿 − 𝑒 −𝑚𝐿
By substituting the C2 into Eq. 6, we obtain
𝜃𝑜 𝑒 𝑚𝐿 − 𝜃𝐿
𝐶1 = 𝜃𝑜 − 𝐶2 ⇒ 𝐶1 = 𝜃𝑜 −
𝑒 𝑚𝐿 − 𝑒 −𝑚𝐿
𝜃𝑜 𝑒 𝑚𝐿 − 𝜃𝐿
𝐶1 = 𝜃𝑜 −
𝑒 𝑚𝐿 − 𝑒 −𝑚𝐿
𝜃𝑜 𝑒 𝑚𝐿 − 𝜃𝑜 𝑒 −𝑚𝐿 − 𝜃𝑜 𝑒 𝑚𝐿 + 𝜃𝐿
𝐶1 =
𝑒 𝑚𝐿 − 𝑒 −𝑚𝐿
𝜃𝐿 − 𝜃𝑜 𝑒 −𝑚𝐿
𝐶1 =
𝑒 𝑚𝐿 − 𝑒 −𝑚𝐿
Now by substituting the values of 𝐶1 and 𝐶2 into the general solution equation (Eq.5), we obtain
𝜃ሺ𝑥ሻ = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
𝜃𝐿 − 𝜃𝑜 𝑒 −𝑚𝐿 𝑚𝑥
𝜃𝑜 𝑒 𝑚𝐿 − 𝜃𝐿
𝜃=( ) 𝑒 + ( ) 𝑒 −𝑚𝑥
𝑒 𝑚𝐿 − 𝑒 −𝑚𝐿 𝑒 𝑚𝐿 − 𝑒 −𝑚𝐿
43 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝐵. 𝐶. 1, 𝑎𝑡 𝑥 = 0, 𝑇 = 𝑇𝑜 ⇒ 𝜃 = 𝑇 − 𝑇∞ = 𝑇𝑜 − 𝑇∞ = 𝜃𝑜
𝜃𝑜 = 𝐶1 𝑒ሺ0ሻ + 𝐶2 𝑒−ሺ0ሻ = 𝐶1 + 𝐶2 ⇒ 𝐶1 = 𝜃𝑜 − 𝐶2
𝐶1 = 𝜃𝑜 − 𝐶2 (6)
By substituting B.C.4 into Eq. 5, we obtain
𝜕𝑇 𝜕𝜃
𝐵. 𝐶. 4, 𝑎𝑡 𝑥 = 𝐿, = =0
𝜕𝑥 𝜕𝑥
𝜃ሺ𝑥ሻ = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
We need to find the first derivative for the temperature distribution
𝜕𝜃
= 𝐶1 𝑒 𝑚𝑥 ሺ𝑚ሻ + 𝐶2 𝑒 −𝑚𝑥 ሺ−𝑚ሻ
𝜕𝑥
𝜕𝜃
= 𝑚𝐶1 𝑒 𝑚𝑥 − 𝑚𝐶2 𝑒 −𝑚𝑥 (13)
𝜕𝑥
Now we can substitute B.C.4 into Eq. 13
𝜕𝜃
| = 0 = 𝑚𝐶1 𝑒 𝑚𝐿 − 𝑚𝐶2 𝑒 −𝑚𝐿
𝜕𝑥 𝑥=𝐿
0 = 𝑚𝐶1 𝑒 𝑚𝐿 − 𝑚𝐶2 𝑒 −𝑚𝐿
(0 = 𝑚𝐶1 𝑒𝑚𝐿 − 𝑚𝐶2 𝑒−𝑚𝐿 ) ÷ 𝑚
0 = 𝐶1 𝑒 𝑚𝐿 − 𝐶2 𝑒 −𝑚𝐿 (14)
𝑆𝑖𝑛𝑐𝑒, 𝐶1 = 𝜃𝑜 − 𝐶2 , 𝑠𝑢𝑏. 𝑖𝑛𝑡𝑜 𝐸𝑞. 14
44 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
0 = ሺ𝜃𝑜 − 𝐶2 ሻ𝑒 𝑚𝐿 − 𝐶2 𝑒 −𝑚𝐿
0 = 𝜃𝑜 𝑒 𝑚𝐿 − 𝐶2 𝑒 𝑚𝐿 − 𝐶2 𝑒 −𝑚𝐿
𝐶2 𝑒 𝑚𝐿 + 𝐶2 𝑒 −𝑚𝐿 = 𝜃𝑜 𝑒 𝑚𝐿
𝜃𝑜 𝑒 𝑚𝐿
𝐶2 ሺ𝑒 𝑚𝐿 +𝑒 −𝑚𝐿 ሻ
= 𝜃𝑜 𝑒 𝑚𝐿
⇒ 𝐶2 = , 𝑠𝑢𝑏. 𝑖𝑛𝑡𝑜 𝐸𝑞. 6
𝑒 𝑚𝐿 + 𝑒 −𝑚𝐿
𝜃𝑜 𝑒 𝑚𝐿
𝐶1 = 𝜃𝑜 − 𝐶2 ⇒ 𝐶1 = 𝜃𝑜 −
𝑒 𝑚𝐿 + 𝑒 −𝑚𝐿
𝜃𝑜 𝑒 𝑚𝐿 + 𝜃𝑜 𝑒 −𝑚𝐿 − 𝜃𝑜 𝑒 𝑚𝐿
𝐶1 =
𝑒 𝑚𝐿 + 𝑒 −𝑚𝐿
𝜃𝑜 𝑒 −𝑚𝐿
𝐶1 =
𝑒 𝑚𝐿 + 𝑒 −𝑚𝐿
Now by substituting the values of 𝐶1 and 𝐶2 in the general solution (Eq.5), we obtain
𝜃 = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
𝜃𝑜 𝑒 −𝑚𝐿 𝑚𝑥
𝜃𝑜 𝑒 𝑚𝐿
𝜃=( )𝑒 + ( 𝑚𝐿 ) 𝑒 −𝑚𝑥
𝑒 𝑚𝐿 + 𝑒 −𝑚𝐿 𝑒 +𝑒 −𝑚𝐿
𝜃𝑜 𝑒 −𝑚𝐿 𝑒 𝑚𝑥 𝜃𝑜 𝑒 𝑚𝐿 𝑒 −𝑚𝑥
𝜃 = ( 𝑚𝐿 ) + ( 𝑚𝐿 )
𝑒 + 𝑒 −𝑚𝐿 𝑒 + 𝑒 −𝑚𝐿
𝜃𝑜 𝑒 −𝑚ሺ𝐿−𝑥ሻ 𝜃𝑜 𝑒 𝑚ሺ𝐿−𝑥ሻ
𝜃 = ( 𝑚𝐿 ) + ( )
𝑒 + 𝑒 −𝑚𝐿 𝑒 𝑚𝐿 + 𝑒 −𝑚𝐿
𝜃𝑜 𝑒 −𝑚ሺ𝐿−𝑥ሻ + 𝜃𝑜 𝑒 𝑚ሺ𝐿−𝑥ሻ
𝜃=
𝑒 𝑚𝐿 + 𝑒 −𝑚𝐿
𝜃𝑜 (𝑒 𝑚ሺ𝐿−𝑥ሻ + 𝑒 −𝑚ሺ𝐿−𝑥ሻ )
𝜃= (15)
𝑒 𝑚𝐿 + 𝑒 −𝑚𝐿
Since
𝑒 𝑚𝐿 + 𝑒 −𝑚𝐿
cosh 𝑚𝐿 =
2
and
𝑒 𝑚ሺ𝐿−𝑥ሻ + 𝑒 −𝑚ሺ𝐿−𝑥ሻ
cosh 𝑚ሺ𝐿 − 𝑥ሻ =
2
45 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝜃𝑜 2 cosh 𝑚ሺ𝐿 − 𝑥ሻ
𝜃=
2 cosh 𝑚𝐿
𝜃 cosh 𝑚ሺ𝐿 − 𝑥ሻ
= (16)
𝜃𝑜 cosh 𝑚𝐿
Note: Eq.16 represents the temperature distribution in the fin. This temperature distribution was
obtained by using B.C.1 and B.C.4.
Now let us see what will be the temperature distribution of fin if we use the 1st and 5 th (Case 4)
boundary conditions.
So let us start with the general solution
𝜃ሺ𝑥ሻ = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥 (5)
By substituting B.C1 into Eq. 5, we obtain
𝐵. 𝐶. 1, 𝑎𝑡 𝑥 = 0, 𝑇 = 𝑇𝑜 ⇒ 𝜃 = 𝑇 − 𝑇∞ = 𝑇𝑜 − 𝑇∞ = 𝜃𝑜
𝜃𝑜 = 𝐶1 𝑒ሺ0ሻ + 𝐶2 𝑒−ሺ0ሻ = 𝐶1 + 𝐶2 ⇒ 𝐶1 = 𝜃𝑜 − 𝐶2
𝐶1 = 𝜃𝑜 − 𝐶2 (6)
By substituting B.C.2 into Eq. 5, we obtain
𝐵. 𝐶. 2, 𝑎𝑡 𝑥 = 𝐿, 𝑇 = 𝑇𝐿 ⇒ 𝜃 = 𝑇 − 𝑇∞ = 𝑇𝐿 − 𝑇∞ = 𝜃𝐿
𝜃 = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
𝜃𝐿 = 𝐶1 𝑒𝑚𝐿 + 𝐶2 𝑒−𝑚𝐿 (17)
By substituting B.C.5 into Eq. 5, we obtain
𝜕𝜃 𝜕𝜃
𝐵. 𝐶. 5, 𝑎𝑡 𝑥 = 𝐿, −𝑘𝐴 | = ℎ𝐴𝜃𝐿 ⇒ −𝑘 | = ℎ𝜃𝐿
𝜕𝑥 𝑥=𝐿 𝜕𝑥 𝑥=𝐿
−𝑘 𝜕𝜃
𝜃𝐿 = | (18)
ℎ 𝜕𝑥 𝑥=𝐿
𝑆𝑖𝑛𝑐𝑒, 𝜃𝐿 = 𝐶1 𝑒 𝑚𝐿 + 𝐶2 𝑒 −𝑚𝐿
Hence, Eq. 17 becomes as follows:
−𝑘 𝜕𝜃
| = 𝐶1 𝑒 𝑚𝐿 + 𝐶2 𝑒 −𝑚𝐿 (19)
ℎ 𝜕𝑥 𝑥=𝐿
Now we need to find the first derivative for the general solution at x=L, and then we can submit
it into Eq. 19
46 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝜃ሺ𝑥ሻ = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
𝜕𝜃
= 𝐶1 𝑒 𝑚𝑥 ሺ𝑚ሻ + 𝐶2 𝑒 −𝑚𝑥 ሺ−𝑚ሻ
𝜕𝑥
𝜕𝜃
= 𝑚𝐶1 𝑒 𝑚𝑥 − 𝑚𝐶2 𝑒 −𝑚𝑥
𝜕𝑥
𝜕𝜃
| = 𝑚𝐶1 𝑒 𝑚𝐿 − 𝑚𝐶2 𝑒 −𝑚𝐿
𝜕𝑥 𝑥=𝐿
By submitting the above equation into Eq. 19, we obtain
−𝑘 𝜕𝜃
| = 𝐶1 𝑒 𝑚𝐿 + 𝐶2 𝑒 −𝑚𝐿 (19)
ℎ 𝜕𝑥 𝑥=𝐿
−𝑘
[𝑚𝐶1 𝑒 𝑚𝐿 − 𝑚𝐶2 𝑒 −𝑚𝐿 ] = 𝐶1 𝑒 𝑚𝐿 + 𝐶2 𝑒 −𝑚𝐿
ℎ
−𝑚𝑘 𝑚𝑘
[ 𝐶1 𝑒 𝑚𝐿 + 𝐶 𝑒 −𝑚𝐿 ] = 𝐶1 𝑒 𝑚𝐿 + 𝐶2 𝑒 −𝑚𝐿
ℎ ℎ 2
𝑚𝑘 𝑚𝑘
𝐶2 𝑒 −𝑚𝐿 − 𝐶2 𝑒 −𝑚𝐿 = 𝐶1 𝑒 𝑚𝐿 + 𝐶 𝑒 𝑚𝐿
ℎ ℎ 1
𝑚𝑘 𝑚𝑘
𝐶2 𝑒 −𝑚𝐿 ( − 1) = 𝐶1 𝑒 𝑚𝐿 (1 + )
ℎ ℎ
𝑚𝑘 − ℎ ℎ + 𝑚𝑘
𝐶2 𝑒 −𝑚𝐿 ( ) = 𝐶1 𝑒 𝑚𝐿 ( )
ℎ ℎ
𝐶2 𝑒 −𝑚𝐿 ሺ𝑚𝑘 − ℎሻ = 𝐶1 𝑒 𝑚𝐿 ሺ𝑚𝑘 + ℎሻ
𝑚𝑘 − ℎ
( ) 𝐶 𝑒 −𝑚𝐿 = 𝐶1 𝑒 𝑚𝐿 (20)
𝑚𝑘 + ℎ 2
𝑚𝑘 − ℎ
𝑙𝑒𝑡 𝑎 = ( ) , 𝑠𝑢𝑏. 𝑖𝑛𝑡𝑜 𝐸𝑞. 20, 𝑤𝑒 𝑜𝑏𝑡𝑎𝑖𝑛
𝑚𝑘 + ℎ
𝑎𝐶2 𝑒 −𝑚𝐿 = 𝐶1 𝑒 𝑚𝐿 (21)
𝑆𝑖𝑛𝑐𝑒, 𝐶1 = 𝜃𝑜 − 𝐶2 , ℎ𝑒𝑛𝑐𝑒 𝐸𝑞. 21 𝑏𝑒𝑐𝑜𝑚𝑒𝑠
𝑎𝐶2 𝑒 −𝑚𝐿 = ሺ𝜃𝑜 − 𝐶2 ሻ𝑒 𝑚𝐿
𝑎𝐶2 𝑒 −𝑚𝐿 = 𝜃𝑜 𝑒 𝑚𝐿 − 𝐶2 𝑒 𝑚𝐿
𝐶2 𝑒 𝑚𝐿 + 𝑎𝐶2 𝑒 −𝑚𝐿 = 𝜃𝑜 𝑒 𝑚𝐿
𝜃𝑜 𝑒 𝑚𝐿
𝐶2 ሺ𝑒 𝑚𝐿 + 𝑎𝑒 −𝑚𝐿 ሻ = 𝜃𝑜 𝑒 𝑚𝐿 ⇒ 𝐶2 =
𝑒 𝑚𝐿 + 𝑎𝑒 −𝑚𝐿
47 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝜃𝑜 𝑒 𝑚𝐿
𝑆𝑖𝑛𝑐𝑒, 𝐶1 = 𝜃𝑜 − 𝐶2 ⇒ 𝐶1 = 𝜃𝑜 −
𝑒 𝑚𝐿 + 𝑎𝑒 −𝑚𝐿
𝜃𝑜 𝑒 𝑚𝐿 + 𝑎𝜃𝑜 𝑒 −𝑚𝐿 − 𝜃𝑜 𝑒 𝑚𝐿
𝐶1 =
𝑒 𝑚𝐿 + 𝑎𝑒 −𝑚𝐿
𝑎𝜃𝑜 𝑒 −𝑚𝐿
𝐶1 =
𝑒 𝑚𝐿 + 𝑎𝑒 −𝑚𝐿
Now by substituting the values of 𝐶1 and 𝐶2 in the general solution equation (Eq. 5), we obtain
𝜃 = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥 (5)
𝑎𝜃𝑜 𝑒 −𝑚𝐿 𝑚𝑥
𝜃𝑜 𝑒 𝑚𝐿
𝜃 = ( 𝑚𝐿 ) 𝑒 + ( ) 𝑒 −𝑚𝑥
𝑒 + 𝑎𝑒 −𝑚𝐿 𝑒 𝑚𝐿 + 𝑎𝑒 −𝑚𝐿
48 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝑒 𝑚ሺ𝐿−𝑥ሻ + 𝑒 −𝑚ሺ𝐿−𝑥ሻ
cosh 𝑚ሺ𝐿 − 𝑥ሻ =
2
𝑒 𝑚ሺ𝐿−𝑥ሻ − 𝑒 −𝑚ሺ𝐿−𝑥ሻ
sinh 𝑚ሺ𝐿 − 𝑥ሻ =
2
𝑒 𝑚𝐿 + 𝑒 −𝑚𝐿
cosh 𝑚𝐿 =
2
𝑒 𝑚𝐿 − 𝑒 −𝑚𝐿
sinh 𝑚𝐿 =
2
Hence, Eq. 23 becomes
𝜃 𝑚𝑘ሺ2 cosh 𝑚ሺ𝐿 − 𝑥ሻሻ + ℎሺ2 sinh 𝑚ሺ𝐿 − 𝑥ሻሻ
=
𝜃𝑜 𝑚𝑘ሺ2 cosh 𝑚𝐿ሻ + ℎሺ2 sinh 𝑚𝐿ሻ
ℎ
𝜃 2𝑚𝑘 [cosh 𝑚ሺ𝐿 − 𝑥ሻ + sinh 𝑚ሺ𝐿 − 𝑥ሻ]
= 𝑚𝑘
𝜃𝑜 ℎ
2𝑚𝑘 [cosh 𝑚𝐿 + sinh 𝑚𝐿]
𝑚𝑘
ℎ
𝜃 cosh 𝑚ሺ𝐿 − 𝑥ሻ + sinh 𝑚ሺ𝐿 − 𝑥ሻ
= 𝑚𝑘 (24)
𝜃𝑜 ℎ
cosh 𝑚𝐿 + sinh 𝑚𝐿
𝑚𝑘
Note: Eq.24 represents the temperature distribution in the fin. This temperature distribution was
obtained by using B.C.1 and B.C.5 (Case 4).
49 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝜃 cosh 𝑚ሺ𝐿 − 𝑥ሻ
= (2)
𝜃𝑜 cosh 𝑚𝐿
Therefore, we need to find the first derivative of the given temperature distribution at x=0
𝜃𝑜
𝜃= cosh 𝑚ሺ𝐿 − 𝑥ሻ
cosh 𝑚𝐿
𝜕𝜃 𝜃𝑜
= sinh 𝑚ሺ𝐿 − 𝑥ሻሺ−𝑚ሻ
𝜕𝑥 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚ሺ𝐿 − 𝑥ሻ
= −𝑚𝜃𝑜
𝜕𝑥 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚ሺ𝐿 − 0ሻ
| = −𝑚𝜃𝑜
𝜕𝑥 𝑥=0 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚𝐿
| = −𝑚𝜃𝑜
𝜕𝑥 𝑥=0 cosh 𝑚𝐿
𝜕𝜃
| = −𝑚𝜃𝑜 tanh 𝑚𝐿 (3)
𝜕𝑥 𝑥=0
By substituting Eq. 3 into Eq. 1, we obtain
𝜕𝜃
𝑞𝑓 = −𝑘𝐴 | (1)
𝜕𝑥 𝑥=0
𝑞𝑓 = −𝑘𝐴ሺ−𝑚𝜃𝑜 tanh 𝑚𝐿ሻ
𝑞𝑓 = 𝑚𝑘𝐴𝜃𝑜 tanh 𝑚𝐿
ℎ𝑝
𝑤ℎ𝑒𝑟𝑒, 𝑚 = √
𝑘𝐴
50 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝜕𝜃
−𝑘𝐴 |
𝑞𝑓 ℎ𝐴ሗሺ𝑇 − 𝑇∞ ሻ 𝜕𝑥 𝑥=0
𝜂𝑓 = = =
𝑞𝑚𝑎𝑥 ℎ𝐴ሗሺ𝑇𝑜 − 𝑇∞ ሻ ℎ𝐴ሗሺ𝑇𝑜 − 𝑇∞ ሻ
𝜕𝜃 𝜕𝜃
−𝑘𝐴 | −𝑘𝐴
|
𝜕𝑥 𝑥=0 𝜕𝑥 𝑥=0
𝜂𝑓 = =
ℎ𝑨ሗሺ𝑇𝑜 − 𝑇∞ ሻ ℎሺ𝑃𝐿ሻሺ𝑇𝑜 − 𝑇∞ ሻ
𝜕𝜃
−𝑘𝐴 |
𝜕𝑥 𝑥=0 (1)
𝜂𝑓 =
ℎሺ𝑃𝐿ሻሺ𝑇𝑜 − 𝑇∞ ሻ
For example, calculate the fin efficiency when the temperature distribution as follows:
𝜃 cosh 𝑚ሺ𝐿 − 𝑥ሻ
= (2)
𝜃𝑜 cosh 𝑚𝐿
Since
𝜕𝜃
−𝑘𝐴 |
𝜕𝑥 𝑥=0 (1)
𝜂𝑓 =
ℎሺ𝑃𝐿ሻሺ𝑇𝑜 − 𝑇∞ ሻ
Therefore, we need to find the first derivative of the given temperature distribution (Eq. 2) at x=0
𝜃𝑜
𝜃= cosh 𝑚ሺ𝐿 − 𝑥ሻ
cosh 𝑚𝐿
𝜕𝜃 𝜃𝑜
= sinh 𝑚ሺ𝐿 − 𝑥ሻሺ−𝑚ሻ
𝜕𝑥 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚ሺ𝐿 − 𝑥ሻ
= −𝑚𝜃𝑜
𝜕𝑥 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚ሺ𝐿 − 0ሻ
| = −𝑚𝜃𝑜
𝜕𝑥 𝑥=0 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚𝐿
| = −𝑚𝜃𝑜
𝜕𝑥 𝑥=0 cosh 𝑚𝐿
𝜕𝜃
| = −𝑚𝜃𝑜 tanh 𝑚𝐿 , sub. into Eq. 1
𝜕𝑥 𝑥=0
𝜕𝜃
−𝑘𝐴 |
𝜕𝑥 𝑥=0 +𝑘𝐴 𝑚𝜃𝑜 tanh 𝑚𝐿
𝜂𝑓 = =
ℎ𝑃𝐿𝜃𝑜 ℎ𝑃𝐿𝜃𝑜
51 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝑘𝐴 𝑚 tanh 𝑚𝐿 𝑚 tanh 𝑚𝐿
𝜂𝑓 = =
ℎ𝑃𝐿 ℎ𝑃 (3)
𝐿
𝑘𝐴
Since
ℎ𝑃 ℎ𝑃
𝑚=√ ⇒ 𝑚2 =
𝐾𝐴 𝐾𝐴
52 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
h=5 W/m2. K,
T∞=290 K
k=50 W/m. K
b=3 mm
To=400 K
L=75 mm
x
Base
𝜕𝑇 𝜕𝜃
𝐵. 𝐶. 4, 𝑎𝑡 𝑥 = 𝐿, | = | =0
𝜕𝑥 𝑥=𝐿 𝜕𝑥 𝑥=𝐿
By substituting the above boundary conditions in the general solution, we obtain
𝜃 cosh 𝑚ሺ𝐿 − 𝑥ሻ
=
𝜃𝑜 cosh 𝑚𝐿
𝜕𝜃
𝑞𝑓 = −𝑘𝐴 | = ℎ𝑨ሗ ሺ𝑇 − 𝑇∞ ሻ
𝜕𝑥 𝑥=0
𝜕𝜃
𝑞𝑓 = −𝑘𝐴 | (1)
𝜕𝑥 𝑥=0
Now we need to find the first derivative for the given temperature distribution at x=0, and then
we will submit in Eq. 1
𝜃 cosh 𝑚ሺ𝐿 − 𝑥ሻ
=
𝜃𝑜 cosh 𝑚𝐿
𝜃𝑜
𝜃= cosh 𝑚ሺ𝐿 − 𝑥ሻ
cosh 𝑚𝐿
53 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝜕𝜃 𝜃𝑜
= sinh 𝑚ሺ𝐿 − 𝑥ሻሺ−𝑚ሻ
𝜕𝑥 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚ሺ𝐿 − 𝑥ሻ
= −𝑚𝜃𝑜
𝜕𝑥 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚ሺ𝐿 − 0ሻ
| = −𝑚𝜃𝑜
𝜕𝑥 𝑥=0 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚𝐿
| = −𝑚𝜃𝑜
𝜕𝑥 𝑥=0 cosh 𝑚𝐿
𝜕𝜃
| = −𝑚𝜃𝑜 tanh 𝑚𝐿 , sub. into Eq. 1
𝜕𝑥 𝑥=0
𝜕𝜃
𝑞𝑓 = −𝑘𝐴 |
𝜕𝑥 𝑥=0
𝑞𝑓 = −𝑘𝐴ሺ−𝑚𝜃𝑜 tanh 𝑚𝐿ሻ
𝑞𝑓 = 𝑚𝑘𝐴𝜃𝑜 tanh 𝑚𝐿
ℎ𝑃
𝑚=√
𝑘𝐴
𝐴 = 𝑏𝑍 = 0.003ሺ1ሻ = 0.003 𝑚2
5 × 2.006
𝑚=√ = 8.18 𝑚−1
50 × 0.003
𝑞𝑓 = 𝑚𝑘𝐴𝜃𝑜 𝐿
54 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
Example 2: A very long copper fin with 2.5 cm in diameter and thermal conductivity of 372
W/m.oC. This fin has one end maintained at 93 oC. The fin exposed t a fluid whose temperature is
38 oC. If the heat transfer coefficient is 3.5 W/m2. oC, how much will be the heat lost by the fin?
Solution:
T∞=38 oC
h=3.5 W/m2. oC
To=93 oC
k= 372 W/m. oC
Assumptions:
• Steady-state conditions.
• One-dimensional conduction along the fin
• Constant thermal conductivity.
• Negligible radiation exchange with surroundings.
• The fin is very long
Since the general solution
𝜃 = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
𝐵. 𝐶. 1, 𝑎𝑡 𝑥 = 0, 𝑇 = 𝑇𝑜 ⇒ 𝜃 = 𝑇 − 𝑇∞ = 𝑇𝑜 − 𝑇∞ = 𝜃𝑜
The fin is very long, and the temperature at the end of the fin is essentially that of the
surrounding fluid. Thus
𝐵. 𝐶. 3, 𝑎𝑡 𝑥 = ∞, 𝑇 = 𝑇∞ ⇒ 𝜃 = 𝑇 − 𝑇∞ = 𝑇∞ − 𝑇∞ = 0
55 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝜃 𝑇 − 𝑇∞ ℎ𝑝
−√ 𝑥
= 𝑒 −𝑚𝑥 = = 𝑒 𝑘𝐴
𝜃𝑜 𝑇𝑜 − 𝑇∞
𝜕𝜃
𝑞𝑓 = −𝑘𝐴 | = ℎ𝑨ሗሺ𝑇 − 𝑇∞ ሻ
𝜕𝑥 𝑥=0
𝜕𝜃
𝑞𝑓 = −𝑘𝐴 |
𝜕𝑥 𝑥=0
Since
𝜃 = 𝜃𝑜 𝑒 −𝑚𝑥
𝜕𝜃
= 𝜃𝑜 𝑒 −𝑚𝑥 ሺ−𝑚ሻ = −𝑚𝜃𝑜 𝑒 −𝑚𝑥
𝜕𝑥
𝜕𝜃
| = −𝑚𝜃𝑜 𝑒 ሺ0ሻ = −𝑚𝜃𝑜
𝜕𝑥 𝑥=0
Therefore
𝜕𝜃
𝑞𝑓 = −𝑘𝐴 |
𝜕𝑥 𝑥=0
𝑞𝑓 = −𝑘𝐴ሺ−𝑚𝜃𝑜 ሻ = +𝑚𝑘𝐴𝜃𝑜
ℎ𝑃
𝑚=√
𝑘𝐴
ℎ𝑃 3.5 × 0.0785
𝑚=√ =√ = 1.22 𝑚−1
𝑘𝐴 372 × 0.0005
56 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
Example 3: A straight rectangular fin with 2.5 cm thick and 15 cm long is constructed of steel
(k=20 W/m. oC) and placed on the outside of a wall maintained at a temperature of 200 oC. The
environmental temperature is 15 oC and heat transfer coefficient of 17 W/m2. oC. If the loose end
of the fin is insulated, calculate the amount of heat loss per unit depth and fin efficiency. Given
that the heat transfer by radiation is negligible.
Solution:
Assumptions:
• Steady-state conditions.
• One-dimensional conduction along the fin
• Constant thermal conductivity.
• Negligible radiation exchange with surroundings.
• The end of the fin is insulated
b=2.5 cm
To=200 oC
L=15 cm
x
Base
𝜕𝑇 𝜕𝜃
𝐵. 𝐶. 4, 𝑎𝑡 𝑥 = 𝐿, | = | =0
𝜕𝑥 𝑥=𝐿 𝜕𝑥 𝑥=𝐿
By substituting the above boundary conditions in the general solution, we obtain
𝜃 cosh 𝑚ሺ𝐿 − 𝑥ሻ
= (2)
𝜃𝑜 cosh 𝑚𝐿
57 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝜕𝜃
Since, 𝑞𝑓 = −𝑘𝐴 | = ℎ𝑨ሗ ሺ𝑇 − 𝑇∞ ሻ
𝜕𝑥 𝑥=0
𝜕𝜃
𝑞𝑓 = −𝑘𝐴 | (3)
𝜕𝑥 𝑥=0
Hence, we need to find the first derivative for the given temperature distribution (Eq.2) at x=0,
and then we will submit into Eq. 3 to calculate the heat transfer rate of the fin.
𝜃 cosh 𝑚ሺ𝐿 − 𝑥ሻ
=
𝜃𝑜 cosh 𝑚𝐿
𝜃𝑜
𝜃= cosh 𝑚ሺ𝐿 − 𝑥ሻ
cosh 𝑚𝐿
𝜕𝜃 𝜃𝑜
= sinh 𝑚ሺ𝐿 − 𝑥ሻሺ−𝑚ሻ
𝜕𝑥 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚ሺ𝐿 − 𝑥ሻ
= −𝑚𝜃𝑜
𝜕𝑥 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚ሺ𝐿 − 0ሻ
| = −𝑚𝜃𝑜
𝜕𝑥 𝑥=0 cosh 𝑚𝐿
𝜕𝜃 sinh 𝑚𝐿
| = −𝑚𝜃𝑜
𝜕𝑥 𝑥=0 cosh 𝑚𝐿
𝜕𝜃
| = −𝑚𝜃𝑜 tanh 𝑚𝐿 (4)
𝜕𝑥 𝑥=0
By substituting Eq. 4 into Eq. 3 we obtain
𝜕𝜃
𝑞𝑓 = −𝑘𝐴 | (3)
𝜕𝑥 𝑥=0
𝑞𝑓 = −𝑘𝐴ሺ−𝑚𝜃𝑜 tanh 𝑚𝐿ሻ
𝑞𝑓 = 𝑚𝑘𝐴𝜃𝑜 tanh 𝑚𝐿
ℎ𝑃
𝑚=√
𝑘𝐴
𝜃𝑜 = 𝑇𝑜 − 𝑇∞ = 200 − 15 = 185 𝐶 °
58 2018-2019
Heat Transfer Lectures
Chemical Engineering Department
University of Technology, Iraq
𝐴 = 𝑏𝑍 = 0.025ሺ1ሻ = 0.025 𝑚2
17 × 2.05
𝑤ℎ𝑒𝑟𝑒, 𝑚 = √ = 8.35 𝑚−1
20 × 0.025
𝑞𝑓 = 𝑚𝑘𝐴𝜃𝑜 tanh 𝑚𝐿
59 2018-2019