Abstract

A numerical study for heat exchanger for spray-assisted low-temperature desalination system is presented for an existing low-temperature desalination unit at Arusha Technical College. This is aimed at recognizing the effect of mass flow and physical parameters like tube layout (diameter and length) on the overall heat transferred and the pressure drop in the shell-and-tube heat exchanger (STHX), as well as the impact of these parameters on the heat transfer coefficient and the overdesign of the STHX. Also, the study provides a suitable mathematical model for the replacement of the current condensation unit which tends to reduce energy consumption by reducing some of the electrical components in the system. A Math CAD model was developed using the Delaware method to obtain the mentioned parameters. The results show that at 0.8 kg/s flow rate a maximum heat transfer coefficient of 23212 W/m2K is achieved in a minimum diameter of 10 mm within a maximum tube length of 1000 mm heat exchanger and the pressure drop seems to be very low in a range of 0.328-0.957 Pa from all configurations. The configuration with 1000 mm tube length and 10 mm diameter performed well on the mass flow of 0.3 kg/s-0.8 kg/s by providing a suitable overall heat transfer coefficient of 2306-2539 W/m2K, while 12.8 is a maximum overdesign coefficient achieved on 0.8 kg/s mass flow. The study results show the possibility of using STHX instead of the current condensation unit in implementing a proposed system layout with the minimum effect of energy consumption.

1. Introduction

The present society demands the efficient utilization of energy and at the same time reducing environmental impact to contribute to sustainable economic development. The increase of freshwater demand is mainly caused by the development of world industries and human population growth, while potable water is one of the key resources for the existence of human beings [1]. Desalination of brackish or seawater is an important approach to solve the water resource dearth. We may categorize the desalination process as thermal and nonthermal based on energy consumption [2].

Spray-assisted low-temperature phase-change desalination technique is the thermal-based method that makes no boundary upon heat and mass transfer mechanism and eliminates the contribution of mechanical energy input [3]. This technology comprises two main units to reach the freshwater, namely, evaporation unit and condensation unit [4]. Fortunately, solar energy can either be used directly or indirectly in supplying thermal energy to an evaporation unit as used in vaporization of liquid into vapor and cause steam formation [5]. However, the heat exchange process occurs in the condensation unit by exchanging steam temperature with that of cooling media (gas or liquid) and condensing the steam into liquid freshwater [6].

Literature shows some of the experimental and numerical studies done by some researchers concerned with the condensation unit in different desalination systems. Alsadaie and Mujtaba [7] through the multistage flash (MSF) desalination developed a dynamic fouling model of the heat exchanger, reaction and diffusion mechanisms were included to investigate the effect of surface temperature and flow velocity, and the fouling results made a clear agreement with the literature on current studies. The cryogenic refrigeration unit embedded with multieffect desalination unit hybrid system was designed by Ghorbani et al. [8] which analyzed the thermos-economic.

In adsorption desalination, Elsayed et al. [9] used an aluminum fumarate metal-organic framework in experimental observation of cooling water desorption, adsorption, and condensation. Olkis et al. [10] analyzed the performance of an adsorption heat transformer through the usage of benchmarking material silica gel in the production of freshwater together with the cooling process as done by Thu et al. [11] and achieved a performance ratio of 0.60 at the inlet temperature of 80°C. Also, silica gel was proposed by Rezk et al. [12] to be employed in solar-assisted adsorption desalination cooling for the performance maximization; the results show an improvement of 70% increase in specific cooling power, 0.961 coefficient of performance, and 6.9 m3/day/ton specific daily water production.

Another study was done by Woo et al. [13] on a silica gel-based adsorption cooling cum desalination system concerned with the performance effect, operating pressure, and brine concentration. Towards eco-friendliness, David et al. [14] recovered the waste heat through the method of evaporative condenser cum thermal desalination. An experimental study was conducted by Rabhi et al. [15] on the external condenser and pin fin absorber from a modified solar still; the evaluation was based on the performance of water production and the thermal behavior.

Chandrakanth et al. [16] computed the tube-side pressure drop and heat transfer area respective in the optimization of the shell and tube condenser process parameters and geometric used in thermal desalination plants. The optimum pressure drop vs. heat transfer area was concluded from the second law analysis to overcome the trade-off in between. Via stoichiometric process in the spray flash desalination system, a water level model for after condenser was constructed by Matsuda et al. [17]; the results show the proper model that could improve the system behavior.

From the above literature, the flash evaporator has been reported from many theoretical and experimental studies as an essential unit in the low-temperature desalination technology. But there is no specific study that discusses the condensation unit on flash evaporation desalination technology. Moreover, currently, there is an existing low-temperature desalination unit at Arusha Technical College (ATC), which uses more energy at the condensation unit. Thus, this paper focuses on a heat exchanger design for spray-assisted low-temperature desalination, whereby the pressure drop, heat transfer coefficient, overall heat transfer coefficient, oversurface, and overdesign are the main concern. The heat exchanger is designed to reduce the level of energy consumption, allowing for a smooth working environment of the desalination unit.

