Open access peer-reviewed chapter

# Direct Contact Heat and Mass Exchanger for Heating, Cooling, Humidification, and Dehumidification

Written By

Marip Kum Ja, Qian Chen, Muhammad Burhan, Doskhan Ybyraiymkul, Muhammad Wakil Shahzad, Raid Alrowais and Kim Choon Ng

Submitted: December 14th, 2021 Reviewed: December 22nd, 2021 Published: February 20th, 2022

DOI: 10.5772/intechopen.102353

From the Edited Volume

## Heat Exchangers

Edited by Laura Castro Gómez, Víctor Manuel Velázquez Flores and Miriam Navarrete Procopio

Chapter metrics overview

View Full Metrics

## Abstract

A direct-contact heat and mass exchanger (DCHME) has many advantages over a traditional surface-type heat exchanger, including a high heat transfer coefficient, simplicity of design, and low OPEX and CAPEX. DCHME has a capability to exchange of both heat and mass between the two fluids in the same process. Hence, DCHMEs are widely used in numerous applications in various industries, including the air conditioning industry for cooling and dehumidification and heating and humidification. Based on their structure, DCHME can be categorized into two groups; two fluids direct contact (TFDC) exchanger and two direct contacts with one non-contact fluid (TDCONF) exchanger. This study developed a mathematical model for these two types of exchangers by using a discretized volume with distributed lumped-parameters method instead of using the conventional log mean enthalpy difference (LMHD) and NTU-effectiveness method. Thus, this model can reflect both heat and mass transfer behavior in every spatially distributed physical system. The objective of this study is to develop a mathematical model to be used as a tool for designing DCHME and to be applied as a sub-function of the model predictive control system to predict the effectiveness and dependent parameters of DCHME under the different load conditions and its various input parameters.

### Keywords

• heat and mass exchanger
• distributed lump model
• humidification and dehumidification

## 1. Introduction

For more than 100 years ago, direct contact heat and mass exchangers (DCHME) have been widely used in various industries, including chemical process plant, food, and beverage industry, geothermal heat recovery, seawater desalination, waste heat recovery, energy storage systems, production of steam generation for the Rankine power cycle, air conditioning and refrigeration industries, and many so forth. DCHME is a device in which the two process streams are flowing and contacting each other to exchange heat and mass between these two streams, which can be gas-solid, gas-liquid, liquid-liquid, liquid-solid, or solid-solid streams. The limitation of DCHME is the contamination of the streams depending on the degree of miscibility. Although it has a limitation, there are many advantages such as no corrode or foul or no degradation of heat transfer performance due to the lack of surfaces, a larger heat transfer surface area, much lower flow resistance compared with surface-type heat exchangers, and less capital and operational cost [1]. Thus, DCHMEs are widely used in air conditioning industry for cooling, heating and humidification, cooling and dehumidification, and cooling and humidification, such devices are swamp cooler or direct contact evaporative cooler, cooling tower, air washer spray chamber, cooling coil of air handling unit (AHU), direct expansion (DX) evaporator coil, indirect evaporative cooler, and M-cycle dew point evaporative cooler. Before developing the numerical model of these devices, the process of each device has required to analysis with the psychometric chart. There are total of eight basic air conditioning processes, which are plotted on a psychometric chart as shown in Figure 1. These processes are the air stream from the initial state O to the state (1) for sensible cooling, (2) for cooling and humidification, (3) for humidification only, (4) for heating and humidification, (5) for sensible heating, (6) for heating and dehumidification, (7) for dehumidification only, and (8) for cooling and dehumidification. In order to achieve the above eight processes, some processes need to transfer heat (sensible heat) only but some need to exchange heat and mass (sensible and latent heat) from the air stream. Process 1, cooling, and process 5, heating, are pure heat transfer processes only, which means removing heat from the air stream for cooling and adding heat to the air stream for heating without changing the moisture of the airflow. Process 3, humidification, and process 7, dehumidification, are pure mass transfer processes without variation in air temperature. The rest processes 2 (cooling and humidification), process 4 (heating and humidification), process 6 (heating and dehumidification), and process 8 (cooling and dehumidification) are both heat and mass (sensible and latent heat) transfer process that is adding or removing of heat and moisture to or from the air stream. For DCHME has a capability to transfer both heat and mass in the same process, almost all the processes except six can be accomplished by using DCHME in a single stage. The typical five DCHMEs mentioned above extensively used in the air conditioning industry are discussed with their basic air conditioning process.

(1) Swamp cooler or direct contact evaporative cooler, shown in Figure 2(a), is a device in which the outside air flows through the medium of cooling pad or fill that is wetted with water by dripping or spraying. The air becomes cool and humidified (process 2) after passing through the medium that results from the absorption of heat from the air by the evaporation of water droplets or water film of the cool pad. For it is a passive cooling device, it has a very efficient electrical energy consumption, and so that, it is a very efficient space cooling in the hot and dry region with the lowest initial cost and operation cost.

(2) Cooling tower, shown in Figure 2(b), is not an air conditioner device, but it is used extensively as a heat sink or heat removal from the condenser of the refrigeration circuit in the air conditioning industry. Likewise, in a swamp cooler, the hot water, carrying the rejection heat from the condenser, can be cool down by evaporating its water droplet or water film from the fill. In the cooling tower, the hot water is cooling down as the air stream is heated up and humidified (process 2) by transferring heat and mass between these two fluids.

(3) Air washer spray chamber, shown in Figure 2(c), is a typical DCHME, which is widely used as an air conditioner and works in synergy with the heating or cooling coil in the air handling unit (AHU). The air washer chamber has the capability to condition the inlet air into three desired outlet air state-points by controlling the spraying water temperature. Heating and humidification, like process 4, can be achieved by spraying hot water into the chamber if its spray temperature is higher than the inlet air dry-bulb temperature. Cooling and dehumidification, like process 8, can be achieved by spraying chilled water which temperature is lower than the dew point temperature of the inlet air of the chamber. The inlet air can also be cooled as per process 2 line without changing the moisture content of air by spraying the chilled water temperature, which is lower than the inlet air dry-bulb temperature but higher than its dew point temperature.

Different structures but similar to another type of DCHME are plate finned tube heat exchanger used as a direct expansion (DX) evaporator coil, cooling coil operating with chilled water, or non-volatile refrigerant, shown in Figure 3(a). Indirect evaporative cooler, shown in Figure 3(b), and M-cycle dew point evaporative cooler, shown in Figure 3(c), is also another type of DCHME. These types of exchangers have three fluid streams in which two streams directly contact each other, and the last stream is separated by the copper tube wall or separator sheet, shown in Figure 3.

(4) Plate finned tube heat exchangers are the most common type of heat exchanger used extensively in the air conditioning industry using as a DX evaporator coil and cooling coil operating with chilled water or ethylene/propylene glycol. For this case, if the dew point temperature of inlet air is lower than the surface temperature of the exchanger, the air stream’s water vapor starts turning into condensate as a form of droplets or film on its surface. Therefore, depending on the dew point and surface temperature, the exchanger, can be thoroughly wetted, fully dry, or partially wetted. The fully dry or dry part of a partially wetted exchanger cannot be presumed as a heat and mass exchanger but as a conventional heat exchanger that can be calculated by its traditional LMDT or NTU-effectiveness method. However, the fully wetted or wet part of the partially wetted exchanger is a direct contact heat and mass exchanger because heat and mass are exchanged between the two fluids of the air stream and the film of water condensate [2], shown in Figure 3(a). The process rendering in the wet region of this exchanger is cooling and dehumidification of the air stream as process 8.

