Next Article in Journal
Degradation of Azo Dye by Ultrasound in Rotating Flow Field
Previous Article in Journal
Capillary Zone Electrophoresis with Light-Emitting Diode-Induced Fluorescence Detection for the Analysis of Monoclonal Antibodies: Detector Optimization through Design of Experiments and Comparison to UV Detection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Study of a Closed-Cycle Evaporation System for Seawater Desalination

1
School of Civil Engineering, Hefei University of Technology, Hefei 230009, China
2
Anhui Tongsu Technology Co., Ltd., Hefei 230000, China
3
Lianyungang Urban Water Works Management Office, Lianyungang 222000, China
4
School of Resources and Environmental Engineering, Hefei University of Technology, Hefei 230009, China
*
Author to whom correspondence should be addressed.
Separations 2023, 10(5), 319; https://doi.org/10.3390/separations10050319
Submission received: 25 April 2023 / Revised: 15 May 2023 / Accepted: 19 May 2023 / Published: 21 May 2023
(This article belongs to the Special Issue Separation Techniques in Sustainable Water Management)

Abstract

:
This paper presents a numerical study of a closed-cycle evaporation system for the desalination of seawater. The system couples the condensing end of a heat pump with a humidifier, where the air is dehumidified in the heat pump evaporator. First, the mechanism of action of the closed-cycle evaporation system was analyzed from the perspective of heat transfer, and the control equations for the heat and mass transfer of the system were investigated. In addition, a mathematical model of the system was developed and validated. The influence of several important parameters of the air and seawater entering the system on the system’s performance under the design conditions was investigated numerically. The parametric analysis showed that the effect of the seawater mass flow rate on the system’s productivity was not significant. As the air mass flow rate increases, the freshwater production rate increases and then decreases. The output ratio (GOR) of the system was estimated and found to be competitive with other reported HDH systems.

1. Introduction

Socio-economic development and increases in population size have led to an increasing shortage of available freshwater [1,2]. As a cost-effective desalination method, seawater desalination is considered as a sustainable freshwater production technology [3,4].
HDH (humidification–dehumidification) systems mainly consist of a heat source, a humidification unit, and a dehumidification unit [5,6]. Conventional HDH systems which use solar energy are limited by solar radiation [7,8]. Heat pump units are increasingly being used in HDH systems. Heat pump units have the ability to function as both a heating source, providing warmth, and a cooling source, providing coolness [9,10,11].
Dul et al. [12] applied a heat pump to an HDH system and a pilot study was carried out which showed a productivity rate of 287.8 L per day. Faegh et al. [13] integrated a heat pump into a humidification–dehumidification system and evaluated the system’s performance, showing that the application of heat pumps to HDH systems is feasible.
Lawal et al. [14] developed a mathematical model of an HDH system coupled with a heat pump in terms of mass and energy balance and performed a theoretical analysis. Saeed et al. [15] developed a numerical model of an HDH system in order to analyze the influence of the main parameters on the production costs of the system. This provides some theoretical reference for the study of the performance of this type of system under different operating conditions. Zhang et al. [16] first developed an HDH system coupled with a heat pump, which was shown to have a GOR of 2.05 and a freshwater production efficiency of 22.26 kg per hour. Based on this, Zhang et al. [17] developed a mathematical analysis model of the system by analyzing the heat and mass transfer of the various components and used the MOPSO algorithm to optimize the hybrid system. Xu et al. [18] coupled solar energy as an auxiliary heat source with a heat-pump-driven HDH system for drinking water desalination. The pilot study found the maximum production efficiency of the system to be 3.45 kg per kWh. In addition, the study found that how the heat pump is connected to the HDH system affects the performance of the coupled system. Lawal [19] et al. analyzed three different HDH cycle layouts from a cost perspective. The results show that the lowest cost for desalinated water is for the HP-HDH air-to-heat system, followed by the HP-HDH water heating cycle, and finally the E-HDH.
The experimental analysis may be limited by the test conditions, thus preventing the analysis of the system performance under various operating conditions and scenarios. Due to cost and time constraints, experiments cannot provide a comprehensive prediction of the system operation law. Mathematical models replace some of the experiments and allow for a practical analysis of the system’s performance. In order to analyze the performance of HDH systems more accurately and to follow the changes in the operating state of the system, a number of scholars have developed different theoretical models of HDH systems from a thermodynamic perspective. He et al. [20] introduced an HDH desalination system featuring a heat pump unit and conducted a performance analysis of the system’s heat and mass transfer under the design conditions. In the system, the refrigerant mainly transfers heat with the seawater, the heated seawater creates heat and mass transfer with the air in the humidifier, and the remaining heat from the seawater is transferred to the refrigerant in the evaporator. In this system, additional dehumidifiers and cooling water are used, and a certain amount of energy is wasted. Kaunga et al. [21] developed a non-linear model for HDH systems, addressing the limitations of conventional models by proposing a non-linear mechanistic model for the humidification–dehumidification–desalination process with higher predictive accuracy. Their analysis showed that an increase in the inlet water temperature had a positive effect on the productivity of the system, while an increase in the inlet water-to-air-flow ratio had a negative effect.
The excellent research work of the aforementioned authors has promoted the development of heat-pump-driven HDH system technology, especially the theoretical research undertaken.
In this paper, a new, compact, closed-cycle HDH desalination system is theoretically investigated by developing a mathematical model. In this system, the humidifier is combined with the condenser of the heat pump, while the dehumidifier is combined with the evaporator of the heat pump, and the system is more compact. In this study, the goal is to examine a new heat-pump-operated HDH system using a mathematical model and assess its theoretical efficiency. To explain the system’s operation, the mass and energy balance equations of each component are applied. The air mass flow and seawater mass flow are the primary factors that affect the system’s effectiveness once the structural parameters are established. The paper analyzes how these parameters impact the system’s performance as well as its internal state. This research is useful for estimating and evaluating the operating conditions of various sorts of heat-pump-operated HDH systems. In addition, the research provides a better way to popularize desalination technology in remote areas where energy and freshwater resources are scarce.

