Numerical Investigation of Small-Scale Adsorption Cooling System Performance Employing Activated Carbon-Ethanol Pair

Adsorber heat exchanger design has great importance in increasing the performance of the adsorption-based cooling system. In this study, a transient two-dimensional axisymmetric Computational Fluid Dynamics (CFD) model has been developed for the performance investigation of finned tube type adsorber using activated carbon and ethanol as the working pair. The operating conditions of the cooling system were 15, 20 and 80 for evaporation, cooling and heating temperatures, respectively. The simulated temperature profiles for different adsorbent thicknesses were validated with those from experimental data measured in our laboratory. Moreover, the error in mass and energy balance were 3% and 7.88%, respectively. Besides, the performance investigation has been performed for cycle time ranging from 600 s to 1400 s. The optimum cycle time was 800 s and the corresponding evaluated specific cooling power (SCP) and coefficient of performance (COP) were found to be 488 W/kg and 0.61, respectively. The developed CFD model will be used for fin height and fin pitch optimization and can be extended to other adsorbent-adsorbate based adsorption cooling system.


Introduction
Refrigeration is an important sector for the growing and sustainable global economy. It contributes significantly to food and drug conservation, thermal comfort, etc. However, greenhouse gas (GHG) emissions from this sector represent 7.8% of total emissions. Conventional refrigeration system employs refrigerants with high global warming potential (GWP) such as Chlorofluorocarbons (CFCs), hydrochlorofluorocarbons (HCFCs) and hydrofluorocarbons (HFCs) gases which produced 37% of the emissions and the remaining 63% was caused indirectly by the energy consumption of cooling and air conditioning systems [1].
In this context, adsorption cooling systems (ACS) are a feasible alternative to vapor compression cooling systems. In fact, they are driven by low-grade waste heat (below 100 • C) or solar thermal energy [2,3]. Moreover, they employ natural and environment-friendly refrigerants such as water [4,5], ethanol [6], methanol [7,8], ammonia [9], CO 2 [10], etc. Therefore, the widespread development be used to optimize the operating conditions and adsorber heat exchanger configurations such as fin pitch, fin height and fin thickness. Moreover, this developed model can be extended for other adsorbent-adsorbate pairs. Figure 1 shows the schematic diagram of a one-bed adsorption cooling system. It generally consists of four components: a condenser, an evaporator, an adsorber/desorber bed filled with the activated carbon powder (ACP) and an expansion valve. An adsorption cycle comprises four processes, which are pre-heating, desorption, pre-cooling and adsorption. In the pre-heating process, the bed is disconnected from the evaporator and hot water is circulated in the adsorption bed to increase the bed pressure. When bed pressure becomes higher than that in the condenser, then it is connected to the condenser and desorption process starts. Ethanol vapor is desorbed from activated carbon bed to the bed because of additional heating and moves to the condenser where it condenses. After the desorption is finished, the bed is disconnected from the condenser and cooling water is circulated in the adsorber bed which is called pre-cooling process. Due to the cooling water supply, the pressure in the bed decreases and becomes lower than that in the evaporator. At this moment, the adsorber bed is connected to the evaporator and evaporated ethanol vapor comes out of the adsorber and is adsorbed in the activated carbon bed. During the adsorption process, cooling water is supplied to the bed continuously to remove the heat of adsorption. from 600 s to 1400 s. The optimum cycle time was reported 800 s and the corresponding evaluated SCP and COP were found to be 488 W/kg and 0.61, respectively. Besides cycle time, our developed CFD model can be used to optimize the operating conditions and adsorber heat exchanger configurations such as fin pitch, fin height and fin thickness. Moreover, this developed model can be extended for other adsorbent-adsorbate pairs. Figure 1 shows the schematic diagram of a one-bed adsorption cooling system. It generally consists of four components: a condenser, an evaporator, an adsorber/desorber bed filled with the activated carbon powder (ACP) and an expansion valve. An adsorption cycle comprises four processes, which are pre-heating, desorption, pre-cooling and adsorption. In the pre-heating process, the bed is disconnected from the evaporator and hot water is circulated in the adsorption bed to increase the bed pressure. When bed pressure becomes higher than that in the condenser, then it is connected to the condenser and desorption process starts. Ethanol vapor is desorbed from activated carbon bed to the bed because of additional heating and moves to the condenser where it condenses. After the desorption is finished, the bed is disconnected from the condenser and cooling water is circulated in the adsorber bed which is called pre-cooling process. Due to the cooling water supply, the pressure in the bed decreases and becomes lower than that in the evaporator. At this moment, the adsorber bed is connected to the evaporator and evaporated ethanol vapor comes out of the adsorber and is adsorbed in the activated carbon bed. During the adsorption process, cooling water is supplied to the bed continuously to remove the heat of adsorption.

