Theoretical Analysis of a Biomass-Driven Single-Effect Absorption Heat Pump for Heating and Cooling Purposes

: Renewable energy exploitation in the building sector can lead to significant energy savings and carbon dioxide emission avoidance. The objective of this study is the detailed investigation of a biomass-driven absorption heat pump for heating and cooling. The heat pump is practically a single-effect absorption chiller operating with the Lithium-bromide/water solution and it has been properly modified for heating production during the winter. This system is a novel one and its combination with a biomass boiler was examined for the first time, especially for covering both heating and cooling needs. For the present study, a typical building in Athens, Greece, with a 400 m 2 floor area is selected to be coupled with the suggested heating/cooling configuration. The analysis was conducted by using TRNSYS software for t he estimation of the building’s thermal loads and with the Engineering Equation Solver for determining the heat pump behavior. According to the results, the yearly biomass consumption is found to be 3.76 tons covering a heating demand of 9136 kWh and cooling demand of 8168 kWh. The seasonal energy cooling performance was found to be 0.751, while the seasonal energy heating performance was at 1.307. Moreover, the proposed configuration was found to have economic and environmental benefits compared to conventional units with an oil boiler and heat pump for cooling. Specifically, the present system leads to 10.8% lower operational costs and 4.8% lower primary energy demand, while there are significant amounts of CO 2 avoidance.


Introduction
Renewable energy utilization is a vital weapon to face the challenges of high electricity prices, global warming, and the increasing energy needs of humanity [1]. The building sector is responsible for about 40% of the worldwide energy demand [2] and thus the utilization of renewables in buildings is a proper way to reach sustainability [3]. Solar irradiation, biomass, wind energy, and geothermal energy are the most representative renewable energy sources that can be exploited in the building sector in various ways. Among them, the use of biomass gives some important advantages [4,5] which make it a very competitive renewable energy source [6]. More specifically, biomass can produce thermal energy of high temperatures with the proper boiler design, and it can be used when there is a demand, something that solves the storage issue that other renewable energies face (e.g., photovoltaics and wind turbines).
In the literature, there are plenty of studies on the utilization of renewable energies in the building sector. Soltero et al. [7] suggested that the use of biomass for district heating purposes is a sustainable choice and they conducted a detailed sustainability analysis for different cities in Spain. Chen et al. [8] found that the incorporation of heat pumps in a district heating system with biomass can enhance the overall performance by about 17%. Zhang et al. [9] investigated a hybrid solar/biomass heating system and they found solar coverage at 63.3% and system exergy efficiency at 16.2%. Nami et al. [10] studied a solar/biomass-fed trigeneration system for heating, cooling, and power production. Thermal photovoltaics, a biomass boiler, power blocks for electricity, heat exchangers for heating, and an absorption chiller for cooling production, are incorporated into the cycle. The system exergy efficiency was found close to 34%. Bellos et al. [11] examined a solar/biomass multigeneration unit for electricity, cooling, and heating production at two temperature levels. They coupled an organic Rankine cycle with a compression heat pump and they found the system exergy efficiency to be 21.8%. Tsimpoukis et al. [12] studied a solar/biomass polygeneration system based on the use of a supercritical CO2 cycle. They found that the total energy efficiency of the system reaches up to 164%, something that is justified by the partial operation of the system as a heat pump. Bellos and Tzivanidis [13] studied a biomass-based polygeneration system based on a novel CO2 supercritical cycle. They calculated the system's energy efficiency at 78.1% and the respective exergy efficiency at 26.3%. Xing and Li [14] studied a biomass/geothermal polygeneration unit for power, heating, cooling, and hydrogen production. The energetic efficiency was found to be 79.5%, while the respective exergetic was at 17.9%. Rezaei et al. [15] studied a biomass/geothermal heating configuration from a financial point of view. They highlighted the need for a detailed economic analysis, as well as the need for using a storage device for storing heat when both heat sources provide heat in the system.
Moreover, it is important to highlight the advantage of thermal renewable energies (e.g., biomass, solar energy, geothermal energy) due to the ability for thermal storage which leads to grid flexibility solution. Advanced techniques have been developed aiming to store thermal energy in a compact and efficient way. The use of phase change materials (PCM) is an important technology that can lead to a significant performance in the case combination with other devices. There are examples that combine PCM with storage tanks [16], finned storage packed tubes [17], heat exchangers [18], batteries [19], etc. Additionally, the PCM can be used in applications for cooling storage [20] which is an important option for the sustainability of future buildings.
The present study investigates a configuration for both heating and cooling in buildings, based on a biomass-powered absorption heat pump operated with LiBr/H2O working substance. Biomass pellets are burned in a boiler to input driving heat to the absorption machine at a temperature of 110 °C . In the summer period, the heat pump evaporator can ensure a cooling temperature of around 5 °C . In cold weather, space heating is supplied at a temperature of 50 °C from the absorber and the condenser, while on very cold days the boiler heats directly the building. Thus, heating and cooling demands during a whole year can be covered by a dual-purposed biomass-driven heat pump, occasionally backed by the biomass boiler, which constitutes a novelty. The energy, economic and environmental benefits of the suggested design were investigated. A residential building with a floor area of 400 m 2 in Athens, Greece, was selected as a potential implementation site. The building thermal loads were analyzed by a purposely developed TRNSYS model [21]. The thermodynamic model of the absorption heat pump employed the Engineering Equation Solver [22].