2. Description

The closed-cycle evaporation system presented in this study was established in Hefei, China [22]. The system consists of humidifiers, two dehumidifiers, and a heat pump unit. The process flow of the system is shown in Figure 1. Figure 2 shows a photograph of the system.The closed-cycle evaporation system contains three media circulations during operation.
There are three media circulations:
  • Seawater: Seawater is circulated in the circulation tank and the humidifier. The seawater is sprayed by a circulation pump onto a heat exchange coil inside the humidifier. The water molecules in the seawater absorb heat to form steam, and the remaining seawater participates in the next circulation.
  • Air: The circulation fan is installed between the first dehumidifier and the second dehumidifier. Through the suction of the circulation fan, a large amount of air enters from the bottom of the humidifier and is heated and humidified. The nearly saturated air leaves the top of the humidifier and enters the two dehumidifiers where it is cooled and dehumidified, resulting in a large amount of freshwater.
  • Refrigerant: The refrigerant is circulated through the humidifier, the second dehumidifier, and the compressor. The heat pump’s condensing end is represented by the heat exchanger coil of the humidifier, where the refrigerant is condensed and gives off heat. Meanwhile, the dehumidifier functions as the evaporating end of the heat pump, wherein the refrigerant gathers heat from the air.

3. Numerical Model

3.1. Sub-Model of the Humidifier

In the sub-model of the humidifier, the heat exchanger coil is simplified. The flow direction of refrigerant is assumed to be vertically downward, parallel to the direction of seawater movement, and the circulating air flows vertically upward. Figure 3 shows the heat and mass transfer analysis unit.
The following assumptions were considered:
(1)
Humidifiers are not affected by external temperature and humidity;
(2)
The nature of seawater is stable;
(3)
The model is in a stable process.
The mass conservation equation for the control body is expressed as:
m k + m c + m w + m a ( 1 + ω ) = m k + d m k + m c + d m c + m w + d m w + m a ( 1 + ω + d ω )
The conservation of energy equation of the control body is expressed as:
m k i k + m c i c + m w c p w T w + m a i a = ( m k + d m k ) i k + ( m c + d m c ) i c + ( m w + d m w ) c p w ( T w + d T w ) + m a ( i a + d i a )
The change in the mass of the refrigerant in the humidifier can be expressed as:
d m c = d m k
The energy change process of the refrigerant in the humidifier can be expressed by the following equation:
d m c ( i k i c ) = d Q c
d Q c = K o ( T k T w ) d A w
The heat transfer coefficient for a coil within the heat exchanger can be determined by utilizing the subsequent formula:
K o = 1 1 α c d w d c + δ λ d w d m + 1 α w
The following equation can be used to compute the mass exchange of air within the humidifier:
d m a = d m w = β a w [ ω w s ω a ] d A a
Expressing the relationship between heat transfer and mass transfer, the Lewis factor plays a crucial role as a parameter in the heat–mass transfer process. The empirical equation of the Lewis factor is:
L e f = 0.865 0.667 ( ω w s + 0.622 ω a + 0.622 1 ) / ln ( ω w s + 0.622 ω a + 0.622 )
By organizing the above equations, the control equation for the mathematical model of the humidifier is:
d m w = β a w [ ω w s ω a ] d A a
d T w d A a = β a w [ L e f ( i a s   i a   ) + ( 1 L e f ) ( ω w s ω a ) i v ] m w c p w + K o ( T k T w ) + c p w T w β a w [ ω w s ω a ] m w c p w
d ω a d A a = - β a w [ ω w s ω a ] m a
d i a d A a = β a w [ L e f ( i a s   i a   ) + ( 1 L e f ) ( ω w s ω a ) i v ] m a
d T a d A a = β a w [ L e f ( i a s   i a   ) + ( 1 L e f ) ( ω w s ω a ) i v ] m a c p m a + ( i f g w o + c p v T a ) [ β a w ( ω w s ω a ) ] m a c p m a

