Evaluation and Optimization of the Annual Performance of a Novel Tri-Generation System Driven by Geothermal Brine in O ﬀ -Design Conditions

: The di ﬀ erence in heating or cooling to power ratio between required demands for district networks and the proposed tri-generation system is the most challenging issue of the system conﬁguration and design. In this work, an adjustable, novel tri-generation system driven by geothermal resources is proposed to supply the thermal energies of a speciﬁc district network depending on ambient temperature in Germany. The tri-generation system is a combination of a modiﬁed absorption refrigeration cycle and a Kalina cycle using NH 3 -H 2 O mixture as a working ﬂuid for the whole tri-generation system. A sensitive analysis of o ﬀ -design conditions is carried out to study the e ﬀ ect of operational parameters on the system performances prior to optimizing its performance. The simulation show that the system is able to cover required heating and cooling demands. The optimization is applied considering the maximum exergy e ﬃ ciency (scenario 1) and minimum total exergy destruction rate (scenario 2). The optimization results show that the maximum mean exergy e ﬃ ciency in scenario 1 is achieved as 44.67% at the expense of 14.52% increase in the total exergy destruction rate in scenario 2. The minimum mean total exergy destruction rate in scenario 2 is calculated as 2980 kW at the expense of 8.32% decrease in the exergy e ﬃ ciency in scenario 1.


Introduction
Due to the non-fluctuating renewable and environmentally friendly nature of geothermal heat sources, it is inevitable to utilize these sources to reduce the consumption of fossil fuels and environmental pollution [1]. For applying geothermal sources to produce power, cooling and heating in different seasons, various multi-generation systems are investigated [2]. It has been observed that co-generation [3] and multi-generation systems [4,5] attract the interest not only of researchers but also of operators and energy suppliers because of their substantial energy conservation [6,7] and economic aspects [8]. Recently, some studies have been carried out to examine co-generation systems for producing power/cooling or power/heating using geothermal heat sources. In addition, some configurations of tri-generation systems have been presented to provide the demand of power, cooling and heating [9,10] with low-grade heat sources like geothermal heat sources. Coskun et al. [11] investigated a tri-generation system driven by geothermal heat source of 173 • C. The system was studied on two different main groups for hot and cold seasons. The exergy efficiency in winter with an ambient temperature of 10.1 • C was calculated as 52.8%. In summer with an ambient temperature of 22.3 • C the exergy efficiency of the system was 49.5%. In addition, in comparison to the single power • to define heating and cooling demand in winter and summer conditions (different ambient temperatures) to evaluate the system performance using operational data of an existing geothermal district heating network in the German Molasse basin [24]. In particular, the geothermal source, mean ambient temperatures and heating/cooling demands are based on this case; • to validate the power part (Kalina cycle) and to verify the cooling part in order to confirm the validity of the simulation model; • to determine the main operational parameters to analyze the system; • to develop and apply an off-design model of the tri-generation system; • to apply the quasi-stationary simulation of a tri-generation system in different weather conditions to fulfill the heating and cooling demands; • to obtain mean annual system performances; • to optimize the system considering two criteria, maximum exergy efficiency and minimum total exergy destruction rate.