The Examined Heating/Cooling System
The first step in this study is the presentation of the single-effect absorption heat pump which is given in Figure 1. This device operates with LiBr/water working pair and it is a single-stage machine. The heat input of high temperature (e.g., 110 °C ) is given in the generator from the boiler and also heat input of low temperature is given from the evaporator. The other devices (condenser and absorber) reject heat from the system to the ambient or the building depending on the operating mode. During the summer, the evaporator absorbs energy from the building and put it in the device, while during the winter the evaporator absorbs heat from the ambient. More specifically, the heat input in the generator makes possible the high-pressure steam evaporation (state point 3) which feeds the condenser where heat is rejected and the outlet (state point 6) is a high-pressure saturated liquid. This quantity is expanded in the throttling valve without losing energy and it becomes a low-pressure liquid/vapor mixture (state point 7) that feeds the evaporator. Inside the evaporator device, saturated steam of low pressure is created (state point 8) with heat absorbance. This quantity goes into the absorber where there is a LiBr/water weak solution of low pressure. In this device, the steam is absorbed by the LiBr crystals, and heat is rejected. The weak solution (state point 1) leaves the absorber, then with a solution pump reaches high pressure (state point 12) and goes to the solution heat exchanger where it is warmed up (state point 2) before the inlet to the generator. On the other side of the circuit, the strong solution leaves the generator (state point 4), then goes to the heat exchanger where it rejects heat to the weak solution and finally it reaches a colder high-pressure strong solution (state point 45). This quantity is expanded without losing energy in a throttling valve (state point 5) and goes into the absorber in order to close the total cyclic process. Figure 2 depicts the suggested unit for heating and cooling with the utilization of biomass. There are three modes of the system, as shown in Figure 2. More specifically, Figure 2a shows the cooling mode where the evaporator of the heat pump feeds by cooling the building, while the biomass boiler feeds the generator with heat input. Thermal oil is used as the heat transfer fluid from the boiler to the generator due to the high-temperature levels. Specifically, the generator temperature is selected at 110 °C in order to be able the operation during the summer period, while heat rejection temperatures in the condenser and the absorber were set at 10 K over the ambient temperature. In this scenario, the evaporator temperature is set at 5 °C in order to be able to feed a cooling system with fan coils.  In the heating mode, there are two cases depending on the ambient temperature. More specifically, the heat pump absorbs heat from the ambient during the winter in order to transfer it, after the upgrade, in the building. However, the water in the cycle cannot go lower than 4 °C , due to freezing reasons. Thus, in low ambient temperatures, the absorption heat pump does not operate and then, the biomass boiler feeds directly the fan coils of the building. Assuming a pinch point temperature of 5 K [23], the critical ambient temperature that separates the heating operating modes is equal to Tam,cr = 9 °C . More specifically, when the Tam ≥ 9 °C , then the heat pump operates and provides heating in the building from the condenser and the absorber. On the other hand, when Tam < 9 °C , the heat pump is not possible to operate and the boiler gives direct heat input in the building. Figure 2b depicts the heating scenario where the condenser and the absorber reject heat to the building to cover its thermal needs, while the evaporator absorbs heat from the ambient. The generator operates at 110 °C in order to give the proper heat input to the system. Thermal oil is used as the heat transfer fluid for transferring thermal energy from the boiler to the generator. Figure 2c shows the system with the operation of the boiler without the use of the heat pump. The boiler gives input directly to the fan coils of the building because the heat pump cannot operate when the ambient temperature is lower than 9 °C . In this scenario, the performance of the system is reduced and thus the climate conditions of the examined region play a significant role in the energy behavior of the system in the heating mode. In other words, in cold climates, the present system is not so effective because (i) the ambient temperature is lower than 9 °C for a great percentage of the year, and (ii) the cooling loads are restricted during the summer. However, this system is an interesting choice for climate conditions such as in the examined study case in Athens, Greece.
In the present work, the boiler has been selected to have a nominal efficiency of 90% [11]. The higher heating value (HHV) of the biomass is 18467 kJ/kg which corresponds to typical pellets [24]. The dry composition of the selected biomass is 49.7% C, 6.0% H, 41.0% O, 1.7% N, 0.0% S, 0.1% Cl, 1.5% Ash [24], while the as received has 15% moisture [11]. Table 1 includes the main information regards the system analysis.