3.2. Sub-Model of the Heat Pump Unit

In this study, the heat pump unit heats the water and cools the air directly. Its energy balance equation is given by:
Q c = Q e + P
Q c = m c ( i c 1 i c 2 )
Q e = m c ( i c 4 i c 3 )
When the evaporating and condensing temperatures are determined, the compressor capacity and the required input power can be estimated using a polynomial function provided by the manufacturer. In the equation fitting method, the compressor cooling capacity and heating capacity can be calculated by the following equations:
Q e = c e , 1 + c e , 2 T k + c e , 3 T e + c e , 4 T k 2 + c e , 5 T k T e + c e , 6 T e 2 + c e , 7 T k 3 + c e , 8 T e T k 2 + c e , 9 T k T e 2 + c e , 10 T e 3
Q c = c c , 1 + c c , 2 T k + c c , 3 T e + c c , 4 T k 2 + c c , 5 T k T e + c c , 6 T e 2 + c c , 7 T k 3 + c c , 8 T e T k 2 + c c , 9 T k T e 2 + c c , 10 T e 3
P = c p , 1 + c p , 2 T k + c p , 3 T e + c p , 4 T k 2 + c p , 5 T k T e + c p , 6 T e 2 + c p , 7 T k 3 + c p , 8 T c T k 2 + c p , 9 T k T e 2 + c p , 10 T e 3

3.3. Sub-Model of the First Dehumidifier

The calculations of the model take into account the following assumptions:
  • Heat radiation losses have been ignored;
  • Heat transfer occurs under steady-state conditions;
  • The thermal resistance due to fins is neglected;
  • The efficiency of the fins under humid conditions is the same as the efficiency of the dry fins;
  • The following equation can be used to represent the heat transfer that occurs within the internal heat exchange channels of the first dehumidifier:
Q 1 d , n = K 1 d , n A 1 d , n Δ T 1 d , n
Q 1 d , n = m a ( i a 2 i a 7 ) = m a ( i a 1 i a 4 ) + m 1 d , n i 1 d , n
ω a 1 = ω a 4
b
The heat exchange in the external heat exchange channels of the first-stage dehumidifier can be expressed as:
Q 1 d , w = K 1 d , w A 1 d , w Δ T 1 d , w
Q 1 d , w = m a ( i a 3 i a 7 ) = m a ( i a 6 i a 5 ) + m 1 d , w i 1 d , w
ω a 5 = ω a 6
To simulate the first dehumidifier, the ε-NTU approach was employed, which brought forth the ideas of the effectiveness of heat exchange and NTU, a measurement for the number of heat transfer units. The effectiveness of a heat exchanger is defined as the proportion of heat transferred by the exchanger to the total heat that could potentially be extracted from the thermal fluid operating at a given starting temperature.
ε can be expressed as:
ε = Q Q max = 1 exp [ ( N T U ) 0.2 { exp [ ( N T U ) 0.78 ] 1 } ]
NTU can be expressed as:
N T U = K A W min = K A m a C a
The individual coefficients within the above equation are calculated as:
K = 1 1 α l η l A A l + 1 α r η r
α = S t c p G
S t = j Pr 2 / 3
Pr = c p μ λ 1 d
ln j = 0.103109 ( ln Re ) 2 1.91091 ( ln Re ) + 3.211

3.4. Sub-Model of the Second Dehumidifier

The calculations of the model take into account the following assumptions:
(1)
Heat radiation losses have been ignored;
(2)
The thermal resistance of the fins is ignored;
(3)
The efficiency of the fins under humid conditions is the same as that of the dry fins;
(4)
The evaporation temperature of the refrigerant is constant.
The heat balance equation for a two-stage dehumidifier is:
Q e = m c ( i c o 2 d i c i 2 d ) + m 2 d i 2 d = m a ( i a 3 i a 4 )

3.5. Validation

In order to verify the mathematical model of the closed-cycle evaporation system, the working conditions in the test were selected and calculated by the mathematical model.
Table 1 displays a comparison between the mathematical model calculations and the experimental data. By analyzing the table, it can be seen that the maximum deviation of temperature is 4.61%, the maximum deviation of moisture content is 3.27%, the maximum error of total productivity is 1.03%, the errors are within 5%, and the errors of the model results and test data are acceptable. Therefore, the calculation platform of the closed-cycle evaporation system is reliable and can be used for theoretical analysis.

3.6. Calculation Process

