PROBLEM 5.

5
KNOWN: Diameter and initial temperature of steel balls cooling in air.
FIND: Time required to cool to a prescribed temperature.
SCHEMATIC:

ASSUMPTIONS: (1) Negligible radiation effects, (2) Constant properties.

ANALYSIS: Applying Eq. 5.10 to a sphere (Lc = ro/3),

hLc h ( ro / 3) 20 W/m ⋅ K (0.002m )

2
Bi = = = = 0.001.
k k 40 W/m ⋅ K

Hence, the temperature of the steel remains approximately uniform during the cooling process, and
the lumped capacitance method may be used. From Eqs. 5.4 and 5.5,

t=
ρ Vcp
ln =
3
( ln
)
Ti − T∞ ρ π D / 6 c p Ti − T∞
hAs T − T∞ hπ D2 T − T∞

7800kg/m 3 ( 0.012m ) 600J/kg ⋅ K 1150 − 325

t= ln
6 × 20 W/m 2 ⋅ K 400 − 325

t = 1122 s = 0.312h <

COMMENTS: Due to the large value of Ti, radiation effects are likely to be significant during the
early portion of the transient. The effect is to shorten the cooling time.
 PROBLEM 5.7
KNOWN: Solid steel sphere (AISI 1010), coated with dielectric layer of prescribed thickness and
thermal conductivity. Coated sphere, initially at uniform temperature, is suddenly quenched in an oil
bath.
FIND: Time required for sphere to reach 140°C.
SCHEMATIC:

(
PROPERTIES: Table A-1, AISI 1010 Steel T = [500 +140 ] C/2 = 320oC ≈ 600K :
o
)
ρ = 7832 kg/m 3 , c = 559 J/kg ⋅ K, k = 48.8 W/m ⋅ K.
ASSUMPTIONS: (1) Steel sphere is space-wise isothermal, (2) Dielectric layer has negligible
thermal capacitance compared to steel sphere, (3) Layer is thin compared to radius of sphere, (4)
Constant properties.
ANALYSIS: The thermal resistance to heat transfer from the sphere is due to the dielectric layer and
the convection coefficient. That is,
l 1 0.002m 1 m2 ⋅ K
R ′′ = + = + = ( 0.050 + 0.0003 ) = 0.0503 ,
k h 0.04 W/m ⋅ K
3300 W/m 2 ⋅ K W
or in terms of an overall coefficient, U = 1/R′′ = 19.88 W/m2 ⋅ K. The effective Biot number is
ULc U ( ro /3 ) 19.88 W/m ⋅ K × ( 0.300/6 ) m
2
Bi e = = = = 0.0204
k k 48.8 W/m ⋅ K
where the characteristic length is Lc = ro/3 for the sphere. Since Bie < 0.1, the lumped capacitance
approach is applicable. Hence, Eq. 5.5 is appropriate with h replaced by U,
ρ c  V  θ i ρ c  V  T ( 0 ) − T∞
t= =
U  As  θ o U  As  T ( t ) − T∞
ln ln .

Substituting numerical values with (V/As ) = ro/3 = D/6,

7832 kg/m 3 × 559 J/kg ⋅ K  0.300m  ( 500 −100 ) C
o
t=  6  ln
19.88 W/m 2 ⋅ K (140 − 100 )
o
C

t = 25,358s = 7.04h. <

COMMENTS: (1) Note from calculation of R′′ that the resistance of the dielectric layer dominates
and therefore nearly all the temperature drop occurs across the layer.
 PROBLEM 5.9
KNOWN: Diameter and radial temperature of AISI 1010 carbon steel shaft. Convection
coefficient and temperature of furnace gases.
FIND: Time required for shaft centerline to reach a prescribed temperature.
SCHEMATIC:

ASSUMPTIONS: (1) One-dimensional, radial conduction, (2) Constant properties.