(5) Indirect evaporative cooler and M-cycle dew point evaporative cooler, shown in Figure 3(b) and (c), are a combination of a direct evaporative cooler and air to the air heat exchanger. These coolers are composed of multiple duo wet and dry channels separated by a thin film sheet. The working air flowing through the wet channel is directly in contact with a thin water film of wick material or water mist sprayed into the wet channel for the evaporation process. The required heat for the evaporation is transferred from the product air flowing in the dry channel through the separator sheet. Thus, the product air temperature is cooling down along process 1 without changing the moisture content. For the indirect evaporative cooler, the working air at the beginning of the wet channel will be cooled and humidified like process 2 due to evaporation, but later part of the wet channel, the product air is heated and humidified due to the heat that transfers from the product air, like process 4. Likewise, indirect evaporative cooler, M-cycle dew point evaporative cooler has a similar process, but there has a concise process 2 at the beginning of the wet channel because the working air is partially taken from the product air [3, 4, 5].

Based on the number of fluids flowing in the exchanger and their arrangement, DCHME for air conditioning can be categorized into two main groups, which are shown in Figure 4. In the first group, the two fluids streams, water, and air are flowing in parallel or counter flow direction and directly contacting each other to exchange heat and mass between these two fluids. This type of exchanger is noted as Type-1 two fluids direct contact (TFDC) heat and mass exchanger. Examples of these exchangers are: (1) air washer chamber, (2) cooling tower, and (3) swamp cooler or direct contact evaporative cooler, shown in Figure 2. A single discretized element of this type of TFDC exchanger is shown in Figure 5.

The second group of DCHME has three kinds of working fluids in which two fluids, water, and air, are directly contacting each other to exchange heat and mass between these two fluids. The third fluid is separated by a tube wall or a thin sheet layer to prevent mass transfer but only allow heat transfer. Examples of this type of heat and mass exchanger are shown in Figure 3, and it can be noted as Type-2, two direct contacts with one non-contact fluid (TDCONF) heat and mass exchanger. This TDCONF exchanger can also be divided into two groups based on the separator type. The first group, Type-2.1 has a separator with an extended surface. An example of an extended surface TDCONF exchanger is the wet region of the plate finned-tube cooling coil unit. Based on the nature of the working fluid, Type-2.1 can be sub-categorized into two groups. The first group, Type-2.1.1 of extended surface TDCONF exchanger, works with non-volatile refrigerants such as chilled water, ethylene/propylene glycol, etc., so that the working fluid temperature is changing continuously along with the coil rows depth. The second group, Type-2.1.2, works with volatile refrigerants such as R134a, R410, etc. Due to the evaporation in the working fluid flow, most of the exchanger has the same surface and fluid temperature profile. Although these two subgroups are the same exchanger, their mathematical models are different due to their different working fluid temperature profiles. Another group, Type-2.2 non-extended surface TDCONF exchangers, are using a plain sheet or tube as a separator, and examples of this type of exchangers are in the direct evaporative cooler and M-cycle dew point evaporative cooler. The mathematical model of this exchanger is the same with Type-2.1.1 of extended surface TDCONF exchanger except for the separator’s area calculation.

Kays and London [6] introduced the definitions of effectiveness and NTU method to use in heat exchanger design in their 1955 publication. London et al. [7] used this method to fit the experimental data of the cooling tower, but this method is not generally consistent with all the other units [8]. The reason for the inconsistency is that the method was based only on the sensible heat transfer process, yet the mass exchanging process is excluded. Berman [9] introduced the log-mean enthalpy method (LMED) to reflect both the heat and mass transfer process of cooling tower design. Several studies used this method to analyze heat and mass exchangers, such as cooling towers and spray chambers [9, 10, 11]. LMED method can calculate the dependent parameters and effectiveness of the exchanger without using iteration process, but this method has some limitations to calculate all the parameters of the spatially distributed system. Thus, discretized volume with distributed lumped-parameters method is an alternative approach to develop the exchanger model with element nodes. This method is widely used to study the behavior of spatially distributed parameters of electrical systems [12], chemical reaction systems [13], heat and mass transfer [14, 15], microwave [16, 17], acoustics [18], and so many other systems. In this study, discretized volume with distributed lumped-parameters model was developed for both types of exchangers, Type-1 TFDC and Type-2 TDCONF, based on the principle of the graphical method mentioned in ASHRAE Fundamentals Handbook [19] and Systems and Equipment Handbook [20]. The derivation of the model for Type-1 TFDC is explained in Section 3, and for Type-2 TDCONF exchanger is described in Section 4. The relation of convective heat transfer with mass transfer coefficients between air and water film, which are related to both models, is explained in Section 2.

## 2. The relation between convective heat and mass transfer coefficient

All types of DCHME mentioned above can be assumed that each exchanger has similar geometry and boundary conditions for heat and mass transfer process so convective mass transfer is analogous to convective heat transfer, and that can be applied for both laminar and turbulent flows [19]. Bird et al. [21] and Incropera and DeWitt [22] defined the analogy Nusselt (Nu) and Sherwood (Sh) number for the calculation of mass transfer coefficient hM[m/s] from the heat transfer coefficient αa[W/m2.K]. Nusselt (Nu) number, related with αa, is a function of Reynold (Re) and Pantanal (Pr) number, expressed in Eq. (1). Sherwood (Sh) number, related with hM, is a function of Reynold (Re) and Schmidt (Sc) number, stated in Eq. (2).

Sh=hMDHDv=aRebSc13E2

where aand bare coefficient numbers, DH[m] is characteristic length of exchanger, ka[W/m.K] is thermal conductivity of air, and Dv[m2/s] is mass diffusivity. Division of Eq. (1) with Eq. (2) are called the Reynold analogy, and it gives the relation of Lewis (Le) number which can be seen in Eqs. (5) and (6).

αahM=PrSc13CpμPrρScμ=Le13LeρCpmE4
αahMρCpm=Le231E5

where μ[N s/m2] is dynamic viscosity, Cpm[J/kg.K] is specific heat capacity of moist air, and ρ[kg/m3] is density of air. For the humidity ratio is taken as a driving force in mass transfer process, mass transfer coefficient should be defied with the term KM = hM ρ[kg/m2.s]. The Lewis relation equation Eq. (5) can be written as follows:

αaKMCpm=Le231E6

Similarly, Chilton and Colburn [23] proposed Chilton-Colburn j-factor analogy using this similarity to relate Nusselt number to friction factor by the analogy. The j-factor analogy has some limitation that is it can be only reliable when the surface conditions are identical. Bedingfield and Drew [24] also proposed the relation equation between heat and mass coefficients. Many others relation equations that can be found in literatures to calculate mass transfer coefficient hMfrom heat transfer coefficient αa, but in this study, Lewis relation will use in this model development due to its simplicity.

## 3. Type-1: two fluids direct contact (TFDC) exchanger

TFDC heat exchangers, such as air washer chamber, cooling tower, and swamp cooler or direct contact evaporative cooler, are the first type of DCHME shown in Figure 2. TFDC has two water and air fluids flowing in parallel or counter, and directly contacting each other to exchange heat and mass between them. Figure 6(a) shows the basic one discretized element of air washer spray chamber and direct contact evaporative cooler. The air stream is directly in contact with the saturated layer of each water droplet or the fill/cool pad, and heat and mass are transferring between them according to the mass transfer coefficient KMwith mass transfer area aMand heat transfer coefficient αawith heat transfer area aH. Figure 6(b) shows the fundamental discretized element of the cooling tower and swamp cooler. Both elements have the same concept of numbering node for three layers parameters; air and its average, water and its average, and saturated layer parameters, shown in Figure 6(c).