The Examined Building
A building with dimensions (20 m × 20 m × 3 m) is examined in the present work. This building is located in Athens, Greece, and it has a floor area of 400 m 2 and a volume of 1200 m 3 . It has four external walls oriented in four directions (south, west, north and east), while the south wall includes a window of 10 m 2 area. The ground slab comes in touch with the ground while the roof is in touch with the ambient air. Additionally, all the external walls are in touch with the external air and there are no shadings from other buildings or obstacles. Table 2 includes the materials of the structural components, as well as the respective material thermal properties. More details regarding the system simulation are included in Table 3. The set point temperature levels were selected at 20 °C in the winter period and 26 °C in the summer period, as typically selected in the residential buildings. The infiltration and natural ventilation rates were selected at 1 air change per hour, while the internal gain of the building for appliances and lighting was set at 7 W/m 2 . The aforementioned parameters were selected to have typical values ac-cording to Greek legislation [25]. Additionally, it was selected to have 7 occupants inside the building and every occupant produces 100 W according to ISO 7730 [26].  The structural materials have the given composition as Table 2 indicates. Using the data of Ref. [25], the heat convection coefficient in the indoor space was selected at hout = 7.69 W/m 2 K and in the outdoor environment at hin = 25 W/m 2 K, which are the usually used values. Therefore, the U-value of the ground was calculated at 0.304 W/m 2 K, of the roof at 0.318 W/m 2 K, and of the external walls at 0.365 W/m 2 K. The window has double-glazing with a U-value of 2.8 W/m 2 K, while the frame has a U-value of 2.27 W/m 2 K and it consists of the 15% of the opening area. The equivalent U-value of the total window was calculated at 2.72 W/m 2 K, the g-value of the glazing was selected at 0.755, while no other shadings were added.

Mathematical Formulation
Basic mathematical equations concerning the present simulation are given in this subsection in order to properly explain the methodology followed. These equations are used in order to simulate properly the examined biomass-based heating/cooling system and the scope of their presentation is a clear description of the followed methodology.

General Equations for the Building and the Boiler
The thermal transmittance of the structural elements (U) is calculated as below [27]: where the internal heat convection coefficient (hin = 7.69 W/m 2 K), the external heat convection coefficient (hout = 25 W/m 2 K), (Li) is the thickness of every material, (ki) is the thermal transmittance of the materials, and (N) is the number of the materials of the examined structural element.
For the window, the (Uwind) is calculated as below [25]: where (A) is the area, "glaz" is the subscript for the glazing, and "fram" is the subscript for the frame of the window. In this work, the frame has 15% of the total area of the window.
The boiler efficiency (ηboiler) can be written as below: where (Qu) is the useful heat production and the biomass heat input is symbolized with (Qb). In the present study, the boiler efficiency is assumed to be 90% [11]. The consumption of the biomass (mb) can be found by using the higher heating value (HHV) as below: In the present work, the boiler feeds the generator and so it can be said: In the cooling mode, the cooling is calculated as below: In the heating mode with the heat pump operation (Tam ≥ 9 °C ), the heating is given as: In the heating mode without the heat pump operation, but with the boiler-only operation (Tam < 9 °C ), the heating is given as:

Equations for the Absorption Heat Pump
The operating working fluid in the single-effect absorption heat pump is the LiBr/water and Figure 1 clearly shows this device, while Figure 2 depicts the different operating modes. The present mathematical modeling is based on some usual assumptions that are made in similar simulations. More specifically [28]: − Every device is assumed to be in steady-state conditions in order to apply the energy balance. − There are no piping pressure losses, and the pressure level changes only in the pumps and in the valves. − The process in the throttling valves is ideal and so the enthalpy is preserved.