PROPERTIES: AISI 1010 carbon steel, Table A.1 ( T = 550 K ) : ρ = 7832 kg / m 3 , k = 51.2
-5 2
W/m⋅K, c = 541 J/kg⋅K, α = 1.21×10 m /s.
ANALYSIS: The Biot number is

hro / 2 100 W/m ⋅ K ( 0.05 m/2)

2
Bi = = = 0.0488.
k 51.2 W/m ⋅ K
Hence, the lumped capacitance method can be applied. From Equation 5.6,
T − T∞   hAs    4h 
= exp  −   t  = exp  − t
Ti − T∞   ρ Vc    ρ cD 
 800 − 1200  4 × 100 W/m 2 ⋅ K
ln   = −0.811 = − t
 300 − 1200  7832 kg/m 3 ( 541 J/kg ⋅ K ) 0.1 m
t = 859 s. <
COMMENTS: To check the validity of the foregoing result, use the one-term approximation to the
series solution. From Equation 5.49c,
To − T∞ −400
=
Ti − T∞ −900
= 0.444 = C1 exp −ς12 Fo ( )
For Bi = hro/k = 0.0976, Table 5.1 yields ς1 = 0.436 and C1 = 1.024. Hence

− ( 0.436)
2
(1.2×10−5 m2 / s) t = ln (0.434) = −0.835
( 0.05 m ) 2
t = 915 s.
The results agree to within 6%. The lumped capacitance method underestimates the actual time,
since the response at the centerline lags that at any other location in the shaft.
 PROBLEM 5.11
KNOWN: Diameter, density, specific heat and thermal conductivity of aluminum spheres used in
packed bed thermal energy storage system. Convection coefficient and inlet gas temperature.
FIND: Time required for sphere to acquire 90% of maximum possible thermal energy and the
corresponding center temperature. Potential advantage of using copper in lieu of aluminum.
SCHEMATIC:

ASSUMPTIONS: (1) Negligible heat transfer to or from a sphere by radiation or conduction due to
contact with other spheres, (2) Constant properties.
ANALYSIS: To determine whether a lumped capacitance analysis can be used, first compute Bi =
2
h(ro/3)/k = 75 W/m ⋅K (0.025m)/150 W/m⋅K = 0.013 < 0.1. Hence, the lumped capacitance
approximation may be made, and a uniform temperature may be assumed to exist in the sphere at any
time. From Eq. 5.8a, achievement of 90% of the maximum possible thermal energy storage
corresponds to
Q
= 0.90 = 1 − exp ( − t / τ t )
ρ cVθi

where τ t = ρ Vc / hA s = ρ Dc / 6h = 2700 kg / m3 × 0.075m × 950 J / kg ⋅ K / 6 × 75 W / m 2 ⋅ K = 427s. Hence,

t = −τ t ln ( 0.1) = 427s × 2.30 = 984s <

From Eq. (5.6), the corresponding temperature at any location in the sphere is

( )
T (984s ) = Tg,i + Ti − Tg,i exp ( −6ht / ρ Dc )

(
T (984s ) = 300°C − 275°C exp −6 × 75 W / m ⋅ K × 984s / 2700 kg / m × 0.075m × 950 J / kg ⋅ K
2 3
)
T (984 ) s = 272.5°C <
3
Obtaining the density and specific heat of copper from Table A-1, we see that (ρc)Cu ≈ 8900 kg/m ×
6 3 6 3
400 J/kg⋅K = 3.56 × 10 J/m ⋅K > (ρc)Al = 2.57 × 10 J/m ⋅K. Hence, for an equivalent sphere
diameter, the copper can store approximately 38% more thermal energy than the aluminum.
COMMENTS: Before the packed bed becomes fully charged, the temperature of the gas decreases
as it passes through the bed. Hence, the time required for a sphere to reach a prescribed state of
thermal energy storage increases with increasing distance from the bed inlet.
 PROBLEM 5.15
