42540_06 - CHAPTER 6 MASS TRANSFER Molecular Diffusion

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: CHAPTER 6 MASS TRANSFER Molecular Diffusion ........................................................................................................................ 6.1 Convection of Mass ......................................................................................................................... 6.5 Simultaneous Heat and Mass Transfer Between Water-Wetted Surfaces and Air.......................... 6.9 Symbols ......................................................................................................................................... 6.13 ASS transfer by either molecular diffusion or convection is the transport of one component of a mixture relative to the motion of the mixture and is the result of a concentration gradient. Mass transfer can occur in liquids and solids as well as gases. For example, water on the wetted slats of a cooling tower evaporates into air in a cooling tower (liquid to gas mass transfer), and water vapor from a food product transfers to the dry air as it dries. A piece of solid CO2 (dry ice) also gets smaller and smaller over time as the CO2 molecules diffuse into air (solid to gas mass transfer). A piece of sugar added to a cup of coffee eventually dissolves and diffuses into the solution, sweetening the coffee, although the sugar molecules are much heavier than the water molecules (solid to liquid mass transfer). Air freshener does not just smell where sprayed, but rather the smell spreads throughout the room. The air freshener (matter) moves from an area of high concentration where sprayed to an area of low concentration far away. In an absorption chiller, lowpressure, low-temperature refrigerant vapor from the evaporator enters the thermal compressor in the absorber section, where the refrigerant vapor is absorbed by the strong absorbent (concentrated solution) and dilutes the solution. In air conditioning, water vapor is added or removed from the air by simultaneous transfer of heat and mass (water vapor) between the airstream and a wetted surface. The wetted surface can be water droplets in an air washer, condensate on the surface of a dehumidifying coil, a spray of liquid absorbent, or wetted surfaces of an evaporative condenser. Equipment performance with these phenomena must be calculated carefully because of simultaneous heat and mass transfer. This chapter addresses mass transfer principles and provides methods of solving a simultaneous heat and mass transfer problem involving air and water vapor, emphasizing air-conditioning processes. The formulations presented can help analyze performance of specific equipment. For discussion of performance of cooling coils, evaporative condensers, cooling towers, and air washers, see Chapters 22, 38, 39, and 40, respectively, of the 2008 ASHRAE Handbook—HVAC Systems and Equipment. M the mixture of gases A and B in a direction that reduces the concentration gradient. Fick’s Law The basic equation for molecular diffusion is Fick’s law. Expressing the concentration of component B of a binary mixture of components A and B in terms of the mass fraction B/ or mole fraction CB /C, Fick’s law is dB J B = – D v --------------------- = – J A dy d CB C * * J B = – CD v ---------------------- = – J A dy (1a) (1b) where = A + B and C = CA + CB. The minus sign indicates that the concentration gradient is negative in the direction of diffusion. The proportionality factor Dv is the mass diffusivity or the diffusion coefficient. The total mass · · flux m B and molar flux m B* are due to the average velocity of the mixture plus the diffusive flux: · mB = Bv dB – D v --------------------dy (2a) d CB C · m B* = C B v * – CD v ---------------------dy (2b) MOLECULAR DIFFUSION Most mass transfer problems can be analyzed by considering diffusion of a gas into a second gas, a liquid, or a solid. In this chapter, the diffusing or dilute component is designated as component B, and the other component as component A. For example, when water vapor diffuses into air, the water vapor is component B and dry air is component A. Properties with subscripts A or B are local properties of that component. Properties without subscripts are local properties of the mixture. The primary mechanism of mass diffusion at ordinary temperature and pressure conditions is molecular diffusion, a result of density gradient. In a binary gas mixture, the presence of a concentration gradient causes transport of matter by molecular diffusion; that is, because of random molecular motion, gas B diffuses through The preparation of this chapter is assigned to TC 1.3, Heat Transfer and Fluid Flow. where v is the mixture’s mass average velocity and v* is the molar average velocity. Bird et al. (1960) present an analysis of Equations (1a) and (1b). Equations (1a) and (1b) are equivalent forms of Fick’s law. The equation used depends on the problem and individual preference. This chapter emphasizes mass analysis rather than molar analysis. However, all results can be converted to the molar form using the relation CB B /MB. Fick’s Law for Dilute Mixtures In many mass diffusion problems, component B is dilute, with a density much smaller than the mixture’s. In this case, Equation (1a) can be written as J B = – Dv dB dy (3) when B << and A . Equation (3) can be used without significant error for water vapor diffusing through air at atmospheric pressure and a temperature less than 80°F. In this case, B < 0.02 , where B is the density of water vapor and is the density of moist air (air and water vapor mixture). The error in JB caused by replacing [d( B/ )/dy] with d B /dy is less than 2%. At temperatures below 140°F where B < 0.10 , Equation (3) can still be used if errors in JB as great as 10% are tolerable. 6.1 6.2 Fick’s Law for Mass Diffusion Through Solids or Stagnant Fluids (Stationary Media) Fick’s law can be simplified for cases of dilute mass diffusion in solids, stagnant liquids, or stagnant gases. In these cases, B << and v 0, which yields the following approximate result: dB · m B = J B = – D v -------dy (4) 2009 ASHRAE Handbook—Fundamentals In the absence of data, use equations developed from (1) theory or (2) theory with constants adjusted from limited experimental data. For binary gas mixtures at low pressure, Dv is inversely proportional to pressure, increases with increasing temperature, and is almost independent of composition for a given gas pair. Bird et al. (1960) present the following equation, developed from kinetic theory and corresponding states arguments, for estimating Dv at pressures less than 0.1pc min: T D v = a ----------------------T cA T cB b Fick’s Law for Ideal Gases with Negligible Temperature Gradient For dilute mass diffusion, Fick’s law can be written in terms of partial pressure gradient instead of concentration gradient. When gas B can be approximated as ideal, B Ru T p B = --------------- = C B R u T MB 11------- + ------MA MB 5 12 (8) pcA p cB T cA T cB --------------------------------------------------------------p where Dv a b T p M = = = = = = diffusion coefficient, ft2/h constant, dimensionless constant, dimensionless absolute temperature, °R pressure, atm relative molecular weight, lbm/lb mol 13 (5) and when the gradient in T is small, Equation (3) can be written as MB Dv d pB J B = – -------------- -------Ru T dy or Dv d pB * J B = – --------- -------R u T dy If v 0, Equation (4) may be written as MB Dv d pB · m B = J B = – -------------- -------Ru T dy or Dv d pB * · m B * = J B = – --------- ---------R u T dy (7b) (7a) (6b) (6a) Subscripts cA and cB refer to the critical states of the two gases. Analysis of experimental data gives the following values of the constants a and b: For nonpolar gas pairs a = 6.518 × 10 4 and b = 1.823 For water vapor with a nonpolar gas a = 8.643 × 10–4 and b = 2.334 In nonpolar gas, intermolecular forces are independent of the relative orientation of molecules, depending only on the separation distance from each other. Air, composed almost entirely of nonpolar gases O2 and N2, is nonpolar. Equation (8) is stated to agree with experimental data at atmospheric pressure to within about 8% (Bird et al. 1960). Mass diffusivity Dv for binary mixtures at low pressure is predictable within about 10% by kinetic theory (Reid et al. 1987). D v = 2.79 where AB D, AB The partial pressure gradient formulation for mass transfer analysis has been used extensively; this is unfortunate because the pressure formulation [Equations (6) and (7)] applies only when one component is dilute, the fluid closely approximates an ideal gas, and the temperature gradient has a negligible effect. The density (or concentration) gradient formulation expressed in Equations (1) to (4) is more general and can be applied to a wider range of mass transfer problems, including cases where neither component is dilute [Equation (1)]. The gases need not be ideal, nor the temperature gradient negligible. Consequently, this chapter emphasizes the density formulation. –5 T 1.5 1110 -------------------------------------- ------- + ------2 MA MB p AB D, AB (9) = characteristic molecular diameter, nm = temperature function, dimensionless Diffusion Coefficient For a binary mixture, the diffusion coefficient Dv is a function of temperature, pressure, and composition. Experimental measurements of Dv for most binary mixtures are limited in range and accuracy. Table 1 gives a few experimental values for diffusion of some gases in air. For more detailed tables, see the Bibliography. Table 1 Gas Ammonia Benzene Carbon dioxide Ethanol Hydrogen Oxygen Water vapor *Gases at 77°F and 14.696 psi. Mass Diffusivities for Gases in Air* Dv , ft2/h 1.08 0.34 0.64 0.46 1.60 0.80 0.99 Dv is in ft2/h, p in atm, and T in °R. If the gas molecules of A and B are considered rigid spheres having diameters A and B [and AB = ( A/2) + ( B/2)], all expressed in nanometers, the dimensionless function D ,AB equals unity. More realistic models for molecules having intermolecular forces of attraction and repulsion lead to values that are functions of temperature. Bird et al. (1960) and Reid et al. (1987) present tabulations of D ,AB . These results show that Dv increases as the 2.0 power of T at low temperatures and as the 1.65 power of T at very high temperatures. The diffusion coefficient of moist air has been calculated for Equation (8) using a simplified intermolecular potential field function for water vapor and air (Mason and Monchick 1965). The following empirical equation is for mass diffusivity of water vapor in air up to 2000°F (Sherwood and Pigford 1952): –4 T 2.5 1.46 10 D v = --------------------------- ----------------T + 441 p (10) Mass Transfer Example 1. Evaluate the diffusion coefficient of CO2 in air at 527.4°R and atmospheric pressure (14.696 psia) using Equation (9). Solution: In Equation (9), Dv is in ft2/h, AB = ( Lennard-Jones energy parameter A B k = of and for each gas are as follows: , nm CO2 Air 0.3996 0.3617 A/2) A 6.3 Air, then, diffuses toward the liquid water interface. Because it cannot be absorbed there, a bulk velocity v of the gas mixture is established in a direction away from the liquid surface, so that the net transport of air is zero (i.e., the air is stagnant): dA · m A = – D v -------- + dy Av +( B B/2), k and the k . Values /k, °R 342 175 =0 (13) The combined values for use in Equation (9) are AB AB = 0.3996/2 + 0.3617/2 = 0.3806 nm The bulk velocity v transports not only air but also water vapor away from the interface. Therefore, the total rate of water vapor diffusion is dB · m B = – D v -------- + dy Bv k= AB 342 175 = 245 at P = 1 atm and T = 527.4°R 1 kT ------- = ------------ = 2.15 0.465 AB (14) 245 ------- = ------------ = 0.465 kT 527.4 From tables for D ,AB at kT/ AB = 2.15 (Bird et al. 1960; Reid et al. 1987), the collision integral D,AB = 1.047. The molecular weights of CO2 and air are 44 and 29, respectively. Substituting these values gives Dv = 2.79 × 10–5[527.41.5/(1 × 0.38062 × 1.047)](1/44 + 1/29)0.5 = 0.533 ft2/h Substituting for the velocity v from Equation (13) and using Equations (11b) and (12) gives D v MB p d A · m B = ----------------- -------dy A Ru T Integration yields D v M B p ln A L A 0 · m B = ------------------ -------------------------------Ru T yL – y0 BL – B0 · m B = – D v PAm ----------------------yL – y0 (15) Diffusion of One Gas Through a Second Stagnant Gas Figure 1 shows diffusion of one gas through a second, stagnant gas. Water vapor diffuses from the liquid surface into surrounding stationary air. It is assumed that local equilibrium exists through the gas mixture, that the gases are ideal, and that the Gibbs-Dalton law is valid, which implies that temperature gradient has a negligible effect. Water vapor diffuses because of concentration gradient, as given by Equation (6a). There is a continuous gas phase, so the mixture pressure p is constant, and the Gibbs-Dalton law yields pA + pB = p = constant or p A B ------- + ------- = --------- = constant Ru T MA MB (11a) (11b) (16a) or (16b) where P Am p -------p AL AL ln A L A 0 -------------------------------AL – A0 (17) The partial pressure gradient of the water vapor causes a partial pressure gradient of the air such that dp B dpA -------- = – -------dy dy or 1 dB 1 dA ------- -------- = – ------- -------MB dy MA dy (11c) (12) PAm is the logarithmic mean density factor of the stagnant air. The pressure distribution for this type of diffusion is illustrated in Figure 2. Stagnant refers to the net behavior of the air; it does not move because bulk flow exactly offsets diffusion. The term PAm in Equation (16b) approximately equals unity for dilute mixtures such as water vapor in air at near-atmospheric conditions. This condition makes it possible to simplify Equations (16) and implies that, for dilute mixtures, the partial pressure distribution curves in Figure 2 are straight lines. Example 2. A vertical tube of 1 in. diameter is partially filled with water so that the distance from the water surface to the open end of the tube is 2.362 in., as shown in Figure 1. Perfectly dried air is blown over the open tube end, and the complete system is at a constant temperature of 59°F. In 200 h of steady operation, 0.00474 lb of water evaporates from the tube. The total pressure of the system is 14.696 psia (1 atm). Using these data, (1) calculate the mass diffusivity of water vapor in air, and (2) compare this experimental result with that from Equation (10). Solution: (1) The mass diffusion flux of water vapor from the water surface is · m B = 0.00474 200 = 0.0000237 lb/h The cross-sectional area of a 1 in. diameter tube is (0.5)2/144 = · 0.005454 ft2. Therefore, m B = 0.004345 lb/ft2 ·h. The partial densities are determined from psychrometric tables. BL AL Fig. 1 Diffusion of Water Vapor Through Stagnant Air = 0; B0 = 0.000801 lb/ft 3 A0 3 3 = 0.0765 lb/ft ; = 0.0752 lb/ft Fig. 1 Diffusion of Water Vapor Through Stagnant Air Because p = pAL = 1 atm, the logarithmic mean density factor [Equation (17)] is 6.4 ln 0.0765 0.0752 P Am = 0.0765 ---------------------------------------------- = 1.009 0.0765 – 0.0752 The mass diffusivity is now computed from Equation (16b) as · –m B y L – y 0 – 0.004345 2.362 D v = -------------------------------------- = ----------------------------------------------------------------P Am B L – B 0 1.009 0 – 0.000801 12 = 1.058 ft 2 2009 ASHRAE Handbook—Fundamentals At steady state, the molar flows of A and B must be equal but opposite: ·* ·* mA + mB = 0 (18) h because the total molar concentration C must stay the same in both chambers if p and T remain constant. Because molar fluxes are the same in both directions, the molar average velocity v * = 0. Thus, Equation (7b) can be used to calculate the molar flux of B (or A): –Dv d pB ·* m B = --------- -------R u T dy or A cs D v p B 0 – p BL ·* m B = -------------- ----------------------L Ru T MB Acs D v p B 0 – p BL · m B = ---------------------- ----------------------L Ru T (19) (2) By Equation (10), with p = 14.696 psi and T = 59 + 460 = 519°R, 519 2 D v = 0.00215 ----------------------- = 0.935 ft h -----------------14.696 519 + 441 2.5 Neglecting the correction factor PAm for this example gives a difference of less than 1% between the calculated experimental and empirically predicted values of Dv. (20) Equimolar Counterdiffusion Figure 3 shows two large chambers, both containing an ideal gas mixture of two components A and B (e.g., air and water vapor) at the same total pressure p and temperature T. The two chambers are connected by a duct of length L and cross-sectional area Acs. Partial pressure pB is higher in the left chamber, and partial pressure pA is higher in the right chamber. The partial pressure differences cause component B to migrate to the right and component A to migrate to the left. or (21) Example 3. One large room is maintained at 70°F (530°R), 1 atm, 80% rh. A 60 ft long duct with cross-sectional area of 1.5 ft2 connects the room to another large room at 70°F, 1 atm, 10% rh. What is the rate of water vapor diffusion between the two rooms? Solution: Let air be component A and water vapor be component B. Equation (21) can be used to calculate the mass flow of water vapor B. Equation (10) can be used to calculate the diffusivity. –4 530 1.46 10 D v = -------------------------- ----------------------530 + 441 1 2.5 = 0.972 ft /h 2 Fig. 2 Pressure Profiles for Diffusion of Water Vapor Through Stagnant Air From a psychrometric table (Table 3, Chapter 1), the saturated vapor pressure at 70°F is 0.363 psi. The vapor pressure difference pB 0 – pBL is p B 0 – p BL = 0.8 – 0.1 0.363 psi = 0.254 psi Then, Equation (21) gives 18 144 1.5 0.972 0.254 · m b = ------------------------------------------------------------------------- = 1.90 1545 530 60 10 – 5 lb m /h Molecular Diffusion in Liquids and Solids Because of the greater density, diffusion is slower in liquids than in gases. No satisfactory molecular theories have been developed for calculating diffusion coefficients. The limited measured values of Dv show that, unlike for gas mixtures at low pressures, the diffusion coefficient for liquids varies appreciably with concentration. Reasoning largely from analogy to the case of one-dimensional diffusion in gases and using Fick’s law as expressed by Equation (4) gives Fig. 2 Pressure Profiles for Diffusion of Water Vapor Through Stagnant Air B1 – B2 · m B = D v ----------------------y2 – y1 (22) Fig. 3 Equimolar Counterdiffusion Fig. 3 Equimolar Counterdiffusion Equation (22) expresses steady-state diffusion of solute B through solvent A in terms of the molal concentration difference of the solute at two locations separated by the distance y = y2 – y1. Bird et al. (1960), Eckert and Drake (1972), Hirschfelder et al. (1954), Reid and Sherwood (1966), Sherwood and Pigford (1952), and Treybal (1980) provide equations and tables for evaluating Dv . Hirschfelder et al. (1954) provide comprehensive treatment of the molecular developments. Diffusion through a solid when the solute is dissolved to form a homogeneous solid solution is known as structure-insensitive diffusion (Treybal 1980). This solid diffusion closely parallels diffusion through fluids, and Equation (22) can be applied to onedimensional steady-state problems. Values of mass diffusivity are generally lower than they are for liquids and vary with temperature. Mass Transfer Diffusion of a gas mixture through a porous medium is common (e.g., diffusion of an air/vapor mixture through porous insulation). Vapor diffuses through the air along the tortuous narrow passages within the porous medium. Mass flux is a function of the vapor pressure gradient and diffusivity, as indicated in Equation (7a). It is also a function of the structure of the pathways within the porous medium and is therefore called structure-sensitive diffusion. All these factors are taken into account in the following version of Equation (7a): · mB = – dpB -------dy (23) where hM = internal mass transfer coefficient, ft/h · m B = mass flux of gas B at interfacial surface, lbm/ft2 ·h 3 Bi = density of gas B at interfacial surface, lbm/ft Bb 6.5 Mass Transfer Coefficients for Internal Flows. Most internal convective mass transfer problems, such as those that occur in channels or in the cores of dehumidification coils, can be solved if an appropriate expression is available to relate the mass transfer flux (to or from the interfacial surface) to the difference between the concentration at the surface and the bulk concentration in the channel, as shown in Figure 5. This formulation leads to the definition of the mass transfer coefficient for internal flows: hM · mB ---------------------Bi – Bb (27) where is called the permeability of the porous medium. Chapter 25 presents this topic in more depth. CONVECTION OF MASS Convection of mass involves the mass transfer mechanisms of molecular diffusion and bulk fluid motion. Fluid motion in the region adjacent to a mass transfer surface may be laminar or turbulent, depending on geometry and flow conditions. 1 u B A cs 1 A cs A A cs uB B dA cs = bulk density of gas B at location x uB u B dA cs = average velocity of gas B at location x, Mass Transfer Coefficient Convective mass transfer is analogous to convective heat transfer where geometry and boundary conditions are similar. The analogy holds for both laminar and turbulent flows and applies to both external and internal flow problems. Mass Transfer Coefficients for External Flows. Most external convective mass transfer problems can be solved with an appropriate formulation that relates the mass transfer flux (to or from an interfacial surface) to the concentration difference across the boundary layer illustrated in Figure 4. This formulation gives rise to the convective mass transfer coefficient, defined as · mB (24) h M ----------------------Bi – B where hM = local external mass transfer coefficient, ft/h 2 · m B = mass flux of gas B from surface, lbm/ft ·h 3 Bi = density of gas B at interface (saturation density), lb m/ft 3 B = density of component B outside boundary layer, lbm/ft fpm Acs = cross-sectional area of channel at station x, ft2 uB = velocity of component B in x direction, fpm 3 B = density distribution of component B at station x, lbm/ft Often, it is easier to obtain the bulk density of gas B from · · m Bo + m B dA A = ----------------------------------u B A cs Bb (28) where · m Bo = mass flow rate of component B at station x = 0, lbm /h A = interfacial area of channel between station x = 0 and station x = x, ft2 Equation (28) can be derived from the preceding definitions. The major problem is the determination of u B . If, however, analysis is restricted to cases where B is dilute and concentration gradients of B in the x direction are negligibly small, u B u . Component B is swept along in the x direction with an average velocity equal to the average velocity of the dilute mixture. If Bi and B are constant over the entire interfacial surface, the mass transfer rate from the surface can be expressed as · mB = hM Bi Analogy Between Convective Heat and Mass Transfer Most expressions for the convective mass transfer coefficient hM are determined from expressions for the convective heat transfer coefficient h. For problems in internal and external flow where mass transfer occurs at the convective surface and where component B is dilute, Bird et al. (1960) and Incropera and DeWitt (1996) found that the Nusselt and Sherwood numbers are defined as follows: Nu = f (X, Y, Z, Pr, Re) (29) (30) – B (25) where h M is the average mass transfer coefficient, defined as hM 1 -A A h m dA (26) Fig. 4 Nomenclature for Convective Mass Transfer from External Surface at Location x Where Surface is Impermeable to Gas A Sh = f (X, Y, Z, Sc, Re) Fig. 5 Nomenclature for Convective Mass Transfer from Internal Surface Impermeable to Gas A Fig. 4 Nomenclature for Convective Mass Transfer from External Surface at Location x Where Surface Is Impermeable to Gas A Fig. 5 Nomenclature for Convective Mass Transfer from Internal Surface Impermeable to Gas A 6.6 and Nu = g Pr, Re Sh = g Sc, Re (31) (32) 2009 ASHRAE Handbook—Fundamentals Sh = 0.664Re Sh = 0.037Re With Nu and Sh known, Sh D v h M = ------------L or Pr h Nu k ------ = ------------- = ----Sc hM Sh D v 1 /3 1/2 4/5 Sc Sc 1/3 1/3 for laminar flow for turbulent flow where f in Equations (29) and (30) indicates a functional relationship among the dimensionless groups shown. The function f is the same in both equations. Similarly, g indicates a functional relationship that is the same in Equations (31) and (32). Pr and Sc are dimensionless Prandtl and Schmidt numbers, respectively, as defined in the Symbols section. The primary restrictions on the analogy are that the surface shapes are the same and that the temperature boundary conditions are analogous to the density distribution boundary conditions for component B when cast in dimensionless form. Several primary factors prevent the analogy from being perfect. In some cases, the Nusselt number was derived for smooth surfaces. Many mass transfer problems involve wavy, droplet-like, or roughened surfaces. Many Nusselt number relations are obtained for constant-temperature surfaces. Sometimes Bi is not constant over the entire surface because of varying saturation conditions and the possibility of surface dryout. In all mass transfer problems, there is some blowing or suction at the surface because of condensation, evaporation, or transpiration of component B. In most cases, this blowing/suction has little effect on the Sherwood number, but the analogy should be examined closely if vi /u > 0.01 or v i u > 0.01, especially if the Reynolds number is large. Example 4. Air at 77°F, 1 atm, and 60% rh flows at 32.8 ft/s (1970 ft/min), as shown in Figure 6. Find the rate of evaporation, rate of heat transfer to the water, and water surface temperature. Solution: Heat transfer to water from air supplies the energy required to evaporate the water. · q = hA t – t s = m h fg = h M A where h hM A · m ts, s t, s Nu k h = ----------L k ----Dv This result is valid for both laminar and turbulent flow. Using this result in the preceding energy balance gives Pr = ----Sc 1 /3 s– k- t – t s ----- -----------------h fg Dv This equation must be solved for s . Then, water surface temperature ts is the saturation temperature corresponding to s. Air properties Sc, Pr, Dv , and k are evaluated at film temperature tf = (t + ts)/2, and hfg is evaluated at ts. Because ts appears in the right side and all the air properties also vary somewhat with ts, iteration is required. Start by guessing ts = 57.2°F (the dew point of the airstream), giving tf = 67.1°F. At these temperatures, values on the right side are found in property tables or calculated as k = 0.01485 Btu/h·ft·°F Pr = 0.709 Dv = 0.9769 ft2/h = 0.01628 ft2/min [from Equation (10)] = 0.07435 lbm/ft3 = 0.04376 lbm/ft·h Sc = / Dv = 0.6025 hfg = 1055.5 Btu/lbm (at 57.2°F) = 8.846 × 10–4 lbm/ft3 (from psychrometric chart at 77°F, 60% rh) ts = 57.2°F (initial guess) Solving yields s = 1.184 × 10–3 lbm/ft3. The corresponding value of ts = 70.9°F. Repeat the process using ts = 70.9°F as the initial guess. The result is s = 0.977 × 10–3 lbm/ft3 and ts = 64.9°F. Continue iterations until s converges to 1.038 × 10–3 lbm/ft3 and ts = 66.9°F. To solve for the rates of evaporation and heat transfer, first calculate the Reynolds number using air properties at tf = (77 + 66.9)/2 = 72.0°F. uL 0.07435 118,200 0.328 Re L = ------------- = ------------------------------------------------------------------- = 65 ,871 0.04376 where L = 0.328 ft, the length of the plate in the direction of flow. Because ReL < 500,000, flow is laminar over the entire length of the plate; therefore, Sh = 0.664Re 1/2 – h fg = convective heat transfer coefficient = convective mass transfer coefficient = 0.328 × 4.92 × 2 = 3.23 ft2 = surface area (both sides) = evaporation rate = temperature and vapor density at water surface = temperature and vapor density of airstream This energy balance can be rearranged to give – h t – ts = ------ -----------------h fg hM s The heat transfer coefficient h is found by first calculating the Nusselt number: Nu = 0.664Re Nu = 0.037Re 1/2 4/5 Sc 1/3 = 144 Pr Pr 1/3 1/3 for laminar flow for turbulent flow Sh D v h M = ------------ = 7.15 fpm L · m = hM A s The mass transfer coefficient hM requires calculation of Sherwood number Sh, obtained using the analogy expressed in Equations (31) and (32): – = 0.00354 lb/min · q = m h fg = 3.74 Btu/min The same value for q would be obtained by calculating the Nusselt number and heat transfer coefficient h and setting q = hA (t – ts). Fig. 6 Water-Saturated Flat Plate in Flowing Airstream The kind of similarity between heat and mass transfer that results in Equations (29) to (32) can also be shown to exist between heat and momentum transfer. Chilton and Colburn (1934) used this similarity to relate Nusselt number to friction factor by the analogy n f Nu j H = ------------------------ = St Pr = ---1– n 2 Re Pr (33) Fig. 6 Water-Saturated Flat Plate in Flowing Airstream Mass Transfer where n = 2/3, St = Nu/(Re Pr) is the Stanton number, and jH is the Chilton-Colburn j-factor for heat transfer. Substituting Sh for Nu and Sc for Pr in Equations (31) and (32) gives the Chilton-Colburn j-factor for mass transfer, jD: n f Sh j D = ------------------------- = St m Sc = ---1– n 2 Re Sc 6.7 Figure 9. Although the dash-dot line in Figure 9 represents the data, it is actually taken from McAdams (1954) as representative of a large collection of data on heat transfer to single cylinders placed transverse to airstreams. To compare these data with friction, it is necessary to distinguish between total drag and skin friction. Because the analogies are based on skin friction, normal pressure drag must be subtracted from the measured total drag. At Re = 1000, skin friction is 12.6% of the total drag; at Re = 31,600, it is only 1.9%. Consequently, values of f /2 at a high Reynolds number, obtained by the difference, are subject to considerable error. In Figure 10, data on evaporation of water into air for single spheres are presented. The solid line, which best represents these data, agrees with the dashed line representing McAdams’ correlation for heat transfer to spheres. These results cannot be compared with friction or momentum transfer because total drag has not been allocated to skin friction and normal pressure drag. Application of these data to air/water-contacting devices such as air washers and spray cooling towers is well substantiated. (34) where Stm = ShPAm/(Re Sc) is the Stanton number for mass transfer. Equations (33) and (34) are called the Chilton-Colburn j -factor analogy. The power of the Chilton-Colburn j-factor analogy is represented in Figures 7 to 10. Figure 7 plots various experimental values of jD from a flat plate with flow parallel to the plate surface. The solid line, which represents the data to near perfection, is actually f /2 from Blasius’ solution of laminar flow on a flat plate (left-hand portion of the solid line) and Goldstein’s solution for a turbulent boundary layer (right-hand portion). The right-hand part also represents McAdams’ (1954) correlation of turbulent flow heat transfer coefficient for a flat plate. A wetted-wall column is a vertical tube in which a thin liquid film adheres to the tube surface and exchanges mass by evaporation or absorption with a gas flowing through the tube. Figure 8 illustrates typical data on vaporization in wetted-wall columns, plotted as jD versus Re. The point spread with variation in / Dv results from Gilliland’s finding of an exponent of 0.56, not 2/3, representing the effect of the Schmidt number. Gilliland’s equation can be written as follows: j D = 0.023Re – 0.17 Fig. 8 Vaporization and Absorption in Wetted-Wall Column --------Dv – 0.56 (35) Similarly, McAdams’ (1954) equation for heat transfer in pipes can be expressed as j H = 0.023Re – 0.20 cp -------k – 0.7 (36) This is represented by the dash-dot curve in Figure 8, which falls below the mass transfer data. The curve f /2, representing friction in smooth tubes, is the upper, solid curve. Data for liquid evaporation from single cylinders into gas streams flowing transversely to the cylinders’ axes are shown in Fig. 7 Mass Transfer from Flat Plate Fig. 8 Vaporization and Absorption in Wetted-Wall Column Fig. 9 Mass Transfer from Single Cylinders in Crossflow Fig. 7 Mass Transfer from Flat Plate Fig. 9 Mass Transfer from Single Cylinders in Crossflow 6.8 When the temperature of the heat exchanger surface in contact with moist air is below the air’s dew-point temperature, vapor condensation occurs. Typically, air dry-bulb temperature and humidity ratio both decrease as air flows through the exchanger. Therefore, sensible and latent heat transfer occur simultaneously. This process is similar to one that occurs in a spray dehumidifier and can be analyzed using the same procedure; however, this is not generally done. Cooling coil analysis and design are complicated by the problem of determining transport coefficients h, hM, and f. It would be convenient if heat transfer and friction data for dry heating coils could be used with the Colburn analogy to obtain the mass transfer coefficients, but this approach is not always reliable, and Guillory and McQuiston (1973) and Helmer (1974) show that the analogy is not consistently true. Figure 11 shows j-factors for a simple parallelplate exchanger for different surface conditions with sensible heat transfer. Mass transfer j-factors and friction factors exhibit the same behavior. Dry-surface j-factors fall below those obtained under dehumidifying conditions with the surface wet. At low Reynolds numbers, the boundary layer grows quickly; the droplets are soon covered and have little effect on the flow field. As the Reynolds number increases, the boundary layer becomes thin and more of the total flow field is exposed to the droplets. Roughness caused by the droplets induces mixing and larger j-factors. The data in Figure 11 cannot be applied to all surfaces, because the length of the flow channel is also an important variable. However, water collecting on the surface is mainly responsible for breakdown of the j-factor analogy. The j-factor analogy is approximately true when surface conditions are identical. Under some conditions, it is possible to obtain a film of condensate on the surface instead of droplets. Guillory and McQuiston (1973) and Helmer (1974) related dry sensible j- and f-factors to those for wetted dehumidifying surfaces. The equality of jH, jD , and f /2 for certain streamlined shapes at low mass transfer rates has experimental verification. For flow past bluff objects, jH and jD are much smaller than f /2, based on total pressure drag. The heat and mass transfer, however, still relate in a useful way by equating jH and jD . Example 5. Using solid cylinders of volatile solids (e.g., naphthalene, camphor, dichlorobenzene) with airflow normal to these cylinders, Bedingfield and Drew (1950) found that the ratio between the heat and mass transfer coefficients could be closely correlated by the following relation: 0.56 2009 ASHRAE Handbook—Fundamentals Solution: For dry air at 70°F and standard pressure, = 0.075 lbm/ft3, = 0.044 lbm/h·ft, k = 0.0149 Btu/h·ft·°F, and cp = 0.240 Btu/lbm·°F. From Equation (10), Dv = 0.973 ft2/h. Therefore, Reda = Sc = u d/ = 0.0749 31 3600 0.300/(12 0.044) = 4750 Pr = cp /k = 0.240 0.044/0.0149 = 0.709 0.973) = 0.604 / Dv = 0.044/(0.0749 From Figure 9 at Reda = 4750, read jH = 0.0088, and jD = 0.0099. From Equations (33) and (34), h = jH cp u / (Pr) 2/3 = 0.0088 0.0749 0.240 31 3600/(0.709)2/3 3600/(0.604)2/3 = 22.2 Btu/h·ft2 ·°F hM = jD u / (Sc)2/3 = 0.0099 = 1550 ft/h h/ hM = 22.2/(0.0749 1550) = 0.191 Btu/lbm ·°F 31 From the Bedingfield-Drew relation, h/ hM = 0.294(0.604)0.56 = 0.222 Btu/lbm ·°F Equations (34) and (35) are called the Reynolds analogy when Pr = Sc = 1. This suggests that h/ hM = cp = 0.240 Btu/lbm ·°F. This close agreement is because the ratio Sc/Pr is 0.604/0.709 or 0.85, so that the exponent of these numbers has little effect on the ratio of the transfer coefficients. h --------- = 0.294 Btu lb m · F hM ----------Dv The extensive developments for calculating heat transfer coefficients can be applied to calculate mass transfer coefficients under similar geometrical and flow conditions using the j-factor analogy. For example, Table 7 of Chapter 4 lists equations for calculating heat transfer coefficients for flow inside and normal to pipes. Each equation can be used for mass transfer coefficient calculations by equating jH and jD and imposing the same restriction to each stated in Table 7 of Chapter 4. Similarly, mass transfer experiments often replace corresponding heat transfer experiments with complex geometries where exact boundary conditions are difficult to model (Sparrow and Ohadi 1987a, 1987b). The j-factor analogy is useful only at low mass transfer rates. As the rate increases, the movement of matter normal to the transfer surface increases the convective velocity. For example, if a gas is blown from many small holes in a flat plate placed parallel to an airstream, the boundary layer thickens, and resistance to both mass and heat transfer increases with increasing blowing rate. Heat transfer data are usually collected at zero or, at least, insignificant mass transfer rates. Therefore, if such data are to be valid for a mass transfer process, the mass transfer rate (i.e., the blowing) must be low. Fig. 11 Sensible Heat Transfer j-Factors for Parallel Plate Exchanger For completely dry air at 70°F flowing at a velocity of 31 fps over a wet-bulb thermometer of diameter d = 0.300 in., determine the heat and mass transfer coefficients from Figure 9 and compare their ratio with the Bedingfield-Drew relation. Fig. 10 Mass Transfer from Single Spheres Fig. 10 Mass Transfer from Single Spheres Fig. 11 Sensible Heat Transfer j-Factors for Parallel Plate Exchanger Mass Transfer The j-factor relationship jH = jD can still be valid at high mass transfer rates, but neither jH nor jD can be represented by data at zero mass transfer conditions. Chapter 24 of Bird et al. (1960) and Eckert and Drake (1972) have detailed information on high mass transfer rates. 6.9 A mass transfer coefficient is defined using W as the driving potential: · mB = KM Wi – W (40) Lewis Relation Heat and mass transfer coefficients are satisfactorily related at the same Reynolds number by equating the Chilton-Colburn j factors. Comparing Equations (33) and (34) gives St Pr n = f /2 = Stm Scn Inserting the definitions of St, Pr, Stm, and Sc gives cp h ----------- ----------k cp u or Dv h---------------- = P Am --------------------hM cp cp k = P Am Dv 23 23 23 where the coefficient KM is in lbm/h·ft2. For dilute mixtures, Ai A ; that is, the partial mass density of dry air changes by only a small percentage between interface and free stream conditions. Therefore, KM · m B = --------Am Bi – (41) hM P Am = ----------------- --------Dv u 23 where Am = mean density of dry air, lbm/ft3. Comparing Equation (41) with Equation (24) shows that KM h M = --------Am (42) The humid specific heat cpm of the airstream is, by definition (Mason and Monchick 1965), (37) or cpm = ( / A cpm = (1 + W ) cp (43a) The quantity /Dv is the Lewis number Le. Its magnitude expresses relative rates of propagation of energy and mass within a system. It is fairly insensitive to temperature variation. For air and water vapor mixtures, the ratio is (0.60/0.71) or 0.845, and (0.845)2/3 is 0.894. At low diffusion rates, where the heat/mass transfer analogy is valid, PAm is essentially unity. Therefore, for air and water vapor mixtures, h/hM cp 1 (38) ) cp (43b) where cpm is in Btu/lbda ·°F. Substituting from Equations (42) and (43b) into Equation (38) gives h Am --------------------------- = 1 K M A cpm h ---------------K M cpm (44) The ratio of the heat transfer coefficient to the mass transfer coefficient equals the specific heat per unit volume of the mixture at constant pressure. This relation [Equation (38)] is usually called the Lewis relation and is nearly true for air and water vapor at low mass transfer rates. It is generally not true for other gas mixtures because the ratio Le of thermal to vapor diffusivity can differ from unity. Agreement between wet-bulb temperature and adiabatic saturation temperature is a direct result of the nearness of the Lewis number to unity for air and water vapor. The Lewis relation is valid in turbulent flow whether or not /Dv equals 1 because eddy diffusion in turbulent flow involves the same mixing action for heat exchange as for mass exchange, and this action overwhelms any molecular diffusion. Deviations from the Lewis relation are, therefore, due to a laminar boundary layer or a laminar sublayer and buffer zone where molecular transport phenomena are the controlling factors. because Am A because of the small change in dry-air density. Using a mass transfer coefficient with the humidity ratio as the driving force, the Lewis relation becomes ratio of heat to mass transfer coefficient equals humid specific heat. For the plate humidifier illustrated in Figure 6, the total heat transfer from liquid to interface is · q = q A + m B h fg (45) Using the definitions of the heat and mass transfer coefficients gives q = h(ti – t ) + KM (Wi – W )hfg Assuming Equation (44) is valid gives q = KM [cpm(ti – t ) + (Wi – W )hfg ] The enthalpy of the air is approximately h = cpa t + Whs (48) (47) (46) SIMULTANEOUS HEAT AND MASS TRANSFER BETWEEN WATER-WETTED SURFACES AND AIR A simplified method used to solve simultaneous heat and mass transfer problems was developed using the Lewis relation, and it gives satisfactory results for most air-conditioning processes. Extrapolation to very high mass transfer rates, where the simple heat-mass transfer analogy is not valid, leads to erroneous results. The enthalpy hs of the water vapor can be expressed by the ideal gas law as hs = cps (t – to) + hfgo (49) Enthalpy Potential The water vapor concentration in the air is the humidity ratio W, defined as W B ----A where the base of enthalpy is taken as saturated water at temperature to . Choosing to = 32°F to correspond with the base of the dry-air enthalpy gives h = (cpa + Wcps)(t – 32) + Whfgo = cpm(t – 32) + Whfgo (50) (39) 6.10 If small changes in the latent heat of vaporization of water with temperature are neglected when comparing Equations (48) and (50), the total heat transfer can be written as q = KM (hi – h ) (51) 2009 ASHRAE Handbook—Fundamentals AH = aH Acs l and AM = aM Acs l (54) The basic equations for the process occurring in the differential length dl can be written for Mass transfer –dGL = Ga dW = KM aM (Wi – W )dl (55) Where the driving potential for heat transfer is temperature difference and the driving potential for mass transfer is mass concentration or partial pressure, the driving potential for simultaneous transfer of heat and mass in an air water/vapor mixture is, to a close approximation, enthalpy. That is, the water evaporation rate, air moisture content increase, and mass transfer rate are all equal. Heat transfer to air Ga cpm dta = ha aH (ti – ta) dl Total energy transfer to air Ga(cpm dta + hfgodW ) = [KM aM (Wi – W )hfg + haaH (ti – ta)] dl (57) Assuming aH = aM and Le = 1, and neglecting small variations in hfg, Equation (57) reduces to Ga dh = KM aM (hi – h) dl (58) (56) Basic Equations for Direct-Contact Equipment Air-conditioning equipment can be classified by whether the air and water used as a cooling or heating fluid are (1) in direct contact or (2) separated by a solid wall. Examples of the former are air washers and cooling towers; an example of the latter is a directexpansion refrigerant (or water) cooling and dehumidifying coil. In both cases, the airstream is in contact with a water surface. Direct contact implies contact directly with the cooling (or heating) fluid. In the dehumidifying coil, contact is direct with condensate removed from the airstream, but is indirect with refrigerant flowing inside the coil tubes. These two cases are treated separately because the surface areas of direct-contact equipment cannot be evaluated. For the direct-contact spray chamber air washer of crosssectional area Acs and length l (Figure 12), the steady mass flow rate of dry air per unit cross-sectional area is · m a A cs = G a (52) The heat and mass transfer areas of spray chambers are assumed to be identical (aH = aM). Where packing materials, such as wood slats or Raschig rings, are used, the two areas may be considerably different because the packing may not be wet uniformly. The validity of the Lewis relation was discussed previously. It is not necessary to account for small changes in latent heat hfg after making the two previous assumptions. Energy balance Ga dh = ±GLcL dtL (59) and the corresponding mass flux of water flowing parallel with the air is · m L A cs = G L where (53) · m a = mass flow rate of air, lb/h Ga = mass flux or flow rate per unit cross-sectional area for air, lb/h·ft2 · m L = mass flow rate of liquid, lb/h GL = mass flux or flow rate per unit cross-sectional area for liquid, lb/h·ft2 A minus sign refers to parallel flow of air and water; a plus refers to counterflow (water flow in the opposite direction from airflow). The water flow rate changes between inlet and outlet as a result of the mass transfer. For exact energy balance, the term (cLtLdGL) should be added to the right side of Equation (59). The percentage change in GL is quite small in usual applications of air-conditioning equipment and, therefore, can be ignored. Heat transfer to water ±GLcL dtL = hLaH (tL – ti) dl (60) Because water is evaporating or condensing, GL changes by an amount dGL in a differential length dl of the chamber. Similar changes occur in temperature, humidity ratio, enthalpy, and other properties. Because evaluating the true surface area in direct-contact equipment is difficult, it is common to work on a unit volume basis. If aH and aM are the areas of heat transfer and mass transfer surface per unit of chamber volume, respectively, the total surface areas for heat and mass transfer are Equations (55) to (60) are the basic relations for solution of simultaneous heat and mass transfer processes in direct-contact airconditioning equipment. To facilitate use of these relations in equipment design or performance, three other equations can be extracted from the above set. Combining Equations (58), (59), and (60) gives hL aH hL h–h ------------i = – -------------- = – ------K M aM KM tL – ti (61) Fig. 12 Air Washer Spray Chamber Equation (61) relates the enthalpy potential for total heat transfer through the gas film to the temperature potential for this same transfer through the liquid film. Physical reasoning leads to the conclusion that this ratio is proportional to the ratio of gas film resistance (1/KM) to liquid film resistance (1/hL). Combining Equations (56), (58), and (44) gives h–h dh ------ = ------------i dt a ta – ti Similarly, combining Equations (55), (56), and (44) gives dW = W – W----------------------i ta – ti dt a (63) (62) Fig. 12 Air Washer Spray Chamber Mass Transfer Equation (63) indicates that, at any cross section in the spray chamber, the instantaneous slope of the air path dW/dta on a psychrometric chart is determined by a straight line connecting the air state with the interface saturation state at that cross section. In Figure 13, state 1 represents the state of the air entering the parallel-flow air washer chamber of Figure 12. The washer is operating as a heating and humidifying apparatus, so the interface saturation state of the water at air inlet is the state designated 1i. Therefore, the initial slope of the air path is along a line directed from state 1 to state 1i. As the air is heated, the water cools and the interface temperature drops. Corresponding air states and interface saturation states are indicated by the letters a, b, c, and d in Figure 13. In each instance, the air path is directed toward the associated interface state. The interface states are derived from Equations (59) and (61). Equation (59) describes how air enthalpy changes with water temperature; Equation (61) describes how the interface saturation state changes to accommodate this change in air and water conditions. The solution for the interface state on the normal psychrometric chart of Figure 13 can be determined either by trial and error from Equations (59) and (61) or by a complex graphical procedure (Kusuda 1957). 6.11 Air heat transfer coefficient per cubic foot of chamber volume ha aH = 72 Btu/h·°F·ft3 Liquid heat transfer coefficient per cubic foot of chamber volume hL aH = 900 Btu/h·°F·ft3 Air volumetric flow rate Q = 6500 cfm · Solution: The air mass flow rate m a = 6500 0.075 = 490 lb/min; the · required spray chamber cross-sectional area is then A = m /G = 490 60/1200 = 24.5 ft2. The mass transfer coefficient is given by the Lewis relation [Equation (44)] as KM aM = (ha aH)/cpm = 72/0.24 = 300 lb/h·ft3 Figure 14 shows the enthalpy/temperature psychrometric chart with the graphical solution for the interface states and the air path through the washer spray chamber. 1. Enter bottom of chart with t a1 of 45°F, and follow up to saturation curve to establish air enthalpy h1 of 17.65 Btu/lb. Extend this enthalpy line to intersect initial air temperature ta1 of 65°F (state 1 of air) and initial water temperature tL1 of 95°F at point A. (Note that the temperature scale is used for both air and water temperatures.) 2. Through point A, construct the energy balance line A-B with a slope of cL GL dh ------- = – -------------- = – 0.7 Ga dt L Point B is determined by intersection with the leaving water temperature tL2 = 75°F. The negative slope here is a consequence of the parallel flow, which results in the air/water mixture’s approaching, but not reaching, the common saturation state s. (Line A-B has no physical significance in representing any air state on the psychrometric chart. It is merely a construction line in the graphical solution.) 3. Through point A, construct the tie-line A-1i having a slope of h – hi hL aH ------------ = – -------------- = – 900 = – 3 -------tL – ti KM aM 300 The intersection of this line with the saturation curve gives the initial interface state 1i at the chamber inlet. [Note how the energy balance cs a a Air Washers Air washers are direct-contact apparatus used to (1) simultaneously change the temperature and humidity content of air passing through the chamber and (2) remove air contaminants such as dust and odors. Adiabatic spray washers, which have no external heating or chilling source, are used to cool and humidify air. Chilled-spray air washers have an external chiller to cool and dehumidify air. Heated-spray air washers, with an external heating source that provides additional energy for water evaporation, are used to humidify and possibly heat air. Example 6. A parallel-flow air washer with the following design conditions is to be designed (see Figure 12). Water temperature at inlet tL1 = 95°F Water temperature at outlet tL2 = 75°F Air temperature at inlet ta1 = 65°F Air wet-bulb at inlet t a1 = 45°F Air mass flow rate per unit area Ga = 1200 lb/h·ft2 Spray ratio GL/Ga = 0.70 Fig. 14 Graphical Solution for Air-State Path in Parallel Flow Air Washer Fig. 13 Air Washer Humidification Process on Psychrometric Chart Fig. 13 Air Washer Humidification Process on Psychrometric Chart Fig. 14 Graphical Solution for Air-State Path in Parallel Flow Air Washer 6.12 line and tie-line, representing Equations (59) and (61), combine for a simple graphical solution on Figure 14 for the interface state.] 4. The initial slope of the air path can now be constructed, according to Equation (62), drawing line 1-a toward the initial interface state 1i. (The length of line 1-a depends on the degree of accuracy required in the solution and the rate at which the slope of the air path changes.) 5. Construct the horizontal line a-M, locating point M on the energybalance line. Draw a new tie-line (slope of –3 as before) from M to ai locating interface state ai . Continue the air path from a to b by directing it toward the new interface state ai . (Note that the change in slope of the air path from 1-a to a-b is quite small, justifying the path incremental lengths used.) 6. Continue in the manner of step 5 until point 2, the final state of air leaving the chamber, is reached. In this example, six steps are used in the graphical construction, with the following results: State tL h ti hi ta 1 95.0 17.65 84.5 49.00 65.0 a 91.0 20.45 82.3 46.25 66.8 b 87.0 23.25 80.1 43.80 68.5 c 83.0 26.05 77.8 41.50 70.0 d 79.0 28.85 75.6 39.10 71.4 2 75.0 31.65 73.2 37.00 72.4 2009 ASHRAE Handbook—Fundamentals spray ratio, it is possible, by trial and error, to determine the proper slope of the tie-line necessary to achieve the measured final air state. The tie-line slope gives the ratio hLaH /KM aM ; KM aM is found from the integral relationship in Example 6 from the known chamber length l. Additional descriptions of air spray washers and general performance criteria are given in Chapter 40 of the 2008 ASHRAE Handbook—HVAC Systems and Equipment. Cooling Towers A cooling tower is a direct-contact heat exchanger in which waste heat picked up by the cooling water from a refrigerator, air conditioner, or industrial process is transferred to atmospheric air by cooling the water. Cooling is achieved by breaking up the water flow to provide a large water surface for air, moving by natural or forced convection through the tower, to contact the water. Cooling towers may be counterflow, crossflow, or a combination of both. The temperature of water leaving the tower and the packing depth needed to achieve the desired leaving water temperature are of primary interest for design. Therefore, the mass and energy balance equations are based on an overall coefficient K, which is based on (1) the enthalpy driving force from h at the bulk water temperature and (2) neglecting the film resistance. Combining Equations (58) and (59) and using the parameters described previously yields G L c L dt = K M a M h i – h d l = G a dh K a dV h – h a = ---------------------------------A cs or Ka V --------- = ·m L t2 t1 The final state of air leaving the washer is ta2 = 72.4°F and h2 = 31.65 Btu/lb (wet-bulb temperature ta2 = 67°F). 7. The final step involves calculating the required length of the spray chamber. From Equation (59), Ga l = -------------K M aM 2 dh ----------------1 hi – h (64) The integral is evaluated graphically by plotting 1/(hi – h) versus h, as shown in Figure 15. Any satisfactory graphical method can be used to evaluate the area under the curve. Simpson’s rule with four equal increments of h equal to 3.5 gives N= 2 dh ----------------1 hi – h h3 y1 + 4 y2 + 2 y3 + 4 y4 + y5 0.0553) c L dt -------------------h – ha (65) N = (3.5/3)[0.0319 + (4 0.0400) + (2 + (4 0.0865) + 0.1870] = 0.975 Chapter 39 of the 2008 ASHRAE Handbook—HVAC Systems and Equipment covers cooling tower design in detail. Therefore, the design length is l = (1200/300)(0.975) = 3.9 ft. Cooling and Dehumidifying Coils When water vapor is condensed out of an airstream onto an extended-surface (finned) cooling coil, the simultaneous heat and mass transfer problem can be solved by the same procedure set forth for direct-contact equipment. The basic equations are the same, except that the true surface area of coil A is known and the problem does not have to be solved on a unit volume basis. Therefore, if, in Equations (55), (56), and (58), aM dl or aH dl is replaced by dA/Acs, these equations become the basic heat, mass, and total energy transfer equations for indirect-contact equipment such as dehumidifying coils. The energy balance shown by Equation (59) remains unchanged. The heat transfer from the interface to the refrigerant now encounters the combined resistances of the condensate film (RL = 1/hL); the metal wall and fins, if any (Rm); and the refrigerant film (Rr = A/hr Ar). If this combined resistance is designated as Ri = RL + Rm + Rr = 1/Ui, Equation (60) becomes, for a coil dehumidifier, · m L c L dt L = U i t L – t i d A (plus sign for counterflow, minus sign for parallel flow). The tie-line slope is then h–h ------------i = tL – ti Ui ------KM (67) (66) This method can also be used to predict performance of existing direct-contact equipment and to determine transfer coefficients when performance data from test runs are available. By knowing the water and air temperatures entering and leaving the chamber and the Fig. 15 Graphical Solution of dh/(hi h) Fig. 15 Graphical Solution of dh/(hi – h) Figure 16 illustrates the graphical solution on a psychrometric chart for the air path through a dehumidifying coil with a constant refrigerant temperature. Because the tie-line slope is infinite in this Mass Transfer Fig. 16 Graphical Solution for Air-State Path in Dehumidifying Coil with Constant Refrigerant Temperature h1= 46.2 Btu/lbda, W1 = 0.023 lbH2O /lbda 6.13 specific volume = 14.3 ft3/lbda h2 = 24.8 Btu/lbda, W2 = 0.010 lbH2O /lbda Enthalpy of the condensate from saturated-water temperature table is hw = hf at 57°F = 25.28 Btu/lb Then, · m a = 425 60 14.3 = 1783 lb/h · m w = 1783 0.023 – 0.010 = 23.18 lb/h · q out = 1783 46.2 – 24.8 – 23.18 25.28 = 37,570 Btu/h So, the air conditioner’s heat transfer and moisture removal rates are 36,500 Btu/h and 23.18 lb/h, respectively. SYMBOLS A = surface area, ft2 a = constant, dimensionless; or surface area per unit volume, ft2/ft3 Acs = cross-sectional area, ft2 b = exponent or constant, dimensionless C = molal concentration of solute in solvent, lb mol/ft2 cL = specific heat of liquid, Btu/lb·°F cp = specific heat at constant pressure, Btu/lb·°F cpm = specific heat of moist air at constant pressure, Btu/lbda ·°F d = diameter, ft Dv = diffusion coefficient (mass diffusivity), ft2/h f = Fanning friction factor, dimensionless G = mass flux, flow rate per unit of cross-sectional area, lbm/h·ft2 gc = gravitational constant, ft·lbm/h2 ·lbf h = enthalpy, Btu/lb; or heat transfer coefficient, Btu/h·ft2 ·°F hfg = enthalpy of vaporization, Btu/lbm hM = mass transfer coefficient, ft/h J = diffusive mass flux, lbm/h·ft2 J* = diffusive molar flux, lb mol/h·ft2 jD = Colburn mass transfer group = Sh/(ReSc1/3), dimensionless jH = Colburn heat transfer group = Nu/(RePr1/3), dimensionless k = thermal conductivity, Btu/h·ft·°F KM = mass transfer coefficient, lb/h·ft2 L = characteristic length, ft l = length, ft Le = Lewis number = /Dv , dimensionless M = relative molecular weight, lbm/lb · m = rate of mass transfer, lb/h · m = mass flux, lb/h·ft2 · m * = molar flux, lb mol/h·ft2 Nu = Nusselt number = hL/k, dimensionless p = pressure, atm or psi PAm = logarithmic mean density factor Pr = Prandtl number = cp /k, dimensionless Q = volumetric flow rate, cfm q = rate of heat transfer, Btu/h q = heat flux per unit area, Btu/h·ft2 Re = Reynolds number = u L / , dimensionless Ri = combined thermal resistance, ft2 ·°F·h/Btu RL = thermal resistance of condensate film, ft2 ·°F·h/Btu Rm = thermal resistance across metal wall and fins, ft2 ·°F·h/Btu Rr = thermal resistance of refrigerant film, ft2 ·°F·h/Btu Ru = universal gas constant = 1545 lbf ·ft/lb mol·°R Sc = Schmidt number = / Dv, dimensionless Sh = Sherwood number = hML/Dv , dimensionless St = Stanton number = h c p u , dimensionless Stm = mass transfer Stanton number = h M PAm u , dimensionless T = absolute temperature, °R t = temperature, °F u = velocity in x direction, fpm Ui = overall conductance from refrigerant to air-water interface for dehumidifying coil, Btu/h·ft2 ·°F V = fluid stream velocity, fpm v = velocity in y direction, fpm Fig. 16 Graphical Solution for Air-State Path in Dehumidifying Coil with Constant Refrigerant Temperature case, the energy balance line is vertical. The corresponding interface states and air states are denoted by the same letter symbols, and the solution follows the same procedure as in Example 6. If the problem is to determine the required coil surface area for a given performance, the area is computed by the following relation: · ma A = ------KM 2 1 dh ----------------hi – h (68) This graphical solution on the psychrometric chart automatically determines whether any part of the coil is dry. Thus, in the example illustrated in Figure 16, entering air at state 1 initially encounters an interface saturation state 1i, clearly below its dew-point temperature td1, so the coil immediately becomes wet. Had the graphical technique resulted in an initial interface state above the dew-point temperature of the entering air, the coil would be initially dry. The air would then follow a constant humidity ratio line (the sloping W = constant lines on the chart) until the interface state reached the air dew-point temperature. Mizushina et al. (1959) developed this method not only for water vapor and air, but also for other vapor/gas mixtures. Chapter 22 of the 2008 ASHRAE Handbook—HVAC Systems and Equipment shows another related method, based on ARI Standard 410, of determining air-cooling and dehumidifying coil performance. Example 7. Air enters an air conditioner at 14.696 psia (1 atm), 86°F, and 85% rh at a rate of 425 cfm and leaves as saturated air at 57°F. Condensed moisture is also removed at 57°F. Calculate the heat transfer and moisture removal rate from the air. Solution: Water mass flow is · · mw = ma W1 – W2 and energy or heat transfer rate is · · · q out = m a h 1 – h 2 – m w h w Properties of air both at inlet and exit states can be determined from the psychrometric chart as follows: 6.14 vi W X,Y,Z x,y,z = = = = velocity normal to mass transfer surface for component i, ft/h humidity ratio, lbw /lbda coordinate direction, dimensionless coordinate direction, ft 2009 ASHRAE Handbook—Fundamentals Chilton, T.H. and A.P. Colburn. 1934. Mass transfer (absorption) coefficients—Prediction from data on heat transfer and fluid friction. Industrial and Engineering Chemistry 26 (November):1183. Eckert, E.R.G. and R.M. Drake, Jr. 1972. Analysis of heat and mass transfer. McGraw-Hill, New York. Guillory, J.L. and F.C. McQuiston. 1973. An experimental investigation of air dehumidification in a parallel plate heat exchanger. ASHRAE Transactions 79(2):146. Helmer, W.A. 1974. Condensing water vapor—Airflow in a parallel plate heat exchanger. Ph.D. dissertation, Purdue University, West Lafayette, IN. Hirschfelder, J.O., C.F. Curtiss, and R.B. Bird. 1954. Molecular theory of gases and liquids. John Wiley & Sons, New York. Incropera, F.P. and D.P. DeWitt. 1996. Fundamentals of heat and mass transfer, 4th ed. John Wiley & Sons, New York. Kusuda, T. 1957. Graphical method simplifies determination of aircoil, wetheat-transfer surface temperature. Refrigerating Engineering 65:41. Mason, E.A. and L. Monchick. 1965. Survey of the equation of state and transport properties of moist gases. In Humidity and moisture, vol. 3. Reinhold, New York. McAdams, W.H. 1954. Heat transmission, 3rd ed. McGraw-Hill, New York. Mizushina, T., N. Hashimoto, and M. Nakajima. 1959. Design of cooler condensers for gas-vapour mixtures. Chemical Engineering Science 9:195. Reid, R.C. and T.K. Sherwood. 1966. The properties of gases and liquids: Their estimation and correlation, 2nd ed. McGraw-Hill, New York. Reid, R.C., J.M. Prausnitz, and B.E. Poling. 1987. The properties of gases and liquids, 4th ed. McGraw-Hill, New York. Sherwood, T.K. and R.L. Pigford. 1952. Absorption and extraction. McGraw-Hill, New York. Sparrow, E.M. and M.M. Ohadi. 1987a. Comparison of turbulent thermal entrance regions for pipe flows with developed velocity and velocity developing from a sharp-edged inlet. ASME Transactions, Journal of Heat Transfer 109:1028-1030. Sparrow, E.M. and M.M. Ohadi. 1987b. Numerical and experimental studies of turbulent flow in a tube. Numerical Heat Transfer 11:461-476. Treybal, R.E. 1980. Mass transfer operations, 3rd ed. McGraw-Hill, New York. Greek = = = = = = = = = = = = = thermal diffusivity = k/ cp, ft2/h Lennard-Jones energy parameter eddy mass diffusivity, ft2/h time parameter, dimensionless absolute (dynamic) viscosity, lbm/ft·h permeability, gr·in/h·ft2 ·in. Hg kinematic viscosity, ft2/h mass density or concentration, lbm/ft3 characteristic molecular diameter, nm time shear stress in the x-y coordinate plane, lbf /ft2 mass fraction, lb/lb temperature function in Equation (9) D i D,AB Subscripts A a Am B c da H i L M m min o s w = = = = = = = = = = = = = = = = gas component of binary mixture air property logarithmic mean more dilute gas component of binary mixture critical state dry-air property or air-side transfer quantity heat transfer quantity air/water interface value liquid mass transfer quantity mean value or metal minimum property evaluated at 0°F water vapor property or transport quantity water vapor property of main fluid stream BIBLIOGRAPHY Bennett, C.O. and J.E. Myers. 1982. Momentum, heat and mass transfer, 3rd ed. McGraw-Hill, New York. DeWitt, D.P. and E.L. Cussler. 1984. Diffusion, mass transfer in fluid systems. Cambridge University Press, U.K. Geankopolis, C.J. 1993. Transport processes and unit operations, 3rd ed. Prentice Hall, Englewood Cliffs, NJ. Kays, W.M. and M.E. Crawford. 1993. Convective heat and mass transfer. McGraw-Hill, New York. Mikielviez, J. and A.M.A. Rageb. 1995. Simple theoretical approach to direct-contact condensation on subcooled liquid film. International Journal of Heat and Mass Transfer 38(3):557. Ohadi, M.M. and E.M. Sparrow. 1989. Heat transfer in a straight tube situated downstream of a bend. International Journal of Heat and Mass Transfer 32(2):201-212. Superscripts * = on molar basis – = average value = wet bulb REFERENCES Bedingfield, G.H., Jr. and T.B. Drew. 1950. Analogy between heat transfer and mass transfer—A psychrometric study. Industrial and Engineering Chemistry 42:1164. Bird, R.B., W.E. Stewart, and E.N. Lightfoot. 1960. Transport phenomena. John Wiley & Sons, New York. ...
View Full Document

Ask a homework question - tutors are online