A balance exists between the evaporation of seawater in the humidifier, the freshwater discharge from the first dehumidifier, and the freshwater discharge from the second dehumidifier. The extent of the evaporation of seawater in the humidifier impacts the heat exchange at the condensing end of the heat pump unit. This, in turn, affects the cooling capability at the evaporating end of the heat pump unit, as well as the freshwater discharge and GOR of the entire system. Therefore, when analyzing the system’s performance using a mathematical model, the first step is to analyze and calculate the humidifier model as the core component.
Only one parameter was changed under any working condition, and the other parameters were kept at their reference values. The reference values and variation range of each key parameter are shown in Table 2. The refrigerant is R134a, whose physical parameters can be found through REFPROP 7.0. In the mathematical model calculation, the characteristics of the seawater are assumed to be the same as those of pure water. In comparison to the power consumed by the heat pump unit, the power utilized by the pump and fan is relatively insignificant. As per the recorded data from the tests, the power consumption of the pump and the fan is approximately 10% of the power consumed by the heat pump unit, and this consumption is considered in the mathematical model.
In actual working conditions, it is difficult for the air state at the outlet of the humidifier to reach just a saturated state. However, if the condensing heat of the refrigerant is not ultimately released, the operational stability of the system will be significantly affected. The purpose of the simulation study is to explore the operation of the closed-cycle low temperature evaporation system under different operating parameters; therefore, in this simulation process, only the humidifier outlet air is not saturated (assuming that the relative humidity of the outlet air is 95%).
The system’s design was implemented by considering different values of the mass flow rate of seawater and air for certain operating conditions, as shown in Figure 4. The main reason for studying the influence of these values on the system’s design is that variations in these parameters can affect the heat transfer rate of the heat pump and therefore the heat transfer efficiency of the overall system.

3.7. Evaluation Indicators

The GOR (gained output ratio) was introduced to evaluate the energy consumption of the closed-cycle evaporation system.
G O R = m f h f W

4. Discussion

The closed-cycle evaporation system has fixed parameters for its components, whereas the system’s performance is primarily influenced by the air mass flow rate and the seawater mass flow rate. In order to assess the system’s productivity and GOR, the impact of these parameters was studied. Moreover, changes in these parameters were analyzed with respect to their effect on the internal air state of the system.

4.1. Seawater Mass Flow

By varying the speed of the seawater pump, the seawater mass flow can be adjusted. Figure 5 shows the effect of seawater mass flow on freshwater production and the GOR. During this calculation, the ambient air condition at the inlet of the radiator fan was (Ta5 = 20 °C, RHa5 = 65%).
Figure 5 demonstrates the small effect of variation in the seawater mass flow on the performance of the closed-cycle evaporation system. The seawater film completely covers the heat exchanger coil during the simulation calculations. The variation in the seawater circulation flow rate has a negligible effect on the individual parameters within the mathematical model.
The charge of sparkling water waft in the closed-cycle evaporation system is extensively impacted by the water evaporation manner in the humidifier. Additionally, a detailed analysis was performed to recognize the alteration within the circumstance of the air gift inside the humidifier.
The change in the seawater mass flow mainly affects the state of the water film in the heat exchanger coil. As can be seen in Figure 6 and Figure 7, there is a slight increase in the temperature and enthalpy of the air inside the humidifier. The increase in seawater mass flow contributes, to a certain extent, to the rate at which the air absorbs heat inside the humidifier. As the seawater mass flow rate increases, the temperature and enthalpy of the air at the humidifier inlet increase by 0.8% and 4%, respectively, and the temperature and enthalpy at the humidifier outlet increase by 0.7% and 1.7%, respectively. The state of the air inside the system is stable during the change in the seawater mass flow rate.

4.2. Air Mass Flow

During this calculation, the ambient air condition on the inlet of the radiator fan is (Ta5 = 20 °C, RHa5 = 65%). Figure 8 illustrates the version of the device’s freshwater along with the flow and the GOR for the different air mass flow values. The effect of air mass flow on freshwater manufacturing rates is illustrated in Figure 8. The increase in air mass flow will increase the performance of the warmth and mass switch inside the humidifier. The increase in air mass flow will increase the humidity distinction and enthalpy difference between the circulating air and the floor of the warmth exchanger coil. As the air mass flow continues to increase, the temperature distinction between the seawater and the refrigerant and the temperature difference between the seawater and the circulating air continues to decrease, the warmth and the mass switch efficiency are affected, and the freshwater flow fee and the GOR generally tend to stabilize.
The heat generated by the condensation of the refrigerant in the humidifier is first transferred to the water film wrapped around the outer wall of the tube, and then the heat is transferred to the air. The increase in the air mass flow rate increases the humidity and enthalpy difference between the surface of the water film and the air, accelerating the efficiency of heat mass transfer between the water film and the air. Figure 9 demonstrates that an increase in the air mass flow rate leads to an increasing temperature and enthalpy of the air at the humidifier inlet. Figure 10 shows that the increase in the air mass flow rate reduces the temperature of the air at the humidifier outlet, but the enthalpy of the air increases. The increase in the air mass flow rate not only increases the efficiency of heat and mass transfer between the air inside the humidifier and the water film outside the heat transfer coil but also reduces the time for heat mass transfer between the air and the water film. The change in air enthalpy is mainly influenced by the latent heat exchange between the air and the water film. The increase in air mass flow promotes the evaporation of the water film, resulting in a higher enthalpy of air leaving the humidifier and a higher share occupied by latent heat transfer, which means that the share occupied by sensible heat transfer is reduced. The change in air temperature is mainly related to the sensible heat exchange between the air and the water film, with the sensible heat transfer occupying a lower share of the heat transfer, resulting in a lower temperature for the air leaving the humidifier. Although the air leaving the humidifier has a lower temperature, its enthalpy is higher, and the heat transfer within the first-level dehumidifier occupies only a smaller share of its heat. When the air with higher enthalpy enters the secondary dehumidifier and releases most of its heat, the air leaving the secondary dehumidifier is saturated. Based on the assumption that the cooling capacity of the compressor in the secondary dehumidifier is constant, the air with the higher enthalpy is relatively hot after it has released heat in the secondary dehumidifier. The main driver of heat transfer within the primary dehumidifier is the temperature difference, but the air leaving the humidifier only accounts for a smaller share of the heat released within the secondary dehumidifier; in addition, according to the control equations for the secondary dehumidifier, an increase in the air mass flow rate increases the heat transfer efficiency of the secondary dehumidifier, ultimately leading to an increase in the temperature and enthalpy of the air entering the humidifier.

