Development of an Integrated Structure for the Tri-Generation of Power, Liquid Carbon Dioxide, and Medium Pressure Steam Using a Molten Carbonate Fuel Cell, a Dual Pressure Linde-Hampson Liquefaction Plant, and a Heat Recovery Steam Generator

: Due to the increase in energy consumption and energy prices, the reduction in fossil fuel resources, and increasing concerns about global warming and environmental issues, it is necessary to develop more efﬁcient energy conversion systems with low environmental impacts. Utilizing fuel cells in the combined process is a method of refrigeration and electricity simultaneous production with a high efﬁciency and low pollution. In this study, a combined process for the tri-generation of electricity, medium pressure steam, and liquid carbon dioxide by utilizing a molten carbonate fuel cell, a dual pressure Linde-Hampson liquefaction plant and a heat recovery steam generator is developed. This combined process produces 65.53 MW of electricity, 27.8 kg/s of medium pressure steam, and 142.9 kg/s of liquid carbon dioxide. One of the methods of long-term energy storage involves the use of a carbon dioxide liquefaction system. Some of the generated electricity is used in industrial and residential areas and the rest is used for storage as liquid carbon dioxide. Liquid carbon dioxide can be used for peak shavings in buildings. The waste heat from the Linde-Hampson liquefaction plant is used to produce the fuel cell inlet steam. Moreover, the exhaust heat of the fuel cell and gas turbine would be used to produce the medium pressure steam. The total efﬁciency of this combined process and the coefﬁcient of performance of the refrigeration plant are 82.21% and 1.866, respectively. The exergy analysis of this combined process reveals that the exergy efﬁciency and the total exergy destruction are 73.18% and 102.7 MW, respectively. The highest rate of exergy destruction in the hybrid process equipment belongs to the fuel cell (37.72%), the HX6 heat exchanger (8.036%), and the HX7 heat exchanger (6.578%). The results of the sensitivity analysis show that an increase in the exit pressure of the V1 valve by 13.33% would result in an increase in the refrigeration energy by 2.151% and a reduction in the refrigeration cycle performance by 9.654%. Moreover, by increasing the inlet fuel to the fuel cell, the thermal efﬁciency of the whole combined process rises by 18.09%, and the whole exergy efﬁciency declines by 12.95%.