CFD Modeling
Performance investigation of finned tube type adsorber employing activated carbon and ethanol pairs was carried out experimentally in our laboratory earlier [25,26]. The experimental apparatus consists of an adsorber, an evaporator and a condenser connected with stainless steel tubes. The temperatures and pressures in all components are measured by K-type thermocouples and pressure gages, respectively. Two copper finned tubes were incorporated inside the adsorber. The finned tube heat exchanger consists of 190 circular fins attached to an annular tube having a length of 700 mm. Fin height, fin pitch and fin thickness were 10, 3.7 and 0.53 mm, respectively.

CFD Modeling
Performance investigation of finned tube type adsorber employing activated carbon and ethanol pairs was carried out experimentally in our laboratory earlier [25,26]. The experimental apparatus consists of an adsorber, an evaporator and a condenser connected with stainless steel tubes. The temperatures and pressures in all components are measured by K-type thermocouples and pressure gages, respectively. Two copper finned tubes were incorporated inside the adsorber. The finned tube heat exchanger consists of 190 circular fins attached to an annular tube having a length of 700 mm. Fin height, fin pitch and fin thickness were 10, 3.7 and 0.53 mm, respectively. Porous activated carbon was filled between two consecutive fins and sealed with a mesh sheet. The annular tube was 1.5 mm thick with an inner diameter of 26 mm. Ethanol was considered as refrigerant. Besides, water was used as a heat transfer fluid and passed through the annular tube. The adsorbent temperatures at 1 and 5 mm thicknesses were recorded experimentally by using K-type thermocouples shown as T 1 and T 5 in Figure 2.
Under similar experimental conditions, the transient CFD simulation was performed using Ansys-Fluent software v.18.1 (ANSYS, Canonsburg, PA, USA). Geometry and meshing were created by using Ansys Design Modeler and Ansys Meshing, respectively. The details of the CFD modeling were presented as follows:

Assumptions
The assumptions used in the simulation are as follows:

•
The adsorbent layer height is equal to the height of the fins.

•
The porous media is considered as homogenous.

•
Darcy model is adopted for flow through porous media. • A thermal equilibrium model is assimilated for porous media i.e., the adsorbent and ethanol vapor are at same temperature.

•
The variation of the heat of adsorption with the ethanol uptake is not considered, therefore, the average value of the heat of adsorption is used. Figure 2 shows the finned tube geometry and computational domain. The finned tube heat exchanger has a symmetric axis; therefore, the geometry is reduced to 2D. Besides, the computational domain used in the simulation is the half space between the two fins as this geometry has several symmetry planes. Porous activated carbon was filled between two consecutive fins and sealed with a mesh sheet. The annular tube was 1.5 mm thick with an inner diameter of 26 mm. Ethanol was considered as refrigerant. Besides, water was used as a heat transfer fluid and passed through the annular tube.

Geometry and Meshing
The adsorbent temperatures at 1 and 5 mm thicknesses were recorded experimentally by using K-type thermocouples shown as T1 and T5 in Figure 2. Under similar experimental conditions, the transient CFD simulation was performed using Ansys-Fluent software v.18.1 (ANSYS, Canonsburg, PA, USA). Geometry and meshing were created by using Ansys Design Modeler and Ansys Meshing, respectively. The details of the CFD modeling were presented as follows:

Assumptions
The assumptions used in the simulation are as follows:


The adsorbent layer height is equal to the height of the fins.  The porous media is considered as homogenous.  Darcy model is adopted for flow through porous media.  A thermal equilibrium model is assimilated for porous media i.e., the adsorbent and ethanol vapor are at same temperature.  The variation of the heat of adsorption with the ethanol uptake is not considered, therefore, the average value of the heat of adsorption is used. Figure 2 shows the finned tube geometry and computational domain. The finned tube heat exchanger has a symmetric axis; therefore, the geometry is reduced to 2D. Besides, the computational domain used in the simulation is the half space between the two fins as this geometry has several symmetry planes.  To perform the simulation, we used fine meshing and divided the computational domain (11.5 mm × 1.85 mm) into 748 elements with 897 nodes. The following mesh quality is reported from Fluent:

Geometry and Meshing
-Minimum orthogonal quality = 0.9721 where orthogonal quality ranges from 0 to 1 and values close to 0 represents low quality. -Maximum ortho skew = 0.0279 where ortho skew ranges from 0 to 1 and values close to 1 correspond to low quality.

Materials
The materials used in the simulations were: (i) copper tubing for tubes and fins, (ii) ethanol as gas phase refrigerant, and (iii) activated carbon powder of type Maxsorb III packed between the fins. The properties of copper were taken from the Fluent database. Real gas properties of ethanol were imported into Fluent from NIST Refprop database (v.9.1, NIST, Gaithersburg, MD, USA). The properties of Maxsorb III were, the particle density ρ p = 464.1 kg·m −3 , adsorbent-adsorbate thermal conductivity k p = 0.2 W/(m·K) [23] and specific heat capacity. The solid adsorbent heat capacity of Maxsorb III was investigated lately and the average value for temperatures ranging from 30 to 80 • C is 0.9 kJ/(kg·K) [27]. Therefore, the particle heat capacity for average instantaneous uptake of 0.7 g/g considering that the adsorbate is in the liquid phase [28] becomes 4.64 kJ/(kg·K) .

Porous Zone Properties
The permeability (α) and inertial loss coefficient (C) of porous media were calculated for the average particle diameter of Maxsorb III D P = 70 µm and adsorbent porosity γ = 0.434 [23] by the following equations [29]. The viscous flow resistance inside the porous media is the reciprocal of permeability of adsorbent: 2 (1)

Boundary Conditions
The boundary conditions applied in this simulation are shown in Figure 3 and the detailed boundary conditions are explained in Sections 3.4.1 and 3.4.2.

Pressure Inlet/Wall/Pressure Outlet Conditions
The interface between activated carbon and refrigerant is set as pressure inlet during adsorption process and the refrigerant inlet pressure was set to 3.85 kPa corresponding to evaporation temperature of 13.1 • C. Similarly, the interface is considered as pressure outlet during the desorption process. In addition, desorption is performed at a pressure of 10.35 kPa corresponding to the condensation temperature of 29.8 • C. Besides, at the time of preheating and precooling processes, the interface is considered as wall condition. The purpose of precooling is to reduce the bed pressure from condensation pressure to that of evaporation pressure. Therefore, when the bed pressure becomes lower than the evaporator's pressure, the precooling process was ended. Similarly, when the bed pressure becomes higher than that of condenser, the preheating process was terminated.

Convection Boundary Conditions
We applied convection boundary condition at the tube inner surface. The water temperature was estimated from the average of inlet and outlet temperature of the tube as presented in Figure 4. The free stream temperature profile was fitted to polynomial functions and implemented into fluent as a user-defined function (UDF). Moreover, Gnielinski correlation [30] is used to calculate the convection heat transfer coefficient for water flow rate of 3 L/min. Then, the heat transfer coefficient is incorporated into fluent as a UDF for the temperature range 20 • C to 80 • C. The heating and cooling processes were applied for 400 s each.

Convection Boundary Conditions
We applied convection boundary condition at the tube inner surface. The water temperature was estimated from the average of inlet and outlet temperature of the tube as presented in Figure 4. The free stream temperature profile was fitted to polynomial functions and implemented into fluent as a user-defined function (UDF). Moreover, Gnielinski correlation [30] is used to calculate the convection heat transfer coefficient for water flow rate of 3 L/min. Then, the heat transfer coefficient is incorporated into fluent as a UDF for the temperature range 20 °C to 80 °C. The heating and cooling processes were applied for 400 s each.

Convection Boundary Conditions
We applied convection boundary condition at the tube inner surface. The water temperature was estimated from the average of inlet and outlet temperature of the tube as presented in Figure 4. The free stream temperature profile was fitted to polynomial functions and implemented into fluent as a user-defined function (UDF). Moreover, Gnielinski correlation [30] is used to calculate the convection heat transfer coefficient for water flow rate of 3 L/min. Then, the heat transfer coefficient is incorporated into fluent as a UDF for the temperature range 20 °C to 80 °C. The heating and cooling processes were applied for 400 s each.