4.3. Performance Comparison between the Current System and Other Systems

Following analysis, the maximum value of the GOR simulated for the system in this study was 4.1. Some of the reported performance parameters for other desalination systems are listed in Table 3. The proposed system is ahead of most desalination systems in terms of freshwater productivity per unit of time and has very competitive production costs.

5. Study Limitations and Challenges

In this paper, the performance of a closed-cycle evaporation system consisting of a specific compressor and fixed-size humidifiers and dehumidifiers was investigated. In the operation of the model system, the study was limited mainly by considering the constant isentropic efficiency of the compressor and the constant saturation temperature of the refrigerant. In future research, variation in the system’s performance under operating conditions with non-fixed size components should be investigated.

6. Conclusions

This paper presents a numerical simulation study of a closed-cycle evaporation system driven by a heat pump. A mathematical model of the system was developed, and its effectiveness confirms that the mathematical model can be used for the theoretical evaluation of closed-cycle evaporation systems. In the context of the mathematical model, the influence of different parameters on the performance of the closed-cycle evaporation system was analyzed. The productivity and GOR of existing systems were evaluated and compared with other systems.
(1)
Variations in the seawater mass flow rate did not significantly improve system productivity, and there was no significant effect on system performance.
(2)
As the air mass flow rate increased, the overall heat and mass transfer efficiency of the system improved significantly, and the system’s productivity improved significantly.
(3)
The productivity of the system was 852~867 kg/h. The GOR of the system is estimated to be 4.16~4.24.
(4)
The system is competitive when compared with other systems.
(5)
The mathematical model developed in this paper provides theoretical support for the optimization of the system. The system is a closed system, and the field of application should not be limited to desalination. The system can be tested and applied in the field of wastewater treatment.

Author Contributions

Conceptualization, J.L. and J.W.; methodology, Y.Z., J.L. and J.W.; software, Y.S., J.L. and J.W.; validation, Y.S., J.L. and J.W.; formal analysis, Y.S., J.L. and J.W.; investigation, Y.Z., J.L. and J.W.; resources, Y.Z., Y.S., J.L. and J.W.; data curation, Y.Z. and J.L.; writing—original draft preparation, Y.Z., J.L. and J.W.; writing—review and editing, J.L. and J.W.; visualization, Y.S., J.L. and J.W.; supervision, Y.S. and J.W.; project administration, Y.S. and J.W.; funding acquisition, Y.S. and J.W.. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

