Solar Technology and District Cooling System in a Hot Climate Regions: Optimal Conﬁguration and Technology Selection

: With the increasing need for cooling and the concerns for pollution due to fossil fuel-based energy use, renewable energy is considered an add-on to cooling technologies. The climatic condition in the Middle East, analyzed in this paper, provides the potential to integrate solar energy with the cooling system. Due to the availability of various solar energy and cooling technologies, multiple conﬁgurations of solar-cooling systems can be considered to satisfy the cooling demand. The research presented in this paper aims to assess and compare these conﬁgurations by considering the energy prices and the installation area. The proposed model is formulated in Mixed-Integer Linear Programming and optimizes the holistic system design and operation. The economic, renewable energy use, and environmental performances of the optimal solution for each conﬁguration are analyzed and compared to the base grid-DCS conﬁguration. Results show that the electricity tariff and the available installation area impact the economic competitiveness of the solar energy integration. When electricity tariff is subsided (low), the conventional grid-based DCS is the most competitive. The PV-DCS conﬁguration is economically competitive among the solar assisted cooling systems, and it can contribute to reducing the environmental impact by 58.3%. The PVT-DCS conﬁguration has the lowest operation cost and the highest environmental performance by decreasing the global warming potential by 89.5%. The T-DCS conﬁguration becomes economically competitive only at high electricity tariffs.


Introduction
In 2016, the total electricity consumption for the cooling required nearly 10% of the electricity consumed globally.From an environmental perspective, the cooling-related CO 2 emissions have tripled since 1990 and reached 1135 MT in 2016; the share of cooling to the total emissions from buildings represents 12% over the same period [1].In the case of Qatar, the consumption for annual cooling demand is more than 14TWh, and it represents about 36% of the annual electricity consumption [2].Out of this consumption, more than 50% is required for cooling during the summer season [3].The district cooling system (DCS) is a centralized plant for the production and circulation of chilled water through a network of insulated pipes laid underground [4].The system comprises a central plant, including the heat rejection system, pumps, chillers, distribution network, and customer interconnection.These systems provide a cooling load to a range of buildings (residential, commercial, industrial, and office).The investigation of Gang et al. [5] concerning the performance of district cooling systems showed that DCS could save 13% of electricity consumption, and its operation cost is lower by 10% than individual cooling at the building level.In Qatar, the Planning and Statistics Authority (PSA) shows that up to 2018, 41 district cooling plants The proposed integrated solar-district cooling system is connected to the electricity grid [12], which allows the exchange of electricity during peak and low solar irradiation periods.The synergy between the generated and consumed energy, respectively, by the solar collector (thermal, electric, or both) and the cooling process (absorption and vapor compression) and the grid connection is considered coherent.The cooling demand is required and varies throughout the day.Similarly, the solar irradiation intensity is not constant.Consequently, the design, components inter-connections, and operational policies should satisfy the cooling demand by adjusting the cooling production.Therefore, there is a need for a systemic optimization approach that leads to an optimal system design and operation.This will allow finding the optimal components selection, including the cooling technologies and operational policies (quantity to be produced and stored during each period) while minimizing the overall costs.The proposed mixed-integer linear programming (MILP) model captures all these aspects and seeks to optimize each configuration design without requiring multi-models to imitate each configuration.
Further discussion in this paper is divided into five sections.In Section 2, the literature review on solar cooling technology is given.Section 3 introduces the problem description and formulation.The experiments and results are given in Section 4. The conclusions, limitations, and future research are given in Section 5.

Literature Review
Research on solar-assisted cooling applications has been initiated for more than four decades (after the energy crisis in the 1970s).During the last decade, the number of research on solar energy for cooling applications has doubled due to the increasing need for cooling and the availability of improved solar cooling technologies.The bibliometric analysis on solar cooling technology shows that the most studied systems are solar thermal cooling systems using the thermal collector and thermally activated chiller (mainly absorption), solar combined power/cooling, and photovoltaic (PV) compression chillers [13].Hybrid technology (photovoltaic-thermal) is also the subject of many types of research [14][15][16][17][18] and its applications [19][20][21].Most current applications are focused on building-level operations, for example, water heating [22], heating and cooling [23], drying [24], and other applications such as desalination [25] and industrial [26].The solar heating and cooling program (IEA SHC task 60) [27][28][29][30][31][32] of the International Energy Agency (IEA) also recommends further investigation concerning the use of solar hybrid technology in cooling applications, the interaction of local electricity use versus grid injection.Therefore, there is an increasing interest in solar hybrid technology in the cooling application, but the applications are lacking for district cooling.
The main approach used in studying solar cooling technology is a simulation using the predefined model components in the tool library.The optimization focuses on the process parameter to improve system efficiency for a given design.Abid et al. [33] analyzed and compared the thermodynamics behavior of different absorption cycles (single, double, triple, and quadruple effect), evacuated solar tube collectors, and hybrid nanofluids.They concluded that the collector efficiency is higher when using hybrid nanofluids.Asadi et al. [34] used the thermo-economic analysis to evaluate an integrated 10 kW singleeffect ammonia-water absorption chiller with diverse solar collectors.They concluded that integrating the evacuated tube collectors with the chillers can provide an economical solution.Weber et al. [35] studied integrated water ammonia absorption chillers and concentrated solar collectors and the performance of the different control and operation strategies.Bellos and Tzivanidis [36] investigated the performance of a solar cooling system that comprises an absorption chiller and nanofluids using plate collectors; they found that nanofluids improve the thermal efficiency by 2.5%, the exegetic performance by 4%, and increase in refrigeration by 0.84%.The authors have also studied the parametric analysis and the optimization of a combined ejector-absorption chiller and parabolic collector system and concluded that the ejector increases the system COP, and the optimized system can achieve a performance enhancement of more than 60% for the respective evaporating and condenser temperature of 12.5 • C and 30 • C [37].Shehadi [38] simulated different absorption chiller components for a range temperature to optimize the chiller COP.The author finds that the optimum COP of the system is 0.776 when the evaporator, condenser (and absorber), and generator temperatures, respectively, are equal to 10 • C, 30 • C, and more than 70 • C. Saleh and Mosa [39] optimize the solar absorption chiller in a hot climate by using a flat plate solar collector with a single-effect water-lithium bromide absorption chiller.Results show that a coating of solar collector is necessary to produce the required operating temperature (70-80 • C) of the chiller.The authors show that their system improved the COP beyond 0.8.Iranmanesh and Mehrabian [40] used a genetic algorithm to optimize an integrated double-effect bromide water absorption chiller and evacuated tube collector system.They concluded that the reduction of auxiliary energy could be achieved by optimizing mass flow rates of hot water passing through the generator and collector.Alobaid et al. [41] reviewed the solar-driven absorption cooling with photovoltaic thermal systems.The authors mentioned that the solar absorption chiller saved 50% of primary energy with electrical PVT efficiency in the range 10-35%, the solar-assisted cooling COP in the range 0.1-0.91, the thermal collector efficiency in the range 0.06-0.64.Alghool et al. [42] developed a MILP model to optimize integrated solar-assisted DCS, the optimal solution that provides the design and operational variables of the system.
The literature review shows that the optimization approach is mainly used for the process parameters rather than the holistic system design and operation.The cost is commonly used when optimizing solar cooling technologies, and the hybrid technology application in district cooling is lacking.Additionally, a model that supports decision-makers in assessing and selecting the appropriate configuration among different alternatives (Table 1) in the context of local energy prices and the area available for installation is also lacking.
Understanding modeling, optimizing, and assessing these configurations becomes necessary for decision-makers.Hence, the proposed MILP model optimizes each configuration design, including the component sizing and related technology, the components interconnection, and the operational policies to satisfy the cooling demand during each period.Furthermore, the performance analysis used in this paper is based on sustainable district cooling guidelines (IEA 2019) [43] by analyzing the environmental benefits of alternative energy production systems, which reduces the emission of greenhouse gasses (GHG) and increase the use of renewable energy with continuity and reliability of energy supply.Hence a set of performance indicators are used to assess the system from the economic, renewable energy use, and environmental perspectives as defined in [33].
The research in this paper aims to: Develop a generalized model to optimize the design and operation of each configuration, and thus without the need for multi-models to imitate each one.Develop and implement the performance indicators of the system based on the numerical results of the proposed model.Compare the diverse configurations performances (economic, renewable energy use, and environmental) by considering the local context (energy price and available installation area).