2. Modification of Existing System Configuration

Figure 1 shows a schematic of the current layout with several electric components, namely, blower, feed pump, filter with a suction fan, cooling water circulation pump, and condenser with a suction fan. These components make a high level of energy intensity. Here, saline water from the source chamber is pumped into the flash chamber via nozzles. The sprayed flash occurs and the small droplets absorb the thermal energy in the flash chamber to form steam. Steam is escaping the chamber by passing between the chamber’s walls, sucked through filter suction fan directly to the filter, then sucked by condenser suction fan and condensed, and finally stored as freshwater in a distillate reservoir. The steam is condensed when the cold water is circulated in the condenser from the water circulator tank via a second feed pump. At the same time, the condenser suction fan sucks the steam from the filter to speed up the process in which much steam is allowed to escape into the atmosphere. Both condenser and filter suction fans increase the amount of energy consumption on the system.

Figure 2 shows a schematic of a proposed layout that remains with the blower and the feed pump only which may alleviate the energy consumption through the implementation of the heat exchanger as a condensation unit. Now, source water from the source chamber is pumped to evacuated tube collector through the heat exchanger then fed to a flash chamber via nozzles. As usual, the sprayed flash occurs and absorbs the thermal energy in the flash chamber to form steam. Steam escapes the chamber by passing at the top, condensed by exchanging heat with cold water from the source at the heat exchanger, and stored as a freshwater. The blower is embedded with a flash chamber to increase temperature and relative humidity [18]. The heat exchanger acts as a condenser and at the same time becomes a preheater for the incoming cold water. The evacuated tube collector (also a preheater) raises further the temperature of the feed water before it is fed into the flash chamber.

3. Heat Exchanger

The heat exchange process takes place between two or more channels at different temperatures when fluids of the channels exchange their temperature by passing through a heat exchanger [19]. The heat exchanger is a piece of equipment that permits heat transfer and often goes together with mass transfer [20]. The equipment has been used in various industries including processing industries, power plants, and many other industries [21]. Among the several heat exchanger models, shell-and-tube heat exchangers (STHX) are the best [22] and widely used [23]. Also, it can be constructed easily in an extensive range of sizes compared to other categories [24]. A model of the shell-and-tube heat exchanger is shown in Figure 3(b).

4. Methodology

Thermal conductivity, density, viscosity, and specific heat capacity of the cold and hot fluid stream were among the factors estimated in the primary stage in the STHX design as can be seen in Table 1, taken from the existing evaporator which appeared in Figure 3(a). Double shell passes were taken into account, and a triangular pitch of 1.25 space between tube to tube centers was considered for the analysis. A single segmental baffle was chosen for ease of maintenance and high thermal features. U-tube (tube sheet) was used to permit a distinction of thermal expansion. Moreover, the outlet temperatures 29 and 28°C ( and ) for steam and water, respectively, were assumed to control the output and easy computation as shown in Figure 3(b).

4.1. Thermal Design

The temperature difference on the tube side or that of shell side was used to define the sensible heat transfer rate of a heat exchanger. The cold and hot fluid stream loads were calculated using equation (1) as described by [25] where and are mass flow rate of steam and water, respectively, is heat transfer rate, and are specific heat capacity of steam and water, respectively, and are outlet and inlet temperature of the water, and and are outlet and inlet temperature of steam, respectively.

Equation (2) is used to calculate the logarithmic mean temperature difference () to use in estimating the true temperature difference in equation (3) by applying a correction factor () in the design of shell and tube exchangers to allow for the departure from the true countercurrent flow:

4.2. Overall Heat Transfer Coefficient Assessment

To attain a preliminary appraisal for the size of the STHX, an approximate value for the overall heat transfer coefficient () is used. According to the working fluid, the range is 1500-4000 W/m2K as stated by [26] where is the heat transfer area, is the heat transfer rate, is the pilot overall heat transfer coefficient, and is the true temperature difference.

4.3. Tube-Side Calculations

From the standard tube layouts, the empirical equation (6) was used to determine bundle diameter. and in the equation are coefficients determined by the tube layout and the number of tube passes from Table 2 for triangular pitch and square pitch, while the number of tubes () was calculated by dividing the total heat transfer area () and the tube outer surface area () as shown in

The heat transfer coefficient () of the hot fluid stream was determined by using equation (7) [28]. It is a function of a tube inner diameter (), tube side Prandtl number (Prt), and tube side Reynolds number (Ret). Also, nondimensional numbers like hot fluid thermal conductivity (), hot fluid viscosity (), and water viscosity () are included in the equation [29].

