Thermoeconomic Optimization of a Hybrid Photovoltaic-Solid Oxide Fuel Cell System for Decentralized Application

: A small-scale, decentralized hybrid system is proposed for autonomous operation in a commercial building (small hotel). The study attempts to provide a potential solution, which will be attractive both in terms of e ﬃ ciency and economics. The proposed conﬁguration consists of the photovoltaic (PV) and solid oxide fuel cell (SOFC) subsystems. The fuel cell subsystem is fueled with natural gas. The SOFC stack model is validated using literature data. A thermoeconomic optimization strategy, based on a genetic algorithm approach, is applied to the developed model to minimize the system lifecycle cost (LCC). Four decision variables are identiﬁed and chosen for the thermoeconomic optimization: temperature at anode inlet, temperature at cathode inlet, temperature at combustor exit, and steam-to-carbon ratio. The total capacity at design conditions is 70 and 137.5 kWe, for the PV and SOFC subsystems, respectively. After the application of the optimization process, the LCC is reduced from 1,203,266 to 1,049,984 USD. This improvement is due to the reduction of fuel consumed by the system, which also results in an increase of the average net electrical e ﬃ ciency from 29.2 to 35.4%. The thermoeconomic optimization of the system increases its future viability and energy market penetration potential.


Introduction
Cogeneration of useful energy in the form of electricity, heating, and cooling has led to the development of combined heat and power (CHP) systems. These systems are available in different capacities, in order to fulfill a range of industrial, commercial, and residential purposes. To improve energy security, availability, and efficiency, the possibility of developing decentralized, autonomous energy systems has been considered in numerous studies [1][2][3][4]. Cogeneration plants usually emphasize on the production of electricity with the highest possible efficiency, while useful heating and cooling can be generated through heat recovery of the flue gas extracted from an electric generator (turbomachinery, fuel cells, etc.). It is also possible and desirable to generate electricity via renewable energy sources (RES), since RES-based systems typically offer production of electricity at zero emissions. A promising RES is solar photovoltaic (PV) technology, especially when applied in areas with high solar radiation [5,6]. However, its main drawback is the mismatch between electricity supply and demand [7,8]. Storage of electricity is possible via batteries, for example, but it can be problematic due to the high capital cost, short lifetime, limited capacity, and energy losses of energy storage devices [9,10].

Configuration of the Hybrid PV-SOFC System
The configuration of the hybrid PV-SOFC system is shown in Figure 1. It is modeled to fully satisfy the annual load profile of a small hotel in Cyprus. The load profile consists of several electric loads, which have been described by the authors in a previous publication [15]. The proposed system generates DC electricity through the SOFC and PV subsystems. When electricity generation from the PV is either unavailable, or inadequate, to fulfill the energy demand of the serviced buildings, the SOFC subsystem is operated. The latter is fueled with natural gas, which is compressed, desulfurized, and preheated through a heat exchanger. Steam is mixed with natural gas, and syngas is generated in the exit of the steam reformer. The hydrogen-rich syngas enters the fuel cell anode. Air is filtered and compressed through an air blower. Air is then preheated and directed to the cathode to generate DC electricity. The fuel cell exhaust mixture is directed to the catalytic combustor. Additional natural gas and air is fed to the combustor, as needed. The combustor generates flue gas, which is used to preheat the natural gas fuel supply (HEx1), generate steam (HEx2), preheat the incoming air flow (HEx3), and generate useful heating through heat recovery (HEx4). The DC/AC inverters convert DC electricity to AC electricity. The main assumptions of the modeling methodology for the study can be found in [15]. Appl. Sci. 2019, 9, x FOR PEER REVIEW 3 of 13

System Modeling
The simulation is accomplished with the development of all component models and their subsequent integration to a total energy system model. The thermoeconomic modeling and optimization is developed in the software Engineering Equation Solver (EES)-Professional version. The component models include the PV and SOFC subsystems. The latter includes the Steam Methane Reformer (SMR) reactor, the fuel cell stack, and auxiliary components (air blower, fuel compressor, water pumps, catalytic combustor, and heat exchangers). All the aforementioned component model equations (except the SOFC stack, which is described below), the overall system balance equations, and cost model equations are identical to those already presented by the authors in a previous publication [15]. In the end of this section, the applied optimization strategy is described in detail.

SOFC Subsystem
The SOFC is of the planar type and its model is presented in detail in this section. The values of the constant parameters in the SOFC subsystem are given in Table 1. The fuel cell effective crosssectional area is set at 100 cm 2 [16]. The fuel utilization factor and the fuel cell operating temperature are fixed at 0.9 and 800 °C, respectively [17][18][19]. The temperature at the fuel preheater exit is fixed at 450 °C [20], while the temperature of the flue gas at the HEx4 exit is fixed at 55 °C [21]. The temperature at the inlet of water pump 1 is fixed at 40 °C [15]. The return and supply temperature values in the hot water storage tank are fixed at 40 and 65 °C, respectively [20].