### 3.1 Energy balance between two fluids of Type-1 model (energy balance line, EBL)

In TFDC exchanger, the water is sprayed into the chamber along the airflow (parallel flow), shown in Figure 5(a), or against the airflow (counter flow), shown in Figure 5(b). Due to the difficulty of area measurement, TFDC exchangers are discretized in length along the fluids flowing, assuming it has same cross sectional area. Figure 5 shows that the system is divided into nine differential volumes numbered by “j,” and the inlet and outlet of the differential volumes are noted by “n.” Since each differential volumes has the same cross sectional area, the differential volume can be changed with differential length “dL.” and mass flowrates of water and air are needed to be changed to mass flux or flow rate per unit cross-sectional area of the exchanger, ACS[m2], for air Ga[kg/(s·m2)], and water GL[kg/(s·m2)].

According to the conservation of energy, the energy transferred from one fluid is equal to the energy gain of the other fluid. Assuming that heat loss from the system and the change of air or water flow rate due to evaporation or condensation are comparatively small and negligible. Hence, the energy equation of two fluids can be expressed as follows:

where, dha[kJ/kg] and dTl[°C] are the enthalpy and temperature difference between outlet and inlet of one discretized element, and Eq. (7) can be rewritten with element node number for the parallel and counter flow as follows:

han+1hanTln+1Tln=GlCplGaE8
han+1hanTlnTln+1=GlCplGaE9

The generalized energy balance line equation with node number for the counter and parallel flow can be described as follows:

han+1hanTln+1Tln=±GlCplGaE10

A minus sign refers to parallel flow of air and water stream. A plus sign refers to counter flow (water flow is in the opposite direction of airflow).

### 3.2 Heat and mass transfer between two fluids of Type-1 model (tie-line slope, TLS)

According to the conservation of mass between the two fluids, change of water amount in liquid stream, dGl, due to evaporation (if saturated humidity ratio, ws[kg/kg], at the saturated layer is greater than humidity ratio, wa, of air stream) or condensation (if ws < wa) is equal to the amount of humidification or dehumidification of air stream, Gadw, in the spray chamber. The above statement can be expressed by equation as follows:

The negative sign means that the increasing of mass in one fluid stream is equal to the losing mass of other fluid stream. Similarly, the change of mass flow rate, massfrom either stream is equal to the mass flux [kg/m2], that is, mass transfer from the air stream to the saturated liquid surface (or) vise visa per contact surface area along the exchanger, ACT[m2]. It can be written in equation as follows:

mass=KMaMwswaavgdLE12

where, KM[kg/(s·m2)] is mass transfer coefficient, that is, mass transfer rate per contact surface area along the exchanger, ACT[m2]. aM[m2/m3] is a mass transferring contact surface area between air and water per unit volume of spray chamber. L[m] is the chamber length. The sensible heat transfer between air stream and saturated liquid surface can be expressed as follows:

GaCpmdTa=αaaHTsTaavgdLE13

where, αa[W/(m2·K)] is convective heat transfer coefficient of air. aH[m2/m3] is heat transfer contact surface area between air and water per unit volume of spray chamber. Cpm[kJ/(kgda·K)] is a specific heat of moist air at constant pressure. Ts[°C] is saturation temperature of water film. Taavg[°C] is the average air temperature of two adjacent nodes, see in Figure 6(c). The total heat transfer between air stream and saturated liquid surface is equal to the sum of sensible heat, CpmdTa, and latent heat, hfgdw. The total heat transfer can be calculated by combining of mass transfer Eq. (12) multiplying with enthalpy of vaporization hfg[kJ/kg] and the sensible heat transfer Eq. (13). The total heat transfer equation can be expressed as follows:

GaCpmdTa+hfgdw=KMaMwswavghfg+αaaHTsTaavgdLE14

Since CpmdTa+hfgdwis equal with enthalpy change of air dha, and a small change of water vaporization heat with the function of temperature is neglected, the Eq. (14) can be rewritten as follows:

The contact surface areas per system volume, aHand aMfor the heat and mass transferring process in spray chambers are identical (aH= aM), however, this hypothetical is not always correct, especially, if the systems are using the packing materials as an extension of contact surface area which are not fully wetted. Assuming Lewis number with the power of 2/3 is similar with 1, see in Eq. (6), and neglecting of a small variations in hfg. McElgin and Wiley [25] simplified Eq. (15) as follows:

The sensible heat transferring from saturated liquid surface to the working fluid stream is equal to the heat gain of the working fluid, the energy balance equations can be depicted as follows:

GlCpldTl=αlaHTsTlavgdLE17

After substitution of Eqs. (16) and (17) into Eq. (7), and rearranged, the Eq. (7) can be written as follows:

αlaHTsTlavgdL=KMaMhshaavgdLE18
hshaavgTsTlavg=αlaHKMaM=αlCpmαaE19

According to the above Eq. (19), the value of, αlCpmαa, is a slope or gradient line of the rise, the enthalpy incensement in vertical over the run, temperature incensement along the horizontal in T-h graph. In other words, the slope is a ratio of the enthalpy difference between air and saturated surface to the temperature difference between working fluid and saturated surface. This slope, the relation of the temperature and enthalpy, can be denoted as tie-line slope, TLSof the Type-1 TFDC exchanger. The TLS, Eq. (19), can be solved by using graphical method, drawing in water-air T-h graph which is proposed in [2, 19, 20, 26], but this method cannot be applied in the numerical model. In order to solve TLS equation numerically, the two unknown parameters, saturated temperature, Ts, and saturated enthalpy, hs, have to reduce one unknown parameter by substitution of their relation equation that can be obtained by fitting of the saturation line on T-h graph, calculated from the psychometric chart equations depicted in ASHRAE code [19], with the third order polynomial regression equation, Eq. (20).

hs=coe1.Ts3+coe2.Ts2+coe3.Ts+coe4E20

Figure 7 shows the third order polynomial regression line with the saturation point, in red circle legend, calculated by psychometric chart equations, and its four related coefficients values, coe1 to coe4. The regression has a very high accuracy as its R-square value 0.9993706, which covering the range of air temperature from 0 to 50°C.

By the substitution of polynomial regression equation Eq. (20) into Eq. (19), the tie-line slope equation become one unknown parameter equation, and that can be rewritten with node number nand jas follows:

coe1Tsj3+coe2Tsj2+coe3+TLSTsj+coe4han+han+12TLSTln+Tln+12=0E21

Newton Raphson iteration method or any other relevant numerical method can be applied to solve out the unknown root value Tsof each node jof the equation, Eq. (21), by using the known value of Tland haof each node n(node numbering of each element is shown in Figure 6).

### 3.3 Dry bulb temperature of each node for Type-1 model

The last unknown parameter of air temperature, Ta[°C], of each node n, can be calculated from the ratio of sensible heat Eq. (13) to total heat and mass transfer Eq. (16), and it gives as follows:

The value of αaKMCpmcan be taken as 1 because the Lewis number with the power of 2/3 is approximately equal to 1, see in Eq. (6). For the TFDC exchanger, this study assumes that the contact surface area per volume of heat transferring process, aH, is identical with the one of mass transferring process, aM. For numerical calculation, the average parameters, Taavgand haavg, are located between two adjacent n nodes, and noted with jnode, shown in Figure 6. After rearranging of the Eq. (22), it gives with node number nand jas follows:

Tan+1Tan=han+1hanhsjhaavgjTsjTan1+han+1han2hsjhaavgjE23

The outlet air temperature of each node, Tan+1, can be calculated from the known inlet temperature, Tan, air enthalpy, ha, average air enthalpy, haavg,vapor saturated air enthalpy, hs, and vapor saturated air temperature, Tsby using Eq. (23).

### 3.4 Contact length calculation for Type-1 TFDEC exchanger model

Calculation of the contact length or depth of packing of the exchanger, L[m], is a primary interest for the designing of the system to achieve a desired out let condition of air or out let water temperature. The length can be calculated by the integration of heat and mass transfer equation, Eq. (16).

L=GaKMaMj=1j=N11hshaavgdhaE24

The integral equation, Eq. (24) can be solved by using Sampson’s rule or Trapezoidal integration method. The integral equation can be rearranged and rewritten with node number nand jas follows:

L=GaKMaMj=1j=N11hsjhaavgj×han+1hanE25

Change of the air mass flow rate, Ga, due to the transfer of water vapor from the water droplet or film to the air by evaporation or condensation is very minimal compared with its mass flowrate, and it can be taken as negligible for all the practical applications. The contact length or exchanger length L, can be calculated from the known air enthalpy, ha, and saturated enthalpy, hs, of each elements.

## 4. Type-2: two direct contact with one non-contact fluids (TDCONF) exchanger

Section 1, discussed about the DCHME that widely used in air conditioning industry, has categorized the exchangers into two groups. Mathematical model of the first type TFDC heat exchanger has been discussed in Section 3. The development of mathematical model for the second type of TDCONF exchanger will be discussed in detail in this section. The examples of TDCONF exchangers are the wet region of plate finned tube heat exchanger for the direct expansion (DX) evaporator coil, cooling coil unit working with chilled water or non-volatile refrigerant, indirect evaporative cooler, and M-cycle dew point evaporative cooler, which are shown in Figure 3.

A discretized volume with distributed lump parameters model of TDCONF, shown in Figure 8, has three fluids, non-contact fluid, contact fluid, and wetted wick layer or thin water film. The product air of indirect evaporative cooler and volatile refrigerant or chilled water of the wetted plate finned tube cooling coil are the examples of non-contact fluid. These fluids are separated by the separator sheet or tube from thin water film and contact fluid air to prevent the transferring of mass and only to allow the transferring of heat between non-contact fluid and thin water film. Thermal resistance due to the conduction of separator sheet or metal tube is Rmw, and due to the convective heat transfer in the non-contact fluid is Rl. But, both heat and mass are transferring between the thin water film and contact fluid air. Heat and mass transfer circuit between from the symmetrical line of contact fluid air to the non-contact fluid with the coefficient of mass transfer, KM, and thermal resistances are shown in Figure 8(c). Each discretized volume of all TDCONF exchanger is composed with three node points locating in line on the four layers of materials. The first node point, Lavgjis at the symmetrical line of non-contact fluid, sjis at the surface of saturated water film, and aavgjis at the symmetrical line of direct contact fluid air. This model assumed that the temperature difference between outer surface of separator sheet or copper tube and saturated surface of water film is negligible. TDCONF exchanger can be divided into three types of different numerical models: (1) Type-2.1.1 extended surface TDCONF exchanger working with non-volatile working fluid, (2) Type-2.1.2 extended surface TDCONF exchanger working with volatile refrigerant, and (3) Type-2.2 non-extended surface TDCONF exchangers, shown in Figure 4.

### 4.1 Type-2.1: Extended surface TDCONF exchanger

In air conditioning process, plate finned tube heat exchange are widely used for cooling and dehumidification process in fan coil unit (FCU) and air handling unit (AHU). If there is no condensation on the outer surface of finned tube heat exchanger, it is a normal heat exchanger and can be calculated by using conventional heat exchanger equation, UA-LMTD or NTU-effectiveness method. However, if the surface temperature of heat exchanger is lower than the dew point temperature of the direct contact fluid air, the condensation process will be taken place and not only heat but also mass will be transferred between the air flow and the water condensate film. Hence, wet region of cooling heat exchanger can be assumed as a direct contact heat and mass exchanger.

### 4.2 Fin efficiency of extended surface TDCONF heat exchanger

In plate finned tube heat exchanger, fin efficiency is one of the critical parameters for the development of exchanger model, and that can be calculated from the physical dimensions of the heat exchanger. Equations (26)(33) are the equations to calculate the physical parameters of heat exchanger; number of fins Nf[−], coil face or frontal area Aa[m2], external exposed prime surface area Ap[m2], external secondary surface area As[m2], internal surface area Ai[m2], total external surface area Ao[m2], the ratio of external to internal surface area Br [−], and coil core surface parameters Fs[−] from its dimensions, width LW[m], height LH[m], depth LD[m], fin gap fg[m], fin thickness ft[m], outside do[m], inside diameter di[m] of tube, longitudinal tube spacing SL[m], and transverse tube spacing ST[m]. This equations are based on heat exchanger with the continuous plate fin and tubes in stagger arrangement. For the other type of heat exchanger, their related equations are stated comprehensively in AHRI standard 410 [27], and that can be applied in the fin efficiency calculation.

Nf=LWfgE26
Aa=LHLW1000000E27
As=NfLHLD500000do2Nt636688E28
Ap=doLWNtdotfNtNf318344E29
Ai=diNtLW318344E30
Ao=As+ApE31
Br=AiAoE32
Fs=AoAfNrE33

Fins can increase heat transfer from a prime surface of heat exchanger. Fin efficiency can be define with the ratio of the actual heat transferred from the fin to the heat that would be transferred if the entire fin were at its root or base temperature [20]. Fin efficiency equation can be written as follows:

where fis the fin efficiency, Tais the temperature of the surrounding environment, and Tfris the temperature at the fin root. Tfis the temperature along the fin. Fin efficiency decreases as the heat transfer coefficient increases because of the increased heat flow. Total heat transfer from a finned tube heat exchanger is sum of heat transfer from finned surface or secondary area Asand un-finned or prime area Ap. It can write in equation as follows:

Q=αpAp+fαsAsTfrTaE35

Assuming the heat transfer coefficients for the finned surface αsand prime surface αpare equal and note as αaair side heat transfer coefficient, Eq. (35) can be rearranged as follows:

s=11fAsAo=Ap+fAsAoE36

where Aois total surface area of (As+Ap). sis fin effectiveness. Schmidt [28] presents empirical equation of fin surface efficiency for circular, rectangular, and hexagonal fins using an equivalent circular fin radius.

f=tanhmproφmproφE37
mp=2αakatfE38
φ=rero11+0.35lnreroE39

where rois the outside tube radius, reis the equivalent circular fin radius, mpis the standard extended surface parameter, φis dimensionless thermal resistance. αais convective heat transfer coefficient of air. kais the thermal conductivity of the fin and tfinis the fin thickness. Plate finned tube heat exchanger for stagger tube array can be considered as the integration of hexagonal fins. For hexagonal fins:

rero=1.27ψβ0.3E40

where Ψand βare defined as; ψ=B/roand β=H/B, βmust be greater than 1. Depending on the ST(the transverse vertical tube spacing) and SL(the longitudinal horizontal tube spacing), shown in Figure 9, Band Hcan be defined as follows:

If SL > ST/2, then B = ST/2. If SL < ST/2, then B = SL.

H=12ST22+SL2E41

However, Schmidt’s empirical equation of fin surface efficiency is limited to the situations of where β > 1.