−
There is no LiBr in the condenser and in the evaporator devices. − There are no thermal losses from the system to the ambient. − The exits of the devices (evaporator, condenser, absorber and generator) are assumed to be saturated state points.
Below, the basic mathematical formulation, based on the aforementioned assumptions, is given: The energy balance in the generator is given as: The energy balance in the absorber is given as: The energy balance in the evaporator is given as: The energy balance in the condenser is given as: The expansions in the throttling valves are assumed to be isenthalpic processes due to the lack of thermal losses. Thus, it can be written as: The work consumption in the solution pump is very small and thus it can be written as: The energy balance in the heat exchanger can be written as follows: The heat exchanger effectiveness of the solution heat exchanger is defined as follows: The total mass flow rate balance in the generator is calculated as follows: The LiBr mass flow rate in the generator is calculated as follows:

Evaluation Indexes
The cooling coefficient performance (COPcool) of the absorption heat pump is given as: The heating coefficient performance (COPheat) of the absorption heat pump is given as: COP heat = Q con +Q abs Q gen (21) The overall system cooling performance index (ηen,cool) is defined as: The overall system heating performance index (ηen,heat) is defined as: For T am ≥ 9 °C: η en,heat = For T am < 9 °C: η en,heat = The exergy efficiency of the unit in the cooling mode (ηex,cool) is given as: The exergy efficiency of the unit in the heating mode (ηex,heat) is given as: The cooling temperature is set at (Tcool = 278.15 K) and the heating temperature at (Theat = 323.15 K). The exergy flow of the biomass rate was selected to be the same as the heat input of the biomass fuel [29].

Followed Methodology
In the present study, the examined building was studied with the TRNSYS software [21] and the heating/cooling loads were calculated during the whole year. The examined location is the city of Athens in Greece (37°59′ N, 23°43′ E) which is an ideal example for simulating a building with significant heating and cooling loads due to the existing climate conditions. Figure 3 shows the ambient temperature distribution in Athens during the year. It is obvious that during the summer there are high-temperature values, while the winter presents relatively low temperatures.  The heat pump model was developed in the Engineering Equation Solver tool [22] which gives the possibility for solving numerous equations together and exploiting the data from the existing libraries for the thermodynamic properties of the working fluids. The final results of the present study were found by inserting the thermal loads of the TRNSYS software into the Engineering Equation Solver tool in order to estimate the energy demands in biomass by assuming 100% satisfaction of the building's thermal loads.
At this point is useful to state that the present modeling of the absorption heat pump is practically the same as for a single-effect absorption chiller operating with LiBr/water. This modeling has also been used in previous studies of the present research team where the respective validation is given. For example, Ref. [30] shows that the deviation of the present model with the data of Ref. [28] is lower than 1%, something that indicates the high validity of the used methodology. Thus, there is not any reason to repeat the validation of the model in the present work. Regarding the TRNSYS software, it is a commercial tool that has been used in numerous studies in the past and it is a validated one according to its company [31]. Additionally, this tool has been used by the authors in the past in various studies, for cooling purposes [32,33], and for heating purposes [34,35], while a verification with the eQUEST software is given in Ref. [36].

Thermal Loads of the Examined Buildings
The first step is the presentation of the thermal loads of the examined building. Figure 4 depicts the instantaneous thermal loads for heating and cooling. It is obvious that in the majority period of the year, there is a need for providing heating or cooling. This result is a reasonable one that is based on the climate conditions of the examined location. The maximum heating load was calculated at 8.35 kW, while the maximum cooling was at 8.65 kW. Figure 5 shows the cumulative values of the thermal loads in terms of [kWh]. Specifically, the yearly demand for heating was calculated at 9136 kWh, while the cooling was at 8168 kWh. The respective specific values were found to be 22.84 kWh/m 2 for heating and 20.42 kWh/m 2 for cooling. The specific value is relatively low, something that is justified by the use of insulation in the building.   Figure 6 and Table 4 show the monthly variation of the monthly thermal loads expressed in [kWh]. Heating demand exists from November up to May, but it is too low in May which can be assumed approximately as zero. The maximum heating need is found in January with 2712 kWh with February to follow with 2193 kWh and the next following value is 1956 kWh for December. Regarding the cooling needs, they appear from May up to October, while they are very low in October. The maximum cooling needs were calculated in July with 2751 kWh, while August follows by 2669 kWh and June with 1412 kWh.