KNOWN: Thickness and properties of furnace wall. Thermal resistance of film on surface of wall
exposed to furnace gases. Initial wall temperature.
FIND: (a) Time required for surface of wall to reach a prescribed temperature, (b) Corresponding
value of film surface temperature.
SCHEMATIC:

ASSUMPTIONS: (1) Constant properties, (2) Negligible film thermal capacitance, (3) Negligible
radiation.
3
PROPERTIES: Carbon steel (given): ρ = 7850 kg/m , c = 430 J/kg⋅K, k = 60 W/m⋅K.
ANALYSIS: The overall coefficient for heat transfer from the surface of the steel to the gas is
−1 −1
−1 1   
+ 10−2 m2 ⋅ K/W 
1
U = ( R′′tot ) =  + R′′f  = = 20 W/m 2 ⋅ K.
h   25 W/m 2 ⋅ K 
Hence,
UL 20 W/m 2 ⋅ K × 0.01 m
Bi = = = 0.0033
k 60 W/m ⋅ K
and the lumped capacitance method can be used.
(a) It follows that
T − T∞
= exp ( − t/ τ t ) = exp ( −t/RC ) = exp ( −Ut/ ρ Lc )
Ti − T∞

ρ Lc T − T∞ 7850 kg/m3 ( 0.01 m ) 430 J/kg ⋅ K 1200 −1300

t=− ln =− ln
U Ti − T∞ 20 W/m 2 ⋅ K 300 − 1300

t = 3886s = 1.08h. <

(b) Performing an energy balance at the outer surface (s,o),
( ) (
h T∞ − Ts,o = Ts,o − Ts,i / R ′′f)
hT∞ + Ts,i / R′′f 25 W/m 2 ⋅ K ×1300 K + 1200 K/10-2 m2 ⋅ K/W
Ts,o = =
h + (1 / R′′f ) ( 25 +100 ) W/m2 ⋅ K
Ts,o = 1220 K. <
COMMENTS: The film increases τ t by increasing Rt but not Ct.
 PROBLEM 5.22
KNOWN: Droplet properties, diameter, velocity and initial and final temperatures.
FIND: Travel distance and rejected thermal energy.
SCHEMATIC:

ASSUMPTIONS: (1) Constant properties, (2) Negligible radiation from space.

3
PROPERTIES: Droplet (given): ρ = 885 kg/m , c = 1900 J/kg⋅K, k = 0.145 W/m⋅K, ε = 0.95.
ANALYSIS: To assess the suitability of applying the lumped capacitance method, use Equation 1.9
to obtain the maximum radiation coefficient, which corresponds to T = Ti.

h r = εσ Ti3 = 0.95 × 5.67 ×10 −8 W/m 2 ⋅ K 4 ( 500 K )3 = 6.73 W/m 2 ⋅ K.

Hence

Bi r =
h r ( ro / 3 )
=
( 6.73 W/m2 ⋅ K )(0.25 × 10−3 m/3)
= 0.0039
k 0.145 W/m ⋅ K
and the lumped capacitance method can be used. From Equation 5.19,
3
L ρc π D / 6
t= =
( )  1 − 1 
V 3ε π D2 σ
( )  T3
 f Ti3 

( 0.1 m/s ) 885 kg/m3 ( 1900 J/kg ⋅ K ) 0.5 × 10−3 m  1 1  1

L=  − 
18 × 0.95 × 5.67 × 10-8 W/m2 ⋅ K4  3003 5003  K3

L = 2.52 m. <
The amount of energy rejected by each droplet is equal to the change in its internal energy.

( 5 ×10−4 m )
3

E i − E f = ρ Vc ( Ti − Tf ) = 885 kg/m π
3 1900 J/kg ⋅ K ( 200 K )
6

E i − E f = 0.022 J. <
COMMENTS: Because some of the radiation emitted by a droplet will be intercepted by other
droplets in the stream, the foregoing analysis overestimates the amount of heat dissipated by radiation
to space.
 PROBLEM 5.24