When calculating the overall area of Type-2.1 extended surface TDCONF exchanger, the fin effectiveness term, s, can used to simplify the thermal resistance equation of heat and mass transfer from air to working fluid.

### 4.3 Type-2.1.1 extended surface TDCONF exchanger operating with non-volatile working fluid

Type-2.1.1 extended surface TDCONF exchanger is widely used as a cooling coil of air conditioning unit that operating with non-volatile working fluid such as chilled water or ethylene/propylene glycol. The air entering to the first part of the cooling coil is start decreasing its temperature but without moisture removing, this is called sensible heat removing part or dry coil region. When the air temperature is decreasing and reaching the dew point temperature of the entering air, moisture in the air stream start condensing on the surface of the coil called sensible and latent heat removing part or wet coil region.

#### 4.3.1 Resistance of heat and mass transfer and dry-wet boundary region

Depending upon the dew point temperature of inlet air and working fluid temperature, the surface of the cooling coil can be fully dry, fully wet, or partially wet. The driving force for heat transfer in dry region is the temperature difference between air temperature Taand working fluid temperature Tl, crossing the thermal resistance Rad, Rmd, and Rl. In the wet region, there are two types of heat and mass transfer processes synergy in series from the air to the working fluid. The first process is the exchange of heat and mass between the air and the saturated surface of water condensate layer due to their potential difference between haand hs. The second process is transferring heat alone from the surface to the working fluid, and their driving force is the temperature difference between Tsand Tl. Figure 10 shows a typical plate finned tube heat exchanger with counter flow process and thermal diagram of dry and wet region with their potential difference and resistances.

The capacity of the sensible heat QDthat transferring in the dry region can be calculated by using the following equation:

QTd=UDAoDΔTm=AoDΔTmRTDE42

where AoD[m2] is total surface area of dry region, RTD[m2.K/W] is total thermal resistance of dry region, and ΔTmis log mean temperature difference between air dry bulb temperature and working fluid temperature. In wet region, total heat capacity,QW, sensible heat and latent heat, transferred from air to saturated surface of condensate water can be calculated by using the following equation:

QTw=UWAoWΔhm=AoWΔhmCpmRaWE43

where AoW[m2] is total surface area of wet region, RaW[m2.K/W] is thermal resistance of air side in wet region, and Δhmis log mean enthalpy difference between air and saturated wetted surface layer. For the calculation of heat transfer from air to working fluid, total overall thermal resistance of finned tube heat exchanger is required to calculate from the heat exchanger parameters and heat transfer coefficient. Overall thermal resistance with clean non-fouled surfaces for dry and wet cooling coils is a combination of three individual thermal resistances: (1) Ra[m2.K/W] convective thermal resistance between air and external surface of the coil, (2) Rm[m2.K/W] total conductance heat transfer of the metal fin, Rf, and tube, Rt, and (3) Rl[m2.K/W] convective thermal resistance between the internal surface of the coil and the working fluid flowing inside the coil. Brown [29] derived thermal resistance of dry region, Rmd,and wet region, Rmw,based on the heat transfer concept, proposed by Ware-Hacha [30], as follows:

Rmw=Rf+Rt=1ssRawCpmm"+BrDiLnDoDi2kt.E45

where m”is dhsdTs. The overall thermal resistance of dry RTdand wet RTwregion can be calculated as follows:

RTw=Raw+Rmw+Rl=1αaw+1ss1αawCpmm"+BrDiLnDoDi2kt+Brαl.E47

where Rad,Raw[m2.K/W] are thermal resistance for convective heat transfer between air and external surface of the dry and wet coil, Rt[m2.K/W] is thermal resistance of the metal tube, and Rl[m2.K/W] is thermal resistance between the internal surface of tube and the working fluid flowing inside the coil. Br [−] is the ratio of total external surface area Aoto internal surface area Ai.

The ratio of wet and dry coil region is influenced by the characteristic of coil such as coil arrangement, refrigerant distribution, coil depth, entering air dew point and dry bulb temperatures, working fluid temperature, and their flowrates. Another parameter, Y[kg °C/kJ], the ratio of fluid temperature rise to air enthalpy drop, is also important parameter to calculate the boundary region. Ycan be driven from the ratio of the heat that gained by working fluid to the heat that lost from the air stream. The equation is as follows:

dTldha=mamlCpl=YE48

The total heat transferring from the air to the saturated water surface in wet region of TDCONF exchanger is similar with TFDC exchanger, Eq. (16). Mass transfer coefficient KM is substituted with Eq. (6), and it can be written as follows:

dqTw=KMhahsdAow=hahsdAowCpmRawE49

Total heat transferred capacity dqTwin wet region also equals with the heat that transfer from the saturated wet surface layer to the working fluid that flowing inside the tube.

dqTw=TsTldAowRmw+RlE50

After equating and rearranging of Eqs. (49) and (50), the following equation will be obtained, and can be denoted as the characteristic of the coil, C[2, 19, 20, 26].

C=Rmw+RlCpmRaw=TsTlhahsE51

For the dry-wet boundary line, the coil characteristic Ccan be expressed as follows:

C=TdpaiTlbhabhsbE52

Enthalpy of air at dry-wet boundary line, hab[kJ/kg] can be derived from Eqs. (48) and (52) as follows:

hab=TdpaiTlo+Yhai+ChsbC+YE53

The exchanger coil can be determined whether it is fully wetted, fully dry, or partially wetted by comparing the values of habwith hai, and hao. If the value habis higher than hai, all the surface of exchanger is fully wetted. If the value of habis higher than haobut less than hai, the surface of exchanger is partially wetted. If the value of habis less than hao, then the surface of exchanger is completely dry. The temperature of air at dry-wet boundary point, Tab[°C], can be calculated from the heat energy equation of the air stream that flowing in dry region.

Tab=TaihaihabCpmE54

The working fluid temperature at dry-wet boundary line, Tlb[°C], cab be calculated from the energy balance equations between the air stream and working fluid in the dry region.

Tlb=TloYCpmTaiTabE55

Air temperature, air enthalpy, and working fluid temperature of dry region can be calculated by using traditional heat exchanger equations. But for the wet region, a discretized volume with distributed lump parameters method, like TFDC model that explained detail in Section 3, can calculate the parameters of air stream, saturated surface, and working fluid.

#### 4.3.2 Energy balance between two fluids of Type-2.1.1 exchanger (energy balance line, EBL)

Energy balance line for the wet region of Type-2.1.1 exchanger can be derived by equating the energy lost from the air stream and the energy gained by the working fluid. EBL equation for each element with the numbering, n, of each node can be expressed as follows:

han+1hanTln+1Tln=±mlCplma=EBLE57

A minus sign refers to a parallel flow between air stream and working fluid flow. A plus sign refers to a counter flow that working fluid is flowing in the opposite direction of airflow.

#### 4.3.3 Heat and mass transfer between two fluids of Type-2.1.1 model (tie-line slope)

Tie-line slope is the ratio of enthalpy difference between air and saturated surface to the temperature difference between working fluid and saturated surface. Tie-line slope, Eq. (59) can be derived from the equating of the sensible heat, transferring from the working fluid to saturated surface, with the total sensible and latent heat that transferring from the saturated surface to the air stream, see in Figure 10. Tie-line slope of Type-2.1.1 extended surface TDCONF exchanger is similar with the one of TFDC exchanger, Eq. (19), but the difference is an addition of thermal resistances of the metal wall and fins in the wet region, Rmw. Hence, Tie-line slope for Type-2.1.1 extended surface TDCONF exchanger can be expressed with the characteristic of coil, C, as follows:

KMhshaavgdA=UslTlavgTsdAE58
hshaavgTsTlavg=UslKM=CpmRawRmw+Rl=1C=TLSE59

The total heat which is transferring from the saturated surface node j, Tsjto the average temperature of air stream, Tlavgjwhich is the average of two adjacent nodes temperature Tlnand similarly from the saturated surface to the air stream are shown in Figure 11. TLSequation, Eq. (59), has two unknown parameters, saturated temperature, Ts, and saturated enthalpy, hs. The third order polynomial regression equation, which is fitted with the saturation 100% RH line, shown in Figure 7, can be used as their relation equation Eq. (20) to the TLSequation Eq. (59), and rearranged with node number as follows:

coe1Tsj3+coe2Tsj2+coe3+TLSTsj+coe4han+han+12TLSTln+Tln+12=0E60

This third order polynomial TLSequation can be solved by using Newton Raphson iteration methods to finds its unknown root parameter, Ts,from the known parameters of tie line slope, enthalpy of air, and working fluid temperature.

#### 4.3.4 Dry bulb temperature of each node for Type-2.1.1 model

Likewise the derivation of air temperature in Section 3.3, the air temperature of each node of TDCONF exchanger can be driven from the ratio of sensible heat that transfer from air to saturated surface to the total heat and mass transfer from the air to saturated surface. The ratio is as follows:

The value of αawKMCpmcan be assumed 1 for the Lewis number to the power of 2/3 is similar with 1, see the derivation in Section 3.2. For numerical calculation, the air temperature of each node, Eq. (61) is needed to be rearranged and rewritten with node number n and j as follows:

Tan+1Tan=han+1hanhsjhavgjTsjTan1+han+1han2hsjhavgjE62

From the above equation outlet air temperature of each node Tan+1can be calculated from the known parameters of air inlet temperature, air enthalpy, and saturated air enthalpy of each node.

#### 4.3.5 Exchanger area calculation of Type-2.1.1

The area of wet region, Aw, can be calculated by integration of the energy balance equation which is the total heat losing from the air stream is equal to the total heat and mass transferring from the air stream to the vapor saturated surface layer.

Mass transfer coefficient, KM, can be substituted with the parameters αawCpm, as the Lewis number to the power of 2/3 is similar with 1, see in Eq. (6). The integral equation Eq. (63) can be solved out with Sampson’s rule or Trapezoidal method, and it can rewrite with node number nand jas follows:

Aow=maCpmαawj=1j=N11hsjhaavgj×han+1hanE64

The wet surface area of TDCONF exchanger can be calculated from the known parameters, air enthalpy haand vapor saturated enthalpy hsof each element.

### 4.4 Type-2.1.2 extended surface TDCONF exchanger operating with volatile working fluid

The exchanger’s geometrical structure of Type-2.1.2 is the same with the one of Type-2.1.1, for both are using plate finned tube heat exchanger having a similar thermal resistance for heat and mass transferring from air to working fluid. But the difference between them is that Type-2.1.1 exchanger is working with non-volatile fluid so that its fluid temperature TLis increasing along the area axis from TLinto TLout, shown in Figure 10, and Type-2.1.2 exchanger is operating with volatile refrigerant so that its working fluid temperature is constant along the coil depth, dTLis 0. As a result, the parameter Ybecome infinity. For coil characteristic Cfor Type-2.1.2 is similar with Type-2.1.1, the margin point, hab, of dry and wet region of Type-2.1.2 exchanger can be calculated from the coil characteristic equation Eq. (65). At boundary point, water saturated surface temperature of dry and wet region is same as dew point temperature of inlet air, Ts=Tdpai. The enthalpy of air at boundary, habcan be expressed as follows:

hab=hsb+TdpaiTlCE65

Similarly with Type-2.1.1, temperature of air at dry-wet boundary point can be calculated from the equation, Eq. (65). As the working fluid temperature is constant in all coil depth, fluid temperature at boundary point TLbis same as with TL. Energy balance line for Type-2.1.2 is vertical for its refrigerant temperature is constant along the coil depth, dTLis 0. As a result, the refrigerant temperature of each node is constant, and the enthalpy of air outlet haocan be calculated from the energy balance equation between the total heat lost from the air stream and the vaporization heat of the refrigerant, and the equation is as follows:

hao=haimlhfgTrma.E67

The equations of tie-line slope, dry bulb air temperature, and contact area for Type-2.1.2 are the same with Eqs. (60), (62), and (64).

### 4.5 Type-2.2: non-extended surface TDCONF heat exchanger

Indirect evaporative cooler and M-cycle dew point evaporative cooler, Type-2.2 exchangers, are passive coolers and very efficient cooler for dry and hot region. The two fluids, a thin water film and working air, are directly contacting each other, and they do not contact with the third fluid, product air, by separating with aluminum or plastic film to prevent the moisture transferring and make sure allowing only heat transfer.

The phenomena of heat and mass transfer and its mathematical model is the same with a wet region of Type-2.1.1 exchanger, except the calculation of separator area. Tube is used as separator in Type-2.1.1 exchanger and flat sheet is used in Type-2.2 exchanger. Thus, total thermal resistance, RT, of non-extended surface TDCONF exchanger can be written with convective heat transfer coefficient of working air, αaw, and of product air, αap. ksepis the thermal conductivity of the separator and tsepis the separator thickness. The total resistance equation can be expressed as follows:

RT=Raw+Rmw+Rap=1αaw+1αawCpmm"+tsepksep+1αapE68

The rest equations of energy balance line (EBL), tie-line slope (TLS), dry bulb air temperature Ta, and contact area for Type-2.2 exchanger are the same with Eqs. (57), (60), (62), and (64).

## 5. Effectiveness of heat and mass exchanger

The effectiveness of heat and mass exchanger are defined in many different ways in many studies. The effectiveness of cooling tower, and direct evaporative cooler based on the temperature is the ratio of the range, the change of water temperature between the inlet, TLin, and the outlet, TLout, to the sum of the range and the approach which is the difference between water outlet temperature and the inlet air wet-bulb temperature, TWBain[31], expressed as follows:

εCoolingTower=RangeRange+Approach=TLinTLoutTLinTWBainE69

Although Jaber and Webb [8] proposed a modified definition of cooling tower effectiveness with enthalpy, temperature based effectiveness, Eq. (69), is widely used in the cooling tower industry and this study also. Saturation effectiveness is a key factor in the determination of direct and indirect evaporative cooler performance. The saturation effectiveness of direct evaporative cooler is the ratio of the dry bulb temperature difference between inlet air, TDBain, and outlet air TDBaout,to the difference between inlet air dry bulb temperature and its wet bulb temperature, TWBain[2], expressed as follows:

εDirectEvapCooler=TDBainTDBaoutTDBainTWBainE70

The saturation effectiveness of indirect evaporative cooler, wet-bulb depression efficiency, is the ratio of the dry bulb temperature difference between the product air inlet, TDBpain, and outlet, TDBpaout,to the difference between the product air inlet and inlet wet bulb temperature of working air, TWBwain[2], expressed as follows:

εIndirectEvapCooler=TDBpainTDBpaoutTDBpainTWBwainE71

The effectiveness of air washer chamber and wet part of cooling coil are based on the enthalpy, and can be defined as the ratio of enthalpy difference between air inlet and out let to the difference between inlet air enthalpy and saturation enthalpy associated with working fluid inlet temperature [8, 32], expressed as follows:

εCoil/Airwasher=hainhaouthainhsWFinE72

