A Decentralized, Hybrid Photovoltaic-Solid Oxide Fuel Cell System for Application to a Commercial Building

: New energy solutions are needed to decrease the currently high electricity costs from conventional electricity-only central power plants in Cyprus. A promising solution is a decentralized, hybrid photovoltaic-solid oxide fuel cell (PV-SOFC) system. In this study a decentralized, hybrid PV-SOFC system is investigated as a solution for useful energy supply to a commercial building (small hotel). An actual load proﬁle and solar/weather data are fed to the system model to determine the thermoeconomic characteristics of the proposed system. The maximum power outputs for the PV and SOFC subsystems are 70 and 152 kWe, respectively. The average net electrical and total efﬁciencies for the SOFC subsystem are 0.303 and 0.700, respectively. Maximum net electrical and total efﬁciencies reach up to 0.375 and 0.756, respectively. The lifecycle cost for the system is 1.24 million USD, with a unit cost of electricity at 0.1057 USD/kWh. In comparison to the conventional case, the unit cost of electricity is about 50% lower, while the reduction in CO 2 emissions is about 36%. The proposed system is capable of power and heat generation at a lower cost, owing to the recent progress in both PV and fuel cell technologies, namely longer lifetime and lower speciﬁc cost.


Introduction
Efforts to increase energy efficiency have intensified over recent years due to the rapid increase of fossil fuel prices and also the need to decrease harmful emissions to the atmosphere [1]. Cogeneration allows the combination of various technologies to improve the fuel efficiency of electricity-only power plants or combined-heat-and-power (CHP) systems [2]. Earlier systems have included combined cycle power plants at large scale (10-100 MWe). In the resulting systems there exists a high level of complexity due to the increased number of parameters. Therefore, there is a need to apply advanced methodologies able to determine optimum solutions. However, this procedure becomes rigorous in case a number of parameters (e.g., thermodynamic, economic and environmental criteria) must be included [3]. Fuel cell technology has been proposed at the kW to the MW scale in a number of proposed systems. In lower temperature proton exchange membrane fuel cells (PEMFCs), CHP systems have been primarily applied at the kW scale, for smaller residential applications, where low-grade heat (recovered from the fuel cell exhaust) is usually adequate to cover residential load profiles, such as space heating and domestic hot water preparation [4,5]. These systems are sometimes operated jointly with vapor compression heat pumps to boost heat generation and/or to provide space cooling [6]. Such systems have been proposed for single-family households and in some cases for multi-unit theoretical predictions. This is done with a detailed thermoeconomic modeling of the components, and further on with their overall integration in the system model. Moreover, through the development of a cost model, a complete thermoeconomic analysis is facilitated to lay out the characteristics of the proposed hybrid system.

System Configuration
The proposed system, shown in Figure 1, was designed to fully fulfill an actual load profile for a commercial building. It includes a natural gas-fueled SOFC subsystem and a solar PV subsystem. The system also includes DC/AC inverters to convert the DC current generated by the PV and the SOFC subsystems to AC electricity prior to distribution to the buildings. In the SOFC subsystem, natural gas (NG) is compressed in the fuel compressor and sulfur is removed with the desulfurizer. The NG is preheated through heat exchanger (HEx) HEx1 before entering the SMR. The endothermic process in the SMR requires external heating, which is generated by a catalytic combustor. The synthesis gas (syngas) at the SMR exit is fed to the fuel cell anode. Air drawn from the atmosphere is filtered and blown to HEx3 for preheating and then fed to the fuel cell cathode. The fuel cell reaction in the SOFC stack generates electricity and also a hot mixture at the fuel cell exit. The hot exhaust mixture is fed to the combustor, along with natural gas from the natural gas supply and air. The flue gas exiting the SMR is used to provide heat for the four heat exchangers (HEx1-HEx4). HEx2 is used to generate steam for the SMR. HEx4 is used to provide low-grade heat externally, i.e., to heat water from recovered heat and supply it to the hot water storage tank. Through the hot water storage tank, hot water is provided to the buildings. At the exit of HEx4, the exhaust flue gas is released to the atmosphere after separation of water through a water separator.
Energies 2018, 11, 3 of 21 with an analysis beyond theoretical predictions. This is done with a detailed thermoeconomic modeling of the components, and further on with their overall integration in the system model. Moreover, through the development of a cost model, a complete thermoeconomic analysis is facilitated to lay out the characteristics of the proposed hybrid system.