KNOWN: Diameter and thermophysical properties of alumina particles. Convection conditions
associated with a two-step heating process.
FIND: (a) Time-in-flight (ti-f) required for complete melting, (b) Validity of assuming negligible
radiation.
SCHEMATIC:

ASSUMPTIONS: (1) Particle behaves as a lumped capacitance, (2) Negligible radiation, (3) Constant
properties.
ANALYSIS: (a) The two-step process involves (i) the time t1 to heat the particle to its melting point and
(ii) the time t2 required to achieve complete melting. Hence, ti-f = t1 + t2, where from Eq. (5.5),
ρ p Vc pθ ρ p Dpcp T −T
t1 = ln i = ln i ∞
hAs θ 6h Tmp − T∞

t1 =
( )
3970 kg m3 50 × 10−6 m 1560 J kg ⋅ K (300 − 10, 000 ) = 4 × 10−4 s
( )
ln
6 30, 000 W m 2 ⋅ K ( 2318 − 10, 000 )
Performing an energy balance for the second step, we obtain
t1 + t 2
∫t1 q conv dt = ∆Est

where qconv = hAs(T∞ - Tmp) and ∆Est = ρpVhsf. Hence,

t2 =
ρpDp h sf
=
(
3970 kg m3 50 × 10−6 m ) × 3.577 ×106 J kg = 5 × 10−4 s
6h (T∞ − Tmp ) (
6 30, 000 W m 2 ⋅ K ) (10, 000 − 2318 ) K

Hence t i − f = 9 × 10−4 s ≈ 1ms <

(b) Contrasting the smallest value of the convection heat flux, q ′′conv,min = h ( T∞ − Tmp ) = 2.3 × 108 W m 2

to the largest radiation flux, q′′rad,max = εσ Tmp

4
( 4
− Tsur )
= 6.5 × 105 W/m2, we conclude that radiation
is, in fact, negligible.

COMMENTS: (1) Since Bi = (hrp/3)/k ≈ 0.05, the lumped capacitance assumption is good. (2) In an
actual application, the droplet should impact the substrate in a superheated condition (T > Tmp), which
would require a slightly larger ti-f.
 PROBLEM 5.25
KNOWN: Diameters, initial temperature and thermophysical properties of WC and Co in composite
particle. Convection coefficient and freestream temperature of plasma gas. Melting point and latent
heat of fusion of Co.
FIND: Times required to reach melting and to achieve complete melting of Co.
SCHEMATIC:

ASSUMPTIONS: (1) Particle is isothermal at any instant, (2) Radiation exchange with surroundings
is negligible, (3) Negligible contact resistance at interface between WC and Co, (4) Constant
properties.
ANALYSIS: From Eq. (5.5), the time required to reach the melting point is
( ρ Vc )tot Ti − T∞
t1 = ln
h π Do2 Tmp − T∞
where the total heat capacity of the composite particle is

( ) / 6 300 J / kg ⋅ K
3
( ρ Vc )tot = ( ρ Vc )c + ( ρ Vc )s = 16, 000 kg / m3 π 1.6 × 10−5 m


( ) ( )   750 J / kg ⋅ K
 3 
3 
+8900 kg / m3 π / 6  2.0 × 10−5 m − 1.6 × 10−5 m
 

( )
= 1.03 × 10−8 + 1.36 ×10−8 J / K = 2.39 ×10−8 J / K

2.39 × 10−8 J / K (300 − 10, 000 ) K = 1.56 × 10−4 s

t1 = ln <
(1770 − 10, 000 ) K
( )( )
2
20, 000 W / m 2 ⋅ K π 2.0 × 10−5 m
The time required to melt the Co may be obtained by applying the first law, Eq. (1.11b) to a control
surface about the particle. It follows that
( )
Ein = hπ Do2 T∞ − Tmp t 2 = ∆Est = ρs (π / 6 ) Do3 − Di3 h sf ( )
( 2 ×10 m ) − (1.6 ×10 m )  2.59 ×10
 3 3
