PROBEM.3 KNOWN: Relation for the Rayleigh number. FIND: Rayleigh number for four fluids for prescribed conditions. ASSUMPTIONS: (1 Perfect gas behavior for specified gases. PROPERTIES: Table A-4, Air (400K, 1 atm: ν =.41 - m /s, α = 38.3 - m /s, β = 1/T = 1/400K =.50-3 K -1 ; Table A-4, Helium (400K, 1 atm: ν = 1 - m /s, α = 5 - m /s, β = 1/T =.50-3 K -1 ; Table A-5, Glycerin (1 C = 85K: ν = 830 - m /s, α = 0.4-7 m /s, β = 0.475-3 K -1 ; Table A-, Water (37 C = 3K, sat. liq.: ν = μ f v f = 5 - N s/m 1.007-3 m 3 /kg = 0.700 - m /s, α = k f v f /c p,f = 0.8 W/m K 1.007-3 m 3 /kg/4178 J/kg K = 0.151 - m /s, β f = 31. - K -1. ANAYSIS: The Rayleigh number, a dimensionless parameter used in free convection analysis, is defined as the product of the Grashof and Prandtl numbers. gβδt 3 μcp gβδt 3 ( νρ cp gβδt 3 Ra Gr Pr = = = ν k ν k να where α = k/ρc p and ν = μ/ρ. The numerical values for the four fluids follow: Air (400K, 1 atm Ra =.8m/s ( 1/400K 30K ( 0.01m 3 /.41 m / s 38.3 m / s = 77 Helium (400K, 1 atm Ra =.8m/s ( 1/ 400K 30K ( 0.01m 3 /1 m / s 5 m / s = 1.5 Glycerin (85K 3 1 3 7 Ra =.8m/s 0.475 K 30K 0.01m / 830 m / s 0.4 m / s = 51 ( ( Water (3K ( ( 3 1 3 Ra =.8m/s 0.3 K 30K 0.01m / 0.700 m / s 0.151 m / s = 1.01. COMMENTS: (1 Note the wide variation in values of Ra for the four fluids. A large value of Ra implies enhanced free convection, however, other properties affect the value of the heat transfer coefficient.
PROBEM.7 KNOWN: arge vertical plate with uniform surface temperature of 130 C suspended in quiescent air at 5 C and atmospheric pressure. FIND: (a Boundary layer thickness at 0.5 m from lower edge, (b Maximum velocity in boundary layer at this location and position of maximum, (c Heat transfer coefficient at this location, (d ocation where boundary layer becomes turbulent. ASSUMPTIONS: (1 Isothermal, vertical surface in an extensive, quiescent medium, ( Boundary layer assumptions valid. PROPERTIES: Table A-4, Air ( 0.030 W/m K, Pr = 0.700. ( = + = Tf Ts T / 350K, 1 atm : ν = 0. - m /s, k = ANAYSIS: (a From the similarity solution results, Fig..4 (see above right, the boundary layer thickness corresponds to a value of η 5. From Eqs..13 and.1, y= ηx( Gr x /4 1/4 (1 m 1 ( ( ( ν 3 3 3 Grx = gβ Ts T x / =.8 130 5 K x / 0. m / s.718 x s = 350K ( 3 1/4 y 5( 0.5m.718 ( 0.5 / 4 1.74 = m = 17.5 mm. (b From the similarity solution shown above, the maximum velocity occurs at η 1 with f η = 0.75. From Eq..15, find ( (3 ν 1/ 1/ 0. m / s 3 u = Gr x f ( η =.718 ( 0.5 0.75 = 0.47 m/s. x 0.5m The maximum velocity occurs at a value of η = 1; using Eq. (3, it follows that this corresponds to a position in the boundary layer given as ymax = 1/ 5 ( 17.5 mm = 3.5 mm. (c From Eq..1, the local heat transfer coefficient at x = 0.5 m is 1/4 3 1/4 x = hxx/k = ( Gr x / 4 g( Pr =.718 ( 0.5 / 4 0.50 = 35.7 h k/x 35.7 0.030 W/m K/0.5 m 4.3 W/m x = x = = K. The value for g(pr is determined from Eq..0 with Pr = 0.700. (d According to Eq..3, the boundary layer becomes turbulent at x c given as Ra Gr Pr x /.718 x,c = x,c c 0.700 = 0.0 m. COMMENTS: Note that β = 1/T f is a suitable approximation for air. ( 1/3
PROBEM.14 KNOWN: Aluminum plate (alloy 04 at an initial uniform temperature of 7 C is suspended in a room where the ambient air and surroundings are at 7 C. FIND: (a Expression for time rate of change of the plate, (b Initial rate of cooling (K/s when plate temperature is 7 C, (c Validity of assuming a uniform plate temperature, (d Decay of plate temperature and the convection and radiation rates during cooldown. ASSUMPTIONS: (1 Plate temperature is uniform, ( Ambient air is quiescent and extensive, (3 Surroundings are large compared to plate. PROPERTIES: Table A.1, Aluminum alloy 04 (T = 500 K: ρ = 770 kg/m 3, k = 18 W/m K, c = 83 J/kg K; Table A.4, Air (T f = 400 K, 1 atm: ν =.41 - m /s, k = 0.0388 W/m K, α = 38.3 - m /s, Pr = 0.0. ANAYSIS: (a From an energy balance on the plate with free convection and radiation exchange, E & out = E & st, we obtain 4 4 4 4 h As ( Ts T As ( Ts Tsur A dt dt stc or h ( Ts T ( Ts Tsur ε σ ρ = = + dt dt ρtc εσ where T s, the plate temperature, is assumed to be uniform at any time. (b To evaluate (dt/dt, estimate h. First, find the Rayleigh number, ( ( ( 3 3.8m s 1 400 K 7 7 K 0.3m Ra 8 = gβ( Ts T να = = 1.308..41 m s 38.3 m s Eq..7 is appropriate; substituting numerical values, find ( 8 1/4 0.70Ra 1/4 0.70 1.308 = 0.8 + = 0.8 + = 55.5 4/ 4/ /1 /1 1+ ( 0.4 Pr 1 + ( 0.4 0.0 h = k = 55.5 0.0338 W m K 0.3m =.5 W m K Continued...
PROBEM.14 (Cont. dt = dt 770 kg m 3 0.015 m 83J kg K 8 4 4 4 4.5 W m K 7 7 K + 0.5 5.7 W m K 500 300 K = 0.0 K s. ( ( ( (c The uniform temperature assumption is justified if the Biot number criterion is satisfied. With c (V/A s = (A s t/a s = (t/ and htot = hconv + hrad, Bi = htot ( t k 0.1. Using the linearized radiation coefficient relation, find ( ( ( ( ( 8 4 3 hrad = εσ Ts + Tsur Ts + Tsur = 0.5 5.7 W m K 500 + 300 500 + 300 K = 3.8 W m K Hence, Bi = (.5 + 3.8 W/m K(0.015 m//18 W/m K = 4.07-4. Since Bi 0.1, the assumption is appropriate. (d The temperature history of the plate was computed by combining the umped Capacitance Model of IHT with the appropriate Correlations and Properties Toolpads. 30 300 50 Temperature, T(C 150 1 70 30 0 5000 000 15000 Time, t(s Heat rate, q(w 00 150 0 50 0 0 3000 000 000 00 15000 Time, t(s Convection heat rate, qconv(w Radiation heat rate, qrad(w Due to the small values of h and h rad, the plate cools slowly and does not reach 30 C until t 14000s = 3.8h. The convection and radiation rates decrease rapidly with increasing t (decreasing T, thereby decelerating the cooling process. COMMENTS: The reduction in the convection rate with increasing time is due to a reduction in the thermal conductivity of air, as well as the values of h and T.
PROBEM.8 KNOWN: Vertical panel with uniform heat flux exposed to ambient air. FIND: Allowable heat flux if maximum temperature is not to exceed a specified value, T max. ASSUMPTIONS: (1 Constant properties, ( Radiative exchange with surroundings negligible. PROPERTIES: Table A-4, Air (T f = (T / + T / = (35.4 + 5 C/ = 30. C = 303K, 1 atm: ν = 1.1 - m /s, k =.5-3 W/m K, α =. - m /s, Pr = 0.707. ANAYSIS: Following the treatment of Section..1 for a vertical plate with uniform heat flux (constant q s, the heat flux can be evaluated as q s = h Δ T/ where Δ T/ = Ts( / T (1, and h is evaluated using an appropriate correlation for a constant temperature vertical plate. From Eq..8, ( 1/5 ΔTx Tx T = 1.15 x / Δ T/ (3 and recognizing that the maximum temperature will occur at the top edge, x =, use Eq. (3 to find o Δ T/ = ( 37 5 C /1.15( 1/1 1/5 =.4 o C or T/ = 35.4 o C. Calculate now the Rayleigh number based upon ΔT /, with T f = (T / + T / = 303K, gβδt 3 Ra = where Δ T =Δ T/ (4 να Ra =.8m / s ( 1/ 303K.4K ( 1m 3 /1.1 m / s. m / s =.07 8. Since Ra, the boundary layer flow is laminar; hence the correlation of Eq..7 is appropriate, 1/4 h 0.70 Ra = = 0.8 + (5 k 4/ /1 1+ ( 0.4 / Pr 0.05 W / m K 8 ( 1/4 /1 ( 4/ h = 0.8 + 0.70.07 / 1 + 0.4 / 0.707 =.38 W / m K. 1m From Eqs. (1 and ( with numerical values for h and ΔT /, find q s =.38W / m K.4 o C = 4.8W / m. COMMENTS: Recognize that radiation exchange with the environment will be significant. 4 4 Assuming Ts T /, Tsur T q = σ T T = W/m. = = and ε = 1, find ( rad s sur
PROBEM.51 KNOWN: Plate, 1m 1m, inclined at 45 from the vertical is exposed to a net radiation heat flux of 300 W/m ; backside of plate is insulated and ambient air is at 0 C. FIND: Temperature plate reaches for the prescribed conditions. ASSUMPTIONS: (1 Net radiation heat flux (300 W/m includes exchange with surroundings, ( Ambient air is quiescent, (3 No heat losses from backside of plate, (4 Steady-state conditions. PROPERTIES: Table A-4, Air (assuming T s = 84 C, T f = (T s + T / = (84 + 0 C/ = 315K, 1 atm: ν = 17.40 - m /s, k = 0.074 W/m K, α = 4.7 - m /s, Pr = 0.705, β = 1/T f. ANAYSIS: From an energy balance on the plate, it follows that q rad = q conv. That is, the net radiation heat flux into the plate is equal to the free convection heat flux to the ambient air. The temperature of the surface can be expressed as Ts = T + q rad /h (1 where h must be evaluated from an appropriate correlation. Since this is the bottom surface of a heated inclined plate, g may be replaced by g cos θ ; hence using Eq..5, find g cosθβ ( Ts T 3.8m / s cos 45 ( 1/ 315K( 84 0 K ( 1m 3 Ra 4.30 = = =. να 17.40 m / s 4.7 m / s Since Ra >, conditions are turbulent and Eq.. is appropriate for estimating 0.387 Ra 1/ 0.85 = + /1 8/7 ( 1+ ( 0.4 / Pr 1/ 0.387 ( 4.30 = 0.85 + 13. 8/7 = /1 1 ( 0.4 / 0.705 + h k / 13. 0.074 W / m K /1m 5. W / m = = = K. (3 Substituting h from Eq. (3 into Eq. (1, the plate temperature is T 0 C 300 W / m / 5. W / m s = + K = 57 C. COMMENTS: Note that the resulting value of T s 57 C is substantially lower than the assumed value of 84 C. The calculation should be repeated with a new estimate of T s, say, 0 C. An alternate approach is to write Eq. ( in terms of T s, the unknown surface temperature and then combine with Eq. (1 to obtain an expression which can be solved, by trial-and-error, for T s.
PROBEM.55 KNOWN: Dimensions and temperature of beer can in refrigerator compartment. FIND: Orientation which maximizes cooling rate. ASSUMPTIONS: (1 End effects are negligible, ( Compartment air is quiescent, (3 Constant properties. PROPERTIES: Table A-4, Air (T f = 88.5K, 1 atm: ν = 14.87 - m /s, k =0.054 W/m K, α = 1.0 - m /s, Pr = 0.71, β = 1/T f = 3.47-3 K -1. ANAYSIS: The ratio of cooling rates may be expressed as qv hv π D ( Ts T h = = v. qh hh π D ( Ts T hh For the vertical surface, find gβ ( T 3 1 s T 3.8m / s 3.47 K ( 3 C Ra 3 3 = = =.5 να 14.87 m / s 1 m / s Ra =.5 0.15 = 8.44, ( 3 ( ( 0.387( 8.44 and using the correlation of Eq.., = 0.85 +.7. /1 8/7 = 1 + 0.4 / 0.71 ( k 0.054 W / m K Hence h = hv = =.7 = 5.03 W / m K. 0.15m gβ ( T For the horizontal surface, find s T RaD = D 3 =.5 ( 0.0 3 = 5.4 5 να 5 0.387( 5.4 and using the correlation of Eq..34, Hence ( 1/ 1/ D = 0.0 + 1.4 /1 8/7 = 1+ 0.55/0.71 k 0.054 W / m K h D = hh = D = 1.4 = 5.18 W / m K. D 0.0m qv 5.03 = = 0.7. qh 5.18 COMMENTS: In view of the uncertainties associated with Eqs.. and.34 and the neglect of end effects, the above result is inconclusive. The cooling rates are approximately the same.
PROBEM.1 KNOWN: Dissipation rate of an immersion heater in a large tank of water. FIND: Surface temperature in water and, if accidentally operated, in air. ASSUMPTIONS: (1 Quiescent ambient fluid, ( Negligible radiative exchange. PROPERTIES: Table A-, Water and Table A-4, Air: T(K k 3 (W/m K ν 7 (μ/ρ,m /s α 7 (k/ρc p,m /s Pr β (K -1 Water 315 34.5 1.531 4.1 400.4 Air 1500 0 400 3500 085.7 ANAYSIS: From the rate equation, the surface temperature, T s, is ( π Ts = T + q/ Dh (1 where h is estimated by an appropriate correlation. Since such a calculation requires knowledge of T s, an iteration procedure is required. Begin by assuming for water that T s = 4 C such that T f = 315K. Calculate the Rayleigh number, 1 ( ( 3 3.8m / s 400.4 K gβδtd 4 0 K 0.0m Ra D = = = 1.804. να.5 7 m / s 1.531 7 m / s Using the Churchill-Chu relation, find 1/ hd 0.387 Ra D D = = 0.0 + k /1 8/7 1 + ( 0.55 / Pr 1/ 0.34 W / m K 0.387( 1.804 h = 0.0 1301W / m + K. 0.01m /1 8/7 = 1 + ( 0.55 / 4.1 Substituting numerical values into Eq. (1, the calculated value for T s in water is T s = 0 C + 550 W / π 0.0m 0.30m 1301W / m K = 4.8 C. Continued ( (3
PROBEM.1 (Cont. Our initial assumption of T s = 4 C is in excellent agreement with the calculated value. With accidental operation in air, the heat transfer coefficient will be nearly a factor of 0 less. Suppose h 5W/m K, then from Eq. (1, T s 30 C. Very likely the heater will burn out. Using air properties at T f 1500K and Eq. (, find Ra D = 1.815. Using Eq..33, n D = C Ra D with C= 0.85 and n = 0.188 from Table.1, find h =.W/m K. Hence, our first estimate for the surface temperature in air was reasonable, Ts 300 C. However, radiation exchange will be the dominant mode, and would reduce the estimate for T s. Generally such heaters could not withstand operating temperatures above 00 C and safe operation in air is not possible.
PROBEM.8 KNOWN: Vertical air vent in front door of dishwasher with prescribed width and height. Spacing between isothermal and insulated surface of 0 mm. FIND: (a Heat loss from the tub surface and (b Effect on heat rate of changing spacing by ± mm. ASSUMPTIONS: (1 Steady-state conditions, ( Vent forms vertical parallel isothermal/adiabatic plates, (3 Ambient air is quiescent. PROPERTIES: Table A-4, (T f = (T s + T / = 31.5K, 1 atm: ν = 17.15 - m /s, α = 4.4 - m /s, k = 7. -3 W/m K, β = 1/T f. ANAYSIS: The vent arrangement forms two vertical plates, one is isothermal, T s, and the other is adiabatic ( q = 0. The heat loss can be estimated from Eq..37 with the correlation of Eq..45 using C 1 = 144 and C =.87 from Table.3: gβ ( Ts T S 3.8m / s ( 1/ 31.5 K( 5 7 K ( 0.00 m 3 RaS = = = 14,88 να 17.15 m / s 4.4 m / s 1/ k C ( 1 C q = A s Ts T + = ( 0.500 0.580 m S 1/ ( RaSS/ ( RaSS/ 1/ 0.07 W / m K C ( 1 C 5 7 K + = 8.8W. 0.00 m 1/ ( RaSS/ ( RaSS/ (b To determine the effect of the spacing at S = 30 and mm, we need only repeat the above calculations with these results S (mm Ra S q (W 1874.1 30 50,585 8.8 Since it would be desirable to minimize heat losses from the tub, based upon these calculations, you would recommend a decrease in the spacing. COMMENTS: For this situation, according to Table.3, the spacing corresponding to the maximum heat transfer rate is S max = (S max /S opt.15(ra S /S 3-1/4 = 14.5 mm. Find q max = 8.5 W. Note that the heat rate is not very sensitive to spacing for these conditions.
PROBEM. KNOWN: Horizontal flat roof and vertical wall sections of same dimensions exposed to identical temperature differences. FIND: (a Ratio of convection heat rate for horizontal section to that of the vertical section and (b Effect of inserting a baffle at the mid-height of the vertical wall section on the convection heat rate. ASSUMPTIONS: (1 Ends of sections and baffle adiabatic, ( Steady-state conditions. PROPERTIES: Table A-4, Air ( ( T = T1+ T /= 77K,1atm : ν = 13.84 - m /s, k = 0.045 W/m K, α = 1.5 - m /s, Pr = 0.713. ANAYSIS: (a The ratio of the convection heat rates is qhor hhorasδt h = = hor. (1 qvert hvertasδt hvert To estimate coefficients, recognizing both sections have the same characteristics length, = 0.1m, with Ra = gβδt 3 /να find ( ( 3 ( (.8 m / s 1/ 77K 18 K 0.1m Ra = = 3.7. 13.84 m / s 1.5 m / s The appropriate correlations for the sections are Eqs..4 and.5 (with H/ = 30, 1/ 3 0.074 1/ 4 0.01 0.3 hor = 0.0 Ra Pr vert = 0.4 Ra Pr ( H /. (3,4 Using Eqs. (3 and (4, the ratio of Eq. (1 becomes, 1/3 0.074 ( ( ( ( ( 1/3 0.074 qhor 0.0 Ra Pr 0.0 3.7 0.713 = = = 1.57. q 1/4 0.01 0.3 1/ 4 vert 0.4 Ra Pr ( H / 0.01 0.3 0.4 3.7 0.713 30 (b The effect of the baffle in the vertical wall section is to reduce H/ from 30 to 15. Using Eq..5, it follows, 0.3 ( H/ 0.3 qbaf hbaf baf 15 = = = 1.3. q h 0.3 = ( H/ 30 That is, the effect of the baffle is to increase the convection heat rate. COMMENTS: (1 Note that the heat rate for the horizontal section is 57% larger than that for the vertical section for the same (T 1 T. This indicates the importance of heat losses from the ceiling or roofs in house construction. ( Recognize that for Eq..5, the Pr > 1 requirement is not completely satisfied. (3 What is the physical explanation for the result of part (b?
PROBEM.111 KNOWN: Vertical array of circuit boards 0.15m high with maximum allowable uniform surface temperature for prescribed ambient air temperature. FIND: Allowable electrical power dissipation per board, q [ W/m ], for these cooling arrangements: (a Free convection only, (b Air flow downward at 0. m/s, (c Air flow upward at 0.3 m/s, and (d Air flow upward or downward at 5 m/s. ASSUMPTIONS: (1 Uniform surface temperature, ( Board horizontal spacing sufficient that boundary layers don t interfere, (3 Ambient air behaves as quiescent medium, (4 Perfect gas behavior. PROPERTIES: Table A-4, Air (T f = (T s + T / 315K, 1 atm: ν = 17.40 - m /s, k = 0.074 W/m K, α = 4.7 - m /s, Pr = 0.705, β = 1/T f. ANAYSIS: (a For free convection only, the allowable electrical power dissipation rate is q = h( ( Ts T (1 where h is estimated using the appropriate correlation for free convection from a vertical plate. Find the Rayleigh number, gβ ΔT 3.8m / s ( 1/ 315K( 0 5 K ( 0.150m 3 Ra = = = 8.551. ( να 17.4 m / s 4.7 m / s Since Ra, the flow is laminar. With Eq..7 find 1/4 1/4 0.70 8.551 h 0.70 Ra = = 0.8 + = 0.8 + = 8.47 (3 k /1 4/ /1 4/ 1 + ( 0.4 / Pr 1 + ( 0.4 / 0.705 ( h = 0.074 W / m K / 0.150m 8.47 = 5.0 W / m K. Hence, the allowable electrical power dissipation rate is, q = 5.0 W / m K ( 0.150m( 0 5 C = 54. W / m. (b With downward velocity V = 0. m/s, the possibility of mixed forced-free convection must be considered. With Re = V/ν, find ( Ra Gr / Re / Re = Pr (4 ( ( ( Gr / Re = 8.551 / 0.705 / 0. m / s 0.150m /17.40 m / s = 0.453. Continued
Since ( PROBEM.111 (Cont. Gr / Re ~ 1, flow is mixed and the average heat transfer coefficient may be found from a correlating equation of the form n = n n F ± N (5 where n = 3 for the vertical plate geometry and the minus sign is appropriate since the natural convection (N flow opposes the forced convection (F flow. For the forced convection flow, Re = 517 and the flow is laminar; using Eq. 7.30, 1/ 1/3 F = 0.4 Re 1/ 1/3 Pr = 0.4( 517 ( 0.705 = 4.50. ( Using = 8.47 from Eq. (3, Eq. (5 now becomes N 3 3 h 3 3 = = ( 4.50 ( 8.47 = 37.7 k 0.074 W / m K h = 37.7.8 W / m = K. 0.150m Substituting for h into the rate equation, Eq. (1, the allowable power dissipation with a downward velocity of 0. m/s is q =.8 W / m K ( 0.150m ( 0 5 C = 7.3 W / m. (c With an upward velocity V = 0.3 m/s, the positive sign of Eq. (5 applies since the N-flow is assisting the F-flow. For forced convection, find Re = V / ν = 0.3m / s 0.150m / ( 17.40 m / s = 58. The flow is again laminar, hence Eq. ( is appropriate. 1/ 1/3 F = 0.4( 58 ( 0.705 = 30.05. From Eq. (5, with the positive sign, and from Eq. (4, N 3 3 3 = ( 30.05 + ( 8.47 or = 3.88 and h =.74 W / m K. From Eq. (1, the allowable power dissipation with an upward velocity of 0.3 m/s is q =.74 W / m K ( 0.150m ( 0 5 C = 70.7 W / m. (d With a forced convection velocity V = 5 m/s, very likely forced convection will dominate. Check by evaluating whether ( 43,3. Hence, Gr / Re 1 where Re = V/ν = 5 m/s 0.150m/(17.40 - m /s = Ra ( ( Gr / Re = / Re = 8.551 / 0.705 / 43,3 = 0.007. Pr The flow is not mixed, but pure forced convection. Using Eq. (, find 1/ 1/3 h = ( 0.074 W / m K / 0.150m 0.4( 43,3 ( 0.705 =.4 W / m K and the allowable dissipation rate is q =.4 W / m K 0.150m 0 5 C = 35 W / m. ( ( COMMENTS: Be sure to compare dissipation rates to see relative importance of mixed flow conditions.