System Modeling
The simulation is accomplished with the development of all component models and their subsequent integration to a total energy system model. The thermoeconomic modeling and optimization is developed in the software Engineering Equation Solver (EES)-Professional version. The component models include the PV and SOFC subsystems. The latter includes the Steam Methane Reformer (SMR) reactor, the fuel cell stack, and auxiliary components (air blower, fuel compressor, water pumps, catalytic combustor, and heat exchangers). All the aforementioned component model equations (except the SOFC stack, which is described below), the overall system balance equations, and cost model equations are identical to those already presented by the authors in a previous publication [15]. In the end of this section, the applied optimization strategy is described in detail.

SOFC Subsystem
The SOFC is of the planar type and its model is presented in detail in this section. The values of the constant parameters in the SOFC subsystem are given in Table 1. The fuel cell effective cross-sectional area is set at 100 cm 2 [16]. The fuel utilization factor and the fuel cell operating temperature are fixed at 0.9 and 800 • C, respectively [17][18][19]. The temperature at the fuel preheater exit is fixed at 450 • C [20], while the temperature of the flue gas at the HEx4 exit is fixed at 55 • C [21]. The temperature at the inlet of water pump 1 is fixed at 40 • C [15]. The return and supply temperature values in the hot water storage tank are fixed at 40 and 65 • C, respectively [20].

Parameter Description Value
A f c Fuel cell effective cross-sectional area 100 cm 2 U f Fuel utilization factor 0.9 T f c Fuel cell operating temperature 800 • C T 6 Temperature at fuel preheater exit 450 • C T 22 Temperature of flue gas at HEx4 exit 55 • C T 25 Temperature at water pump 1 inlet 40 • C T 28 Return temperature in hot water storage tank 40 • C T 29 Supply temperature in hot water storage tank 65 • C The SOFC stack model includes the direct internal reforming (DIR) and fuel cell reaction processes. The DIR process takes place at the catalyst surface [22] and it is equivalent to the steam reforming process [15]. A 0-dimensional model is considered, to allow coupling to the other system components. Therefore, changes in thermodynamic and electrochemical properties are considered negligible along the fuel cell. The SOFC temperature and pressure are fixed at a constant value throughout the fuel cell, because the outlet radial temperature gradients are negligible. Inlet air provides oxygen for the fuel cell reaction and cooling for the fuel cell stack (N 2 is considered inert both at cathode and anode).
The Gibbs free energy, the reversible voltage, and the open circuit voltage, respectively, are defined as [23]: The activation losses in the anode and cathode can be determined as follows [17]: The Ohmic losses are determined as follows [17]: where σ ele , σ ano , σ cat and σ int are the ionic conductivity of the electrolyte, the electronic conductivity of the anode, the electronic conductivity of the cathode, and the electronic conductivity of the interconnect, respectively. These are determined as follows: The concentration losses are determined at the anode and cathode, respectively [17]: Then the cell voltage can be calculated: The voltage, the current, and the power are determined as follows, respectively: .
The hydrogen consumption molar flow rate is determined as follows: . n con,h2 = I f c · n cell n e · F An energy balance on the fuel cell stack allows determination of the amount of oxygen that needs to be supplied to the cathode:

Optimization Strategy
The optimization strategy aims to minimize the LCC and is applied to the system model, assuming operation at design (full load) conditions throughout the annual load profile and the total system lifetime [24]. Specifically, the objective function of the optimization problem is defined as follows: The decision variables for the optimization are chosen based on realistic ranges, where they can be varied to optimize the objective function [25,26]. Four parameters have been identified as the decision variables of the optimization problem: the temperature at anode inlet (T 8 ), the temperature at cathode inlet (T 13 ), the temperature at combustor exit (T 17 ), and the steam-to-carbon ratio (SC). The minimum values for the bounds of T 8 and T 13 are fixed at 700 • C, because the temperature difference between the fuel cell and the air/fuel inlet must not exceed 100 • C, to avoid thermal stresses in the SOFC [21]. The steam-to-carbon ratio can vary between 2.5-4.0, which is a realistic operating range for SOFC systems [27]. The minimum, maximum, and initial values of the decision variables for the optimization process in the current study are tabulated in Table 2.

Model Validation
The SOFC stack model is validated by variation of the current density from the design value of 1000 A/m 2 to part-load conditions, as shown in the polarization curve of Figure 2. The simulation shows good agreement with typical polarization curves for SOFC models found in the literature [27]. The PV model has been validated by the authors in a previous publication [15].
Appl. Sci. 2019, 9, x FOR PEER REVIEW 6 of 13 in the SOFC [21]. The steam-to-carbon ratio can vary between 2.5-4.0, which is a realistic operating range for SOFC systems [27]. The minimum, maximum, and initial values of the decision variables for the optimization process in the current study are tabulated in Table 2.