System Configuration
The proposed system, shown in Figure 1, was designed to fully fulfill an actual load profile for a commercial building. It includes a natural gas-fueled SOFC subsystem and a solar PV subsystem. The system also includes DC/AC inverters to convert the DC current generated by the PV and the SOFC subsystems to AC electricity prior to distribution to the buildings. In the SOFC subsystem, natural gas (NG) is compressed in the fuel compressor and sulfur is removed with the desulfurizer. The NG is preheated through heat exchanger (HEx) HEx1 before entering the SMR. The endothermic process in the SMR requires external heating, which is generated by a catalytic combustor. The synthesis gas (syngas) at the SMR exit is fed to the fuel cell anode. Air drawn from the atmosphere is filtered and blown to HEx3 for preheating and then fed to the fuel cell cathode. The fuel cell reaction in the SOFC stack generates electricity and also a hot mixture at the fuel cell exit. The hot exhaust mixture is fed to the combustor, along with natural gas from the natural gas supply and air. The flue gas exiting the SMR is used to provide heat for the four heat exchangers (HEx1-HEx4). HEx2 is used to generate steam for the SMR. HEx4 is used to provide low-grade heat externally, i.e., to heat water from recovered heat and supply it to the hot water storage tank. Through the hot water storage tank, hot water is provided to the buildings. At the exit of HEx4, the exhaust flue gas is released to the atmosphere after separation of water through a water separator.   The main assumptions for the current study are the following: 1.
The proposed system operates in complete autonomy, i.e., it is not connected to a central power grid (no import/export of electricity).

2.
The maximum PV power output is set equal to the minimum electric load in the load profile to ensure no power is wasted. In turn, this value is used to size the SOFC subsystem. The system model is modeled in such a way to ensure that the proposed system is capable of completely covering the building load profile at all times, throughout the year. 3.
Heat losses are considered in the three main components of the SOFC subsystem, namely: The SMR reactor, the SOFC stack, and the catalytic combustor. Additionally, pressure losses are considered in every component of the SOFC subsystem.

4.
Additional heating (space heating and domestic hot water) is provided through natural gas-fired boilers, while space cooling is provided through electric vapor-compression heat pumps. This equipment is already available in the buildings and therefore its associated capital cost is not considered in the thermoeconomic modeling for this study.

5.
The hourly solar and ambient temperature data used in the simulation of the PV subsystem are based on a Typical Meteorological Year-TMY2 for Nicosia, Cyprus [27]. 6.
The consumption data system is applied for a small hotel with load profile data extracted from [28]. The load profile includes the following loads (all in an electrical energy basis): Fans, interior equipment, interior lights, space cooling, space heating, and domestic hot water. The load profile is shown in Figure 2. The main assumptions for the current study are the following: 1. The proposed system operates in complete autonomy, i.e., it is not connected to a central power grid (no import/export of electricity). 2. The maximum PV power output is set equal to the minimum electric load in the load profile to ensure no power is wasted. In turn, this value is used to size the SOFC subsystem. The system model is modeled in such a way to ensure that the proposed system is capable of completely covering the building load profile at all times, throughout the year. 3. Heat losses are considered in the three main components of the SOFC subsystem, namely: The SMR reactor, the SOFC stack, and the catalytic combustor. Additionally, pressure losses are considered in every component of the SOFC subsystem. 4. Additional heating (space heating and domestic hot water) is provided through natural gas-fired boilers, while space cooling is provided through electric vapor-compression heat pumps. This equipment is already available in the buildings and therefore its associated capital cost is not considered in the thermoeconomic modeling for this study. 5. The hourly solar and ambient temperature data used in the simulation of the PV subsystem are based on a Typical Meteorological Year-TMY2 for Nicosia, Cyprus [27]. 6. The consumption data system is applied for a small hotel with load profile data extracted from [28]. The load profile includes the following loads (all in an electrical energy basis): Fans, interior equipment, interior lights, space cooling, space heating, and domestic hot water. The load profile is shown in Figure 2.

System Modeling
The modeling of the components of the proposed hybrid system was based on first principles to accurately represent the coupling and operation of the system as realistically as possible. After modeling each component, the components were coupled together to form the SOFC subsystem. Subsequently, simulation of the PV subsystem generates PV data for the simulation of the overall system model. Additionally, a cost model was developed for the economic analysis of the proposed system. It includes all necessary cost functions and inputs needed for the calculation of capital costs, fuel cost, lifecycle cost and unit cost of electricity. The modeling of the system was developed with