−5 −5
8900 kg / m
3
(π / 6 )  5
J / kg
−5
t2 = = 2.28 × 10 s <
(20, 000 W / m ⋅ K )π (2 × 10 m ) (10, 000 − 1770 ) K
2
2 −5

COMMENTS: (1) The largest value of the radiation coefficient corresponds to hr = εσ (Tmp + Tsur)
(T2 2
mp + Tsur . ) For the maximum possible value of ε = 1 and Tsur = 300K, hr = 378 W/m ⋅K << h =
2

2
20,000 W/m ⋅K. Hence, the assumption of negligible radiation exchange is excellent. (2) Despite the
large value of h, the small values of Do and Di and the large thermal conductivities (~ 40 W/m⋅K and
70 W/m⋅K for WC and Co, respectively) render the lumped capacitance approximation a good one.
(3) A detailed treatment of plasma heating of a composite powder particle is provided by Demetriou,
Lavine and Ghoniem (Proc. 5th ASME/JSME Joint Thermal Engineering Conf., March, 1999).
 PROBLEM 5.49
KNOWN: A long cylinder, initially at a uniform temperature, is suddenly quenched in a large oil bath.
FIND: (a) Time required for the surface to reach 500 K, (b) Effect of convection coefficient on surface
temperature history.
SCHEMATIC:

ASSUMPTIONS: (1) One-dimensional radial conduction, (2) Constant properties, (3) Fo > 0.2.
ANALYSIS: (a) Check first whether lumped capacitance method is applicable. For h = 50 W/m2⋅K,

hLc h ( ro 2 ) 50 W m ⋅ K ( 0.015 m / 2 )
2
Bic = = = = 0.221 .
k k 1.7 W m ⋅ K
Since Bic > 0.1, method is not suited. Using the approximate series solution for the infinite cylinder,

( ) ( ) ( )
θ * r* , Fo = C1 exp −ζ12 Fo × J 0 ζ1r* (1)

Solving for Fo and setting r* = 1, find

1  θ* 
Fo = − ln  
ζ12  C1J 0 (ζ1 ) 
T ( ro , t o ) − T∞ (500 − 350 ) K
where θ * = (1, Fo ) = = = 0.231 .
Ti − T∞ (1000 − 350 ) K
From Table 5.1, with Bi = 0.441, find ζ1 = 0.8882 rad and C1 = 1.1019. From Table B.4, find J0(ζ1) =
0.8121. Substituting numerical values into Eq. (2),
1
Fo = − ln [0.231 1.1019 × 0.8121] = 1.72 .
(0.8882 )2
From the definition of the Fourier number, Fo = α t ro2 , and α = k/ρc,

r2 ρc
t = Fo o = Fo ⋅ ro2
α k
t = 1.72 ( 0.015 m ) × 400 kg m3 × 1600 J kg ⋅ K 1.7 W m ⋅ K = 145s .
2
<
(b) Using the IHT Transient Conduction Model for a Cylinder, the following surface temperature
histories were obtained.
Continued...
 PROBLEM 5.49 (Cont.)

1000

Surface temperature, T(K)

900

800

700

600

500

400

300
0 50 100 150 200 250 300

Time, t(s)

h = 250 W/m^2.K
h = 50 W/m^2.K

Increasing the convection coefficient by a factor of 5 has a significant effect on the surface temperature,
greatly accelerating its approach to the oil temperature. However, even with h = 250 W/m2⋅K, Bi = 1.1
and the convection resistance remains significant. Hence, in the interest of accelerated cooling,
additional benefit could be achieved by further increasing the value of h.

COMMENTS: For Part (a), note that, since Fo = 1.72 > 0.2, the approximate series solution is
appropriate.
 PROBLEM 5.51
KNOWN: Sapphire rod, initially at a uniform temperature of 800K is suddenly cooled by a convection
process; after 35s, the rod is wrapped in insulation.
FIND: Temperature rod reaches after a long time following the insulation wrap.
SCHEMATIC:

ASSUMPTIONS: (1) One-dimensional radial conduction, (2) Constant properties, (3) No heat losses
from the rod when insulation is applied.
3
PROPERTIES: Table A-2, Aluminum oxide, sapphire (550K): ρ = 3970 kg/m , c = 1068 J/kg⋅K, k =
-5 2
22.3 W/m⋅K, α = 5.259×10 m /s.

ANALYSIS: First calculate the Biot number with Lc = ro/2,

h Lc h ( ro / 2 ) 1600 W/m ⋅ K ( 0.020 m/2)
2
Bi = = = = 0.72.
k k 22.3 W/m ⋅ K
Since Bi > 0.1, the rod cannot be approximated as a lumped capacitance system. The temperature
distribution during the cooling process, 0 ≤ t ≤ 35s, and for the time following the application of
insulation, t > 35s, will appear as

Eventually (t → ∞), the temperature of the rod will be uniform at T ( ∞ ) . To find T (∞ ) , write the
conservation of energy requirement for the rod on a time interval basis, E in − E out = ∆ E ≡ E final − E initial .

Using the nomenclature of Section 5.5.3 and basing energy relative to T∞, the energy balance becomes
−Q = ρ cV ( T ( ∞ ) − T∞ ) − Qo

where Qo = ρcV(Ti - T∞). Dividing through by Qo and solving for T (∞ ) , find

T ( ∞ ) = T∞ + ( Ti − T∞ )(1 − Q/Qo ) .
From the Groeber chart, Figure D.6, with
hr 1600 W/m 2 ⋅ K × 0.020m
Bi = o = = 1.43
k 22.3 W/m ⋅ K

( ) (
Bi 2Fo = Bi 2 α t/ro2 = (1.43 ) 2 5.259 × 10-6 m 2 /s × 35s/ ( 0.020m ) 2 = 0.95. )
find Q/Qo ≈ 0.57. Hence,
T ( ∞ ) = 300K + ( 800 − 300 ) K (1-0.57 ) = 515 K. <
COMMENTS: From use of Figures D.4 and D.5, find T(0,35s) = 525K and T(ro,35s) = 423K.
 PROBLEM 5.54
KNOWN: Diameter and initial temperature of roller bearings. Temperature of oil bath and
convection coefficient. Final centerline temperature. Number of bearings processed per hour.
FIND: Time required to reach centerline temperature. Cooling load.
SCHEMATIC:

ASSUMPTIONS: (1) One-dimensional, radial conduction in rod, (2) Constant properties.

( ) 3
PROPERTIES: Table A.1, St. St. 304 T = 548 K : ρ=7900 kg/m , k = 19.0 W/m⋅K, cp = 546
-6 2
J/kg⋅K, α = 4.40 × 10 m /s.

ANALYSIS: With Bi = h (ro/2)/k = 0.658, the lumped capacitance method can not be used. From
the one-term approximation of Eq. 5.49 c for the centerline temperature,
T −T
θ o∗ = o ∞ =
50 − 30
Ti − T∞ 500 − 30 ( )
= 0.0426 = C1 exp −ζ12 Fo = 1.1382 exp  − ( 0.9287 ) Fo 


2


where, for Bi = hro/k = 1.316, C1 = 1.1382 and ζ1 = 0.9287 from Table 5.1.

Fo = −n ( 0.0374 ) / 0.863 = 3.81

t f = Fo ro2 / α = 3.81( 0.05 m ) / 4.40 × 10−6 = 2162s = 36 min

2
<
From Eqs. 5.44 and 5.51, the energy extracted from a single rod is

 2θ ∗ 
Q = ρ cV ( Ti − T∞ ) 1 − o J1 (ζ1 )
 ζ1 

With J1 (0.9287) = 0.416 from Table B.4,

 0.0852 × 0.416 
Q = 7900 kg / m3 × 546 J / kg ⋅ K π ( 0.05m ) 1m  470 K 1 −
2
 = 1.53 × 107 J
   0.9287 