Governing Equations
The main mathematical equations used in fluent are mass, momentum and energy conservation equations. To consider the adsorption phenomena, mass and energy source terms are added to these conservation equations through user-defined functions (UDFs).

Mass Conservation Equation in Porous Media
The mass conservation equation for porous media is given by: where the mass source term (S m ) denote the adsorbed ethanol gas inside the porous media and is specified as follows: The term dq dt represents the adsorption rate due to the diffusion of ethanol vapor through the Maxsorb III micro-pores which will be stated in Section 3.5.4.

Momentum Conservation Equation
The momentum conservation equation for porous media is expressed by: where the momentum source term (F) is estimated for the viscous and inertial losses in porous media.
The viscous resistance (inverse of permeability: 1/α) and inertial losses coefficients were defined previously in Section 3.3.2:

Energy Conservation Equation
The energy conservation equation for porous media considering thermal equilibrium condition is given by: S h is the heat source term corresponding to the heat released by adsorption process and is given by: Here, Q st = 1002 kJ/kg is the average heat of adsorption for Maxsorb III-ethanol pair [31].

Adsorption Characteristics
The linear driving force (LDF) model [32] is used in this study to predict the adsorption kinetics ( dq dt ) and it is expressed by: where k and q * are the diffusion time constant and equilibrium uptake, respectively. The diffusion time constant is defined by the Arrhenius equation (Equation (10)) and the equilibrium uptake was calculated by Dubinin-Astakhov (D-A) adsorption isotherm equation (Equation (11)): A = 0.2415 s −1 , E a = 225 kJ/kg, q s = 1.2 kg/kg, E = 139.5 kJ/kg and n = 1.8 denote the pre-exponential factor, activation energy, saturated uptake, characteristic energy and heterogeneity parameter, respectively [33].
In D-A equation, P s denotes the saturated pressure (bar) calculated by the Antoine equation: where

Performance Investigation
The performance of adsorption cooling system was investigated using Equations (13) and (14) to calculate specific cooling power (SCP) and coefficient of performance (COP), respectively: . m ads dt t cycle ·m ac (13) where:

Mass and Energy Balance Equation
The mass balance in the adsorber/desorber bed implies that the mass of refrigerant adsorbed onto the adsorbent is equal to that desorbed from the bed as described in Equation (16) The energy balance in the adsorption cooling system indicates that the heat released by the condenser and adsorber bed is equal to that adsorbed by the evaporator and desorber bed as expressed by the following equation: Figure 5 shows the simulated pressure change vs. experimental pressure profile in the adsorber/desorber bed for two cycles. The experimental pressure was not constant during adsorption and desorption processes as it was affected by the temperature fluctuation in condenser and evaporator, respectively and this was not considered in this study. In the simulated pressure profile, adsorption and desorption occur at constant pressure condition as seen from Figure 5. Besides, the simulation can obtain the cyclic steady-state pressure within two cycles. Good agreement has been found between experimental and simulated pressure profiles at the beginning of pre-heating and pre-cooling processes. adsorber/desorber bed for two cycles. The experimental pressure was not constant during adsorption and desorption processes as it was affected by the temperature fluctuation in condenser and evaporator, respectively and this was not considered in this study. In the simulated pressure profile, adsorption and desorption occur at constant pressure condition as seen from Figure 5. Besides, the simulation can obtain the cyclic steady-state pressure within two cycles. Good agreement has been found between experimental and simulated pressure profiles at the beginning of pre-heating and pre-cooling processes.

Temperature Profiles
Experimental temperature change vs. simulated temperature profile at 1 and 5 mm adsorbent thicknesses (T1 and T5) were presented in Figure 6. Good agreement was found between the experimental data and the simulated results. From Figure 6, simulated temperature profile perfectly matched with that of experimental data at the end of desorption and in pre-cooling process. However, a little deviation is found during the adsorption process and at the starting of desorption. One possible reason for deviation at adsorption process is that Gnielinski correlation under-estimated the heat transfer coefficient from the adsorber bed to the heat transfer water.