Problem Description and Formulation
It is worth noting that The TRNSYS is one of the tools that can be used to simulate the configurations mentioned in this paper once the design is fixed.It can help consolidate the design through specific analysis.However, modeling becomes necessary to optimize the design of each configuration and select the appropriate components, including the technology.Such a step is necessary to compare further the different configurations based on the economic and environmental aspects and perform other analyses such as the effect of the electricity tariff and the reserved installation area.
Figure 1 shows the proposed integrated system diagram.It is composed mainly of (i) solar collector representing a given technology (photovoltaic, thermal, and hybrid), (ii) vapor compression chiller, (iii) absorption chiller, (iv) chilled water storage tank, (v) hot water storage tank, (vi) auxiliary boiler.The system is connected to the electricity grid.It is worth noting that the inclusion of hot water and chilled water thermal energy storage and the auxiliary boiler is not mandatory in the system.However, their inclusion contributes to lowering the total system cost.

Problem Formulation
The proposed model assumes that the cooling load is deterministic, the solar collector efficiency is constant, the system operates in a steady-state, and the cool-losses are negligible.The efficiencies of the auxiliary boiler, vapor compression, and absorption chiller are assumed as constant.It is worth noting that chiller performance depends on the pa- Figure 1 also shows the energy flow, hot water, chilled water, gas, and electricity between the different system components.Based on the used collector technology (photovoltaic, thermal, hybrid), the collected solar irradiance is transformed into electricity, hot water, or both.The generated electricity drives the vapor compression chiller.In the case of extra-generation during the peak sunlight, the surplus electricity is sold back to the grid.If the generated electricity by the solar collector is insufficient, the grid feeds the system with the needed electricity.The hot water generated feeds the absorption chiller.The surplus heat generated during the peak irradiance will be stored in a hot water storage tank to be used by the absorption chiller during limited irradiance periods.Suppose the generated hot water by the solar collector plus the inventory level of the hot water storage tank is insufficient to run the absorption chiller, or the hot water temperature is under the operating temperature of the absorption chiller (80 • C).In that case, the auxiliary boiler is activated to provide the requested amount of heat.The chilled water produced by the vapor-compression and absorption chillers feeds the building to satisfy the cooling demand.Extra chilled water is produced and stored in the chilled water storage tank during peak irradiance.
The integration of these components forms an interconnected system with intertwined behaviors.The proposed model captures all these interdependencies and provides an optimal solution based on an approximated annual hourly cooling demand and solar irradiance (8760 h).