System Modeling
The modeling of the components of the proposed hybrid system was based on first principles to accurately represent the coupling and operation of the system as realistically as possible. After modeling each component, the components were coupled together to form the SOFC subsystem. Subsequently, simulation of the PV subsystem generates PV data for the simulation of the overall system model. Additionally, a cost model was developed for the economic analysis of the proposed system. It includes all necessary cost functions and inputs needed for the calculation of capital costs, fuel cost, lifecycle cost and unit cost of electricity. The modeling of the system was developed with the software Engineering Equation Solver (EES)-Professional version. Hourly simulation data were generated for a complete year, i.e., 8760 hourly segments.

Photovoltaic Subsystem
The PV subsystem was based on the Hay-Davies-Klucher-Reindl (HDKR) modeling methodology [29], i.e., the total incident solar radiation on a tilted surface is calculated with a consideration of both the ground-reflected and the beam effects: In the PV array, the temperature was calculated with the relation (the effect of wind speed is considered negligible): The array's maximum power point efficiency is: The PV array's electricity output is:

SOFC Subsystem
The SOFC subsystem includes the fuel processing subsystem with the fuel pre-reformer (SMR reactor), four heat exchangers, SOFC stack and actuators. For the configuration shown in Figure 1, the inputs are given in Table 1. The fuel utilization factor was set at 0.92, and fuel cell temperature was set at 750 • C [10]. The temperature of fuel at the fuel preheater exit, the temperature of the reformate at the SMR reactor exit/anode inlet, and the temperature of the flue gas exiting the catalytic combustor were set at 450, 650 and 1005 • C, respectively [11]. HEx4 flue gas exit temperature was set at 55 • C because it must be at 25 • C above the dew point of the combustion product gases [10]. The steam-to-carbon ratio was set at 2.5, which although it is a relatively low value, the SOFC can treat CO as fuel [30], and therefore CO content does not need to be significantly reduced prior to anode inlet. HEx4 flue gas exit temperature 55 • C T 25 Water pump 1 inlet temperature 40 • C T 28 Hot water storage tank return temperature 40 • C T 29 Hot water storage tank supply temperature 65 • C SC Steam-to-carbon ratio 2.5 Energies 2018, 11, 3512 6 of 20

SMR Reactor
An SMR reactor configuration was assumed for the pre-reformer. The SMR reactor model is based on chemical equilibrium [10][11][12]. Two chemical reactions were included: SMR reaction (methane-steam), and water gas shift (WGS) reaction (carbon monoxide-steam) [31]. Since the SMR reaction is endothermic, heat must be supplied by an external source [32] (in this case from the catalytic combustor).
For the SMR reaction, i.e., CH 4 +H 2 O CO 2 +3H 2 , the overall change in Gibbs free energy is: The equilibrium constant at the given temperature for the SMR reaction is: For the WGS reaction, CO + H 2 O CO 2 +H 2 , the overall change in Gibbs free energy is: The equilibrium constant at the given temperature for the WGS reaction is: The molar flow output is defined as: . n re f ,out = .
The equilibrium constants for the aforementioned reactions are [33]: K wgs · y re f ,out,CO · y re f ,out,H 2 O = y re f ,out,CO 2 · y re f ,out,H 2 .
A molar flow rate balance for each species can be applied at the reformer inlet and outlet: . .
. n re f ,out,CO = X smr − X wgs .
The flue gas temperature at exit is calculated through an energy balance in the reformer: Energies 2018, 11, 3512 7 of 20

SOFC Stack
The SOFC stack model includes both the fuel cell reaction and direct internal reforming processes. For the latter, the reforming process takes place at the surface of the catalysts (anode), where hydrogen gas is mixed with steam before entering the anode [34]. Internal reforming is identical to the SMR reactor modeling equations, described in Section 3.2.1.
For the fuel cell reaction, the open circuit voltage is modeled as follows [10][11][12]: The reversible voltage is: The Gibbs free energy is: The activation losses are based on the Butler-Volmer equation, defined for the anode and cathode, respectively, to determine the current density: where i 0,ano and i 0,cat are the exchange current densities for the anode and cathode, respectively: The activation overvoltage is determined as the sum of anode and cathode losses: Concentration losses are the gradual losses due to the reactant depletion in the catalyst layer, and they are defined as the difference between the Nernst potential at the catalyst layer and the bulk flow at both anode and cathode [35]. The limiting current densities for hydrogen and oxygen species are defined as follows, respectively: where C H 2 ,0 and C O 2 ,0 are the concentration of species for hydrogen and oxygen, respectively. The concentration losses are [35]: The Ohmic losses are defined as the product of current density and Ohmic resistance: Based on the above definitions, the cell voltage can be defined as follows [34]: Fuel cell stack voltage, current, and power are defined as follows, respectively: .
The molar flow rate of oxygen at the inlet of the cathode can be calculated through an energy balance: .