Temperature Profiles
Experimental temperature change vs. simulated temperature profile at 1 and 5 mm adsorbent thicknesses (T 1 and T 5 ) were presented in Figure 6. Good agreement was found between the experimental data and the simulated results. From Figure 6, simulated temperature profile perfectly matched with that of experimental data at the end of desorption and in pre-cooling process. However, a little deviation is found during the adsorption process and at the starting of desorption. One possible reason for deviation at adsorption process is that Gnielinski correlation under-estimated the heat transfer coefficient from the adsorber bed to the heat transfer water.   Figure 7 provides the simulated temperature profile at different adsorbent thicknesses from the tube outer surface such as 0, 1, 5 and 8 mm. During pre-heating and pre-cooling processes, the change of temperature is almost same in all these points due to high heat transfer rate. However, during desorption and adsorption process, the change of temperature is noticeable at 0, 1 and 5 mm thickness.    Figure 7 provides the simulated temperature profile at different adsorbent thicknesses from the tube outer surface such as 0, 1, 5 and 8 mm. During pre-heating and pre-cooling processes, the change of temperature is almost same in all these points due to high heat transfer rate. However, during desorption and adsorption process, the change of temperature is noticeable at 0, 1 and 5 mm thickness. Figure 7 provides the simulated temperature profile at different adsorbent thicknesses from the tube outer surface such as 0, 1, 5 and 8 mm. During pre-heating and pre-cooling processes, the change of temperature is almost same in all these points due to high heat transfer rate. However, during desorption and adsorption process, the change of temperature is noticeable at 0, 1 and 5 mm thickness.

Adsorption Characteristics
The profiles of equilibrium uptake ( * ) and instantaneous uptake ( ) in the adsorber/desorber bed are displayed in Figure 8. The difference between maximum and minimum values of equilibrium uptake ∆ * is 0.51 kg/kg whereas the change in case of instant uptake ∆ is 0.314 kg/kg. This means that ethanol adsorption onto activated carbon is only 61.6% of its capacity within the 400 s adsorption time.

Adsorption Characteristics
The profiles of equilibrium uptake (q * ) and instantaneous uptake (q) in the adsorber/desorber bed are displayed in Figure 8. The difference between maximum and minimum values of equilibrium uptake ∆q * is 0.51 kg/kg whereas the change in case of instant uptake ∆q is 0.314 kg/kg. This means that ethanol adsorption onto activated carbon is only 61.6% of its capacity within the 400 s adsorption time.

Energy and Mass Balance
The total amount of heat removed from adsorber and that added to desorber were calculated by integrating the heat transfer profile as shown in Figure 9. The values of heat released from condenser and adsorber bed and those adsorbed by desorber and evaporator are used to check the energy balance as shown in Equation (17) and the error in energy balance is found to be 7.9%.
The simulated ethanol flow rate to/from the adsorber/desorber bed is presented in Figure 10. The error in mass balance was 3% which is estimated by Equation (16).
One of the main reasons for mass balance error is that the time for adsorption and desorption was not exactly same in the case of the experiment. The time difference during the experiment was around 12 s. This mass balance error also affects the energy balances. Besides, during adsorption at 20 °C temperature, the heat transfer coefficient was underestimated by Gnielinski correlation that is

Energy and Mass Balance
The total amount of heat removed from adsorber and that added to desorber were calculated by integrating the heat transfer profile as shown in Figure 9. The values of heat released from condenser and adsorber bed and those adsorbed by desorber and evaporator are used to check the energy balance as shown in Equation (17) and the error in energy balance is found to be 7.9%.
The simulated ethanol flow rate to/from the adsorber/desorber bed is presented in Figure 10. The error in mass balance was 3% which is estimated by Equation (16).
One of the main reasons for mass balance error is that the time for adsorption and desorption was not exactly same in the case of the experiment. The time difference during the experiment was around 12 s. This mass balance error also affects the energy balances. Besides, during adsorption at 20 • C temperature, the heat transfer coefficient was underestimated by Gnielinski correlation that is also visible in Figure 6. Therefore, heat input was higher than the heat removed from the bed.

Energy and Mass Balance
The total amount of heat removed from adsorber and that added to desorber were calculated by integrating the heat transfer profile as shown in Figure 9. The values of heat released from condenser and adsorber bed and those adsorbed by desorber and evaporator are used to check the energy balance as shown in Equation (17) and the error in energy balance is found to be 7.9%.
The simulated ethanol flow rate to/from the adsorber/desorber bed is presented in Figure 10. The error in mass balance was 3% which is estimated by Equation (16).
One of the main reasons for mass balance error is that the time for adsorption and desorption was not exactly same in the case of the experiment. The time difference during the experiment was around 12 s. This mass balance error also affects the energy balances. Besides, during adsorption at 20 °C temperature, the heat transfer coefficient was underestimated by Gnielinski correlation that is also visible in Figure 6. Therefore, heat input was higher than the heat removed from the bed.

