Dynamic Evaluation of Desiccant Dehumidiﬁcation Evaporative Cooling Options for Greenhouse Air-Conditioning Application in Multan (Pakistan)

: This study provides insights into the feasibility of a desiccant dehumidiﬁcation-based Maisotsenko cycle evaporative cooling (M-DAC) system for greenhouse air-conditioning application. Conventional cooling techniques include direct evaporative cooling, refrigeration systems, and pas-sive/active ventilation. which are commonly used in Pakistan; however, they are either not feasible due to their energy cost, or they cannot efﬁciently provide an optimum microclimate depending on the regions, the growing seasons, and the crop being cultivated. The M-DAC system was therefore proposed and evaluated as an alternative solution for air conditioning to achieve optimum levels of vapor pressure deﬁcit (VPD) for greenhouse crop production. The objective of this study was to investigate the thermodynamic performance of the proposed system from the viewpoints of the temperature gradient, relative humidity level, VPD, and dehumidiﬁcation gradient. Results showed that the standalone desiccant air-conditioning (DAC) system created maximum dehumidiﬁcation gradient (i.e., 16.8 g/kg) and maximum temperature gradient (i.e., 8.4 ◦ C) at 24.3 g/kg and 38.6 ◦ C ambient air conditions, respectively. The DAC coupled with a heat exchanger (DAC+HX) created a temperature gradient nearly equal to ambient air conditions, which is not in the optimal range for greenhouse growing conditions. Analysis of the M-DAC system showed that a maximum air temperature gradient, i.e., 21.9 ◦ C at 39.2 ◦ C ambient air condition, can be achieved, and is considered optimal for most greenhouse crops. Results were validated with two microclimate models (OptDeg and Cft) by taking into account the optimality of VPD at different growth stages of tomato plants. This study suggests that the M-DAC system is a feasible method to be considered as an efﬁcient solution for greenhouse air-conditioning under the climate conditions of Multan (Pakistan). which lies within the required optimum temperature conditions of the air-conditioning. Additionally, the M-DAC created a dehumidiﬁcation gradient 16.8 g/kg) at 24 g/kg ambient conditions, which also lies inside the range of optimum humidity conditions inside the greenhouse. Moreover, the vapor pressure deﬁcit (VPD) proﬁle of the greenhouse environment related to the growth of plants indicated that a VPD of 0.8–1.2 kPa (for leaf temperature equal to surrounding air temperature) is most suitable for healthy transpiration from the plant leaves at late vegetation and early ﬂowering stages. These results were validated with OptDeg and Cft-ratio microclimate models, which showed that the M-DAC system maintained the VPD of the product air inside the healthy transpiration zone throughout most of the summer months (i.e., May to August). Moreover, the comfort ratio model results of the proposed M-DAC system were compared with experimental results from a tomato greenhouse, indicating that the M-DAC system could potentially achieve the required comfort level compared to a pad-and-fan type air-conditioning system. Thus, the study indicates the M-DAC system is a feasible option for greenhouse air-conditioning for the climatic conditions of Multan (Pakistan).


Introduction
Closed-field production of crops and vegetables by means of greenhouses with different covering materials depends on the efficiency of the climate control system to provide optimum growth conditions with a reasonable energy cost. The extensive solar radiation during the daytime causes a greenhouse effect, resulting in high air temperature and relative humidity that can significantly exceed the optimum conditions. This can lead to a significant loss of yield and fruit quality in the absence of a proper ventilation and dehumidification system. Additionally, the plants and vegetables are more prone to pest/fungus/disease attack in higher relative humidity conditions than normal, which may also affect the production level of the greenhouse. To reduce the impact of suboptimal microclimate conditions on production, passive and active ventilation, in addition to direct evaporative cooling systems by means of pad-and-fans or swamp coolers, are conventionally used by greenhouse growers in Pakistan. The problem with these solutions is that they do not allow control over the relative humidity level of the greenhouse environment, hence leaving the plants and vegetables exposed to pests, fungus, and disease attacks. Therefore, a more appropriate air-conditioning (AC) solution that is cost-efficient and can be implemented in large scale commercial production is required. In this regard, desiccant dehumidification-based indirect evaporative cooling air-conditioning systems are gaining research attention that needs to be evaluated for AC application in different closed-field plant production environments.
The focus of this study is the city of Multan, Pakistan, which lies in a warm desert climate where direct evaporative cooling methods are used in combination with passive and active ventilation methods inside agricultural greenhouses. For the purpose of this research, a desiccant dehumidification-based Maisotsenko cycle evaporative cooling (M-DAC) system is proposed and evaluated as an alternative solution for controlling humidity levels to achieve the optimum vapor pressure deficit (VPD). The objective was to provide a systematic thermodynamic analysis of the proposed M-DAC system from the viewpoints of temperature gradient, relative humidity levels, optimum VPD, and dehumidification gradient. A review of the relevant published literature is summarized in Section 2 to identify the most recent achievements and the gaps. Modeling of the M-DAC and dynamic evaluation by means of the optimality degree and comfort ratio models are described in Section 3. The proposed solution is validated and discussed in Section 4 based on the thermodynamic analysis and the microclimate models. Section 5 concludes the study and highlights the potential aspects for further research and development.