4.4. Shell-Side Equations

The shell diameter is calculated as follows:

In the above equations, is a shell diameter and is a clearance between the diameter of the shell and a tube bundle diameter , where the values of and are presented in Table 3. and are empirical coefficients that relate to the head type designed in the STHX.

The heat transfer coefficient () on the shell side was determined from equation (10), where is a cold fluid thermal conductivity, is a shell equivalent diameter, is a Colburn factor, Prs is a shell side Prandtl number, and and are cold fluid and water viscosity, respectively.

4.5. Overall Coefficient Calculations

In this, the required overall heat transfer coefficient was calculated using

The clean overall heat transfer coefficient was calculated using

If , the following step was considered.

Then, the design overall heat transfer coefficient was determined through

If , the next step was worthy to be carried out.

Overdesign is computed in equation (14). The final design safety margin is provided by overdesign in which the required fouling compensation is represented. Oversurface deals with exchanger surface area and depends on the wall and film resistances and fouling allowance; it can be obtained through equation (15).

4.6. Pressure Drop Calculations
4.6.1. Tube-Side Pressure Drop

A pressure drop () of the hot fluid stream was computed as shown in equation (16), where is a Darcy friction factor, is the specific gravity in the tube side, is the tube length, is the hot fluid mass velocity, is the number of tube passes, and is the tube internal diameter.

The tube entrance, exit, and return losses can be computed through

The pressure loss () in the tube side nozzles is considered through equation (18), where is a tube side nozzle mass velocity and is the number of tube passes.

Equation (19) was used to evaluate the total pressure drop () across the hot fluid stream on the tube side as presented below:

4.6.2. Shell-Side Pressure Drop

Equation (20) was used to find the initial shell side pressure losses (), where is a specific gravity in the shell side, is the friction factor, is the equivalent diameter, is the number of the baffle, and is a shell diameter.

Also, the pressure loss in the shell side nozzles () was considered through equation (21), where is a shell side nozzle mass velocity and is the number of tube passes.

Then, the net pressure drop () across the cold fluid stream on the shell side should be computed by summing the nozzle pressure loss and initial shell side pressure losses as per

5. Results and Discussion

5.1. Effect on the Pressure Drop and the Heat Transfer Coefficient

Figures 4(a)4(c) show the pressure drop together with heat transfer coefficient variations within the STHX with the length of 600 mm, 800 mm, and 1000 mm for three different tube outer diameters of 10 mm, 12 mm, and 14 mm with 1 mm thick for each, respectively, through a given range of mass flow of 0.1-0.8 kg/s as regards the shell side. From Figures 4(a)4(c), it can be seen that both the heat transfer coefficient and the pressure drop increase proportionally to the mass flow rate by considering the nine configurations of heat exchangers. Also, confirm that the heat transfer coefficient is much greater at any point in the range of mass flow rate compared to the pressure drop. Hence, it is noted that the diameter of 10 mm shows the greatest heat transfer coefficient of 18162, 20881, and 23212 W/m2K and the maximum pressure drop of 0.328, 0.593, and 0.957 Pa, while the diameter of 14 mm shows a minimum pressure drop of 0.117, 0.223, and 0.370 Pa and a lesser heat transfer coefficient of 13265, 15326, and 17107 W/m2K in the length of 600 mm, 800 mm, and 1000 mm, respectively.

Figure 4(d) illustrates the variation of all minimum pressure drop (blue) which appears in the tube with a diameter of 14 mm versus the cold fluid stream mass flow rate in the shell side for the three different tube length exchangers. The results show that the tube configuration with a minimum length (600 mm) has a minimum pressure drop of 0.117 Pa corresponding to the rest of the configurations. Figure 4(d) also shows the results of the variation of the maximum heat transfer coefficient (black) which appears in the tube with a diameter of 10 mm versus cold stream mass flow rate in the shell side for the three different tube length exchangers; thus, the tube configuration with a maximum length of 1000 mm portrays maximum heat transfer coefficient of 23212 W/m2K compared to the rest of the configurations.

5.2. Effect on the Overall Heat Transfer Coefficient

Figures 5(a)5(c)) express the effect of mass flow rate on the overall heat transfer coefficients of design heat transfer coefficient () and clean heat transfer coefficient () within the STHX with the length of 600 mm, 800 mm, and 1000 mm for three different tube outer diameters of 10 mm, 12 mm, and 14 mm, respectively, with 1 mm thick for each. The results for both and coefficients demonstrate that when the cold fluid flows over the whole surface of the shell were varied, and increase constantly at any rate regardless of the difference in the tube configuration throughout the study. On the other hand, is shown to increase twice or more than when there is an increase in mass flow in the shell side. This is due to the significant effect of the fouling factor of the working fluids. However, the increase in mass flow and sustaining the temperatures keep and rising proportionally; the design coefficient is used rather than the clean coefficient based on the interest of overdesign. The highest and smallest values of appeared in 10 and 14 mm tube diameter, respectively.