Problem Formulation
The proposed model assumes that the cooling load is deterministic, the solar collector efficiency is constant, the system operates in a steady-state, and the cool-losses are negligible.The efficiencies of the auxiliary boiler, vapor compression, and absorption chiller are assumed as constant.It is worth noting that chiller performance depends on the parameters such as the partial load ratio and the weather conditions.Integrating these aspects in the model requires additional decision variables and constraints.Such an addition increases problem complexity in terms of total numbers of binary variables, non-integers variables, and constraints, making it unlikely that the real life-size problem of the solar energy-DCS integrated system could be optimally solved.Therefore, the developed model is based on a trade-off between the model breadth (number of components, operational decisions, and time horizon length) and model depth (degree of fidelity of component model in terms of efficiency) that aims at providing a solution in a reasonable computational time.Hence, it can support decision-makers to ensure the feasibility of such an integrated system and select the appropriate configuration in the context of energy prices and available installation areas for the solar collector.The detailed proposed model nomenclature is given below.
Figure 2 shows the organization of the proposed model is similar to [44].The model is fed with solar irradiance, cooling demand, energy tariff, and component attributes (capacity, cost, efficiency).The core of the optimization process is composed of the structure and the behavior.The structure refers to the design, including the definition block (components) and internal block reflecting the interconnection between the components (Figure 1).The behavior (operation) refers to the system activities related to energy (E, G, H, C) production, consumption, and storage (CI, HI).The model's output is the optimal design and operation that minimizes the objective function, which represents the configuration's economic performance.The renewable energy use and environmental performance are obtained after the post-treatment of the numerical value of the optimal solution for each configuration.
The optimal operation seeks to find the quantity of chilled water (C), hot water (H), and electricity (E) to be produced, consumed, and stored (CI, HI) for each period to satisfy the cooling demand.It is worth noting that every component has attributes (constant parameters in the proposed model) such as capacity, efficiency, and costs (acquisition and maintenance).Some of these parameters are used in the model constraints to define the feasible region.Others (cost parameters) are used in the objective function to evaluate the solution's fitness.The objective function (Equation ( 1)) minimizes the total annualized cost and includes three terms: the first term is for the capital cost for all implemented equipment; the second term is the operating cost of the running equipment (boiler and vapor compression chiller) that is required energy from outside the system, and the third term is for the maintenance cost.
The investment cost (Equation ( 2)) includes all installed equipment (solar collector, vapor compression chiller, absorption chiller, hot water storage tank, chilled water storage tank, and auxiliary boiler) [45].The annualized investment cost is obtained by summing the product of each selected equipment's fixed cost () by its amount to be installed and the capital recovery factor (CRF = ).The area to be installed for the solar collector (thermal, photovoltaic, or hybrid) is represented by the variable (x).For the remaining, each piece of equipment is represented by its design binary variable ( ,  , , , ).Each technology is represented by a set of elements and identified by its index in the summation (Equation (2)).For example, index (i) represents various vapor compression chillers candidates, and each one is characterized by its attributes such as capacity, acquisition and maintenance costs, and efficiency.The optimal design seeks to find the required components and their interactions.Each component is presented by a design variable (binary or non-integer).For example, if the vapor compression chiller (i) is selected, the related design variable ( y i ) takes the value one and zero otherwise, which is similar for the absorption chiller (z j ), hot water tank (u k ), chiller water tank (v m ), and the auxiliary boiler (w n ).However, the design variable (x) represents the area to be installed for solar collector technology.If the variable (x > 0), the considered configuration (1 or 2 or 3) is economically competitive than the base configuration (grid-DCS).Otherwise, when (x = 0), the base configuration is optimal.
The optimal operation seeks to find the quantity of chilled water (C), hot water (H), and electricity (E) to be produced, consumed, and stored (CI, HI) for each period to satisfy the cooling demand.It is worth noting that every component has attributes (constant parameters in the proposed model) such as capacity, efficiency, and costs (acquisition and maintenance).Some of these parameters are used in the model constraints to define the feasible region.Others (cost parameters) are used in the objective function to evaluate the solution's fitness.
The objective function (Equation ( 1)) minimizes the total annualized cost and includes three terms: the first term is for the capital cost for all implemented equipment; the second term is the operating cost of the running equipment (boiler and vapor compression chiller) that is required energy from outside the system, and the third term is for the maintenance cost. Minimize The investment cost (Equation ( 2)) includes all installed equipment (solar collector, vapor compression chiller, absorption chiller, hot water storage tank, chilled water storage tank, and auxiliary boiler) [45].The annualized investment cost is obtained by summing the product of each selected equipment's fixed cost (Fc) by its amount to be installed and the capital recovery factor (CRF = r(1+r) n (1+r) n −1 ).The area to be installed for the solar collector (thermal, photovoltaic, or hybrid) is represented by the variable (x).For the remaining, each piece of equipment is represented by its design binary variable (y, z, u, v, w).Each technology is represented by a set of elements and identified by its index in the summation (Equation (2)).For example, index (i) represents various vapor compression chillers can-Energies 2022, 15, 2657 7 of 24 didates, and each one is characterized by its attributes such as capacity, acquisition and maintenance costs, and efficiency.
The operating cost (Equation ( 3)) represents the yearly (8760 h) sum of the running cost for all implemented equipment in the district cooling and requiring energy from outside the system, which are the boiler and the vapor compression chiller [46].The first term in Equation ( 3) is the cost for running the auxiliary boiler.It is equal to the product of the gas price (P Gas ) and the amount of gas consumed by the boiler (G Bo ).The second term in Equation ( 3) represents the cost related to the net electricity of the system.The net electricity cost is equal to the product of the electricity cost (P Elec ) and the difference between the amount of electricity purchased (E Grid→Vc ) and sold to the grid (E Col→Grid ).If the purchase and selling costs are different, the feed-in tariff coefficient (α) is used as a correction factor.If the solar panel electricity yield is more than the need of the electric chiller, the extra generated electricity is sold to the grid and will be counted as an income.Consequently, the overall system cost will be decreased.Contrary, if the solar collector-generated electricity is insufficient, the system will purchase electricity from the grid, and the overall system cost increases.
The annual maintenance cost equals the sum of the annual cost of maintaining all the system components.The maintenance cost of each element is equal to the product of its maintenance cost (Fm) and the installed amount.The installed amount for the vapor compression, absorption chiller, hot and chilled water tanks, and the auxiliary boiler is equal to its design binary variable (y, z, u, v, w).For the solar panel is equal to the installed area (x).
The constraints for the model are classified into three blocks.The definition block describes the system in terms of structure (components).The internal block describes the system layout or components interconnection and translates the energy balance.The activities block describes the system's elements behavior in terms of consumed, generated, and stored energy.

Definition Block
Constraints ( 5)-(9) define the system structure in terms of components.Six components are integrated into the proposed system: solar collector, vapor compression chiller, absorption chiller, hot water storage tank, chilled water storage tank, and auxiliary boiler.
Constraint (5) ensures that the total installed area of the solar collector does not exceed the available area (A).The available area for installing solar collectors depends on the building type.In the urban region, government institutions like universities and hospitals cover a great area and generally have more available space than commercial and residential buildings.
Energies 2022, 15, 2657 8 of 24 Constraint (6) ensures that two technologies of the chilling process (vapor compression and absorption) are a candidate to produce the chilled water to satisfy the cooling demand, and at least one technology should be selected.Therefore, the system structure may contain an absorption chiller, electric chiller, or both.
Constraints ( 7)-( 9) ensure that selecting the hot water storage tank, chilled water storage tank, and the auxiliary boiler is optional.Their selection contributes to decreasing the total cost.

