HMT Combined Sir Asad
HMT Combined Sir Asad
ME-315
Recommended Book
Heat and Mass Transfer
by Younus Cengel and Afshin J. Ghajar, 5th Edition, McGraw Hill
Face EFGH is at
Face ABCD is at temperature 𝑇2
temperature 𝑇1 C F
H
A
𝑻𝟏 > 𝑻𝟐
D E
• Heat will flow from face ABCD to EFGH since ABCD is at higher temperature than EFGH.
• The area perpendicular to the direction of heat is ABCD or EFGH.
• The heat flux will be amount of heat transfered divided by the area ABCD.
Asad Akhter Naqvi
The wall of an industrial furnace is constructed from 0.15𝑚 thick fireclay brick having a thermal
conductivity of 1.7 𝑊/𝑚 𝐾 . Measurements made during steady-state operation reveal
temperatures of 1400 and 1150 𝐾 at the inner and outer surfaces, respectively. What is the rate
of heat loss through a wall that is 0.5 𝑚 × 1.2 𝑚 on a side?
∆𝑥 = 0.15𝑚 ∆𝑻
𝑘 = 1.7𝑊/𝑚 𝐾 𝒒 = −𝒌𝑨 ∆𝒙
𝑇2 = 1400𝐾
𝑇1 = 1150𝐾
𝐴 = 0.5𝑚 × 1.2𝑚
𝟏𝟒𝟎𝟎 − 𝟏𝟏𝟓𝟎
𝒒 = −𝟏. 𝟕 × 𝟎. 𝟓 × 𝟏. 𝟐 ×
𝟎. 𝟏𝟓
𝒒 = −𝟏𝟕𝟎𝟎 𝑾 Asad Akhter Naqvi
Convection
• Convection heat transfer, which occurs between a fluid in motion and a bounding surface when the two are at
different temperatures.
• Regardless of the nature of the convection heat transfer process, the appropriate rate equation is of the form:
𝒒𝒙 ′′ = 𝒉 𝑻𝒔 − 𝑻∞
• where , the convective heat flux (W/m2 ), is proportional to the difference between the surface and fluid
temperatures, 𝑇𝑠 and 𝑇∞ , respectively.
• This expression is known as Newton's law of cooling, and the parameter h (𝑊/𝑚2 𝐾) is termed the convection
heat transfer coefficient.
Asad Akhter Naqvi
Radiation
• Thermal radiation is energy emitted by matter that is at a nonzero temperature.
• The energy of the radiation field is transported by electromagnetic waves (or
alternatively, photons).
• While the transfer of energy by conduction or convection requires the presence of
a material medium, radiation does not.
• In fact, radiation transfer occurs most efficiently in a vacuum.
• the rate at which energy is released per unit area (𝑊/𝑚2 ) is termed the surface
emissive power, E.
• There is an upper limit to the emissive power, which is prescribed by the Stefan
Boltzmann law.
𝑬𝒃 = 𝝈𝑻𝒔 𝟒
• where 𝑇𝑠 is the absolute temperature (K) of the surface and 𝜎 is the Stefan
Boltzmann constant . Such a surface is called an ideal radiator or blackbody.
• The heat flux emitted by a real surface is less than that of a blackbody at the same
temperature and is given by
𝑬𝒃 = 𝝐𝝈𝑻𝒔 𝟒
Asad Akhter Naqvi
• The amount of heat transfer by radiation is given by
𝑻𝒔𝒖𝒓𝒓
𝑻𝒔
5
2
Some important points
➢ The cubical freezer compartment is shown in the fig.
➢ The bottom side which is denoted by 6 is perfectly 3
insulated. 1 4
➢ It means that remaining 5 sides (1 to 5) are
responsible of heat loss from the freezer
compartment.
6
Asad Akhter Naqvi
The heat transfer through the freezer compartment is given by Fourier’s law of
heat transfer.
𝑻𝟏 − 𝑻𝟐
𝒒 = 𝒌𝑨
∆𝒙 𝟑𝟓 + 𝟏𝟎
𝑞 = 500 𝑊 𝟓𝟎𝟎 = 𝟎. 𝟎𝟑𝟎 × 𝟐𝟎
𝒌 = 0.030 𝑊/𝑚 𝐾 ∆𝒙
𝑻𝟏 = 𝑶𝒖𝒕𝒔𝒊𝒅𝒆 𝒕𝒆𝒎𝒑𝒆𝒓𝒂𝒕𝒖𝒓𝒆 = 35𝑜 𝐶
𝑻𝟐 = 𝑰𝒏𝒔𝒊𝒅𝒆 𝒕𝒆𝒎𝒑𝒆𝒓𝒂𝒕𝒖𝒓𝒆 = −10𝑜
𝑨 = 𝑨𝒓𝒆𝒂 𝒐𝒇 𝒕𝒉𝒆 𝒄𝒐𝒎𝒑𝒂𝒓𝒕𝒎𝒆𝒏𝒕
Here Total area will be 5 × 2 × 2 = 20 𝑚2
∆𝒙 = 𝟎. 𝟎𝟓𝟒 𝒎
∆𝒙 = 𝟓𝟒 𝒎𝒎
Asad Akhter Naqvi
1-15. The 5 𝑚𝑚 thick bottom of a 200 𝑚𝑚 diameter pan may be made from
aluminum (𝑘 = 240 𝑊/𝑚 𝐾) or copper (𝑘 = 390 𝑊/𝑚 𝐾). When used to boil
water, the surface of the bottom exposed to the water is nominally at 110𝑜 𝐶. If
heat is transferred from the stove to the pan at a rate of 600 𝑊, what is the
temperature of the surface in contact with the stove for each of the two materials?
𝑇2 = 110𝑜 𝐶 𝐷 = 200 𝑚𝑚
∆𝑥 = 5 𝑚𝑚
𝑇1 =? 𝑞 = 600 𝑊
𝑇𝑠 = 90𝑜 𝐶
𝒒
𝑻∞ = 𝟐𝟓𝒐 𝑪
𝑽 = 6.25 × 10−5 ℎ2
𝑳 = 𝟐𝟎 𝒎𝒎 𝑻∞ = 𝟐𝟓𝒐 𝑪
𝑇 = 75𝑜 𝐶
𝑆
Asad Akhter Naqvi
The anemometer is heated by resistance heating. It means that heat dissipated is equal to the
input electrical power
𝒒 = 𝑷 = 𝑰𝑽
𝑰 = 𝟎. 𝟏 𝑨 𝒉 = 𝟑𝟏𝟖 𝑾/𝒎𝟐 𝑲
𝑽=𝟓𝑽
𝑨 = 𝝅𝑫𝑳 = 𝝅 × 𝟎. 𝟓 × 𝟏𝟎−𝟑 × 𝟎. 𝟎𝟐 = 𝟎. 𝟎𝟎𝟎𝟎𝟑𝟏𝟒𝟐 𝒎𝟐
𝑞𝑟𝑎𝑑
𝑞𝑐𝑜𝑛𝑣
𝑇∞ = 25𝑜 𝐶
𝑇𝑆 = 150𝑜 𝐶
• However, we would also find that, for equal values of 𝐴 , ∆𝑥 and ∆𝑇, the value of 𝑞𝑥 would be smaller for the
plastic than for the metal.
• This suggests that the proportionality may be converted to an equality by introducing a coefficient that is a
measure of the material behavior.
∆𝑻
𝒒𝒙 = 𝒌𝑨
∆𝒙
• where 𝑘, the thermal conductivity (𝑊/𝑚 𝐾) is an important property of the material.
• Evaluating this expression in the limit as ∆𝑥 → 0, we obtain for the heat rate
𝒅𝑻
𝒒𝒙 = −𝒌𝑨
𝒅𝒙 Asad Akhter Naqvi
The Conduction Rate Equation (cont.)
• Heat flux can be given as
′′
𝒒𝒙 𝒅𝑻
𝒒𝒙 = = −𝒌
𝑨 𝒅𝒙
• The minus sign is necessary because heat is always transferred in the direction of
decreasing temperature.
• Above equation, implies that the heat flux is a directional quantity.
• In particular, the direction of 𝒒𝒙 ′′ is normal to the cross-sectional area 𝑨.
• More generally, the direction of heat flow will always be normal to a surface of
constant temperature, called an isothermal surface.
• Note that the isothermal surfaces are planes normal to the 𝑥-direction.
• Recognizing that the heat flux is a vector quantity, we can write a more general
statement of the conduction rate equation (Fourier's law ) as follows:
′′
𝝏𝑻 𝝏𝑻 𝝏𝑻
𝒒 = −𝒌𝜵𝑻 = −𝒌 𝒊 +𝒋 +𝒌
𝝏𝒙 𝝏𝒚 𝝏𝒛
• where 𝜵 is the three-dimensional del operator and T(x, y, z) is the scalar temperature
field.
Asad Akhter Naqvi
The Conduction Rate Equation (cont.)
• The heat flux vector can be resolved into components such that, in Cartesian coordinates, the general expression
for 𝒒′′ is
′′ ′′ ′′ ′′
𝒒 = 𝒊𝒒 𝒙
+ 𝒋𝒒 𝒚
+ 𝒌𝒒 𝒛
• Here
′′
𝝏𝑻
𝒒 𝒙 = −𝒌
𝝏𝒙
′′
𝝏𝑻
𝒒 𝒚 = −𝒌
𝝏𝒚
′′
𝝏𝑻
𝒒 𝒛 = −𝒌
𝝏𝒛
Asad Akhter Naqvi
The Thermal Properties of Matter
Thermal Conductivity
• From Fourier’s law, the thermal conductivity associated with conduction in the 𝑥-direction is defined as:
𝒒′′ 𝒙
𝒌𝒙 = −
𝝏𝑻
𝝏𝒙
• Similar definitions are associated with thermal conductivities in the y- and z-directions (𝒌𝒚 , 𝒌𝒛 ),
• But for an isotropic material the thermal conductivity is independent of the direction of transfer,
𝒌𝒙 = 𝒌𝒚 = 𝒌𝒛 = 𝒌
Differential control volume, 𝒅𝒙, 𝒅𝒚, 𝒅𝒛, for conduction analysis in Cartesian coordinates.
Asad Akhter Naqvi
The Heat Diffusion Equation (cont.)
• Consider a homogeneous medium within which there is no bulk motion (advection) and the temperature
distribution T(x, y, z) is expressed in Cartesian coordinates.
• In the absence of motion (or with uniform motion), there are no changes in mechanical energy and no work being
done on the system.
• If there are temperature gradients, conduction heat transfer will occur across each of the control surfaces.
• For an infinitesimally small (differential) control volume, 𝒅𝒙, 𝒅𝒚, 𝒅𝒛, the conduction heat rates perpendicular to
each of the control surfaces at the x-, y-, and z-coordinate locations are indicated by the terms 𝒒𝒙 , 𝒒𝒚 and 𝒒𝒛
respectively.
• The conduction heat rates at the opposite surfaces can then be expressed as a Taylor series expansion where,
neglecting higher-order terms,
𝝏𝑻
𝒒𝒙 = −𝒌𝒅𝒚𝒅𝒛
𝝏𝒙 𝑧 𝑑𝑧
𝝏𝑻 𝑦 𝑑𝑧
𝒒𝒚 = −𝒌𝒅𝒙𝒅𝒛 𝑑𝑥 𝑞𝑥 𝑑𝑧
𝝏𝒚 𝑥
𝝏𝑻 𝑑𝑦
𝑑𝑦
𝒒𝒛 = −𝒌𝒅𝒙𝒅𝒚 𝒒𝒚 𝑑𝑥
𝝏𝒛
Asad Akhter Naqvi
The Heat Diffusion Equation (cont.)
𝝏 𝝏𝑻 𝝏 𝝏𝑻 𝝏 𝝏𝑻 𝝏𝑻
− −𝒌𝒅𝒚𝒅𝒛 𝒅𝒙 − −𝒌𝒅𝒛𝒅𝒚 𝒅𝒚 − −𝒌𝒅𝒙𝒅𝒚 𝒅𝒛 + 𝒒𝒅𝒙𝒅𝒚𝒅𝒛 = 𝝆𝑪𝑷 𝒅𝒙𝒅𝒚𝒅𝒛
𝝏𝒙 𝝏𝒙 𝝏𝒚 𝝏𝒚 𝝏𝒛 𝝏𝒛 𝝏𝒕
𝝏 𝝏𝑻 𝝏 𝝏𝑻 𝝏 𝝏𝑻 𝝏𝑻
𝒌 + 𝒌 + (𝒌 ) + 𝒒 = 𝝆𝑪𝑷
𝝏𝒙 𝝏𝒙 𝝏𝒚 𝝏𝒚 𝝏𝒛 𝝏𝒛 𝝏𝒕
This is the general form of the heat diffusion equation in Cartesian coordinates.
𝝏 𝝏𝑻 𝝏 𝝏𝑻 𝝏 𝝏𝑻
𝒌 + 𝒌 + (𝒌 ) + 𝒒 = 𝟎
𝝏𝒙 𝝏𝒙 𝝏𝒚 𝝏𝒚 𝝏𝒛 𝝏𝒛
• Moreover, if the heat transfer is one-dimensional (e.g., in the x-direction) and there is no energy generation.
𝒅 𝒅𝑻
𝒌 =𝟎
𝒅𝒙 𝒅𝒙
• Under steady-state, one-dimensional conditions with no energy generation, the heat flux is a constant in the
direction of transfer. ′′
𝒅𝒒 𝒙
=𝟎
𝒅𝒙 Naqvi
Asad Akhter
The Heat Diffusion Equation in cylinder coordinate
𝝏𝑻 𝟏 𝝏𝑻 𝝏𝑻
𝒒′′ = −𝒌𝜵𝑻 = −𝒌 𝒊 +𝒋 +𝒌
𝝏𝒓 𝒓 𝝏∅ 𝝏𝒛
• Where,
′′ 𝝏𝑻 ′′
𝟏 𝝏𝑻 ′′
𝝏𝑻
𝒒𝒓 = −𝒌 𝒒∅ = −𝒌 𝒒𝒛 = −𝒌
𝝏𝒓 𝒓 𝝏∅ 𝝏𝒛
• General form of the heat equation is:
𝝏𝑻 𝟏 𝝏𝑻 𝟏 𝝏𝑻
𝒒′′ = −𝒌𝜵𝑻 = −𝒌 𝒊 +𝒋 +𝒌
• Where, 𝝏𝒓 𝒓 𝝏𝜽 𝒓𝒔𝒊𝒏𝜽 𝝏∅
′′ 𝝏𝑻 ′′
𝟏 𝝏𝑻 ′′
𝒌 𝝏𝑻
𝒒𝒓 = −𝒌 𝒒𝜽 = −𝒌 𝒒∅ =−
𝝏𝒓 𝒓 𝝏𝜽 𝒓𝒔𝒊𝒏𝜽 𝝏∅
• General form of the heat equation is:
𝐴 𝑥 = (1 − 𝑥)
𝑑𝑇 𝟐𝟎
= 300(−2 − 3𝑥 2 ) 𝒌=
𝑑𝑥 𝟏 − 𝒙 (𝟐 + 𝟑𝒙𝟐 )
𝑑𝑇
Substituting the values of 𝑇 , 𝐴 and in eq. 1
𝑑𝑥
2-11. Consider steady-state conditions for one-dimensional conduction in a plane
wall having a thermal conductivity 𝑘 = 50𝑊/𝑚 𝐾 and a thickness 𝐿 = 0.25 𝑚,
with no internal heat generation. Determine the heat flux and the unknown
quantity for each case and sketch the temperature distribution, indicating the
direction of the heat flux.
Given data
𝐶𝑎𝑠𝑒 𝑇1 𝑇2 𝑑𝑇/𝑑𝑥
𝐶𝑎𝑠𝑒 𝑇1 𝑇2 𝑑𝑇/𝑑𝑥 𝑞
1 50 −20
1 50 −20 −280 14000
2 −30 −10
2 −30 −10 80 −4000
3 70 160
4 40 −80 3 70 110 160 −8000
5 30 200 4 60 40 −80 4000
5 −20 30 200 −10000
𝐿 = 0.25 𝑚
Using Fourier’s law of heat transfer
′′
𝒅𝑻
𝒒 = −𝒌
𝒅𝒙
𝒅𝑻 𝑻𝟐 − 𝑻𝟏
=
𝒅𝒙 𝑳
2-26. At a given instant of time, the temperature distribution within an infinite
homogeneous body is given by the function:
𝝏 𝝏𝑻 𝝏 𝝏𝑻 𝝏 𝝏𝑻 𝝏𝑻
𝒌 + 𝒌 + 𝒌 ሶ
+ 𝒒 = 𝝆𝑪𝑷
𝝏𝒙 𝝏𝒙 𝝏𝒚 𝝏𝒚 𝝏𝒛 𝝏𝒛 𝝏𝒕
𝜕𝑇
𝑘 2 − 4 + 2 = 𝜌𝐶𝑃
𝜕𝑡
𝝏𝑻
=𝟎 There is no region where Temperature changes with time
𝝏𝒕
2-30. The steady-state temperature distribution in a one dimensional wall of
thermal conductivity 50 𝑊/𝑚 𝐾 and thickness 50 𝑚𝑚 is observed to be 𝑇 𝐶 =
𝑎 + 𝑏𝑥 2 , where 𝑎 = 200𝑜 𝐶, 𝑏 = −2000𝑜 𝐶/𝑚2 , and 𝑥 is in meters.
(a) What is the heat generation rate in the wall?
(b) Determine the heat fluxes at the two wall faces. In what manner are these heat
fluxes related to the heat generation rate?
Using heat diffusion equation
𝝏 𝝏𝑻 𝝏 𝝏𝑻 𝝏 𝝏𝑻 𝝏𝑻
Face 1 Face 2 𝒌 + 𝒌 + 𝒌 + 𝒒ሶ = 𝝆𝑪𝑷
𝝏𝒙 𝝏𝒙 𝝏𝒚 𝝏𝒚 𝝏𝒛 𝝏𝒛 𝝏𝒕
Since there is only one dimensional flow and no energy storage so,
𝒅 𝒅𝑻
𝒌 + 𝒒ሶ = 𝟎
𝒅𝒙 𝒅𝒙
50 𝑚𝑚 𝒅𝟐 𝑻
𝒌 + 𝒒ሶ = 𝟎
𝑥 𝒅𝒙𝟐
Heat Flux will be given by using Fourier’s law of Heat transfer
𝑻=𝒂+ 𝒃𝒙𝟐
𝒅𝑻
𝒒′′ = −𝒌
𝑻 = 𝟐𝟎𝟎 − 𝟐𝟎𝟎𝟎𝒙𝟐 𝒅𝒙
𝟐
𝒅 𝑻 𝒒′′ = −𝒌(−𝟒𝟎𝟎𝟎𝒙)
= −𝟒𝟎𝟎𝟎
𝒅𝒙𝟐
At face 1 𝒙 = 𝟎
Now we have,
𝒒′′ = 𝟎
𝒅𝟐 𝑻
𝒒ሶ = −𝒌
𝒅𝒙𝟐 At face 2 𝒙 = 𝟓𝟎 𝒎𝒎
𝒅𝑻
𝒒𝒊𝒏 𝒒𝒐𝒖𝒕 𝒒 = −𝒌𝑨
𝒅𝒙
𝑬𝒔𝒕
𝒅𝑻
𝒒𝒊𝒏 = −𝒌𝑨 (𝒙 = 𝟎)
𝒅𝒙
0.3 𝑚 𝒅𝑻
𝑥 𝒒𝒐𝒖𝒕 = −𝒌𝑨 (𝒙 = 𝑳)
𝒅𝒙
The temperature distribution is Now we are going to find the temperature at both sides
𝑻 𝒙 = 𝒂 + 𝒃𝒙 + 𝒄𝒙𝟐 𝑻 𝒙 = 𝟎 = 𝟐𝟎𝟎𝒐 𝑪
𝑇2 = 142.7𝑜 𝐶
𝑻 𝒙 = 𝟐𝟎𝟎 − 𝟐𝟎𝟎𝒙 + 𝟑𝟎𝒙𝟐 𝑻 𝒙 = 𝟎. 𝟑 = 𝟏𝟒𝟐. 𝟕𝒐 𝑪
𝑇1 = 200𝑜 𝐶
𝒅𝑻
= −𝟐𝟎𝟎 + 𝟔𝟎𝒙
𝒅𝒙
𝒅𝑻
The heat gained by the fluid is 𝒒𝒐𝒖𝒕
𝒒𝒊𝒏 = −𝒌𝑨 (𝒙 = 𝟎) given by Newton’s law of cooling
𝒅𝒙
𝒒𝒊𝒏 = −𝟏 × 𝟏 × (−𝟐𝟎𝟎 + 𝟔𝟎 × 𝟎) 𝒒𝒐𝒖𝒕 = 𝒉𝑨(𝑻𝒔 − 𝑻∞ )
𝑇 𝑟=0 =𝑎
𝑇 𝑟 = 0.025 𝑚 = 𝑎 + 0.00625𝑏
𝒒′ (𝒓 = 𝟎) = 𝟎 𝝏𝑻
= 𝟓𝟔. 𝟖𝟐 𝑲/𝒔
𝝏𝒕
𝒒′ 𝒓 = 𝟎. 𝟎𝟐𝟓 = −𝟑𝟎 × 𝟐𝝅 × 𝟎. 𝟎𝟐𝟓 × −𝟖. 𝟑𝟑 × 𝟏𝟎𝟓 × 𝟎. 𝟎𝟐𝟓
𝑇∞ = 20𝑜 𝐶 (𝑥 = −𝐿) 𝑥 (𝑥 = 𝐿)
𝑇∞ = 20𝑜 𝐶
ℎ(𝑥=−𝐿) =? (𝑥 = 0) ℎ(𝑥=𝐿) =?
Part a
𝑻 𝒙 = 𝒂 + 𝒃𝒙 + 𝒄𝒙𝟐
𝑎 = 82.0𝑜 𝐶
𝑏 = −210𝑜 𝐶/𝑚
𝑐 = −2 × 104 𝐶/𝑚2
𝑷𝒂𝒓𝒕 𝒃 𝑷𝒂𝒓𝒕 𝒄
Using One Dimensional Heat Diffusion ′′
𝒅𝑻
𝒒 = −𝒌
in Rectangular Coordinate system 𝒅𝒙
𝒅 𝒅𝑻 𝒒ሶ
+ =𝟎 𝐴𝑡 𝑥 = −𝐿 = −0.02m
𝒅𝒙 𝒅𝒙 𝒌
𝒅𝑻 𝒅𝑻
𝒒𝒄𝒐𝒏𝒗 ′′ 𝒙 = −𝑳 = 𝒒′′ (𝒙 = −𝑳) 𝒒′′ = −𝒌 = −210 − 4 × 104 𝑥
𝒅𝒙 𝒅𝒙
𝒒′′ 𝒙 = −𝑳 = 𝐡(𝑻𝒔 − 𝑻∞ )
𝒒′′ = −𝒌(−𝟐𝟏𝟎 − 𝟒 × 𝟏𝟎𝟒 𝒙)
𝐴𝑡 𝑥 = −𝐿 = −0.02m
The location where heat flux is zero (𝒒′′ = 𝟎)
𝑻𝒔 = 𝟕𝟖. 𝟐𝒐 𝑪
𝒉 = 𝟏𝟎𝟏. 𝟒 𝐖/𝒎𝟐 𝑲
2-34 One-dimensional, steady-state conduction with uniform internal energy generation occurs
in a plane wall with a thickness of 50 𝑚𝑚 and a constant thermal conductivity of 5 𝑊/𝑚 𝐾. For
these conditions, the temperature distribution has the form 𝑇 𝑥 = 𝑎 + 𝑏𝑥 + 𝑐𝑥 2 . The surface
at 𝑥 = 0 has a temperature of 𝑇 0 = 𝑇𝑜 = 120𝐶 and experiences convection with a fluid for
which 𝑇 = 20𝑜 𝐶 and ℎ = 500 𝑊/𝑚2 𝐾. The surface at 𝑥 = 𝐿 is well insulated.
a. Applying an overall energy balance to the wall, calculate the volumetric energy generation
rate .
b. Determine the coefficients a, b, and c by applying the boundary conditions to the prescribed
temperature distribution. Use the results to calculate and plot the temperature distribution.
c. Consider conditions for which the convection coefficient is halved, but the volumetric energy
generation rate remains unchanged. Determine the new values of a, b, and c, and use the
results to plot the temperature distribution. Hint: recognize that 𝑇(0) is no longer 120𝐶.
d. Under conditions for which the volumetric energy generation rate is doubled, and the
convection coefficient remains unchanged (ℎ = 500 𝑊/𝑚2 𝐾), determine the new values of
a, b, and c and plot the corresponding temperature distribution.
𝑻∞ = 𝟐𝟎𝒐 𝑪
𝒒𝒊𝒏 𝒒𝒐𝒖𝒕
𝒙
Part a
500 20 − 120
𝑇 𝑥 = 𝑎 + 𝑏𝑥 + 𝑐𝑥 2 𝑞ሶ = −
0.05
Apply Energy Balance Equation
𝐸𝑖𝑛 − 𝐸𝑜𝑢𝑡 + 𝐸𝑔𝑒𝑛 = 𝐸𝑠𝑡 𝒒ሶ = 𝟏 × 𝟏𝟎𝟔 𝑾/𝒎𝟑
Steady State conditions are mentioned Part b
and at 𝑥 = 𝐿 wall is perfectly insulated so
𝑇 𝑥 = 𝑎 + 𝑏𝑥 + 𝑐𝑥 2
𝐸𝑜𝑢𝑡 = 0
𝑇 𝑥 = 0 = 𝑇𝑜 = 𝑎 + 𝑏(0) + 𝑐(0)
𝐸𝑠𝑡 = 0
𝒅𝑻
𝐸𝑖𝑛 + 𝐸𝑔𝑒𝑛 = 0 𝒂 = 𝟏𝟐𝟎𝐨 𝐂 𝒒𝒐𝒖𝒕 = −𝒌𝑨 𝒙=𝑳 =𝟎
𝒅𝒙
𝐸𝑖𝑛 = −𝐸𝑔𝑒𝑛
𝒅𝑻 𝒅𝑻
𝒒𝒊𝒏 = 𝒉𝑨 𝑻∞ − 𝑻𝑺 = −𝒌𝑨 (𝒙 = 𝟎) = 𝒃 + 𝟐𝒄𝑳 = 𝟎
At 𝑥 = 0, heat transfer takes place 𝒅𝒙 𝒅𝒙
between fluid and the heated surface
𝒅𝑻 𝒃 𝟏 × 𝟏𝟎𝟒
𝐸𝐸𝑖𝑛 =
= 𝑞𝐴𝐿
ሶ
𝑞 = ℎ𝐴(𝑇∞ − 𝑇𝑠 ) 𝒅𝒙
= 𝒃 + 𝟐𝒄𝒙 𝒄=− =
𝟐𝑳 𝟐 × 𝟎. 𝟎𝟓
𝑔𝑒𝑛
𝑞𝐴𝐿
ሶ = −ℎ𝐴(𝑇∞ − 𝑇𝑠 ) 𝟓𝟎𝟎 × 𝟐𝟎 − 𝟏𝟐𝟎 = −𝟓(𝒃 + 𝟐𝒄 × 𝟎)
ℎ 𝑇∞ − 𝑇𝑠 𝟓𝒐
𝑞ሶ = − 𝒃 = 𝟏 × 𝟏𝟎 𝑪/𝒎 𝟒𝒐 𝒄 = −𝟏 × 𝟏𝟎 𝑪/𝒎𝟐
𝐿
Part c
Now convection coefficient is halved but
volumetric heat generation remains same.
𝐸𝑖𝑛 = −𝐸𝑔𝑒𝑛
ℎ 𝑇∞ − 𝑇𝑠 = −𝑞𝐿
ሶ
250 × 20 − 𝑇𝑠 = −1 × 106 × 0.05
𝑇𝑠 = 𝑇𝑜 = 220𝑜 𝐶 𝑑𝑇
𝑞𝑜𝑢𝑡 = −𝑘𝐴 𝑥=𝐿 =0
𝑇 𝑥 = 0 = 𝑇𝑜 = 𝑎 + 𝑏(0) + 𝑐(0) 𝑑𝑥
𝑑𝑇
𝒂 = 𝟐𝟐𝟎𝒐 𝑪 = 𝑏 + 2𝑐𝐿 = 0
𝑑𝑥
𝑑𝑇
𝑞𝑖𝑛 = ℎ𝐴 𝑇∞ − 𝑇𝑆 = −𝑘𝐴 (𝑥 = 0) 𝑏 1 × 104
𝑑𝑥 𝑐=− =
𝑑𝑇 2𝐿 2 × 0.05
= 𝑏 + 2𝑐𝑥
𝑑𝑥
250 × 20 − 220 = −5(𝑏 + 2𝑐 × 0) 𝟓𝒐
𝒄 = −𝟏 × 𝟏𝟎 𝑪/𝒎𝟐
𝟒𝒐
𝒃 = 𝟏 × 𝟏𝟎 𝑪/𝒎
Part d 𝑑𝑇
𝑞𝑜𝑢𝑡 = −𝑘𝐴 𝑥=𝐿 =0
Now convection coefficient remains same 𝑑𝑥
but volumetric heat generation gets
𝑑𝑇
doubled. = 𝑏 + 2𝑐𝐿 = 0
𝑑𝑥
𝐸𝑖𝑛 = −𝐸𝑔𝑒𝑛
ℎ 𝑇∞ − 𝑇𝑠 = −𝑞𝐿
ሶ 𝑏 1 × 104
𝑐=− =
2𝐿 2 × 0.05
500 × 20 − 𝑇𝑠 = −2 × 106 × 0.05
𝑇𝑠 = 𝑇𝑜 = 220𝑜 𝐶
𝟓𝒐
𝑇 𝑥 = 0 = 𝑇𝑜 = 𝑎 + 𝑏(0) + 𝑐(0) 𝒄 = −𝟏 × 𝟏𝟎 𝑪/𝒎𝟐
𝒂 = 𝟐𝟐𝟎𝒐 𝑪
𝑑𝑇
𝑞𝑖𝑛 = ℎ𝐴 𝑇∞ − 𝑇𝑆 = −𝑘𝐴 (𝑥 = 0)
𝑑𝑥
𝑑𝑇
= 𝑏 + 2𝑐𝑥
𝑑𝑥
250 × 20 − 220 = −5(𝑏 + 2𝑐 × 0)
𝟒𝒐
𝒃 = 𝟏 × 𝟏𝟎 𝑪/𝒎
2-46. A steam pipe is wrapped with insulation of inner and outer radii 𝑟𝑖 and 𝑟𝑜 ,
respectively. At a particular instant the temperature distribution in the insulation is
known to be of the form
𝒓
𝑻 𝒓 = 𝑪𝟏 𝐥𝐧 + 𝑪𝟐
𝒓𝒐
Are conditions steady-state or transient? How do the heat flux and heat rate vary
with radius?
Using Heat diffusion equation in cylindrical
coordinate system
𝝏𝑻
=𝟎
𝝏𝒕
1 𝜕 𝜕𝑇 𝜕𝑇 Eq. 1
𝑘𝑟 = 𝜌𝐶𝑃
𝑟 𝜕𝑟 𝜕𝑟 𝜕𝑡
There is no change in temperature with
respect to time so the conditions are steady.
Since there is no internal heat generation
and temperature distribution is only in
radial direction.
𝑇 𝑟 = 𝐶1 ln 𝑟 − 𝐶1 𝑙𝑛𝑟𝑜 + 𝐶2
𝜕𝑇 𝐶1
=
𝜕𝑟 𝑟
𝑘 𝜕 𝐶1 𝜕𝑇
𝑟× = 𝜌𝐶𝑃
𝑟 𝜕𝑟 𝑟 𝜕𝑡
𝑘 𝜕 𝜕𝑇
𝐶1 = 𝜌𝐶𝑃
𝑟 𝜕𝑟 𝜕𝑡
2-49. Two-dimensional, steady-state conduction occurs in a hollow cylindrical solid
of thermal conductivity 𝑘 = 16 𝑊/𝑚 𝐾, outer radius 𝑟0 = 1𝑚 and overall length
2𝑧𝑜 = 5𝑚, where the origin of the coordinate system is located at the midpoint of
the center line. The inner surface of the cylinder is insulated, and the temperature
distribution within the cylinder has the form 𝑇 𝑟, 𝑧 = 𝑎 + 𝑏𝑟2 + 𝑐𝑙𝑛𝑟 + 𝑑𝑧 2 ,
where 𝑎 = 20𝑜 𝐶, 𝑏 = 150𝑜 𝐶/𝑚2 , 𝑐 = 12𝑜 𝐶, 𝑑 = 300𝑜 𝐶/𝑚2 and r and z are in
meters.
a. Determine the inner radius 𝑟𝑖 of the cylinder.
b. Obtain an expression for the volumetric rate of heat generation (𝑞ሶ 𝑊/𝑚 ሶ 3 ),
c. Determine the axial distribution of the heat flux at the outer surface 𝑞𝑟 ′′ (𝑟𝑜 , 𝑧).
What is the heat rate at the outer surface? Is it into or out of the cylinder?
d. Determine the radial distribution of the heat flux at the end faces of the
cylinder 𝑞𝑧 ′′ (𝑟, +𝑧𝑜 ) and 𝑞𝑧 ′′ 𝑟, −𝑧𝑜 . What are the corresponding heat rates?
Are they into or out of the cylinder?
e. Verify that your results are consistent with an overall energy balance on the
cylinder.
𝒓𝒊
𝒓𝒐
𝝏𝑻
𝒒𝒓 ′′ = −𝒌
𝝏𝒓
𝑻 𝒓, 𝒛 = −𝟐𝟎 + 𝟏𝟓𝟎𝒓𝟐 − 𝟏𝟐𝒍𝒏𝒓 − 𝟑𝟎𝟎𝒛𝟐 Steady State conditions are given and temperature gradient is only
in radial and axial direction so
𝝏𝑻 𝟏𝟐
= 𝟑𝟎𝟎𝒓 − 𝟏 𝝏 𝝏𝑻 𝝏 𝝏𝑻
𝝏𝒓 𝒓 𝒌𝒓 + 𝒌 + 𝒒ሶ = 𝟎
𝒓 𝝏𝒓 𝝏𝒓 𝝏𝒛 𝝏𝒛
𝑨𝒕 𝒓 = 𝒓𝒊 𝒒𝒓 ′′ = 𝟎 𝝏𝑻 𝟏𝟐 𝝏𝑻
= 𝟑𝟎𝟎𝒓 − = −𝟔𝟎𝟎𝒛
𝝏𝒓 𝒓 𝝏𝒛
𝝏𝑻
𝟎 = −𝒌 (𝒓 = 𝒓𝒊 )
𝝏𝒓 𝒌 𝝏 𝟏𝟐 𝝏
𝒓 𝟑𝟎𝟎𝒓 − +𝒌 −𝟔𝟎𝟎𝒛 + 𝒒ሶ = 𝟎
𝟏𝟐 𝒓 𝝏𝒓 𝒓 𝝏𝒛
𝟑𝟎𝟎𝒓𝒊 − =𝟎
𝒓𝒊
𝒌 𝝏 𝟐
𝝏
𝟑𝟎𝟎𝒓 − 𝟏𝟐 + 𝒌 −𝟔𝟎𝟎𝒛 + 𝒒ሶ = 𝟎
𝒓𝒊 = 𝟎. 𝟐 𝒎 𝒓 𝝏𝒓 𝝏𝒛
𝒌
𝒓
𝟔𝟎𝟎𝒓 − 𝒌(𝟔𝟎𝟎) + 𝒒ሶ = 𝟎 𝒒ሶ = 𝟎
Part c
𝝏𝑻
𝒒𝒓 ′′ = −𝒌 (𝒓 = 𝒓𝒐 )
𝝏𝒓
𝝏𝑻 𝟏𝟐
= 𝟑𝟎𝟎𝒓 −
𝝏𝒓 𝒓
𝑨𝒕 𝒓 = 𝒓𝒐 = 𝟏𝒎
𝝏𝑻
= 𝟐𝟖𝟖
𝝏𝒓
𝝏𝑻 𝒓𝒐
= −𝟔𝟎𝟎𝒛
𝝏𝒛
𝑨𝒕 𝒛 = +𝒛𝒐 = 𝟐. 𝟓𝒎
The area perpendicular to the axial heat transfer = 𝝅𝒓𝟐𝒐 − 𝝅𝒓𝟐𝒊 = 𝟑. 𝟎𝟏𝟓 𝒎𝟐
𝑻∞,𝟏 𝑻∞,𝟐
𝑻𝑺,𝟐
𝑻𝑺,𝟏
𝒉𝟏 𝒉𝟐
𝒙
The Plane Wall
• A plane wall separates two fluids of different temperatures.
• The temperature distribution in the wall can be determined by solving the
heat equation with the proper boundary conditions.
• For steady state with no internal heat generation and for one dimensional
system: 𝒅 𝒅𝑻
𝒌 =𝟎
𝒅𝒙 𝒅𝒙
• If the thermal conductivity of the wall material is assumed to be constant,
the equation may be integrated twice to obtain the general solution.
𝑻 𝒙 = 𝑪𝟏 𝒙 + 𝑪𝟐
• To obtain the constants of integration, 𝑪𝟏 and 𝑪𝟐 , boundary conditions must
be introduced.
• We choose to apply conditions of the first kind at 𝑥 = 0 and 𝑥 = 𝐿, in which
case.
• Applying the condition at 𝑥 = 0 to the general solution, it follows that
𝑻𝒔,𝟏 = 𝑪𝟐
The Plane Wall (Cont.)
• Similarly, at 𝑥 = 𝐿, 𝑻𝒔,𝟐 = 𝑪𝟏 𝑳 + 𝑪𝟐 = 𝑪𝟏 𝑳 + 𝑻𝒔,𝟏
𝑻𝒔,𝟐 − 𝑻𝒔,𝟏
𝑪𝟏 =
𝑳
• Substituting into the general solution, the temperature distribution is then
𝑻𝒔,𝟐 −𝑻𝒔,𝟏
𝑻 𝒙 = 𝒙+ 𝑻𝒔,𝟏
𝑳
• Now that we have the temperature distribution, we may use Fourier’s law, to determine the conduction heat
transfer rate.
𝒅𝑻 𝑻𝒔,𝟏 −𝑻𝒔,𝟐
𝒒𝒙 = −𝒌𝑨 𝒅𝒙 = 𝒌𝑨 𝑳
• Note that 𝑨 is the area of the wall normal to the direction of heat transfer and, for the plane wall, it is a constant
independent of 𝑥.
• The heat flux is then 𝒒𝒙 𝑻𝒔,𝟏 − 𝑻𝒔,𝟐
𝒒𝒙 ′′ = =𝒌
𝑨 𝑳
• The above equations show that both heat transfer and heat flux are independent of 𝑥.
Thermal Resistance
• An analogy exists between the diffusion of heat and electrical charge.
• Just as an electrical resistance is associated with the conduction of electricity, a thermal resistance may be associated with the
conduction of heat.
• Defining resistance as the ratio of a driving potential to the corresponding transfer rate.
• It follows from that the thermal resistance for conduction in a plane wall is
𝑻𝒔,𝟏 −𝑻𝒔,𝟐
𝒒𝒙 = 𝒌𝑨
𝑳
• For electrical conduction in the same system, Ohm’s law provides an electrical resistance of the form
• A thermal resistance may also be associated with heat transfer by convection at a surface.
• From Newton’s law of cooling,
𝒒 = 𝒉𝑨(𝑻𝒔 − 𝑻∞ )
In terms of the overall temperature difference, 𝑻∞,𝟏 −𝑻∞,𝟐 , and the total thermal resistance, 𝑅𝑡𝑜𝑡 , the heat transfer rate may also
be expressed as
𝑻∞,𝟏 − 𝑻∞,𝟐 𝟏 𝑳 𝟏
𝒒𝒙 = 𝑹𝒕𝒐𝒕 = + +
𝑹 𝒉𝟏 𝑨 𝒌𝑨 𝒉𝟐 𝑨
𝒕𝒐𝒕
The Composite Wall 𝑻∞,𝟏 > 𝑻∞,𝟐
𝑻∞,𝟏
𝑻𝑺,𝟏
𝒉𝟏
𝑻𝟐
𝑻𝟑
𝑻𝑺,𝟒
𝑘𝐴 𝑘𝐵 𝑘𝐶 𝑻∞,𝟒
𝐿𝐴 𝐿𝐵 𝐿𝐶
𝒉𝟒
𝒙
The Composite Wall
➢ The one dimensional heat transfer rate for this system may be expressed as
𝑻∞,𝟏 − 𝑻∞,𝟒
𝒒𝒙 =
σ 𝑹𝒕
➢ where 𝑻∞,𝟏 − 𝑻∞,𝟒 is the overall temperature difference, and the summation
includes all thermal resistances. Hence
➢ Alternatively, the heat transfer rate can be related to the temperature difference
and resistance associated with each element. For example,
𝑇∞,1 − 𝑇𝑠,1 𝑇𝑠,1 − 𝑇𝑠,2 𝑇𝑠,2 − 𝑇𝑠,3 𝑇𝑠,3 − 𝑇𝑠,4 𝑇𝑠,2 − 𝑇∞,2
𝑞𝑥 = = = = =
1 𝐿𝐴 𝐿𝐵 𝐿𝑐 1
ℎ1 𝐴 𝑘𝐴 𝐴 𝑘𝐵 𝐴 𝑘𝐶 𝐴 ℎ4 𝐴
➢ With composite systems, it is often convenient to work with an overall heat
transfer coefficient U, which is defined by an expression analogous to Newton’s
law of cooling. Accordingly,
𝒒𝒙 = 𝑼𝑨∆𝑻
➢ ∆𝑻 is the overall temperature difference. The overall heat transfer coefficient U
is related to the total thermal resistance,
The Composite Wall
𝑇∞,1 − 𝑇∞,4 𝟏
𝑞𝑥 =
σ 𝑅𝑡 𝑞𝑥 = 𝑈𝐴∆𝑇 𝑼𝑨 =
σ 𝑹𝒕
𝑻𝒔,𝟏
𝑻∞,𝟐
𝑻∞,𝟏
𝒒 𝑟2
𝑻𝒔,𝟐
𝑟1
𝑳
𝑻∞,𝟏
𝒉𝟏
The Cylinder (Cont.)
➢ For steady-state conditions with no heat generation, the appropriate form of the heat diffusion equation.
1 𝑑 𝑑𝑇
𝑘𝑟 =0
𝑟 𝑑𝑟 𝑑𝑟
➢ Integrating twice to get the general distribution of temperature.
𝑑𝑇 𝑑𝑇 𝐶1
𝑘𝑟
𝑑𝑟
=C = 𝑇 𝑟 = 𝐶1 𝑙𝑛𝑟 + 𝐶2
𝑑𝑟 𝑟
➢ Applying boundary conditions.
𝐴𝑡 𝑟 = 𝑟1 , 𝑇 = 𝑇𝑠,1 and 𝐴𝑡 𝑟 = 𝑟2 , 𝑇 = 𝑇𝑠,2
𝑇𝑠,1 = 𝐶1 𝑙𝑛𝑟1 + 𝐶2 1
𝑇𝑠,2 = 𝐶1 𝑙𝑛𝑟2 + 𝐶2 2
Solving for 𝐶1 𝑎𝑛𝑑 𝐶2 and substituting into the general solution, we then obtain
The Cylinder (Cont.)
➢ The temperature distribution is known and we can use Fourier’s law to determine the heat
transfer rate
➢ Alternatively, Consider a long cylindrical layer (such as a circular pipe) of inner radius 𝑟1 and
outer radius 𝑟2 , length 𝐿, and average thermal conductivity 𝑘
➢ The two surfaces of the cylindrical layer are maintained at constant temperatures 𝑇𝑠,1 and 𝑇𝑠,2 .
➢ There is no heat generation in the layer and the thermal conductivity is constant.
➢ For one-dimensional heat conduction through the cylindrical layer, we have T(r). 𝒒
➢ Then Fourier’s law of heat conduction for heat transfer through the cylindrical layer can be
expressed as
𝑑𝑇
𝑞 = −𝑘𝐴
𝑑𝑟
➢ The area perpendicular to the direction of heat transfer will be 𝐀 = 𝟐𝝅𝒓𝑳
𝑑𝑇
𝑞 = −𝑘 × 2𝜋𝑟𝐿 ×
𝑑𝑟
The Cylinder (Cont.)
𝑟2 𝑇2
𝑞 𝑞
𝑑𝑟 = −𝑘 × 𝑑𝑇 න 𝑑𝑟 = −𝑘 න 𝑑𝑇
2𝜋𝑟𝐿 𝑟1 2𝜋𝑟𝐿 𝑇1
𝑇1 − 𝑇2
𝑞= 𝑟2
ln ൗ𝑟1 ൘
2𝜋𝐿𝑘
➢ From this result it is evident that, for radial conduction in a cylindrical wall, the thermal resistance is of the form.
𝑟2
ln ൗ𝑟1
𝑅𝑡,𝑐𝑜𝑛𝑑 =
2𝜋𝐿𝑘
➢ This definition is arbitrary, and the overall coefficient may also be defined in terms of A4 or
any of the intermediate areas. Note that
Critical Radius of Insulation
➢ Consider a cylindrical pipe of outer radius 𝑟1 whose outer surface temperature 𝑇1 is
maintained constant
➢ The pipe is now insulated with a material whose thermal conductivity is 𝑘 and outer
radius is 𝑟2 .
➢ Heat is lost from the pipe to the surrounding medium at temperature 𝑇∞ , with a
convection heat transfer coefficient h.
➢ The total thermal resistance per unit length of tube is then
𝑟
ln 𝑟 1
′
𝑅𝑡𝑜𝑡 = 𝑖
+ 𝑇1 > 𝑇∞
2𝜋𝑘 2𝜋𝑟ℎ
➢ The rate of heat transfer per unit length of tube is
𝑇1 − 𝑇∞
𝑟 𝑟𝑖 𝑇∞
𝑞′ = ′
𝑅𝑡𝑜𝑡
➢ An optimum insulation thickness would be associated with the value of r that 𝑇1
minimized 𝑞 ′ or maximized 𝑅𝑡𝑜𝑡′ 𝒒
➢ Since this result is always positive, it follows that 𝑟 = 𝑘/ℎ is the insulation radius for which the
total resistance is a minimum, not a maximum
➢ From the above result it makes more sense to think in terms of a critical insulation radius
𝑘
𝑟𝑐𝑟,𝑐𝑦𝑙 =
ℎ
Critical Radius of Insulation
𝑘
𝑟𝑐𝑟,𝑐𝑦𝑙 =
ℎ
➢ Note that the critical radius of insulation depends on the thermal conductivity of the insulation k
and the external convection heat transfer coefficient h.
➢ The rate of heat transfer from the cylinder increases with the addition of insulation for 𝑟 < 𝑟𝑐𝑟 ,
reaches a maximum when 𝑟 = 𝑟𝑐𝑟 , and starts to decrease for 𝑟 > 𝑟𝑐𝑟 .
➢ Thus, insulating the pipe may actually increase the rate of heat transfer from the pipe instead of
decreasing it when 𝑟 < 𝑟𝑐𝑟 .
The Sphere
➢ For steady-state, one-dimensional conditions with no heat generation. The appropriate form of
Fourier’s law is
4𝜋𝑘(𝑇𝑠,1 − 𝑇𝑠,2 )
𝑞𝑟 =
1ൗ − (1ൗ )
𝑟1 𝑟2
➢ Remembering that the thermal resistance is defined as the temperature difference divided by the
heat transfer rate
1 1 1
𝑅𝑡,𝑐𝑜𝑛𝑑 = −
4𝜋𝑘 𝑟1 𝑟2
➢ Critical radius of insulation for a spherical shell
𝑘
𝑟𝑐𝑟,𝑠𝑝ℎ =
2ℎ
Lecture 9
Problems Related to Thermal Network
3-3. The rear window of an automobile is defogged by passing warm air over its inner surface.
If the warm air is at 𝑇∞,𝑖 = 40𝑜 𝐶 and the corresponding convection coefficient is ℎ𝑖 =
30 𝑊/𝑚2 𝐾, what are the inner and outer surface temperatures of 4 𝑚𝑚 thick window glass, if
the outside ambient air temperature is 𝑇∞,2 = −10𝑜 𝐶 and the associated convection coefficient
is ℎ𝑜 = 65 𝑊/𝑚2 𝐾?
𝑇𝑠,2
𝑇𝑠,1
4 𝑚𝑚
𝑇∞,𝑖 𝑇𝑠,1 𝑇𝑠,2 𝑇∞,𝑜
1 𝐿 1
ℎ𝑖 𝑘 ℎ𝑜
Heat flux will be Same heat flux will flow through each resistances so,
𝑇∞,𝑖− 𝑇∞,𝑜 𝑇∞,𝑖− 𝑇𝑠,1 𝑇𝑠,1− 𝑇𝑠,2
𝑞 ′′ = ′′ 𝑞 ′′ = 𝑞 ′′ =
𝑅𝑡𝑜𝑡 1 𝐿
′′
1 𝐿 1 ℎ𝑖 𝑘
𝑅𝑡𝑜𝑡 = + +
ℎ𝑖 𝑘 ℎ𝑜
1 0.004 1
𝑇𝑠,1 = 𝑇∞,𝑖 − 𝑞 ′′ /ℎ𝑖 𝑇𝑠,2 = 𝑇𝑠,1 − 𝑞′′ × 𝐿/𝑘
′′
𝑅𝑡𝑜𝑡 = + +
30 1.4 65
𝑇𝑠,1 = 40 − 970.8/30 𝑇𝑠,2 = 7.64 − 970.8 × 0.004/1.4
′′
𝑅𝑡𝑜𝑡 = 0.0515 𝑚2 𝐾/𝑊
𝒒′′𝟏 𝒒′′𝟐
𝒒′′𝒐
Amount of heat flux through substrate will be Amount of heat flux through film will be Total heat flux will be the sum of heat flux through
film and substrate
𝑇𝑜 − 𝑇1 𝑇𝑜 − 𝑇∞
𝑞′′1 = 𝑞′′2 =
𝐿𝑠 𝐿𝑓 1 𝑞′′𝑜 = 𝑞′′1 + 𝑞′′2
𝑘𝑠 +
𝑘𝑓 ℎ
60 − 30 60 − 20 𝑞′′𝑜 = 1500 + 1333.33
𝑞′′1 = 𝑞′′2 =
0.001/0.05 0.00025 1
+
0.025 50 𝑞′′𝑜 = 2833.33 𝑊/𝑚2
𝑞′′1 = 1500 𝑊/𝑚2 𝑞′′2 = 1333.33 𝑊/𝑚2
3-7. The walls of a refrigerator are typically constructed by sandwiching a layer of insulation
between sheet metal panels. Consider a wall made from fiberglass insulation of thermal
conductivity 𝑘𝑖 = 0.046 𝑊/𝑚 𝐾 and thickness 𝐿𝑖 = 50 𝑚𝑚 and steel panels, each of thermal
conductivity 𝑘𝑝 = 60 𝑊/𝑚 𝐾 and thickness 𝐿𝑃 = 3 𝑚𝑚. If the wall separates refrigerated air at
𝑇∞,𝑖 = 4𝑜 𝐶 from ambient air at 𝑇∞,𝑜 = 25𝑜 𝐶, what is the heat gain per unit surface area?
Coefficients associated with natural convection at the inner and outer surfaces may be
approximated as ℎ𝑖 = ℎ𝑜 = 5 𝑊/𝑚2 𝐾.
𝑇∞,𝑖 = 4𝑜 𝐶
𝑇∞,𝑜 = 25𝑜 𝐶
50 𝑚𝑚
3 𝑚𝑚 3 𝑚𝑚
𝑇∞,𝑖 𝑇∞,𝑜
1/ℎ𝑖 𝐿𝑠 𝐿𝑃 𝐿𝑠 1/ℎ𝑜
ൗ𝑘 ൗ𝑘 ൗ𝑘
𝑠 𝑃 𝑠
′′
𝑇∞,𝑖− 𝑇∞,𝑜
𝑞 = ′′
𝑅𝑡𝑜𝑡
′′
𝑅𝑡𝑜𝑡 =
1 𝐿𝑠
ℎ𝑖 𝑠
𝐿
𝑃
𝐿
+ ൗ𝑘 + 𝑃ൗ𝑘 + 𝑠ൗ𝑘 +
𝑠
1
ℎ𝑜
𝑞 ′′ = 14.12 𝑊/𝑚2
1 0.003 1
′′
𝑅𝑡𝑜𝑡 = + ൗ60 + 0.05ൗ0.046 + 0.003ൗ0.046 +
5 5
′′
𝑅𝑡𝑜𝑡 = 1.48 𝑚2 𝐾/𝑊
25 − 4
𝑞 ′′ =
1.552
3-26. A composite wall separates combustion gases at 2600𝑜 𝐶 from a liquid coolant at 100𝑜 𝐶,
with gas- and liquid-side convection coefficients of 50 and 1000𝑊/𝑚2 𝐾 . The wall is composed
of a 10 𝑚𝑚 thick layer of beryllium oxide on the gas side and a 20 𝑚𝑚 thick slab of stainless
steel (𝐴𝐼𝑆𝐼 304) on the liquid side. The contact resistance between the oxide and the steel is
0.05 𝑚2 𝑊/𝐾. What is the heat loss per unit surface area of the composite? Sketch the
temperature distribution from the gas to the liquid.
10 𝑚𝑚 20 𝑚𝑚
𝐿𝐵 𝐿𝑠
ൗ𝑘 𝑅′′𝑐𝑜𝑛𝑡 ൗ𝑘
𝑇∞,𝑖 𝐵 𝑠
𝑇∞,𝑜
𝑇∞,𝑖 = 2600𝑜 𝐶
𝑇𝑠,1 = 1906𝑜 𝐶
Similarly
𝑇𝑐,1 = 1889𝑜 𝐶
𝑇𝑐,1 = 𝑇𝑠,1 − 𝑞′′ × 𝐿𝐵 /𝑘𝐵
𝑇𝑐,1 = 1906 − 34672 × 0.01/21.5
𝑇𝑐,1 = 1889𝑜 𝐶
𝑇∞,0 = 100𝑜 𝐶
𝑇𝑐,1 = 155.4𝑜 𝐶
Similarly
𝑇𝑠,2 = 134𝑜 𝐶
𝑇𝑐,2 = 𝑇𝑐,1 − 𝑞′′ × 𝑅′′𝑐𝑜𝑛
𝑇𝑐,2 = 1889 − 34672 × 0.05
𝑇𝑐,2 = 155.4𝑜 𝐶
Lecture 10
Problem Related to Thermal Network (Cont.)
𝑞 ′′
𝑞 𝑇∞,𝑖 = 400 𝐾
′′
𝑇∞,𝑜 − 𝑇∞,𝑖 𝑇𝑠,2 = 𝑇∞,𝑖 + ൘ℎ
𝑞 = 𝑖 ℎ𝑖 = 500 𝑊/𝑚2 𝐾
𝑅′′𝑡𝑜𝑡𝑎𝑙
1 𝐿𝑖𝑛 1 𝑇𝑠,2 = 400 + 406250ൗ500 5 𝑚𝑚
𝑅′′𝑡𝑜𝑡𝑎𝑙 = + +
ℎ𝑜 𝑘𝑖𝑛 ℎ𝑖
𝑇𝑠,2 = 1212 𝐾
1 0.005 1
𝑅′′𝑡𝑜𝑡𝑎𝑙 = + +
1000 25 500
′′
𝐿𝑖𝑛
𝑇𝑠,1 = 𝑇𝑠,𝑖 + 𝑞
𝑅′′𝑡𝑜𝑡𝑎𝑙 = 0.0032 𝑚2 𝐾/𝑊 𝑘𝑖𝑛
𝑇𝑠,2 = 1212 + 406250 0.005ൗ25
1700 − 400
𝑞 ′′ = 𝑇𝑠,2 = 1293 𝐾
0.0032
𝑟 𝑇∞ = −10𝑜 𝐶
1/ℎ × 2𝜋𝑟𝑜 𝑅′𝑡,𝑐 ln( 𝑜ൗ𝑟 )
𝑖 ℎ = 100 𝑊/𝑚2 𝐾
2𝜋𝑘
𝒒𝒐 𝑟𝑖 = 25 𝑚𝑚
Heat transfer from the heater to the Heat transfer from the heater to the
inside wall of the cylinder will be fluid will be
𝑇𝑜 − 𝑇∞
𝑇𝑜 − 𝑇𝑖 𝑞′2 =
𝑞′1 = 𝑟 1
ln( 𝑜ൗ𝑟 ) ℎ × 2𝜋 × 𝑟𝑜
𝑖 Total heat transfer will be
𝑅′𝑡,𝑐 +
2𝜋𝑘
25 − (−10)
25 − 5 𝑞′2 = 𝒒′𝒐 = 𝒒′𝟏 + 𝒒′𝟐
1
𝑞′1 =
ln(0.075ൗ0.025) 1 × 2𝜋 × 0.075
𝒒′𝒐 = 𝟕𝟐𝟕 + 1649
0.01 + 2𝜋 × 10
𝒒′𝟐 = 𝟏𝟔𝟒𝟗 𝑾/𝒎 𝒒′𝒐 = 𝟐𝟑𝟕𝟑 𝑾/𝒎
𝒒′𝟏 = 𝟕𝟐𝟕 𝑾/𝒎
Asad Akhter Naqvi
3-51. Superheated steam at 575𝑜 𝐶 is routed from a boiler to the turbine of an electric power
plant through steel tubes (𝑘 = 35 𝑊/𝑚 𝐾) of 300 𝑚𝑚 inner diameter and 30 𝑚𝑚 wall
thickness. To reduce heat loss to the surroundings and to maintain a safe-to-touch outer surface
temperature, a layer of calcium silicate insulation (𝑘 = 0.10 𝑊/𝑚 𝐾) is applied to the tubes,
while degradation of the insulation is reduced by wrapping it in a thin sheet of aluminum having
an emissivity of ε = 0.20. The air and wall temperatures of the power plant are 27𝑜 𝐶.
a. Assuming that the inner surface temperature of a steel tube corresponds to that of the steam
and the convection coefficient outside the aluminum sheet is 6 𝑊/𝑚2 𝐾, what is the
minimum insulation thickness needed to ensure that the temperature of the aluminum does
not exceed 50𝑜 𝐶? What is the corresponding heat loss per meter of tube length?
b. Explore the effect of the insulation thickness on the temperature of the aluminum and the
heat loss per unit tube length.
𝑟2 = 0.18𝑚
𝑇∞ = 27𝑜 𝐶
ℎ = 6 𝑊/𝑚2 𝐾
𝑇𝑖 = 575𝑜 𝐶
575 − 50
= 2𝜋𝑟3 (6 50 − 27 + 0.2 × 5.67 × 10−8 (3234 − 3004 )
0.18ൗ
0.15) 𝑟
ln( + ln( 3ൗ0.18)/2𝜋(0.10)
2𝜋 30
𝒓𝟑 = 𝟎. 𝟑𝟗𝟒 𝐦 = 𝟑𝟗𝟒 𝐦𝐦
𝑡𝑖𝑛 = 𝑟3 − 𝑟2
𝑇𝑖 − 𝑇𝑜
𝑞′𝑐𝑜𝑛𝑑 = 𝑟2
ൗ𝑟 ) 𝑟
1
ln( + ln( 3ൗ𝑟2 )/2𝜋𝑘𝑖𝑛
2𝜋𝑘𝑠𝑡
575 − 50
𝑞′𝑐𝑜𝑛𝑑 =
(0.18ൗ0.15)
ln + ln(0.394ൗ0.18)/2𝜋(0.10)
2𝜋 30
𝒒𝑩 ln(𝑟2ൗ𝑟 ) 1
1
𝑟2 = 100 𝑚𝑚
𝜋𝑘𝐵 ℎ × 𝜋𝑟2
𝑇𝑠,𝑖 = 500 𝐾
𝒒𝑩
𝑻𝒔,𝒐(𝑩)
𝒌𝑩 = 𝟎. 𝟐𝟓 𝑾/𝒎𝑲
ln 0.1ൗ0.05 ln 0.1ൗ0.05
𝑇𝑠,𝑜(𝐴) = 500 − 842 𝑇𝑠,𝑜(𝐵) = 500 − 198
𝜋(2) 𝜋(0.25)
𝑻∞ 𝑻∞
𝑻𝒔 𝑻𝒔
𝑞𝑥 = 𝑞𝑥+𝑑𝑥 + 𝑑𝑞𝑐𝑜𝑛𝑣 1
➢ From Fourier’s law we know that
𝑑𝑇
𝑞𝑥 = −𝑘𝐴𝑐
𝑑𝑥
➢ where 𝐴𝑐 is the cross-sectional area, which may vary with x.
Since the conduction heat rate at 𝑥 + 𝑑𝑥 may be expressed as
𝑑𝑞𝑥
𝑞𝑥+𝑑𝑥 = 𝑞𝑥 + 𝑑𝑥
𝑑𝑥
𝑑𝑇 𝑑 𝑑𝑇
𝑞𝑥+𝑑𝑥 = −𝑘𝐴𝑐 −𝑘 𝐴𝑐 𝑑𝑥
𝑑𝑥 𝑑𝑥 𝑑𝑥 Asad Akhter Naqvi
Fin Equation
➢ The convection heat transfer rate may be expressed as
𝑑𝑞𝑐𝑜𝑛𝑣 = ℎ𝑑𝐴𝑠 (𝑇 − 𝑇∞ )
➢ where 𝑑𝐴𝑠 is the surface area of the differential element.
➢ Eq. 1 can be written as
𝑑𝑇 𝑑𝑇 𝑑 𝑑𝑇
−𝑘𝐴𝑐 = −𝑘𝐴𝑐 −𝑘 𝐴𝑐 𝑑𝑥 + ℎ𝑑𝐴𝑠 (𝑇 − 𝑇∞ )
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑑𝑥
𝑑 𝑑𝑇 ℎ 𝑑𝐴𝑠
𝐴𝑐 − 𝑇 − 𝑇∞ = 0
𝑑𝑥 𝑑𝑥 𝑘 𝑑𝑥
𝑑2 𝑇 𝑑𝐴𝑐 𝑑𝑇 ℎ 𝑑𝐴𝑠
𝐴𝑐 2 + − 𝑇 − 𝑇∞ = 0
𝑑𝑥 𝑑𝑥 𝑑𝑥 𝑘 𝑑𝑥
𝑑2 𝑇 ℎ𝑃
− 𝑇 − 𝑇∞ = 0
𝑑𝑥 2 𝑘𝐴𝑐
Defining
𝜽 𝒙 = 𝑻 𝒙 − 𝑻∞
𝑑𝜽 𝑑𝑇
➢ since 𝑻∞ is a constant, = 𝑥
𝑑𝑥 𝑑𝑥
𝑑2𝑇 ℎ𝑃
− 𝑇 − 𝑇∞ = 0
𝑑𝑥 2 𝑘𝐴𝑐
𝑑2𝜃 ℎ𝑃
− 𝑚2 𝜃 = 0 𝑚2 =
𝑑𝑥 2 𝑘𝐴𝑐
𝜃 0 = 𝜃𝑏 = 𝑇𝑏 − 𝑇∞
➢ Consider case A, where at the tip the amount of conduction
is equal to the convection to the air
𝑞𝑐𝑜𝑛𝑑 = 𝑞𝑐𝑜𝑛𝑣
𝑑𝑇
ℎ𝐴𝑐 𝑇 𝐿 − 𝑇∞ = −𝑘𝐴𝑐 ቤ𝑥 = 𝐿
𝑑𝑥
𝑑𝜃
ℎ𝜃 𝐿 = −𝑘 ቤ𝑥 = 𝐿 3
𝑑𝑥
➢ Now applying boundary conditions
𝜃𝑏 = 𝐶1 + 𝐶2 𝐴 𝒒𝒄𝒐𝒏𝒗
𝒒𝒄𝒐𝒏𝒅
𝜃 𝑥 = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥
𝑑𝜃
ቤ 𝑥 = 𝐿 = 𝑚(𝐶1 𝑒 𝑚𝐿 − 𝐶2 𝑒 −𝑚𝐿 )
𝑑𝑥
𝒒𝒇 = 𝒉𝑷𝒌𝑨𝒄 𝜽𝒃 𝐭𝐚𝐧𝐡 𝒎𝑳
Asad Akhter Naqvi
Fins of Uniform Cross-Sectional Area
➢ Assume Case C, where fin is very long it means
𝑳 → ∞ 𝒕𝒉𝒆𝒏 𝜽𝑳 → 𝟎
𝜽
= 𝒆−𝒎𝒙
𝜽𝒃
𝒒𝒇 = 𝒉𝑷𝒌𝑨𝒄 𝜽𝒃
𝑻∞
𝑻∞
𝑞𝑓𝑖𝑛,𝑚𝑎𝑥 = ℎ𝐴𝑠 𝜃𝑏 𝒙
𝑇(𝑥)
𝑞𝑓 𝑞𝑓
𝜂𝑓 = =
𝑞𝑚𝑎𝑥 ℎ𝐴𝑠 𝜃𝑏
𝐷𝑖 = 1 𝑚
𝒒𝒉𝒆𝒎𝒊 𝑻 = 𝟒𝟎𝟎 𝑲
𝒒𝒉𝒆𝒎𝒊
𝑇∞ = 300 𝐾
ℎ = 10𝑊/𝑚2 𝐾
𝒒𝒄𝒚𝒍𝒊
𝑫𝒐 = 𝑫𝒊 + 𝟐𝒕
𝑫𝒐 = 𝟏 + 𝟐(𝟎. 𝟎𝟐)
𝑫𝒐 = 𝟏. 𝟎𝟒 𝒎
1
𝑅ℎ𝑒𝑚𝑖
𝑇𝑖 ℎ × 2𝜋𝑟𝑜2
𝒒 𝑇𝑜
1
𝑅ℎ𝑒𝑚𝑖 ℎ × 2𝜋𝑟𝑜2
𝑞 = 𝑞𝑐𝑦𝑙 + 2𝑞ℎ𝑒𝑚𝑖𝑠𝑝ℎ𝑒𝑟𝑒
𝒒 = 𝟖𝟔𝟔𝟏 𝑾
𝑞 = 𝑞𝑐𝑦𝑙 + 𝑞𝑠𝑝ℎ𝑒𝑟𝑒
𝑇𝑖 − 𝑇𝑜 𝑇𝑖 − 𝑇𝑜
𝑞= 𝑟 +
ln 𝑜ൗ𝑟𝑖 1 1 1 1
1 − +
2𝜋𝑘𝐿
+
ℎ × 2𝜋𝑟𝑜 𝐿
2𝜋𝑘 𝐷𝑖 𝐷𝑜 ℎ × 4𝜋𝑟𝑜2
𝑻𝒐 = 𝟐𝟎𝒐 𝑪
𝒉 = 𝟔 𝑾/𝒎𝟐 𝑲
𝒅𝒊 = 𝟑𝒎
𝑻𝒊 = −𝟔𝟎𝒐 𝑪
1 1 1 1
−
2𝜋𝑘 𝑑𝑖 𝑑𝑜 ℎ × 4𝜋𝑟𝑜2 𝑘
𝑇∞ − 𝑇𝑖
𝑞=
1 1 1 1
− +
2𝜋𝑘 𝑑𝑖 𝑑𝑜 ℎ × 4𝜋𝑟02 0 − (−60)
612.4 =
1 1 1
20 − (−60) 2𝜋 × 0.06 3 − 𝑑𝑥
𝑞=
1 1 1 1
− +
2𝜋 × 0.06 3 3.5 6 × 4𝜋 × 1.752
𝑑𝑥 = 3.37 m
𝒒 = 𝟔𝟏𝟐. 𝟒 𝑾
𝑟𝑥 = 1.685 m
𝑇𝑜 − 𝑇𝑖
𝑞=
1 1 1
−
2𝜋𝑘 𝑑𝑖 𝑑𝑥
Asad Akhter Naqvi
3-129. A long, circular aluminum rod is attached at one end to a heated wall and
transfers heat by convection to a cold fluid.
a. If the diameter of the rod is tripled, by how much would the rate of heat
removal change?
b. If a copper rod of the same diameter is used in place of the aluminum, by how
much would the rate of heat removal change?
3𝐷
𝒒𝟑𝑫
= 𝟓. 𝟐
𝒒𝑫
Asad Akhter Naqvi
3-130. A brass rod 100 𝑚𝑚 long and 5 𝑚𝑚 in diameter extends horizontally from a
casting at 200𝑜 𝐶. The rod is in an air environment with 𝑇∞ = 20𝑜 𝐶 and ℎ = 30 𝑊/𝑚2 𝐾.
What is the temperature of the rod 25, 50, and 100 𝑚𝑚 from the casting?
𝐿 = 100 𝑚𝑚
𝑇∞ = 20𝑜 𝐶 ℎ = 30 𝑊/𝑚2 𝐾
𝐷 = 5 𝑚𝑚
𝑇𝑏 = 200𝑜 𝐶
𝜃 𝑐𝑜𝑠ℎ 𝑚 𝐿 − 𝑥 + ( ℎΤ𝑚𝑘) sinh 𝑚(𝐿 − 𝑥) 𝜃(𝑥) 𝑐𝑜𝑠ℎ 13.43 0.1 − 𝑥 + ( 30Τ13.43 × 133) sinh 13.43(0.1 − 𝑥)
= =
𝜃𝑏 𝑐𝑜𝑠ℎ 𝑚𝐿 + ( ℎΤ𝑚𝑘) 𝑠𝑖𝑛ℎ 𝑚𝐿 180 𝑐𝑜𝑠ℎ(13.43 × 0.1) + ( 30Τ13.43 × 133) 𝑠𝑖𝑛ℎ(13.43 × 0.1)
ℎ𝑃
𝑚= 𝜃(𝑥) = 𝑇(𝑥) − 𝑇∞
𝑘𝐴𝑐
𝑃 = 𝜋𝐷 = 0.005𝜋
200 + 20
𝑇ത = = 110𝑜 𝐶 𝑥 cosh 𝑚(𝐿 − 𝑥) sinh 𝑚(𝐿 − 𝑥) 𝜃 𝑇 (𝑜 𝐶)
2
At 𝑇ത = 110𝑜 𝐶
0.025 1.55 1.19 136.5 156.5
𝑘𝑏𝑟𝑎𝑠𝑠 = 133 𝑊/𝑚𝐾
0.05 1.24 0.725 108.9 128.9
30 × 0.005𝜋 0.1 1 0 87 107
𝑚= 𝜋
133 × ( 4 × 0.0052 )
𝑚 = 13.43 𝑚−1
𝜃𝑏 = 𝑇𝑏 − 𝑇∞ = 200 − 20 = 180𝑜 𝐶
𝐿 = 10 𝑚𝑚
𝑡 = 1 𝑚𝑚
𝑤≫𝑡
Asad Akhter Naqvi
Case A
sinh 𝑚𝐿 + (ℎΤ𝑚𝑘) cosh 𝑚𝐿
𝑞𝑓 = ℎ𝑃𝑘𝐴𝑐 𝜃𝑏
cosh 𝑚𝐿 + (ℎΤ𝑚𝑘) sinh 𝑚𝐿
𝑃 = 2 𝑤 + 𝑡 = 2𝑤 𝑤
𝑀= ℎ𝑃𝑘𝐴𝑐 𝜃𝑏 = ℎ × 2𝑤 × 𝑘(𝑤𝑡)𝜃𝑏 = 𝑤 2ℎ𝑘𝑡𝜃𝑏 𝑡
𝟏𝟒𝟒. 𝟔
𝜺𝒇 = = 𝟏𝟗. 𝟑
𝟏𝟎𝟎 × 𝟎. 𝟎𝟎𝟏 × 𝟕𝟓 Asad Akhter Naqvi
Asad Akhter Naqvi
Lecture 13
Transient Conduction
𝑇∞ < 𝑇𝑖
𝑞𝑜𝑢𝑡
𝑇∞
𝑇(𝑡)
𝜌𝑉𝑐 𝑑𝜃
= −𝜃
ℎ𝐴𝑠 𝑑𝑡
➢ Separating variables and integrating from the initial condition,
for which 𝑡 = 0 and 𝑇 0 = 𝑇𝑖 , we then obtain
𝑄 = −𝐸𝑠𝑡
𝑇∞
➢ Rearranging, we then obtain
ℎ𝐴𝑠 𝑡
= 𝐵𝑖. 𝐹𝑜
𝜌𝑉𝑐
𝛼𝑡
𝐹𝑜 = 2
𝐿𝑐
𝜃 𝑇 − 𝑇∞
= = exp −𝐵𝑖. 𝐹𝑜
𝜃𝑖 𝑇𝑖 − 𝑇∞
Asad Akhter Naqvi
5-6. Steel balls 12 𝑚𝑚 in diameter are annealed by heating to 1150 𝐾 and then slowly cooling
to 400 𝐾 in an air environment for which 𝑇∞ = 325 𝐾 and ℎ = 20 𝑊/𝑚2 𝐾. Assuming the
properties of the steel to be 𝑘 = 40 𝑊/𝑚 𝐾, ρ = 7800 𝑘𝑔/𝑚3 , and 𝑐 = 600 𝐽/𝑘𝑔 𝐾, estimate
the time required for the cooling process.
𝑻
𝑇∞ = 325 𝐾 1150 𝐾
𝐷 = 12 𝑚𝑚
400 𝐾
𝑇𝑖 = 1150 𝐾
𝒕
𝑡
𝐵𝑖 =
ℎ𝐿𝑐 20 × 0.002
= 𝒕 = 𝟏𝟗 𝒎𝒊𝒏
𝑘 40
𝐵𝑖 = 0.001 ≪ 1
Lumped system
𝑇 − 𝑇∞ ℎ𝐴𝑠
= exp − 𝑡
𝑇𝑖 − 𝑇∞ 𝜌𝑉𝑐
𝑇 − 𝑇∞ ℎ
= exp − 𝑡
𝑇𝑖 − 𝑇∞ 𝜌𝐿𝑐 𝑐
𝑇 − 𝑇∞ ℎ𝐴𝑠
= exp − 𝑡
𝑇𝑖 − 𝑇∞ 𝜌𝑉𝑐
𝑇 − 𝑇∞ ℎ
= exp − 𝑡
𝑇𝑖 − 𝑇∞ 𝜌𝐿𝑐 𝑐
𝑻𝒔 > 𝑻𝒔𝒖𝒓𝒓
𝑻𝒔
Vacuum
Asad Akhter Naqvi
Introduction
➢ Consider a hot object that is suspended in an evacuated chamber whose walls
are at room temperature.
➢ The hot object will eventually cool down and reach thermal equilibrium with its
surroundings.
➢ That is, it will lose heat until its temperature reaches the temperature of the
walls of the chamber.
➢ Heat transfer between the object and the chamber could not have taken place by
conduction or convection, because these two mechanisms cannot occur in a
vacuum.
➢ Therefore, heat transfer must have occurred through another mechanism that
involves the emission of the internal energy of the object. This mechanism is
radiation.
➢ We associate thermal radiation with the rate at which energy is emitted by
matter as a result of its finite temperature.
Incident radiation
Reflected
Transmitted
Asad Akhter Naqvi
Radiation Heat Fluxes
Irradiation
We define,
➢ Reflectivity (𝝆) as the fraction of the irradiation that is reflected.
➢ Absorptivity (𝜶) as the fraction of the irradiation that is absorbed.
➢ Transmissivity (𝝉) as the fraction of the irradiation that is transmitted.
𝝆+𝜶+𝝉=𝟏
A medium that experiences no transmission (𝝉 = 𝟎) is opaque, in which
case
𝝆+𝜶=𝟏
Radiosity
➢ It accounts for all the radiant energy leaving the surface.
➢ For an opaque surface, it includes emission and the reflected portion of
the irradiation.
➢ It is therefore expressed as
𝑱 = 𝑬 + 𝑮𝒓𝒆𝒇 = 𝑬 + 𝝆𝑮
𝒒′′𝒓𝒂𝒅 = 𝑱 − 𝑮 𝒒′′𝒓𝒂𝒅 = 𝑬 + 𝝆𝑮 − 𝑮
𝑟
𝑠
𝜃
𝒔 = 𝒓𝜽
𝐬
𝜽=
𝐫
𝑟
𝑟 𝑟
𝑟 𝑟
𝜃
𝑑𝜃
𝑑𝜃
𝑑∅
∅
𝒅𝑺 = 𝒓 𝐬𝐢𝐧 𝜽 𝒅∅ × 𝒓𝒅𝜽
𝒅𝑺 = 𝒓𝟐 𝐬𝐢𝐧 𝜽 𝒅∅𝒅𝜽 𝑑𝜔 =
𝑑𝑆
Asad Akhter Naqvi 𝑟2
Solid angle
𝑑𝑆
𝑑𝜔 = 2
𝑟
𝑑𝑆 = 𝑟 2 sin 𝜃 𝑑∅𝑑𝜃
𝑑𝜔 = 𝑠𝑖𝑛 𝜃 𝑑∅𝑑𝜃
2π π/2
ω = න dω = න න sin θ d∅dθ
0 0
ω = 2π steradians
=𝜋
ℎ = 6.63 × 10−34 𝐽. 𝑠
𝑘𝐵 = 1.381 × 10−23 𝐽/𝐾
𝑐𝑜 = 2.99 × 108 𝑚/𝑠
➢ Its spectral emissive power is
∞
𝑪𝟏
𝑬𝒃 = න 𝒅𝝀
𝟎 𝝀𝟓 𝐞𝐱𝐩(𝑪𝟐 Τ𝝀𝑻) − 𝟏
➢ Performing the integration, it may be shown that
𝑬𝒃 = 𝛔𝑻𝟒𝒔
➢ They may also be used to obtain the fraction of the radiation between any
two wavelengths 𝜆1 and 𝜆2
Table 12.2 is used to determine the fraction of the total black body emissions.
5 10 15 20 ∞
𝐸 = න 0𝑑𝜆 + න 100𝑑𝜆 + න 200𝑑𝜆 + න 100𝑑𝜆 + න 0𝑑𝜆
0 5 10 15 20
Part c
𝐸 = 0 + 100 10 − 5 + 200 15 − 10 + 100 20 − 15 + 0
𝐼𝑒 = 637𝑊/𝑚2 𝑠𝑟
Asad Akhter Naqvi
Lecture 16
Radiation Heat Transfer
𝑑𝐴𝑗 cos 𝜃𝑗
𝑑𝜔𝑗−𝑖 =
𝑅2
𝐼𝑒+𝑟,𝑖 cos 𝜃𝑖 cos 𝜃𝑗 𝑑𝐴𝑖 𝑑𝐴𝑗
𝑑𝑞𝑖→𝑗 =
𝑅2
➢ Assuming that surface i emits and reflects diffusely
𝑱=𝑬
𝐸𝑖 𝐽𝑖
𝐼𝑖 = = Asad Akhter Naqvi
𝜋 𝜋
The View Factor
𝐼𝑒+𝑟,𝑖 cos 𝜃𝑖 cos 𝜃𝑗 𝑑𝐴𝑖 𝑑𝐴𝑗 𝐽𝑖 cos 𝜃𝑖 cos 𝜃𝑗 𝑑𝐴𝑖 𝑑𝐴𝑗
𝑑𝑞𝑖→𝑗 = 𝑑𝑞𝑖→𝑗 =
𝑅2 𝜋𝑅 2
➢ The total rate at which radiation leaves surface i and is intercepted by
j may then be obtained by integrating over the two surfaces.
➢ From the definition of the view factor as the fraction of the radiation
that leaves 𝐴𝑖 and is intercepted by 𝐴𝑗
𝑞𝑖→𝑗
𝐹𝑖𝑗 =
𝐴𝑖 𝐽𝑖
𝑨𝒊 𝑭𝒊𝒋 = 𝑨𝒋 𝑭𝒋𝒊
➢ This expression, termed as reciprocity relation.
➢ It is useful in determining one view factor from knowledge
of the other.
𝐹𝑖𝑗 = 1
𝑗=1
4
➢ This is known as the summation rule.
3
➢ For example, applying the summation rule to surface 1 of a 5
six-surface enclosure yields
6 2
6
𝐹1𝑗 = 𝐹11 + 𝐹12 + 𝐹13 + 𝐹14 + 𝐹15 + 𝐹16 = 1
𝑗=1
𝐴1 𝐹12 = 𝐴2 𝐹21
𝐴1
𝐹21 =
𝐴2
𝐴1
𝐹22 =1−
𝐴2
𝑻𝒋 , 𝑨𝒋 , 𝜺𝒋
𝑱𝒊 𝒒𝒊 𝑮𝒊
𝒒𝒊 = 𝑨𝒊 ( 𝑱𝒊 −𝑮𝒊 )
𝑻𝒊 , 𝑨𝒊 , 𝜺𝒊
Asad Akhter Naqvi
Radiation Exchange Between Opaque, Diffuse, Gray
Surfaces in an Enclosure
➢ In general, radiation may leave an opaque surface due to both
reflection and emission, and on reaching a second opaque surface,
experience reflection as well as absorption.
➢ Analyzing radiation exchange in an enclosure may be simplified by
making certain assumptions.
➢ Each surface of the enclosure is assumed to be isothermal and to be
characterized by a uniform radiosity and a uniform irradiation.
➢ The surfaces are also assumed to be opaque 𝜏 = 0 .
➢ Surfaces have emissivity, absorptivity, and reflectivity that are
independent of direction (the surfaces are diffuse) and independent of
wavelength (the surfaces are gray).
➢ The emissivity is equal to the absorptivity 𝜀 = 𝛼 .
𝐴𝑖 𝐺𝑖 = 𝐹𝑖𝑗 𝐴𝑖 𝐽𝑗 𝐺𝑖 = 𝐹𝑖𝑗 𝐽𝑗 𝑱𝟒
𝑗=1 𝑗=1
➢ By comparing
𝑵
𝑬𝒃𝒊 − 𝑱𝒊 (𝑱𝒊 − 𝑱𝒋 )
=
(𝟏 − 𝜺𝒊 )Τ𝜺𝒊 𝑨𝒊 (𝑭𝒊𝒋 𝑨𝒋 )−𝟏
𝒋=𝟏
➢ The rate of radiation transfer (current flow) to i through its surface resistance
must equal the net rate of radiation transfer (current flows) from i to all other
surfaces through the corresponding geometrical resistances.
𝑬𝒃𝟏 − 𝑱𝟏 𝑱𝟐 − 𝑬𝒃𝟐
𝒒𝟏 = 𝒒𝟐 = 𝐴1 , 𝑇1 , 𝜀1
𝑬𝒃𝟏 (𝟏 − 𝜺𝟏 )Τ𝜺𝟏 𝑨𝟏 𝒒𝟏𝟐 (𝟏 − 𝜺𝟐 )Τ𝜺𝟐 𝑨𝟐
𝑬𝒃𝟏 − 𝑬𝒃𝟐
𝒒𝟏𝟐 =
𝑹𝟏 + 𝑹𝟏𝟐 + 𝑹𝟐
𝐴1 , 𝑇1 , 𝜀1 𝐴2 , 𝑇2 , 𝜀2
𝜀3,1 𝜀3,2
𝐴3 , 𝑇3
Asad Akhter Naqvi
Radiation Shields
➢ Radiation heat transfer between two surfaces can be reduced greatly by
inserting a thin, high-reflectivity (low-emissivity) sheet of material between
the two surfaces.
➢ Such highly reflective thin plates or shells are called radiation shields.
➢ Radiation heat transfer between two large parallel plates of emissivity
𝜀1 𝑎𝑛𝑑 𝜀2 maintained at uniform temperatures 𝑇1 and 𝑇2 .
𝐴𝜎(𝑇14 − 𝑇24 )
𝑞12 =
1 1
𝜀1 𝜀2 − 1
+
➢ To reduce radiation heat transfer, radiation shield will be added.
➢ With the radiation shield, additional resistances are present, and the heat
transfer rate is reduced.
𝐴𝜎(𝑇14 − 𝑇24 ) 𝐴2 , 𝑇2 , 𝜀2
𝑞12 = 𝐴1 , 𝑇1 , 𝜀1
1 1 1 − 𝜀3,1 1 − 𝜀3,2
𝜀1 𝜀2 + 𝜀3,1 + 𝜀3,2
+
➢ Note that the resistances associated with the radiation shield become very
large when the emissivity 𝜀3,1 and 𝜀3,2 are very small.
Asad Akhter Naqvi
Radiation Shields
𝐴𝜎(𝑇14 − 𝑇24 )
𝑞12 =
1 1 1 − 𝜀3,1 1 − 𝜀3,2
𝜀1 𝜀2 + 𝜀3,1 + 𝜀3,2
+
(𝟏 − 𝜺𝟏 )Τ𝜺𝟏 𝑨𝟏 (𝑭𝟏𝟑 𝑨𝟏 )−𝟏 (𝟏 − 𝜺𝟑,𝟏 )Τ𝜺𝟑,𝟏 𝑨𝟏 (𝟏 − 𝜺𝟑,𝟐 )Τ𝜺𝟑,𝟐 𝑨𝟏 (𝑭𝟑𝟐 𝑨𝟐 )−𝟏 (𝟏 − 𝜺𝟐 )Τ𝜺𝟐 𝑨𝟏
𝑬𝒃𝟏 𝑱𝟏 𝑱𝟑,𝟏 𝑬𝒃𝟑 𝑱𝟑,𝟐 𝑱𝟐 𝑬𝒃𝟐
𝐴2 𝐹12 = 1.0
𝐴1
𝐹21 = 𝐹
𝐴2 12
𝐴1 𝐴1 = 2𝑅 × 𝐿
3
𝐴2 = × 2𝜋𝑅 × 𝐿
4
2𝑅 × 𝐿
𝐹21 = (1.0)
3
4 × 2𝜋𝑅 × 𝐿
4
𝐹21 =
3𝜋
𝐹21 = 0.424
Look carefully
𝐴2 𝐹11 = 0
By Symmetry
𝐴1
𝐹12 = 0.5 𝐹13 = 0.5
𝑨𝟏 𝑨𝟑 𝑨𝟑
𝑨𝟐
𝑨𝟐
𝑭𝟐(𝟏,𝟑) = 𝟎. 𝟐𝟓
𝑭𝟐𝟑 = 𝟎. 𝟏𝟗
𝐴2
𝐹12 = × 𝐹21
𝐴1
0.3 × 0.5
𝐹12 = × 0.06
0.2 × 0.5
𝑭𝟏𝟐 = 𝟎. 𝟎𝟗
𝐴2
𝐴3 𝐿 = 0.3 𝑚
𝐴1
𝜋 2
𝐴1 = 𝐴2 = 𝐷 = 0.07069 𝑚2 (𝟏 − 𝜺𝒊 )Τ𝜺𝒊 𝑨𝒊
4 (𝑭𝒊𝒋+𝟏 𝑨𝒊 )−𝟏
𝐴3 = 𝜋𝐷𝐿 = 0.2827 𝑚2
𝑵
𝑬𝒃𝒊 − 𝑱𝒊 (𝑱𝒊 − 𝑱𝒋 )
=
(𝟏 − 𝜺𝒊 )Τ𝜺𝒊 𝑨𝒊 (𝑭𝒊𝒋 𝑨𝒊 )−𝟏
𝒋=𝟏
𝐿 0.3
= =2
𝑟1 0.15
𝑟2 0.15
= = 0.5
𝐿 0.3
𝐹12 = 0.17
𝐹11 = 0 (𝟏 − 𝜺𝟏 )Τ𝜺𝟏 𝑨𝟏
(𝑭𝟏𝟑 𝑨𝟏 )−𝟏
𝐹13 = 1 − 𝐹12
𝐹12 = 0.17
𝐹21 = 0.17
𝐹22 = 0
𝟑𝟓𝟒𝟒 − 𝑱𝟐 (𝑱𝟐 − 𝑱𝟏 ) (𝑱𝟐 − 𝑱𝟑 )
𝐹23 = 1 − 𝐹21 = +
(𝟏 − 𝟎. 𝟓)Τ𝟎. 𝟓 × 𝟎. 𝟎𝟕𝟎𝟔𝟗 (𝟎. 𝟏𝟕 × 𝟎. 𝟎𝟕𝟎𝟔𝟗)−𝟏 (𝟎. 𝟖𝟑 × 𝟎. 𝟎𝟕𝟎𝟔𝟗)−𝟏
𝐹23 = 1 − 0.17 = 0.83
−𝟎. 𝟏𝟕𝑱𝟏 + 𝟐𝑱𝟐 − 𝟎. 𝟖𝟐𝟖𝑱𝟑 = 𝟑𝟓𝟒𝟒 2
𝐹31 = 0.207
Similarly
𝐹32 = 0.207
𝒖∞ 𝒖
2
𝜕2𝑇 𝜕𝑢 𝜕𝑇 𝜕𝑇 𝜕𝑢 𝜕𝑣
𝑘 𝑑𝑥 𝑑𝑦 + 𝜇 𝑑𝑥𝑑𝑦 = 𝜌𝐶𝑝 𝑢 +𝑣 + 𝑇 + 𝑑𝑥𝑑𝑦
𝜕𝑦 2 𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑦
2
𝜕2𝑇 𝜕𝑢 𝜕𝑇 𝜕𝑇 𝜕𝑢 𝜕𝑣
𝑘 +𝜇 = 𝜌𝐶𝑝 𝑢 +𝑣 + =0
𝜕𝑦 2 𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑦
𝜕𝑇 𝜕𝑇 𝜕2 𝑇 𝜇 𝜕𝑢 2
𝑢 𝜕𝑥 + 𝑣 𝜕𝑦 = 𝛼 𝜕𝑦 2 + 𝜌𝐶𝑝 𝜕𝑦
For low velocity flows the viscous work term may be neglected
𝜕𝑇 𝜕𝑇 𝜕2𝑇
𝑢 +𝑣 =𝛼 2
𝜕𝑥 𝜕𝑦 𝜕𝑦
The differential form of the boundary layer energy equation
Asad Akhter Naqvi
Local and Average Convection Coefficients
➢ A fluid of velocity 𝑉 and temperature 𝑇∞ flows over a surface of
arbitrary shape and of area 𝐴𝑠 .
➢ The surface is presumed to be at a uniform temperature, 𝑇𝑠 ..
➢ We know that convection heat transfer will occur if 𝑇𝑠 ≠ 𝑇∞ .
➢ the surface heat flux and convection heat transfer coefficient both
vary along the surface.
➢ The total heat transfer rate q may be obtained by integrating the
local flux over the entire surface.
➢ where n is a positive exponent. Hence for a gas ; for a liquid metal ; for an oil .
Asad Akhter Naqvi
Dimensionless Parameters
Nusselt Number
➢ Ratio of convection to pure conduction heat transfer.
𝑞𝑐𝑜𝑛𝑣
𝑁𝑢 =
𝑞𝑐𝑜𝑛𝑑
➢ This parameter is equal to the dimensionless temperature gradient at the surface, and it
provides a measure of the convection heat transfer occurring at the surface.
➢ The Nusselt number is to the thermal boundary layer what the friction coefficient is to the
velocity boundary layer.
𝐻 𝛿
𝛿𝑡
a d
𝑑𝑥
𝐻
1. Rate at which energy advected into the left face = 𝜌𝐶𝑝 0 𝑢𝑇𝑑𝑦
𝐻 𝑑 𝐻
2. Rate at which energy advected out of the right face = 𝜌𝐶𝑝 0 𝑢𝑇𝑑𝑦 + 𝑑𝑥 0 𝑢𝑇𝑑𝑦 𝑑𝑥
𝑑 𝐻
ሶ 𝐶𝑝 𝑇∞ = (𝑚ሶ 𝑐𝑑 − 𝑚𝑎𝑏
3. Rate at which energy advected from the top = 𝑚𝑏𝑐 ሶ )𝐶𝑝 𝑇∞ = 𝜌𝐶𝑝 𝑇∞ 0 𝑢 𝑑𝑦 𝑑𝑥
𝑑𝑥
𝜕𝑇
4. Rate at which energy is conducted into the bottom face= −𝑘 𝑑𝑥 𝜕𝑦ቚ 𝑦 = 0)
Asad Akhter Naqvi
Integral Analysis of the boundary layer
Applying Energy balance
1 + 3 + 4 = (2)
𝐻
𝑑 𝜕𝑇
න 𝑇 − 𝑇 𝑢𝑑𝑦 = 𝛼 ቤ (𝑦 = 0) 1
𝑑𝑥 0 ∞ 𝜕𝑦
𝑘
𝛼=
𝜌𝐶𝑝
We have learned
3
𝑢 3 𝑦 1 𝑦 3 𝑇 − 𝑇𝑠 3 𝑦 1 𝑦
= − 2 = − 3
𝑢∞ 2 𝛿 2 𝛿 𝑇∞ − 𝑇𝑠 2 𝛿𝑡 2 𝛿𝑡
𝜕𝑇 3 𝑇∞ − 𝑇𝑠
ቤ (𝑦 = 0) = 4
𝜕𝑦 2 𝛿𝑡
• Substitute eqs. 2, 3, and 4 in eq. 1 • Define 𝜁 = 𝛿𝑡ൗ𝛿 • Integrate from 0 𝑡𝑜 𝛿𝑡
𝑑 3 2 3 4 3𝛼(𝑇∞ − 𝑇𝑠 )
𝑇∞ − 𝑇𝑠 𝑢∞ 𝛿 𝜁 − 𝜁 =
𝑑𝑥 20 280 2𝛿𝑡 Asad Akhter Naqvi
Integral Analysis of the boundary layer
𝑑 3 2 3 4 3𝛼(𝑇∞ − 𝑇𝑠 )
𝑇∞ − 𝑇𝑠 𝑢∞ 𝛿 𝜁 − 𝜁 =
𝑑𝑥 20 280 2𝛿𝑡
Since 𝜁 < 1 , it follows that 𝜁 4 < 𝜁 2 and the corresponding term may be neglected. Thus
3 𝑑 2 3𝛼 𝑢∞ 𝑑 2 𝛼
𝑢∞ 𝛿𝜁 = 𝛿𝜁 =
20 𝑑𝑥 2𝛿𝑡 10 𝑑𝑥 𝛿𝑡
𝑢∞
2𝛿𝜁
𝑑𝜁
+ 𝜁 2 𝑑𝛿
=
𝛼 𝑢∞ 2 2
𝑑𝜁 3
𝑑𝛿
10 𝑑𝑥 𝑑𝑥 δζ 2𝛿 𝜁 + δ𝜁 =𝛼 5
10 𝑑𝑥 𝑑𝑥
Also
1/2 1/2 280 𝜈
𝛿 4.64 280 𝜈 𝛿2 = 6
= 1/2 = 13 𝑢∞ 𝑥
𝑥 𝑅𝑒 13 𝑢∞ 𝑥
𝑥
Differentiating eq. 6
𝑑𝛿 140 𝜈
𝛿 = 7
𝑑𝑥 13 𝑢∞
Substituting eq. 6 and 7 in eq. 5 Asad Akhter Naqvi
Integral Analysis of the boundary layer
𝑑𝜁 13 𝛼 3
4 𝑑(𝜁 3 ) 13 𝛼
𝜁 3 + 4𝑥𝜁 2 = 𝜁 + 𝑥 = 8
𝑑𝑥 14 𝜈 3 𝑑𝑥 14 𝜈
4
Dividing eq. 8 by
3
3 −1/4 3 𝑑(𝜁 3 ) 39 𝛼 −1/4 𝑑(𝑥 3/4 𝜁 3 ) 39 𝛼 −1/4 9
𝑥 𝜁 +𝑥 3/4
= 𝑥 = 𝑥
4 𝑑𝑥 56 𝜈 𝑑𝑥 56 𝜈
Integrating eq. 9
13 𝛼 3/4
𝑥 3/4 𝜁 3 = 𝑥 +𝐶 𝑢∞ 𝑇∞
14 𝜈
13 1 𝑇𝑠 = 𝑇∞ 𝑇𝑠 > 𝑇∞
𝜁3 = + 𝐶𝑥 −3/4 10
14 𝑃𝑟
At 𝑥𝑜 𝑥=𝐿
𝑥 = 𝑥𝑜 𝜁=0
𝑥
13 1 3/4
𝐶=− 𝑥
14 𝑃𝑟 𝑜 Asad Akhter Naqvi
Integral Analysis of the boundary layer
Eq. 10 will be
3/4 3/4 1/3
13 1 𝑥 1 1 𝑥
𝜁= 1− 11
𝜁3 = 1− 1.025 𝑃𝑟 −1/3 𝑥𝑜
14 𝑃𝑟 𝑥𝑜
We have learned
𝜕𝑇
−𝑘 ฬ (𝑦 = 0) 𝜕𝑇 3 𝑇∞ − 𝑇𝑠
𝜕𝑦 ቤ (𝑦 = 0) =
ℎ= 𝜕𝑦 2 𝛿𝑡
(𝑇𝑠 − 𝑇∞ )
3𝑘 3 𝑘
ℎ= ℎ=
2 𝛿𝑡 2 𝛿𝜁
3/4 −1/3
𝑘 1/2 𝑥
ℎ = 0.332 𝑅𝑒𝑥 𝑃𝑟1/3 1 −
𝑥 𝑥𝑜
ℎ𝑥
𝑁𝑢𝑥 =
𝑘
3/4 −1/3
1/2 𝑥
𝑁𝑢𝑥 = 0.332𝑅𝑒𝑥 𝑃𝑟1/3 1−
𝑥𝑜
This equation is used to determine the Nusselt number for laminar flow with constant surface temperature
The friction coefficient will be given by The convection heat transfer coefficient will be given by
𝜏𝑠 𝜕𝑇
𝐶𝑓 = 2 −𝑘 ฬ (𝑦 = 0)
𝜌𝑢∞ /2 𝜕𝑦
ℎ=
(𝑇𝑠 − 𝑇∞ )
𝜕𝑇
= 𝐸 + 2𝐹𝑦 − 3𝐺𝑦 2
𝜕𝑢 𝜕𝑦
= 𝐴 + 2𝐵𝑦 − 3𝐶𝑦 2
𝜕𝑦
At 𝑦 = 0
At 𝑦 = 0
𝜕𝑢 𝜕𝑇
(𝑦 = 0) = 𝐴 =𝐸
𝜕𝑦 𝜕𝑦
𝜏𝑠 = 𝐴𝜇 −𝑘𝐸
ℎ=
(𝐷 − 𝑇∞ )
2𝐴𝜇
𝐶𝑓 = 2 Asad Akhter Naqvi
𝜌𝑢∞
6-2. Parallel flow of atmospheric air over a flat plate of length 𝐿 = 3 𝑚 is disrupted by an array of stationary rods placed in the
flow path over the plate.
Laboratory measurements of the local convection coefficient at the surface of the plate are made for a prescribed value of 𝑉 and
𝑇𝑠 > 𝑇. The results are correlated by an expression of the form ℎ𝑥 = 0.7 + 13.6𝑥 − 3.4𝑥 2 , where ℎ𝑥 has units of 𝑊/𝑚2 𝐾 and 𝑥
is in meters. Evaluate the average convection coefficient ℎ𝐿 ℎ𝐿 Τℎ𝐿 for the entire plate and the ratio ℎ𝐿 Τℎ𝐿 at the trailing edge.
For 𝐿 = 3𝑚
ℎ𝑥 = 0.7 + 13.6(3) − 3.4(3)2
𝒉𝒙 = 𝟏𝟎. 𝟗 𝑾/𝒎𝟐 𝑲
1 𝐿
ℎത 𝐿 = න ℎ𝑥 𝑑𝑥
𝐿 0
1 3
ℎത 𝐿 = න 0.7 + 13.6𝑥 − 3.4𝑥 2 𝑑𝑥
3 0
ഥ 𝑳 = 𝟏𝟎. 𝟗 𝑾/𝒎𝟐 𝑲
𝒉
ഥ 𝑳 𝟏𝟎. 𝟗
𝒉
= =𝟏
𝒉𝒙 𝟏𝟎. 𝟗
𝑇∞ = 300 𝐾 𝑇∞ = 300 𝐾
𝑇𝑠 = 400𝐾 𝑇𝑠 = 400𝐾
𝐿𝑐 = 1 𝑚 𝐿𝑐 = 5 𝑚
Assume 𝑞
𝑇𝑠 = 𝑇∞ +
VL ℎ𝐴
𝑇𝑠 = 45𝑜 𝐶 𝑅𝑒 =
𝜈 30 × 10−3
𝑇𝑠 + 𝑇∞ 𝑇𝑠 = 25 +
𝑇ത = = 35𝑜 𝐶 10 × 0.12 107 × 0.004 × 0.004
2 𝑅𝑒 =
16.69 × 10−6
At 𝑇𝑠 = 42.5𝑜 𝐶
𝑅𝑒 = 71.9 × 103
𝑇ത = 35𝑜 𝐶
ℎ𝑥 𝑥
𝑁𝑢𝑥 = = 0.332𝑅𝑒 1/2 𝑃𝑟1/3
𝑘
1 𝑥
ℎ𝑥 = න ℎ𝑥 𝑑𝑥
𝑥 0
1/2 𝑥
𝑘 𝑢∞ 1
ℎ𝑥 = 0.332 𝑃𝑟1/3 න 𝑑𝑥
𝑥 𝜈 0 𝑥 1/2
ℎ𝑥 = 2ℎ𝑥
𝛿 = 0.37𝑥 𝑅𝑒 −1/5
➢ Comparing these results with those for the laminar boundary layer,
we see that turbulent boundary layer growth is much more rapid.
➢ The local Nusselt number for turbulent flow is
4/5
𝑁𝑢𝑥 = 0.0296𝑅𝑒𝑥 𝑃𝑟1/3
➢ Enhanced mixing causes the turbulent boundary layer to grow more
rapidly than the laminar boundary layer and to have larger friction
and convection coefficients.
𝑇∞ = 29𝑜 𝐶
𝑢∞ = 5 𝑚/𝑠 𝑇𝑠 = 75𝑜 𝐶
0.5 𝑚
1𝑚
At 𝑇𝑓 = 325 𝐾
𝑁𝑢𝐿 = 308
Asad Akhter Naqvi
𝑇∞ + 𝑇𝑠
𝑇𝑓 = = 52𝑜 𝐶 = 325 𝐾 𝑇∞ = 29𝑜 𝐶
2 𝑇𝑠 = 75𝑜 𝐶
𝑢∞ = 5 𝑚/𝑠
At 𝑇𝑓 = 325 𝐾
At trailing edge
1/2
𝑁𝑢𝑥 = 0.332𝑅𝑒𝑥 𝑃𝑟1/3
4.64𝑥
𝛿=
𝑅𝑒𝑥
4.64(1)
𝛿=
2.71 × 105
𝛿 = 8.91 × 10−3 𝑚
Asad Akhter Naqvi
Air at 15𝑜 𝐶flows at 1.8 𝑚/𝑠 over a 0.6 𝑚 long heating panel. The panel is intended to supply
420 𝑊/𝑚2 to the air but the surface can sustain only about 110𝑜 𝐶 without being damaged. Is it
safe? What is the average temperature of the plate?
𝑢∞ = 1.8 𝑚/𝑠
𝑇∞ = 15𝑜 𝐶
0.6 𝑚
𝑢∞ = 60 𝑚/𝑠
𝑇∞ = 25𝑜 𝐶
𝑇𝑠 = 230𝑜 𝐶
50 𝑚𝑚
𝑞 = ℎ𝐴(𝑇𝑠 − 𝑇∞ ) ℎ𝐿
𝑁𝑢𝑥 =
𝑘
𝑁𝑢𝑥 = 1671
Asad Akhter Naqvi
7-43. The cover plate of a flat-plate solar collector is at 15𝑜 𝐶, while ambient air at 10𝑜 𝐶 is in
parallel flow over the plate, with 𝑢∞ = 2 𝑚/𝑠.
Flow is laminar
𝑞 = 39 W
Asad Akhter Naqvi
Lecture 22
External Flow over the cylinders.
➢ As the fluid decelerates, the velocity gradient at the surface, 𝜕𝑢 Τ𝜕𝑦ȁ(𝑦 = 0),
eventually becomes zero.
➢ At this location, termed the separation point, fluid near the surface lacks
sufficient momentum to overcome the pressure gradient, and continued
downstream movement is impossible.
➢ This is a condition for which the boundary layer detaches from the surface,
and a wake is formed in the downstream region.
➢ Flow in this region is characterized by vortex formation and is highly
irregular.
➢ The separation point is the location for which 𝜕𝑢 Τ𝜕𝑦ȁ(𝑦 = 0).
➢ If 𝑅𝑒𝐷 ≤ 2 × 105 , the boundary layer remains laminar, and separation occurs
at θ = 80𝑜
➢ However, if 𝑅𝑒𝐷 ≥ 2 × 105 , boundary layer transition occurs, and separation
is delayed to θ = 140𝑜
𝑉 = 15 𝑚/𝑠
𝑇∞ = 25𝑜 𝐶
𝑇𝑠 = 100𝑜 𝐶
𝑇𝑠 +𝑇∞
𝑇𝑓 = = 335 𝐾 ℎത = 88.2 𝑊/𝑚2 𝐾
2
At 335 𝐾 ത
𝑞 = ℎ𝐴(𝑇𝑠 − 𝑇∞ )
𝑘 = 0.0288 𝑊/𝑚𝐾
𝜈 = 19.31 × 10−6 𝑚2 /𝑠 𝑞 = ℎത × (𝜋𝐷𝐿) × (𝑇𝑠 − 𝑇∞ )
𝑃𝑟 = 0.7
𝐶 = 0.193 𝑚 = 0.618 𝑞
𝑉𝐷 𝑞′ = = ℎത × (𝜋𝐷) × (𝑇𝑠 − 𝑇∞ )
𝑅𝑒 = = 19420 𝐿
𝜈 𝑁𝑢𝐷 = 0.193 × 194200.618 × 0.71Τ3
ത
ℎ𝐷 ത
ℎ𝐷
𝑁𝑢𝐷 = = 𝐶𝑅𝑒𝐷𝑚 𝑃𝑟 1Τ3 𝑁𝑢𝐷 =
𝑘 𝑘 𝑞′ = 520 𝑊/𝑚
𝑁𝑢𝐷 𝑘
ℎത =
𝐷
Asad Akhter Naqvi
External Flow over a sphere
➢ Boundary layer effects associated with flow over a sphere are much like those
for the circular cylinder, with transition and separation playing prominent
roles.
➢ Numerous heat transfer correlations have been proposed, and Whitaker [15]
recommends an expression of the form
2
𝑢 𝑟
=2 1−
𝑢𝑚 𝑟𝑜
The Entrance Region
➢ When dealing with internal flows, it is important to be cognizant of the extent
of the entry region, which depends on whether the flow is laminar or
turbulent.
➢ The Reynolds number for flow in a circular tube is defined as
𝜌𝑢𝑚 𝐷 𝑢𝑚 𝐷
𝑅𝑒𝐷 = =
𝜇 𝜈
➢ In a fully developed flow, the critical Reynolds number corresponding to the
onset of turbulence is
𝑅𝑒𝐷,𝑐 ≈ 2300
➢ For laminar flow the hydrodynamic entry length may be obtained from an
expression of the form
𝑥𝑓𝑑,ℎ
≈ 0.05𝑅𝑒𝐷
𝐷 𝑙𝑎𝑚
➢ This expression is based on the presumption that fluid enters the tube from a
rounded converging nozzle and is hence characterized by a nearly uniform
velocity profile at the entrance.
➢ There is no satisfactory general expression for the entry length in turbulent
flow
➢ It is approximately independent of Reynolds number and that,
Thermal Considerations
➢ If fluid enters the tube at a uniform temperature T(r, 0) that is
less than the surface temperature
➢ Convection heat transfer occurs and a thermal boundary
layer begins to develop.
➢ The shape of the fully developed temperature profile
𝑇(𝑟, 𝑥) differs according to whether a uniform surface
temperature or heat flux is maintained.
➢ For laminar flow the thermal entry length may be expressed as
𝑥𝑓𝑑,𝑡
≈ 0.05𝑅𝑒𝐷 𝑃𝑟
𝐷 𝑙𝑎𝑚
➢ Comparing Equations of hydrodynamic and thermal
conditions it is evident that, if Pr > 1, the hydrodynamic
boundary layer develops more rapidly than the thermal
boundary layer , while the inverse is true for Pr < 1.
𝑑𝑇𝑚 𝑞 ′′ 𝑠 𝑃 𝑃
= = ℎ(𝑇𝑠 − 𝑇𝑚 )
𝑑𝑥 𝑚𝐶𝑝 𝑚𝐶𝑝
The Energy Balance
Constant Surface Heat Flux
𝑞𝑐𝑜𝑛𝑣 = 𝑞 ′′ 𝑠 ( 𝑃. 𝐿)
➢ For constant 𝑞 ′′ 𝑠
𝑑𝑇𝑚 𝑞 ′′ 𝑠 𝑃
= ≠ 𝑓(𝑥)
𝑑𝑥 𝑚𝐶𝑝
𝑞′′ 𝑠 𝑃
𝑇𝑚 𝑥 = 𝑇𝑚,𝑖 + 𝑥
𝑚𝐶𝑝
Accordingly, the mean temperature varies linearly with x along the tube
we also expect the temperature difference (𝑇𝑠 − 𝑇𝑚 ) to vary with x,
The Energy Balance
Constant Surface Temperature
➢ Results for the total heat transfer rate and the axial distribution of the mean
temperature are entirely different for the constant surface temperature
condition.
➢ Defining ∆𝑇 = 𝑇𝑠 − 𝑇𝑚 ,
𝑑𝑇𝑚 𝑑(∆𝑇) 𝑃
=− = ℎ∆𝑇
𝑑𝑥 𝑑𝑥 𝑚𝐶𝑝
➢ Separating variables and integrating from the tube inlet to the outlet
∆𝑇𝑜 𝐿
𝑑(∆𝑇) 𝑃 ∆𝑇𝑜 𝑃𝐿 1 𝐿
න =− න ℎ𝑑𝑥 ln =− න ℎ𝑑𝑥
∆𝑇𝑖 ∆𝑇 𝑚𝐶𝑝 0 ∆𝑇𝑖 𝑚𝐶𝑝 𝐿 0
∆𝑇𝑜 𝑃𝐿 ℎ𝐿 𝐴𝑆
ln =− ℎ 𝑚𝐶𝑝 = −
∆𝑇𝑖 𝑚𝐶𝑝 𝐿 ∆𝑇
ln ∆𝑇𝑜
𝑖
𝑞𝑐𝑜𝑛𝑣 = ℎ𝐿 𝐴𝑆 ∆𝑇𝑙𝑚
∆𝑇𝑜 − ∆𝑇𝑖
∆𝑇𝑙𝑚 =
∆𝑇
ln ∆𝑇𝑜
𝑖
Lecture 24
Internal Flows (Cont.)
Temperature Profile
Assumptions
Energy Balance
𝑑𝑥
4 − 3 = 1 − (2) 𝑥 𝑟
𝜕𝑇 𝜕𝑇 𝜕𝑇 𝜕2 𝑇
𝜌𝐶𝑝 𝑢 𝑇 + 𝑑𝑥 2𝜋𝑟𝑑𝑟 − 𝜌𝐶𝑝 𝑢𝑇 2𝜋𝑟𝑑𝑟 = −𝑘 2𝜋𝑟𝑑𝑥 + 𝑘 2𝜋 𝑟 + 𝑑𝑟 𝑑𝑥 + 𝑑𝑟
𝜕𝑥 𝜕𝑟 𝜕𝑟 𝜕𝑟 2
𝜕𝑇 𝜕2𝑇 𝜕𝑇 𝜕2𝑇 2
𝜌𝐶𝑝 𝑢 2𝜋𝑟𝑑𝑟𝑑𝑥 = 2𝜋𝑘𝑑𝑥 𝑟𝑑𝑟 + dr + 2 𝑑𝑟
𝜕𝑥 𝜕𝑟 2 𝜕𝑟 𝜕𝑟
2
Since 𝑑𝑟 = 0, the above expression will be
Temperature Profile
𝜕𝑇 𝜕2𝑇 𝜕𝑇
𝜌𝐶𝑝 𝑢 2𝜋𝑟𝑑𝑟𝑑𝑥 = 2𝜋𝑘𝑑𝑥 𝑟𝑑𝑟 + dr
𝜕𝑥 𝜕𝑟 2 𝜕𝑟
𝜌𝐶𝑝 𝜕𝑇 𝜕2𝑇 𝜕𝑇 𝜌𝐶𝑝 𝜕𝑇 𝜕 𝜕𝑇
𝑢𝑟 = 𝑟 + 𝑢𝑟 = 𝑟
𝑘 𝜕𝑥 𝜕𝑟 2 𝜕𝑟 𝑘 𝜕𝑥 𝜕𝑟 𝜕𝑟
𝑢 𝑟
2 𝑢𝑐
𝑢𝑟 𝜕𝑇 𝜕 𝜕𝑇 =2 1− 𝑢𝑚 =
= 𝑟 𝑢𝑚 𝑟𝑜 2
𝛼 𝜕𝑥 𝜕𝑟 𝜕𝑟
𝜕 𝜕𝑇 𝑟 2 𝑟 𝜕𝑇 (1)
𝑟 𝜕𝑟 = 𝑢𝑐 1 −
𝜕𝑟 𝑟𝑜 𝛼 𝜕𝑥
Integrating eq. 1
𝜕𝑇 𝑢𝑐 𝑟 2 𝑟 4 𝜕𝑇 (2)
𝑟 = − + 𝐶1
𝜕𝑟 𝛼 2 4𝑟𝑜2 𝜕𝑥
Integrating eq. 2
𝑢𝑐 𝑟 2 𝑟 4 𝜕𝑇
𝑇= − + 𝐶1 𝑙𝑛𝑟 + 𝐶2
𝛼 4 16𝑟𝑜2 𝜕𝑥
Temperature Profile
𝑢𝑐 𝑟 2 𝑟 4 𝜕𝑇
𝑇= − + 𝐶1 𝑙𝑛𝑟 + 𝐶2
𝛼 4 16𝑟𝑜2 𝜕𝑥
Boundary Conditions
𝜕𝑇
= 0 𝑎𝑡 𝑟 = 0
𝜕𝑟
𝑇 = 𝑇𝑐 𝑎𝑡 𝑟 = 0
𝐶1 = 0 𝐶2 = 𝑇𝑐
𝑢𝑐 𝑟 2 𝑟 4 𝜕𝑇
𝑇 = 𝑇𝑐 + −
𝛼 4 16𝑟𝑜2 𝜕𝑥
The Mean Temperature
𝑟𝑜
𝑞 = 𝑚𝐶𝑝 𝑇𝑚 = න 𝐶𝑝 𝑇𝑑𝑚
0
𝑟𝑜
𝑚𝐶𝑝 𝑇𝑚 = න 𝜌𝑢𝐶𝑝 𝑇𝑑𝐴𝑐 (3)
0
For circular pipes
𝐴𝑐 = 𝜋𝑟 2 𝑑𝐴𝑐 = 2𝜋𝑟𝑑𝑟
𝜌𝑢𝑐
𝑚 = 𝜌𝑢𝑚 𝜋𝑟𝑜2 = (𝜋𝑟𝑜2 )
2
𝑢𝑐 𝑟 2 𝑟 4 𝜕𝑇
𝑇 = 𝑇𝑐 + −
𝛼 4 16𝑟𝑜2 𝜕𝑥
2
𝑢 𝑟
= 1−
𝑢𝑐 𝑟𝑜
At 𝑟 = 𝑟𝑜 𝑇𝑚.𝑖 + 𝑇𝑚,𝑜
𝑇𝑚 =
2
𝜕𝑇 𝑢𝑐 𝑟𝑜 𝜕𝑇
𝑟 = 𝑟𝑜 = Constant Surface temperature Condition
𝜕𝑟 4𝛼 𝜕𝑥
24 𝑘 48 𝑘
ℎ= =
11 𝑟𝑜 11 𝐷
𝑁𝑢𝐷 = 4.364
The Entry Region
➢ The results of the preceding section are valid only when both the
velocity and temperature profiles are fully developed,
➢ However, two different entry length are also associated
➢ The simplest solution is for the thermal entry length problem, and it is
based on assuming that thermal conditions develop in the presence of a
fully developed velocity profile.
➢ the combined (thermal and velocity) entry length problem corresponds
to the case for which the temperature and velocity profiles develop
simultaneously.
➢ It would never be the case that thermal conditions are fully developed
and hydrodynamic conditions are developing.
➢ Since the temperature distribution depends on the velocity distribution,
as long as the velocity is still changing, thermal conditions cannot be
fully developed.
The Entry Region
For thermal entry length problem
0.0668(𝐷 Τ𝐿)𝑅𝑒𝐷 𝑃𝑟
𝑁𝑢𝐷 = 3.66 +
1 + 0.04 (𝐷 Τ𝐿)𝑅𝑒𝐷 𝑃𝑟 2/3
1Τ3 0.14
𝑅𝑒𝐷 𝑃𝑟 𝜇
𝑁𝑢𝐷 = 1.86
𝐿Τ𝐷 𝜇𝑠
Turbulent Flow in Circular Tubes
➢ For fully developed (hydrodynamically and thermally) turbulent flow in a
smooth circular tube, the local Nusselt number may be obtained from the
Dittus. Boelter equation
Τ
𝑁𝑢𝐷 = 0.023𝑅𝑒𝐷4 5 𝑃𝑟 𝑛
➢ where 𝑛 = 0.4 for heating (𝑇𝑠 > 𝑇𝑚 ) and 0.3 for cooling (𝑇𝑠 < 𝑇𝑚 ).
➢ These equations have been confirmed experimentally for the range of
conditions
➢ Since entry lengths for turbulent flow are typically short, 10 ≤ (𝑥𝑓𝑑 /𝐷) ≤ 60,
it is often reasonable to assume that the average Nusselt number for the entire
tube is equal to the value associated with the fully developed region, 𝑁𝑢𝐷 ≈
𝑁𝑢𝐷,𝑓𝑑 .
Turbulent Flow in Circular Tubes
➢ For flows characterized by large property variations, the following equation, due
to Sieder and Tate
0.14
Τ 𝜇
𝑁𝑢𝐷 = 0.027𝑅𝑒𝐷4 5 𝑃𝑟1/3
𝜇𝑠
𝑄ሶ = 3.25 × 10−6 𝑚3 /𝑠
𝑇𝑏𝑜
𝑇𝑏𝑖 = 20𝑜 𝐶
𝑇𝑠
𝐷 = 1.4 𝑐𝑚
𝐿 =3𝑚
Taking Properties at 𝑇𝑏𝑖 = 20𝑜 𝐶 𝑞 = 𝑚𝐶
ሶ 𝑝 (𝑇𝑏𝑜 − 𝑇𝑏𝑖 )
𝑞
𝑇𝑏𝑜 = + 𝑇𝑏𝑖
From table A 5 𝑚𝐶
ሶ 𝑝
𝐶𝑃 = 2382 𝐽/𝑘𝑔𝐾 𝑞
𝑇𝑏𝑜 = + 𝑇𝑏𝑖
𝜌 = 1116 𝑘𝑔/𝑚3 𝜌𝑄𝐶𝑝
ℎ=
Nuk 𝑇𝑠𝑜 = 81.9𝑜 𝐶
D
Checking either flow is developed or not
4.364 × 0.249
𝑥ℎ = 0.05𝑅𝑒𝐷 𝐷 ℎ=
0.014
𝑥ℎ = 0.011 𝑚
ℎ = 77.62 𝑊/𝑚2 𝐾
𝑥𝑡 = 0.05𝑅𝑒𝐷 𝐷𝑃𝑟
𝑥𝑡 = 2.303 𝑚
𝑞 ′′ = ℎ(𝑇𝑠𝑜 − 𝑇𝑏𝑜 )
𝑚ሶ = 0.5 𝑘𝑔/𝑠
𝑇𝑏𝑜
𝑇𝑏𝑖 = 20𝑜 𝐶
𝑇𝑠
𝐷 = 50 𝑚𝑚
𝐿 = 25 𝑚
Taking Properties at 𝑇𝑏𝑖 = 20𝑜 𝐶 Checking flow nature
ℎ𝐷 𝑞 = 15564 𝑊
𝑁𝑢 =
𝑘
𝑁𝑢𝑘
ℎത = = 32.61 𝑊/𝑚2 𝐾
𝐷
𝑇𝑠 − 𝑇𝑏𝑖 𝜋𝐷𝐿
ln = ℎ
𝑇𝑠 − 𝑇𝑏𝑜 𝑚𝐶𝑝 𝐿
Engine oil flows through a 25 𝑚𝑚 diameter tube at a rate of 0.5 𝑘𝑔/𝑠. The oil enters the tube at a
temperature of 25𝑜 𝐶, while the tube surface temperature is maintained at 100𝑜 𝐶
Determine the oil outlet temperature for a 5 − 𝑚 and for a 100 − 𝑚 long tube. For each case,
compare the log mean temperature difference to the arithmetic mean temperature difference
𝑚ሶ = 0.5 𝑘𝑔/𝑠
𝑇𝑏𝑜
𝑇𝑏𝑖 = 20𝑜 𝐶
𝑇𝑠
𝐷 = 25 𝑚𝑚
𝐿 =5𝑚
Assume 𝑚ሶ = 𝜌𝐴𝑢𝑚 Checking either flow is developed or not
𝑇𝑚𝑜 = 85𝑜 𝐶 𝑚ሶ 𝑥ℎ = 0.05𝑅𝑒𝐷 𝐷
𝑢𝑚 =
𝜌𝐴 𝑥ℎ = 0.380 𝑚
𝑇𝑚𝑖 + 𝑇𝑚𝑜
ത
𝑇= 0.5
2 𝑢𝑚 = 𝑥𝑡 = 0.05𝑅𝑒𝐷 𝐷𝑃𝑟
𝜋
865.8 × 4 (0.025)2 𝑥𝑡 = 458.3 𝑚
𝑇ത = 328𝐾 ≈ 330 𝐾
𝑢𝑚 = 1.176 𝑚/𝑠
From table A 5 𝑥ℎ is less than 5 m but 𝑥𝑡 is greater than 5 m
𝐶𝑃 = 2035 𝐽/𝑘𝑔𝐾 Checking flow nature It means flow is thermally developing.
𝜌 = 865.8 𝑘𝑔/𝑚3
𝑘 = 0.141 𝑊/𝑚𝐾 𝑢𝑚 𝐷 1.176 × 0.025
𝑅𝑒𝐷 = =
𝜈 = 96.6 × 10−6 𝑚2 /𝑠 𝜈 96.6 × 10−6
𝑃𝑟 = 1205
𝑅𝑒𝐷 = 304.3 < 2300
Flow is laminar
0.0668(𝐷 Τ𝐿)𝑅𝑒𝐷 𝑃𝑟 ∆𝑇𝑜 − ∆𝑇𝑖
𝑁𝑢𝐷 = 3.66 + ∆𝑇𝑙𝑚 =
1 + 0.04 (𝐷 Τ𝐿)𝑅𝑒𝐷 𝑃𝑟 2/3 ∆𝑇
ln ∆𝑇𝑜
𝑖
0.0668(0.025Τ5) × 304.3 × 1250
𝑁𝑢𝐷 = 3.66 +
1 + 0.04 (0.025Τ5) × 304.3 × 1250 2/3
𝑁𝑢𝐷 = 21.16
(100 − 28.4) − (100 − 25)
ℎ𝐷 ∆𝑇𝑙𝑚 =
𝑁𝑢 = (100 − 28.4)
𝑘 ln
(100 − 25)
𝑁𝑢𝑘
ത
ℎ= = 119 𝑊/𝑚2 𝐾 ∆𝑇𝑙𝑚 = 73.3𝑜 𝐶
𝐷
The outlet temperature will be determined by ∆𝑇𝑜 + ∆𝑇𝑖
∆𝑇𝑎𝑚 =
2
𝑇𝑠 − 𝑇𝑏𝑖 𝜋𝐷𝐿
ln = ℎ
𝑇𝑠 − 𝑇𝑏𝑜 𝑚𝐶𝑝 𝐿 ∆𝑇𝑎𝑚 = 73.3𝑜 𝐶
𝑇𝑏𝑜 = 28.4𝑜 𝐶
Water is to be heated from 15°C to 65°C as it flows through a 3-cm-internal diameter 5-m-long
tube (Fig. 8–29). The tube is equipped with an electric resistance heater that provides uniform
heating throughout the surface of the tube. The outer surface of the heater is well insulated, so that
in steady operation all the heat generated in the heater is transferred to the water in the tube. If the
system is to provide hot water at a rate of 10 L/min, determine the power rating of the resistance
heater. Also, estimate the inner surface temperature of the pipe at the exit.
Average Film temperature 𝑞 = 𝑚𝐶
ሶ 𝑝 (𝑇𝑏𝑜 − 𝑇𝑏𝑖 )
Τ
𝑁𝑢𝐷 = 0.023𝑅𝑒𝐷4 5 𝑃𝑟 0.4
𝑁𝑢𝐷 = 0.023(10750)0.8 4.320.4
𝑁𝑢𝐷 = 69.3
𝑁𝑢𝑘
ℎത = = 1457 𝑊/𝑚2 𝐾
𝐷
𝑞
𝑇𝑠𝑜 = + 𝑇𝑏𝑜
ℎ(𝜋𝐷𝐿)
𝑇𝑠𝑜 = 115𝑜 𝐶
Lecture 26
Heat Transfer due to free convection
Introduction
• In free convection fluid motion is due to buoyancy forces within the fluid, while in forced
convection it is externally imposed.
• Buoyancy is due to the combined presence of a fluid density gradient and a body force that is
proportional to density.
• In practice, the body force is usually gravitational, although it may be a centrifugal force in
rotating fluid machinery or a Coriolis force in atmospheric and oceanic rotational motions.
• The fluid will move due to density gradient.
• The density gradient is due to a temperature gradient and the body force is gravitational.
Thermal Stratification
𝑇2 𝑇1
𝑇1 𝑇2
𝑇1 > 𝑇2 𝑇1 < 𝑇2
Heat transfer will be due to convection. Heat transfer will be due to conduction only
Buoyancy Forces
➢ In a gravitational field, there is a net force that pushes upward a light fluid placed in a heavier fluid.
➢ The upward force exerted by a fluid on a body completely or partially immersed in it is called the buoyancy force.
➢ In heat transfer studies, the primary variable is temperature, and it is desirable to express the net buoyancy force in terms of
temperature differences.
➢ The property that provides that information is the volume expansion coefficient 𝜷 ,
𝛿𝑚. 𝑎𝑥 = 𝐹𝑥 (1)
𝛿𝑚 = 𝜌(𝑑𝑥. 𝑑𝑦. 1)
➢ The acceleration in the x-direction
➢ The forces acting on the differential volume element in the vertical direction are the pressure
forces acting on the top and bottom surfaces, the shear stresses acting on the side surfaces
and the force of gravity acting on the entire volume element.
𝜕𝑢
𝜏=𝜇
𝜕𝑦
Equation of motion
(2)
➢ The x-momentum equation in the quiescent fluid outside the boundary layer can be
obtained from the relation above as a special case by setting 𝑢 = 0. It give
➢ For a given x the pressure in the boundary layer is equal to the pressure in the quiescent
fluid. 𝑃 = 𝑃 𝑥 = 𝑃∞ (𝑥)
𝜕𝑃 𝜕𝑃∞
= = −𝜌∞ 𝑔
𝜕𝑥 𝜕𝑥
➢ Eq. 2 will be
(3)
➢ The last term represents the net upward force per unit volume of the fluid. This is the force
that initiates and sustains convection currents.
➢ From the definition of volumetric expansion coefficient
𝜌∞ − 𝜌 = 𝜌𝛽(𝑇 − 𝑇∞ )
➢ Eq. 3 will be
𝜕𝑢 𝜕𝑢 𝜕2𝑢
𝑢 +𝑣 = 𝜈 2 + 𝑔𝛽(𝑇 − 𝑇∞ ) (4)
𝜕𝑥 𝜕𝑦 𝜕𝑦
The Grashof Number
➢ Assume L be the characteristic length and we are introducing nondimensional parameters
𝑥 𝑦 𝑢 𝑣 𝑇 − 𝑇∞
𝑥∗ = 𝑦∗ = 𝑢∗ = 𝑣∗ = 𝑇∗ =
𝐿 𝐿 𝑢𝑜 𝑢𝑜 𝑇𝑠 − 𝑇∞
➢ Substitute in eq. 4
𝑔𝛽(𝑇𝑠 − 𝑇∞ )𝐿3
𝐺𝑟𝐿 =
𝜈2
➢ Grashof number plays the same role in free convection that the Reynolds number plays
in forced convection.
➢ Reynolds number provides a measure of the ratio of the inertial to viscous forces acting
on a fluid element.
➢ In contrast, the Grashof number is a measure of the ratio of the buoyancy forces to the
viscous forces acting on the fluid.
➢ In free convection
𝑁𝑢𝐿 = 𝑓(𝐺𝑟𝐿 , 𝑃𝑟)
Combined Free and Forced Convection
➢ Richardson Number
𝐺𝑟
𝑅𝑖 = 2
𝑅𝑒
➢ The combined effects of free and forced convection must be considered when 𝑅𝑖 ≈ 1
➢ If 𝑅𝑖 ≪ 1 then free convection effects will be neglected
➢ If 𝑅𝑖 ≫ 1 then forced convection effects will be neglected
𝑛 = 3 (𝑇𝑦𝑝𝑖𝑐𝑎𝑙𝑙𝑦)
➢ The effect of buoyancy on heat transfer in a forced flow is strongly influenced by the
direction of the buoyancy force relative to that of the flow.
➢ Three special cases that have been studied extensively correspond to buoyancy-
induced and forced motions having the same direction (assisting flow), opposite
directions (opposing flow), and perpendicular directions (transverse flow).
➢ Upward and downward forced motions over a hot vertical plate are examples of
assisting and opposing flows, respectively.
➢ Examples of transverse flow include horizontal motion over a hot cylinder, sphere, or
horizontal plate.
Governing Equations
Continuity equation
Momentum equation
Energy Equation
Empirical Correlations: External Free Convection Flows
➢ Rayleigh number: 𝑅𝑎𝑥 = 𝐺𝑟𝑥 𝑃𝑟
ത
ℎ𝐿
𝑁𝑢𝐿 = = 𝐶𝑅𝑎𝐿𝑛
𝑘
The Vertical Plate
𝐴𝑠
𝐿=
𝑃
Upper surface heated or lower cooled Lower surface heated or upper cooled
Τ Τ4
𝑁𝑢𝐿 = 0.54𝑅𝑎𝐿1 4 104 ≤ 𝑅𝑎𝐿 ≤ 107 𝑁𝑢𝐿 = 0.27𝑅𝑎𝐿1 105 ≤ 𝑅𝑎𝐿 ≤ 1010
Τ
𝑁𝑢𝐿 = 0.15𝑅𝑎𝐿1 3 107 ≤ 𝑅𝑎𝐿 ≤ 1011
Spheres
Lecture 27
Problem related to free convection
Air flows through a 15 𝑚 long rectangular heating duct that is 0.75 𝑚 wide and
0.3 𝑚 high maintains the outer duct surface at 45𝑜 𝐶. If the duct is uninsulated and
exposed to still air 15𝑜 𝐶, what is the heat loss from the duct?
𝑇∞ = 15𝑜 𝐶
𝑇𝑠 = 45𝑜 𝐶
0.3 𝑚
0.75 𝑚
The average film temperature is: 𝑅𝑎𝐿 = 𝐺𝑟𝐿 𝑃𝑟
𝑅𝑎𝐿 = 7.09 × 107 (Laminar)
15 + 45
𝑇𝑓 = = 30𝑜 𝐶 = 303 𝐾 104 ≤ 𝑅𝑎 ≤ 109
2
For vertical side
𝜈 = 16.2 × 10−6 𝑚2 /𝑠
𝑘 = 0.0265 𝑊/𝑚𝐾
𝑃𝑟 = 0.71
1
𝛽 = = 0.0033 /𝐾 𝑁𝑢𝐿 = 55.13
𝑇
ℎ𝐿
𝑁𝑢 =
For vertical side 𝐿 = 0.3 𝑚 𝑘
𝑁𝑢𝑘
𝑔𝛽(𝑇𝑠 − 𝑇∞ )𝐿3 ℎ=
𝐿
𝐺𝑟𝐿 =
𝜈2
ℎ = 4.87 𝑊/𝑚2 𝐾
9.8 × 0.033(45 − 15) × 0.33 𝑞𝑠𝑖𝑑𝑒 = 2(ℎ𝐴 𝑇𝑠 − 𝑇∞ )
𝐺𝑟𝐿 =
(16.2 × 10−6 )2
𝑞𝑠𝑖𝑑𝑒 = 2(4.87 × 15 × 0.3 × 30)
𝐺𝑟𝐿 = 9.99 × 107
𝑞𝑠𝑖𝑑𝑒 = 1315 𝑊
For horizontal side
Τ4
𝑅𝑎𝐿 = 𝐺𝑟𝐿 𝑃𝑟 𝑁𝑢𝐿 = 0.27𝑅𝑎𝐿1
ℎ𝐿 = 2.096 𝑊/𝑚2 𝐾
Upper surface heated
𝑞𝑏𝑜𝑡𝑡𝑜𝑚 = ℎ𝐴 𝑇𝑠 − 𝑇∞
Τ
𝑁𝑢𝐿 = 0.15𝑅𝑎𝐿1 3
𝑞𝑏𝑜𝑡𝑡𝑜𝑚 = 2.096 × 15 × 0.75 × 30
𝑁𝑢𝐿 = 73.89
𝑞𝑏𝑜𝑡𝑡𝑜𝑚 = 707 𝑊
A horizontal rod 5 𝑚𝑚 in diameter is immersed in water maintained at 18𝑜 𝐶. If the
rod surface temperature is 56𝑜 𝐶, estimate the free convection heat transfer rate per
unit length of the rod.
5 𝑚𝑚
𝑔𝛽(𝑇𝑠 − 𝑇∞ )𝐷3
The average film temperature is: 𝐺𝑟𝐿 =
𝜈2
18 + 56
𝑇𝑓 = = 37𝑜 𝐶 = 310 𝐾 9.8 × 361.9 × 10−6 (56 − 18) × 0.0053
2 𝐺𝑟𝐿 =
(16.2 × 10−6 )2
𝜈 = 6.99 × 10−7 𝑚2 /𝑠
𝐺𝑟𝐿 = 34478.9
𝑘 = 0.268 𝑊/𝑚𝐾
𝑃𝑟 = 4.62 𝑅𝑎𝐿 = 𝐺𝑟𝐿 𝑃𝑟
𝛽 = 361.9 × 10−6 /𝐾
𝑅𝑎𝐿 = 34478.9 × 4.62
𝑅𝑎𝐿 = 159292
𝑞 = ℎ𝐴(𝑇𝑠 − 𝑇∞ )
𝑞 = ℎ × (𝜋𝐷𝐿) × (𝑇𝑠 − 𝑇∞ )
𝑁𝑢 = 10.40
′
𝑞
𝑞 = = ℎ × (𝜋𝐷) × (𝑇𝑠 − 𝑇∞ )
ℎ𝐷 𝐿
𝑁𝑢 =
𝑘
𝑞 ′ = 1306 × (𝜋 × 0.005) × (56 − 18)
𝑁𝑢𝑘
ℎ=
𝐷
10.40 × 0.628
𝑞′ = 780 𝑊/𝑚
ℎ=
0.005
ℎ = 1306 𝑊/𝑚2 𝐾
Airflow through a long, 0.2 𝑚 square air conditioning duct maintains the outer duct
surface temperature at 10𝑜 𝐶. If the horizontal duct is uninsulated and exposed to air
at 35𝑜 𝐶 in the crawlspace beneath a home, what is the heat gain per unit length of
the duct?
𝑇∞ = 35𝑜 𝐶
𝑇𝑠 = 10𝑜 𝐶
0.2 𝑚
0.2 𝑚
The average film temperature is:
𝑞𝑠𝑖𝑑𝑒 = 2 × ℎ𝐴(𝑇𝑠 − 𝑇∞ )
𝑅𝑎𝐿 = 𝐺𝑟𝐿 𝑃𝑟
10 + 35
𝑇𝑓 = = 22.5𝑜 𝐶 = 295.5 𝐾
2 𝑞′𝑠𝑖𝑑𝑒 = 2 × 4.92 × 0.2 × (35 − 10)
𝑅𝑎𝐿 = 25.61 × 106 × 0.71
𝜈 = 15.89 × 10−6 𝑚2 /𝑠
𝑅𝑎𝐿 = 18.11 × 106 𝑞′𝑠𝑖𝑑𝑒 = 49.2 𝑊/𝑚
𝑘 = 0.0265 𝑊/𝑚𝐾
𝑃𝑟 = 0.71
1
𝛽 = = 0.0033 /𝐾
𝑇
ℎ = 4.92 𝑊/𝑚2 𝐾
𝐺𝑟𝐿 = 25.61 × 106
𝑁𝑢𝑘
ℎ=
For horizontal side 𝐿
𝐴𝑠 ℎ𝐿 = 5.51 𝑊/𝑚2 𝐾
𝐿= 𝑞𝑡𝑜𝑡𝑎𝑙 = 𝑞𝑡𝑜𝑝 + 𝑞𝑏𝑜𝑡𝑡𝑜𝑚 + 𝑞𝑠𝑖𝑑𝑒
𝑃
𝑞𝑡𝑜𝑝 = ℎ𝐴 𝑇𝑠 − 𝑇∞
𝐿 = 0.1 𝑚 𝑞𝑡𝑜𝑡𝑎𝑙 = 90.5 𝑊/𝑚
𝑞′𝑡𝑜𝑝 = 5.5 × 0.2 × 25
𝑔𝛽(𝑇𝑠 − 𝑇∞ )𝐿3
𝐺𝑟𝐿 = 𝑞𝑡𝑜𝑝 = 27.55 𝑊/𝑚
𝜈2
𝐺𝑟𝐿 = 3.202 × 106 Bottom surface heated
Τ4
𝑅𝑎𝐿 = 𝐺𝑟𝐿 𝑃𝑟 𝑁𝑢𝐿 = 0.27𝑅𝑎𝐿1
ℎ𝐿 = 2.75 𝑊/𝑚2 𝐾
Upper surface heated
𝑞𝑏𝑜𝑡𝑡𝑜𝑚 = ℎ𝐴 𝑇𝑠 − 𝑇∞
Τ
𝑁𝑢𝐿 = 0.54𝑅𝑎𝐿1 4
𝑞′𝑏𝑜𝑡𝑡𝑜𝑚 = 2.75 × 0.2 × 25
𝑁𝑢𝐿 = 20.96
𝑞𝑏𝑜𝑡𝑡𝑜𝑚 = 13.75 𝑊/𝑚
Lecture 28
Boiling and Condensation
Boiling
➢ When evaporation occurs at a solid liquid interface, it is termed boiling.
➢ The process occurs when the temperature of the surface 𝑇𝑠 exceeds the
saturation temperature 𝑇𝑠𝑎𝑡 corresponding to the liquid pressure.
➢ Heat is transferred from the solid surface to the liquid, and the appropriate
form of Newton’s law of cooling is
➢ It is important to acknowledge, however, that above equation applies only to clean, unfinned surfaces.
➢ During normal heat exchanger operation, surfaces are often subject to fouling by fluid impurities, rust formation, or other
reactions between the fluid and the wall material.
➢ This effect can be treated by introducing an additional thermal resistance in equation 1, termed the fouling factor, 𝑅𝑓 .
➢ where i is the fluid enthalpy. The subscripts h and c refer to the hot and cold
fluids, whereas the subscripts i and o designate the fluid inlet and outlet
conditions.
➢ If the fluids are not undergoing a phase change and constant specific heats are
assumed, these expressions reduce to
since ∆𝑇 varies with position in the heat exchanger, it is necessary to work with a
rate equation of the form
The Parallel-Flow Heat Exchanger
➢ The energy balances and the subsequent analysis are subject to the following assumptions.
➢ where 𝑪𝒉 and 𝑪𝒄 are the hot and cold fluid heat capacity rates, respectively
➢ Heat transfer across the surface area 𝑑𝐴 may also be expressed as
𝒅𝒒 = 𝑼∆𝑻𝒅𝑨 5
∆𝑻 = 𝑻𝒉 − 𝑻𝒄
𝒅(∆𝑻) = 𝒅𝑻𝒉 − 𝒅𝑻𝒄
➢ From eq. 3 and 4
𝟏 𝟏
𝒅 ∆𝑻 = −𝒅𝒒 +
𝑪𝒉 𝑪𝒄
The Parallel-Flow Heat Exchanger
𝟏 𝟏 𝟏 𝟏
𝑑 ∆𝑇 = −𝒅𝒒 + 𝑑 ∆𝑇 = −𝑈∆𝑇𝑑𝐴 +
𝑪𝒉 𝑪𝒄 𝑪𝒉 𝑪𝒄
∆𝑇2 2
𝑑 ∆𝑇 𝟏 𝟏 ∆𝑇2 1 1
න =𝑈 + න 𝑑𝐴 ln = −𝑈𝐴 + 6
∆𝑇1 ∆𝑇 𝑪𝒉 𝑪𝒄 1 ∆𝑇1 𝐶ℎ 𝐶𝑐
𝑞 = 𝑈𝐴∆𝑇𝑙𝑚
The Counterflow Heat Exchanger
➢ Note that, for the same inlet and outlet temperatures, the log mean temperature difference for counterflow exceeds that for
parallel flow
➢ Hence the surface area required to effect a prescribed heat transfer rate q is smaller for the counterflow than for the parallel-
flow arrangement, assuming the same value of U.
➢ Also note that 𝑇𝑐,𝑜 can exceed 𝑇ℎ,𝑜 for counterflow but not for parallel flow.
Special Operating Conditions
➢ If the hot fluid has a heat capacity rate, much larger than that of the cold fluid then the temperature of the hot fluid remains
approximately constant throughout the heat exchanger, while the temperature of the cold fluid increases.
➢ The same condition is achieved if the hot fluid is a condensing vapor.
➢ Condensation occurs at constant temperature, and, for all practical purposes, 𝐶ℎ → ∞.
➢ In an evaporator or a boiler, it is the cold fluid that experiences a change in phase and remains at a nearly uniform temperature
(𝐶𝑐 → ∞).
➢ The third special case involves a counterflow heat exchanger for which the heat capacity rates are equal (𝐶ℎ = 𝐶𝑐 ).
Lecture 30
Heat Exchanger (Cont.)
Heat Exchanger Analysis: The Effectiveness–NTU Method
➢ To define the effectiveness of a heat exchanger, we must first determine the maximum possible heat transfer rate, 𝑞𝑚𝑎𝑥 , for the
exchanger.
➢ Consider a counterflow heat exchanger of infinite length and also 𝐶𝑐 < 𝐶ℎ .
➢ The cold fluid would then experience the larger temperature change, and since 𝐿 → ∞, it would be heated to the inlet
temperature of the hot fluid (𝑇𝑐,𝑜 = 𝑇ℎ,𝑖 )
𝒒𝒎𝒂𝒙 = 𝑪𝒄 (𝑻𝒉,𝒊 − 𝑻𝒄,𝒊 )
➢ If 𝐶ℎ < 𝐶𝑐 the hot fluid would experience the larger temperature change and would be cooled to the inlet temperature of the cold
fluid (𝑇ℎ,𝑜 = 𝑇𝑐,𝑖 )
𝒒𝒎𝒂𝒙 = 𝑪𝒉 (𝑻𝒉,𝒊 − 𝑻𝒄,𝒊 )
➢ From the foregoing results we are then prompted to write the general expression
𝒒𝒎𝒂𝒙 = 𝑪𝒎𝒊𝒏 (𝑻𝒉,𝒊 − 𝑻𝒄,𝒊 )
➢ The effectiveness, 𝜀 , is the ratio of the actual heat transfer rate for a heat exchanger to the maximum possible heat transfer rate:
𝒒
𝜺=
𝒒𝒎𝒂𝒙
𝑪𝒉 (𝑻𝒉,𝒊 − 𝑻𝒉,𝒐 ) 𝑪𝒄 (𝑻𝒄,𝒐 − 𝑻𝒄,𝒊 )
𝜺= 𝜺=
𝑪𝒎𝒊𝒏 (𝑻𝒉,𝒊 − 𝑻𝒄,𝒊 ) 𝑪𝒎𝒊𝒏 (𝑻𝒉,𝒊 − 𝑻𝒄,𝒊 )
➢ The number of transfer units (NTU) is a dimensionless parameter that is widely used for heat exchanger analysis and is defined
as
𝑈𝐴
𝑁𝑇𝑈 =
𝐶𝑚𝑖𝑛
𝐶𝑚𝑖𝑛
𝐶𝑟 =
𝐶𝑚𝑎𝑥
A shell-and-tube exchanger (two shells, four tube passes) is used to heat 10,000 𝑘𝑔/ℎ of
pressurized water from 35𝑜 to 120𝑜 𝐶 with 5000 𝑘𝑔/ℎ pressurized water entering the exchanger
at 300 𝑜 𝐶. If the overall heat transfer coefficient is 1500 𝑊/𝑚2 𝐾, determine the required heat
exchanger area.
𝑚ሶ 𝑐 = 10000 𝑘𝑔Τℎ𝑟 = 2.77 𝑘𝑔/𝑠 𝑞 = 𝑚ሶ 𝑐 𝐶𝑝,𝑐 (𝑇𝑐,𝑜 − 𝑇𝑐,𝑖 )
∆𝑇𝑙𝑚 = 143𝑜 𝐶
𝑞 = 𝑈𝐴∆𝑇𝑙𝑚
𝑞
𝐴=
𝑈 ∆𝑇𝑙𝑚
A counterflow, concentric tube heat exchanger is used to cool the lubricating oil for a large
industrial gas turbine engine. The flow rate of cooling water through the inner tube (𝐷𝑖 = 25 𝑚𝑚)
is 0.2 𝑘𝑔/𝑠, while the flow rate of oil through the outer annulus (𝐷𝑜 = 45 𝑚𝑚) is 0.1 𝑘𝑔/𝑠. The
oil and water enter at temperatures of 100𝑜 and 30𝑜 𝐶, respectively. How long must the tube be
made if the outlet temperature of the oil is to be 60𝑜 𝐶?
At
𝑚ሶ 𝑐 = 0.2 𝑘𝑔/𝑠 30 + 40
𝑇𝑐,𝑚 = = 35𝑜 𝐶 = 308𝐾
𝑚ሶ ℎ = 0.1 𝑘𝑔/𝑠 2
𝐶𝑃,𝑐 = 4187 𝐽/𝑘𝑔𝐾
𝑇ℎ,𝑖 = 100𝑜 𝐶 𝑇ℎ,𝑜 = 60𝑜 𝐶
𝑜 𝜇 = 725 × 10−6 𝑁𝑠/𝑚2
𝑇𝑐,𝑖 = 30 𝐶
𝑘 = 0.625 𝑊/𝑚𝐾
At
𝑃𝑟 = 4.85
100 + 60 Amount of heat lost by the hot fluid
𝑇ℎ,𝑚 = = 80𝑜 𝐶 = 353𝐾
2
𝑞 = 𝑚ሶ ℎ 𝐶𝑝,ℎ (𝑇ℎ,𝑖 − 𝑇ℎ,𝑜 )
𝐶𝑃,ℎ = 2131 𝐽/𝑘𝑔𝐾
𝑞 = 0.1 × 2131 × (100 − 60)
𝜇 = 3.25 × 10−3 𝑁𝑠/𝑚2
𝑞 = 8524 𝑊
𝑘 = 0.138 𝑊/𝑚𝐾 Same amount of heat will be gained by cold fluid
Assume 𝑇𝑐,𝑜 = 40𝑜 𝐶 𝑞 = 𝑚ሶ 𝑐 𝐶𝑝,𝑐 (𝑇𝑐,𝑜 − 𝑇𝑐,𝑖 )
𝑇𝑐,𝑜 = 40.1𝑜 𝐶
For flow through tube
𝑞 = 𝑈𝐴∆𝑇𝑙𝑚
4𝑚ሶ 𝑐
𝑅𝑒𝐷 =
𝜋𝐷𝑖 𝜇
(𝑇ℎ,𝑖 − 𝑇𝑐,𝑜 ) − (𝑇ℎ,𝑜 − 𝑇𝑐,𝑖 )
∆𝑇𝑙𝑚 = 4 × 0.2
(𝑇 − 𝑇𝑐,𝑜 ) 𝑅𝑒𝐷 =
ln ℎ,𝑖 𝜋 × 0.025 × 725 × 10−6
(𝑇ℎ,𝑜 − 𝑇𝑐,𝑖 )
𝑅𝑒𝐷 = 14050 > 2300
(100 − 40.1) − (60 − 30) Turbulent flow
∆𝑇𝑙𝑚 =
(100 − 40.1) Using Dittus Boelter equation and water is heated
ln
(60 − 30)
𝑁𝑢 = 0.023𝑅𝑒 4/5 𝑃𝑟 0.4
∆𝑇𝑙𝑚 = 43.2𝑜 𝐶
𝑁𝑢 = 90
The over all heat transfer coefficient will be ℎ𝑖 𝐷𝑖
𝑁𝑢 =
𝑘
1 1 1
= + ℎ𝑖 = 2250 𝑊/𝑚2 𝐾
𝑈 ℎ𝑖 ℎ𝑜
For outer tube flow (annulus) 1 1 1
= +
𝜌𝑢𝑚 𝐷ℎ 𝑈 ℎ𝑖 ℎ𝑜
𝑅𝑒 = 𝐴
𝜇 1 1 1 𝐿=
= + 𝜋𝐷𝑖
𝜌 (𝐷𝑜 − 𝐷𝑖 ) 𝑚ሶ ℎ 𝑈 2250 38.847
𝑅𝑒 = × 𝜋
𝜇 𝜌 × 4 × (𝐷𝑜2 − 𝐷𝑖4 ) 𝐿 = 65.9 𝑚
2
𝑈 = 38.187 𝑊/𝑚 𝐾
𝑅𝑒 = 56 < 2300
𝑞 = 𝑈𝐴∆𝑇𝑙𝑚
Flow is laminar 𝑞
𝐴=
𝑈∆𝑇𝑙𝑚
𝐷𝑖
= 0.56 8524
𝐷𝑜 𝐴=
38.187 × 43.2
Using table 8-2
𝐴 = 5.177 𝑚2
ℎ𝑜 𝐷ℎ ℎ𝑜 (𝐷𝑜 −𝐷𝑖 )
𝑁𝑢 = = = 5.63 𝐴 = 𝜋𝐷𝑖 𝐿
𝑘 𝑘
ℎ𝑜 = 38.847 𝑊/𝑚2 𝐾
Hot exhaust gases, which enter a finned-tube, cross-flow heat exchanger at 300𝑜 𝐶 and leave at
100𝑜 𝐶, are used to heat pressurized water at a flow rate of 1 𝑘𝑔/𝑠 from 35𝑜 to 125𝑜 𝐶. The
overall heat transfer coefficient based on the gas-side surface area is 𝑈ℎ = 100 𝑊/𝑚2 𝐾.
Determine the required gas-side surface area 𝐴ℎ using the 𝑁𝑇𝑈 method.
The amount of heat lost by hot fluid is same as that of gained by cold fluid
35 + 125
𝑇ത𝑐 = = 80𝑜 𝐶
2
𝑞ℎ = 𝑞𝑐
At
𝑇ത𝑐 = 80𝑜 𝐶 𝐶ℎ 𝑇ℎ,𝑖 − 𝑇ℎ,𝑜 = 𝐶𝑐 (𝑇𝑐,𝑜 − 𝑇𝑐,𝑖 )
𝐶ℎ = 1889 𝑊/𝐾
𝐶𝑚𝑖𝑛
𝑞𝑚𝑎𝑥 = 𝐶𝑚𝑖𝑛 (𝑇ℎ,𝑖 − 𝑇𝑐,𝑖 ) 𝜀 = 𝑓 𝑁𝑇𝑈,
𝐶𝑚𝑎𝑥
𝑞𝑚𝑎𝑥 = 1889(300 − 35)
➢ Combining Equations
➢ The total molar transfer rate for an entire surface, 𝑁𝐴 (𝑘𝑚𝑜𝑙/𝑠), may then be expressed as
𝑁𝐴 = ℎ𝑚 𝐴𝑠 (𝐶𝐴,𝑠 − 𝐶𝐴,∞ ) 3
➢ where the average and local mass transfer convection coefficients are related by an equation of the form
➢ A Schmidt number of near unity (Sc = 1) indicates that momentum and mass
transfer by diffusion are comparable, and velocity and concentration boundary
layers almost coincide with each other.
Lewis number
➢ It is the ratio of thermal diffusivity to the mass diffusivity.
➢ It represent the relative magnitudes of heat and mass diffusion in the thermal
and concentration boundary layers.
Dimensionless parameters in Mass Transfer
Sherwood number:
➢ It is the dimensionless concentration gradient at the surface.
➢ Nusselt number to Heat Transfer as Sherwood number to mass
transfer.
Stanton number
➢ Sometimes it is more convenient to express the heat and mass
transfer coefficients in terms of the dimensionless Stanton
number as
Dimensionless parameters in Mass Transfer
➢ For a given geometry, the average Nusselt number in forced convection
depends on the Reynolds and Prandtl numbers, whereas the average
Sherwood number depends on the Reynolds and Schmidt numbers.
5 × 500
𝑅𝑒 = 𝑆ℎ = 0.037𝑅𝑒 0.8 𝑆𝑐 1/3
15.89 × 10−6
𝑆ℎ = 1.13 × 105
𝑅𝑒 = 157.3 × 106 > 5 × 105
Flow is turbulent The convective mass transfer will be determined by
15.89 × 10−6 ℎ𝑚 𝐿
𝑥𝑐𝑟 = 5 × 105 × 𝑆ℎ =
5 𝐷𝐴𝐵
𝑥𝑐𝑟 = 1.589 𝑚
Assume turbulent flow completely ℎ𝑚 = 5.87 × 10−3 𝑚/𝑠
The convective mass transfer will be determined by Concentration at air side
𝑛𝐴 = ℎ𝑚 𝐴(𝜌𝐴,𝑠 − 𝜌𝐴,∞ )
𝐶𝐴,∞ = 1.415 × 10−4 𝑘𝑚𝑜𝑙/𝑚3
ℎ𝑚 𝐿 𝑛𝐴 = 0.029 × 1(0.436 − 0)
𝑆ℎ =
𝐷𝐴𝐵
𝑛𝐴 = 1.262 × 10−3 𝑘𝑔/𝑠
ℎ𝑚 = 0.029 𝑚/𝑠
𝑛𝐴 = 4.55 𝑘𝑔/ℎ𝑟
The density at air side : 𝜌𝐴,∞ = 0
1
𝜌𝐴,𝑆 =
𝑣𝑔
1
𝜌𝐴,𝑆 = = 0.436 𝑘𝑔/𝑚3
22.93
A stream of air at 150 KPa and 300 K is flowing over a top surface of a thin flat
sheet of solid naphthalene of length 0.2 m with velocity of 20 m/s. If diffusivity of
naphthalene in air at a given condition is 4 × 10−6 𝑚2 /𝑠 and concentration of
naphthalene is 1 × 10−5 𝐾𝑚𝑜𝑙 Τ𝑚3 . Find the rate of loss of naphthalene from the
surface per unit width.
𝜈 = 15.89 × 10−6 𝑚2 Τ𝑠
Air at 150 KPa The Schmidt number will be determined by
𝜈𝑎𝑖𝑟
101325 𝑆𝑐 =
𝜈 = 15.89 × 10−6 × 𝐷
150 × 103 𝑆𝑐 = 2.6825
𝜈 = 10.73 × 10−6 𝑚2 Τ𝑠 The Sherwood number will be determined by
ℎ𝑚 𝐿
𝑆ℎ =
𝐷𝐴𝐵
ℎ𝑚 = 0.01127 𝑚/𝑠
𝐶𝐴,∞ = 0