Performance Assessment of a Hybrid System with Hydrogen Storage and Fuel Cell for Cogeneration in Buildings

: The search for new fuels to supersede fossil fuels has been intensiﬁed these recent decades. Among these fuels, hydrogen has attracted much interest due to its advantages, mainly cleanliness and availability. It can be produced from various raw materials (e.g., water, biomass) using many resources, mainly water electrolysis and natural gas reforming. However, water electrolysis combined with renewable energy sources is the cleanest way to produce hydrogen while reducing greenhouse gases. Besides, hydrogen can be used by fuel cells for producing both electrical and thermal energy. The aim of this work was towards e ﬃ cient integration of this system into energy e ﬃ cient buildings. The system is comprised of a photovoltaic system, hydrogen electrolyzer, and proton exchange membrane fuel cell operating as a cogeneration system to provide the building with both electricity and thermal energy. The system’s modeling, simulations, and experimentations were ﬁrst conducted over a short-run period to assess the system’s performance. Reported results show the models’ accuracy in analyzing the system’s performance. We then used the developed models for long-run testing of the hybrid system. Accordingly, the system’s electrical e ﬃ ciency was almost 32%. Its overall e ﬃ ciency reached 64.5% when taking into account both produced electricity and thermal energy.


Introduction
The deployment of renewable energy sources (RES) as green energy producers, such as solar energy, wind, geothermal, and biomass, as well as the search for new sources of fuel (e.g., hydrogen and biogas) have been intensified these last decades in order to supersede fossil products (e.g., coal, oil, natural gas) [1][2][3][4][5]. The latter currently produce more than 70% of the global electricity needs [6]. Therefore, Morocco as a sunny country, with an average daily global radiation of about 5 kWh/m 2 [7], is putting more emphasis on solar energy. It aims to increase the proportion of RES from 42% in 2020 to 52% in 2030 [8] in order to meet the increase of electricity consumption. This will reduce its dependence on fossil products, which is currently around 90% [6]. Thus, many research and development (R&D) institutions (e.g., IRESEN-Institut de Recherche en Energie Solaire et Energies Nouvelles, MASEN-Moroccan Agency for Solar Energy) have been built, and many RES plants (e.g., Noor in Ouarzazate) have been deployed [9]. Furthermore, this North African country was involved, through MASEN, in the DESERTEC project, which had the intention of deploying solar their state (e.g., steady state or dynamic conditions) and spatial dimensions (e.g., 0-D, 1-D, 2-D, or 3-D): (i) analytical models that are based on approximations and assumptions; (ii) semi-empirical models, which use experiments to validate theoretical models; and (iii) mechanistic models that are based on physics and electrochemistry laws. Concerning hydrogen storage, numerous models have been developed in order to study the behavior of its different types [44,48,50,51]. In addition, numerous analytical models [52], mechanistic models [49,[52][53][54], and empirical models [48,52,[54][55][56][57] have been developed for the purpose of assessing the performance and the behavior of the different types of fuel cells and estimating and predicting their current/voltage relationship.
The work presented in this article is a part of two projects. The first project, named PROPRE.MA (PROductivité Photovoltaïque à l'échelle REgionale dans tout le Maroc), is funded by IRESEN and managed by the Faculty of Sciences Semlalia Marrakech (FSSM) [58][59][60]. It aims at mapping the productivity and the efficiency of three different technologies of PV panels (monocrystalline silicon, polycrystalline silicon, and amorphous silicon) in 21 different locations in Morocco. The second project, named MIGRID, is an ongoing project aiming at deploying a micro-grid (MG) system comprised of PV panels in order to generate electricity from solar irradiation; wind turbines to generate electricity from the wind; batteries for storing the overproduction of electricity in order to use it when there is no production; converters (e.g., DC/DC converters, AC/DC converters, and DC/AC converters); a hydrogen electrolyzer to produce hydrogen using the overproduction of our RES systems; hydrogen tanks to store the produced hydrogen; a fuel cell that uses the stored hydrogen and the air's oxygen for supplying the building with both valuable electricity and thermal energy when needed; and finally, the electric grid in order to feed electricity to the building when there is no production, i.e., when there is no hydrogen left and the batteries are empty. Mainly, the main aim of the MIGRID project is to test some approaches (e.g., dimensioning and control strategies) and validate the components' models using experiments [39,61,62]. Therefore, the hybrid system that we have considered in this article is comprised of a PV system, an electrolyzer for hydrogen production, and a PEM fuel cell that is principally operating as a CHP system in order to provide the building with both valuable electrical and thermal energy.
The aim of the work presented in this paper is towards efficient integration of a hybrid system into buildings together with its performance assessment. Thus, the main contributions in this article are twofold: (i) modeling, simulating, and experimenting the hybrid system and showing the accuracy of the developed models, and (ii) using the validated models for studying the behavior of the hybrid system over one year and assessing its performance when operating as a CHP system for supplying the building with both electrical and thermal energy.
The remainder of this article is structured as follows. The proposed hybrid system is described with more detail in Section 2. Section 3 introduces the modeling of the MG system's components. The obtained simulation and experimental results are presented in Section 4. Finally, the main conclusions and perspectives are given in Section 5. Figure 1 depicts the architecture of the proposed hybrid system. It is comprised of PV panels, converters, a hydrogen electrolyzer, hydrogen tanks, and a fuel cell. Actually, as a part of the PROPRE.MA project, the PV system used in this experiment is set up in the Faculty of Sciences El Jadida (latitude 33.22 • N, longitude 8.48 • W, altitude 24 m, and 2.1 km away from the Atlantic Ocean) in Morocco, more precisely on the roof of the physics department building. It is comprised of PV panels in order to generate electricity from solar irradiation and an inverter in order to convert the direct current (DC) generated by the PV panels into an alternative current (AC) [63]. Concerning the PV panels, we installed 8 polycrystalline silicon PV modules, which were manufactured by Solar World, type Sun module plus. Each module had a rated power of about 255 W, an open circuit voltage (V OC ) equal to 37.6 V, and a short circuit current (Isc) of about 8.81 A, as presented in Table 1. In this experiment, we chose to work with the polycrystalline silicon PV modules due to their cost-performance ratio  Regarding the inverter that is manufactured by SMA, type Sunny boy 2000HF-30, it has already been reported with more detail by the authors in [59]. It is noteworthy that we chose the city of El Jadida because of its weather conditions (e.g., moderate temperature), and also the PV system was already set up in this city in the frame of the project. Concerning the data acquisition, an acquisition platform comprised of many sensors and electronic cards was installed in order to obtain the PV system's data as well as the weather data. Regarding the PV system's data, it was measured using current and voltage sensors, which were integrated into the inverter. The chosen sampling time in this experiment was 5 min, and the data could be accessed via a Bluetooth modem or transferred by means of an Ethernet or Internet network [60,64]. As for the weather data, a weather station was deployed near the PV panels in order to collect some meteorological data, such as wind speed and direction, solar irradiation, and the ambient and PV module's temperature. Mainly, it was comprised of several sensors: i) a temperature sensor for measuring the ambient temperature; ii) a Pt100 temperature sensor for capturing the PV module's temperature; iii) an anemometer for determining the wind speed; iv) a wind vane in order to measure the wind direction; v) two small PVs (one tilted by 30° and the other horizontal), with a rated power of about 20 W each, in order to measure the inplane and horizontal solar irradiances, which are proportional to their short circuit current; and vi) a PC-duino and Raspberry pi in order to gather, process, and store the measured data [63,64]. Furthermore, an electrolyzer was used to produce hydrogen, which was stored and then used when needed. The used PEM electrolyzer, which had an operative power consumption of about 800 W and  Regarding the inverter that is manufactured by SMA, type Sunny boy 2000HF-30, it has already been reported with more detail by the authors in [59]. It is noteworthy that we chose the city of El Jadida because of its weather conditions (e.g., moderate temperature), and also the PV system was already set up in this city in the frame of the project. Concerning the data acquisition, an acquisition platform comprised of many sensors and electronic cards was installed in order to obtain the PV system's data as well as the weather data. Regarding the PV system's data, it was measured using current and voltage sensors, which were integrated into the inverter. The chosen sampling time in this experiment was 5 min, and the data could be accessed via a Bluetooth modem or transferred by means of an Ethernet or Internet network [60,64]. As for the weather data, a weather station was deployed near the PV panels in order to collect some meteorological data, such as wind speed and direction, solar irradiation, and the ambient and PV module's temperature. Mainly, it was comprised of several sensors: (i) a temperature sensor for measuring the ambient temperature; (ii) a Pt100 temperature sensor for capturing the PV module's temperature; (iii) an anemometer for determining the wind speed; (iv) a wind vane in order to measure the wind direction; (v) two small PVs (one tilted by 30 • and the other horizontal), with a rated power of about 20 W each, in order to measure the in-plane and horizontal solar irradiances, which are proportional to their short circuit current; and (vi) a PC-duino and Raspberry pi in order to gather, process, and store the measured data [63,64]. Furthermore, an electrolyzer was used to produce hydrogen, which was stored and then used when needed. The used PEM electrolyzer, which had an operative power consumption of about 800 W and produced hydrogen with a purity of about 99.999% (see Table 2 for more characteristics), is manufactured by Heliocentris, type HG 100. It was used in order to produce hydrogen using deionized water and electricity. The produced hydrogen was stored in hydrogen tanks so that it could be used later by the fuel cell for producing both valuable electrical and thermal energy according to the building's need. The PEM fuel cell used in this experiment, which had a rated power of about 1.2 kW, an output voltage in the range of 20-36 VDC, and a rated current of around 65 ADC (see Table 3 for more features), is manufactured by Heliocentris, type Nexa 1200. Concerning the DC produced electricity, an inverter was used in order to convert it into an AC for feeding the building's AC appliances. Besides, several current and voltage sensors were deployed in order to measure the power in its different branches. In addition, the hydrogen flow and temperature of the fuel cell's stack and hydrogen tanks were measured using an integrated flow meter and temperature sensors. In fact, the fuel cell was purchased together with these sensors, which were fully embedded into its acquisition system. The measured data were then stored and visualized using a software tool, which was provided by the same company that manufactured the fuel cell. Regarding the produced thermal energy (Q), it could be used for either heating the building in the winter time using an air to air heat exchanger or cooling it during the summer period using a refrigeration process. It is computed using Equation (1): where ρ air is the density of the air (kg/m 3 ), V air being the volume of the heated air (m 3 ), C p,air represents the specific heat capacity of the air (kJ/kg.K), and T f inal and T initial are the final and initial temperatures of the air, respectively (K) [65,66]. In order to capture most of the produced thermal energy by the fuel cell, we deployed a box at the fuel cell's heat exhaust. Firstly, the used fuel cell system, which is presented in Figure 2, is constituted of hydrogen tanks, a 1.2 kW PEM fuel cell (see Table 3), converters (DC/DC converter and inverter), batteries to provide the computer used for data visualization (e.g., fuel cell's current and voltage, hydrogen flow rate, stack temperature) with the needed electricity, and a DC load.