Background and Literature Review
Evaporative cooling (EC) systems can potentially be classified into direct evaporative (DEC), indirect evaporative (IEC), and Maisotsenko cycle evaporative (MEC) systems. The MEC system is widely studied by the research community for heating, ventilation, and air-conditioning (HVAC) system applications. Pandelidis and Anisimov [1] numerically studied and compared eight different configurations of a standalone Maisotsenko-cycle evaporative heat and mass exchanger. Caliskan et al. [2] studied energy and exergy analyses of the standalone MEC system for building air-conditioning. Zhan et al. [3] studied counter flow standalone MEC system for building air-conditioning application. Caliskan et al. [4] compared traditional vapor compression air-conditioning systems with the MEC system based on energy and exergy analyses. Chua et al. [5] reviewed environment-friendly airconditioning options for building applications. Pandelidis et al. [6] numerically simulated the MEC heat exchanger for different air-conditioning applications. Rogdakis et al. [7] numerically simulated and experimentally validated the MEC system for building airconditioning application for the climatic conditions of Greece. Cui et al. [8] studied the MEC system for precooling of ambient air as an energy-saving technique for hot and humid climatic conditions. Riangvilaikul and Kumar [9] numerically investigated the MEC Energies 2021, 14, 1097 3 of 21 system for different inlet conditions. Sultan [10] studied the MEC system for greenhouse air-conditioning and ventilation applications. Anisimov et al. [11] theoretically and experimentally analyzed the standalone MEC system for air-conditioning and found the MEC system to be an efficient way of indirect evaporative cooling. In addition, results from the study concluded that inlet air velocity impacted the performance of the system. Zhao et al. [12] numerically investigated the performance of the MEC system for climatic conditions of the United Kingdom, which resulted in wet-bulb effectiveness of 1.3. Zube and Gillan [13] experimentally studied a commercial type MEC air-conditioning system and evaluated heat mass transfer parameters inside a heat exchanger for the first time. Weerts [14] reported the power-saving potential of commercially available standalone MEC system compared to the vapor compression air-conditioning (VCAC) system. Maisotsenko and Treyger [15] also reported the energy-saving potential of the standalone MEC system compared to standalone VCAC system. Similarly, [16][17][18] also studied the standalone MEC system in-depth and concluded that the MEC system can potentially be an environmentally friendly air-conditioning option compared to the VCAC system.
The standalone desiccant air-conditioning (DAC) system has been extensively studied in the literature for various applications. Generally, the DAC system better suits regions with higher relative humidity in the ambient air. The DAC system was studied for the possible application of vehicular air-conditioning [19]. Additionally, the authors used a continuous input and output modulator (i.e., proportional-derivative (PD) controller) to increase the efficiency of the system and concluded that the DAC system coupled with the PD controller was energy conservative compared to conventional AC systems in vehicles [19]. Moreover, the DAC system was successfully installed in a wet market in Hong Kong. The authors numerically simulated the energy loads of a hypothetical wet-market using EnergyPlus for the climatic conditions of Hong Kong and compared results with the actual wet market [20]. The authors concluded that the DAC could potentially be used as a replacement for conventional air-conditioning systems in the wet markets. A solar-assisted open cycle DAC system was experimentally tested for grain storage in Melbourne, Australia [21]. The authors concluded that a 5.85 m 3 solar collector area was sufficient for cooling of up to 200 tonnes of grain. The coefficient of performance (COP) was a function of surrounding conditions and air mass flow rate. A maximum COP of 86.2 was observed under humid conditions for an air mass flow rate of 0.019 kg/s [21]. The standalone DAC system was studied for a greenhouse air temperature and humidity control system [22]. Moreover, the authors studied the greenhouse air-conditioning from a vapor pressure deficit (VPD) point of view. The standalone DAC system with dual desiccant wheels was studied for residual waste heat recovery in marine ships [23]. The authors concluded that the DAC system was 33.4% more power-efficient compared to traditional air-conditioning systems for marine ships. An EC-coupled DAC system was studied in two modes (i.e., ventilation and recirculation modes) [24]. The authors concluded that conventional ventilation mode is more efficient when humidity ratio is below~10.9 g/kg, whereas the recirculation mode of the air-conditioning system is more efficient for areas where humidity ratio is above~10.9 g/kg, consuming a higher amount of input energy. Enteria et al. [25] studied the analyses of the first and second law of thermodynamics for a solar thermal electric desiccant air-conditioning system [25]. A simulated solar-operated evaporative cooling-assisted solid desiccant air-conditioning system was developed in a tropical region of Malaysia [26]. The authors used the TRNSYS simulation environment to simulate the desiccant air-conditioning system for climatic conditions of Malaysia. The authors concluded that two-stage ventilation mode desiccant air-conditioning system is better for tropical conditions such as Malaysia, with a temperature gradient of 17.6 • C at 30 • C ambient temperature [26]. Another study investigated a dual-wheel cooling system [27]. Results from the study indicate that higher regeneration temperature (RT) produces higher total COP at both air mass flow rates (i.e., 100 and 200 m 3 /h) [27]. Moreover, the performance and energy-saving potential of a two-stage solar-driven rotary wheel DAC system was compared with a conventional vapor compression air-conditioning (VCAC) system for building air-conditioning [28]. Building energy performance simulation (BEPS) was conducted for Berlin and Shanghai. The DAC system was found be to energy-saving and more efficient for both of the cities compared to VCAC system. However, in the case of Shanghai, the RT of the DAC system was 30 • C higher than that of Berlin [28]. Figure 1 shows the schematic diagram and working principle of a conventional vapor compression air-conditioning system, and solar thermal liquid desiccant dehumidification-based airconditioning system. Moreover, the authors studied desiccant dehumidification systems and evaporative cooling systems in detail for various applications, including greenhouse air-conditioning [29][30][31][32], building air-conditioning [33], livestock air-conditioning [34], and agricultural product storage [35][36][37][38][39][40][41][42][43][44].
Song and Sobhani [45] studied solar-assisted desiccant air-conditioning system coupled with a Maisotsenko evaporative cooling system and phase change material for building air-conditioning application. The authors concluded that the maximum total COP of 0.404 was achieved using a solar-operated desiccant air-conditioning system coupled with a Maisotsenko evaporative cooler using a phase change material for air-conditioning of an official building in Iran. In another study, a Maisotsenko cycle-assisted desiccant air-conditioning system was studied for the climatic conditions of Japan [46]. Pandelidis et al. [47] numerically studied the Maisotsenko cycle-assisted desiccant air-conditioning system with two modifications, i.e., bypass and bypass with rotary sensible heat exchanger. The authors concluded that all the configurations successfully managed the sensible load of the environment while being highly dependent on the humidity ratio of ambient air. The addition of a rotary sensible heat exchanger in system configuration 3 increased the cooling capacity up to~2.35 kW at 40 • C ambient temperature generating~17.25 • C supply air temperature at the same ambient conditions [47]. Table 1 shows the summary of different desiccant-based Maisotsenko cycle evaporative cooling systems for air-conditioning applications. From Table 1 and the reviewed literature, it is evident that there is a gap in the literature regarding OptDeg and Cft plant thermal comfort indices for the proposed M-DAC system. Moreover, this study aims to use a low and easily achievable regeneration temperature (i.e., 50 • C) coupled with different comfort levels of plant growth stages, which is lacking in previous studies. The study area (i.e., Multan) lies in a warm desert climate of the Köppen-Geiger climatic classification ( Figure 2). Therefore, this justifies the need for an air-conditioning system for greenhouses. Usually, pad-type direct evaporative cooling or natural/forced ventilation systems are used for greenhouse air-conditioning. These systems are not feasible throughout the year and fail to achieve the required optimum temperature and relative humidity conditions inside a greenhouse. Moreover, the DEC system fails to handle humidity inside the greenhouse environment, which is crucial to the vapor pressure deficit. The desiccant-based hybrid evaporative cooling systems have not been explored for greenhouse air-conditioning applications in Pakistan.    bug Rhino v6 for the climatic conditions of Multan (Pakistan). A typical meteorological EnergyPlus weather file was used in Ladybug Rhino to represent the hourly weather conditions. In Figure 3, the light fill area shows the average hourly variation, whereas the lines show the average daily variation in climatic parameters. From the viewpoint of the literature reviewed above, pad-type direct evaporative cooling systems are only suitable for controlling the temperature, and fail to achieve the optimum humidity level inside the greenhouse. To overcome this issue, an indirect evaporative cooling system could be used, which provides a sensibly cooled product air with no latent transfer from the wet to dry channels of the IEC system. However, due to the low wet bulb effectiveness of the system, it cannot be considered a suitable option for achieving the required temperature level inside the greenhouse. Additionally, higher levels of humidity inside the greenhouse can potentially result in pest, fungus, and disease attacks on the plants and vegetables. Therefore, in this study, a desiccant dehumidification-based Maisotsenko cycle evaporative cooling air-conditioning (M-DAC) system was thermodynamically analyzed for the climatic conditions of Multan (Pakistan). The performance of the system was analyzed for optimum greenhouse air-conditioning conditions from the viewpoints of temperature gradient, relative humidity level, and wet bulb effectiveness of the systems. The study area (i.e., Multan) lies in a warm desert climate of the Köppen-Geiger climatic classification ( Figure 2). Therefore, this justifies the need for an air-conditioning system for greenhouses. Usually, pad-type direct evaporative cooling or natural/forced ventilation systems are used for greenhouse air-conditioning. These systems are not feasible throughout the year and fail to achieve the required optimum temperature and relative humidity conditions inside a greenhouse. Moreover, the DEC system fails to handle humidity inside the greenhouse environment, which is crucial to the vapor pressure deficit. The desiccant-based hybrid evaporative cooling systems have not been explored for greenhouse air-conditioning applications in Pakistan.  Figure 3 shows the ambient climatic conditions of Multan (Pakistan) and the greenhouse optimum temperature and relative humidity zone. Figure 3 was made using Ladybug Rhino v6 for the climatic conditions of Multan (Pakistan). A typical meteorological EnergyPlus weather file was used in Ladybug Rhino to represent the hourly weather conditions. In Figure 3, the light fill area shows the average hourly variation, whereas the potentially result in pest, fungus, and disease attacks on the plants and vegetables. Therefore, in this study, a desiccant dehumidification-based Maisotsenko cycle evaporative cooling air-conditioning (M-DAC) system was thermodynamically analyzed for the climatic conditions of Multan (Pakistan). The performance of the system was analyzed for optimum greenhouse air-conditioning conditions from the viewpoints of temperature gradient, relative humidity level, and wet bulb effectiveness of the systems.