Auxiliary Components
The auxiliary components are the actuators (air blower, fuel compressors and two water pumps), the catalytic combustor and the four heat exchangers. The actuators were modeled using fundamental equations, while the catalytic combustor model was based on an energy balance of products and reactants. The modeling of the heat exchangers was based on the Logarithmic Mean Temperature Difference (LMTD) method.

Overall System
The proposed hybrid system includes two prime movers for the generation of electrical energy. Additionally, in the case of the SOFC subsystem, heat is generated and recovered for external use in the buildings to satisfy the heating loads. Therefore, an algorithm must be included in the code of the system model to relate fuel cell power output, PV power output and power demand. Additionally, since the system is non-grid connected, it must be ensured that no excess power is generated from the PV subsystem.
If P load > P pvs then P f cs = P load − P pvs P pv,exc = 0 Else If P load < P pvs then P f cs = 0 P pv,exc = P pvs − P load Else P f cs = 0 P pv,exc = 0 EndIf EndIf (37) The inverter power losses for the PV subsystem are calculated as follows: . P loss,inv,pv = .
The net electrical power output for the SOFC subsystem is defined as follows: .
The net electrical efficiency for the SOFC subsystem can be based on the lower heating value (LHV) or the higher heating value (HHV), respectively [36]: The thermal efficiency is the ratio of recovered heat rate actually used to cover the building heating loads (fully or partly) to the chemical energy rate of the fuel consumed by the SOFC subsystem: The total SOFC subsystem efficiency is the sum of SOFC subsystem net electrical efficiency and thermal efficiency: The thermal-to-electric ratio is the ratio of recovered heat rate to net electrical power output: When the recovered heat from the SOFC subsystem is inadequate to cover the heating loads, additional heat must be generated externally: The total load profile electrical energy requirement is the sum of electricity required to operate the fans, the interior lights, the interior equipment, the space cooling, and supplementary heating: P load = P f an + P light + P equip + P cool + P heat,net .

Cost Model
A cost model was developed to determine the economic performance of the proposed hybrid system, based on the methodology found in [37]. The modeling equations are shown in Table 2, while the values of the constant parameters are given in Table 3. The specific cost of the PV array was set at 2.00 USD/W, which is based on approximate values given in [38]. The specific cost of the SOFC subsystem and the power subsystem were approximated from values given in [39], and they were set at 2.00 and 1.00 USD/W, respectively. The specific cost of the power subsystem included the two DC/AC inverters and the power conditioning components. The cost of fuel (i.e., natural gas) was set at 7.19 USD/MMBTU, which is the current cost in the European Union (EU) [40]. The lifetime was set at 20 years for the system (i.e., PV arrays and power subsystem) [19,37] and 5 years for the SOFC subsystem [39], with a fuel cell operation factor set at 0.50, since the fuel cell is operated for approximately 50% of the time. The cost of the hot water storage tank was based on an approximation from values given in [41]. The values for the remaining parameters were taken from [37]. Table 2. Cost modeling equations for the proposed hybrid system.

C loan
Cost of the loan (USD) Total worth of capital (USD) Tax paid on property (USD) Operation, maintenance and insurance cost (USD) Life cycle cost (USD) LCC = C twc + C prop + C omi + C tc f − D dep + D cred + D salv c el Unit cost of electricity (USD/kWh) P cs,li f e = N lt · P load,yr c el = LCC/P cs,li f e

Results and Discussion
In this section the system model is validated with available literature data. Then, the performance of the proposed hybrid system is presented in detail. Finally, the proposed hybrid system is compared with conventional and alternative system configurations to analyze and investigate its competitiveness in regard to key thermoeconomic parameters.

Validation
For the validation of the SOFC stack, relevant literature data from [42] were used. As shown in Figure 3, the literature data compare well against the simulation data generated by the system model, with only a small deviation in the results. The PV subsystem was validated in a previous publication by some of the authors [43].

Results and Discussion
In this section the system model is validated with available literature data. Then, the performance of the proposed hybrid system is presented in detail. Finally, the proposed hybrid system is compared with conventional and alternative system configurations to analyze and investigate its competitiveness in regard to key thermoeconomic parameters.