Description of the System
The schematic diagram of the proposed tri-generation system is shown in Figure 1. The system is able to simultaneously generate heating, power and cooling for various purposes in different seasons. The system consists of three parts: the power section which produces electricity, the adaptable cooling part which is applied as absorption refrigeration in summer days and heat pump in winter days (by switching the evaporator) and the domestic heat exchanger (DHE) is used to supply heat demand. As shown in Figure 1, the geothermal hot water (state 21) is divided into two parts; one part provides the required heat input to satisfy the heat demand (state 24) via the DHE. The other part is coupled to the power section by the boiler (BOI) (state 22). The required heat for the generator (GEN) is supplied by stream 26 which is the combination of the two returning geothermal water streams coming from DHE (state 25) and BOI (state 23). The description of the tri-generation system is started with the ammonia-water solution in the power section. A saturated liquid ammonia-water solution (state 4) is pumped through pump 2 (PU2) (state 8) to the high pressure of the system. The stream is delivered to separator (SEP) after heating in the boiler via geothermal water (state 9). In the SEP, the liquid-vapor stream is separated to saturated liquid (state 10) and saturated vapor phases (state 11). The stream 12 is expanded via the turbine (TUR) and is mixed with the stream leaving the solution heat exchanger (SHE) via expansion valve 2 (EX2). This stream (stream 15) is led through CON 2 and then to absorber (ABS) (stream 17) through PU3. The ammonia-water solution at the outlet of the ABS (state 1) is forced to the medium pressure level (state 2) by PU1 and is led to SHE. This stream at the outlet of the SHE (stream 3) which is heated in the SHE via saturated liquid (state 11) from SEP. The stream 3 goes to GEN to heat with the geothermal water following to divide into two saturated liquid and vapor streams (stream 4 and 5). The vapor stream coming from the GEN (state 5) is cooled through rectifier (REC) and again the stream is divided into liquid and vapor phase to provide a strong ammonia solution (stream 7). The pure ammonia is saturated in the CON 1 (stream 18) and is led through the evaporator (EVA) via EX 1 . In the EVA, the pure ammonia solution is heated by absorbing heat and can produce a cooling demand in the summer condition (stream [39][40]. The stream coming from EVA (stream 20) leads to ABS and is mixed with stream 17 via PU3 and the proposed system is completed. The heat sink streams (streams 31-36) are used as coolant and the water-ethylene glycol streams 39 and 40 are utilized to cover the cooling demand. The water streams 29 and 30 are used to produce heating demand. Moreover, it is worthwhile to indicate that in the proposed system, two different evaporators are applied in summer and winter. In summer days, when the tri-generation mode is required, the evaporator for cooling is applied to the system. In winter days, when no cooling demand is required, the geothermal fluid coming from GEN (stream 27) is led through the evaporator for heating. In this regard, the useful productions in the tri-generation system are outlined as follows: heating via domestic heat exchanger power generation via turbine cooling via evaporator

Assumption
For a process simulation and a detailed analysis of the proposed system, each component of the cogeneration cycle is considered as a control volume and a mathematical model is established by applying the conservation of mass and energy balances as well as the exergy principle. Simulation is carried out by the software Engineering Equation Solver (EES) [25]. The fluid properties of ammoniawater solution and water are estimated in EES by implementing correlations of Ibrahim [26] and Wagner [27], respectively. In addition, the thermodynamic and transport properties of waterethylene glycol are applied in EES by adapting the correlations of Fujuita [28] and Lee [29].
Some assumptions are made to simplify the mathematical model [30]: -Pressure and heat losses in the heat exchangers and pipelines are neglected. - The isotropic efficiency for pumps is considered as constant. - The changes of kinetic and potential energies are negligible. - The state of leaving streams of absorber, condenser, generator and rectifier is saturated. - In the cold season when heating demand is required, the extracted stream of geothermal fluid from the generator is used in the evaporator to supply the needed heat (stream [27][28] in the evaporator in the absorption heat pump cycle. In the summer days, the stream 27 goes to reinjection well and the cooling demand is produced in the evaporator for cooling.
The basic assumption and assumed boundary conditions for the simulation are summarized in Table 1.

Assumption
For a process simulation and a detailed analysis of the proposed system, each component of the cogeneration cycle is considered as a control volume and a mathematical model is established by applying the conservation of mass and energy balances as well as the exergy principle. Simulation is carried out by the software Engineering Equation Solver (EES) [25]. The fluid properties of ammonia-water solution and water are estimated in EES by implementing correlations of Ibrahim [26] and Wagner [27], respectively. In addition, the thermodynamic and transport properties of water-ethylene glycol are applied in EES by adapting the correlations of Fujuita [28] and Lee [29].
Some assumptions are made to simplify the mathematical model [30]: -Pressure and heat losses in the heat exchangers and pipelines are neglected. - The isotropic efficiency for pumps is considered as constant. - The changes of kinetic and potential energies are negligible. - The state of leaving streams of absorber, condenser, generator and rectifier is saturated. - In the cold season when heating demand is required, the extracted stream of geothermal fluid from the generator is used in the evaporator to supply the needed heat (stream [27][28] in the evaporator in the absorption heat pump cycle. In the summer days, the stream 27 goes to reinjection well and the cooling demand is produced in the evaporator for cooling.
The basic assumption and assumed boundary conditions for the simulation are summarized in Table 1. Table 1. Input data and boundary conditions used for design the system and process the simulation.

Off-Design Model
The operating conditions of the system are influenced by the performance of the heat exchangers and turbine [33,34]. In this regard, the different steps of off-design calculation are indicated in Figure 2. As shown, the first step in off-design process calculation is related to obtain the heat transfer capacity (UA) of heat exchangers in the design mode which is outlined in Table 1. The heat transfer rate ( . Q) in the heat exchangers is expressed as [35]: where U is the overall heat transfer coefficient, the area of the heat exchanger is A and ∆T LMTD represents the logarithmic mean temperature difference between the hot and cold streams in the heat exchanger. In the next step, the off-design behavior of the heat exchangers is considered according to the method presented by Manente [34]. In off-design conditions, the variations in mass flow rates from the design conditions bring about a change in the UA. The heat exchange capacity is modified considering the actual mass flow rate: The heat exchanger configuration determines the exponent m. The exponent m is set to 0.44 for the air-cooled condenser, 0.66 for the generator, solution heat exchanger and evaporators, 0.15 for the boiler [36]. After obtaining the off-design characteristics such as UA for the considered heat exchangers, calculation of the turbine isentropic efficiency is the next step. The isentropic efficiency of the turbine is considered as a function of enthalpy drop and volumetric flow rate at the outlet of the turbine. The off-design behavior of the isentropic efficiency of the turbine depends on two coefficients r h and r v according to [37]. The coefficients depend on the outlet volume flow rate and the enthalpy drop in the turbine: where h is the enthalpy and exponents design means the calculated value at design [37]: Appl. Sci. 2020, 10, 6532 6 of 20 r v = ([(−0.21r VT + 1.117) × r VT − 2.533] × r VT + 2, 588) × r VT + 0.038 (6) where .
V shows the volumetric flow rate at the outlet of the turbine. Therefore, the isentropic efficiency of the turbine in off-design can be calculated as [32]: Appl. Sci. 2020, 10 where shows the volumetric flow rate at the outlet of the turbine. Therefore, the isentropic efficiency of the turbine in off-design can be calculated as [32]: It is important to note that the off-design values of , and can take a value higher than the design mode values. In addition, the UA value and turbine isentropic efficiency, the thermodynamic properties of different states have been calculated. The whole off-design case is calculated in an iterative process which ends when the simulation shows convergence. The main criteria of convergence are related to the equality of the thermodynamic properties such as temperature, pressure and ammonia concentration in the outlet of different heat exchangers and the inlet of other heat exchangers. The convergence is obtained within 1%. The proposed system is designed under operational parameters under base case which is presented in Table 1. The design case results are calculated for the heating demand of 11.66 MW in the ambient temperature of −5.67 °C. It is important to note that the off-design values of η s,tur and r VT can take a value higher than the design mode values. In addition, the UA value and turbine isentropic efficiency, the thermodynamic properties of different states have been calculated. The whole off-design case is calculated in an iterative process which ends when the simulation shows convergence. The main criteria of convergence are related to the equality of the thermodynamic properties such as temperature, pressure and ammonia concentration in the outlet of different heat exchangers and the inlet of other heat exchangers. The convergence is obtained within 1%. The proposed system is designed under operational parameters under base case which is presented in Table 1. The design case results are calculated for the heating demand of 11.66 MW in the ambient temperature of −5.67 • C.
The heat exchange capacities (UA) and mass flow rates of all heat exchangers in the system for design case are shown in Table 2. In addition, in design conditions the volumetric flow rate at the turbine outlet is calculated to be 3.376 m 3 /s.

Heating and Cooling Demand Profiles
In order to analyze the proposed cycle, real heating data of an existing geothermal district heating network in the German Molasse Basin [24] and estimated cooling demand are applied. The VDI 4655 method [38] is applied to estimate the annual system performance. In this method, 10 typical days are analyzed. This procedure is distinguished from seasonal, i.e., summer (S), winter (W) and transition (T), and from the typical days, i.e., workdays (W) and Sundays (S). In this category, Saturdays are selected as the workday. Furthermore, fine (F) and cloudy (C) days are differentiated to reach acceptable results. However, in the summer days, no distinction between fine and cloudy days (X) is made. These 10 typical days categories are; WWF, WWC, WSF, WSC, TWC, TSC, TSF, TWF, SWX, SSX.
According to the test reference years (TRY), Germany is divided into various climate zones. As shown in Table 3, according to this climate zones in the mentioned category "VDI 4655", the number of typical days during the year changes. In this work, for the purpose of the annual simulation, the 10 typical days are simulated for TRY13 which is associated to the southern German Molasse Basin. The annual performances of the system are calculated by weighting these typical days considering their frequency. In order to conduct a quasi-stationary simulation, the heating and cooling demands for each typical day category are required. The profiles of corresponding heating demand for different typical days categories are developed depending on operational data of a real geothermal heat plant in the German Molasse Basin. The corresponding mean values of heating demands and the ambient temperatures for different typical days categories are illustrated in Table 4, respectively [18]. It is seen that the ambient temperature takes the lowest values in winter days and the highest values occur in the summer days. It is clear that the ambient temperature in cloudy days is basically lower than the fine days. For heating demand values which are listed in Table 4, it is observed that the heating demand has higher values in winter days in comparison to the summer and transition days. In addition, the heating demand in cloudy days, especially in the afternoon is considerably higher than the one in fine days. Furthermore, comparing the heat demand on Sundays and workdays shows that the heating demand is generally lower on Sundays and has a small variation. The cooling demand is not available for specific real application; therefore, it is calculated with the following assumption. The cooling demand for an ambient temperature higher than 13 • C in four typical days according to some available real data in public buildings is assumed to be 0.33-0.85 times of the corresponding heating demand which is depending on the ambient temperature. It is noteworthy to note that this value is based on the ambient temperature in which with increasing of the ambient temperature, the corresponding value increases. Referring to Table 4, it is shown that the cooling demand takes higher values on workdays in comparison to Sundays. In addition, the cooling demand on all 4 typical days in the morning is lower than the one in the afternoon.

Exergetic Evaluation
In order to evaluate the cycle performance, the second law efficiency and exergy destruction rate are considered as main evaluation criteria. In general, in the absence of magnetic, electrical, nuclear and surface tension effects as well as ignoring the kinetic and potential exergies, the exergy rate consists of two physical and chemical exergy rates as [39]: Physical exergy can be obtained from [33]: and the chemical exergy for ammonia-water solution is expressed as [40]: .
The subscript 0 in Equation (3) indicates the dead state and T 0 shows its temperature. The standard chemical exergy for ammonia and water is set to 336.684 and 45 kJ/kmol [40].
For a detailed analysis, the exergy destruction rate is calculated. The exergy destruction rate is defined as the difference between the product exergy ( . E P ), loss exergy ( . E L ) and fuel exergy ( . E F ) for each component of the proposed system [41]: For parametric analysis, the energy and exergy efficiencies are calculated to be [42]: It is evident that the energy or thermal efficiency is the ratio of the net available products (net product power, heating and cooling rates) to the rate of thermal energy, which is absorbed from the geothermal heat source. In addition, the exergy efficiency is the ratio of the net available products exergy (net product power, heating and cooling exergies) to the rate of input exergy from the geothermal heat source. In Equations (12) and (13) . W net is the net output power of the system, . Q cooling is the cooling capacity in the summer condition, the heating capacity and the total heat input are . Q heating and . Q in , respectively. . E cooling is the generated exergy rate for cooling purpose used in the summer days, the heating exergy and the total input exergy are . E heating and . E in , respectively. These parameters are defined as: . . . . . .
The total exergy destruction rate of the system is calculated to be: on winter days: .
on summer days: .

Results and Discussion
The system is analyzed in off-design conditions by considering mass, energy and exergy balances as well as off-design model equations for each component as a control volume. In this regard, by adapting the geothermal mass flow division ratio (r geo = . m 24 / . m 21 ), the heating demand is totally covered by the geothermal resource. In addition, the sensitive analysis is applied to each typical day category displayed in Table 4 to find the operational parameters prior to optimization.

Validation and Verification
In Table 5, the real power plant parameters given in Ogriseck [16] are compared to the conducted simulations. The represented results display a good agreement between the present work and real data of the existing Kalina power plant with a deviation of less than 0.9%. The percentage of relative deviation between available data in the literature and obtained data from simulation in the present work is calculated as follows: The cooling system (absorption refrigeration system) is verified using available data reported in literature [43]. In the absorption system, the verification results are obtained with some main design parameters such as absorber, condenser and generator temperatures and ammonia concentration. Figure 3 represents the data obtained in the present work and the data by Adewusi [43] for the absorption refrigeration ammonia-water cycle. The represented results display a good agreement between the present work and literature with a deviation of less than 1%.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 11 of 20 The cooling system (absorption refrigeration system) is verified using available data reported in literature [43]. In the absorption system, the verification results are obtained with some main design parameters such as absorber, condenser and generator temperatures and ammonia concentration. Figure 3 represents the data obtained in the present work and the data by Adewusi [43] for the absorption refrigeration ammonia-water cycle. The represented results display a good agreement between the present work and literature with a deviation of less than 1%.

Parametric Analysis
A parametric study is applied to investigate the influences of decision parameters such as the turbine inlet pressure (TIP), the generator temperature (T ), the temperature differences at the Absorber (δT ) and condensers exit (δT = T − T , δT = T − T ) on the system

Parametric Analysis
A parametric study is applied to investigate the influences of decision parameters such as the turbine inlet pressure (TIP), the generator temperature (T gen ), the temperature differences at the Absorber (δT ABS ) and condensers exit (δT CON1 = T 18 − T 35 , δT CON2 = T 16 − T 37 ) on the system performance prior to optimization. It is noteworthy to mention that, in parametric analysis for off-design conditions only one parameter, among the ones outlined in Table 1, is changed and the other parameters are calculated in off-design conditions. The ammonia concentration in the power section is kept constant as 0.5 kg(NH 3 )/kg(solution).
The variations of exergy efficiency as well as total exergy destruction rate with the turbine inlet pressure are shown in Figure 4a in off-design conditions. As shown, the exergy efficiency and . E D,tot take a local maximum and minimum values, respectively. The trends of these variations can be justified by representing the variation in power production, input exergy and loss exergy in some components with the TIP. Due to the reduction in reinjection geothermal temperature, the . E in is increased with an increase in TIP. In addition, changes in TIP has no effect on the produced heating capacity and the heating exergy of the system. On the other side, the power production is maximized locally at specific values of TIP. With an increase in TIP, the turbine mass flow rate decreases, however the turbine specific work has been optimized locally. As presented in the Figure 4b, the contrary effect of these parameters occurs in a way that leads to a local optimum for the produced net power output and the exergy efficiency. Furthermore, with an increase in TIP, the loss exergy rate in the ABS, REC, CON 1 , CON 2 is increased. According to Equation (21), the total exergy destruction rate takes locally minimum values.  The influence of the generator temperature on the exergy efficiency and total exergy destruction rate in off-design mode is presented in Figure 5a. It is seen that as the T increases, the exergy efficiency is locally maximized as well as the E , is increased. With an increase in T , the E and W are decreased as shown in Figure 5b. The reduction rate of these two parameters occur in such a way that leads to a local maximum value for the exergy efficiency. The loss exergy which occurs in ABS, REC, CON1, CON2, is decreased by increasing T . In this regard, according to Equation (21) the total exergy destruction rate decreases with an increase in T .  The influence of the generator temperature on the exergy efficiency and total exergy destruction rate in off-design mode is presented in Figure 5a. It is seen that as the T gen increases, the exergy efficiency is locally maximized as well as the . E D,tot is increased. With an increase in T gen , the . E in and . W net are decreased as shown in Figure 5b. The reduction rate of these two parameters occur in such a way that leads to a local maximum value for the exergy efficiency. The loss exergy which occurs in ABS, REC, CON 1 , CON 2 , is decreased by increasing T gen . In this regard, according to Equation (21) the total exergy destruction rate decreases with an increase in T gen . rate in off-design mode is presented in Figure 5a. It is seen that as the T increases, the exergy efficiency is locally maximized as well as the E , is increased. With an increase in T , the E and W are decreased as shown in Figure 5b. The reduction rate of these two parameters occur in such a way that leads to a local maximum value for the exergy efficiency. The loss exergy which occurs in ABS, REC, CON1, CON2, is decreased by increasing T . In this regard, according to Equation (21) the total exergy destruction rate decreases with an increase in T . The influence of the temperature differences at the exit of the condensers on exergy efficiency and total exergy destruction rate in off-design mode is depicted in Figure 6. It is seen that with an increase in δT , the η and E , are increased and decreased, respectively. In addition the η is decreased and E , increased as the δT changes. The net produced power is decreased with an increase in both δT and δT . However, the input exergy to the system is kept constant in changing δT and is decreased in variation δT . The influence of the temperature differences at the exit of the condensers on exergy efficiency and total exergy destruction rate in off-design mode is depicted in Figure 6. It is seen that with an increase in δT CON1 , the η exergy and . E D,tot are increased and decreased, respectively. In addition the η exergy is decreased and . E D,tot increased as the δT CON2 changes. The net produced power is decreased with an increase in both δT CON1 and δT CON2 . However, the input exergy to the system is kept constant in changing δT CON2 and is decreased in variation δT CON1 .
Appl. Sci. 2020, 10, x FOR PEER REVIEW 13 of 20 Figure 6. The variation of exergy efficiency and the total exergy destruction rate with the temperature differences at the condensers exit.

Discussion of Parametric Analysis
The general behavior of the system on the selected parameters may be summarized as follows:  Mean total ficiency (%) Figure 6. The variation of exergy efficiency and the total exergy destruction rate with the temperature differences at the condensers exit.

Discussion of Parametric Analysis
The general behavior of the system on the selected parameters may be summarized as follows: -There is a local optimum value in the exergy efficiency at the specific value of TIP = 16.3 bar and T gen = 51.5 • C in the typical day of WWC. In the other typical day categories, the exergy efficiency has local optimum at the different values of TIP and T gen in comparison to typical day of WWC. - The total exergy destruction rate in minimized locally at the specific values of TIP = 16.53 bar in the typical day of WWC. This value of TIP is changed in the other typical day categories in order to minimize locally the total exergy destruction rate. - The general behavior of operational parameters for all typical days categories is almost the same of the typical day of WWC as shown in Figures 5-7. -In the hot days; TSF, TWF, SWX and SSX, with an increase in the δT CON1 the exergy efficiency decreases which is in reverse for other category days.

-
The values of δT CON1 and δT CON1 are limited by the ambient temperature and take different range in different typical day categories.
-The general behavior of operational parameters for all typical days categories is almost the same of the typical day of WWC as shown in In the hot days; TSF, TWF, SWX and SSX, with an increase in the δT the exergy efficiency decreases which is in reverse for other category days. - The values of δT and δT are limited by the ambient temperature and take different range in different typical day categories.

The Optimized Quasi-Steady State Results For Different Typical Days
The quasi-steady state study under off-design conditions is calculated for the proposed system to examine the system performance. Then the system is optimized based on two main criteria, exergy efficiency and total exergy destruction rate. In the first scenario, the optimization is performed for the maximum exergy efficiency. The minimizing of the total exergy destruction rate is considered in the second scenario. According to the sensitive analysis, parameters such as: the turbine inlet pressure, the generator temperature, the temperature differences at the Absorber and condensers affect the system performance. The possible range of these parameters for optimization purpose is considered as follows: maximize η exergy and minimize . E D,tot (TIP, T gen , δT CON1 , δT CON2 , δT ABS ) 11 ≤ TIP (bar) ≤ 32 It is noteworthy to mention that in all annual calculations for 10 typical days, the ammonia concentration in the power section is kept constant as 0.5 kg(NH 3 )/kg(solution). The other operational parameters (pinch point temperature difference in the boiler, hot outlet-cold inlet temperature difference in the boiler, heat sink outlet temperature at the absorber, rectifier and condensers, reinjection geothermal hot water temperature) are calculated in off-design conditions for each 10 typical days category. The optimization results are presented in Figure 7 for maximum mean exergy efficiency and in Figure 8 for minimum mean total exergy destruction rate. In addition, the values of different operational parameters as well as the other parameters, which are obtained to optimize the system, are underlined in Table A1 for winter conditions and Table A2 for summer conditions in Appendix A.
as seen in the Figures 7 and 8, the exergy efficiency in the winter days is varied between 38.79−53.21%. In the summer days, it varies between 25.2−36.16%. In addition, in the transient conditions, it is in the range of 36.08−41.16%. Furthermore, in Ref. [17] the ambient temperature was taken constant values of in summer and winter days 22.3 °C and 10.1 °C, respectively. However, in the present study it varies in different typical day categories according to an existing geothermal district heating network in the German Molasse basin, which makes the investigation more practical.   Referring to Tables A1 and A2 in the Appendix, the turbine isotropic efficiency in the off-design mode in different typical day categories varies between 56−85%, however it took almost in the previous works constant value of 85%. It is important to mention that the turbine isotropic efficiency In order to maximize the exergy efficiency based on yearly calculation in the proposed system, the δT CON2 and δT ABS take their highest values, the δT CON1 is considered its lowest value and the TIP and T gen are selected as the variable parameters. In the first scenario, the mean annual exergy and energy efficiencies are calculated to be 44.67% and 35.47%, respectively due to the weighting according to the frequency of typical days. The mean total exergy destruction rate in the first scenario is calculated to be 3412.83 kW. As indicated in Figure 7, due to the high values of heating demand in winter days in comparison to the transition and summer days, the exergy efficiency in winter days is higher. In the weather category of WWC, the higher mean exergy efficiency (53.21%) happens. In addition, the mean exergy efficiency in summer days takes the lowest value because of the lower values of heating demand. The lowest decrease of the mean exergy efficiency takes place at SSX due to the lowest value of required heating demand. The cooling demand in summer days is higher; however, the cooling exergy takes lower values in comparison to heating exergy. Besides, in the first scenario the highest and lowest mean energy efficiency occurs in the WWF and SSX with values of 53.12% and 19.07%, respectively. The lower mean energy efficiency in hot days happens due to the lower values of heat input to the system. The higher reinjection geothermal temperature brings about reduction in heat input. In the typical day of WSF, in the first scenario the system is run under part load of 90% with net produced power of 1517 kW. However, in the typical day of TSF with the net produced power of 986.1 kW the system is run under part load of 58%.
In minimizing the total exergy destruction rate, the δT CON1 , δT ABS and δT CON2 select their lowest values. The TIP and T gen are considered as the variable parameters. In this regard, in the second scenario the annual mean exergy destruction rate is 2980.32 kW. In addition, the annual mean exergy and energy efficiencies in the second scenario are calculated to be 40.95% and 42.59%, respectively. It is observed that the highest mean exergy destruction rate occurs in the typical day TWC of 3747 kW. The lowest mean exergy achieves a value of 2433 kW in the typical day TSF. In the winter days, the total exergy destruction rate in the work days is lower than the one in the Sundays, however, in the transient and summer days this value in the workdays is higher than the one in the Sundays. In addition, a comparison between cloudy and fine days in the winter days shows that in the fine days the total exergy destruction rate takes higher value than the one in the cloudy days in the winter conditions. In the transient condition, this value in cloudy days in lower than the one in the fine days. In this regard, the performances of tri-generation system have been changed with variation in the weather conditions. The system works under part load of 35.35% and 70.32% in the typical day of WWF and SWX, respectively.
Comparing two optimization cases, the reduction in the total exergy destruction rate in the second scenario is achieved at the expense of 8.32% decrease in the exergy efficiency. In addition, the higher mean exergy efficiency in the first scenario is obtained at the expense of 14.52% increase in the total exergy destruction rate. In this regard, the optimization results for minimizing the total exergy destruction rate bring about a better performance enhancement.
In detail discussion, the on-design exergy efficiency in [17] driven by geothermal heat source of 173 • C was 52.8% in winter days, however it was calculated to be 49.5% in summer days. However as seen in the Figures 7 and 8, the exergy efficiency in the winter days is varied between 38.79−53.21%. In the summer days, it varies between 25.2−36.16%. In addition, in the transient conditions, it is in the range of 36.08−41.16%. Furthermore, in Ref. [17] the ambient temperature was taken constant values of in summer and winter days 22.3 • C and 10.1 • C, respectively. However, in the present study it varies in different typical day categories according to an existing geothermal district heating network in the German Molasse basin, which makes the investigation more practical.
Referring to Tables A1 and A2 in the Appendix A, the turbine isotropic efficiency in the off-design mode in different typical day categories varies between 56−85%, however it took almost in the previous works constant value of 85%. It is important to mention that the turbine isotropic efficiency in the off-design calculation in the second scenario is almost 5−20% lower than the one in the first scenario. In addition, the pinch point difference temperature in the boiler in previous works always took a constant value of 10 • C. However, as shown in Tables A1 and A2, it varies between 14.22−25.47 • C for maximizing the exergy efficiency. It is obvious that the pinch point difference temperature in the boiler takes higher values in order to minimizing the total destruction rate in comparison to maximizing the exergy efficiency. On the other hand, the temperature difference of the heat sink streams (δT air ) in the previous studies almost was taken 10 • C. These values in the present work and referring to Tables A1 and A2, are calculating in the off-design conditions and vary in the range of 3.5−19.2 • C. It is shown that the heat sink temperature difference in the REC (δT air,REC ) is decreases with a reduction in the heating demand rate. Additionally, these values in the second scenario are higher than the ones in the first scenario. The heat sink temperature difference of the ABS (δT air,ABS ) and CON 1 (δT air,CON1 ) in maximizing the exergy efficiency is almost higher than the one in minimizing the total destruction rate. However, the heat sink temperature difference in the CON 2 (δT air,CON2 ) in the first scenario is lower than the one in the second one.

Conclusions
A novel tri-generation system using ammonia-water solution as working fluid driven by geothermal heat source is studied with the aim to fulfill corresponding heating and cooling demand in winter and summer conditions in Germany. The sensitive analyses are performed to study the influence of various optional parameters on the system performance. In addition, the result of parametric analysis for 10 different typical days is investigated in off-design modes.
The following main conclusions can be drawn according to the corresponding issues mentioned in the introduction:

•
The required heating demand is covered by the geothermal mass flow division ratio in both winter and summer conditions. This value is decreased with a reduction in heating demand.

•
There is a local optimum for the exergy efficiency and the total destruction rate at specific values of Turbine inlet pressure.

•
With an increase in generator temperature, there is a local maximum in exergy efficiency, however, the total exergy destruction rate is decreased.

•
Under optimization condition, the maximum annual exergy efficiency is obtained as 44.67% at the expense of 14.52% increase in the total exergy destruction rate in comparison to second scenario.

•
The ooptimization results show the minimum mean total exergy destruction rate of the system is calculated as 2980 kW at the expense of 8.32% decrease in the exergy efficiency in the first scenario.

•
Considering both exergy efficiency and total exergy destruction rate in off-design condition, the best enhancement through the optimization process is achieved on the typical day WWC with 53.21% of exergy efficiency and 2570 kW of total exergy destruction rate.

•
Comparing the optimization results for maximum exergy efficiency and minimum total exergy destruction rate depicts a better performance improvement obtained in the second scenario.
In future studies, designing the heat exchangers will be applied to compare the off-design results. In addition, the annual economic analysis will be considered, as well.
Author Contributions: All authors contributed to this work by collaboration. M.A.K. is the main author of this manuscript. F.H. had the initial idea and assisted in the conceptual design of the study as well as in the writing of the manuscript. D.B. supervised the whole project. All authors have read and agreed to the published version of the manuscript.

Funding:
The funding from the Bavarian State Ministry for Education, Science and the Arts in the framework of the project "Geothermie-Allianz Bayern" is gratefully acknowledged. This publication was funded by the German Research Foundation (DFG) and the University of Bayreuth in the funding programme Open Access Publishing.

Conflicts of Interest:
The authors declare no conflict of interest.