1.5 Unsteady Heat Conduction

Download as pdf or txt
Download as pdf or txt
You are on page 1of 6

This energy must be conducted to the surface, and therefore,

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.

For Bi << 1, it is reasonable to assume a uniform temperature distribution across a solid


at any time during a transient process.
Founer Modulus - It is also a dimensionless number and is define as

Fo= α t/L2 (3.2)


Where, L is the characteristic length of the body, a is the thermal diffusivity, and t is the
time

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.

Fig. 1.18 (a) Nomenclature for Biot Modulus

1.7.3. Lumped Capacity System-Necessary Physical Assumptions

We know that a temperature gradient must exist in a material if heat energy is to be


conducted into or out of the body. When Bi < 0.1, it is assumed that the internal thermal
resistance of the body is very small in comparison with the external resistance and the transfer of
heat energy is primarily controlled by the convective heat transfer at the surface. That is, the
temperature within the body is approximately uniform. This idealized assumption is possible, if

(a) The physical size of the body is very small,

(b) The thermal conductivity of the material is very large, and

(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.

7.4. An Expression for Evaluating the Temperature Variation in a Solid Using


Lumped Capacity Analysis

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

during a time dt the body during that time dt

or, hAs ( T − T∞ ) dt = - pCVdT

Where As is the surface area of the body, V is the volume of the body and C is the
specific heat capacity.

or, (hA/ ρ CV) dt = - dT /( T − T∞ )

with the initial condition being: at t = 0, T = Ts

The solution is: ( T − T∞ )/( Ts − T∞ ) = exp(-hA / ρ CV)t (3.3)

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

Bi × Fo = (hL/k) × ( α t/L2) = (hL/k) (k/ ρ C) (t/L2) = (hA/ ρ CV) t,


Where, V / A is the characteristic length L.

And, the solution describing the temperature variation of the object with respect to time
is given by

( T − T∞ )/ ( Ts − T∞ ) = exp (-Bi· Fo) (3.4)

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.

Solution: Characteristic length, L = VIA = 4/3 π r3/4 π r2 = r/3 = 5 × 10-3/3m

Bi = hL/k = 25 × 5 × 10-3/ (3 × 48) = 8.68 × 10-4<< 0.1,

Since the internal resistance is negligible, we make use of lumped capacity analysis: Eq.
(3.4),

( T − T∞ ) / ( Ts − T∞ )=exp(-Bi Fo) ; (150 35) / (750 35) = 0.16084

∴ Bi × Fo = 1827; Fo = 1.827/ (8. 68 × 10-4) 2.1× 103

or, α t/ L2 = k/ ( ρ CL2) t = 2100 and t = 568 = 0.158 hour

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

= −ρCV ( Ts − T∞ ) exp ( − hAt / ρCV ) − 1 (3.5)

= -7800 × 600 × (4/3) π (5 × 10-3)3 (750-35) (0.16084 - 1)

= 1.47 × 103 J = 1.47 kJ.

If we allow the time’t’ to go to infinity, we would have a situation that corresponds to


steady state in the new environment. The change in internal energy will be U0 - U∞ =

[ ρ CV( Ts − T∞ ) exp (- ∞ ) - 1] = [ ρ CV( Ts − T∞ ].

We can also compute the instantaneous heal transfer rate at any time.

or. Q = - ρ VCdT/dt = - ρ VCd/dt [ T∞ + ( Ts − T∞ ) exp (-hAt/ ρ CV)]


= hA( Ts − T∞ )[exp(-hAt/ ρ CV)) and for t = 60s,

Q = 25 × 4 × 3.142 (5 × 10-3)2(750 35) [exp (-25 × 3 × 60/5 × 10-3 × 7800 × 600)]

= 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.

Solution: Characteristic length. L = V/A = π r2L/2 π r(r+ L) = rL/2(r + L)

= 5 × 10-2 × 30 × 10-2/2 (2 (5 +30) × 10-2)

= 2.143 × 10-2 m.

Bi = hL/k = 100 × 2.143 × 10-2/40 = 0.0536 << 0.1

Fo = α t/L2 = (k/ ρ C) × (t/L2)

= 40 × t/ (7600 × 600 × [2.143 × 1 0-2)2] = 191 × 10-2 t

and ( T − T∞ )/ ( Ts − T∞ ) = exp (-Bi .Fo)

or, (850 - 1300) / (50 - 1300) = 0.36 = exp (- Bi Fo)

∴ Bi Fo= 102

and Fo = 19.06 and t = 19.06/( 1.91 × 1 0 -2) = 16.63 min

(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.

Solution: Characteristic length, L= [2 × 3 × 4 /2(2 × 3 + 2 × 4 + 3 × 4)] × 10-2

= 4.6 × 10-3m

Bi = hL/k = 50 × 4.6 × 10-3/180 = 1.278 × 10-3 << 0.1


Fo = α t/L2 =10-4 × 180 / (4.6×10 -3)2 = 850

exp (-Bi Fo) = exp (-1.278 × 10-3×'850) = 0.337

(T - T ∞ ) (Ts - T ∞ ) - (T - 30)/(300 - 30) = 0.337

∴ 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,

Bi = hL/k = 100 × 0.5 × 10 -3/ (2 × 370) = 6.76 × 10-5<< 0.1

Fo = α t/L2 = (k/ ρC ) × (t/L2)

= [370t/ (8800 × 381 × (0.25 × 10-3)2] = 1760t

exp (-Bi Fo) = (T - T ∞ )(Ts -T ∞ )

= (90 - 35)/(150 - 35) = 0.478

Bi Fo = 0.738 = 6.76 × 10-5 × 1760 t; ∴ t = 6.2s

when h = 40 W/ m2K, Bi = 2.7 × 10 -5 and 2.7 × 10 -5 ×1760 t = 0.738;

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.

Solution: hA/ ρCV = hA/ mC = h × 40 × 10-4 /(0.1 × 350) = 1.143×10-4h

and, hAt / ρCV = 1.143 × 10-4 h × 100 = 1.143 × 10 -2h

(T - T ∞ )/ (Ts - T ∞ ) = (40 - 25) / (100 - 25) = 0.2

∴ exp (-1.143 × 10-2h) = 0.2 or, 1.143 × 10-2h = 1.6094, and h = 140W/m2K.

You might also like