Internal Block
The constraints ( 10)-( 18) describe the interconnection between the system components (the system layout).Hence, they describe the energy balance (arrow in Figure 1).
Constraint (10) imposes the solar collector (Col) thermal energy balance (H).It shows that the hot water produced (H) by the solar collector (col) is delivered to the hot water storage tank (Col → Hwt ) and the absorption chiller ( Col → Ab ).Constraint (11) imposes the electric energy balance (E) for the solar collector (Col), ensuring that the produced electricity (E) equals the amount delivered to the grid ( Col → Grid ) and vapor compression chiller ( Col → Vc ).
Constraint (12) imposes the energy balance for the hot water storage tank (Hwt).It imposes that the sum of hot water received from the solar collector ( Col → Hwt ) and the inventory level (H I) at the period (t − 1) equals the sum delivered to the absorption chiller ( Hwt → Ab ) and the inventory level (H I) at the period t.
Like constraints ( 12) and ( 13) imposes the energy balance for the chilled water storage tank (Cwt).It imposes that the sum of chilled water (C) received from the absorption ( Ab → Cwt ) and vapor compression chillers ( Vc → Cwt ) plus the inventory level (CI) at the period (t − 1) equals the amount delivered to the building ( Cwt → Bld ) plus the inventory level (CI) at the period t.
Constraints ( 14) and ( 15) represent the absorption chiller's (Ab) hot water (H) and chilled water (C) balances.The left-hand side of constraint (14) reflects the needed hot water to operate the absorption chiller as the sum of hot water provided by the solar collector ( Col → Ab ), the hot water storage tank ( Hwt → Ab ), and the auxiliary boiler ( Bo → Ab ).However, constraint (15) defines that the produced chilled water by the absorption chiller equals the quantity delivered to the building ( b → Bld ) and the chilled water storage tank ( Ab → Cwt ).
Constraints ( 16) and ( 17) impose the vapor compression chiller's (Vc) electricity (E) and chilled water (C) balances.The left-hand side of constraint ( 16) reflects the needed electricity to operate the vapor compression chiller, and the right-hand side represents the electricity delivered from the solar collector ( Col → Vc ) and the grid ( Grid → Vc ).Constraint (17) represents the balance of chilled water (C); the chilled water produced equals the amount delivered to the building ( Vc → Bld ) and the chilled water storage tank ( Vc → Cwt ).
Constraint (18) imposes the thermal energy balance (C) for the building cooling demand satisfaction (d).It enforces that the cooling demand is assumed by the amount of chilled water delivered from the vapor compression chiller ( Vc → Bld ), the absorption chiller ( Ab → Bld ), and the chilled water storage tank ( Cwt → Bld ).

Activities Block
The constraints ( 19)-( 28) describe the behavior of each system-components in terms of activities, which are energy generation (electricity (E), cooling (C), heat (H)), consumption (electricity and heat), and storage (chilled (CI) and hot water (HI)).
Constraints (19) and (20) ensure that the amount of electricity (E) and the amount of hot water (H) produced by the solar collector (Col) do not exceed the capacity of the installed area (x).The solar irradiance, electric and thermal efficiency of the solar collector are respectively S, η elec , and η ther .
Constraints ( 21)- (22) ensure that the inventory level of hot water (H I) and chilled water (CI) does not exceed respectively, the installed capacity (Q) of the hot water tank (Hwt) and chilled water tank (Cwt).
Constraint (23) ensures that heat production (H) of the auxiliary boiler (Bo) does not exceed its installed capacity (Q).Constraint (24) ensures that the relationship between the heat generated (H) and the required amount of gas (G) in the boiler can be expressed as a linear expression using the efficiency (η Bo ).The efficiency is considered constant during heat production.
Energies 2022, 15, 2657 10 of 24 Constraint (25) ensures that the amount of chilled water (C) produced by the installed absorption chiller (Ab) does not exceed the chiller capacity (Q).Constraint (26) expresses the linear relationship between the cooling energy and the required heat in the absorption chiller by using the energy efficiency ratio (EER) which is considered as a constant.
Similar to constraints (25)- (27) ensures that the amount of chilled water (C) produced by the installed vapor compression chiller (Vc) does not exceed the chiller capacity (Q).The energy efficiency ratio (EER) of the vapor compression chiller defines the linear relationship between the cooling energy (C) produced and the required electricity (E). Linearization: The right-hand side of constraints ( 24), ( 26) and ( 28) contain a product of binary and non-integers variables w n G Bo (n,t) , z j H Ab (j,t) , and y i E Vc (i,t) , which can be linearized by performing a variable change and using the BigM method [42].

Assessment Indicators
The performance indicators for integrated solar systems and their definition are introduced in SHC Task60\Report D1 [24].They are classified into three sub-categories: economic, energy, and environmental.These indicators allow the assessment of the system over a certain period.The numerical results obtained from the proposed model allow the calculation of these indicators and the assessment of configurations over one running year (8760 h).

Economic indicators:
The model objective function represents the economic indicator and includes the annualized capital, operational, and maintenance costs.

Energy indicators:
The energy indicators assess the system based on renewable energy use and saving.
-The solar electrical fraction is expressed as the ratio of the yield of electricity and electric chiller consumption: The solar thermal fraction is expressed as the yield heat and the absorption chiller The renewable energy fraction is expressed as the ratio of the sum of the heat and electricity yield and the sum of absorption and electric chillers consumption: .
-The system's final energy saving for gas is expressed as the ratio of the heat yield and the efficiency of the selected boiler.If the optimal system design does not include boiler, an efficiency reference of 0.85 is considered: The system's final energy-saving for electricity equals the electricity yield:

Environment indicators:
The environmental indicators assess the system based on the non-renewable use and its impact on global warming.

-
The non-renewable primary energy consumption (kWh oil-equivalents) is expressed as the sum of the primary energy of the gas and net electricity consumption.The primary energy consumption for gas is equal to the product of the quantity consumed and its cumulated energy demand (CED = 1.06).The primary energy consumption for electricity is equal to the product of the net quantity consumed and the cumulated energy ).