Figure 5(d) illustrates which appears through all-tube configurations. The result shows that the tube configuration with a maximum length (1000 mm) has maximum compared to the other configurations. Moreover, it can be seen that as the length of the tube increases with a decrease in diameter the maximum is achieved. To provide the required rate of heat transfer, the value of should be greater than or equal to the value of the required coefficient, . The results show that the suitable mass flow of 0.3-0.8 kg/s with a range of 2306-2539 W/m2K can be selected for the proper design of the exchanger as shown in Figure 6(d).

5.3. Effect on the Oversurface and Overdesign

Figures 6(a)6(c) demonstrate the oversurface () and overdesign () changes when designing the STHX with the length of 600 mm, 800 mm, and 1000 mm for three different tube outer diameters of 10 mm, 12 mm, and 14 mm with 1 mm thick for each, respectively, through a given range of 0.1-0.8 kg/s. Since is much smaller than , appeared to increase in the range of 19.2-463.3 as the mass flow increases. were shown to increase slightly from negative value on 0.1-0.3 kg/s to a positive value on 0.4-0.8 kg/s in a range of -33.9-12.8.

Since engineers do recognize that there will be uncertainties on the data provided and that there may be times when the feedstock will not exactly match up to what was originally specified, a certain amount of conservatism will be required just to achieve satisfactory performance despite unforeseen circumstances. Figure 6(d) clarifies the possibility of designing the suitable STHX by displaying the greatest values of which mostly appeared in the smallest tube diameter (10 mm) to all three tube length configurations with a total of nine samples. However, is a negative value at a low mass flow rate of 0.1-0.3 kg/s specifically to the configurations of length 600 mm and 800 mm. Conversely, as mass flow increases from 0.4 to 0.8 kg/s, becomes a positive value. Hereafter, the configuration with tube length 1000 mm performed well even from the mass flow of 0.3 kg/s compared to other configurations.

6. Conclusion

A heat exchanger model for a low-temperature desalination system was established through the influence tube length and diameter, mass flow rate on pressure drop, heat transfer coefficient, overall heat transfer coefficient, oversurface, and overdesign. After making a deep analysis of these heat transfer parameters, the following summarized conclusions were obtained. (i)Both the heat transfer coefficient and the pressure drop increase proportionally to the mass flow rate among all nine configurations of heat exchangers. But the pressure drop is noticed to be very low 0.328-0.957 Pa for all studied configurations that will lower the pumping power(ii)The clean overall heat transfer coefficient increases twice or more to that of the design coefficient when there is an increase in mass flow in the shell side due to the effect of the fouling factor(iii)The maximum design coefficient is achieved by increasing the tube length with a decrease in diameter. Also, the mass flow of 0.3-0.8 kg/s is suitable for the proper design of exchanger in the present study(iv)The configuration with the 1000 mm length and 10 mm diameter which provides high heat transfer is recommended to work with a maximum flow rate of 0.8 kg/s to achieve a maximum heat transfer coefficient of 23212 W/m2K, while 12.8 is a maximum overdesign coefficient achieved on 0.8 kg/s mass flow(v)Therefore, the implementation of the proposed layout will ensure less energy consumption in the system. Hence, an experimental study is recommended to be done to compare the results obtained by the present mathematical model

Nomenclature

:Area (m2)
Cp:Specific heat capacity (kJ/kg)
, :Diameter (m)
:Mass flow rate (kg/s)
:Pass number, number
:Correction factor
:Darcy friction factor
:Mass velocity (kg/s m2)
:Heat transfer coefficient (W/m2K)
:Thermal conductivity (W/m K)
:Heat load (kW)
:Temperature (°C)
:Specific gravity
:Overall heat transfer coefficient (W/m2K)
Pr:Prandtl number
Re:Reynolds number.
Greek Letters
:Viscosity (kg/m s)
:Density (kg/m3)
:Pressure drop (Pa)
:Temperature difference (K).
Subscripts
6:Hot fluid inlet
7:Hot fluid outlet
2:Cold fluid inlet
3:Cold fluid outlet
b:Baffle
c:Clearance, clean
d:Design
e:Equivalent
f:Liquid, cold fluid
t:Tube
g:Vapor, gas, steam
w:Water
i:Initial, inner
r:Return
req:Required
lm:Logarithmic mean
n:Nozzles
m:Material, mean
s:Shell
o:Outer.
Abbreviations
STHX:Shell-and-tube heat exchanger.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

There are no conflicts to declare.

Acknowledgments

The authors would like to acknowledge the Arusha Technical College for allowing performing this study.