System Results
The performance of the boiler-driven absorption heat pump is given in Figure 7. This figure shows its performance for the different operating conditions and more specifically there are three operating modes, two for heating and one for cooling. In this simplistic presentation, it is assumed that the cooling mode starts for ambient temperatures over 20 °C and the cooling performance has a decreasing rate with the ambient temperature increase. The range of the cooling performance is from 0.62 to 0.78. It is significant to state that the increase in the ambient temperature over 35 °C leads to an important reduction in the performance because the system faces difficulties to produce cooling on very hot days. The heating modes are separated by the critical ambient temperature of 9 °C . In lower ambient temperatures, only the boiler operates with an efficiency of 90%, while in higher temperatures the heat pump operates. The performance of the heat pump ranges from 1.49 to 1.63 and it has an increasing trend with the ambient temperature increase.  Figure 8 shows the instantaneous demand for biomass heat and Figure 9 shows the respective instantaneous demand in fuel mass flow rate. There is a great variation in the demand due to (i) the variation of the load, (ii) the variation of the efficiency of the absorption heat pump and ii) the different operating modes of the system. The maximum demand for the heating mode was found to be 9.28 kW with a respective fuel mass flow rate of around 0.5 g/s. In the cooling mode, the maximum demand is 12.83 kW with a respective fuel demand of around 0.7 g/s. The difference between the maximum demands is justified by the reduced efficiency of the system during the cooling period.

Heating period Cooling period
Figures 9 and 10 depict the cumulative energy demand curves for the fuel energy and the fuel mass consumption. In total, the heat input demand from the biomass was calculated at 17859 kWh and it is separated into two quantities one for heating and one for cooling. The fuel energy demand for cooling was calculated at 10870 kWh (61% of the total demand), while the heating demand is 6989 kWh (39% of the total demand), significantly lower due to the higher heating performance. It is useful to state that the respective mass of the fuel, according to Figure 11, for the heating production is 1644 kg, while for cooling is significantly higher at 2119 kg, while the total consumption during the year is 3763 kg. More specifically, Figure 12 makes clear that the efficiency during the winter is higher compared to the summer period. It was calculated that the mean heating energy performance is 1.307, while the respective cooling efficiency is 0.751.

Time (h)
Heating period Cooling period Total Figure 11. The cumulative biomass consumption for heating, cooling, and total.  Figure 13 shows that the heating performance is characterized by higher exergy efficiency values compared to the cooling period. It is important to state that the mean heating mode exergy efficiency is 14.48%, while for the cooling mode is only 6.31%, respectively. The main reason for this result is the reduced energy performance during the summer which leads to higher fuel consumption in this period.  Table 5 includes all the aforementioned data. Moreover, this table indicates that the cooling period is 2627 h and the heating period is 3041 h. The heating period is separated into two sub-periods, one of 1126 h where the system operates with only the boiler (very cold conditions with an ambient temperature lower than 9 °C ) and 1915 h where the heat pump operates. The period the absorption heat pump operates during the winter is 63% of the time that there is a heat demand. At this point, it is critical to comment that the aforementioned percentage is dependent on the climate conditions of the examined location and it is a critical parameter that determines the mean efficiency of the heating mode. Also, it is important to state that it influences the mean heating performance in a great way. Heating period Cooling period Figure 12 shows that the system operates during the heating mode with an efficiency of around 1.55 and 0.90 depending on the ambient temperature. So, the mean value of around 1.31 is found after the integration during the winter period and a critical parameter in the final result is the percentage of the heating mode operation with the absorption heat pump. A simplistic mathematical analysis would state that [63%·1.55 + 37%·0.9 ≈ 1.31]. This calculation validates the analysis that has been presented regarding the variation of the modes during the heating period. At this point, it is useful to separate the heating input in the system in order to determine the contribution of the ambient to heat production. The yearly heating load is calculated at 9136 kWh and the useful heat from biomass that inserts into the system (in the heat pump or in the building directly) is found to be 6290 kWh. Consequently, the heat input from the ambient to the evaporator of the heat pump is around 2846 kWh. Therefore, it can be said that the ambient gives the 31.1% of the heating demand, which is a significant amount that indicates the existence of fuel savings.
Another interesting conclusion is that the heating mode is always more efficient than the cooling mode, as Figure 7 indicates. Thus, it is critical to state that the present system is an ideal choice for locations with high both heating and cooling loads. Moreover, this system has to be avoided in very cold climates because in these locations the heating heat pump does not operate easily, and also there are no cooling loads. Additionally, in extremely hot climates, this system is not an optimum choice because there is a lack of heating demand, and also the cooling demand increases due to the high ambient temperatures, the fact that reduces the performance. To conclude, climate conditions such as Athens (Greece), which is a typically Mediterranean city, are ideal for applying the present configuration.