mcMass flow rate of the refrigerant liquid, kg/s
mkMass flow rate of the refrigerant vapor, kg/s
ikSpecific enthalpy of the refrigerant vapor, J/kg
icSpecific enthalpy of the refrigerant liquid, J/kg
QcCompressor heat exchange, W
QeCompressor cooling capacity, W
K0Heat transfer coefficient of the heat exchanger coils, W/(m2·K)
TkTemperature of the refrigerant in the tube, K
AwSurface area of the heat exchanger coils, m2
αcHeat transfer coefficient of the refrigerant, W/(m2·K)
αwHeat transfer coefficient of the water film, W/(m2·K)
ΒawMass transfer coefficient
dcInner diameter of the heat exchanger coils, m
dwHeat exchanger coils’ outer diameter, m
dmAverage value of the inner–outer diameter of the heat exchanger coils, m
ΔThickness of the heat exchanger coils, m
ΛThermal conductivity of the heat exchanger coils, W/(m·K)
a1State of the air point 1
a2State of the air point 2
a3State of the air point 3
a4State of the air point 4
a5State of the air point 5
a6State of the air point 6
a7State of the air point 7
ic1Enthalpy of the refrigerant entering the humidifier, kJ/kg
ic2Enthalpy of the refrigerant flowing out of the humidifier, kJ/kg
ic3Enthalpy of the refrigerant entering the second dehumidifier, kJ/kg
ic4Enthalpy of refrigerant flowing out of the second dehumidifier, kJ/kg
PPower consumption of the heat pump unit, W
cc,1~cc,10Constants
ce,1~ce,10Constants
cp,1~cp,10Constants
TeEvaporation temperature, K
TkCondensation temperature, K
Q1d,nHeat exchange in the internal heat exchange channel of the first dehumidifier, W
A1d,nHeat transfer area of the heat transfer channel inside the first dehumidifier, m2
∆T1d,nTemperature difference in the heat exchange channels inside the first dehumidifier, K
K1d,nHeat transfer coefficient of the heat transfer channels inside the first-stage dehumidifier, W/(m2·K)
m1d,nMass of freshwater flowing out of the heat exchange channel inside the first dehumidifier, kg
i1d,nEnthalpy of freshwater flowing out of the heat transfer channel inside the first dehumidifier, J/kg
Q1d,wHeat exchange in the external heat exchange channel of the first dehumidifier, W
A1d,wHeat transfer area of the external heat transfer channel of the first dehumidifier, m2
T1d,wTemperature difference between the external heat exchange channels of the first dehumidifier, K
K1d, wAverage heat transfer coefficient of the external heat transfer channel of the first dehumidifier, W/(m2·K)
m1d, wMass of freshwater flowing out of the external heat exchange channel of the first dehumidifier, kg
i1d, wEnthalpy of freshwater flowing out of the external heat exchange channel of the first dehumidifier, J/kg
CSpecific heat of the air
αlHeat transfer coefficient between the cold fluid and the heat exchange surface
αrHeat transfer coefficient between the thermal fluid and the heat exchange surface
ηlCold fluid channel efficiency
ηrThermal fluid path efficiency
AlHeat transfer area of cold fluid channels
ArHeat transfer area of thermal fluid channels
cpSpecific heat at a constant pressure, kJ/(kg ·K)
λ1dThermal conductivity, W/(kg·K)
ΜPower viscosity, kg/(m·s)
GMass flow rate, kg/(m2·s)
JHeat transfer factor
ReReynolds number
m2dThe mass of freshwater flowing out of the second dehumidifier, kg
i2dEnthalpy of freshwater flowing out of the second dehumidifier, J/kg
WPower consumption of the system, W
MfMass flow rate of freshwater, kg