Performance Investigation
The performance of adsorption cooling system (ACS) has been evaluated in terms of specific cooling power (SCP) and coefficient of performance (COP) by using Equations (13) and (14), correspondingly. The estimated SCP and COP for 800 s cycle time are found 488 W/kgac and 0.61, respectively. This performance is also relevant to the reported performance for activated carbon-ethanol pair in other studies previously. The specific cooling capacity for activated-carbon ethanol pair considering two different particle size and three different domains has been reported between the range of 424 and 710 W/kgac by Sourav et al. [24]. Besides, the performance of activated carbon fiber-ethanol based adsorption cooling system was found 200 W/kg for SCP and 0.65 for COP previously [17]. Figure 11 shows the performance variation with cycle time. The COP increases with the increase of cycle time. However, the SCP reaches its maximum at a cycle time of 800 s. With

Performance Investigation
The performance of adsorption cooling system (ACS) has been evaluated in terms of specific cooling power (SCP) and coefficient of performance (COP) by using Equations (13) and (14), correspondingly. The estimated SCP and COP for 800 s cycle time are found 488 W/kg ac and 0.61, respectively. This performance is also relevant to the reported performance for activated carbon-ethanol pair in other studies previously. The specific cooling capacity for activated-carbon ethanol pair considering two different particle size and three different domains has been reported between the range of 424 and 710 W/kg ac by Sourav et al. [24]. Besides, the performance of activated carbon fiber-ethanol based adsorption cooling system was found 200 W/kg for SCP and 0.65 for COP previously [17]. Figure 11 shows the performance variation with cycle time. The COP increases with the increase of cycle time. However, the SCP reaches its maximum at a cycle time of 800 s. With adsorption time increase, the cooling effect increases as the total mass adsorbed increases. However, the flow rate of adsorbed mass increases sharply at starting of adsorption then decreases as shown previously in Figure 10. Therefore the cooling power increases with cycle time increase and reaches its maximum at optimum cycle time of 800 s. Then, for higher adsorption/desorption time, the cooling power decreases. The COP increases with adsorption/desorption time increase which as the increase in cooling effect generated is higher than the increase of heat consumed to desorb the refrigerant from activated carbon.

Conclusions
CFD modeling of a finned tube adsorber employing activated carbon-ethanol as the adsorbent-adsorbate pair has been performed. The simulated pressure and temperature profiles were compared with experimental data of an adsorption cooling system tested in our laboratory for evaporation temperature of 15 °C, cooling temperature of 20 °C, the heating temperature of 80 °C and total cycle time of 800 s. Good agreement was found between experimental and simulation results. The performance of the system was investigated numerically for different adsorption/desorption phase times ranging from 300 to 700 s. The optimal cycle time for pre-heating and adsorption phases as well as pre-cooling and adsorption phases corresponded to 400 s. At this condition, the cooling power reached 488 W per kilogram of adsorbent and the COP was 0.61. This CFD model will allow performance optimization of the system by optimizing the operating conditions as well as fin dimensions.
Author Contributions: M.K. and S.J. have equal contributions in CFD modeling and paper writing. T.M. has provided supervision and guidance in the simulation. B.B.S. and S.K. delivered overall guidance and were involved in the discussion.

Conflicts of Interest:
The founding sponsors had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Conclusions
CFD modeling of a finned tube adsorber employing activated carbon-ethanol as the adsorbent-adsorbate pair has been performed. The simulated pressure and temperature profiles were compared with experimental data of an adsorption cooling system tested in our laboratory for evaporation temperature of 15 • C, cooling temperature of 20 • C, the heating temperature of 80 • C and total cycle time of 800 s. Good agreement was found between experimental and simulation results. The performance of the system was investigated numerically for different adsorption/desorption phase times ranging from 300 to 700 s. The optimal cycle time for pre-heating and adsorption phases as well as pre-cooling and adsorption phases corresponded to 400 s. At this condition, the cooling power reached 488 W per kilogram of adsorbent and the COP was 0.61. This CFD model will allow performance optimization of the system by optimizing the operating conditions as well as fin dimensions.