Modelling M-DAC System
A general working principle of the M-DAC system is first shown in Figure 4, followed by a schematic representation in Figure 5 that demonstrates the application of the proposed system for greenhouse air-conditioning for a possible higher yield. Taking into account the climatic conditions of Multan (Pakistan), the performance of a silica gel-based desiccant wheel for M-DAC system was investigated using Equations (1)-(4), as described in the literature [52,53].
where F 1,ip and F 2,ip represent the combined potential as a function of the humidity ratio in kg/kg and temperature in • C of the solid desiccant-based air-conditioning system. ηF 1 and ηF 2 represent the efficiencies of the system correlating to the combined potentials. ηF 1 and ηF 2 have typical values of 0.05 and 0.95 for a high-efficiency silica gel-based solid desiccant wheel air-conditioning system. The coefficients A, B, and C are presented in Table 2.

Modelling M-DAC System
A general working principle of the M-DAC system is first shown in Figure 4, followed by a schematic representation in Figure 5 that demonstrates the application of the proposed system for greenhouse air-conditioning for a possible higher yield. Taking into account the climatic conditions of Multan (Pakistan), the performance of a silica gel-based desiccant wheel for M-DAC system was investigated using Equations (1)-(4), as described in the literature [52,53].

Modelling M-DAC System
A general working principle of the M-DAC system is first shown in Figure 4, followed by a schematic representation in Figure 5 that demonstrates the application of the proposed system for greenhouse air-conditioning for a possible higher yield. Taking into account the climatic conditions of Multan (Pakistan), the performance of a silica gel-based desiccant wheel for M-DAC system was investigated using Equations (1)-(4), as described in the literature [52,53].   In comparison with the conventional vapor compression air-conditioning systems, the desiccant dehumidification-assisted evaporative cooling air-conditioning system consumes less energy due to the absence of compressor work and is able to operate at lower regenera-Energies 2021, 14, 1097 9 of 21 tion temperatures (i.e., 50-90 • C), which are easily attainable using solar thermal, biogas, and waste heat sources. The Maisotsenko cycle evaporative cooling system was modeled using the model of [32]. Equation (5) shows the performance model of the Maisotsenko cycle evaporative cooling system.  (1)-(4)).

