1. Introduction
Water is the world’s most precious resource yet is becoming increasingly scarce due to the competing demands of the world’s burgeoning population. Finding solutions to the water crisis that the world is currently facing is a significant challenge, as is the need for greatly improved global equitable and sustainable water resources management. To date, water desalination is the most viable solution, but it has a large ecological footprint, as it is energy-intensive and environmentally taxing. Heightened water requirements for domestic water supply, sanitation, and irrigation purposes, are the leading causes of stress on the global availability of water [
1]. Such water-fueled stress has serious social, economic, and developmental implications. The United Nations (UN) reported that 2.2 billion people lack access to safely managed drinking water, more than 3.6 billion people do not have adequate sanitation, and more than 2.3 billion people do not have basic hygiene facilities [
2]. Adopted by the UN Member States in 2015 at the UN Summit, the Sustainable Development Goals (SDGs) of the 2030 Agenda for Sustainable Development seek to reduce global poverty and achieve enhanced worldwide prosperity. SDG 6 relates specifically to providing the sustainable management of clean water and sanitation for all. Such is the importance of water; it is the common theme in almost all the SDGs. It is essential to food security (SDG 2); health and sanitation (SDGs 3 and 6); energy, employment, and industries (SDGs 7, 8, and 9); climate, and the protection of ecosystems (SDGs 13, 14, and 15) [
3]. The Middle East and North Africa (MENA) region is the planet’s most arid terrain mainly due to harsh weather conditions, increasing dependence on imported water bottles, and seawater desalination [
4]. Annual freshwater withdrawals refer to total water withdrawals from renewable aquifers used in all sectors such as agriculture, industry, and municipal or domestic. It is noteworthy that most MENA countries fall in the extremely high-stress region where the withdrawal is more than 100% of the renewable sources, and seawater desalination is involved in covering the shortage [
5], see
Figure 1. With an ever-increasing population and rising water demands, it is predicted to confront severe water shortages and even warmer climates [
3]. Of particular concern in this region are the water-intensive methods used in agriculture and the energy-intensive method in desalination [
6].
The world will continue to face significant challenges where the water supply is limited, and the quality is variable unless a sustainable water source can be found and managed efficiently. The sole source of accessible water on the planet is moisture in the air [
7], and it is thought that approximately 13 sextillions (13 × 10
21) L of water exist in the atmosphere at any given time [
8]. To date, harvesting water from the atmosphere occurs using three processes and is a promising solution to the world’s water crisis [
9]. One such method involves desiccant materials that act as absorbents of atmospheric moisture and from which water vapor is then released. This research presents the results of experiments conducted to assess and optimize the behavior of Calcium Chloride desiccant as a means of efficient Atmospheric Water Harvesting (AWH).
The following section,
Section 2, contains a literature review on atmospheric water generation using hygroscopic salts;
Section 3 outlines the materials and methodology used in the different experiments;
Section 4 presents the results and discussion and highlights the main findings of this study. In
Section 5, we present an empirical mathematical formula. Conservation governing equations have been developed and solved numerically in
Section 6. Several scenarios were simulated and discussed in
Section 7. Finally, conclusions are drawn in
Section 8.
2. Literature Review
In recent years, AWH has gained popularity due to the rising water shortage problem, stressed water systems, issues in public water infrastructure, and natural disasters [
4]. As mentioned previously, in recognizing the increasing challenge of water shortage, in 2015, the UN established SDG 6 (among 17 SDGs) to “Ensure availability and sustainable management of water and sanitation for all” [
2]. In 2018, the UN launched the Water Action Decade 2018–2028 to make significant progress toward sustainable development [
10]. Because the quantity of fresh water available on Earth is just 2.5% compared to 96.54% of salty water in oceans, seas, and bays [
7], significant efforts have been devoted to harvesting water from clouds, fog, or atmospheric water vapor, particularly in dry regions [
8,
9,
10,
11,
12]. Two main methods of atmospheric water generators have been used for atmospheric harvesting water: passive and active methods [
13,
14]. With passive techniques, water may accumulate naturally without needing an external power source, relying only on natural temperature variations. The most employed passive method is a fog catcher, which generates a reasonable amount of water [
8,
15,
16]. However, they only work in a limited number of geographical areas and are restricted to specific climate conditions such as high humidity. Besides, fog collection technology only contributes a little to alleviating water shortages.
In contrast, active methods used for AWH are not restricted to a specific location and are less affected by site or climate conditions. Still, they require an external source of energy, such as solar energy. Atmospheric water generators can use solar energy (refrigeration-based) or solar thermal (sorption-based) energy. The amount of water generated depends on relative humidity and ambient temperature.
Bar [
17] patented a three-stage process technology that extracts water from the air. A solid desiccant absorbs humidity, and a heat pump heats it to desorb water vapor. Then the vapor turns to water by passive condensation. Such a system is environmentally friendly due to moderate low-cost heating (solar or waste heat). With a low energy consumption of 100–150 kcal/L, the system generates water at a rate of 1000 L/day. Extracting water by cooling air to the dew point produces low water output in passive techniques and consumes high energy with active techniques. Alternative methods that capture water molecules regardless of the dew point are more energy-efficient, especially at low relative humidity [
11]. Membrane-assisted humidity harvesting is an active method that depends on mechanical energy (vacuum pump) to push water molecules through a selective membrane and then condense it to produce water. With the help of a low-pressure sweep stream, the membrane-assisted device harvests an average amount of 3.36 L/day [
18].
Because of the abundance of solar energy, solar-driven hygroscopic water harvesting is increasingly scrutinized, as researchers are working on improving structural designs, sorbents’ capacity, and thermal management [
13,
14,
19]. Many challenges are still observed with these types of Atmospheric Water Generators (AWG), such as increasing efficiency, scalability, and cost reduction. This implies ensuring that research should focus on the microscopic mechanism of the process, develop new hygroscopic agents, modify solar absorbers, and design new condensation systems. Talaat et al. [
20] tested a double-cycle absorbent-based water harvester where the cycles alternate between day and night. At night, the apparatus exposes the Calcium Chloride (CaCl
2) solution to natural humid air. In the daytime, the system is covered with a double-faced transparent surface, forcing the water vapor to desorb due to high temperature, which then condenses on the cover sides and is collected inside a bottle. This system is affected by weather and ambient conditions, generating a maximum amount of water of 0.63 L/day. Sleiti et al. [
21] recommended multiple layers of stacked high-quality sorbent to improve the water output of the AWG unit. A layer of 25 mm thickness of silica gel generated around 0.8 L/day, and the unit consisted of a water sorbent bed, condenser, radiant heat flux, and reflector. Due to different experimental conditions such as ambient temperature, relative humidity, silica gel thickness, surface area to volume ratio, and radiant heat flux, the amount of water recovered varies. Increasing the radiant flux will reach peak efficiency and then drop; the relative humidity is proportional to water capture.
Solar-driven hygroscopic AWGs’ water production depends on external environmental and internal device factors. Understanding the interaction between external and internal factors such as temperature, pressure, wind speed, relative humidity, sorbent characteristics, types of solar collectors, and condenser design is the key to manufacturing an effective scaled-up product in the future. These parameters may only affect the production rate or maximum sorbent capacity or affect both rate and capacity. In the case of desiccants’ absorption, maximum water absorption is related to their inherent characteristics and the surrounding environment. Changing the chemical properties of the desiccant or making a physical modification to the system, such as using a porous medium as a host, will increase the water uptake.
Several studies demonstrate that some changes may increase the water uptake but decrease the absorption rate, as in the work of Zhang et al. [
22]. They used porous silica gel as a host for CaCl
2, but the equilibrium process took 35 h, while using carbon nanotubes only took 4 h. Moreover, a composite of 85% CaCl
2 and 15% Lithium Bromide (LiBr) was found to absorb three times its weight, while CaCl
2 alone only absorbed twice its own weight [
23]. Experiments of Elashmawy et al. [
24] clearly demonstrated the effect of increasing the temperature of the saturated sorbent by adding a solar parabolic concentrator to the tubular solar-still with a CaCl
2 trough inside and increasing its temperature to 98 °C, thereby tripling water production to 2.4 L/m
2/day. One of the advantages of using desiccants instead of cooling to dew point is that the cooling method loses its reliability below 30% relative humidity. On the other hand, most desiccant-based AWGs succeeded in generating a reasonable amount of water in low relative humidity conditions.
Sorbent characteristics can differ between absorption and desorption; one such example is the absorption rate of Lithium Chloride (LiCl), which is higher than CaCl
2 but lower in the desorption phase [
25]. In some ways, the carrying capacity factors (temperature, relative humidity) may also affect the process rate. However, on the other hand, factors such as surface area, wind velocity, and liquid desiccant flow are only related to how fast/slow the sorbent will reach equilibrium. Liquid desiccants have excellent flexibility with respect to their surface area, and can adopt many host surfaces, such as porous silica gel and carbon nanotubes, as mentioned earlier [
22]. Elashmawy et al. [
26] used a black cotton cloth as a host for their experiments and collected 1.06 L/m
2/day using five shelves, to increase surface area. This can be compared to a 1 kg flat layer of saw wood with a 10–40 μm distance hosting a 60% concentration of CaCl
2, which only collected around 0.18 L/m
2/day [
27]. Furthermore, the sand hosting 37% CaCl
2 concentration solution was able to generate 0.11 L/m
2/day with the help of a solar Scheffler reflector [
28]. Since the sorbent creates a low vapor pressure area around its surface, water molecules in the air will move to fill it; therefore, circulating more water molecules over the surface of the sorbent will speed up the saturation. In the city of Hail (northern Saudi Arabia), a forced airflow was injected in solar-still tubes containing black cotton cloth impregnated with CaCl
2 solution via small fans at different speeds. The highest absorption amount was from the sample subjected to the fastest fan (4 m/s) and collected 0.47 L/m
2/day [
24]. It was thought that the absorption rate of CaCl
2 solution fluctuated due to the wind speed in Dhahran, the eastern region of Saudi Arabia. This was verified by an experiment that demonstrated that the high flow of the desiccant solution achieved the highest absorption rate of 2.11 L/m
2 at night. The opposite happened with regard to desorption during the daytime, where the lower flow rate caused a greater amount of collected water [
29].
Before utilizing desiccants to generate freshwater from the air, researchers focused on the desiccants’ industrial and air conditioning properties [
16,
30,
31]. Properties such as solubility, solution vapor pressure, density, surface tension, viscosity, and thermal conductivity are essential to understanding the behavior of the desiccant solution at various concentrations and in certain environments. Thus, despite recent studies on AWG, utilizing desiccants in the AWG field requires conducting further detailed experimental investigations. Moreover, from the aspect of mathematical modeling, few studies have been conducted to study atmospheric water harvesting in recent years. Milani et al. [
32] developed a model for a complete dehumidification system using TRNSYS software. Sibie et al.’s [
33] modeling and simulation focused on Copper Chloride (CuCl
2), Copper Sulfate (CuSO
4), and Magnesium Sulfate (MgSO
4) anhydrous salts and introduced an analytical solution for the mathematical model. Furthermore, Alkinani et al. [
34] developed a CuCl
2 water vapor absorption mathematical model. The experiments in this paper aim to discover CaCl
2 desiccant behavior under AWH-related conditions and provide a mathematical model describing the AWH behavior, thereby enhancing the absorption and desorption processes that can be applied to various liquid desiccants.
5. Proposed Formulas
The exponential decay equation is commonly used in modeling chemical reactions, electrostatics, heat transfer, and optics. By reviewing the experimental data, the desiccant appears to follow an exponential decay function in both forms, decay for desorption and growth for absorption. CaCl
2 and similar desiccants may also follow the curve of the exponential decay/growth function shown in
Figure 10. Internal factors can be thought of as the type of desiccant, solubility, surface area, and depth. External factors are assumed to be temperature, relative humidity, and wind velocity (airflow).
On considering relevant internal and external factors affecting both absorption and desorption, the formula that describes this study’s experimental data on CaCl
2 can be expressed as the following:
where
where
equals the mass ratio of water-to-CaCl
2 after elapsed time
t and
is its initial value,
b is a variable related to the absorption/desorption rate due to temperature. The ratio of ambient temperature and solution temperature is
τ.
T is the temperature in Celsius and
RH is the relative humidity. The computation of the saturation water vapor pressure,
, at a given temperature, has been adopted from the work of Huang [
36]. His formula is less complicated than previous ones, with a lower error.
is the water vapor partial pressure in kPa.
Finally, is the maximum amount of water ratio the salt can absorb in a certain ambient condition, in other words, the vapor–liquid dynamic equilibrium points at a particular condition. The equation is derived based on the observation of the linear relationship between and water vapor pressure.
The equation’s constants
are optimized using the Generalized Reduced Gradient (GRG) algorithm [
37]. The convergence is set to 0.0001, the start set to multi, and the derivative is set to central. The target was to minimize the average absolute error between measured saturation points and the
equation output.
The mean absolute error (MAE) was found to be 6.9%. As the temperature of the solution increases, the values deviate from the actual readings due to crystallization (see
Figure 11). The constants
generated and found by the optimization are presented in
Table 2.
It appears that the proposed formula describes the absorption process better than the desorption process. The formula predicted the water absorbed during the day (see
Figure 12) more accurately than the desorbed amount when the solution is heated up to 45 °C (see
Figure 13). This implies that for calculating
, the initial water content is required.
The mean absolute error for the formula in the absorption process was 1.9% after fixing b to equal 0.05 (see
Figure 14a). The same formula gave a MAE of 3.56% in the desorption process after b was set at 0.294 (see
Figure 14b).
6. Mathematical Modeling
Since the sample of the desiccant solution is contained in a small cylindrical beaker, the only important axis is the depth. However, the distribution of the solution concentration concerning radial and azimuthal axes is uniform, assuming isothermal walls. Moreover, the model represents all the expected physics in these phenomena. It contains a phase-change term, coupled with the empirical equation, in addition to the energy equation, which is necessary for the desorption case.
The CaCl
2 moisture absorption depends on water vapor sorption at the phase boundary between the air and the solution and CaCl
2 diffusion in the water. However, it is a slow process compared to the water vapor diffusion in the air. By neglecting the diffusion of CaCl
2 in water and considering the case of a thin layer solution, solving the calculations will be easier. A one-dimensional mathematical model in the z-direction is therefore proposed. A system of two transport (mass conservation) equations can solve the sorption phenomena; however, to solve the desorption process, a heat equation is needed:
where
is the water vapor density [kg/m
3],
is the water vapour diffusion coefficient in the solution [m
2/s],
is the density of water-CaCl
2 mixture [kg/m
3] and
and
are the solution-specific heat and thermal conductivity, respectively. The transfer term
represents the mass transfer between vapor and liquid equation equations. We suggest the following definition of
,
In order to close the system of equations, it has been linked to the proposed empirical equations,
where
are defined in the previous section (
Section 5),
and
are defined in
Table 2,
is the ideal gas constant 8.3136 [
] and
is the molar mass of water 0.018 [
].
The initial conditions for Equations (7)–(9) are considered uniform and expressed in Equations (14)–(16), but the boundary conditions are defined only for the water vapor and solution temperature for Equations (7) and (9) in the Equations (17)–(20).
Boundary conditions:
Upper Boundary | |
|
Lower Boundary | |
|
where
is the water vapor density related to the relative humidity at the surface of the solution and
is the mass ratio of water-to-CaCl
2 related to the initial concentration of the solution.
is the solution temperature and usually equals ambient temperature,
is the temperature of the controllable heat pad (set to high temperatures to mimic solar heating for desorption),
is the water vapor density related to the ambient relative humidity of the surrounding airflow [kg/m
3] and
is the depth of the solution [m].
The above highly nonlinear parabolic partial differential equation is solved numerically using an efficient algorithm. The governing Equations (7)–(13) are solved along with their initial and boundary conditions (14–20). The Galerkin method handles spatial discretization, while an adaptive time step is used with time integration. The model is solved by using the “pdepe” function in MATLAB software that can solve a system of parabolic and elliptic partial differential equations. The pdepe is capable of solving the differential-algebraic equations that arise when the PDE contains elliptic equations by using the ode15 s solver, and for handling Jacobians with a specified sparsity pattern.
Figure 15 shows the modeling results for both absorption/desorption processes.
The actual readings from the thin layer CaCl
2 solution are compared with the mathematical model in
Figure 16. In the absorption process, the initial conditions for
and
were set to 17%, 0.011 kg/m
3 (around 48% relative humidity), and 26.6 °C, respectively. The boundary condition for
only is set as 0.011 kg/m
3, since we assumed the temperature of the solution remains constant, and the MAE was determined to be 3.13%. In the desorption process, the initial conditions for
were set to 89%, 0.011 kg/m
3, and 25.6 °C, respectively. The boundary conditions for
and
set as 0.011 kg/m
3 and 45 °C, respectively. Since the solution is at a high temperature, and the model only considers the average temperature and relative humidity of the air around the solution, regardless of the forced convection applied to air at the surface of the solution, this results in a less accurate
θmax equation than the absorption model, so the model’s MAE was calculated to be 7.32%.
7. Simulation Scenarios
Various simulation scenarios of water absorption and desorption were performed using the proposed mathematical model to approximate the amount of water that could potentially be generated in the city of Jeddah, KSA.
Table 3 shows the average values for day/night temperatures, relative humidity, and nighttime hours in Jeddah, used as data inputs for the model. The solution temperature is assumed to match the temperature during the daytime and the condensation temperature is the same as the ambient temperature at night.
The initial amount of water in the solution changed in each simulation until we ended with the same amount, to complete the absorption/desorption cycle.
Figure 17 shows the simulation results of one day (one cycle) in each season. Assuming all desorbed water is condensed to water, the maximum amount that can be collected is in Fall, which equates to the highest relative humidity average in
Table 3, and around 37% of the CaCl
2 mass. This means that an amount of 370 mL of water can be collected from 1 kg of CaCl
2/m
2 day at a depth of between 1.6–2.0 mm.
Changing the amount of desiccant will need more depth per m
2 and it will reduce the absorbed water ratio, as shown in
Figure 18. Since the initial mass ratio is fixed, each increase in the depth means doubling the amount of CaCl
2 in the solution, yet the absorbed amount is less than double.
The simulation shows a linear relationship between collected water and depth of the solution (see
Figure 19). A total of 1800 g of CaCl
2 can absorb 2160 mL of water from an area of 3 m
2 at 1 mm depth, but the same 1800 g of CaCl
2 can absorb only 630 mL if the area is 1 m
2 at 3 mm depth. This shows the importance of finding a solution to increase the absorption and desorption in deep containers. Adopting a thin layer approach for designing future AWG devices will consume vast areas and may complicate the maintenance process.