Model Validation
The SOFC stack model is validated by variation of the current density from the design value of 1000 A/m 2 to part-load conditions, as shown in the polarization curve of Figure 2. The simulation shows good agreement with typical polarization curves for SOFC models found in the literature [27]. The PV model has been validated by the authors in a previous publication [15].

Thermoeconomic Optimization: Application to the Hybrid PV-SOFC System
The variation of the temperature at anode inlet, the temperature at cathode inlet, the temperature at combustor exit, and the steam-to-carbon ratio throughout the optimization process are shown graphically in terms of number of generations in Figures 3-6, respectively. The evolution of the objective function LCC 1 is shown graphically in terms of the number of generations in Figure 7. 1 The values of cost parameters are given in US dollars (USD).

Thermoeconomic Optimization: Application to the Hybrid PV-SOFC System
The variation of the temperature at anode inlet, the temperature at cathode inlet, the temperature at combustor exit, and the steam-to-carbon ratio throughout the optimization process are shown graphically in terms of number of generations in Figures 3-6, respectively. The evolution of the objective function LCC 1 is shown graphically in terms of the number of generations in Figure 7.         Table 3 includes the resulting values of the four decision variables and the objective function after the application of the optimization process at design conditions. Table 4 includes the thermophysical parameter values, corresponding to the nodes shown in Figure 1, for operation at nominal conditions (after the application of the optimization process). Table 3. Optimum values of the decision variables and the objective function after the application of the optimization process at design conditions.

Parameter Description
Optimum value 8 T   Table 3 includes the resulting values of the four decision variables and the objective function after the application of the optimization process at design conditions. Table 4 includes the thermophysical parameter values, corresponding to the nodes shown in Figure 1, for operation at nominal conditions (after the application of the optimization process). Table 3. Optimum values of the decision variables and the objective function after the application of the optimization process at design conditions.

Parameter Description
Optimum Value

Discussion
The results of the thermoeconomic optimization applied to the proposed hybrid system suggest that significant room for improvement of the objective function is available. The optimization strategy has successfully decreased the LCC to a more competitive value, in comparison to the pre-optimization performance of the system model. To find the actual effect of the applied optimization methodology on the system, the optimum values of the decision variables are applied to the simulation model for the actual varying annual load profile. Table 5 includes all the key economic and efficiency parameter values before and after the application of the thermoeconomic optimization process, when the optimum decision variables are applied on the actual load profile. The LCC is reduced from 1,203,266 USD to 1,049,984 USD, i.e., a 12.7% improvement. The reason for this significant cost reduction is almost entirely due to the minimization of the total cost of fuel, which decreases by 17.3%. The effect on the other cost parameters is almost insignificant (0.5%). The fuel reduction also has a positive effect on the net electrical efficiency, which increases by 6.2% (average value). The maximum and minimum values for the net electrical of the hybrid system are 38.2 and 30.6%, respectively. The overall system improvement, both in terms of economics and thermodynamics, is due to the better use of resources (i.e., fuel consumption). In other words, the system can now operate more efficiently because of the reduction of exhausted flue gas from the system. Specifically, at design conditions, the molar flow rate of the flue gas exhausted ( . n 23 ) is reduced from 0.0365 to 0.0295 kmol/s. Finally, the application of the optimization strategy results in an important decrease of fuel consumption and carbon emissions. The mass flow rate of methane per year is reduced from 31.7 to 26.2 kg/s, while the mass flow rate of CO 2 per year is reduced from 67.3 to 58.0 kg/s. Table 5. Key economic and efficiency parameter values for the hybrid PV-SOFC system after the application of the thermoeconomic optimization process (actual load profile).

Parameter Description Value
Before After Improvement

Conclusions
A small-scale PV-SOFC system was proposed for autonomous operation in a small hotel. The total capacity at design conditions was 70 and 137.5 kWe (kW of electric power), for the PV and SOFC subsystems, respectively. A thermoeconomic optimization strategy was applied to the developed simulation model to minimize the lifecycle cost (LCC) to a more competitive value, which could improve its future viability and energy market penetration. After the application of the optimization process, the LCC was reduced from 1,203,266 to 1,049,984 USD; an improvement of 12.7%. This improvement is due to the reduction of fuel consumption, which also improves average net electrical efficiency from 29.2% to 35.4%. The results accomplished by this study suggest a significant potential for the introduction of the proposed configuration in the energy infrastructure for decentralized, autonomous energy systems.