## 6. Results and discussion

This study has developed four types of numerical models for DCHME widely used in the air-conditioning industry. The first model explained in Section 3 is for Type-1 TFDC exchanger: air washer chamber, cooling tower, swamp cooler, or direct contact evaporative cooler, shown in Figures 2, 5, and 6. The second model derived in Section 5.2 is for Type-2.1.1 extended surface TDCONF exchanger working with non-volatile refrigerant (2.1.1. Extended Surface-Non Vol). An example of these exchanges is the wet region of plate finned tube heat exchanger used for cooling and dehumidification working with chilled water or ethylene/propylene glycol, shown in Figures 3(a), 8(b) and (c), and 10. The third model derived in Section 5.3 is for Type-2.1.2 extended surface TDCONF exchanger working with volatile refrigerant (2.1.2. Extended Surface-Vol). Examples of these exchangers are the wet region of DX-coil for cooling and dehumidification working with R134a, R410, etc., shown in Figures 3(a) and 8(b) and (c). Finally, the fourth model, derived detail in Section 6, is for Type-2.2 non-extended surface TDCONF exchanger (2.2. Non-Extended Surface), which are indirect evaporative cooler and M-cycle dew point evaporative cooler, shown in Figures 3(b), 8(a) and (c), and 12.

There are two types of problems with known (given) parameters and parameters to calculate for DCHME, which are listed in Table 1. The first type of problem, Case 1a, and 1b, is the case for designing DCHME to solve the length or area of the exchanger to achieve the desired fluid outlet temp or air outlet condition. The second type of problem, Case 2, is the case for the model predictive control system to predict the effectiveness and output parameter of a given DCHME under the load’s variation and different operational parameters. Model 1. TFDC and 2.1.1. Extended Surface-Non Volatile working fluids have three problem cases, but model 2.1.2. Extended Surface-Vol has only two problems excluded Case 1b because there is no temperature difference between the inlet and outlet of the volatile working fluid. Similarly, model Type-2.2. Non-Extended Surface has no Case 1a because the product air operating as a working fluid is the only parameter needed to design and develop the control model of the exchanger. The models with different types of problem cases are shown in Figure 13.

Problem typesKnown parametersParameters to calculate
Case 1a
[Parallel/Counter]
Calculation of length or area of exchanger to get the desired working fluid outlet temp.
1. Both inlets [Tai, wi| TLi]
2. Both flowrates [Ga, GL] TFDC
[ma, mL] TDCONF
3. Working fluid outlet [TLo]
4. Heat transfer coef. [αa aH, αL aH] TFDC
[αa, αL] TDCONF
1. Air outlet condition [Tao, wo]
2. Saturate line [Ts, hs]
3. Length [L] or area [Aw]
Case 1b
[Parallel/Counter]
Calculation of length or area of exchanger to get the desired air outlet condition.
1. Both inlets [Tai,wi| Tli]
2. Both flowrates [Ga, Gl]
3. Air outlet condition [Tao, wb, ha]
4. AW heat transfer coef. [haaH, hlaH]
PFTHEX thermal resistance
1. Working fluid outlet [Tlo]
2. Saturate line [Ts, hs]
3. Length [L] or area [Aw]
Case 2
[Parallel/Counter]
Prediction of the effectiveness and the working fluid outlet temp or air outlet condition of exchanger with known length or area.
1. Both inlets [Tai, wi| Tli]
2. Both flowrates [Ga, Gl]
3. Length [L]
4. AW heat transfer coef. [haaH, hlaH]
PFTHEX thermal resistance
1. Both outlets [Tao, wo| Tlo]
2. Saturate line [Ts, hs]

### Table 1.

List of problem types with known and calculation parameters for DCHME.

Figure 14 shows the simulation result of temperature and enthalpy of the air, water, and saturated layer of each element of Type-1 TFDC exchanger, which are cooling tower or air washer chamber with Case 1 problem, calculation of exchanger length or area. The numbering of nodes iand jfor exchanger is shown in Figure 5, in which water is spraying into an (a) parallel or counter flow with airflow direction. The process is heating and humidifying airflow in the air washer chamber or cooling the water in the cooling tower. Figure 14(a) shows the result of parallel flow exchanger where the energy balance line is inclined to the left side due to its negative slope per Eq. (10), and similarly, the positive slope for counter flow, thus the line inclined to the right side as shown in Figure 14(b). Both counter and parallel flow have the same negative tie line slope value, see Eq. (19), both tie lines are inclined into the left side. The outlet air temperature, green circle, of counter flow is 1.5°C higher than parallel flow under all the same condition. Thus, the counter flow has higher efficiency than the parallel flow.

Figure 15 shows the simulation result of temperature and enthalpy of the air, water, and saturated layer of each element of Type-1 TFDC exchanger, air washer chamber, Type-2.1.1 exchanger, wet region of cooling coil working with non-volatile fluid, for the process of air cooling and dehumidification. The behavior of the energy balance line and tie line is the same as the heating and humidification process, but the waterline, blue circle, is left side of the saturation line because the working fluid temperature is lower than the saturation temperature, red circle, at each node. At the same time, the water line is on the right side of the saturation line because the working fluid temperature is higher than the saturation temperature in the air heating and humidification process, see Figure 14.

Figure 16 shows the effectiveness of the exchanger based on the water stream, Eq. (69), and based on the airflow stream Eq. (72) for cooling tower and air washer chamber under the same giving length with different liquid air ratio and inlet air enthalpy. At a low liquid-air flow ratio, the effectiveness of cooling tower based on water temperature stream is higher because the range, the temperature difference between water inlet and outlet, is higher due to the low water flow rate. However, for the air washer based on air stream, the effectiveness, Eq. (72), is lower than the higher liquid-air flow ratio because the temperature difference between the air inlet and outlet temperature is higher due to the high airflow rate. Therefore, this Type-1 exchanger with the Case 2 problem model can predict the effectiveness and outlet condition of water or air stream under the several of liquid-air flow ratio and inlet air condition. Thus this model can be applied as a sub-function of the model predictive control system.

Figure 17 shows the simulation result of the temperature and enthalpy of each element. Figure 17(a) depicts the result of the model, Type-2.1.1 extended surface working with the non-volatile working fluid, and Figure 17(b) shows the result of the model, Type-2.1.2 extended surface working with volatile refrigerant for the process of cooling and dehumidification. Both models run at 6°C of working fluid temperature with the exact dimension of the plane finned tube heat exchanger and the same flow rate, flow rate ratio, and condition of the inlet air. Depending on the nature of the volatilization of the working fluid, the temperature of the non-volatile working fluid changes from inlet 6°C to outlet 10°C. The outlet air temperature is dropped from 31 to 24.1°C, as shown in Figure 17(a), whereas, in the model of Type-2.1.2, the refrigerant temperature is constant at 6°C, but the outlet air temperature is dropped from 31 to 17°C, as shown in Figure 17(b). Thus, the Type-2.1.2 exchanger has a better performance than the Type-2.1.1 exchanger.

Figure 18 shows the simulation result of model Type-2.1.1 exchanger running with the same plate finned tube heat exchanger and operating with the same volatile refrigerant at 5°C with the same refrigerant-air flow ratio. The result clearly shows that the model can estimate the dry and wet area and predict the outlet air temperature and moisture removing rate under different air inlet conditions (same inlet temperature but different humidity ratio). Thus, these models are suitable for use as a sub-function of the model predictive control system.

## 7. Conclusion