Numerical Results
The proposed MILP model is coded and implemented using AIMMS 4.82.6.10 software, solved with CPLEX 20.1 branch and bound algorithm on Intel(R) Core(TM) i7-1165G7 CPU@2.80GHz with 48 GB RAM.The MILP model is generalized to represent four system configurations: each configuration represents an integrated energy source for DCS.The model parameters setting of the solar collector efficiency η ther and η elec allow the selection of the desired system.Therefore, the proposed model imitates integrating all solar technology with a district cooling system in addition to the grid-powered DCS without the need for multi-models.Table 2 illustrates the setting parameters for each configuration.For example, to optimize configuration (2), which represents a photovoltaic integrated district cooling system, the thermal efficiency parameter of the model should be set to zero.Similarly, to optimize the base configuration (Grid-DCS), the thermal and electric efficiency should be set to zero.The optimal selection of the cooling technology is performed (i) via the presence of constraint (5), where both technologies (absorption and electric chillers) are a candidate for the system design without forcing a specific one, (ii) the objective function, which includes the investment, maintenance, and operation costs.This latter considers the energy costs (electricity and gas) for running the system.It is worth noting that the integrated system is connected to the grid, and both cooling technologies are the candidates.Therefore, the optimal solution of configurations 1, 2, and 3 may lead to configuration 4, which refers to us as the base system (Grid-DCS), which means that the integration of the solar technology is not competitive compared to the Grid-DCS.
Figure 3 summarizes the structure of the data used in the model.The solar collector cost and efficiency are collected from IEA. SHC TASK 60 report [47] and solar collector providers.The remaining system components (absorption chiller, vapor compression chiller, hot water, chiller water TES, and auxiliary boiler) data such as capacity, fixed cost, coefficient of performance, and the cooling load estimation are taken from the data paper [48].The electricity cost is collected from the energy provider KAHRAMAA [49].The electricity tariff in Qatar is subsidized, and there is a plan to remove these subsidies by employing the time of use or dynamic pricing schemes [3].Currently, the electricity price (P Elec ) depends on the customer category (productive farms, residential, industrial, hotel, commercial, and government).The average tariff is $0.055 /kWh, and the feed-in-tariff coefficient (α) is set to one.The gas cost fluctuates over time; therefore, an average tariff of $0.017 kWh is considered [50].The life span (n) considered of the study system is 25 years, and the interest rate (r) is 6%.The maximum allowed area (A) for the solar collector installation is set to 6000 m 2 .The solar irradiance represents the typical meteorological year 2016 (TMY) of Doha/Qatar.The peak cooling load considered is 10,000 kW.The objective function is obtained by running the system for one year, equivalent to an 8760 hourly period.
Table 3 details the solar collector used parameters values.The detailed characteristics of chillers, storage tanks, and auxiliary data can be found in [42].The set size column in Table 3 represents the number of elements for each set.
In research, hybrid solar technology (PVT) is less investigated in district cooling applications than the others.Therefore, configuration (1) detailed optimal solution is presented for illustration in the following section.

Integration of Hybrid Solar Collector and DCS
Figure 4 shows the cooling load and the solar irradiance patterns for a typical summer day.The cooling demand is over the whole day with an important load from 11 AM to 10 PM.The solar irradiance is present from 6 AM to 6 PM, with a peak at 1 PM.The storage becomes essential to profit from extra energy (hot water and electricity) during the peak irradiance.The connection to the grid allows for selling the surplus generated electricity and purchasing when needed.
The optimal design of the integrated system is analyzed in terms of design, cost, electricity flow, heat flow, and cooling flow.The detailed energy and environmental characteristics of the system are calculated based on the numerical results of the model.The implementation of these indicators allowed us to further compare the different configurations from an energy and environmental perspective.
Table 4 shows the optimal design for the configuration (1).The optimal installed area is equal to the maximum allowed (6000 m 2 ).Both cooling technologies (vapor compression and absorption) are selected, with capacities of 5300 and 1454 kW.Two storage tanks for hot and chilled water are selected with the same capacity of 63,000 kW.However, the auxiliary boiler is not selected in the optimal design.The investment and maintenance cost of configuration (1) is higher by 416% compared to the base configuration.However, Table 3 details the solar collector used parameters values.The detailed characteristics of chillers, storage tanks, and auxiliary data can be found in [42].The set size column in Table 3 represents the number of elements for each set.In research, hybrid solar technology (PVT) is less investigated in district cooling applications than the others.Therefore, configuration (1) detailed optimal solution is presented for illustration in the following section.

Integration of Hybrid Solar Collector and DCS
Figure 4 shows the cooling load and the solar irradiance patterns for a typical summer day.The cooling demand is over the whole day with an important load from 11 AM to 10 PM.The solar irradiance is present from 6 AM to 6 PM, with a peak at 1 PM.The storage becomes essential to profit from extra energy (hot water and electricity) during the peak irradiance.The connection to the grid allows for selling the surplus generated electricity and purchasing when needed.