Validation
For the validation of the SOFC stack, relevant literature data from [42] were used. As shown in Figure 3, the literature data compare well against the simulation data generated by the system model, with only a small deviation in the results. The PV subsystem was validated in a previous publication by some of the authors [43].

Performance Characteristics of the Proposed Hybrid System
The proposed hybrid PV-SOFC system was sized in accordance with the requirements of the assumptions defined in Section 2. Based on these assumptions, the PV maximum power output is 70 kWe, while the SOFC maximum power output at full-load (i.e., design conditions) is 152 kWe. The average annual net electrical efficiency of the SOFC subsystem is 0.303, while total efficiency is 0.700. Maximum net electrical and total efficiencies can reach up to 0.375 and 0.756, respectively. In terms of annual useful energy generation, the electricity output (actual electricity delivered to the buildings) of the PV and SOFC subsystems is 135.9 and 451.2 MWh, respectively. The SOFC subsystem also provides 694.5 MWh of heating through heat recovery of the flue gas exhaust by the SOFC subsystem. This amount can almost completely cover the heating needs of the buildings, with only 7.2 MWh needed to be generated in addition. The current density at design conditions is 1228 A/m 2 . A summary of the performance characteristics of the proposed system is given in Table 4. For an illustration of the performance of the system, Table 5 includes the values for the thermophysical parameters of the system at full load conditions for the SOFC subsystem.

Performance Characteristics of the Proposed Hybrid System
The proposed hybrid PV-SOFC system was sized in accordance with the requirements of the assumptions defined in Section 2. Based on these assumptions, the PV maximum power output is 70 kWe, while the SOFC maximum power output at full-load (i.e., design conditions) is 152 kWe. The average annual net electrical efficiency of the SOFC subsystem is 0.303, while total efficiency is 0.700. Maximum net electrical and total efficiencies can reach up to 0.375 and 0.756, respectively. In terms of annual useful energy generation, the electricity output (actual electricity delivered to the buildings) of the PV and SOFC subsystems is 135.9 and 451.2 MWh, respectively. The SOFC subsystem also provides 694.5 MWh of heating through heat recovery of the flue gas exhaust by the SOFC subsystem. This amount can almost completely cover the heating needs of the buildings, with only 7.2 MWh needed to be generated in addition. The current density at design conditions is 1228 A/m 2 . A summary of the performance characteristics of the proposed system is given in Table 4. For an illustration of the performance of the system, Table 5 includes the values for the thermophysical parameters of the system at full load conditions for the SOFC subsystem. Overall heat transfer coefficient of HEx4 8396 W/m 2 ·K Table 5. Values for the thermophysical parameters of the proposed PV-SOFC hybrid system at full load conditions for the SOFC subsystem.

System Cost Analysis
A cost analysis of the proposed hybrid system is given in Table 6. In terms of capital cost, the highest cost is allocated for the purchase of the SOFC subsystem at 607,540 USD, while the cost for the PV arrays is 140,132 USD. The cost of the power subsystem is also significant at 221,951 USD, which means that it constitutes about 1 4 of the total system cost. The total cost of fuel for the operation of the system during its lifetime is estimated at 891,735 USD. Although natural gas prices are constantly fluctuating, it is not expected that this cost estimation will be significantly altered in the near future for the EU market, based on a statistical analysis of the prices for the last 10 years [40]. The lifecycle cost for the system is 1,241,369 USD, with a unit cost of electricity at 0.1057 USD/kWh. For a quantitative investigation of the possible merits of the proposed hybrid PV-SOFC system, four different case configurations were analyzed, in terms of thermoeconomic performance: (A) Power supply from a central power grid (conventional case), (B) power supply from a central power grid assisted with an on-site PV system, (C) power (and heat) supply from a decentralized SOFC system, and (D) power (and heat) supply from the proposed hybrid PV-SOFC system. A schematic representation of the four cases is given in Figure 4. The four cases can be compared in terms of two parameters: Unit cost of electricity and CO 2 emissions. The results from this comparison are shown in Figures 5 and 6, respectively. As observed, the proposed system outperforms all other configurations, in terms of both the unit cost of electricity and CO 2 emissions. In particular, in comparison to case A, the unit cost of electricity is about 50% lower (0.2128 vs. 0.1057 USD/kWh), while the reduction in CO 2 emissions is about 36% (673 vs. 428 g(CO 2 )/kWh).
In comparison to cases B and C, the additional capital cost for purchasing the SOFC subsystem and the PV subsystem, respectively, is well reasoned by the reduction in fuel consumption, hence on the unit cost of electricity (0.1700 USD/kWh (case B) and 0.1265 USD/kWh (case C)). Similarly, in terms of CO 2 emissions, the proposed system manages to significantly reduce emissions. For case B, power generation remains heavily dependent on inefficient central power grid supply, and therefore CO 2 emission generation remains high. For case C, CO 2 emissions are even higher than case B, because power (and heat) generation is completely dependent on the SOFC system. On an annual basis, the fuel consumption is 154,530 kg of natural gas for case C, compared to a reduced consumption of 115,848 kg for the proposed system in case D. In terms of lifecycle cost, for case C this is 1,468,209 USD, i.e., 226,840 USD higher than the equivalent cost for the proposed system in case D.
CO2 emission generation remains high. For case C, CO2 emissions are even higher than case B, because power (and heat) generation is completely dependent on the SOFC system. On an annual basis, the fuel consumption is 154,530 kg of natural gas for case C, compared to a reduced consumption of 115,848 kg for the proposed system in case D. In terms of lifecycle cost, for case C this is 1,468,209 USD, i.e., 226,840 USD higher than the equivalent cost for the proposed system in case D.