The nominal cooling load is

N Q 10 × 1.53 × 107 J
q= = = 70,800 W = 7.08 kW <
tf 2162s
COMMENTS: For a centerline temperature of 50°C, Eq. 5.49b yields a surface temperature of

T ( ro , t ) = T∞ + ( Ti − T∞ )θ o∗ Jo (ζ1 ) = 30°C + 470°C × 0.0426 × 0.795 = 45.9°C

 PROBLEM 5.58
KNOWN: A ball bearing is suddenly immersed in a molten salt bath; heat treatment to harden occurs
at locations with T > 1000K.
FIND: Time required to harden outer layer of 1mm.
SCHEMATIC:

ASSUMPTIONS: (1) One-dimensional radial conduction, (2) Constant properties, (3) Fo ≥ 0.2.
ANALYSIS: Since any location within the ball whose temperature exceeds 1000K will be hardened,
the problem is to find the time when the location r = 9mm reaches 1000K. Then a 1mm outer layer
will be hardened. Begin by finding the Biot number.
hr 5000 W/m2 ⋅ K ( 0.020m/2 )
Bi = o = = 1.00.
k 50 W/m ⋅ K
Using the one-term approximate solution for a sphere, find

( )
 
ln θ ∗ / C1 sin ζ1r∗
1 1
Fo = − .
2
ζ1  ζ1r ∗ 
From Table 5.1 with Bi = 1.00, for the sphere find ζ 1 = 15708
. rad and C1 = 1.2732. With r*
= r/ro = (9mm/10mm) = 0.9, substitute numerical values.
−1  (1000 − 1300 ) K 1 
Fo = ln  /1.2732 sin (1.5708 × 0.9 rad ) = 0.441.
(1.5708)2  (300 − 1300 ) K 1.5708 × 0.9 
From the definition of the Fourier number with α = k/ρc,
2
ro2 2 ρc  0.020m  kg J
t = Fo = Fo ⋅ ro = 0.441×   7800 3 × 500 / 50 W/m ⋅ K = 3.4s. <
α k  2  m kg ⋅ K
COMMENTS: (1) Note the very short time required to harden the ball. At this time it can be easily
shown the center temperature is T(0,3.4s) = 871 K.
-1
(2) The Heisler charts can also be used. From Fig. D.8, with Bi = 1.0 and r/ro = 0.9, read θ/θo =
0.69(±0.03). Since
θ = T − T∞ = 1000 − 1300 = −300K θ i = Ti − T∞ = −1000K
it follows that
θ θ θ θo θ θ
= 0.30. Since = ⋅ , then = 0.69 o
θi θi θo θi θi θi
and θ o / θ i = 0.30 / 0.69 = 0.43 ( ±0.02 ).
-1
From Fig. D.7 at θo/θi=0.43, Bi =1.0, read Fo = 0.45 (±0.03) and t = 3.5 (±0.2)s. Note the use of
tolerances associated with reading the charts to ±5%.
 PROBLEM 5.79
KNOWN: Initial temperature, density and specific heat of a material. Convection coefficient and
temperature of air flow. Time for embedded thermocouple to reach a prescribed temperature.
FIND: Thermal conductivity of material.
SCHEMATIC:

ASSUMPTIONS: (1) One-dimensional conduction in x, (2) Sample behaves as a semi-infinite

modium, (3) Constant properties.
ANALYSIS: The thermal response of the sample is given by Case 3, Eq. 5.60,

T ( x, t ) − Ti  x    hx h 2α t     x h α t 
= erfc   −  exp  +    erfc  + 
 k  2 αt k  
T∞ − Ti  2 α t    k 2    

where, for x = 0.01m at t = 300 s, [T(x,t) – Ti]/(T∞ - Ti) = 0.533. The foregoing equation must be
solved iteratively for k, with α = k/ρcp. The result is