Introduction
Fuel cells are one of the suitable alternatives to conventional power generation technologies from fossil fuels. They have been of interest to many researchers in recent years due to their low pollution and high performance [1]. Molten carbonate fuel cells (MCFC), as well as solid oxide fuel cells, are high-temperature fuel cells. The exhaust gases from them can be used in devices such as gas turbines and heat recovery steam generators (HRSG) [2]. Various studies on the modeling of MCFCs have been performed. The complexity of the studied models ranges from simple and brief models to more complex 3-dimensional models of fuel cells [3]. Haghighat et al. [4] analyzed a hybrid MCFC and a gas turbine Carbon dioxide has better physical properties than air for storage since it has a critical temperature and pressure of 31 • C and 37.37 MPa, respectively. Therefore, liquid CO 2 can be used as a working fluid in energy storage structures [16,17]. Aliyon et al. [18] investigated CO 2 liquefaction units by several methods, such as the utilization of the vapor compression refrigeration cycle, the water-ammonia absorption refrigeration unit, and the Linde-Hampson refrigeration cycles. The simulation results revealed that the performance coefficient was the lowest (0.537) for the water-ammonia absorption refrigeration unit and was the highest (2.617) for the compression refrigeration cycle. Xu et al. [16] developed two combined power generation systems for peak hours by liquefying carbon dioxide and air. Solar collectors were used to heat the CO 2 liquefaction unit. The round trip and exergy efficiencies of the CO 2 liquefaction unit were 45.35% and 67.2%, respectively. Shakib et al. [19] utilized the waste heat from a gas turbine to produce steam in the HRSG system. This steam was implemented in the multi-stage flash distillation. Then they optimized the main parameters of the developed hybrid system by the genetic algorithm and particle swarm optimization. With these optimization algorithms, the exergy efficiency of the developed integrated structure increased from 3.5% to 4%. Moreover, Hosseini et al. [20] performed exergy, economic, and risk analyses on a freshwater production integrated structure using a gas turbine to produce steam in the HRSG system and the multistage flash distillation. Akrami et al. [21] performed exergy and exergoeconomic analyses of three power and CO 2 generation integrated structures using a new hybrid cycle of biomass gasification, a gas turbine, a MCFC, an organic Rankine unit, and a CO 2 separation and liquefaction unit. The results of the exergy assessment of the three hybrid structures revealed that the maximum value of exergy efficiency of the entire integrated unit and its specific power consumption parameters were 43.95% and 854.9 kgCO 2 /MWh, respectively. Ebrahimi-Moghadam et al. [22] investigated a combined structure for the cogeneration of power and heat in a residential complex using the Kalina power cycle and a gas turbine. The simulation results showed that the thermal efficiency, exergy efficiency, and period of return were 69.43%, 37.90%, and 3.34 years, respectively. Shokouhi Tabrizi et al. [23] applied a CO 2 recompression Brayton unit and solar collectors in a city gas station to preheat the natural gas and generate power. The results showed that the exergy efficiency and period of return were 41% and 4 years, respectively. Ebrahimi-Moghadam et al. [24] investigated a hybrid configuration for the tri-generation of power, cooling, and heat by the Kalina unit and the ejector refrigeration plant and a regenerative turbine system. The simulation results revealed that the thermal, electrical, and exergy efficiencies were 36.55%, 37.87%, and 38.85%, respectively. Ghorbani et al. [25] investigated a new liquid CO 2 cryogenic energy storage system for peak shaving by a post-combustion structure, a Linde-Hampson liquefaction unit, an organic Rankine, and gas turbine plants. Pinch, exergy, and sensitivity analyses were used to evaluate the developed integrated structure. Ebrahimi-Moghadam et al. [26] proposed a novel integrated structure for the cogeneration of the power and heat required by the gas pressure drop station. The simulation results revealed that the thermal and exergy efficiencies were 42.51%, and 38.17%, respectively. Ebrahimi-Moghadam et al. [27] evaluated a novel hybrid system for the trigeneration of power, cold, and heat using a gas turbine, a Kalina power production unit, and ejector systems. The combined process produced 1025.9 kW of power, 1642.3 kW of heat, and 304.9 kW of cooling.
There are some studies in the literature that suggest useful guidelines to enhance the efficiency of cogeneration systems. In fact, the principal purpose of the reported investigations is to optimize the necessary thermal energy, generated electricity, and economic prices of the process and finally to decrease the negative environmental impacts. To the best of our knowledge, no systematic investigation has been conducted in relation to the use of the waste heat and electricity generated in fuel cell systems in the CO 2 liquefaction industry. In this study, an integrated structure for the tri-generation of electricity, steam, and liquid carbon dioxide is developed using a molten carbonate fuel cell, a dual pressure Linde-Hampson liquefaction plant, and a heat recovery steam generator. Energy, exergy, and sensitivity assessments are employed to investigate the thermodynamic laws and the main parameters of the combined system. After compression and preheating, the liquid CO 2 can be used for peak shaving in industrial and residential areas.

