Incroperas Principles of Heat and Mass Transfer 8e
Incroperas Principles of Heat and Mass Transfer 8e
Incroperas Principles of Heat and Mass Transfer 8e
Incropera’s Principles of Heat and Mass Transfer, 8/E (Pusan National University)
PROBLEM 1.1
PROBLEM 1.2
PROBLEM 1.3
PROBLEM 1.4
ASSUMPTIONS: (1) One-dimensional conduction in the x-direction, (2) Steady-state conditions, (3) Constant
properties.
ANALYSIS: From Fourier’s law, if q′′x and k are each constant it is evident that the gradient, dT dx = − q′′x k , is a
constant, and hence the temperature distribution is linear. The heat flux must be constant under one-dimensional,
steady-state conditions; and k is approximately constant if it depends only weakly on temperature. The heat flux and
heat rate when the outside wall temperature is T2 = -15°C are
q′′x =
−k
dT
= k
T1 − T2
=1W m ⋅ K
25 C − −15 C
=
(
133.3 W m 2 . (1)
)
dx L 0.30 m
2500
Heat loss, qx (W)
1500
500
-500
-1500
-20 -10 0 10 20 30 40
Ambient
Outside air temperature, T2 (C)
surface
For the concrete wall, k = 1 W/m⋅K, the heat loss varies linearly from +2667 W to -867 W and is zero when the inside
and outer surface temperatures are the same. The magnitude of the heat rate increases with increasing thermal
conductivity.
COMMENTS: Without steady-state conditions and constant k, the temperature distribution in a plane wall would not
be linear.
PROBLEM 1.5
PROBLEM 1.6
PROBLEM 1.7
PROBLEM 1.8
PROBLEM 1.9
PROBLEM 1.10
PROBLEM 1.11
PROBLEM 1.12
PROBLEM 1.13
PROBLEM 1.14
PROBLEM 1.15
KNOWN: Power required to maintain the surface temperature of a long, 25-mm diameter cylinder with an imbedded
electrical heater for different air velocities.
FIND: (a) Determine the convection coefficient for each of the air velocity conditions and display the results
graphically, and (b) Assuming that the convection coefficient depends upon air velocity as h = CVn, determine the
parameters C and n.
SCHEMATIC:
V(m/s) 1 2 4 8 12
Pe¢(W/m) 450 658 983 1507 1963
h (W/m2⋅K) 22.0 32.2 48.1 73.8 96.1
ASSUMPTIONS: (1) Temperature is uniform over the cylinder surface, (2) Negligible radiation exchange between
the cylinder surface and the surroundings, (3) Steady-state conditions.
ANALYSIS: (a) From an overall energy balance on the cylinder, the power dissipated by the electrical heater is
transferred by convection to the air stream. Using Newton’s law of cooling on a per unit length basis,
=Pe′ h (π D )( Ts − T∞ )
where Pe′ is the electrical power dissipated per unit length of the cylinder. For the V = 1 m/s condition, using the data
from the table above, find
=h 450 W m π × 0.025 m 300 − 40 = (
C 22.0 W m 2⋅K ) <
Repeating the calculations, find the convection coefficients for the remaining conditions which are tabulated above
and plotted below. Note that h is not linear with respect to the air velocity.
(b) To determine the (C,n) parameters, we plotted h vs. V on log-log coordinates. Choosing C = 22.12 W/m2⋅K(s/m)n,
assuring a match at V = 1, we can readily find the exponent n from the slope of the h vs. V curve. From the trials with
n = 0.8, 0.6 and 0.5, we recognize that n = 0.6 is a reasonable choice. Hence, C = 22.12 and n = 0.6. <
100
Coefficient, h (W/m^2.K)
100
Coefficient, h (W/m^2.K)
80
80 60
40
60
40 20
20
10
0 2 4 6 8 10 12 1 2 4 6 8 10
Air velocity, V (m/s) Air velocity, V (m/s)
PROBLEM 1.16
PROBLEM 1.17
PROBLEM 1.18
PROBLEM 1.19
KNOWN: Width, input power and efficiency of a transmission. Temperature and convection coefficient associated
with air flow over the casing.
FIND: Surface temperature of casing. Thermal convection resistance.
SCHEMATIC:
ASSUMPTIONS: (1) Steady state, (2) Uniform convection coefficient and surface temperature, (3) Negligible
radiation.
ANALYSIS: From Newton’s law of cooling,
q hAs ( Ts − T∞=
= ) 6 hW 2 ( Ts − T∞ )
where the output power is hPi and the heat rate is
Hence,
q 7833 W
Ts = T∞ + = 30°C + = 102.5°C <
6 hW 2 6 × 200 W / m 2 ⋅ K × ( 0.3m )
2
∆T / q x = ( Ts − T∞ ) / q x =
R t,conv = (102.5 − 30 ) K / 7833 W =
0.00926 K/W <
COMMENTS: (1) There will, in fact, be considerable variability of the local convection coefficient over the
transmission case and the prescribed value represents an average over the surface. (2) The convection thermal
resistance could equivalently be calculated from Rt,conv = 1/hA.
PROBLEM 1.20
PROBLEM 1.21
KNOWN: Length, diameter and calibration of a hot wire anemometer. Temperature of air stream.
Current, voltage drop and surface temperature of wire for a particular application.
FIND: Air velocity
SCHEMATIC:
ASSUMPTIONS: (1) Steady-state conditions, (2) Negligible heat transfer from the wire by natural
convection or radiation.
ANALYSIS: If all of the electric energy is transferred by convection to the air, the following equality must
be satisfied
= hA ( Ts − T∞ )
= EI
Pelec
π DL =
where A = 3.14 ×10−5 m 2 .
π ( 0.0005m × 0.02m ) =
Hence,
EI 5V × 0.1A
= = = 318 W/m 2 ⋅ K
( )
h
A ( Ts − T∞ ) 3.14 ×10−5m 2 50 C
( )
2
V = 6.25 ×10−5 h 2 = 6.25 ×10−5 318 W/m 2 ⋅ K = 6.3 m/s <
COMMENTS: The convection coefficient is sufficiently large to render buoyancy (natural convection) and
radiation effects negligible.
PROBLEM 1.22
KNOWN: Chip width and maximum allowable temperature. Coolant conditions.
FIND: Maximum allowable chip power for air and liquid coolants.
SCHEMATIC:
ASSUMPTIONS: (1) Steady-state conditions, (2) Negligible heat transfer from sides and bottom, (3) Chip
is at a uniform temperature (isothermal), (4) Negligible heat transfer by radiation in air.
ANALYSIS: All of the electrical power dissipated in the chip is transferred by convection to the coolant.
Hence,
P=q
and from Newton’s law of cooling,
2
P = hA(T - T∞) = h W (T - T∞).
In air,
2 2
Pmax = 250 W/m ⋅K(0.005 m) (85 - 15) ° C = 0.44 W. <
In the dielectric liquid
2 2
Pmax = 3000 W/m ⋅K(0.005 m) (85-15) ° C = 5.25 W. <
COMMENTS: Relative to liquids, air is a poor heat transfer fluid. Hence, in air the chip can dissipate far
less energy than in the dielectric liquid.
PROBLEM 1.23
PROBLEM 1.24
PROBLEM 1.25
PROBLEM 1.26
KNOWN: Air and wall temperatures of a room. Surface temperature, convection coefficient and
emissivity of a person in the room.
FIND: Basis for difference in comfort level between summer and winter. Ratio of thermal convection
resistance to thermal radiation resistance in summer and winter.
SCHEMATIC:
Air
“
qconv q“rad
Tsur = 27oC (summer)
o
= 20 C Tsur = -14oC (winter)
h = 2 W/m2-K
TS = 32oC
ε = 0.9
However, the heat flux due to radiation will differ, with values of
Summer: (4 4
) −8 2 4 4 4 4
q ′′rad =εs Ts − Tsur =0.9 × 5.67 × 10 W/m ⋅ K 305 − 300 K =28.3 W/m
2
( )
Winter: (4 4
) −8 2 4 4 4 4
q ′′rad =εs Ts − Tsur =0.9 × 5.67 × 10 W/m ⋅ K 305 − 287 K =95.4 W/m
2
( )
There is a significant difference between winter and summer radiation fluxes, and the chilled condition is
attributable to the effect of the colder walls on radiation. <
From Eq. 1.11, the thermal resistance due to convection is
R t,conv = ( Ts − T∞ ) / q′′conv A
∆T / q conv =
Continued…
R t,conv
=
( Ts − T∞ ) / q′′conv
R t,rad ( Ts − Tsur ) / q′′rad
=
Summer:
R t,conv
=
( 32 − 20 ) K / 24 W/m 2 2.83 <
R t,rad ( 32 − 27 ) K / 28.3 W/m 2
=
Winter:
R t,conv ( 32 − 20 ) K / 24 W/m 2 2.65
= <
R t,rad ( 32 − 14 ) K / 95.4 W/m 2
2
COMMENTS: (1) For a representative surface area of A = 1.5 m , the heat losses are qconv = 36 W,
qrad(summer) = 42.5 W and qrad(winter) = 143.1 W. The winter time radiation loss is significant and if
maintained over a 24 h period would amount to 2,950 kcal. (2) The convection resistance is larger than the
radiation resistance but they are the same order of magnitude. There isn’t much difference in the resistances
between summer and winter conditions; the main difference is the larger temperature difference through
which radiation occurs in the winter as compared to summer.
PROBLEM 1.27
PROBLEM 1.28
KNOWN: Spherical shaped instrumentation package with prescribed surface emissivity within a large space-
simulation chamber having walls at 77 K.
FIND: Acceptable power dissipation for operating the package surface temperature in the range Ts = 40 to 85°C.
Show graphically the effect of emissivity variations for 0.2 and 0.3.
SCHEMATIC:
ASSUMPTIONS: (1) Uniform surface temperature, (2) Chamber walls are large compared to the spherical package,
and (3) Steady-state conditions.
ANALYSIS: From an overall energy balance on the package, the internal power dissipation Pe will be transferred by
radiation exchange between the package and the chamber walls. From Eq. 1.7,
q rad = Pe = εAs σ Ts4 - Tsur
4
( )
For the condition when Ts = 40°C, with As = πD2 the power dissipation will be
( )
Pe = 0.25 π × 0.102 m 2 × 5.67 ×10-8 W m 2 ⋅ K 4 × ( 40 + 273) - 77 4 K 4 = 4.3 W
4
<
Repeating this calculation for the range 40 ≤ Ts ≤ 85°C, we can obtain the power dissipation as a function of surface
temperature for the ε = 0.25 condition. Similarly, with 0.2 or 0.3, the family of curves shown below has been
obtained.
10
Power dissipation, Pe (W)
2
40 50 60 70 80 90
COMMENTS: (1) As expected, the internal power dissipation increases with increasing emissivity and surface
temperature. Because the radiation rate equation is non-linear with respect to temperature, the power dissipation will
likewise not be linear with surface temperature.
(2) What is the maximum power dissipation that is possible if the surface temperature is not to exceed 85°C? What
kind of a coating should be applied to the instrument package in order to approach this limiting condition?
PROBLEM 1.29
PROBLEM 1.30
PROBLEM 1.31
PROBLEM 1.32
PROBLEM 1.33
PROBLEM 1.34
PROBLEM 1.35
KNOWN: Hot and cold reservoir temperatures of an internally reversible refrigerator. Thermal resistances between
refrigerator and hot and cold reservoirs.
FIND: Expressions for modified Coefficient of Performance and power input of refrigerator.
SCHEMATIC:
qout
∙ Rt,h
Rt,c
qin
ANALYSIS: Heat is transferred from the low temperature reservoir (the refrigerated space) at Tc to the refrigerator
unit, through the resistance Rt,c, with Tc > Tc,i. Heat is rejected from the refrigerator unit to the higher temperature
reservoir (the surroundings), through the resistance Rt,h, with Th,i > Th. The heat input and output rates can be
expressed in a manner analogous to Equations 1.18a and 1.18b.
q=
in (Tc − Tc ,i ) / Rt ,c (1)
q=
out (Th ,i − Th ) / Rt ,h (2)
Equations (1) and (2) can be solved for the internal temperatures, to yield
1 + COPm
Th ,i =
Th + qout Rt ,h =
Th + qin Rt ,h (3)
COPm
Tc ,=
i Tc − qin Rt ,c (4)
1 + COPm
qout = qin (5)
COPm
using the definition of COPm given in the problem statement. The modified Coefficient of Performance can then be
expressed as
Continued…
Tc ,i Tc − qin Rt ,c
=
COPm =
Th ,i − Tc ,i 1 + COPm
Th + qin Rt ,h − Tc + qin Rt ,c
COPm
Tc − qin Rtot
COPm = <
Th − Tc + qin Rtot
qin T −T + q R
=W = qin h c in tot <
COPm Tc − qin Rtot
COMMENTS: As qin or Rtot goes to zero, the Coefficient of Performance approaches the maximum Carnot value,
COPm = COPC = Tc /(Th - Tc).
PROBLEM 1.36
KNOWN: Hot and cold reservoir temperatures of an internally reversible refrigerator. Thermal resistances between
refrigerator and hot and cold reservoirs under clean and dusty conditions. Desired cooling rate.
FIND: Modified Coefficient of Performance and power input of refrigerator under clean and dusty conditions.
SCHEMATIC:
= 25°C
qout
qin = 750 W
qin
= 5°C
ASSUMPTIONS: (1) Refrigerator is internally reversible, (2) Steady-state operation, (3) Cold side thermal resistance
does not degrade over time.
ANALYSIS: According to Problem 1.35, the modified Coefficient of Performance and power input are given by
Tc − qin Rtot
COPm = (1)
Th − Tc + qin Rtot
T −T + q R
W = qin h c in tot (2)
Tc − qin Rtot
Under new, clean conditions, with Rtot,n = Rh,n + Rc,n = 0.09 K/W, we find
Under dusty, conditions, with Rtot,d = Rh,d + Rc,n = 0.15 K/W, we find
Continued...
COMMENTS: (1) The cooling rates and power input values are time-averaged quantities. Since the refrigerator
does not run constantly, the instantaneous power requirements would be higher than calculated. (2) In practice, when
the condenser coils become dusty the power input does not adjust to maintain the cooling rate. Rather, the
refrigerator’s duty cycle would increase. (3) The ideal Carnot Coefficient of Performance is COPC = Tc/(Th – Tc) = 14
and the corresponding power input is 54 W. (4) This refrigerator’s energy efficiency is poor. Less power would be
consumed by more thoroughly insulating the refrigerator, and designing the refrigerator to minimize heat gain upon
opening its door, in order to reduce the cooling rate, qin.
PROBLEM 1.37
KNOWN: Width, surface emissivity and maximum allowable temperature of an electronic chip. Temperature of air
and surroundings. Convection coefficient.
2 1/4
FIND: (a) Maximum power dissipation for free convection with h(W/m ⋅K) = 4.2(T - T∞) , (b) Maximum power
2
dissipation for forced convection with h = 250 W/m ⋅K.
SCHEMATIC:
ASSUMPTIONS: (1) Steady-state conditions, (2) Radiation exchange between a small surface and a large enclosure,
(3) Negligible heat transfer from sides of chip or from back of chip by conduction through the substrate.
ANALYSIS: Subject to the foregoing assumptions, electric power dissipation by the chip must be balanced by
convection and radiation heat transfer from the chip. Hence, from Eq. (1.10),
where
(a) If heat transfer is by natural convection,
<
(b) If heat transfer is by forced convection,
<
COMMENTS: Clearly, radiation and natural convection are inefficient mechanisms for transferring heat from the
2
chip. For Ts = 85°C and T∞ = 25°C, the natural convection coefficient is 11.7 W/m ⋅K. Even for forced convection
2
with h = 250 W/m ⋅K, the power dissipation is well below that associated with many of today’s processors. To
provide acceptable cooling, it is often necessary to attach the chip to a highly conducting substrate and to thereby
provide an additional heat transfer mechanism due to conduction from the back surface.
PROBLEM 1.38
KNOWN: Width, input power, and efficiency of a transmission. Temperature and convection coefficient for air flow
over the casing. Emissivity of casing and temperature of surroundings.
FIND: Surface temperature of casing. Resistances due to convection and radiation.
SCHEMATIC:
ASSUMPTIONS: (1) Steady state, (2) Uniform convection coefficient and surface temperature, (3) Radiation
exchange with large surroundings.
ANALYSIS: Heat transfer from the case must balance heat dissipation in the transmission, which may be expressed
as q = Pi – Po = Pi (1 - h) = 150 hp × 746 W/hp × 0.07 = 7833 W. Heat transfer from the case is by convection and
radiation, in which case
=
(
q As h ( Ts − T∞ ) + εs Ts4 − Tsur
4
)
2
where As = 6 W . Hence,
7833 W
2
( )
6 ( 0.30 m ) 200 W / m 2 ⋅ K ( Ts − 303K ) + 0.8 × 5.67 × 10−8 W / m 2 ⋅ K 4 Ts4 − 3034 K 4
Ts ≈ 373 K =
100°C <
The thermal resistances can be found from Eq. 1.11, that is, Rt = ∆T/q. For convection, the relevant temperature
difference is Ts − T∞ .
Continued …
COMMENTS: (1) For Ts ≈ 373 K, qconv ≈ 7,560 W and qrad ≈ 270 W, in which case heat transfer is dominated by
convection. This can also be seen by the fact that the resistance to radiation heat transfer is much larger than the
resistance to convection heat transfer. (2) If radiation is neglected, the corresponding surface temperature is Ts =
102.5°C.
PROBLEM 1.39
PROBLEM 1.40
KNOWN: Surface areas, convection heat transfer coefficient, surface emissivity of gear box and generator.
Temperature of nacelle. Electric power generated by the wind turbine and generator efficiency.
qconv,gb qconv,gen
qrad,gb qrad,gen
ASSUMPTIONS: (1) Steady-state conditions, (2) Interior of nacelle can be treated as large surroundings, (3)
Negligible heat transfer between the gear box and the generator.
ANALYSIS: Heat is generated within both the gear box and the generator. The mechanical work into the generator
can be determined from the electrical power, P = 2.5 × 106 W, and the efficiency of the generator as
Wgen,in ==
P / hgen 2.5 × 106 W / 0.95 =
2.63 × 106 W
The generator surface temperature may be found by using a numerical solver, or by trial-and-error, yielding
Heat is also generated by the gear box. The heat generated in the gear box may be determined from knowledge of the
heat generated cumulatively by the gear box and the generator, which is provided in Example 3.1 and is q = qgen + qgb
= 0.33 × 106 W. Hence, qgb = q - qgen = 0.33 × 106W - 0.13 × 106W = 0.20 × 106W and
PROBLEM 1.41
KNOWN: Radial distribution of heat dissipation in a cylindrical container of radioactive wastes. Surface
convection conditions.
FIND: Total energy generation rate and surface temperature.
SCHEMATIC:
ASSUMPTIONS: (1) Steady-state conditions, (2) Negligible temperature drop across thin container wall.
ANALYSIS: The rate of energy generation is
o r
E g = qdV=q ∫ 1- ( r/ro ) 2π rLdr
∫
2
o
0
= (
E g 2π Lq o ro2 / 2 − ro2 / 4 )
or per unit length,
πq r 2
E ′g = o o . <
2
Performing an energy balance for a control surface about the container yields, at an instant,
E ′g − E ′out = 0
and substituting for the convection heat rate per unit length,
π q o ro2
= h ( 2π ro )( Ts − T∞ )
2
q o ro
Ts = T∞ + . <
4h
COMMENTS: The temperature within the radioactive wastes increases with decreasing r from Ts at ro to a
maximum value at the centerline.
PROBLEM 1.42
KNOWN: Thickness and initial temperature of an aluminum plate whose thermal environment is changed.
FIND: (a) Initial rate of temperature change, (b) Steady-state temperature of plate.
SCHEMATIC:
ASSUMPTIONS: (1) Negligible end effects, (2) Uniform plate temperature at any instant, (3) Constant properties,
(4) Adiabatic bottom surface, (5) Negligible radiation from surroundings, (6) No internal heat generation.
ANALYSIS: (a) Applying an energy balance, Eq. 1.12c, at an instant of time to a control volume about the plate,
E in − E out =
E st , it follows for a unit surface area.
( ) ( ) ( )
αSGS 1m 2 − E 1m 2 − q′′conv 1m 2 =( d dt )( McT ) =ρ 1m 2 × L c ( dT dt ) . ( )
Rearranging and substituting from Eqs. 1.3 and 1.5, we obtain
=
dT dt (1 ρ Lc ) αSGS − εs Ti4 − h ( Ti − T∞ ) .
αSGS = εs T 4 + h ( T − T∞ ) (2)
PROBLEM 1.43
KNOWN: Blood inlet and outlet temperatures and flow rate. Dimensions of tubing.
FIND: Required rate of heat addition and estimate of kinetic and potential energy changes.
SCHEMATIC:
ASSUMPTIONS: (1) Steady-state conditions, (2) Incompressible liquid with negligible kinetic and potential energy
changes, (3) Blood has properties of water.
PROPERTIES: Table A.6, Water ( T ≈ 300 K): cp,f = 4179 J/ kg∙ K, ρf = 1/vf = 997 kg/m3.
p (Tout - Tin )
q = mc
where
m = 997 kg/m3 × 200 m/min × 10-6 m3 /m
= ρf ∀ 60 s/min = 3.32 × 10-3 kg/s
Thus
q = 3.32 × 10-3 kg/s × 4179 J/kg ⋅ K × (37o C - 10o C) = 375 W <
The velocity in the tube is given by
/A = 200 m/min × 10-6 m3 /m
V=∀ (60 s/min × 6.4 × 10-3 m × 1.6 × 10-3 m) = 0.33 m/s
c
PROBLEM 1.44
PROBLEM 1.45
PROBLEM 1.46
PROBLEM 1.47
PROBLEM 1.48
KNOWN: Silicon wafer positioned in furnace with top and bottom surfaces exposed to hot and cool zones,
respectively.
FIND: (a) Initial rate of change of the wafer temperature corresponding to the wafer temperature Tw,i = 300 K, and
(b) Steady-state temperature reached if the wafer remains in this position. How significant is convection for this
situation? Sketch how you’d expect the wafer temperature to vary as a function of vertical distance.
SCHEMATIC:
ASSUMPTIONS: (1) Wafer temperature is uniform, (2) Transient conditions when wafer is initially positioned, (3)
Hot and cool zones have uniform temperatures, (3) Radiation exchange is between small surface (wafer) and large
enclosure (chamber, hot or cold zone), and (4) Negligible heat losses from wafer to mounting pin holder.
ANALYSIS: The energy balance on the wafer illustrated in the schematic above includes convection from the upper
(u) and lower (l) surfaces with the ambient gas, radiation exchange with the hot- and cool-zone (chamber)
surroundings, and the rate of energy storage term for the transient condition.
E ′′in − E ′′out =
E ′′st
dT
q′′rad,h + q′′rad,c − q′′conv,u − q′′conv,l =
ρ cd w
dt
(
εs Tsur,h
4
)
− Tw4 + εs Tsur,c
4
( )
− Tw4 − h u ( Tw − T∞ ) − h l ( Tw − T∞ ) =
ρ cd
d Tw
dt
(a) For the initial condition, the time rate of temperature change of the wafer is determined using the energy balance
above with = Tw T= w,i 300 K,
( ) (
0.65 × 5.67 × 10−8 W / m 2 ⋅ K 4 1500 4 − 3004 K 4 + 0.65 × 5.67 × 10−8 W / m 2 ⋅ K 4 3304 − 3004 K 4 )
−8 W / m 2 ⋅ K ( 300 − 700 ) K − 4 W / m 2 ⋅ K ( 300 − 700 ) K =
( d Tw / dt )i = 104 K / s <
(b) For the steady-state condition, the energy storage term is zero, and the energy balance can be solved for the steady-
state wafer temperature, Tw = Tw,ss .
Continued …..
(
0.65s 15004 − Tw,ss
4
) (
K 4 + 0.65s 3304 − Tw,ss
4
K4 )
−8 W / m 2 ⋅ K ( Tw,ss − 700 ) K − 4 W / m 2 ⋅ K ( Tw,ss − 700 ) K =
0
If the wafer were elevated above the present operating position, its temperature would increase, since the lower
surface would begin to experience radiant exchange with progressively more of the hot zone chamber. Conversely, by
lowering the wafer, the upper surface would experience less radiant exchange with the hot zone chamber, and its
temperature would decrease. The temperature-distance trend might appear as shown in the sketch.
PROBLEM 1.49
PROBLEM 1.50
PROBLEM 1.51
PROBLEM 1.52
KNOWN: Inner surface heating and new environmental conditions associated with a spherical shell of prescribed
dimensions and material.
FIND: (a) Governing equation for variation of wall temperature with time. Initial rate of temperature change, (b)
Steady-state wall temperature, (c) Effect of convection coefficient on canister temperature.
SCHEMATIC:
ASSUMPTIONS: (1) Negligible temperature gradients in wall, (2) Constant properties, (3) Uniform, time-
independent heat flux at inner surface.
PROPERTIES: Table A.1, Stainless Steel, AISI 302: ρ = 8055 kg/m3, c p = 535 J/kg⋅K.
− E E st . Identifying
ANALYSIS: (a) Performing an energy balance on the shell at an instant of time, E in out =
relevant processes and solving for dT/dt,
( ) ( )
q′′i 4π ri2 − h 4π ro2 ( T − T=
4
(
∞ ) ρ π ro − ri c p
3
3
)
3 dT
dt
dT 3 q′′i ri2 − hro2 ( T − T∞ ) .
=
p (
dt ρ c r 3 − r 3 o i )
Continued …..
(c) Parametric calculations were performed using the IHT First Law Model for an Isothermal Hollow Sphere. As
shown below, there is a sharp increase in temperature with decreasing values of h < 1000 W/m2⋅K. For T > 380 K,
boiling will occur at the canister surface, and for T > 410 K a condition known as film boiling (Chapter 10) will occur.
The condition corresponds to a precipitous reduction in h and increase in T.
1000
900
800
Temperature, T(K)
700
600
500
400
300
100 400 800 2000 6000 10000
Although the canister remains well below the melting point of stainless steel for h = 100 W/m2⋅K, boiling should be
avoided, in which case the convection coefficient should be maintained at h > 1000 W/m2⋅K.
COMMENTS: The governing equation of part (a) is a first order, nonhomogenous differential equation with constant
coefficients. Its solution is θ= (S/R ) (1 − e− Rt ) + θi e− Rt , where θ ≡ T − T∞ , S ≡ 3q′′i ri2 / ρ cp ( ro3 − ri3 ) ,
R ( )
3hro2 /ρ c p ro3 − ri3 . Note results for t → ∞ and for S = 0.
PROBLEM 1.53
PROBLEM 1.54
PROBLEM 1.55
PROBLEM 1.56
PROBLEM 1.57
KNOWN: Daily thermal energy generation, surface area, temperature of the environment, and heat transfer
coefficient.
FIND: (a) Skin temperature when the temperature of the environment is 20ºC, and (b) Rate of perspiration to maintain
skin temperature of 33ºC when the temperature of the environment is 33ºC.
SCHEMATIC:
T∞
h = 3 W/m2⋅K E out
Air
E g
ASSUMPTIONS: (1) Steady-state conditions, (2) Thermal energy is generated at a constant rate throughout the day,
(3) Air and surrounding walls are at same temperature, (4) Skin temperature is uniform, (5) Bathing suit has no effect
on heat loss from body, (6) Heat loss is by convection and radiation to the environment, and by perspiration in Part 2.
(Heat loss due to respiration, excretion of waste, etc., is negligible.), (7) Large surroundings.
PROPERTIES: Table A.11, skin: ε = 0.95, Table A.6, water (306 K): ρ = 994 kg/m3, hfg = 2421 kJ/kg.
ANALYSIS:
(a) The rate of energy generation is:
E g =
2000 × 103 cal/day (0.239 cal/J × 86,400 s/day ) = 96.9 W
= q = 96.9 W. Energy outflow is due to convection and net radiation from the surface to the environment,
Thus E out
Equations 1.3a and 1.7, respectively.
(b) If the skin temperature is 33ºC when the temperature of the environment is 33ºC, there will be no heat loss due to
convection or radiation. Thus, all the energy generated must be removed due to perspiration:
E out = mh
fg
We find:
This is the perspiration rate in mass per unit time. The volumetric rate is:
=
∀ = 4.0 × 10−5 kg/s / 994 kg/m3 = 4.0 × 10−8 m3 /s = 4.0 × 10−5 /s
m/ρ <
COMMENTS: (1) In Part 1, heat losses due to convection and radiation are 32.4 W and
60.4 W, respectively. Thus, it would not have been reasonable to neglect radiation. Care must be taken to include
radiation when the heat transfer coefficient is small, even if the problem statement does not give any indication of its
importance. (2) The rate of thermal energy generation is not constant throughout the day; it adjusts to maintain a
constant core temperature. Thus, the energy generation rate may decrease when the temperature of the environment
goes up, or increase (for example, by shivering) when the temperature of the environment is low. (3) The skin
temperature is not uniform over the entire body. For example, the extremities are usually cooler. Skin temperature
also adjusts in response to changes in the environment. As the temperature of the environment increases, more blood
flow will be directed near the surface of the skin to increase its temperature, thereby increasing heat loss. (4) If the
perspiration rate found in Part 2 was maintained for eight hours, the person would lose 1.2 liters of liquid. This
demonstrates the importance of consuming sufficient amounts of liquid in warm weather.
PROBLEM 1.58
KNOWN: Thermal conductivity, thickness and temperature difference across a sheet of rigid extruded insulation.
Cold wall temperature, surroundings temperature, ambient temperature and emissivity.
FIND: (a) The value of the convection heat transfer coefficient on the cold wall side in units of W/m2⋅°C or W/m2⋅K,
and, (b) The cold wall surface temperature for emissivities over the range 0.05 ≤ ε ≤ 0.95 for a hot wall temperature of
T1 = 30 °C.
SCHEMATIC:
T2 = 20°C = 293 K
TT11 ==30°C
32°C
L=25 mm
ASSUMPTIONS: (1) One-dimensional conduction in the x-direction, (2) Steady-state conditions, (c) Constant
properties, (4) Large surroundings.
ANALYSIS:
(a) An energy balance on the control surface shown in the schematic yields
E in = E out or q cond = q conv + q rad
Substituting from Fourier’s law, Newton’s law of cooling, and Eq. 1.7 yields
T1 - T2
k = h(T2 - T∞ ) + εσ(T24 - Tsur
4
) (1)
L
1 T -T
or h= [k 1 2 - εσ(T24 - Tsur
4
)]
(T2 - T∞ ) L
Substituting values,
1 W (30 - 20)o C W
h= [0.029 × - 0.95 × 5.67 × 10-8 2 4 (2934 - 3204 ) K 4 ]
(20 - 5)o C m⋅K 0.02 m m ⋅K
W
h = 12.2 <
m2 ⋅ K
W
Since temperature differences of 1 K and 1°C are equal, h = 12.2 . <
m ⋅°C
2
Continued...
(b) Equation (1) may be solved iteratively to find T2 for any emissivity ε. IHT was used for this purpose, yielding the
following.
295
290
285
280
0 0.2 0.4 0.6 0.8 1
Emissivity
COMMENTS: (1) Note that as the wall emissivity increases, the surface temperature increases since the surroundings
temperature is relatively hot. (2) The IHT code used in part (b) is shown below. (3) It is a good habit to work in
temperature units of kelvins when radiation heat transfer is included in the solution of the problem.
//Problem 1.56
//Equation (1) is
PROBLEM 1.59
PROBLEM 1.60
PROBLEM 1.61
PROBLEM 1.62
PROBLEM 1.63