Comparison with a Conventional System
At this point, it is important to compare the present system with another conventional one in order to investigate the sustainability of the suggested unit. The conventional system includes an oil boiler of 85% and an electric heat pump for cooling with COPcool = 2.8, while the lower heating value of the oil was selected at the typical value of 42000 kJ/kg [25]. The cost of the oil was estimated at 2 €/Lt, according to the recent economic trends, the electricity price at 0.2 €/kWh, the biomass cost at 0.65 €/kg, and the oil density selected at 0.86 kg/Lt. The primary energy factors for oil, electricity, and biomass were selected at 1.1, 2.9, and 1.0, respectively [25]. The specific CO2 emissions coefficients were selected at 0.264 kgCO2/kWhth for oil [25] and 0.487 kgCO2/kWhel for electricity [37].
According to the results, the yearly operating cost of the present system is 2446 €, while the conventional system leads to 2726 €. Therefore, there is a decrease of 10.3% in the yearly operating costs. The primary energy consumption of the present system is 19303 kWh per year, while the conventional system consumes 20284 kWh. Therefore, there are 4.8% savings in primary energy during the year. Additionally, the yearly CO2 avoidance with the exploitation of biomass in the present systems was found to be 4258 kgCO2 per year. The aforementioned values show that the suggested system has significant advantages in terms of economic gain, primary energy, and CO2 emissions avoidance. Figure 14 depicts the cumulative primary energy demand during the year for the present system and the conventional one. It is clear that the absorption heat pump is significantly more efficient during the heating period, while it is less efficient during the winter period. However, the final primary energy demand is lower for the absorption heat pump and thus it is the preferable system energetically. Moreover, Figure 15 illustrates the cumulative operating costs of the two examined cases. The absorption heat pump is again a more viable choice during the winter, while during the summer the conventional electric heat pump is the best choice. However, the overall year period makes clear that the absorption heat pump leads to lower operating costs and thus it is the optimum choice financially.  At this point, some limitations of the suggested system have to be discussed. Biomass is an environmentally friendly fuel sometimes its use is associated with a smell which makes some users avoid it. However, its use is safe and there is not any danger. Moreover, the biomass boilers in the households need a relatively often cleaning from the ash, something that has to be taken into consideration. Additionally, another aspect that has to be taken into account is the availability of biomass, when it is needed in huge quantities, something that can be a possible issue if its utilization is dramatically expanded.

Conclusions
Biomass exploitation for covering the thermal loads of the buildings, both for heating and cooling, consists of a sustainable solution with multiple benefits. The present study suggested a novel system based on an absorption heat pump operating with the LiBr/water working pair. As a reference, a case study of a building with a 400 m 2 floor area is examined in the location of Athens, Greece. The examined building presents specific heating and cooling loads at 22.84 kWh/m 2 and 20.42 kWh/m 2 , respectively. The maximum heating loads were found in January and the maximum cooling load in July.
According to the simulation results, the total biomass energy demand was found to be a total of 17859 kWh and it is separated at 10870 kWh for cooling (61% of the total demand) and 6989 kWh for heating (39% of the total demand). Moreover, the yearly energy performance for heating is found to be 1.307, while the yearly energy performance for cooling is calculated at 0.751. The exergy efficiency of the heating period was found to be 14.48% and for the cooling period at 6.31%. Additionally, it is interesting to state that the cooling period of the system is 2627 h and the heating period is 3041 h. The heating period is separated into two parts; 63% with heat pump operation and 37% with the only boiler operation. It is important to highlight that the present system was found to be more efficient compared to another conventional system which operates with an oil boiler and electric heat pump. More specifically, a 10.3% decrease in the yearly operating cost was calculated, as well as 4.8% savings in primary energy were found. The CO2 avoidance was estimated at 4258 kgCO2 per year.
In the future, the examined system can be tested in different climate conditions, and also it can be evaluated financially in detail. Additionally, different working pairs can be tested and for example, the use of water/NH3 is a choice that has to be investigated. Lastly, the hybridization of the system with solar thermal systems, as well as the design of the system for domestic hot water production are extra ideas that have the potential for sustainable designs.
Author Contributions: E.B.: conceptualization, methodology, software, validation, formal analysis, investigation, writing-original draft preparation, writing-review and editing; P.L.: software, writing-original draft preparation; C.T.: writing-original draft preparation, supervision. All authors have read and agreed to the published version of the manuscript.

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