Process Description
Fuel cells are a competitive alternative to conventional power generation technologies from fossil fuels due to their low pollution and high efficiency, especially in small-scale and scattered power generation [28]. The MCFCs are kind of high-temperature fuel cell which have many advantages but have received less attention. The high operating temperature of this type of fuel cell makes it compatible with different fuels. On the other hand, the high heat generated by the fuel cell can be recovered for heating applications [29]. In MCFCs, hydrogen is delivered to the anode and the cathode receives carbon dioxide and air to produce heat and electricity due to chemical reactions within the fuel cell. Since the MCFC operates at high temperatures, different hydrocarbon fuels can be used. Before the injection of hydrocarbon fuels, a reformer can be used to convert hydrocarbons to hydrogen using water vapor. Due to the endothermic reaction of the reformer, the heat required for this reaction is provided by electrochemical reactions inside the fuel cell. In the developed MCFC unit, the chemical energy of the natural gas is converted to power. The block diagram for the production of liquid CO 2 , electricity, and medium pressure steam using a molten carbonate fuel cell, a Linde-Hampson liquefaction plant, and a HSRG is illustrated in Figure 1. This combined configuration produces 65.53 MW of electricity, 27.80 kg/s of medium pressure steam, and 142.9 kg/s of liquid carbon dioxide.  The natural gas (stream A29) at a 25 • C temperature, 2.5 bar pressure and a 4.202 kg/s flow rate is mixed with the A28 steam stream and enters the reformer through stream A31 after it has been preheated in the HX5 heat exchanger by the return stream from the reformer. The required heat for generating the A28 vapor stream is provided by the hightemperature outflow from the C2 compressor located in the dual pressure Linde-Hampson liquefaction plant. In the reformer, an approximate amount of the methane of stream A31 reacts to produce the hydrogen required for the anode pole reaction. The outflow from the reformer (stream A32) enters the anode pole after transferring heat to the reformer inflow in the HX5 heat exchanger and produces electrons during the reaction with the carbonate ion. The process flow diagram of the combined configuration for the tri-generation of liquid carbon dioxide, electricity, and medium pressure steam using the molten carbonate fuel cell, a liquefaction plant, and a HSRG is presented in Figure 2.   The outflow from the anode is mixed with the A4 return flow and enters the catalytic burner through stream A36 so that its entire CO is converted to CO 2 . Air is compressed in the C2 compressor and then it is mixed with the outflow from the catalytic burner and enters the cathode. In the cathode, carbonate ions are produced due to a combination of carbon dioxide and oxygen to maintain the electrolyte concentration. The carbonate ion is carried from the cathode to the anode via the electrolyte. In the integrated structure, the carbonate ion is separated from the cathode outflow with a splitter and enters the anode pole. Part of the fuel cell outflow is used as a return flow to be combined with the outflow from the anode to enter the catalytic burner, and the rest enters the T1 gas turbine through stream A47 and produces 40.02 MW of electricity.
The A48 outflow of the gas turbine enters the economizer and the evaporator with a temperature of 461.4 • C to produce 27.8 kg/s of steam at the temperature of 182.7 • C and 10 bar pressure. The outflow from the evaporator then enters the environment at a temperature of 92.22 • C. The reactions that occur at the MCFC are listed below.
η an = 2.27 × 10 −9 e ( E act,an Note that the reformer is applied to convert natural gas to the required hydrogen for the reactions of the molten carbonate fuel cell, so that based on the reactions with steam, CH 4 is converted to CO 2 and H 2 [6].

•
Steam reforming reaction: • Water-gas shift reaction: • The overall reaction: Rxn Heat: 1.647 × 10 5 kJ/kmol. Information on the structure of the MCFC system is given in [6]. Stream A1, which contains 97 mol % of carbon dioxide, enters the C2 compressor, and its pressure rises by 19.70 bar. Then stream A2 enters the HX4 heat exchanger with a temperature of 281.6 • C and provides the required heat to the steam entering the reformer. The discharged A3 stream is in a two-phase state and by the D2 flash drum, most of the carbon dioxide impurity is extracted in liquid form. The carbon dioxide flow is mixed with the A25 return stream and enters the C3 compressor through stream A6 at a temperature of 19.73 • C, 19.70 bar pressure, and a 179.6 kg/s flow rate, and its pressure increases by 39 bar.
The outflow of the C3 compressor goes to the HX1 heat exchanger and its temperature is reduced by 32 • C before it enters the D3 flash drum. The gas stream extracted from the top of the D3 flash drum is mixed with the A20 return stream and enters the C4 compressor through stream A12 and its pressure increases by 100 bar. The outflow of the C4 compressor enters the HX2 heat exchanger so that its temperature is reduced by 35 • C before it enters the HX3 heat exchanger and its temperature is reduced by the return cold flow. Stream A15 coming out of the heat exchanger goes to the V1 valve and with a pressure drop of 39 bar, its temperature decreases by 3.323 • C and stream A16 enters the D4 flash drum. The gas flow coming out of the D4 flash drum is used to chill the HX3 heat exchanger and then as a return flow, it is mixed with the carbon dioxide which has entered the second part. The liquid carbon dioxide outflow of the D4 flash drum goes to the V2 pressure reducing valve and its pressure decreases to 19.70 bar. With the decrease in the pressure, the carbon dioxide flow would be in a two-phase state. Stream A22 enters the D5 flash drum and the gas stream coming out of the top of the D5 flash drum is used to cool the HX3 heat exchanger and as a return flow, it blends with the CO 2 which has entered the first part. The liquid CO 2 flow from the bottom of the D5 flash drum is stored as the final product. More theory and background information on the Linde-Hampson liquefaction plant can be seen in [18].

Exergy Analysis
Exergy analysis by combining the first and second laws of thermodynamics has become one of the most important instruments for the quantitative and qualitative study of energy consumption in processes. The exergy of an energy source refers to the maximum work achievable when the process of the thermodynamic state identifies its source as dead [34]. The dead state is actually the same ambient conditions that are typically considered to be 25 • C and 1 atm pressure. The total exergy rate of a stream can be examined as the sum of the physical exergy and chemical exergy [35]: The physical exergy of a stream can be determined as follows [36,37]: where γ is the activity coefficient of component i, which can be more or less than one, and e 0 i refers to the standard chemical exergy of any compound in the stream. The exergy balance relation for the control volume of the system is obtained from Equation (17)

Results and Discussion
A fuel cell is an electrochemical device that directly turns chemical energy into electrical energy. This conversion continues as long as fuel is injected into the fuel cell [39]. In fact, by eliminating the intermediate stages of energy conversion, fuel cells provide more energy conversion efficiency. In thermal power plants, chemical energy is first turned into heat and then into mechanical energy and finally, the mechanical energy is converted to electrical energy. Each of these stages causes an efficiency loss. Fuel cells do not store energy like a battery; rather, in a fuel cell, a state of energy is converted to another state in a manner that ensures that no material is consumed within the cell. Furthermore, the battery density is lower than the fuel cell, and the process of charging the battery is more sophisticated than filling the fuel tank of the fuel cell. In batteries, the potential of electrochemical conversions decreases after several charges, while in fuel cells there is no such limit. One method of storing energy during off-peak hours is to convert some of the electrical energy into liquid carbon dioxide. On the other hand, during the peak hours of energy consumption, power can be generated by cooling and compressing the liquid carbon dioxide to be used at peak hours. In this paper, molten carbonate fuel cells, the Linde-Hampson liquefaction plant, and a HSRG were used to produce liquid carbon dioxide, electricity, and medium pressure steam. The results of the energy, exergy, and sensitivity examinations of the combined process are presented as follows.

Energy Analysis
This hybrid structure of a molten carbonate fuel cell, a dual pressure Linde-Hampson liquefaction plant, and a HSRG produces 65.53 MW of electricity, 27.80 kg/s of medium pressure steam, and 142.9 kg/s of liquid CO 2 . The simulation of the hybrid structure was performed using the HYSYS V10 software and MATLAB V10 programming. Based on the obtained results from the literature, the Peng-Robinson equation of state was used for the simulation. The thermodynamic conditions of the streams in the integrated structure of liquid carbon dioxide, electricity, and medium pressure steam production using a molten carbonate fuel cell, a Linde-Hampson liquefaction plant, and a HSRG, are presented in Table 1. The validation of the several parameters of the Linde-Hampson liquefaction cycle in this paper with reference [18] is illustrated in Figure 3. The results of the simulation analysis with the HYSYS software show that the fuel cell performance coefficient, total cooling duty, and refrigeration energy differ from the evaluated values reported in the reference article by 2.07%, 0.604%, and 2.022%, respectively. This indicates that the simulation results are acceptable.  Figure 4 shows the voltage variations in the current density change of the current density in the MCFC for this integrated structure and in [6]. The temperature and pressure of the fuel cell in the validation were considered to be 650 °C and 1 bar, respectively. The results of the validation show that the simulated fuel cell in this study is in an appropriate  Figure 4 shows the voltage variations in the current density change of the current density in the MCFC for this integrated structure and in [6]. The temperature and pressure of the fuel cell in the validation were considered to be 650 • C and 1 bar, respectively. The results of the validation show that the simulated fuel cell in this study is in an appropriate agreement with the one reported in [6]. To achieve the optimal design, it is necessary to evaluate the influence of each variable in the fuel cell system on its performance by sensitivity analysis.
Sustainability 2021, 13, x FOR PEER REVIEW 11 of Figure 4. The voltage variation with the current density change of the current density in the MCFC for this integrated structure and in [6].
As the current density increases, the cell voltage decreases due to irreversible loss The power density rises with the increment in current density till it reaches a maximu value and then it starts to decrease. As the current density increases, the cell efficien decreases due to the irreversible losses, and as a result, less of the generated heat is co verted into electricity by the electrochemical reaction, and the rest is released as wa heat. This phenomenon leads to a reduction in power production and therefore, more is needed to remove the excessive heat from the cell to cool it.  Figure 5 shows the variations in power density, voltage, and efficiency of the fuel cell with the variation of current density in the operating conditions of the fuel cell used in the combined system. The design temperature and pressure of the fuel cell were 650 • C and 2.5 bar, respectively. The rate of voltage drop due to mass transfer and activation increased with the increase in current density.
As the current density increases, the cell voltage decreases due to irreversible losses. The power density rises with the increment in current density till it reaches a maximum value and then it starts to decrease. As the current density increases, the cell efficiency decreases due to the irreversible losses, and as a result, less of the generated heat is converted into electricity by the electrochemical reaction, and the rest is released as waste heat. This phenomenon leads to a reduction in power production and therefore, more air is needed to remove the excessive heat from the cell to cool it.
The inlet airflow to the cathode has an important impact on the regulation of the cell temperature, the size of the fuel cell unit, and the electricity consumption in the gears. The operational specifications of the fuel cell utilized in the hybrid process of producing liquid carbon dioxide, electricity, and medium pressure steam are presented in Table 2. The current density, fuel utilization, and cell voltage in the molten carbonate fuel cell were 0.2000 A·cm −2 , 0.6800, and 0.7917 V, respectively. The results of the simulation of the hybrid process for the tri-generation of electricity, medium pressure steam, and liquid CO 2 are presented in Table 3. The inlet airflow to the cathode has an important impact on the regulation of the cell temperature, the size of the fuel cell unit, and the electricity consumption in the gears. The operational specifications of the fuel cell utilized in the hybrid process of producing liquid carbon dioxide, electricity, and medium pressure steam are presented in Table 2. The current density, fuel utilization, and cell voltage in the molten carbonate fuel cell were 0.2000 A·cm⁻ 2 , 0.6800, and 0.7917 V, respectively. The results of the simulation of the hybrid process for the tri-generation of electricity, medium pressure steam, and liquid CO2 are presented in Table 3.   The MCFC cycle and the net overall thermal and overall thermal efficiencies of the developed combined configuration are 52.13%, 63.43%, and 82.21%, respectively. The simulation results reveal that the coefficient of performance and refrigeration energy in the dual pressure Linde-Hampson liquefaction plant is 1.866 and 26.54 MW, respectively.

Exergy Analysis
The exergy assessment of the combined process was performed by the relationship between the HYSYS package and MATLAB programming. The output of the exergy analysis of the processes including the physical, chemical, and total exergy of the streams are listed in Table 4.
The equations used for the exergy analysis of the integrated structure are presented in Appendix A. The input and output exergy, the exergy degradation, and the efficiency of the implemented devices in the combined process for the tri-generation of liquid carbon dioxide, electricity, and medium pressure steam are presented in Table 5.
For the exergy analysis, the performance of each part of the equipment should be evaluated from the irreversibility and exergy yield points of view. The simulation results show that the exergy yield of the valves is lower than other devices, while they had less irreversibility. Among the expansion valves, the V2 valve with 67.43% efficiency has the highest exergy yield and the V1 valve with 53.17% efficiency has the lowest exergy yield. Among the heat exchangers, the HX5 exchanger with 80.85% efficiency has the lowest exergy efficiency and the HX4 exchanger with 96.95% efficiency has the highest exergy efficiency. Among the heat exchangers, the HX6 exchanger has the highest exergy destruction with a magnitude of 8.258 MW, and the HX3 exchanger has the lowest exergy destruction with a magnitude of 0.2251 MW. Among the compressors of the developed hybrid process, the C3 compressor with 78.71% efficiency has the lowest exergy efficiency and the C2 compressor with 85.87% efficiency has the highest exergy yield. Moreover, among the compressors, the C2 compressor has the highest exergy destruction with 4.651 MW and the C3 compressor has the lowest exergy destruction with 1.818 MW. The contribution of exergy destruction of each piece of equipment in the combined process is depicted in Figure 6 separately. The results show that compared to other equipment, the MCFC, heat exchangers, and compressors utilized in the combined process have the most exergy destruction with a share of 37.72%, 18.71%, and 12.26%, respectively. The total exergy efficiency and exergy destruction of the entire combined process is 73.18% and 102.7 MW, respectively.

Sensitivity Analysis
With the sensitivity study, important and influential parameters of the combined system of liquid carbon dioxide, electricity, and medium pressure steam production using a molten carbonate fuel cell, a Linde-Hampson double-effect refrigeration system, and a HSRG were identified and investigated and their behavior in regard to these changes was studied. Figure 7 shows the variation of the refrigeration cycle performance coefficient and the refrigeration cooling duty with the reduction in the outlet pressure of the V1 expansion valve. The results of the sensitivity analysis show that when the outlet pressure of the V1 expansion valve was increased from 39 bar to 45 bar, the amount of steam coming out of the D4 flash drum (in other words, the return flow of carbon dioxide) decreased, and as a result, the refrigeration cooling duty increased by 0.7165% and the refrigeration cycle performance coefficient was reduced by 9.654%.

Sensitivity Analysis
With the sensitivity study, important and influential parameters of the combined system of liquid carbon dioxide, electricity, and medium pressure steam production using a molten carbonate fuel cell, a Linde-Hampson double-effect refrigeration system, and a HSRG were identified and investigated and their behavior in regard to these changes was studied. Figure 7 shows the variation of the refrigeration cycle performance coefficient and the refrigeration cooling duty with the reduction in the outlet pressure of the V1 expansion valve. The results of the sensitivity analysis show that when the outlet pressure of the V1 expansion valve was increased from 39 bar to 45 bar, the amount of steam coming out of the D4 flash drum (in other words, the return flow of carbon dioxide) decreased, and as a result, the refrigeration cooling duty increased by 0.7165% and the refrigeration cycle performance coefficient was reduced by 9.654%. The changes in the refrigeration load and the refrigeration energy with the reduction in the outlet pressure of the V1 expansion valve are presented in Figure 8. The results of the sensitivity analysis show that when the outlet pressure of the V1 expansion valve was increased from 39 bar to 45 bar, the amount of refrigeration energy in the refrigeration cycle increased by 2.151%. Moreover, by increasing the outlet pressure of the V1 expansion valve to 41 bar, the amount of the refrigeration load decreased and then increased to 45 bar.  The changes in the refrigeration load and the refrigeration energy with the reduction in the outlet pressure of the V1 expansion valve are presented in Figure 8. The results of the sensitivity analysis show that when the outlet pressure of the V1 expansion valve was increased from 39 bar to 45 bar, the amount of refrigeration energy in the refrigeration cycle increased by 2.151%. Moreover, by increasing the outlet pressure of the V1 expansion valve to 41 bar, the amount of the refrigeration load decreased and then increased to 45 bar.  Figure 9 shows the changes in electricity produced by the fuel cell and electrical efficiency with the reduction in fuel input to the fuel cell. The output illustrates that when fuel input to the fuel cell increased, the electricity generated by the fuel cell and electrical efficiency increased from 89.68 MW and 29.06% to 134.5 MW and 37.75%, respectively.  Figure 9 shows the changes in electricity produced by the fuel cell and electrical efficiency with the reduction in fuel input to the fuel cell. The output illustrates that when fuel input to the fuel cell increased, the electricity generated by the fuel cell and electrical efficiency increased from 89.68 MW and 29.06% to 134.5 MW and 37.75%, respectively. The variation of the net overall thermal efficiency and the overall exergy efficiency with respect to the reduction in fuel input to the fuel cell is depicted in Figure 10. The output illustrates that when the fuel input to the fuel cell increased, the net overall thermal efficiency in the integrated structure increased by 18.09%, and the overall exergy efficiency decreased by 12.95%.

Conclusions
The storage of liquid carbon dioxide is a method of storing energy for long periods which is portable and can be used during peak hours. Providing the required electricity for the liquefaction cycle during off-peak hours is a concern of many researchers in this field. Using molten carbonate fuel cells as a source of high efficiency and low pollution power generation is one of the ways to provide electricity and heat to store energy during peak hours. In this research, a combined process for the production of electricity, medium pressure steam, and energy storage in the form of liquid carbon dioxide was developed. For co-production, a combination of a molten carbonate fuel cell, a dual pressure Linde-Hampson liquefaction plant and a heat recovery steam generator was used. The results of the energy and exergy analyses were as follows.
In this integrated structure, 4.202 kg/s of natural gas and 144 kg/s of carbon dioxide with impurities were the main input feed. The primary products of this integrated structure were 65.53 MW of electricity, 27.80 kg/s of medium pressure steam, and 142.9 kg/s of liquid CO 2 . The thermal efficiency of the fuel cell, the total thermal efficiency of the entire integrated structure, and the performance coefficient of the liquefaction plant were 46.94%, 82.21%, and 1.866, respectively. The results of the CO 2 liquefaction plant validation showed that the performance coefficient of the developed refrigeration plant and the reference were 1.87 and 1.83, respectively. The low relative error of the compared parameters depends on the version of HYSYS software used for the simulation.
The output of the exergy investigation showed that fuel cells, heat exchangers, and compressors had the most exergy destruction with values of of 37.72%, 18.71%, and 12.26, respectively. Moreover, the heat exchangers had the highest exergy yield and the valves had the lowest exergy yield. The exergy efficiency and destruction of equipment results showed that the equipment with high exergy destruction shares, such as heat exchangers, also had high exergy efficiencies. Moreover, the equipment with low exergy efficiencies such as the expansion valves had a low exergy destruction share. Considering the above showed that the developed simulation had a high efficiency of 73.18%.
The results obtained from the sensitivity analysis showed that a 50% increase in fuel input to the fuel cell led to an increase in the electrical efficiency of the fuel cell and an increase in the thermal efficiency of the entire integrated structure by 37.75% and 18.09%, respectively. Moreover, the exergy yield of the entire hybrid process was reduced by 12.95%. Furthermore, by increasing the outlet pressure of the V1 expansion valve from 39 bar to 45 bar, the total cooling duty rate increased from 116.4 MW to 117 MW, and the refrigeration cycle performance coefficient decreased from 1.866 to 1.770.
As a suggestion for future works, the development of the power generation cycle in peak hours for network use could be mentioned. Advanced exergy and advanced exergoeconomic analyses could also be used to investigate the developed co-generation cycle.

Conflicts of Interest:
The author declares no conflict of interest. Table A1. Equations needed for the exergy analysis of the equipment employed in the developed integrated process [40][41][42][43].

Equipment Exergy Destruction Exergy Efficiency
Heat Exchangers m∆h cold