Operation cost
321,286 129,363 Figure 5 shows the system electricity flow.The generated electricity feeds the electric chiller during the daytime, and the excess is transmitted to the grid.During the absence of solar irradiance, the grid supplies the chiller.The installed area covers 45.07% of the annual DCS electricity consumption.The annual generated electricity equals 1,930,012 kWh, which is considered energy saved.The optimal design of the integrated system is analyzed in terms of design, cost, electricity flow, heat flow, and cooling flow.The detailed energy and environmental characteristics of the system are calculated based on the numerical results of the model.The implementation of these indicators allowed us to further compare the different configurations from an energy and environmental perspective.
Table 4 shows the optimal design for the configuration (1).The optimal installed area is equal to the maximum allowed (6000 m 2 ).Both cooling technologies (vapor compression and absorption) are selected, with capacities of 5300 and 1454 kW.Two storage tanks for hot and chilled water are selected with the same capacity of 63,000 kW.However, the auxiliary boiler is not selected in the optimal design.The investment and maintenance cost of configuration (1) is higher by 416% compared to the base configuration.However, the total cost is lower by 5.25%.This is due to the low operation cost (129,363$) in comparison with the base configuration (321,286$), which is equivalent to decreasing the energy bill by more than 59%. Figure 5 shows the system electricity flow.The generated electricity feeds the electric chiller during the daytime, and the excess is transmitted to the grid.During the absence of solar irradiance, the grid supplies the chiller.The installed area covers 45.07% of the annual  Figure 6 shows the system heat flow.The generated solar heat feeds the absorption chiller during the peak solar irradiance.The excess heat is stored in the hot water tank; during the absence of solar irradiance, the stored hot water operates the absorption chiller.In contrast to electricity, the heat generated covers 100% of the absorption chiller need and saves the equivalent of 9,055,221 kWh of natural gas.Considering the electricity and the heat generated, renewable energy covers 80.36% of the total energy to operate the district cooling.The yearly primary energy consumption and the global warming potential are 6,797,412 kWh oil-equivalent and 1,232,472 kg CO2-equivalent, which contributes to decreasing the environmental impact of the Grid-DCS configuration by 59.74% (Table 5).
Figure 7 shows the system cooling flow.Absorption and electric chiller contribute to satisfying the cooling load demand.The absorption chiller is running all day hours with its total capacity.However, the electric chiller production is adjusted to supply the complement.The surplus production during the morning hours (5 a.m. to 9 a.m.) is stored in the chilled water tank and used during the peak hours (Figure 7).The absorption chiller ensures 21.15% of the annual cooling demand, and the electric chiller production supplies the remaining amount.Table 5 summarizes the performance of configuration 1 and the base configuration.Figure 6 shows the system heat flow.The generated solar heat feeds the absorption chiller during the peak solar irradiance.The excess heat is stored in the hot water tank; during the absence of solar irradiance, the stored hot water operates the absorption chiller.In contrast to electricity, the heat generated covers 100% of the absorption chiller need and saves the equivalent of 9,055,221 kWh of natural gas.Considering the electricity and the heat generated, renewable energy covers 80.36% of the total energy to operate the district cooling.The yearly primary energy consumption and the global warming potential are 6,797,412 kWh oil-equivalent and 1,232,472 kg CO 2 -equivalent, which contributes to decreasing the environmental impact of the Grid-DCS configuration by 59.74% (Table 5).Figure 6 shows the system heat flow.The generated solar heat feeds the absorption chiller during the peak solar irradiance.The excess heat is stored in the hot water tank during the absence of solar irradiance, the stored hot water operates the absorption chiller In contrast to electricity, the heat generated covers 100% of the absorption chiller need and saves the equivalent of 9,055,221 kWh of natural gas.Considering the electricity and th heat generated, renewable energy covers 80.36% of the total energy to operate the distric cooling.The yearly primary energy consumption and the global warming potential ar 6,797,412 kWh oil-equivalent and 1,232,472 kg CO2-equivalent, which contributes to de creasing the environmental impact of the Grid-DCS configuration by 59.74% (Table 5).
Figure 7 shows the system cooling flow.Absorption and electric chiller contribute to satisfying the cooling load demand.The absorption chiller is running all day hours with its total capacity.However, the electric chiller production is adjusted to supply the com plement.The surplus production during the morning hours (5 a.m. to 9 a.m.) is stored in the chilled water tank and used during the peak hours (Figure 7).The absorption chille ensures 21.15% of the annual cooling demand, and the electric chiller production supplie the remaining amount.Table 5 summarizes the performance of configuration 1 and th base configuration.Figure 7 shows the system cooling flow.Absorption and electric chiller contribute to satisfying the cooling load demand.The absorption chiller is running all day hours with its total capacity.However, the electric chiller production is adjusted to supply the complement.The surplus production during the morning hours (5 a.m. to 9 a.m.) is stored in the chilled water tank and used during the peak hours (Figure 7).The absorption chiller ensures 21.15% of the annual cooling demand, and the electric chiller production supplies the remaining amount.Table 5 summarizes the performance of configuration 1 and the base configuration.
Energies 2022, 15, x FOR PEER REVIEW