Coefficient Value
where T in represents the inlet temperature of the process air in • C and w in represents the inlet humidity ratio of the process air in kg/kg. Equation (5) models the temperature of product air produced from the Maisotsenko cycle evaporative cooling system. The model was developed by [32] using experimental data of a developed Maisotsenko cycle based evaporative cooling system [11]. The sensible heat exchanger (HX) used in this study was modeled using the heat exchanger standard equation (Equation (6)) provided in the literature by The American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) [54].
where T 3 represents the outlet temperature of the sensibly cooled process air resulting from the heat exchanger, T 2 represents the outlet temperature of heated and dehumidified air resulted from the desiccant air-conditioning system, ε HX represents the sensible heat exchange efficiency of the heat exchanger (assumed as 0.9), and T 1 represents the dry bulb temperature of the ambient air inlet into the desiccant air-conditioning system. Figure 5 shows the working principle and schematic of the sensible heat exchanger. Regeneration of the desiccant is done usually through regeneration air at higher than ambient temperature. To heat the regeneration air, the air passes through the heat exchanger which sensibly heats the air to a certain limit. For further heating of the regeneration air before entering the desiccant, a heating unit/source is used. This heating source could be a heating electric coil, solar heater, or biogas-operated thermal heat source. Waste heat from a condenser or any other source in industrial uses can also be used to regenerate the desiccant material and absorb its moisture. For this study, a regeneration temperature of 50 • C was used. This relatively lower temperature is achievable through a solar thermal system, biogas-operated thermal heat source, or waste heat. A flat plate solar collector and heat exchanger were used for this study to provide the required regeneration temperature (i.e., 50 • C). Figure 5 shows the working principle and schematic of a solar thermal heating source.