This study has developed a mathematical model based on a discretized volume with distributed lumped-parameters method for two fluid direct contact (TFDC) exchangers and two direct contacts with one non-contact fluid (TDCONF) exchanger. Based on the flow system and structure, this study has developed four models; Type-1 TFDC exchanger model (air washer chamber, cooling tower, and swamp cooler or direct contact evaporative cooler), Type-2.1.1 extended surface TDCONF exchanger model working with non-volatile refrigerant (wet region of plate finned tube heat exchanger cooling coil working with chilled water or ethylene/propylene glycol), Type-2.1.2 extended surface TDCONF exchanger working with volatile refrigerant (wet region of plate finned tube DX-coil), and Type-2.2 non-extended surface TDCONF exchanger (indirect evaporative cooler and M-cycle dew point evaporative cooler). From the simulation result, these models can reflect both heat and mass transfer behavior in every spatially distributed physical system. Moreover, they can predict well the effectiveness and dependent parameters of DCHME under the different load conditions and its various input parameters. Hence, these models can be a valuable tool for designing the exchangers mentioned above and can be applied as a sub-function of the model predictive control system.

## References

1. 1. Jacobs HR. Direct contact heat exchangers. DOI:10.1615/AtoZ.d.direct_contact_heat_exchangers
2. 2. 2020 ASHRAE Handbook—HVAC Systems and Equipment (SI). Atlanta, GA: ASHRAE; 2020
3. 3. Lin J, Wang RZ, Kumja M, Bui TD, Chua KJ. Modelling and experimental investigation of the cross-flow dew point evaporative cooler with and without dehumidification. Applied Thermal Engineering. 2017;121:1-13
4. 4. Lin J, Thu K, Bui TD, Wang RZ, Ng KC, Kumja M, et al. Unsteady-state analysis of a counter-flow dew point evaporative cooling system. Energy. 2016;113:172-185
5. 5. Lin J, Wang RZ, Kumja M, Bui TD, Chua KJ. Multivariate scaling and dimensional analysis of the counter-flow dew point evaporative cooler. Energy Conversion and Management. 2017;150:172-187
6. 6. Kays WM, London AL. Compact Heat Exchangers. New York: McGraw-Hill Companies; 1964
7. 7. London AL, Mason WF, Boelter LMK. Performance characteristics of a mechanically induced draft, counter flow, packed cooling tower. Transactions of ASME. 1940;62:41-50
8. 8. Jaber H and Webb BL. Design of cooling towers by the effectiveness-NTU method. Journal of Heat Transfer. November 1989;111(4):837-843
9. 9. Berman LD. Chapter 2. In: Sawistowski H, editor. Evaporative Cooling of Circulating Water. 2nd ed. New York: Pergamon Press; 1961. pp. 94-99. [translated from Russian by R. Hardbottle]
10. 10. Webb RL. A critical review of cooling tower design methods. In: Shah RK, Subbarao EC, Mashelkar RA, editors. Heal Transfer Equipment Design. Washington, DC: Hemisphere Pub. Corp.; 1988. pp. 547-558
11. 11. Whillier A. A fresh look at the performance of cooling towers. ASHRAE Transactions. 1976;82:269-282
12. 12. Agarwal A, Lang J. Course Materials for 6.002 Circuits and Electronics. MIT OpenCourseWare (PDF). Boston: Massachusetts Institute of Technology; 2007
13. 13. Groppia G, Bellolie A, Tronconip E, Forzatti P. A comparison of lumped and distributed models of monolith catalytic combustors. Chemical Engineering Science. 1995;50(17):2705-2715
14. 14. Ja MK, Koyama S. Numerical and experimental investigations on a two-bed mode adsorption chiller, ICR 2007. Refrigeration creates the future. In: Proceedings of the 22nd IIR International Congress of Refrigeration; 21–26 August 2007; Beijing, China
15. 15. Ja MK, Choo FH, Li B, Chakraborty A, Dass ETM, Zhao K, et al. Two-dimensional numerical analysis of membrane-based heat and mass cross-flow exchanger. Heat Transfer Engineering. 2017;38(4):438-445
16. 16. Kumja M, Ng KC, Yap C, Yanagi H, Koyama S, Saha BB, et al. Modeling of a novel desorption cycle by dielectric heating. Modern Physics Letters B. 2011;23(3):425-428
17. 17. Kum Ja M. Performance of microwave-activated adsorption cycle: Theory and experiments [Ph.D. thesis]. 2010. Available from:http://scholarbank.nus.edu.sg/handle/10635/18640
18. 18. Dimitra Psychogiou, Roberto Gómez-García, Raúl Loeches-Sánchez, Dimitrios Peroulis. Hybrid acoustic-wave-lumped-element resonators (AWLRs) for high-Q bandpass filters with quasi-elliptic frequency response. IEEE Transactions on Microwave Theory and Techniques. JULY 2015;63(7):447-449
19. 19. 2001 ASHRAE Fundamentals Handbook (SI), Section 5.9. Simultaneous Heat and Mass Transfer Between Water-Wetted Surfaces and Air. Atlanta, GA: ASHRAE; 2001
20. 20. ASHRAE Systems and Equipment Handbook (SI). Atlanta, GA: ASHRAE; 2016
21. 21. Bird RB, Stewart WE, Lightfoot EN. Transport Phenomena. New York: John Wiley and Sons; 1960
22. 22. Incropera FP, DeWitt DP. Fundamentals of Heat and Mass Transfer. 4th ed. New York: John Wiley and Sons; 1996
23. 23. Chilton TH, Colburn AP. Mass transfer (absorption) coefficients—Prediction from data on heat transfer and fluid friction. Industrial and Engineering Chemistry. 1934;26:1183
24. 24. Bedingfield GH Jr, Drew TB. Analogy between heat transfer and mass transfer—A psychometric study. Industrial and Engineering Chemistry. 1950;42:1164
25. 25. McElgin J and Wiley DC. Calculation of coil surface areas for air cooling and dehumidification. Heating, Piping, and Air Conditioning. March 1940;12(3):195
26. 26. Kusuda T. Graphical method simplifies determination of air coil, wet heat-transfer surface temperature. Refrigerating Engineering. 1957;65:41
27. 27. 2001 Standard for Forced-Circulation Air-Cooling and Air-Heating Coils, AHRI Standard 410-2001. Arlington, VA: Air-Conditioning, Heating and Refrigeration Institute; 2001
28. 28. Schmidt TE. La Production Calorifique des Surfaces Munies d’ailettes, Annexe, Du bulletin De L’Institut International Du Froid, Annexe G-5. 1945
29. 29. Brown G. Theory of Moist Air Heat Exchangers. vol. 77. Stockholm, Sweden: Transactions of the Royal Institute of Technology; 1945. pp. 12-15
30. 30. Ware CD, Hacha TH. Heat Transfer From Humid Air to Fin and Tube Extended Surface Cooling Coils, ASME Paper No. 60-HT-17. New York, NY, USA: American Society of Mechanical Engineers; 1960
31. 31. Cheremisinoff NP, Cheremisinoff PN. Cooling Towers: Selection, Design, and Practice. Ann Arbor, MI: Ann Arbor Science Publishers; 1981
32. 32. Nellis G, Klein S. Heat Transfer. Cambridge: Cambridge University Press; 2008

Written By

Marip Kum Ja, Qian Chen, Muhammad Burhan, Doskhan Ybyraiymkul, Muhammad Wakil Shahzad, Raid Alrowais and Kim Choon Ng

Submitted: December 14th, 2021 Reviewed: December 22nd, 2021 Published: February 20th, 2022