k = 0.45 W / m ⋅ K <
-7 2
with α = 4.30 × 10 m /s.
COMMENTS: The solution may be effected by inserting the Transient Conduction/Semi-infinite
Solid/Surface Conduction Model of IHT into the work space and applying the IHT Solver. However,
the ability to obtain a converged solution depends strongly on the initial guesses for k and α.
 PROBLEM 5.82
KNOWN: Thickness, initial temperature and thermophysical properties of concrete firewall.
Incident radiant flux and duration of radiant heating. Maximum allowable surface temperatures at the
end of heating.
FIND: If maximum allowable temperatures are exceeded.
SCHEMATIC:

ASSUMPTIONS: (1) One-dimensional conduction in wall, (2) Validity of semi-infinite medium

approximation, (3) Negligible convection and radiative exchange with the surroundings at the
irradiated surface, (4) Negligible heat transfer from the back surface, (5) Constant properties.
ANALYSIS: The thermal response of the wall is described by Eq. (5.60)

2 q′′o (α t / π )  − x 2  q′′ x
1/ 2
 x 
T ( x, t ) = Ti + exp   − o erfc  
k  4α t  k  2 αt 
 

where, α = k / ρ c p = 6.92 × 10 −7 m 2 / s and for t = 30 min = 1800s, 2q′′o (α t / π )1/ 2 / k = 284.5 K. Hence,
at x = 0,

T ( 0,30 min ) = 25°C + 284.5°C = 309.5°C < 325°C <

( )
At x = 0.25m, − x 2 / 4α t = −12.54, q ′′o x / k = 1, 786K, and x / 2 (α t )1/ 2 = 3.54. Hence,

( )
T ( 0.25m, 30 min ) = 25°C + 284.5°C 3.58 × 10−6 − 1786°C × ( ~ 0 ) ≈ 25°C <
Both requirements are met.
COMMENTS: The foregoing analysis is conservative since heat transfer at the irradiated surface
due to convection and net radiation exchange with the environment have been neglected. If the
emissivity of the surface and the temperature of the surroundings are assumed to be ε = 1 and Tsur =
298K, radiation exchange at Ts = 309.5°C would be q′′rad = εσ Ts4 − Tsur
4
( )
= 6, 080 W / m 2 ⋅ K,
which is significant (~ 60% of the prescribed radiation).
 PROBLEM 5.84
KNOWN: Initial temperatures, properties, and thickness of two plates, each insulated on one
surface.
FIND: Temperature on insulated surface of one plate at a prescribed time after they are pressed
together.
SCHEMATIC:

ASSUMPTIONS: (1) One-dimensional conduction, (2) Constant properties, (3) Negligible

contact resistance.
3
PROPERTIES: Stainless steel (given): ρ = 8000 kg/m , c = 500 J/kg⋅K, k = 15 W/m⋅K.
ANALYSIS: At the instant that contact is made, the plates behave as semi-infinite slabs and, since
the (ρkc) product is the same for the two plates, Equation 5.63 yields a surface temperature of
Ts = 350K.

The interface will remain at this temperature, even after thermal effects penetrate to the insulated
surfaces. The transient response of the hot wall may therefore be calculated from Equations 5.40
and 5.41. At the insulated surface (x* = 0), Equation 5.41 yields
To − Ts
Ti − Ts (
= C1 exp −ζ12 Fo )
where, in principle, h → ∞ and T∞ → Ts . From Equation 5.39c, Bi → ∞ yields ζ 1 = 1.5707, and
from Equation 5.39b
4sinζ 1
C1 = = 1.273
2ζ1 + sin ( 2ζ1 )

αt 3.75 ×10 −6 m 2 / s ( 60s )

Also, Fo = = = 0.563.
L2 ( 0.02 m )2
Hence,
To − 350
400 − 350 (
= 1.273exp −1.5707 2 × 0.563 = 0.318 )
To = 365.9 K. <
COMMENTS: Since Fo > 0.2, the one-term approximation is appropriate.