Hybrid System's Description
Sustainability 2020, 12, x FOR PEER REVIEW 6 of 20 visualization (e.g., fuel cell's current and voltage, hydrogen flow rate, stack temperature) with the needed electricity, and a DC load. Regarding the deployed box, the volume of the air trapped into it was approximately equal to 0.0478 m 3 . Besides, an IoT/Big Data platform, which consisted of three temperature sensors (DS18B20), an Arduino Nano board, and Raspberry pi, was developed and deployed for the aim of measuring the temperature of the air trapped in the box [61]. Then, the collected data was sent by means of the Raspberry pi to a cluster, where it was processed, stored, and visualized. It could also be visualized and extracted using a web application that we had already developed. Finally, the temperature data was used in order to compute the amount of produced thermal energy using Equation (1).

Figure 2.
The deployed fuel cell system as well as the box for capturing the released heat by the fuel cell.

System's Modeling
In this section, the system's components (e.g., PV panels, converters, hydrogen electrolyzer, hydrogen tanks, and fuel cell) have been modeled and then connected in order to build the system's model. Then, the simulation was conducted under MATLAB using similar experimental conditions (e.g., solar irradiation, ambient temperature). Figure 3 presents the flowchart diagram of the system's model. It shows mainly the different phases followed for computing the components' production, i.e., power, hydrogen volume, and thermal energy. Besides, the remainder of this section is devoted to the description and development of the models of the system's components. Regarding the deployed box, the volume of the air trapped into it was approximately equal to 0.0478 m 3 . Besides, an IoT/Big Data platform, which consisted of three temperature sensors (DS18B20), an Arduino Nano board, and Raspberry pi, was developed and deployed for the aim of measuring the temperature of the air trapped in the box [61]. Then, the collected data was sent by means of the Raspberry pi to a cluster, where it was processed, stored, and visualized. It could also be visualized and extracted using a web application that we had already developed. Finally, the temperature data was used in order to compute the amount of produced thermal energy using Equation (1).

System's Modeling
In this section, the system's components (e.g., PV panels, converters, hydrogen electrolyzer, hydrogen tanks, and fuel cell) have been modeled and then connected in order to build the system's model. Then, the simulation was conducted under MATLAB using similar experimental conditions (e.g., solar irradiation, ambient temperature). Figure 3 presents the flowchart diagram of the system's model. It shows mainly the different phases followed for computing the components' production, i.e., power, hydrogen volume, and thermal energy. Besides, the remainder of this section is devoted to the description and development of the models of the system's components.

PV System's Model
PV solar panels have been extensively used in order to produce electricity from solar irradiation. Many PV technologies have been developed these last decades, such as crystalline silicon (c-Si) solar cells, thin film cells, and multi-junction cells. In fact, the most used and manufactured ones are the crystalline silicon (c-Si) solar cells. These latter are classified into three categories, namely, monocrystalline silicon (mc-Si), polycrystalline silicon (pc-Si), and amorphous silicon (a-Si). The polycrystalline silicon PV modules are the most used due to their cost-performance ratio compared to the monocrystalline and amorphous ones. Therefore, their modeling and performance assessment is extremely needed in order to evaluate their behavior and performance under various conditions (e.g., solar irradiation, temperature). Thus, based on the Gow and Manning model [35,36,39] and the Shockley diode equation, a mathematical model of the solar PV cells' equivalent electrical circuit has been developed. This model consists of a current source parallel to a diode, a series resistance , and a shunt resistance ℎ , as shown in Figure 4. From the electrical circuit of the PV cell, which is illustrated in Figure 4, the equation of the PV cell's -characteristic is expressed using Equation (2): where ℎ is the photocurrent, being the diode current, and ℎ is the shunt current [35,36,39]. Using the Shockley diode equation as well as the expressions of the photocurrent and the shunt current presented in our previous work [39], the -characteristic of the PV module can be given by Equation (3):

PV System's Model
PV solar panels have been extensively used in order to produce electricity from solar irradiation. Many PV technologies have been developed these last decades, such as crystalline silicon (c-Si) solar cells, thin film cells, and multi-junction cells. In fact, the most used and manufactured ones are the crystalline silicon (c-Si) solar cells. These latter are classified into three categories, namely, monocrystalline silicon (mc-Si), polycrystalline silicon (pc-Si), and amorphous silicon (a-Si). The polycrystalline silicon PV modules are the most used due to their cost-performance ratio compared to the monocrystalline and amorphous ones. Therefore, their modeling and performance assessment is extremely needed in order to evaluate their behavior and performance under various conditions (e.g., solar irradiation, temperature). Thus, based on the Gow and Manning model [35,36,39] and the Shockley diode equation, a mathematical model of the solar PV cells' equivalent electrical circuit has been developed. This model consists of a current source parallel to a diode, a series resistance R s , and a shunt resistance R sh , as shown in Figure 4.

PV System's Model
PV solar panels have been extensively used in order to produce electricity from solar irradiation. Many PV technologies have been developed these last decades, such as crystalline silicon (c-Si) solar cells, thin film cells, and multi-junction cells. In fact, the most used and manufactured ones are the crystalline silicon (c-Si) solar cells. These latter are classified into three categories, namely, monocrystalline silicon (mc-Si), polycrystalline silicon (pc-Si), and amorphous silicon (a-Si). The polycrystalline silicon PV modules are the most used due to their cost-performance ratio compared to the monocrystalline and amorphous ones. Therefore, their modeling and performance assessment is extremely needed in order to evaluate their behavior and performance under various conditions (e.g., solar irradiation, temperature). Thus, based on the Gow and Manning model [35,36,39] and the Shockley diode equation, a mathematical model of the solar PV cells' equivalent electrical circuit has been developed. This model consists of a current source parallel to a diode, a series resistance , and a shunt resistance ℎ , as shown in Figure 4. From the electrical circuit of the PV cell, which is illustrated in Figure 4, the equation of the PV cell's -characteristic is expressed using Equation (2): where ℎ is the photocurrent, being the diode current, and ℎ is the shunt current [35,36,39]. Using the Shockley diode equation as well as the expressions of the photocurrent and the shunt current presented in our previous work [39], the -characteristic of the PV module can be given by Equation (3): From the electrical circuit of the PV cell, which is illustrated in Figure 4, the equation of the PV cell's I PV -V PV characteristic is expressed using Equation (2): where I ph is the photocurrent, I d being the diode current, and I sh is the shunt current [35,36,39].
Using the Shockley diode equation as well as the expressions of the photocurrent and the shunt current presented in our previous work [39], the I PV -V PV characteristic of the PV module can be given by Equation (3): where N p is the number of parallel cells in a module, I s being the diode saturation current (A), q stands for the electron charge (1.6 × 10 −19 C), N s represents the number of series cells in a module, n is the diode ideality factor, K being the Boltzmann constant (1.38 × 10 −23 J/K), and T represents the module's temperature (K). Regarding the converter (e.g., DC/DC converter and DC/AC converter) models, they have been reported and described with more detail in our previous work, presented in [39]. They were used in order to connect the different components, convert and adapt the variable voltages of the different energy sources to a constant one that could be consumed by the loads, and control the energy flow provided by the different energy sources based on the loads' consumption.

Hydrogen Electrolyzer Model
Hydrogen can be produced, as stated above, from various raw materials (e.g., water, biomass, fossil fuels) using many resources, mainly water electrolysis and natural gas reforming [7,[22][23][24]. Since the main purpose of using hydrogen is to supersede fossil fuels, water electrolysis, especially when using the electricity produced by RES (e.g., solar energy), remains the cleanest method for hydrogen generation from non-fossil fuels while reducing greenhouse gases due to its multiple advantages, such as cleanliness, availability of water, and the high purity of produced hydrogen [24]. Actually, based on the electrolyte and the operating temperature, hydrogen electrolyzers are categorized into three types, namely, the alkaline electrolyzer, the PEM electrolyzer, and the solid oxide electrolyzer [24,26,28]. Besides, the PEM electrolyzer is the most used, especially in buildings, due to its various benefits compared to the other electrolyzers, such as a compact system design, rapid system response, high dynamic operation, high efficiency, and high purity of produced hydrogen [26].
The PEM electrolyzer is constituted of a cathode, an anode, and a solid polymer membrane (electrolyte). It is an electrochemical device that is used for producing both hydrogen and oxygen from electricity and deionized water. In fact, deionized water is inserted into the anode electrode, where it is separated into a hydrogen proton (H + ) and oxygen (O 2 ). The hydrogen proton is then transported to the cathode electrode through the electrolyte. In the cathode, where a DC current is applied, the hydrogen proton consumes electrons for producing hydrogen gas (H 2 ) [24]. Equations (4) and (5) present the reactions that take place at the PEM electrolyzer's electrodes (anode and cathode, respectively).
Regarding the modeling of hydrogen electrolyzers, several models, which are based especially on their polarization curve as well as the volume of the produced hydrogen, have been proposed in order to simulate their behavior under different conditions and compute the produced hydrogen's volume [45]. The assessment of their performances is also based on the evaluation of the temperature and the pressure of their electrodes and electrolytes. However, these two parameters could be neglected for the sake of simplicity. In fact, the authors in [45] reviewed four models used for estimating the electrolyzers' behavior as well as computing the produced hydrogen's volume, while indicating the degree of complexity of each model. Then, the four models were compared to an electrolyzer's real behavior. The aim was to validate and assess the accuracy of these models. Finally, the authors in [45] found that the responses of the four models were too close to the real response of the device, with slight errors. For all these reasons, and for the sake of simplicity, the simplest model, which depends on the applied current, was used in order to estimate the behavior of the electrolyzer and the produced Sustainability 2020, 12, 4832 9 of 21 hydrogen's volume. This equivalent circuit model is comprised of a DC voltage source in series with a resistor, which symbolizes the internal electrical losses [45]. Actually, the produced hydrogen's volume (V H 2 in m 3 ) is calculated using Equation (6): where n c is the number of electrolyzer cells, z stands for the number of electrons (2 for hydrogen and 4 for oxygen), R is the gas constant (8.3145 J/mol·K), P represents the working pressure of the electrolyzer (Pa), I ez being the applied current to one of the electrolyzer's cells (A), F stands for the Faraday constant (96485.33 C/mol), T represents the operating temperature of the electrolyzer (K), and finally t is the operating time of the electrolyzer (second).
The produced hydrogen has to be stored and then used when it is needed by the building. In fact, since hydrogen has a small density at ambient conditions (1 kg of hydrogen is equivalent to 11 m 3 [14]), its storage necessitates the increase of its density and thus the diminution of its volume. Actually, based on its application, storage size, temperature, and pressure, different methods of reversibly storing hydrogen exist: (i) its compression into high pressure gas cylinders, (ii) its storage as a liquid, and (iii) its storage in a solid state in the metal and non-metal hydrides [14]. For instance, the high pressure gas cylinders are the most used systems for storing hydrogen, with a maximum pressure of 20 MPa [14]. This is attributed to the fact that they are currently the best known and understood methods [67]. In this case, the hydrogen is in its gaseous state, with a density of 0.089886 kg/m 3 at 0 • C and 1 bar [14]. On another note, hydrogen gas, which is not an ideal gas, could be described by the Van der Waals equation (Equation (7)): where P represents the gas pressure (Pa), n being the number of moles (mol), R stands for the gas constant that is equal to 8.314 J/K·mol, T is the gas absolute temperature (K), V represents the gas volume (m 3 ), b being the occupied volume by hydrogen molecules that is equal to 2.661 × 10 −5 m 3 /mol, and finally a stands for the dipole interaction or repulsion constant, which is equal to 2.476 × 10 −2 m 6 ·Pa/mol 2 [14]. As for liquid hydrogen, it can be stored in cryogenic tanks at −253 • C, which is found in a small zone between its triple point and its critical point, and ambient pressure [14]. Since liquefying hydrogen implies an increase of its volumetric energy density (70.8 kg/m 3 ), a reduction of storage volume compared to the high pressure gas cylinders is achieved [67]. However, the liquefaction of hydrogen needs too much energy to be performed, which makes it a very expensive process. The expensive cost of this method is also caused by the fact that the cryogenic tanks have to be thermally insulated in order to avoid the boiling of hydrogen (−252.879 • C) [14]. Finally, hydrogen can be stored in its solid state, with a density equal to 70.6 kg/m 3 at a low temperature (−262 • C) [64] in metal-hydrides and non-metal hydrides. For the metal-hydrides, the metals, the intermetallic compounds, and their alloys combine with hydrogen in order to form solid metallic hydrides [14][15][16]. The metal-hydrides are considered to be an excellent alternative to hydrogen storage due to their high volumetric energy density (70.6 kg/m 3 ); however, they have a low specific energy (kWh/kg), which leads to heavier storage tanks. In addition, the high cost of metals makes their use on small-scale applications limited [67]. For these reasons, the non-metal hydrides have attracted much interest these last years due to their small weights [67]. They are a combination between some reactive non-metals, mainly, carbon, boron, nitrogen and oxygen, and hydrogen [67]. Alas, they still have the issue of not binding with hydrogen as strongly as the metal-hydrides [67].
In this work, we have chosen to store it in its gaseous state using pressurized tanks. The pressurized tanks used in this work have a recommended charging pressure of 10 bar, a maximum storage temperature of about 50 • C, and a storage capacity of 500 l. Besides, since the electrolyzer works on the same pressure as the hydrogen tanks (10 bar), there is no need for a dedicated compressor, and then V H 2 is also the volume of hydrogen stored in the pressurized tanks.

Fuel Cell's Model
Fuel cells have been extensively used these last decades in both vehicles and buildings for producing valuable electricity, thermal energy, and water from the combination of hydrogen and oxygen. This is attributed to their various benefits, mainly, cleanliness, modularity, low noise level, and higher efficiency, which could achieve 85-90% [13,22]. Actually, based on the operating temperature and the used fuel and electrolyte, fuel cells are categorized into different types, mainly, alkaline fuel cells (AFCs), phosphoric acid fuel cells (PAFCs), sulphuric acid fuel cells (SAFCs), proton exchange membrane fuel cells (PEMFCs), molten carbonate fuel cells (MCFCs), direct methanol fuel cells (DMFCs), solid polymer fuel cells (SPFCs), and solid oxide fuel cells (SOFCs) [24,30,31]. In fact, the most suitable fuel cells for micro-cogeneration applications are the low temperature PEMFCs working around 80 • C, and the high temperature SOFCs operating in the range of 600-1000 • C. However, the PEMFC, which is comprised of two platinum catalyzed porous electrodes and separated by a solid polymeric membrane electrolyte, is still the most used for micro-cogeneration applications. This is attributed to its multiple benefits, mainly, high power density, quick start-up and load following, less corrosion, and a low operating temperature [22,32].
Like the PEM electrolyzer, the PEMFC is comprised of a cathode, an anode, and a solid polymer membrane (electrolyte). It is an electrochemical device that generates electricity, thermal energy, and water from the exothermic reaction between hydrogen and oxygen. Actually, hydrogen, which is delivered to the anode, is divided into electrons (e − ) and hydrogen protons (H + ). The latter goes across the solid polymer membrane in order to achieve the cathode. In this latter, the oxygen (oxidant) has been provided, and thus it is combined with hydrogen protons in order to form water (H 2 O) [29]. The reactions that take place at the PEMFC's anode and cathode can be expressed by Equations (8) and (9), respectively: Concerning the modeling of PEMFCs, several models (i.e., analytical [52], mechanistic [52,54], and empirical [52,54]) and software tools [68] have been developed and reported in the literature for estimating and predicting their current/voltage relationship, and assessing their behavior under different circumstances. In our simulations, and for the sake of simplicity, we used the model developed and presented in [57] for the purpose of assessing the performance of our PEMFC and, therefore, estimating the amount of produced electrical energy. The used model is comprised of a controlled voltage source in series with a resistor, as presented in Figure 5. Fuel cells have been extensively used these last decades in both vehicles and buildings for producing valuable electricity, thermal energy, and water from the combination of hydrogen and oxygen. This is attributed to their various benefits, mainly, cleanliness, modularity, low noise level, and higher efficiency, which could achieve 85-90% [13,22]. Actually, based on the operating temperature and the used fuel and electrolyte, fuel cells are categorized into different types, mainly, alkaline fuel cells (AFCs), phosphoric acid fuel cells (PAFCs), sulphuric acid fuel cells (SAFCs), proton exchange membrane fuel cells (PEMFCs), molten carbonate fuel cells (MCFCs), direct methanol fuel cells (DMFCs), solid polymer fuel cells (SPFCs), and solid oxide fuel cells (SOFCs) [24,30,31]. In fact, the most suitable fuel cells for micro-cogeneration applications are the low temperature PEMFCs working around 80 °C, and the high temperature SOFCs operating in the range of 600-1000 °C. However, the PEMFC, which is comprised of two platinum catalyzed porous electrodes and separated by a solid polymeric membrane electrolyte, is still the most used for microcogeneration applications. This is attributed to its multiple benefits, mainly, high power density, quick start-up and load following, less corrosion, and a low operating temperature [22,32].
Like the PEM electrolyzer, the PEMFC is comprised of a cathode, an anode, and a solid polymer membrane (electrolyte). It is an electrochemical device that generates electricity, thermal energy, and water from the exothermic reaction between hydrogen and oxygen. Actually, hydrogen, which is delivered to the anode, is divided into electrons (e − ) and hydrogen protons (H + ). The latter goes across the solid polymer membrane in order to achieve the cathode. In this latter, the oxygen (oxidant) has been provided, and thus it is combined with hydrogen protons in order to form water (H2O) [29]. The reactions that take place at the PEMFC's anode and cathode can be expressed by Equations (8) and (9), respectively: Concerning the modeling of PEMFCs, several models (i.e., analytical [52], mechanistic [52,54], and empirical [52,54]) and software tools [68] have been developed and reported in the literature for estimating and predicting their current/voltage relationship, and assessing their behavior under different circumstances. In our simulations, and for the sake of simplicity, we used the model developed and presented in [57] for the purpose of assessing the performance of our PEMFC and, therefore, estimating the amount of produced electrical energy. The used model is comprised of a controlled voltage source in series with a resistor, as presented in Figure 5. A diode was also used in order to prevent the current flowing into the fuel cell's stack. From the electrical circuit of the fuel cell, which is presented in Figure 5, the fuel cell's stack voltage ( ) is computed using Equation (10): where is the fuel cell's stack current (A), being the internal resistance (Ω), and stands for the controlled voltage source (V). The detailed description of this model as well as the computation A diode was also used in order to prevent the current flowing into the fuel cell's stack. From the electrical circuit of the fuel cell, which is presented in Figure 5, the fuel cell's stack voltage (V FC ) is computed using Equation (10): where I FC is the fuel cell's stack current (A), R being the internal resistance (Ω), and E stands for the controlled voltage source (V). The detailed description of this model as well as the computation of the controlled voltage source (E) are presented in [57]. Besides, the amount of generated thermal energy is computed using the enthalpy change ∆H r of the reaction, which is equal to −285.84 kJ/mol of H 2 , especially when the generated water is in a liquid state. The enthalpy change of this exothermic reaction is expressed by Equation (11): where ∆G represents the Gibbs free energy, T being the stack's temperature, and ∆S stands for the entropy [65]. The reaction's enthalpy change is divided into three fractions: (i) the amount of energy responsible for the water vaporization (∆H vap = 44.01 kJ/mol), (ii) the Gibbs free energy that represents the maximum usable produced energy (i.e., electrical and thermal energy) (∆G = −237.2 kJ/mol of H 2 ), and finally (iii) the amount of generated heat (T∆S = −4.63 kJ/mol of H 2 ) [22]. However, only a small fraction of the Gibbs free energy could be used to supply the building with both electrical and thermal energy, since most of it is used for increasing the hydrogen gas temperature as well as the temperature of the fuel cell, mainly the frame, the electrodes, and the electrolyte. The amount of the energy, which is lost being used in increasing the temperature of the fuel cell as well as the hydrogen gas, is calculated using Equation (12): where m is the mass of the fuel cell's components (kg), C p stands for the specific heat capacity of the fuel cell's components (kJ/kg.K), and T f inal and T initial represent the final and the initial temperature of the fuel cell's components, respectively (K).

Simulation and Experimental Results
This section focuses on the simulations and experimentations of the deployed hybrid system. This was performed in two phases. First, technical integration and testing was carried out in order to validate the components of the system. Besides the technical integration of the components including equipment and software in preparation of the experiments, real testing scenarios were established and tested for validating the system's models over a short-run period using both simulations and experimentations. In the second phase, virtual test runs were carried out to test the performance of the hybrid system for day-to-day use, and during one year of simulations.

Results of Short-Term Simulations and Experimentations of the Hybrid System
Regarding the PV system's model, which is presented in Section 3, it was implemented into MATLAB. Then, it was simulated under the same experimental conditions, namely, the solar irradiation and the ambient temperature, in order to assess the performance of the system. In fact, an experiment was conducted over 24 h in order to collect the data to be used for simulating and validating the system's model. For instance, Figure 6a depicts the in-plane solar irradiation and the ambient temperature during 24 h. These data were measured using our deployed weather station (more details in [63,64]), and were used in order to simulate the system's model. In fact, we noticed that the maximum solar irradiation exceeded 900 W/m 2 and the ambient temperature was moderate. This is attributed to the freshness brought by the Atlantic Ocean. Besides, Figure 6b presents the simulation and the experimental results of the power generated by the PV system during 24 h. We observe that there is no production during the night since there is no solar irradiation. However, after sunrise, the PV production rises until achieving its maximum, which is almost 1.8 kW at midday, and then declines until sunset. Moreover, we observe that the simulation and experimental results match very well. Furthermore, the simulation and experimental produced energies were computed and found equal to 13.3 kWh and 13.44 kWh, respectively. Hence, the relative error between the simulation and experimental results was found equal to 1.04%. This latter shows the accuracy of the system's model and its good approximation of the real-world scenarios. Concerning the fuel cell's modeling, we used the model described in Figure 5 in order to assess its performance and determine the amount of produced electrical and thermal energy [57]. Once the model was implemented into MATLAB/Simulink, many simulations were carried out under the same experimental conditions. The stack's temperature (Figure 7a) that achieves 50 °C ; the 2 flow rate at the entrance of the fuel cell, which is approximately 5.92 nl/min; and the load's power (Figure 7b) were measured during the experiment and used as inputs to our model. Then, the simulation and experimental results were reported in order to assess the fuel cell's performance and determine the amount of produced electrical and thermal energy. Figure 7b depicts the simulation and experimental results of the fuel cell system as well as the load's power. We can obviously observe that the fuel cell system's power is always higher than the load's power in both the simulation and experiment. The difference between the experimental fuel cell system's power and the load's power, which is around 20 W, is consumed by the auxiliary components (e.g., valves, fans), which are extremely needed for the fuel cell's operation. Besides, the difference between the simulated and experimental power of the fuel cell is attributed to the fact that its electrical efficiency was taken in the simulation equal to 50%; however, this efficiency varied around 50% in the experiment according to the load's consumption. On another note, the difference between the reference electric power ( ) given by the manufacturer for a 50% efficiency (600 W) and the simulated fuel cell system's power does not achieve 10%, as depicted in Figure 7b. This shows the accuracy of the fuel cell's model.
Regarding the amount of generated thermal energy, it was estimated using the reaction's enthalpy change. In fact, we used the hydrogen flow rate, which was converted into mol/min using the hydrogen's density and molar mass, and the system's efficiency, which was approximately equal to 50%, in order to compute the total amount of generated thermal energy. This latter was found approximately equal to 33.58 KJ/min. Accordingly, the total amount of the generated heat throughout the entire experiment was 380.46 kJ. Alas, most of this energy was used for increasing the fuel cell's temperature (e.g., frame, electrodes, and electrolyte) and the hydrogen gas temperature. The amount of energy used for increasing the fuel cell's temperature was found approximately equal to 379.14 kJ. It was computed by Equation (12) using the mass of the fuel cell, the stack's temperature, and the average specific heat capacity of the fuel cell's components (0.881 KJ/Kg.K for the frame and 1.32 KJ/Kg.K for the electrodes and electrolyte). Moreover, the quantity of energy consumed for preheating hydrogen was found to be 0.0268 kJ. It was computed by Equation (12) using the hydrogen flow rate that was converted into kg/min, the specific heat capacity of hydrogen (10.140 kJ/kg.K), and Concerning the fuel cell's modeling, we used the model described in Figure 5 in order to assess its performance and determine the amount of produced electrical and thermal energy [57]. Once the model was implemented into MATLAB/Simulink, many simulations were carried out under the same experimental conditions. The stack's temperature (Figure 7a) that achieves 50 • C; the H 2 flow rate at the entrance of the fuel cell, which is approximately 5.92 nl/min; and the load's power (Figure 7b) were measured during the experiment and used as inputs to our model. Then, the simulation and experimental results were reported in order to assess the fuel cell's performance and determine the amount of produced electrical and thermal energy. Figure 7b depicts the simulation and experimental results of the fuel cell system as well as the load's power. We can obviously observe that the fuel cell system's power is always higher than the load's power in both the simulation and experiment. The difference between the experimental fuel cell system's power and the load's power, which is around 20 W, is consumed by the auxiliary components (e.g., valves, fans), which are extremely needed for the fuel cell's operation. Besides, the difference between the simulated and experimental power of the fuel cell is attributed to the fact that its electrical efficiency was taken in the simulation equal to 50%; however, this efficiency varied around 50% in the experiment according to the load's consumption. On another note, the difference between the reference electric power (P R ) given by the manufacturer for a 50% efficiency (600 W) and the simulated fuel cell system's power does not achieve 10%, as depicted in Figure 7b. This shows the accuracy of the fuel cell's model.
Regarding the amount of generated thermal energy, it was estimated using the reaction's enthalpy change. In fact, we used the hydrogen flow rate, which was converted into mol/min using the hydrogen's density and molar mass, and the system's efficiency, which was approximately equal to 50%, in order to compute the total amount of generated thermal energy. This latter was found approximately equal to 33.58 KJ/min. Accordingly, the total amount of the generated heat throughout the entire experiment was 380.46 kJ. Alas, most of this energy was used for increasing the fuel cell's temperature (e.g., frame, electrodes, and electrolyte) and the hydrogen gas temperature. The amount of energy used for increasing the fuel cell's temperature was found approximately equal to 379.14 kJ.
It was computed by Equation (12) using the mass of the fuel cell, the stack's temperature, and the average specific heat capacity of the fuel cell's components (0.881 KJ/Kg.K for the frame and 1.32 KJ/Kg.K for the electrodes and electrolyte). Moreover, the quantity of energy consumed for preheating hydrogen was found to be 0.0268 kJ. It was computed by Equation (12) using the hydrogen flow rate that was converted into kg/min, the specific heat capacity of hydrogen (10.140 kJ/kg.K), and the difference of hydrogen's temperature, which did not exceed 0.5 • C. Hence, the amount of usable energy, which could be released outside the fuel cell, was approximately 1.2930 kJ. For the aim of comparing the theoretical and experimental results and thus validating this computation methodology, the temperature of the air trapped into the box (Figure 2) was measured. Figure 8 presents the variation of the temperature of the air trapped into the box (the average value of the three sensors). Then, for computing the heat released by the fuel cell, we assumed that the air's density did not change with temperature. So, the released heat was computed by Equation (1) using the density of the air (1.2 kg/m 3 ), the volume of the air trapped in the box (0.0478 m 3 ), the specific heat capacity of the air (1.005 kJ/kg.K), and the temperature variation throughout the experiment (K). The amount of released heat was found equal to 1.265 kJ. Therefore, we concluded that the theoretical and experimental results were approximately similar, with a slight error (a relative error of about 2.2%). This error could be provoked by many reasons, mainly the measurement noises and the lack of stability of the temperature sensors, the box was not totally insulated, and not taking into consideration the variation of some factors with temperature such as air density. After validating the system's components models, they were connected and used in order to simulate and assess the performance of the proposed hybrid system (Figure 1) during one year. It is For the aim of comparing the theoretical and experimental results and thus validating this computation methodology, the temperature of the air trapped into the box (Figure 2) was measured. Figure 8 presents the variation of the temperature of the air trapped into the box (the average value of the three sensors). Then, for computing the heat released by the fuel cell, we assumed that the air's density did not change with temperature. So, the released heat was computed by Equation (1) using the density of the air (1.2 kg/m 3 ), the volume of the air trapped in the box (0.0478 m 3 ), the specific heat capacity of the air (1.005 kJ/kg.K), and the temperature variation throughout the experiment (K). The amount of released heat was found equal to 1.265 kJ. Therefore, we concluded that the theoretical and experimental results were approximately similar, with a slight error (a relative error of about 2.2%). This error could be provoked by many reasons, mainly the measurement noises and the lack of stability of the temperature sensors, the box was not totally insulated, and not taking into consideration the variation of some factors with temperature such as air density. After validating the system's components models, they were connected and used in order to simulate and assess the performance of the proposed hybrid system (Figure 1) during one year. It is worth noting that the electrolyzer was not modeled, but we used the manufacturer data, since it was not operating well during our experiment. capacity of the air (1.005 kJ/kg.K), and the temperature variation throughout the experiment (K). The amount of released heat was found equal to 1.265 kJ. Therefore, we concluded that the theoretical and experimental results were approximately similar, with a slight error (a relative error of about 2.2%). This error could be provoked by many reasons, mainly the measurement noises and the lack of stability of the temperature sensors, the box was not totally insulated, and not taking into consideration the variation of some factors with temperature such as air density. After validating the system's components models, they were connected and used in order to simulate and assess the performance of the proposed hybrid system (Figure 1) during one year. It is worth noting that the electrolyzer was not modeled, but we used the manufacturer data, since it was not operating well during our experiment.

Results of Long-Term Simulations of the Hybrid System
The proposed hybrid system was simulated during one year in order to evaluate its behavior and assess its performance when supplying the building with both valuable electricity and thermal energy. Concerning the PV system, its evaluation was based on the computation of some important performance indicators, namely, the total produced energy (E AC ), the system's efficiency (η system ), and the performance ratio (PR). Thus, the PV model was simulated under MATLAB and experimented under the same conditions, which were the ambient temperature and the solar irradiation presented in Figure 9. In fact, we have also presented the experimental results of the PV system since we already had related data. Therefore, the simulation and experimental results of the monthly energy generated by the PV system during one year were computed using Equation (13) and compared in order to assess the performance of the PV system. In addition, the system's efficiency was calculated using Equation (14) for the purpose of assessing its performance [59].
where E AC,m and E AC,d are the total monthly and daily energy produced by the PV system (kWh/kWp), respectively, N being the number of days of a month, G t,d stands for the monthly average daily solar irradiation received by the surface of the PV modules (kWh/m 2 ), and S PV represents the PV system's surface that is equal to 13.41 m 2 . Figure 10 depicts the simulated and experimental monthly energy generated by 1 kWp of the PV system (kWh/kWp) as well as the PV system's efficiency (%) over the whole period. We can obviously observe that the produced energy by the PV system varies between a minimum value of about 120 kWh/kWp, which corresponds to the month of November (winter period), and a maximum value of around 180 kWh/kWp, which corresponds to the month of April (spring period). Moreover, the annual produced energy by 1 kWp of the PV system is almost 1828 kWh/kWp per year, with an annual in-plane solar irradiation of about 2353 kWh/m 2 . On another note, the annual error between the simulation and the experimental monthly produced energy was around 3.9%. Thus, we conclude that the PV system's model gives a good approximation of the real-world scenarios with a slight error. Besides, we noticed that the average system's efficiency decreased from 15% (the manufacturer's efficiency) to 12.5% (the experimental efficiency), which was principally provoked by the weather conditions. ɳ , = , × * 100 (14) where , and , are the total monthly and daily energy produced by the PV system (kWh/kWp), respectively, being the number of days of a month, , stands for the monthly average daily solar irradiation received by the surface of the PV modules (kWh/m²), and represents the PV system's surface that is equal to 13.41 m².  Figure 10 depicts the simulated and experimental monthly energy generated by 1 kWp of the PV system (kWh/kWp) as well as the PV system's efficiency (%) over the whole period. We can obviously observe that the produced energy by the PV system varies between a minimum value of about 120 kWh/kWp, which corresponds to the month of November (winter period), and a maximum value of around 180 kWh/kWp, which corresponds to the month of April (spring period). Moreover, the annual produced energy by 1 kWp of the PV system is almost 1828 kWh/kWp per year, with an annual in-plane solar irradiation of about 2353 kWh/m². On another note, the annual error between the simulation and the experimental monthly produced energy was around 3.9%. Thus, we conclude that the PV system's model gives a good approximation of the real-world scenarios with a slight error. Besides, we noticed that the average system's efficiency decreased from 15% (the manufacturer's efficiency) to 12.5% (the experimental efficiency), which was principally provoked by the weather conditions. This decrease in the PV system's efficiency could also be noticed in other performance indicators, especially the performance ratio ( ), which is the ratio between the actual and the theoretically produced energy. It is an indicator that gives detailed information concerning the overall losses in the PV system, and it depends on several parameters (e.g., the quality of the cables, the efficiency of the converters). It is calculated using Equation (15): where is the energy produced by the PV system (kWh), stands for the solar radiation on the PV modules at STC (standard test conditions, which are 1000 W/m² solar irradiation and 25 °C temperature), represents the in-plane solar irradiation (kWh/m²), and , being the rated power of the installed PV system at STC. We found that the performance ratio was approximately equal to 79%. Accordingly, its high value discloses the good behavior and performance of the PV system during the whole year.
Then, the produced energy by the PV system could be either injected into the traditional electricity grid or used directly by an electrolyzer in order to produce hydrogen. For instance, the injected electricity into the electric grid could be used to produce hydrogen when needed in any other location. This electricity could be also used locally or stored in the form of hydrogen to be used when needed. The particularity of the first case is that the two processes (i.e., electricity production using PV panels and hydrogen generation by means of electrolyzer) could be carried out in different locations, which are connected through the traditional electricity grid. In this manner, electricity This decrease in the PV system's efficiency could also be noticed in other performance indicators, especially the performance ratio (PR), which is the ratio between the actual and the theoretically produced energy. It is an indicator that gives detailed information concerning the overall losses in the PV system, and it depends on several parameters (e.g., the quality of the cables, the efficiency of the converters). It is calculated using Equation (15): where E AC is the energy produced by the PV system (kWh), G STC stands for the solar radiation on the PV modules at STC (standard test conditions, which are 1000 W/m 2 solar irradiation and 25 • C temperature), G m represents the in-plane solar irradiation (kWh/m 2 ), and P PV, STC being the rated power of the installed PV system at STC. We found that the performance ratio was approximately equal to 79%. Accordingly, its high value discloses the good behavior and performance of the PV system during the whole year. Then, the produced energy by the PV system could be either injected into the traditional electricity grid or used directly by an electrolyzer in order to produce hydrogen. For instance, the injected electricity into the electric grid could be used to produce hydrogen when needed in any other location. This electricity could be also used locally or stored in the form of hydrogen to be used when needed. The particularity of the first case is that the two processes (i.e., electricity production using PV panels and hydrogen generation by means of electrolyzer) could be carried out in different locations, which are connected through the traditional electricity grid. In this manner, electricity would be transported instead of hydrogen, which has various benefits (e.g., security, less losses, less expensive). However, the only drawback of this process is that we would have more losses (e.g., the Moroccan electric power transmission and distribution losses are almost 15% [69]). Therefore, we will study the second case in this article, and thus all the produced energy by the PV system was directly converted into hydrogen. Hence, based on the electricity produced by the PV panels ( Figure 10); the electrolyzer's efficiency, which was around 63% [70]; and the heating value of hydrogen [71], the volume of produced hydrogen was determined. The volume of produced hydrogen in each month of the year is presented in Figure 11. Besides, the annual volume of hydrogen produced by 1 kWp of the PV system was almost 393 Nm 3 if the hydrogen was compressed under 10 bar. Finally, it is worth noting that the experimental results are unfortunately not presented here since the electrolyzer was not operating well, and we kept only the simulated results. This produced and stored hydrogen will be used by the above-mentioned fuel cell in order to provide the building with both valuable electrical and thermal energy. Therefore, we have computed the monthly produced electricity by the fuel cell, which is depicted in Figure 12, using the fuel cell's electrical efficiency as well as the hydrogen's heating value. As a result, an electrical energy as well as an equivalent amount of thermal energy of about 589 kWh/year could be generated. The produced electricity could be used in order to provide the building's appliances with the needed electricity, while the produced thermal energy could be used for either heating the building during the winter time using an air to air exchanger or cooling it during the summer period using a refrigeration process. This produced and stored hydrogen will be used by the above-mentioned fuel cell in order to provide the building with both valuable electrical and thermal energy. Therefore, we have computed the monthly produced electricity by the fuel cell, which is depicted in Figure 12, using the fuel cell's electrical efficiency as well as the hydrogen's heating value. As a result, an electrical energy as well as an equivalent amount of thermal energy of about 589 kWh/year could be generated. The produced electricity could be used in order to provide the building's appliances with the needed electricity, while the produced thermal energy could be used for either heating the building during the winter time using an air to air exchanger or cooling it during the summer period using a refrigeration process.
Accordingly, the electrical efficiency of the whole system was about 32%. However, its overall efficiency reached almost 64.5% when taking into account both produced electricity and thermal energy. These results show the effectiveness of the studied system when providing the building with both valuable electrical and thermal energy. electrical efficiency as well as the hydrogen's heating value. As a result, an electrical energy as well as an equivalent amount of thermal energy of about 589 kWh/year could be generated. The produced electricity could be used in order to provide the building's appliances with the needed electricity, while the produced thermal energy could be used for either heating the building during the winter time using an air to air exchanger or cooling it during the summer period using a refrigeration process. Accordingly, the electrical efficiency of the whole system was about 32%. However, its overall efficiency reached almost 64.5% when taking into account both produced electricity and thermal energy. These results show the effectiveness of the studied system when providing the building with both valuable electrical and thermal energy.

Conclusions and Perspectives
The aim of this work was towards efficient integration of a hybrid system into buildings together with its performance assessment. Therefore, we have focused on the modeling, simulations,

Conclusions and Perspectives
The aim of this work was towards efficient integration of a hybrid system into buildings together with its performance assessment. Therefore, we have focused on the modeling, simulations, experimentations, and performance assessment of the proposed hybrid system. The latter is comprised of a photovoltaic system, hydrogen electrolyzer, and PEM fuel cell operating as a CHP system to provide the building with both electricity and thermal energy. Several simulations and experiments were carried out first in order to show the accuracy of the developed models and assess the system's performance. Reported results showed the accuracy of the developed models in analyzing the system's performance. The validated models were afterwards used for long-run simulations (one year), and results were reported to show the effectiveness of the hybrid system. Accordingly, the system's electrical efficiency was almost 32%. However, its overall efficiency reached 64.5% when taking into account both produced electricity and thermal energy. These results show the effectiveness of the studied hybrid system when providing the building with both valuable electricity and thermal energy. It is noteworthy that the proposed system and computation methodology could be scaled up according to the system's size, and could be simulated and experimented under different conditions of any location. Our ongoing work puts more emphasis on the performance assessment of the proposed system when including the use of the thermal energy (i.e., the efficiency of the air to air exchanger (heating) and the refrigeration process (cooling)).