1.5 Unsteady Heat Conduction
1.5 Unsteady Heat Conduction
1.5 Unsteady Heat Conduction
Q = -kA(dT / dr) r = R
or, h ( Ts − T∞ ) = -k (dT/dr) r=R ≈ -k (Tc-Ts)/R
Where, Tc is the temperature at the axis of the cylinder
By rearranging, (Ts - Tc) / ( Ts − T∞ ) h/Rk (3.1)
The term, hR/k, IS called the 'BlOT MODULUS'. It is a dimensionless number and is
the ratio of internal heat flow resistance to external heat flow resistance and plays a fundamental
role in transient conduction problems involving surface convection effects. I t provides a
measure 0 f the temperature drop in the solid relative to the temperature difference between the
surface and the fluid.
The Fourier modulus measures the magnitude of the rate of conduction relative to the
change in temperature, i.e., the unsteady effect. If Fo << 1, the change in temperature will be
experienced by a region very close to the surface.
Fig. 1.18 Effect of Biot Modulus on steady state temperature distribution in a plane wall
with surface convection.
(c) The convective heat transfer coefficient at the surface is very small and there is a
large temperature difference across the fluid layer at the interface.
Let us consider a small metallic object which has been suddenly immersed in a fluid
during a heat treatment operation. By applying the first law of
Heat flowing out of the body = Decrease in the internal thermal energy of
Where As is the surface area of the body, V is the volume of the body and C is the
specific heat capacity.
Fig. 3.2 depicts the cooling of a body (temperature distribution time) using lumped
thermal capacity system. The temperature history is seen to be an exponential decay.
We can express
And, the solution describing the temperature variation of the object with respect to time
is given by
Example 1.19 Steel balls 10 mm in dia (k = 48 W/mK), (C = 600 J/kgK) are cooled in
air at temperature 35°C from an initial temperature of 750°C. Calculate the time required for the
temperature to drop to 150°C when h = 25 W/m2K and density p = 7800 kg/m3.
Since the internal resistance is negligible, we make use of lumped capacity analysis: Eq.
(3.4),
We can also compute the change in the internal energy of the object as:
1 1
U 0 − U t = − ∫ ρCVdT = ∫ ρCV ( Ts − T∞ ) ( −hA / ρCV ) exp t ( −hAt / ρCV ) dt
0 0
We can also compute the instantaneous heal transfer rate at any time.
= 4.63 W.
Example 1.20 A cylindrical steel ingot (dia 10 cm. length 30 cm, k = 40 W mK. ρ =
7600 kg/m3, C = 600 J/kgK) is to be heated in a furnace from 50°C to 850°C. The temperature
inside the furnace is 1300oC and the surface heat transfer coefficient is 100 W/m2K. Calculate
the time required.
= 2.143 × 10-2 m.
∴ Bi Fo= 102
(The length of the ingot is 30 cm and it must be removed from the furnace after a period
of 16.63 min. therefore, the speed of the ingot would be 0.3/16.63 = 1.8 × 10-2 m/min.)
Example 1.21 A block of aluminum (2cm × 3cm × 4cm, k = 180 W/mK, α = 10 -4m2/s)
initially at 300oC is cooled in air at 30oC. Calculate the temperature of the block after 3 min.
Take h = 50W/ m2K.
= 4.6 × 10-3m
∴ T= 121.1°C.
Example 1.22 A copper wire 1 mm in dia initially at 150°C is suddenly dipped into
water at 35°C. Calculate the time required to cool to a temperature of 90°C if h = 100 W/ m2K.
What would be the time required if h = 40 W/m2K. (for copper; k = 370 W/mK, ρ = 8800
kg/m3. C = 381 J/kgK.
Solution: The characteristic length for a long cylindrical object can be approximated as
r/2. As such,
or, t = 15.53s.
Example 1.23 A metallic rod (mass 0.1 kg, C = 350 J/kgK, dia 12.5 mm, surface area
40cm2) is initially at 100°C. It is cooled in air at 25°C. If the temperature drops to 40°C in 100
seconds, estimate the surface heat transfer coefficient.
∴ exp (-1.143 × 10-2h) = 0.2 or, 1.143 × 10-2h = 1.6094, and h = 140W/m2K.