Optimality Degree and Comfort Ratio Model
Dynamic assessment of the proposed M-DAC system for greenhouse crop production application was carried out by means of two microclimate models, namely, OptDeg [55] and Cft-ratio [56] that take into account VPD measurements as inputs. The output of the first model, optimality degree of VPD, denoted by Opt(VPD) GS, (Light) = α : VPD → [0, 1] is a quantitative value between 0 and 1 that represents how close a VPD reading inside the greenhouse is to the optimum value (0 ≤ α ≤ 1) at a specific growth stage (GS) and light condition. Graphical representation of the membership functions for defining Opt(VPD) under different conditions is provided in Figure 6. In this model, a VPD measurement in the greenhouse is first calculated from Equation (7) at the time t m,n (where m and n refer to a specific minute and date of the growing stage), and is then mapped to a value between 0 and 1 that quantifies the optimality degree of the microclimate for greenhouse production. In other words, an optimality-degree equal to 1 refers to a high potential yield with marketable value and high-quality fruit.
VPD t m,n = f T t m,n , RH t m,n , τ t m,n = C g 10 (7.5τ tm,n )/(237.3+τ tm,n ) 1000 − C g 10 (7.5T tm,n )/(237.3+T tm,n ) 1000 RH t m,n 100 Energies 2021, 14, x FOR PEER REVIEW 11 of 23 The OptDeg model determines the variation of a single VPD measurement from the optimal references in an instant time. Although this model can provide an overview of the performance of the M-DAC system in sample time, a different model called comfort ratio, denoted by Cft-ratio, that incorporates a sample time frame, such as 24 h was used. The comfort ratio of VPD, denoted by ( , , ) = , represents the percent of vapor pressure deficit at a specific time and growth stage, which lies between the reference conditions of vapor pressure deficit related to . Ideal microclimate conditions inside a greenhouse could potentially be defined as ( , , 1) = 1. Self-adjusted optimality degree is denoted by to define the reference conditions required for performance evaluation of the proposed M-DAC system or control for microclimate evaluation. The set of Simulink blocks that were implemented to determine the Cft-ratio of the VPD values for different generated from the M-DAC system are shown in Figure 7. For this research, we selected three levels of as 0, 0.5, and 1, representing failure, marginal, and ideal production. It is worth mentioning that OptDeg and Cft-ratio models were derived with a prime focus on improved fruit quality and yield. The derivations at the back of the models (which can be found in published literature) for greenhouse tomato production were condensed for the sake of simplicity in this study.  It should be noted that there are various functions for calculating VPD that take into accounts different assumptions and constants; however, that presented in Equation (7) is widely accepted and used by meteorologists and commercial greenhouse growers [34]. The first part of Equation (7) represents the vapor pressure of the saturated environment (i.e., the leaf of the plant) in kPa, C g denotes a coefficient that is equal to 610.7, and τ represents the surface temperature of the leaf, which is usually measured by means of an infrared gun for such experiments. For the sake of simplicity, this parameter was assumed to be equal to the temperature of the air inside the greenhouse at both ambient and evaporative assisted desiccant air-conditioning system supply air temperature. In fact, assuming leaf temperature equal to the surrounding air temperature is normal practice of greenhouse growers. The second part of Equation (7) represents the vapor pressure of the air inside the greenhouse in kPa, T represents the temperature of the air inside the greenhouse in • C, and RH represents the relative humidity of the air inside the greenhouse in (%). Because VPD is the function of temperature of the leaf surface, air temperature, and relative humidity, any slight variation in one of the parameters from the optimal values results in restricted growth of the plants, and can have significant consequences on reducing the overall yield and quality. The OptDeg model determines the variation of a single VPD measurement from the optimal references in an instant time. Although this model can provide an overview of the performance of the M-DAC system in sample time, a different model called comfort ratio, denoted by Cft-ratio, that incorporates a sample time frame, such as 24 h was used. The comfort ratio of VPD, denoted by C f t(VPD, t, α s ) GS = β, represents the percent of vapor pressure deficit at a specific time t and growth stage, which lies between the reference conditions of vapor pressure deficit related to α s . Ideal microclimate conditions inside a greenhouse could potentially be defined as C f t(VPD, t, 1) = 1. Self-adjusted optimality degree is denoted by α s to define the reference conditions required for performance evaluation of the proposed M-DAC system or control for microclimate evaluation. The set of Simulink blocks that were implemented to determine the Cft-ratio of the VPD values for different α s generated from the M-DAC system are shown in Figure 7. For this research, we selected three levels of α s as 0, 0.5, and 1, representing failure, marginal, and ideal production. It is worth mentioning that OptDeg and Cft-ratio models were derived with a prime focus on improved fruit quality and yield. The derivations at the back of the models (which can be found in published literature) for greenhouse tomato production were condensed for the sake of simplicity in this study.  Figure 8 shows the annual profile of the climatic ambient air conditions of Multan (Pakistan). The daily, hourly, and monthly variation of dry-bulb temperature shows that cooling inside the greenhouse environment is mostly required in summer months (i.e., May to August). However, the applicability of standalone pad-type direct evaporative cooling systems is limited during the monsoon (heavy rainfall season, i.e., July to August,  Figure 8 shows the annual profile of the climatic ambient air conditions of Multan (Pakistan). The daily, hourly, and monthly variation of dry-bulb temperature shows that cooling inside the greenhouse environment is mostly required in summer months (i.e., May to August). However, the applicability of standalone pad-type direct evaporative cooling systems is limited during the monsoon (heavy rainfall season, i.e., July to August, shown in Figure 8, humidity ratio (HR)) which results in excess relative humidity in the air. Consequently, a desiccant dehumidification-based Maisotsenko cycle evaporative cooling system was thermodynamically analyzed for the climatic conditions of Multan. Accordingly, results from the study are presented in Figures 8-13, which show the annual representation of the performance of the proposed system and the profile of the vapor pressure deficit inside the greenhouse. air. Consequently, a desiccant dehumidification-based Maisotsenko cycle evaporative cooling system was thermodynamically analyzed for the climatic conditions of Multan. Accordingly, results from the study are presented in Figures 8-13, which show the annual representation of the performance of the proposed system and the profile of the vapor pressure deficit inside the greenhouse. The standalone DAC system created a maximum temperature gradient, i.e., 8.4 °C more at ambient conditions of 38.6 °C with a humidity ratio of 21.6 g/kg at ambient conditions of 25 g/kg, which lies well outside the required greenhouse optimum temperature and humidity conditions (Figure 9a). It is worth mentioning that the temperature gradient is defined as the difference between the dry bulb temperature of the ambient air and the dry bulb temperature of the product air at the system outlet node. The studied standalone DAC system created unsuitable thermal and humidity conditions inside the greenhouse throughout the year due to a relatively higher temperature and very dehumidified air (Figure 9a). According to Figure 9b, the standalone DAC system coupled with a sensible heat exchanger (DAC+HX) (working efficiency assumed to be 0.9) created a maximum temperature gradient, i.e., 0.77 °C more at ambient conditions of 39.2 °C with a humidity ratio of 16.8 g/kg at ambient conditions of 24.3 g/kg, which lies slightly outside the required greenhouse optimum temperature and humidity conditions (Figure 9b), and is very close to ambient air conditions. The studied DAC+HX system created very close thermal and humidity conditions to ambient conditions inside the greenhouse throughout the year. Figure 9c shows the psychrometric profile of temperature and relative humidity of the desiccant dehumidification-based Maisotsenko cycle evaporative cooling (M-DAC) system for the climatic conditions of Multan (Pakistan). According to Figure 9c, most of the hourly points of the M-DAC system lie inside the required thermal and humidity conditions of the greenhouse environment. The M-DAC system created a maximum temper- The standalone DAC system created a maximum temperature gradient, i.e., 8.4 • C more at ambient conditions of 38.6 • C with a humidity ratio of 21.6 g/kg at ambient conditions of 25 g/kg, which lies well outside the required greenhouse optimum temperature and humidity conditions (Figure 9a). It is worth mentioning that the temperature gradient is defined as the difference between the dry bulb temperature of the ambient air and the dry bulb temperature of the product air at the system outlet node. The studied standalone DAC system created unsuitable thermal and humidity conditions inside the greenhouse throughout the year due to a relatively higher temperature and very dehumidified air (Figure 9a). According to Figure 9b, the standalone DAC system coupled with a sensible heat exchanger (DAC+HX) (working efficiency assumed to be 0.9) created a maximum temperature gradient, i.e., 0.77 • C more at ambient conditions of 39.2 • C with a humidity ratio of 16.8 g/kg at ambient conditions of 24.3 g/kg, which lies slightly outside the required greenhouse optimum temperature and humidity conditions (Figure 9b), and is very close to ambient air conditions. The studied DAC+HX system created very close thermal and humidity conditions to ambient conditions inside the greenhouse throughout the year. Figure 9c shows the psychrometric profile of temperature and relative humidity of the desiccant dehumidification-based Maisotsenko cycle evaporative cooling (M-DAC) system for the climatic conditions of Multan (Pakistan). According to Figure 9c, most of the hourly points of the M-DAC system lie inside the required thermal and humidity conditions of the greenhouse environment. The M-DAC system created a maximum temperature gradient i.e., 21.98 • C. According to Figure 9c, the M-DAC system coupled with a sensible heat exchanger (working efficiency assumed to be 0.9) created a maximum temperature gradient i.e., 21.9 • C more at ambient conditions of 39.2 • C with a humidity ratio of 16.8 g/kg at ambient conditions of 19.6 g/kg, which lies well inside the required thermal and humidity conditions inside the greenhouse environment (Figure 9c).       Figure 11 shows the temperature profile of the studied systems against the ambient air conditions of Multan (Pakistan) for greenhouse air-conditioning application. The ambient air conditions lie outside the required temperature conditions (marked by the dotted red line, Figure 11) from April to September (i.e., typical summer months) of Multan. Airconditioning is required inside the greenhouse environment during the summer months. Only M-DAC system created a uniform temperature gradient which lies well inside the required temperature conditions (marked in dotted red line, Figure 11) of the greenhouse environment. Figure 12 shows the temperature gradient of the DAC, DAC+HX, and M-DAC system in the summer months (i.e., May to August) for the greenhouse air-conditioning for the climatic conditions of Multan (Pakistan). The DAC system created higher than ambient temperature conditions throughout the summer months, which is unsuitable for greenhouse air-conditioning requirements ( Figure 12). The DAC + HX system created slightly higher/very close to ambient temperature conditions throughout the summer months, which is also not suitable for greenhouse air-conditioning requirements ( Figure  12). The M-DAC system created a temperature gradient which is highly suitable for the required temperature conditions throughout the summer months for greenhouse air-conditioning application (Figure 12).  Figure 11 shows the temperature profile of the studied systems against the ambient air conditions of Multan (Pakistan) for greenhouse air-conditioning application. The ambient air conditions lie outside the required temperature conditions (marked by the dotted red line, Figure 11) from April to September (i.e., typical summer months) of Multan. Air-conditioning is required inside the greenhouse environment during the summer months. Only M-DAC system created a uniform temperature gradient which lies well inside the required temperature conditions (marked in dotted red line, Figure 11) of the greenhouse environment. Figure 12 shows the temperature gradient of the DAC, DAC+HX, and M-DAC system in the summer months (i.e., May to August) for the greenhouse airconditioning for the climatic conditions of Multan (Pakistan). The DAC system created higher than ambient temperature conditions throughout the summer months, which is unsuitable for greenhouse air-conditioning requirements ( Figure 12). The DAC + HX system created slightly higher/very close to ambient temperature conditions throughout the summer months, which is also not suitable for greenhouse air-conditioning requirements ( Figure 12). The M-DAC system created a temperature gradient which is highly suitable for the required temperature conditions throughout the summer months for greenhouse air-conditioning application ( Figure 12). Figure 13 shows the dehumidification profile of the desiccant dehumidification unit of the proposed system for the summer months (i.e., May to August) for the climatic conditions of Multan (Pakistan). The desiccant dehumidification unit created a maximum dehumidification gradient of the ambient air, i.e., 6.61 g/kg at ambient conditions of 14.14 g/kg. It is worth mentioning that the dehumidification gradient is defined as the difference between the ambient air humidity ratio and the humidity ratio of the product air exiting the outlet of the proposed system. Compared to current proposed passive cooling techniques for greenhouse air-conditioning (i.e., pad-and-fan type systems for greenhouse air-conditioning systems), the proposed M-DAC system has the ability to manipulate the moisture inside the greenhouse environment, which is a key factor in defining the vapor pressure deficit inside the greenhouse. In turn, this defines the growth rate and ultimately results in better yield and fruit quality. Moreover, the M-DAC system can easily achieve the required temperature conditions inside a greenhouse compared to passive cooling techniques, i.e., pad-and-fan type AC systems.   Figure 13 shows the dehumidification profile of the desiccant dehumidification unit of the proposed system for the summer months (i.e., May to August) for the climatic conditions of Multan (Pakistan). The desiccant dehumidification unit created a maximum dehumidification gradient of the ambient air, i.e., 6.61 g/kg at ambient conditions of 14.14 g/kg. It is worth mentioning that the dehumidification gradient is defined as the difference between the ambient air humidity ratio and the humidity ratio of the product air exiting the outlet of the proposed system. Compared to current proposed passive cooling techniques for greenhouse air-conditioning (i.e., pad-and-fan type systems for greenhouse airconditioning systems), the proposed M-DAC system has the ability to manipulate the moisture inside the greenhouse environment, which is a key factor in defining the vapor pressure deficit inside the greenhouse. In turn, this defines the growth rate and ultimately results in better yield and fruit quality. Moreover, the M-DAC system can easily achieve   Figure 13 shows the dehumidification profile of the desiccant dehumidification unit of the proposed system for the summer months (i.e., May to August) for the climatic conditions of Multan (Pakistan). The desiccant dehumidification unit created a maximum dehumidification gradient of the ambient air, i.e., 6.61 g/kg at ambient conditions of 14.14 g/kg. It is worth mentioning that the dehumidification gradient is defined as the difference between the ambient air humidity ratio and the humidity ratio of the product air exiting the outlet of the proposed system. Compared to current proposed passive cooling techniques for greenhouse air-conditioning (i.e., pad-and-fan type systems for greenhouse airconditioning systems), the proposed M-DAC system has the ability to manipulate the moisture inside the greenhouse environment, which is a key factor in defining the vapor pressure deficit inside the greenhouse. In turn, this defines the growth rate and ultimately results in better yield and fruit quality. Moreover, the M-DAC system can easily achieve

Validation Based on Optimality Degree and Comfort Ratio
Validation of the performance of the proposed M-DC system with respect to the required optimal influential parameters of the greenhouse microclimate (i.e., temperature of the air, inside RH, and vapor pressure deficit) is shown in the plots of Figure 14 for 24 h. Results clearly show that the output of the M-DAC system was very close to the requirement microclimate data generated by the OptDeg model for the vegetative to mature fruiting growth stages. Moreover, the VPD plot implies that during the mid-day hours (i.e., 13:00-17:00) when the air temperature was at the highest peak and relative humidity had the lowest values, the VPD increased to 2 kPa, which is consistent with the output of the M-DAC system. Results of the VPD data generated from the OptDeg model for actual greenhouse tomato cultivation in 130 days versus the simulated VPD data generated from the M-DAC system are plotted in Figure 15. The first derivative of the VPD data sets, plotted in Figure 15a, shows that the optimum VPD values from the OptDeg model and the VPD values from the M-DAC system are not significantly different. This finding was also verified using the one-way analysis of variance (ANOVA) test for each of the five growth stages. Figure 15b shows the deviation of the VPD data of the M-DAC system from 100% optimal condition (α = 1), and compares that with actual VPD data of a commercial greenhouse that operated with ventilation and a pad-and-fan evaporative cooling system.

Validation Based on Optimality Degree and Comfort Ratio
Validation of the performance of the proposed M-DC system with respect to the required optimal influential parameters of the greenhouse microclimate (i.e., temperature of the air, inside RH, and vapor pressure deficit) is shown in the plots of Figure 14 for 24 h. Results clearly show that the output of the M-DAC system was very close to the requirement microclimate data generated by the OptDeg model for the vegetative to mature fruiting growth stages. Moreover, the VPD plot implies that during the mid-day hours (i.e., 13:00-17:00) when the air temperature was at the highest peak and relative humidity had the lowest values, the VPD increased to 2 kPa, which is consistent with the output of the M-DAC system. Results of the VPD data generated from the OptDeg model for actual greenhouse tomato cultivation in 130 days versus the simulated VPD data generated from the M-DAC system are plotted in Figure 15. The first derivative of the VPD data sets, plotted in Figure 15a, shows that the optimum VPD values from the OptDeg model and the VPD values from the M-DAC system are not significantly different. This finding was also verified using the one-way analysis of variance (ANOVA) test for each of the five growth stages. Figure 15b shows the deviation of the VPD data of the M-DAC system from 100% optimal condition ( = 1), and compares that with actual VPD data of a commercial greenhouse that operated with ventilation and a pad-and-fan evaporative cooling system. A more in-depth result for comparing the performance of the M-DAC system and evaporative pad-and-fan system in a random day of cultivation under fruit formation growth stage is provided by means of the Cft-ratio model in the plots of Figure 16. In Figure 16, the reference conditions related to the optimal, marginal, and failure vapor pressure deficit (i.e., α s = 1, 0.5, and 0, respectively) are presented with a green solid line, blue solid line, and red solid line, respectively. The Cft-ratio Simulink model then determined the percentage of VPD data of each case, shown in Figure 16 that were inside each reference border. From the graphical presentation, it is evident that the VPD data from the M-DAC system resulted in a higher comfort ratio than those from the evaporative pad-and-fan system. In detail, VPD values resulting from the M-DAC system never crossed the marginal reference border, whereas for the conventional pad-and-fan case, during the hours of 13:30-14:30, VPD exceeded this border and increased to 3 kPa, which is considered a failure value for greenhouse production. Sample results of the comfort ratio model for comparing the performance of the evaporative system and M-DAC system in providing optimum vapor pressure deficit for greenhouse production are given in Table 3 for two different growth stages of the tomato plant. A more in-depth result for comparing the performance of the M-DAC system and evaporative pad-and-fan system in a random day of cultivation under fruit formation growth stage is provided by means of the Cft-ratio model in the plots of Figure 16. In Figure 16, the reference conditions related to the optimal, marginal, and failure vapor pressure deficit (i.e., = 1,0.5, and 0, respectively) are presented with a green solid line, A more in-depth result for comparing the performance of the M-DAC system and evaporative pad-and-fan system in a random day of cultivation under fruit formation growth stage is provided by means of the Cft-ratio model in the plots of Figure 16. In Figure 16, the reference conditions related to the optimal, marginal, and failure vapor pressure deficit (i.e., = 1,0.5, and 0, respectively) are presented with a green solid line, pad-and-fan system. In detail, VPD values resulting from the M-DAC system never crossed the marginal reference border, whereas for the conventional pad-and-fan case, during the hours of 13:30-14:30, VPD exceeded this border and increased to 3 kPa, which is considered a failure value for greenhouse production. Sample results of the comfort ratio model for comparing the performance of the evaporative system and M-DAC system in providing optimum vapor pressure deficit for greenhouse production are given in Table 3 for two different growth stages of the tomato plant. Figure 16. Comparison between comfort ratios of VPD for 24 h resulting from the proposed M-DAC system and VPD data from a conventional evaporative (pad-and-fan) cooled greenhouse.

Conclusions
The present study aimed to investigate the applicability of desiccant dehumidificationbased Maisotsenko cycle evaporative cooling (M-DAC) system for greenhouse air-conditioning application for the climatic conditions of Multan (Pakistan). The study area (i.e., Multan) lies in a warm desert climate of the Köppen-Geiger climatic classification. In this regard, air-conditioning is required in greenhouses to optimize the temperature and humidity conditions, which could potentially increase the production level. Therefore, the present study proposed a desiccant dehumidification-based Maisotsenko cycle evaporative cooling system (M-DAC) for greenhouse air-conditioning application. The standalone DAC system created a (maximum) dehumidification gradient (i.e., 16.8 g/kg) at 24.3 g/kg ambient air conditions, and a (maximum) temperature gradient (i.e., 8.4 • C) at 38.6 • C ambient air conditions. The DAC+HX system created a temperature gradient roughly equal to the ambient air conditions, which are unsuitable for greenhouse air-conditioning. The M-DAC system created a (maximum) temperature gradient (i.e., 21.9 • C) at 39.2 • C ambient air conditions, which lies well within the required optimum temperature conditions of the greenhouse air-conditioning. Additionally, the M-DAC created a dehumidification gradient (i.e., 16.8 g/kg) at 24 g/kg ambient conditions, which also lies inside the range of optimum humidity conditions inside the greenhouse. Moreover, the vapor pressure deficit (VPD) profile of the greenhouse environment related to the growth of plants indicated that a VPD of 0.8-1.2 kPa (for leaf temperature equal to surrounding air temperature) is most suitable for healthy transpiration from the plant leaves at late vegetation and early flowering stages. These results were validated with OptDeg and Cft-ratio microclimate models, which showed that the M-DAC system maintained the VPD of the product air inside the healthy transpiration zone throughout most of the summer months (i.e., May to August). Moreover, the comfort ratio model results of the proposed M-DAC system were compared with experimental results from a tomato greenhouse, indicating that the M-DAC system could potentially achieve the required comfort level compared to a pad-and-fan type air-conditioning system. Thus, the study indicates the M-DAC system is a feasible option for greenhouse air-conditioning for the climatic conditions of Multan (Pakistan).