Conclusions
In this study a decentralized, hybrid PV-SOFC system is proposed for the fulfillment of a load profile for a commercial building (small hotel) in Cyprus. The system components are modeled in detail to allow a realistic simulation of the operation of the system. An actual load profile and solar/weather data are fed to the system model to determine the thermoeconomic characteristics of Figure 6. Comparison of the four cases under study in terms of CO 2 emissions (in g (CO 2 )/kWh).

Conclusions
In this study a decentralized, hybrid PV-SOFC system is proposed for the fulfillment of a load profile for a commercial building (small hotel) in Cyprus. The system components are modeled in detail to allow a realistic simulation of the operation of the system. An actual load profile and solar/weather data are fed to the system model to determine the thermoeconomic characteristics of the proposed system. The system is sized based on the requirements of the load profile, with maximum power outputs for the PV and SOFC subsystems at 70 and 152 kWe, respectively. The system operates efficiently throughout the whole year for a transient load profile. The average net electrical and total efficiencies for the SOFC subsystem are 0.303 and 0.700, respectively. Maximum net electrical and total efficiencies reach up to 0.375 and 0.756, respectively. The total contribution of the two subsystems on a yearly basis for the fulfillment of the load profile is at 135.9 and 451.2 MWh for the PV and the SOFC subsystems, respectively. Application of the proposed hybrid system is favored over conventional power generation with electricity-only central power stations for technical and economic reasons. The proposed system can operate more efficiently in terms of net electrical efficiency (especially at part-load operation over a heat engine-based power generator), and, more importantly, it can take advantage of the heat recovery capability of the SOFC subsystem. Additionally, fuel consumption is reduced significantly, primarily because of the integration of the PV subsystem, and also due to the elimination of transmission and distribution losses.
The cost analysis of the proposed system shows that in terms of capital cost, the highest cost is for the purchase of the SOFC subsystem (607,540 USD), while the cost for the PV arrays is 140,132 USD. The cost of the power subsystem, which is usually underestimated, is also significant at 221,951 USD. The total cost of fuel for the operation of the system during its lifetime is estimated at 891,735 USD. Although natural gas prices are constantly fluctuating, it is not expected that this cost estimation will be significantly altered in the near future for the EU market, based on a statistical review of the prices for the last 10 years [33]. The lifecycle cost for the system is 1,241,369 USD, with a unit cost of electricity at 0.1057 USD/kWh. The proposed system outperforms conventional and other system configurations, in terms of both the unit cost of electricity and CO 2 emissions. In comparison to the conventional case, the unit cost of electricity is about 50% lower (0.2128 vs. 0.1057 USD/kWh), while the reduction in CO 2 emissions is about 36% (673 vs. 428 g(CO 2 )/kWh). The additional capital cost for purchasing the PV and the SOFC subsystems is well reasoned by the reduction in fuel consumption, hence on the unit cost of electricity. Similarly, in terms of CO 2 emissions, the proposed system manages to significantly reduce emissions, because power generation is independent of the inefficient central power grid supply. Additionally, the integration of the PV subsystem allows a significant reduction in power generation from the SOFC subsystem during solar energy availability.

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