1. Introduction
Air conditioning (AC) is one of the primary human activities that contributes to heavy power consumption, particularly during the summer months. According to data provided by International Energy Agency (IEA, Paris, France), the global power consumed by air-conditioners and electric fans in building sector is around 20% of the total power consumption [
1]. That, in turn, would exacerbate the environmental pollution since the key sources of electricity production worldwide are still fossil fuels (i.e., coal, oil, and natural gas). The traditional vapor-based compression cooling systems use ammonia [
2], carbon dioxide [
3], hydrofluorocarbons (HCFCs) [
4], and recently hydrates [
5,
6] as refrigerants. These refrigerants, however, pose certain environmental rinks. They can contribute to the global warming both directly via the release of refrigerant into the atmosphere and indirectly by the release of carbon oxides during the combustion of fossil fuels in the power plants that used to generate the electricity needed to keep the cooling equipment running. Solar adsorption cooling systems are found to be very efficient alternatives for dealing with these challenges. They neither use environment-damaging refrigerants, nor require fossil-fuel based electricity. These cooling systems can be driven by low grade energy like solar energy which is more abundant and more evenly distributed in nature than other energy sources. Moreover, the good solar radiation conditions in summer, which closely corresponds to the periods when the peak demand for cooling is expected, enables the adsorption cooling to be a reliable technology for sustainable energy utilization. This technology can also aid to climate-change mitigation with making progress toward the United Nations sustainable development goals of reaching net-zero carbon emissions by the second half of the twenty-first century [
7,
8]. The main objective for utilizing the adsorption phenomenon in cooling systems is to improve the evaporation rate of the coolant fluid in use. This improves the system’s cooling efficiency as the potential for removing more heat from the cooling space increases. Therefore, the cooling efficiency of adsorptive systems depends on the evaporation rate, which in turn, largely depends on the thermophysical properties of the adsorbent-adsorbate working pairs and the effective design of the porous adsorbent bed to provide the necessary heat and mass transfer augmentation. These influencing factors have been a subject of interest for numerous studies, and several cooling, refrigeration, and desalination applications have been explored using different working pairs with different bed designs [
9,
10,
11]. The adsorption cycles with silica gel/water or zeolite/water as working pairs are found to be more functional than others for the use with solar panels or flat plate collectors as they can be powered by quite low, near ambient temperatures [
12,
13,
14]. The two-chamber design for adsorption/desorption processes is usually applied to enable the adsorption chillers to operate alternatively so that intermittence-free cooling cycle is performed [
15].
Many theoretical and technological advancements were conducted on water adsorbents like silica gel and zeolite to improve their utilization in the adsorption cooling systems. Miyazaki and Akisawa [
16] examined the operating parameters of silica gel–water adsorption chiller based on single-bed design and observe that although increasing the size of adsorbent bed can decrease the coefficient of performance (COP) due to the shorter cycle time, it can also increase the specific cooling power (SCP). Wang et al. [
17] measures how the pollution by the impurity ions or solid particles would increase the adsorption deterioration of silica gel, and found that soaking in acidic solution along with washing via distilled water would be a very good solution to restore the silica gel adsorption capacity. Effect of grain size and grain layers on the adsorption performance of silica gel was numerically examined by Chakraborty et al. [
18]. It is found that COP and SCP can be both improved if smaller grain size is applied. Freni et al. [
19] tested a novel composite of silica modified by calcium nitrate with water sorbent (SWS-8L) and the resulting COP based on experiments was 0.18–0.31 for cycle time of 10 min. Wang et al. [
20] developed a lumped-parameter model for multi-bed silica adsorption chiller and showed that adopting a multi-bed design can enable the chiller to utilize heat source at a temperature as low as 65 °C more effectively than its mature two-bed technology. Manila et al. [
21] proved through numerical simulation that the two-bed design with smaller cylinder’s depth and optimum silica gel grain diameter of 0.8 mm is preferable for faster vapor uptake and lower bed pressure drop.
The good match between the high sunshine and the peak air-conditioning demand during summer represents the most feasible utilization of solar energy. Therefore, the feasibility of integrating solar collectors with adsorption chillers have been extensively investigated in literature. To enable more economic and technical feasibility of dual-stage absorption chillers in the building sector, Chemisana et al. [
22] tested integrating of highly effective solar collectors (evacuated tube collectors) with adsorptive chillers to achieve driving temperatures as high as 150 °C and found that this arrangement can provide a substantial reduction in the collector area requirements. Habib et al. [
23] evaluated by numerical simulation the single-stage and dual-stage operation modes of a four-bed silica gel-water adsorption chiller driven by evacuated tube collectors in Durgapur, India and found that the operation mode needs to be switched to dual-stage mode when the heat source temperature is below (60 °C). Drosou et al. [
24] discussed the feasibility of parabolic trough collectors for use with solar adsorptive cooling system to cover the cooling needs of typical office building in Greece, and proved that the use of parabolic trough can bring higher COP within less space requirements compared to traditional flat plate collectors. Wang et al. [
25] studied the effect of the mass transfer on the performance of activated carbon-methanol adsorption chiller and demonstrated an improvement of 35.9% compared with the average COP of the natural mass transfer adsorption refrigeration system when the input radiation energy was not less than 14.7 MJ during a refrigeration cycle.
The performance of adsorption cooling systems powered by solar energy has been the subject of numerous research studies [
26,
27,
28,
29,
30,
31,
32]. However, there were no reported studies adopting their performance under unequal adsorption/desorption phases times with long switching times. For an equal state, the time of adsorption is equal to the time of desorption. If there is a significant difference between the two phases (adsorption/desorption) times it is called unequal state. Under this circumstance of unequal state, the adsorption capacity rate is not equal to the desorption capacity rate. The difference of both states can lead to the distinct difference for the value of the system performance. However, in realistic cooling systems driven by solar energy, several factors have a strong effect on the adsorption/desorption phases times, especially for the systems that are operated under short cycle. For example, varying heating/cooling water temperatures or their mass flow rates, as a result of equipment obsolescence and weather changes, are commonly arising in the real cooling process. Based on these characteristics, this paper presents a dynamic model for an adsorption cooling system driven by solar evacuated tubes collectors under unequal states conditions. The existence of varying adsorption/desorption times with long switching times for a single bed cooling system was investigated experimentally in our earlier work [
26]. Moreover, the minimum adsorption time for the present system is fixed, which indicates the shortest adsorption time at which the system still in service of producing cooling effect. Thus, the influence of the unequal state needs to be considered and cannot be neglected for obtaining the best model results. During these analyses, the roles of adsorption velocity, adsorption capacity, overall heat transfer coefficient, and the performance of the cooling system are also elucidated.
2. System Description
A schematic diagram of the main arrangements of the adsorption cooling system under investigation is illustrated in
Figure 1. This system has been built and implemented by Lattieff et al. [
33] for solar-powered adsorption air-conditioning application under the climatic conditions of Iraq. In the present study, the same system has been retested under the unequal adsorption/desorption operating time conditions. The principal innovation of this concept is the capability to run the cooling system under different operation conditions with variable adsorption/desorption intervals. The system consists of three components, namely: the supplier of hot water which uses a set of solar evacuated-tube collectors as heat source, the supplier of cool water which is obtained from the municipality network, and the adsorption chiller equipment. The chiller which is considered the main component consists of four smaller parts: adsorber, evaporator, condenser, and set of valves that are operating on an intermittent basis according to whether the active mode is adsorption or desorption. The adsorber is fabricated as a shell-and-tube heat exchanger having diameters 17-mm ID and 19-mm OD with total heat-transfer area for the adsorber of 0.35 m
2. The tube is made of copper embedded in silica-gel adsorbent type A as the shell of adsorber heat exchanger. The adsorbent was embedded in an aluminum mesh for enhancing the rates of heat and mass transfer inside the adsorber. The condenser is fabricated as 1-m coil tube inserted in a metal outer shell currying 2-L volumetric capacity of stagnant water. The evaporator is a small cylindrical vessel containing 0.25 L of distilled water as coolant located at the bottom of the chiller equipment.
A solar water heater consisting of two evacuated tube collectors, which were made up of parallel rows of transparent glass tubes and manufactured by gtc-solar co.ltd./China, was used to provide the required heating load. Each collector consists of a set of 20 evacuated tubes of 1.25 m in length and 0.025 m in diameter. A cylindrical storage tank with a capacity of 0.12 m3 was used for collecting the hot water from the solar heater. A hydraulic pump was used to transport by forced circulation the heat-transfer fluid to the collectors. A rotameter (Model DH4500/LZS) with accuracy of ±4% and measuring range of 20–170 L/min was used for water flow rate measurements. The scenario begins when the silica-gel bed is linked firstly to the evaporator for achieving the first half of the cycle, i.e., the adsorption mode. Once the evaporation comes to end up at the evaporator side, the silica gel in the bed side gets saturated with water vapor that is being drawn from the evaporator. Since adsorption is classified as an exothermic process, the existence of an external cooling is necessary for cooling the silica gel down so that a continuous uptake of vapor is preserved. In the second half of the cycle, the bed is linked to the condenser to condensate the water vapor by releasing its heat content to the coolant flowing in the condenser. In between these two processes, there is a process when water vapor circulation is stopped by closing all coolant valves. This process is known as the preheating/precooling process. The proposed operating time for each process in the cycle is: 10-min pre-heating, 20-min desorption, 20-min pre-cooling, and 10-min adsorption.
3. Mathematical Modelling
In order to evaluate the efficiency of the adsorption cooling system under the conditions of unequal adsorption/desorption times, a mathematical model of the adsorption/desorption (A/D) cycle was developed and validated. This model employs an interconnected cycle of the evaporator, condenser, and adsorber components to analyze the efficiency of the cooling system. The hot water flowrate is assumed to be constant during the desorption half cycle, and the cooling water flowrate undergoes a fluctuation during the next adsorption half cycle. In this model, the chiller system is assumed to be operating with equal mode (10 min pre-heating + 20 min desorption + 10 min pre-cooling + 20 min adsorption) to demonstrate the effectiveness of unexpected physical changes on the system performance. In evaluating the model being developed, the mass and energy balances of silica gel-water as the absorbent-adsorbed pair is used under the relevant temperature and pressure operating conditions. The lumped-parameter mathematical model of the adsorption chiller is employed to predict the temperature and coolant concentration variation in the bed. Then, the influence of shifting the cycle time from unequal mode to equal one on the performance of the chiller is discussed from the viewpoint of the COP and SCP.
The energy balances that govern the temperature distribution in the adsorber bed during the adsorption and desorption modes are expressed as
where the left-hand side represents the change rate of sensible heat in the adsorber bed components, including the adsorbent material, adsorbate material and metal tubes. The right-hand side includes the rates of heat input and heat release during adsorption. The outlet temperature of the hot water
from adsorption bed and cold water from desorper
are estimated using the log mean temperature difference (LMTD) method as
where Δ
Ha,
,
, and
represents the heat of adsorption, silica-gel heat capacity, heat capacity of adsorbed water, heat capacity of metallic part in the A/D bed, respectively.
In condenser, the energy balance is expressed as
The outlet temperature of cooling water from the condenser is described by
In the evaporator, the energy balance is expressed as
The outlet temperature of chilled water from the evaporator is described by
The cooling effect in solar-powered A/D cycles is primarily attributed to how high is the uptake rate of the water vapor particles by the adsorbent materials. Therefore, the water adsorption/desorption rate (x) is considered the major important parameter for evaluating the efficiency of adsorbent (silica-gel). The water adsorption/desorption rate is estimated using Dubinin-Astakhov model as
where
Pb is the saturation pressure of water vapor in the adsorption bed. It is estimated using the correlation:
Km is referred to as the overall mass transfer coefficient being defined as
Here
Do,
Rp,
Ea, and R denote the surface diffusion factor of the adsorbate in absorbent particles, the radius of absorbent particles, the surface diffusion activation energy, and the universal gas constant. For silica gel-water, these parameters take the following numerical values: 2.54 × 10
−4 m
2/s, 3.5 × 10
−4 m, 4.20 × 10
4 J/mol, and 8.314 J/mol.K, respectively. The linear driving force (LDF) model is employed for the estimation of the water uptake velocity by the adsorbent particles as
The heat brought from the heat source during desorption mode (
Qin) is estimated as:
The heat released by working fluid in the evaporator section (
Qevap) is estimated as:
Finally, the coefficient of performance (COP) and the specific cooling power (SCP) for one cycle can be estimated as:
The parameters used in the simulation model with their numerical values are listed in
Table 1.
4. Results and Discussion
The performance of the intermittent cooling system with different adsorption/desorption operating times was investigated in terms of adsorption velocity, adsorption capacity, overall heat transfer coefficient, COP, and SCP. After validation of the model developed with the experimental findings, a parametric analysis of the system performance under both the equal and unequal phases times was conducted and reported. The obtained results were calculated by considering the inlet mass flow rate of 50 L/min with total cycle time of 60 min being divided into 30 min for adsorption and 30 min for desorption. The flow rate of 50 L/min was selected as this rate yields the largest bed temperature lift which causes a highest cooling capacity. Increasing the cycle time more than 30 min has an insignificant effect on the cooling capacity. Therefore, adsorption/desorption time of 30 min was chosen to run this system. The performance of this proposed system in terms of (SCP and COP) was evaluated by using dynamic model to calculate the adsorption time variation that simulates the actual operation of the real system (single bed).
Figure 2 depicts the variation of the bed temperature with the cycle time for both cases, equal and unequal adsorption/desorption phases times. A coincide agreement can be observed in the first 30 min during desorption phase, while during the adsorption phase a slight deviation can be noticed at the last 10 min when the unequal adsorption phase is starting. This is because the unequal state has a long pre-cooling time (20 min) and the short adsorption time (10 min) that makes the adsorber pressure is much lower than the evaporator pressure. Moreover, starting the unequal adsorption process resulted in a larger amount of adsorbed water vapor and this is quite crucial to affect the bed temperature. Identical trend for both states confirm the validity of the model.
Figure 3 shows that the desorption velocity, 0.00012 kg water/kg silica.s, is highest at the start of the desorption process and becomes very slow after 20 min of operation. As per this figure, the adsorption velocity for both states rise sharply at the starting of the adsorption phase, but after that they move to decrease gradually with the progress of time. However, the unequal state has a maximum adsorption velocity value of about 0.0002 kg water/kg silica.s, which is noticeably higher than that of the equal state, 0.00008 kg water/kg siliaca.s. The reason of this difference is due to the variance in pressure between the evaporator and the adsorber as previously mentioned in
Figure 2 and this difference is reduced to a very small value. It can be concluded that the faster the adsorption velocity is, the lower the cycle time becomes. The figure also indicates that the system can be driven more stably and efficiently by the equal state. This information is very useful when managing a cooling system being tailored to operate at the optimum adsorption velocity. This result is comparable to that founded by Yang PZ [
34] who reported that the adsorption/desorption will not completely accomplish and the cooling potential of the adsorber cannot exert deeply if the adsorption time is too short.
A complete cycle of adsorption/desorption for the silica gel-water adsorber is shown in
Figure 4. As seen in this figure, desorption process is very close in concentration change to adsorption process, and many studies assume that both are the same in adsorption refrigeration process. During the desorption phase, the minimum value reached is 0.06 kg/kg after the first half-cycle of 30 min, and during the second half-cycle of 60 min, the maximum is 0.12 and 0.124 kg/kg for unequal and equal cycles, respectively. It is clear that changes in adsorption rate with time are not simultaneous. This means that the adsorption rate cannot be entirely reached when the unequal adsorption time reaches its end value. Using these results, it was determined that 4% of silica gel was not effectively used in unequal cycle compared with equal cycle. It could be observed that the equal adsorption rate deviates from the unequal one and this deviation does not seem very high at the adsorption completion point. Without doubt, the unequal state is influenced by the real process of cooling system and should not be neglected on designing cycle adsorption chiller, especially when the cycle is short because in real short cycle, the adsorbed water could not be sufficient due to the non-equilibrium adsorption process [
35].
Figure 5 aims to determine the amount of exchanged water vapor involved in one adsorption/desorption cycle. It can be noticed that the bed is able to exchange a large amount of 0.064 kg/kg.s under equal condition and 0.06 kg/kg.s under unequal condition, so it is typical for air conditioning applications. The main reason for such differences between the equal and the unequal simulated data presented in
Figure 5, would be the pre-assumption that the governing parameters like density, specific heat, and viscosity of the working fluid are pressure and temperature independent. Another reason would be the adoption of a simplified mathematical model for the heat-and-mass transfer simulation. It is worth to mention that the results of this work reveal low values compared with the other previous works [
19,
36,
37].
It can be realized that by far the heat-and-mass transfer problems in adsorptive cooling systems are not fully solvable. The major obstacle that prevents the employing of the adsorption cooling for widespread use is the low heat-and-mass transfer in the adsorption cooling system. Therefore, enhancing the heat-and-mass transfer in the adsorber to obtain faster adsorption/desorption rates is the key for boosting a superior adsorption cooling efficiency. In
Figure 6, we have reported the overall heat transfer model results; the figure depicts the equal and unequal times evaluation of predicated overall heat transfer coefficient. The comparison between these values shows a significant difference. This difference may be due to the amount of the adsorbed water vapor is also not stable and is changing with time. Yet, the lack of an efficient thermal interaction between the adsorbent and the metal part of the adsorber results in an undesirable steep temperature gradient at the contact surface. As a matter of fact, stimulating the adsorbent to be in a good direct contact with the metal is not easy in practice, particularly when granular adsorbent is utilized. Therefore, several alternative designs of adsorptive cooling systems were proposed to mainly improve the internal heat and mass transfer [
38,
39,
40,
41,
42]. In the majority of these systems, a constant
Ub value is assumed for implementing their models, but they couldn’t fully achieve the complete convergence. This work proves that the use of variable heat transfer coefficient gives more accurate and realistic results with closer convergence between the model results and the experimental data. It is clear that the heat transfer coefficient of the equal state is constant because of the semi-stable adsorption velocity as shown in
Figure 3. Unequal heat transfer coefficient results show a similar trend to that of adsorption velocity trend depicted in
Figure 3. It indicates the effect of the adsorption process on the heat transport rate between the adsorbent and its heat exchanger.
Variation of the desorbed water in the condenser with cycle time is predicted in
Figure 7. It is clear that the condensate water rises sharply due the flow of the vapor from the desorber, then after it decreases due to the decrease in the amount of water vapor in the condenser. During adsorption phase, it can be seen that the condensate water is of about 0.6 kg for equal state and 0.57 kg for unequal state. There is a difference of 5% in condensate water between the two states. The results demonstrated that the condenser performance, for given operating conditions, is not strongly dependent on the adsorption steps.
Figure 8 shows the variation of evaporating temperature with the cycle time. A good agreement can be observed between the unequal model results and the experimental data. Results of the equal model show a different behavior with a lower evaporating temperature. In comparison to the unequal cycle, the proposed equal cycle is capable of producing chilled water at a temperature below that of the unequal cycle. For the unequal cycle, producing a cooling effect such as that achieved via the proposed equal cycle is almost impossible in practice. However, the proposed cycle mode enables generating a cooling effect and cold chilled water in either the adsorption or desorption mode. Thus, in each mode of the cycle, there is at least one bed is linked to the evaporator so that the cooling effect can be interminably produced. Additionally, it can be seen that during a period of 40 min of equal state, the evaporating temperature increases from 9 to 12.2 °C and to 14 °C for unequal state when the pre-cooling phase extended to 50 min. This difference in temperature does not seem very big (1.6 °C) although there is a difference in time of adsorption, 10 min, between the two cycles. Furthermore, to keep the temperature of refrigerant inside the evaporator as low as possible, the evaporator was insulated tightly. Additionally, the results showed that equal adsorption capacity decreased the evaporator temperature to 6.6 °C, while the equal adsorption capacity resulted in a lower evaporator temperature of 5.8 °C. This indicates that the system is working with a good statement for both capacities and a high vacuum is generated due to the high performance of the adsorption process.
Figure 9 shows that the equal cycle produces a higher SCP as compared to that of the unequal cycle. According to this figure the SCP rises sharply at the beginning of the adsorption phase for both cycles. Afterward, this SCP decreases gradually with the progress in the adsorption process. The possible reason is that at the starting of adsorption, the silica gel has the lowest adsorption capacity, so water can be adsorbed by silica gel easily. As the process of adsorption goes, adsorption amount declines, and so does the production of the cooling effect. As seen in this figure, the SCP peak value of equal cycle is 185 W/kg, which is higher than that of unequal cycle, 120 W/kg, there is 35% an increment in SCP value for the equal cycle, although the adsorption velocity for the former cycle is lower than the latter cycle. In contrast, the adsorption capacity rate (x) and the adsorption exchanged rate (∆x) are higher for the equal cycle. It can be noted that the adsorption velocity has a lower effect on the SCP than the amount of the adsorbed refrigerant. However, increasing the SCP leads to increase the cooling effect inside the evaporator, and this increases the performance of the system. Without doubt, the unequal cycle model should be accounted and not be neglected on designing cycle adsorption cooling system, especially when the cycle is varied because in real system, the adsorbed refrigerant could not be sufficient due to the non-equilibrium adsorption process.
Figure 10 has shown that the COP results of the equal adsorption process is similar in behavior to the unequal process and the reason is simply because of the high rate of the adsorption process. Those two figures give the dynamics of the system performance from the first time of adsorption phase up to the end of the phase. However, adsorption does not occur during the precooling phase while a decline in pressure exists, and the COP and SCP values are almost zero. Meanwhile, the COP value appears to almost increase linearly with the adsorption time as the cycle moves to the adsorption phase because the adsorption process is enhanced causing a higher cooling effect for a given heat input. The maximum COP value reached is 0.3 during the equal phase, and it is 0.15 during the unequal phase. Nonetheless, adequate performance can be achieved even for unequal cycle. In particular, the COP of 0.15 and SCP of 120 W/kg were estimated for the unequal cycle, demonstrating that it is still an acceptable value to work on and this system can be driven by a varied adsorption time. In general, there is a maximum of the COP at the equal cycle, but less adsorption velocity. Consequently, a strategy was adopted to assess the effect of the adsorption time. The different results are obtained for the COP and SCP: equaling the adsorption and desorption times leads to increase the system performance. For the two beds, equal time of adsorption and desorption phases is commonly used in cooling cycle but they are not always the same in the real system. A comparison with the results calculated by Angelo, et al. (2012) who employed similar conditions to those of our simulation (T
des = 90 °C, T
cw = 30 °C, T
e = 10 °C) shows that the peak value of the SCP is about 390 W/kg after 10 min, while for the present work, the peak value is 120 W/kg after 20 min. This difference in values is due to the different properties of the silica gel used in the two experiments. Angelo et al. used a “modified silica gel” which has an adsorption capacity rate exchange of ∆x > 0.1 kg water/kg silica gel, whereas in our work a silica gel of ∆x = 0.06 kg water/kg silica gel is used. Additionally, the maximum value of the COP (0.15) of our work is lower than that of Angelo’s experiment (0.41). The reason is that the larger amount of heat is required to raise the larger amount of silica gel for our experiment (10 kg).