Scenarios and Configurations Assessment
In Qatar, the electricity price ( ) depends on the customer category (pro farms, residential, industrial, hotel, commercial, and government).The governm courages investment in sectors like agriculture, industry, and tourism by subsidiz tricity with different degrees.Moreover, the applied tariff for residential and gov institutions is different.The electricity tariff is in the range ($0.02 kWh, $0.09 kW In the urban region, the available installation area for the solar collector is not iden all buildings.It may depend on the area of the building, the presence of parking,

Scenarios and Configurations Assessment
In Qatar, the electricity price (P Elec ) depends on the customer category (productive farms, residential, industrial, hotel, commercial, and government).The government encourages investment in sectors like agriculture, industry, and tourism by subsidizing electricity with different degrees.Moreover, the applied tariff for residential and government institutions is different.The electricity tariff is in the range ($0.02 kWh, $0.09 kWh) [43].In the urban region, the available installation area for the solar collector is not identical for all buildings.It may depend on the area of the building, the presence of parking, and the unused area.To reflect this local context, 9 scenarios are designed (Figure 8).For the electricity, low, medium, and high tariffs represent $0.02 kWh, $0.055 kWh, and $0.09 kWh.Low, medium, and high availability for the installation area is considered consecutively representing 3000 m 2 , 6000 m 2 , and 9000 m 2 .For example, scenario 5 represents medium area availability and energy tariff.This is true for configuration 3 at a high electricity tariff.Therefore, the available installation area may impact the competitiveness of the solar cooling configuration.

Base Configuration Performance
The Grid-DCS is considered base configuration; the economic, energy, and environmental performance are as follows.Concerning the economic performance, the investment and maintenance cost remains unchanged for all scenarios since the design is composed only of an electric chiller and chilled water tank.However, the operation cost varies and depends on the electricity tariff, impacting the total cost.The base configuration is powered exclusively by grid electricity.Hence, the final energy saving (renewable energy) is considered null.The electricity generation in the power plant is mainly based on fossil fuels, and their combustion generates CO2 emissions.The yearly electricity consumption of the Grid-DCS contributes to the consumption of 16,882,130 kWh oil-equivalents of primary energy and the emissions of 3,060,981 kg CO2-equivalents (Figure 10).

Optimal Configuration (Economic Assessment)
All scenarios are simulated for the three configurations and the base Grid-DCS by varying the electricity tariff (P Elec ) and the available installation area (A). Figure 9 shows the optimal configuration based on the economic assessment (objective function).At low electricity tariffs (Scenarios 1-3), the base configuration (Grid-DCS) is more competitive than configurations 1-3.This is true for configuration 3 at a high electricity tariff.Therefore, the available installation area may impact the competitiveness of the solar cooling configuration.

Base Configuration Performance
The Grid-DCS is considered base configuration; the economic, energy, and environmental performance are as follows.Concerning the economic performance, the investment and maintenance cost remains unchanged for all scenarios since the design is composed only of an electric chiller and chilled water tank.However, the operation cost varies and depends on the electricity tariff, impacting the total cost.The base configuration is powered exclusively by grid electricity.Hence, the final energy saving (renewable energy) is considered null.The electricity generation in the power plant is mainly based on fossil fuels, and their combustion generates CO2 emissions.The yearly electricity consumption When the electricity tariff increases (medium or high), the solar cooling configurations become more competitive than the base.However, configuration 3 (T-DCS) is competitive only at high tariffs.Configuration 1 is competitive at medium electricity tariff when the available area is medium and above.This is true for configuration 3 at a high electricity tariff.Therefore, the available installation area may impact the competitiveness of the solar cooling configuration.At a low electricity tariff (Scenarios 1-3), The Grid-DCS is the optimal configuration, and the economic performance is detailed in Figure 10.For the remaining Scenarios (Figure 11), the investment and maintenance cost of the Grid-DCS system is the lowest due to the absence of an absorption chiller, hot water tank, and solar collector in the system design.However, the operation cost is the highest among the other configurations due to the grid electricity consumption.The total cost of the PV system is the lowest.The investment and maintenance cost of the PV system is less by 40% to 60% than PVT and T systems due to the absence of absorption chiller and hot water tank in the design, and superior by 60% to 180% than the Grid-DCS configuration mainly due to the PV panel installation cost.The operation cost of the PVT system is the lowest due to the heat and electricity cogeneration.The total cost of the PVT system is slightly superior to the PV system and is in the range of 2% to 4% in the case of high electricity tariff and 19% to 31% for the medium tariff.However, the operation cost of the PVT-DCS configuration is lower by 59% to 89% compared to the Grid-DCS configuration, when electricity tariff is medium, and by 30% to 90% when electricity tariff is high.The T system is more competitive than the base configuration at high electricity tariffs and medium or high available areas.The total cost is almost similar to the base configuration (−2%).This is explained by compensating for the operation cost-saving and the investment cost.At a low electricity tariff (Scenarios 1-3), The Grid-DCS is the optimal configuration, and the economic performance is detailed in Figure 10.For the remaining Scenarios (Figure 11), the investment and maintenance cost of the Grid-DCS system is the lowest due to the absence of an absorption chiller, hot water tank, and solar collector in the system design.However, the operation cost is the highest among the other configurations due to the grid electricity consumption.The total cost of the PV system is the lowest.The investment and maintenance cost of the PV system is less by 40% to 60% than PVT and T systems due to the absence of absorption chiller and hot water tank in the design, and superior by 60% to 180% than the Grid-DCS configuration mainly due to the PV panel installation cost.The operation cost of the PVT system is the lowest due to the heat and electricity cogeneration.The total cost of the PVT system is slightly superior to the PV system and is in the range of 2% to 4% in the case of high electricity tariff and 19% to 31% for the medium tariff.However, the operation cost of the PVT-DCS configuration is lower by 59% to 89% compared to the Grid-DCS configuration, when electricity tariff is medium, and by 30% to 90% when electricity tariff is high.The T system is more competitive than the base configuration at high electricity tariffs and medium or high available areas.The total cost is almost similar to the base configuration (−2%).This is explained by compensating for the operation cost-saving and the investment cost.At low electricity tariff (Scenarios 1-3), The Grid-DCS is the optimal configuration based on the economic performance, and the energy performance is detailed in Figure 10. Figure 12 shows renewable energy use saving performance for the remaining Scenarios (4-9).The equivalent energy saved (gas or electricity) depends on the renewable energy generated by the installed area.The gas equivalent saved is proportional to the produced heat and depends on the solar radiation, thermal efficiency, and the installed area.The electricity saved is equal to the electricity generated and depends on the solar radiation, electric efficiency, and the installed area.The hybrid solar technology generates both heat and electricity, which is not the case for photovoltaic and thermal.Consequently, the renewable fraction is superior to PV and T systems.The maximum renewable energy fraction achieved by the PVT, PV, and T systems are respectively 94.70%, 58.30%, and 79.26% for the high installation area scenarios.The minimum fraction for PVT and PV systems are, respectively, 54.82% and 19.43% for the low installation area scenarios.The minimum for the T system is 63.23% for the medium installation area scenario (Figure 13).At low electricity tariff (Scenarios 1-3), The Grid-DCS is the optimal configuration based on the economic performance, and the energy performance is detailed in Figure 10. Figure 12 shows renewable energy use saving performance for the remaining Scenarios (4-9).The equivalent energy saved (gas or electricity) depends on the renewable energy generated by the installed area.The gas equivalent saved is proportional to the produced heat and depends on the solar radiation, thermal efficiency, and the installed area.The electricity saved is equal to the electricity generated and depends on the solar radiation, electric efficiency, and the installed area.The hybrid solar technology generates both heat and electricity, which is not the case for photovoltaic and thermal.Consequently, the renewable fraction is superior to PV and T systems.The maximum renewable energy fraction achieved by the PVT, PV, and T systems are respectively 94.70%, 58.30%, and 79.26% for the high installation area scenarios.The minimum fraction for PVT and PV systems are, respectively, 54.82% and 19.43% for the low installation area scenarios.The minimum for the T system is 63.23% for the medium installation area scenario (Figure 13).At low electricity tariff (Scenarios 1-3), The Grid-DCS is the optimal configuration, and the environmental performance is detailed in Figure 10.At low electricity tariff (Scenarios 1-3), The Grid-DCS is the optimal configuration, and the environmental performance is detailed in Figure 10.At low electricity tariff (Scenarios 1-3), The Grid-DCS is the optimal configuration, and the environmental performance is detailed in Figure 10.
Figure 14 shows the environmental performance for the remaining Scenarios (4-9).The Grid-DCS system has the highest primary energy use and global warming potential.Unlike, the hybrid technology has the lowest primary energy use and global warming potential.For high available installation areas, the PVT system reduces the environmental impact by 89.5% compared to the Grid-DCS system, by 59.74% in medium installation area scenarios, and by 30.33% in the case of low installation area scenarios.The PV system reduces the environmental impact by 58.3%, 38.87%, 19.43% when the installation area is high, medium, and low.The T system reduces the environmental impact by 43.76% and 25.71% when the installation area is high and medium.
Energies 2022, 15, x FOR PEER REVIEW 20 of 24 Figure 14 shows the environmental performance for the remaining Scenarios (4-9).The Grid-DCS system has the highest primary energy use and global warming potential.Unlike, the hybrid technology has the lowest primary energy use and global warming potential.For high available installation areas, the PVT system reduces the environmental impact by 89.5% compared to the Grid-DCS system, by 59.74% in medium installation area scenarios, and by 30.33% in the case of low installation area scenarios.The PV system reduces the environmental impact by 58.3%, 38.87%, 19.43% when the installation area is high, medium, and low.The T system reduces the environmental impact by 43.76% and 25.71% when the installation area is high and medium.

Conclusions, Limitations, and Future Research
In the hot climatic regions, the abundance in solar irradiation and high cooling load patterns form a synergy that can be exploited for an ideal sustainable integrated system.Various solar and cooling technologies exist, implying different possible configurations that satisfy the cooling demand.Therefore, this paper focused on optimizing and assessing solar cooling systems in comparison to the conventional Grid-based district cooling system.The optimal solution leads to the best sizing of the system components, including the technology selection and the optimal production policy in terms of energy generation, hot and chilled water production, and storage.

Conclusions, Limitations, and Future Research
In the hot climatic regions, the abundance in solar irradiation and high cooling load patterns form a synergy that can be exploited for an ideal sustainable integrated system.Various solar and cooling technologies exist, implying different possible configurations that satisfy the cooling demand.Therefore, this paper focused on optimizing and assessing solar cooling systems in comparison to the conventional Grid-based district cooling system.The optimal solution leads to the best sizing of the system components, including the technology selection and the optimal production policy in terms of energy generation, hot and chilled water production, and storage.
The assessment and the comparison of the four configurations is performed based on the economic, renewable energy use, and environmental performances.The electricity tariff and available installation area shape the local building context.Therefore, nine scenarios are simulated by varying these parameters.The analysis made in this paper shows that: - The electricity tariff and the available installation area impact the integration of solar energy.When electricity tariff is medium or high, the more installed the area is, the better it is for the total cost, renewable energy use, and the environmental impact.-At a low electricity tariff rate (subsidized), the conventional Grid-based DCS remains the most competitive among the other configurations.However, the environmental performance is the worst.The annual assessment of 10,000 kW district cooling system powered with grid-electricity shows that 16,882,130 kWh oil-equivalent of primary energy is consumed, which impacts the environment by 3,060,981 kg CO 2 -equivalent.- The PV-DCS configuration is the most competitive among solar-assisted cooling systems.The annual total cost is less by 2% to 4% than PVT-DCS configuration when electricity tariff is high, and less by 19% to 31% when electricity tariff is medium.It reduces the environmental impact compared to the Grid-DCS configuration by 58.3%, 38.87%, 19.43% when the installation area is high, medium, and low.- The PVT-DCS configuration has the lowest operation cost among solar-assisted cooling systems and is lower by 59% to 89% than the Grid-DCS configuration, when electricity tariff is medium, and by 30% to 90% when electricity tariff is high.The PVT system also has the best environmental performance by decreasing the primary energy use and the global warming potential by 89.5% compared to the Grid-DCS system, by 59.74% in medium installation area scenarios, and by 30.33% in the case of low installation area scenarios.- The T-DCS configuration is more economically competitive than the conventional Grid-DCS configuration only at high electricity tariffs and medium to high available installation areas.The total system cost is slightly less than the Grid-DCS configuration (−2%).Despite the high renewable energy fraction in the range 63.23% and 79.26%, respectively, for medium and high installation areas, the environmental performance is relatively lower than PVT and PV systems.The T system reduces the environmental impact by 43.76% and 25.71% when the installation area is high, and medium compared to the Grid-DCS configuration.

Limitations and Future Research
The proposed model has some limitations.The model assumes that the chiller efficiency is constant during the operation.In general, the efficiency outcomes can be varied.Therefore, the energy consumption pattern can change, exhibiting nonlinear effects in terms of energy requirements and the cooling load supply.However, this would increase the complexity of the model.Therefore, further research can focus on changing the solution approach by decomposing the problem and developing a heuristic to solve the sub-problems.
It should be further noted that the exploitation of solar energy is based on solar irradiation and weather conditions, and these factors are stochastic by nature.Therefore, the model can be extended by considering the probabilistic distribution of irradiance and its impact on electricity and thermal energy generation.
The current model does not consider heating and cooling losses.Therefore, further research can extend to heat losses to establish a more efficient configuration.
The proposed model focuses on the system's conceptual design, which is one of the earlier stages of the system life cycle development.Further research can focus on the system measurement and the predictive control aspects when the system development reaches advanced stages, such as preliminary and final designs.

Figure 1 .
Figure 1.Integrated solar collector and mix-technologies DCS Diagram.

Figure 2 .
Figure 2. Framework of the developed model.

Figure 2 .
Figure 2. Framework of the developed model.

Figure 3 .
Figure 3. Structure of used data.

Figure 3 .
Figure 3. Structure of used data.

Figure 4 .Figure 4 .
Figure 4. Cooling load and solar for a typical summer day.

Figure 5 .
Figure 5. System's electricity flow for a typical August day.

Figure 5 .
Figure 5. System's electricity flow for a typical August day.

Figure 5 .
Figure 5. System's electricity flow for a typical August day.
PE sys = ∑ t CED Gas G Bo t + ∑ t CED Elec (E Grid→Vc The global warming potential (Kg CO 2 -equivalents) is expressed as the sum of the global warming potential for each nonrenewable energy consumed.Similar to the calculation of the PE sys in the GWP sys , the CED is replaced by the GWP factor for each type of non-renewable energy.(GWP Gas = 0.228 and GWP Elec = 0.524).GWP sys = ∑ t GWP Gas G Bo t + ∑ t GWP Elec (E Grid→Vc t − E Col→Grid t

Table 3 .
Solar collectors' parameters and set size.

Table 3 .
Solar collectors' parameters and set size.

Table 4 .
Optimal design characteristics and costs for base configuration and configuration 1.
The annual generated electricity equals 1,930,012 kWh, which is considered energy saved.

Table 5 .
Annual performances of configuration (1) and base configuration.

Table 5 .
Annual performances of configuration (1) and base configuration.