References

  1. Ullah, I.; Rasul, M. Recent Developments in Solar Thermal Desalination Technologies: A Review. Energies 2018, 12, 119. [Google Scholar] [CrossRef]
  2. Dastbaz, A.; Karimi-Sabet, J.; Ahadi, H.; Amini, Y. Preparation and characterization of novel modified PVDF-HFP/GO/ODS composite hollow fiber membrane for Caspian Sea water desalination. Desalination 2017, 424, 62–73. [Google Scholar] [CrossRef]
  3. Chowdhury, S.; Rahaman, M.; Mazumder, M.A.J. Global mapping of seawater desalination research: A bibliometric analysis of research trends from 1980–2022. Qual. Quant. 2023, 1–22. [Google Scholar] [CrossRef]
  4. Liu, J.; Sun, Y.; Yv, S.; Wang, J.; Hu, K. Experimental study on a Closed-Cycle Humidification and Dehumidification System for Treating Wastewater Containing High Concentrations of Inorganic Salts and Organic Matter. Processes 2021, 9, 671. [Google Scholar] [CrossRef]
  5. Yıldırım, C.; Solmuş, İ. A parametric study on a humidification–dehumidification (HDH) desalination unit powered by solar air and water heaters. Energy Convers. Manag. 2014, 86, 568–575. [Google Scholar] [CrossRef]
  6. Kang, H.; Wang, T.; Zheng, H. Comparative analysis of regenerative and air-extraction multi-stage humidification–dehumidification desalination system using pinch technology. Desalination 2016, 385, 158–166. [Google Scholar] [CrossRef]
  7. Zamen, M.; Amidpour, M.; Soufari, S.M. Cost optimization of a solar humidification–dehumidification desalination unit using mathematical programming. Desalination 2009, 239, 92–99. [Google Scholar] [CrossRef]
  8. Sathyamurthy, R.; El-Agouz, S.A.; Nagarajan, P.K.; Subramani, J.; Arunkumar, T.; Mageshbabu, D.; Madhu, B.; Bharathwaaj, R.; Prakash, N. A review of integrating solar collectors to solar still. Renew. Sustain. Energy Rev. 2017, 77, 1069–1097. [Google Scholar] [CrossRef]
  9. Xu, H.; Dai, Y.J. Parameter analysis and optimization of a two-stage solar assisted heat pump desalination system based on humidification-dehumidification process. Sol. Energy 2019, 187, 185–198. [Google Scholar] [CrossRef]
  10. He, W.F.; Han, D.; Ji, C. Investigation on humidification dehumidification desalination system coupled with heat pump. Desalination 2018, 436, 152–160. [Google Scholar] [CrossRef]
  11. Wu, G.; Zheng, H.F.; Wang, F.; Chang, Z.H. Parametric study of a tandem desalination system based on humidification-dehumidification process with 3-stage heat recovery. Appl. Therm. Eng. 2017, 112, 190–200. [Google Scholar] [CrossRef]
  12. Lawal, D.U.; Antar, M.A.; Khalifa, A.; Zubair, S.M.; Al-Sulaiman, F. Expermental investigation of heat pump driven humidification-dehumidification desalination system for water desalination and space conditioning. Desalination 2020, 475, 114199. [Google Scholar] [CrossRef]
  13. Faegh, M.; Shafii, M.B. Thermal performance assessment of an evaporative condenser-based combined heat pump and humidification-dehumidification desalination system. Desalination 2020, 496, 114733. [Google Scholar] [CrossRef]
  14. Lawal, D.; Antar, M.; Khalifa, A.; Zubair, S.; Al-Sulaiman, F. Humidification-dehumidification desalination system operated by a heat pump, published online in Energy Conversion and Management. Energy Convers. Manag. 2018, 161, 128–140. [Google Scholar] [CrossRef]
  15. Saeed, D.; Abhijit, D.; Aliakbar, A. Performance analysis of a heat pump driven humidification-dehumidification desalination system. Desalination 2018, 445, 95–104. [Google Scholar]
  16. Zhang, Y.; Zhu, C.; Zhang, H.; Zheng, W.; You, S.; Zhen, Y. Experimental study of a humidification-dehumidification desalination system with heat pump unit. Desalination 2018, 442, 108–117. [Google Scholar] [CrossRef]
  17. Zhang, Y.; Zhang, H.; Zheng, W.; You, S.; Wang, Y. Optimal operating conditions of a hybrid humidification-dehumidification and heat pump desalination system with multi-objective particle swarm algorithm. Desalination 2019, 468, 114076. [Google Scholar] [CrossRef]
  18. Xu, H.; Zhao, Y.; Jia, T.; Dai, Y.J. Experimental investigation on a solar assisted heat pump desalination system with humidification-dehumidification. Desalination 2018, 437, 89–99. [Google Scholar] [CrossRef]
  19. Lawal, D.U.; Zubair, S.M.; Antar, M.A. Exergo-economic analysis of humidification-dehumidification (HDH) desalination systems driven by heat pump (HP). Desalination 2018, 443, 11–25. [Google Scholar] [CrossRef]
  20. He, W.F.; Yang, H.X.; Han, D. Thermodynamic analysis of a novel humidification dehumidification desalination system driven by heat pump. Energy Procedia 2019, 158, 6030–6037. [Google Scholar] [CrossRef]
  21. Kaunga, D.; Patel, R.; Mujtaba, I.M. Humidification-dehumidification desalination process: Performance evaluation and improvement through experimental and numerical methods. Therm. Sci. Eng. Prog. 2022, 27, 101159. [Google Scholar] [CrossRef]
  22. Liu, J.; Sun, Y.; Yv, S.; Wang, J.; Hu, K. Design and Experimental Study on a New Closed-Cycle Desalination System Based on Ambient Temperature. Processes 2020, 8, 1131. [Google Scholar] [CrossRef]
  23. Liu, Z.H.; Guan, H.Y.; Wang, G.S. Performance optimization study on an integrated solar desalination system with multi-stage evaporation/heat recovery processes. Energy 2014, 76, 1001–1010. [Google Scholar] [CrossRef]
  24. Behnam, P.; Shafii, M.B. Examination of a solar desalination system equipped with an air bubble column humidifier, evacuated duct collectors and thermosyphon heat pipes. Desalination 2016, 397, 30–37. [Google Scholar] [CrossRef]
  25. Rahimi-Ahar, Z.; Hatamipour, M.S.; Ghalavand, Y. Experimental investigation of a solar vacuum humidification-dehumidification (VHDH) desalination system. Desalination 2018, 437, 73–80. [Google Scholar] [CrossRef]
  26. Zubair, M.I.; Al-Sulaiman, F.A.; Antar, M.A.; Al-Dini, S.A.; Ibrahim, N.I. Performance and cost assessment of solar driven humidification dehumidification desalination system. Energy Convers. Manag. 2017, 132, 28–39. [Google Scholar] [CrossRef]
  27. Shafii, M.B.; Jafargholi, H.; Faegh, M. Experimental investigation of heat recovery in a humidification-dehumidification desalination system via a heat pump. Desalination 2018, 437, 81–88. [Google Scholar] [CrossRef]
Figure 1. Process flow chart of the closed cycle evaporation system [22].
Figure 1. Process flow chart of the closed cycle evaporation system [22].
Separations 10 00319 g001
Figure 2. Experimental set-up.
Figure 2. Experimental set-up.
Separations 10 00319 g002
Figure 3. Analysis unit for the humidifiers.
Figure 3. Analysis unit for the humidifiers.
Separations 10 00319 g003
Figure 4. Mathematical model flow chart.
Figure 4. Mathematical model flow chart.
Separations 10 00319 g004
Figure 5. Impact of seawater mass flow on freshwater productivity and the GOR.
Figure 5. Impact of seawater mass flow on freshwater productivity and the GOR.
Separations 10 00319 g005
Figure 6. Impact of seawater mass flow on the state of the air at the humidifier inlet.
Figure 6. Impact of seawater mass flow on the state of the air at the humidifier inlet.
Separations 10 00319 g006
Figure 7. Impact of seawater mass flow on the state of the air at the humidifier outlet.
Figure 7. Impact of seawater mass flow on the state of the air at the humidifier outlet.
Separations 10 00319 g007
Figure 8. Impact of the air mass flow on the system’s performance.
Figure 8. Impact of the air mass flow on the system’s performance.
Separations 10 00319 g008
Figure 9. Impact of the air mass flow on the state of the air at the inlet of the humidifier.
Figure 9. Impact of the air mass flow on the state of the air at the inlet of the humidifier.
Separations 10 00319 g009
Figure 10. Impact of the air mass flow on the state of the air at the outlet of the humidifier.
Figure 10. Impact of the air mass flow on the state of the air at the outlet of the humidifier.
Separations 10 00319 g010
Table 1. Comparison of the mathematical model simulation results and the experimental results of the closed cycle evaporation system.
Table 1. Comparison of the mathematical model simulation results and the experimental results of the closed cycle evaporation system.
ParametersUnit Condition 1Condition 2Condition 3Condition 4
Ta2°C 51.7051.3051.5051.90
ωa2g/kg 86.4879.3987.0085.51
mwkg/s 8.898.899.457.78
Makg/s 3.503.803.503.50
Ta1°CExp.43.4042.3044.0042.50
°CNum.42.9041.0544.1042.30
%Deviation1.152.960.230.47
ωa1g/kgExp.15.2514.3015.3014.70
g/kgNum.14.8014.2014.8014.20
%Deviation2.950.703.273.40
ia1kJ/kgExp.83.7580.1084.6081.46
kJ/kgNum.82.1078.6083.4079.19
%Deviation2.011.911.442.87
Ta4°CExp.20.6020.2020.7020.10
°CNum.20.2719.2720.6719.90
%Deviation1.604.610.161.00
mfg/sExp.249.30245.18250.81247.84
g/sNum.250.88247.72252.70249.59
%Deviation0.631.030.750.70
Table 2. Parameters under the condition of a fixed structure and variable working conditions.
Table 2. Parameters under the condition of a fixed structure and variable working conditions.
mw (kg/s)ma (kg/s)Ta5 (°C)φa5 (%)φa2 (%)
Reference values93.6206595
Interval7.5~10.52.8~4.2
Table 3. Comparison of different systems.
Table 3. Comparison of different systems.
ResearchersHeat SourceMaximum ProductivityGOR
Xu H. et al. [9]Solar-assisted heat pump20.54 kg/h2.42
Wu G. et al. [11]Solar182 kg/h2.65
Liu Z.H. et al. [23]Solar2.138 kg/(h∙m2)2.5
Behnam P. et al. [24]Solar6.275 kg/(day·m2)-
Rahimi-Ahar Z. et al. [25]Solar1.07 kg/(h·m2)3.43
Zubair M.I. et al. [26]Solar19,445 kg/year2.6
Lawal D.U. et al. [14]Heat pump287.8 kg/day4.07
He W.F. et al. [19]Heat pump82.12 kg/h5.14
Shafii M.B. et al. [27]Heat pump2.79 kg/h2.08
Current studyHeat pump852~867 kg/h4.16~4.24
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, J.; Sun, Y.; Zhang, Y.; Wang, J. Theoretical Study of a Closed-Cycle Evaporation System for Seawater Desalination. Separations 2023, 10, 319. https://doi.org/10.3390/separations10050319

AMA Style

Liu J, Sun Y, Zhang Y, Wang J. Theoretical Study of a Closed-Cycle Evaporation System for Seawater Desalination. Separations. 2023; 10(5):319. https://doi.org/10.3390/separations10050319

Chicago/Turabian Style

Liu, Jun, Yong Sun, Yizhu Zhang, and Jiaquan Wang. 2023. "Theoretical Study of a Closed-Cycle Evaporation System for Seawater Desalination" Separations 10, no. 5: 319. https://doi.org/10.3390/separations10050319

APA Style

Liu, J., Sun, Y., Zhang, Y., & Wang, J. (2023). Theoretical Study of a Closed-Cycle Evaporation System for Seawater Desalination. Separations, 10(5), 319. https://doi.org/10.3390/separations10050319

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop