Modeling Analysis of a Polygeneration Plant Using a CeO2/Ce2O3 Chemical Looping

In the current context of complexity between climate change, environmental sustainability, resource scarcity, and geopolitical aspects of energy resources, a polygenerative system with a circular approach is considered to generate energy (thermal, electrical, and fuel), contributing to the control of CO2 emissions. A plant for the multiple productions of electrical energy, thermal heat, DME, syngas, and methanol is discussed and analyzed, integrating a chemical cycle for CO2/H2O splitting driven using concentrated solar energy and biomethane. Two-stage chemical looping is the central part of the plant, operating with the CeO2/Ce2O3 redox couple and operating at 1.2 bar and 900 °C. The system is coupled to biomethane reforming. The chemical loop generates fuel for the plant’s secondary units: a DME synthesis and distillation unit and a solid oxide fuel cell (SOFC). The DME synthesis and distillation unit are integrated with a biomethane reforming reactor powered by concentrated solar energy to produce syngas at 800 °C. The technical feasibility in terms of performance is presented in this paper, both with and without solar irradiation, with the following results, respectively: overall efficiencies of 62.56% and 59.08%, electricity production of 6.17 MWe and 28.96 MWe, and heat production of 111.97 MWt and 35.82 MWt. The fuel production, which occurs only at high irradiance, is 0.71 kg/s methanol, 6.18 kg/s DME, and 19.68 kg/s for the syngas. The increase in plant productivity is studied by decoupling the operation of the chemical looping with a biomethane reformer from intermittent solar energy using the heat from the SOFC unit.


Introduction
Greenhouse gas emissions, climate change, and carbon scarcity are serious problems that require a strong focus. The sixth assessment report by the IPCC has been published, in which the role of anthropogenic activities is updated [1,2]. To reduce this impact on the environment, it was decided to limit the increase in global average temperature to 1.5 • C above pre-industrial levels [3]. The power generation sector plays a key role in reducing greenhouse gas emissions [4,5]. The use of solar energy in power generation could be a key tool to limit the use of fossil fuels. Solar radiation is highly available, globally distributed, clean, and renewable. However, the intermittent nature and recurring mismatch between consumers and direct supply are the main factors limiting its expansion. These problems can be overcome by non-fossil solar fuels [6,7], which store solar energy [8]. In this paper, solar energy is stored using a mixture of carbon monoxide and hydrogen. This gas mixture can be used as a direct fuel in a gas turbine [9,10] and an internal combustion engine [11,12] or in a fuel cell to generate electricity [12]. Furthermore, the conversion of syngas into fuels and other chemicals makes it possible to produce carbon-neutral materials [13][14][15]. The production of dimethyl ether from syngas is becoming an increasingly attractive route. The scientific community regards DME as the fuel that could eliminate the dependence on fossil fuels [16]. This is due to its physical properties, which are very similar to those 1.
Reduction reaction (1): the concentrated solar energy provides thermal heat. The metal oxide MOox undergoes thermal reduction and releases oxygen (MOred is obtained).

2.
Oxidation reaction (2): the reduced metal oxide MOred reacts with water and carbon dioxide to achieve both the reoxidation of metal and the reduction of H 2 O and CO 2 with the subsequent production of syngas. becoming an increasingly attractive route. The scientific community regards DME as the fuel that could eliminate the dependence on fossil fuels [16]. This is due to its physical properties, which are very similar to those of liquefied petroleum gas (LPG, e.g., propane and butane), with a low environmental and health impact [17]. Moreover, DME distillation also produces methanol (CH3OH) after dimethyl ether synthesis from syngas [18]. The main methods for producing solar syngas are photobiological or photochemical conversion, thermochemical conversion, and electrochemical conversion. The first route occurs at relatively low temperatures by directly using solar energy from photons to power the process [19,20]. Electrochemical conversion processes consist of the coelectrolysis of CO2 and H2O in an electrolyzer that supplies electricity from photovoltaic systems or solar thermal power plants [21,22]. Thermochemical conversion processes use concentrating solar collectors (CSP) to concentrate solar radiation to a point where a receiver reactor is located [23,24]. The simplest thermochemical pathway for syngas production would be the direct one-step solar thermal dissociation of water and carbon dioxide (thermolysis). However, this route is not practical, because an extremely high temperature of more than 2200 °C is required and cannot be achieved, due to the limitations imposed by the reactor materials, the thermal losses, and the possible products recombination to avoid explosive mixtures [25][26][27][28]. One way to lower the reduction temperature is to use redox cycles [20,29]. Two phases characterize the basic and simplest thermochemical cycles [30][31][32], with overall efficiencies of around 40-50% [33,34]. Considering a metal oxide (MO) with different oxidation states, the two steps are shown in Figure 1: 1. Reduction reaction (1): the concentrated solar energy provides thermal heat. The metal oxide MOox undergoes thermal reduction and releases oxygen (MOred is obtained). 2. Oxidation reaction (2): the reduced metal oxide MOred reacts with water and carbon dioxide to achieve both the reoxidation of metal and the reduction of H2O and CO2 with the subsequent production of syngas.
(Reduction reaction) 2CeO 2 + CH 4 → CO + 2H 2  Various metal-redox pairs have been investigated for the TSC processes. Ceria is currently considered the most attractive metal-redox pair for various reasons [35][36][37][38][39]: high crystallographic stability during extended thermal cycling (non-volatile cycle), high oxygen release and storage capacities (>100 mmol min −1 g −1 ), fast oxygen exchange rates, reversible shift between Ce 4+ and Ce 3+ oxidation states, and fast kinetics during thermochemical cycles as compared to other non-volatile metal oxides. However, the reduction reaction of pure ceria occurs only at very high temperatures (above 1400 °C) and very low pressures (between 100 and 200 mbar) [38]. The introduction of biomethane as a reducing agent (see Figure 2) could solve these problems [40][41][42][43] while preserving sustainability criteria, allowing the realization of an isothermal and isobaric cycle according to reactions (1) and (2).
Polygeneration systems could represent an efficient and cost-effective solution for the integration of solar thermochemical conversion processes [47]. These systems have gained increasing attention in recent years due to their high overall efficiency and the possibility to use waste heat streams to generate useful heat or to obtain additional useful power [33,48]. As a result, these solutions generate multiple products such as electricity, cooling, heat, freshwater, and chemicals from one or more input streams in a single system. Several examples of polygeneration systems integrated with solar thermochemical conversion units can be found in the literature. Kaniyal et al. [49] performed a polygeneration system analysis based on hybrid solar gasification of coal at atmospheric pressure to produce liquid fuels and electricity. To solve the problem of intermittent solar energy, they proposed the pressurized storage of upgraded syngas and oxygen that can be used during periods of low solar radiation. Liu et al. [50] conducted a Other thermochemical processes to produce syngas could be [25]: steam or dry natural gas/biomethane reforming [44], and biomass or coal gasification [45]. Biomethane could be produced by different separation mechanisms for upgrading the biogas extracted from the organic fraction of MSW (e.g., pressure swing adsorption, cryogenic separation, and absorption techniques [46]).
Polygeneration systems could represent an efficient and cost-effective solution for the integration of solar thermochemical conversion processes [47]. These systems have gained increasing attention in recent years due to their high overall efficiency and the possibility to use waste heat streams to generate useful heat or to obtain additional useful power [33,48]. As a result, these solutions generate multiple products such as electricity, cooling, heat, freshwater, and chemicals from one or more input streams in a single system. Several examples of polygeneration systems integrated with solar thermochemical conversion units can be found in the literature. Kaniyal et al. [49] performed a polygeneration system analysis based on hybrid solar gasification of coal at atmospheric pressure to produce liquid fuels and electricity. To solve the problem of intermittent solar energy, they proposed the pressurized storage of upgraded syngas and oxygen that can be used during periods of low solar radiation. Liu et al. [50] conducted a thermodynamic study on a lignite polygeneration system powered by solar energy. They demonstrated that the use of solar energy in lignite drying pyrolysis and gasification can save 63.84% of lignite energy compared to conventional methods. Kong et al. studied a polygeneration plant based on an isothermal thermochemical cycle fed with concentrated solar energy and methane to produce electricity and methanol simultaneously. They found that the fossil fuel consumption to produce one unit mass of methanol was about 22 GJ/ton with an optimum efficiency of more than 44%. Bai et al. [51] carried out a techno-economic analysis of a polygeneration system based on solar thermal gasification of biomass to produce methanol and electricity. They concluded that the energy and exergy efficiencies could reach values of 56.09% and 54.86%, respectively.
All literature studies are characterized by the production of CO 2 streams, making the system not exactly carbon neutral. These gas streams could be emitted into the atmosphere without proper and accurate carbon capture and sequestration processes.
As far as we know, there are no studies on zero-emissions polygeneration plants at an industrial scale. In this paper, a polygeneration system with the solar-driven thermochemical cycle is simulated, analyzed, and discussed.

Case Study: Polygeneration Plant
The plant is built with four units:

1.
Chemical looping unit. This section is the driving force of the plant and produces the fuel for the other units. The CeO 2 /Ce 2 O 3 redox couple associated with biomethane reforming is chosen ( Figure 2). The cycle is unpressurized and isothermal at 900 • C and 1.2 bar. The reduction reaction takes place within a solar receiver via a solar tower, while the oxidation reaction takes place in a non-solar reactor. These two reactors operate simultaneously when solar energy is available. Under steady-state conditions, the reduction reactor is fed with 0.59 kmol/s of ceria with a particle size of 0.4-1 mm; 2.
Fixed SOFC unit. The secondary unit is fed by the syngas produced in the CL reduction reactor. Electricity is generated in this unit; 3. DME synthesis and distillation unit integrated with biomethane reforming reactor fed by concentrated solar energy for CCU unit [52]. This is another secondary unit fed with syngas. The syngas is produced in the oxidation reactor of CL. DME, methanol, and syngas are produced from this unit; 4.
Steam generation plant. There are five steam generators with heat recovery in the plant that use the hot flows of the system to generate steam. This steam is used to heat other cold process flows and to generate thermal energy to supply other consumers (e.g., a district heating network).
The chemical cycle can operate properly at sufficient irradiance to heat the reduction reactor at 900 • C, while at lower irradiance levels, the chemical cycle turns off. A comparison is made with a similar non-solar polygeneration system published in the literature [53]. An attempt to increase the system's productivity by decoupling its operation from solar energy fluctuations is made. The results of the chemical cycle integrated into the system are compared with experimental [54] and modeling studies [41,55].
The proposed plant produces thermal energy, electricity, syngas, methanol, and dimethyl ether (DME), in two operating conditions: (1) Operation on a day with a clear sky and irradiance sufficient for the CL operation; (2) Operation on a day with low solar irradiance.
In the first operating condition (see Figure 3), all components of the plant can work properly. The reduction reactor of the CL is oversized compared to the SOFC. An excess of syngas is obtained from the reduction reactor. This gas mixture is stored and used during periods of low solar irradiation. This storage system is fundamental to couple the solar syngas with the correct operation of the SOFC system. This system is very sensitive to load variations and transients due to its high operating temperature.
The uniqueness of the SOFC system's operation is provided in the second operating condition plant. This system is fed by properly stored syngas (see Figure 4). When solar energy is not present, the chemical cycle and consequently the synthesis of DME, the distillation unit, and the reforming reactor do not come into operation. Therefore, the anode fluxes of the SOFC system are stored for reuse when the CL system is operational. solar syngas with the correct operation of the SOFC system. This system is very sensitive to load variations and transients due to its high operating temperature. The uniqueness of the SOFC system's operation is provided in the second operating condition plant. This system is fed by properly stored syngas (see Figure 4). When solar energy is not present, the chemical cycle and consequently the synthesis of DME, the distillation unit, and the reforming reactor do not come into operation. Therefore, the anode fluxes of the SOFC system are stored for reuse when the CL system is operational.  These different operating conditions are defined by considering the seasonal average daily temperatures of the reduction reactor. The experimental data for the modeling analysis are taken from the dish facility placed at the Turin-Energy Center [56]. Looking at the seasonal daily average temperature of the dish system, it is found that the temperature is higher/equal to 900 °C (see Figure 5  These different operating conditions are defined by considering the seasonal average daily temperatures of the reduction reactor. The experimental data for the modeling analysis are taken from the dish facility placed at the Turin-Energy Center [56]. Looking at the seasonal daily average temperature of the dish system, it is found that the temperature is higher/equal to 900 • C (see Figure 5):  These different operating conditions are defined by considering the seasonal average daily temperatures of the reduction reactor. The experimental data for the modeling analysis are taken from the dish facility placed at the Turin-Energy Center [56]. Looking at the seasonal daily average temperature of the dish system, it is found that the temperature is higher/equal to 900 °C (see Figure 5  In summary, the chemical cycle can operate around 21% of the year. The SOFC system is designed to operate continuously. In summary, the chemical cycle can operate around 21% of the year. The SOFC system is designed to operate continuously.

Modeling and Plant Simulation in Aspen Plus
In this study, the simulation tool used is Aspen Plus ® from AspenTech (v10, Bedford, MA, USA). In both plant schemes (Figures 3 and 4), the plant is modeled mainly with built-in components and assumes chemical equilibrium in each component. The DME reactor is characterized by a kinetic approach. The main plant components are: -Turbines, valves, coolers, distributors, mixers, heaters, compressors, heat recovery steam generators, flash units, and cyclones; -RGIBBS reactor units are used to model the chemical looping oxidation and reduction reactors, the post-combustion unit, the SOFC anode, and the biomethane reforming reactor; -Separator (coupled with a heat exchanger) for the cathode of the SOFC; -RADFRAC columns used to model the DME distillation unit; -RPLUG reactor coupled to the kinetic model Langmuir-Hinshelwood-Hougen-Watson (LHHW) is used to simulate the DME synthesis reactor with a catalytic behavior. In this unit, it is used for the thermodynamic properties of the Soave-Redlich-Kwong (SRK) equation of state (EOS). A similar approach was followed by Graaf et al. [58], where the SRK equation of state was used to model the chemical equilibrium of the methanol and water gas shift (WGS) reaction. This model is applied to binary components [59].
In the model, the flows used are conventional (H 2 O, CO 2 , H 2 , CO, CH 4 , N 2 , O 2 , CH 3 OH, and CH 3 OCH 3 ) and fixed flows (oxygen carriers of chemical looping, CeO 2 and Ce 2 O 3 ). The Barin equation is used for the evaluation of the fixed flow [60]. For conventional flows, the Peng-Robinson-Boston-Mathias (PR-BM) properties method is used. This approach is widely used in the literature for the treatment of hydrocarbons, such as refining in petrochemical processes [61,62]. Table 1 lists the main assumptions used. Table 1. Main assumptions and hypotheses adopted in the process simulation.

Oxidation and reduction reactors
Model: RGIBBS, no heat losses

Compressors, pumps, and turbines
η is,comp = 0.9 η mech,comp = 0.98 η is, tur = 0.9 η mech,tur = 0.98 η is,pump = 0.9 η driver,pump = 0.9 The constant-temperature and -pressure operating conditions for the CL process increase the efficiency by reducing losses associated with heat flows and high system pressures. Concerning the reduction reaction, the operating temperature and the amount of methane to be fed are chosen to consider the following three requirements (see Figure 6):

SOFC
- The temperature should be compatible with a solar system; - The complete reduction of the oxygen carrier to Ce 2 O 3 must take place with a limited amount of biomethane; - The desired syngas concentration has a molar ratio of H 2 /CO around 2. This molar ratio is more suitable to feed fuel cell systems. H 2 diffuses quicker in the anode than CO and the diffusion overvoltage is reduced [63]. The flow mixture between the biomethane and CeO 2 ratio varies directly with syngas production. Increasing the biomethane content increases the CO flow rate and, thus, decreases the H 2 /CO ratio.
The H 2 O/CO 2 molar ratio is ≈1.13 for the oxidation reaction to obtain: - The complete re-oxidation of Ce 2 O 3 to CeO 2 ; -An H 2 /CO molar ratio around 1. This is an optimal value for the downstream DME production in the proper reactor (see Figure 7).   The inlet streams of the reduction reactor (RED-CL) are constituted by biomethane (METH-1) and CeO2 (CER-OXY). Particle size is defined from the particle distribution reported in Supplementary Table S1 [64]. At the outlet of this reactor, the flow PROD-1 is separated in cyclone CYC-1. Two different streams are produced. The first stream (Ce2O3 (CER-RED)) is sent to the oxidation reactor, while the second stream is the syngas fraction (SYN-SOF1). This stream is compressed and cooled to be stored in the AISI316L tank [65]. This fraction is stored at 10 bar and 800 °C, following the mechanical and physical properties of AISI316L. However, in Aspen Plus, there is no model for the storage system simulation. Thus, downstream COMP-SYR and SRG1, there is a splitter (SPL-SYNS),    The inlet streams of the reduction reactor (RED-CL) are constituted by biomethane (METH-1) and CeO2 (CER-OXY). Particle size is defined from the particle distribution reported in Supplementary Table S1 [64]. At the outlet of this reactor, the flow PROD-1 is separated in cyclone CYC-1. Two different streams are produced. The first stream (Ce2O3 (CER-RED)) is sent to the oxidation reactor, while the second stream is the syngas fraction (SYN-SOF1). This stream is compressed and cooled to be stored in the AISI316L tank [65]. This fraction is stored at 10 bar and 800 °C, following the mechanical and physical properties of AISI316L. However, in Aspen Plus, there is no model for the storage system simulation. Thus, downstream COMP-SYR and SRG1, there is a splitter (SPL-SYNS), The inlet streams of the reduction reactor (RED-CL) are constituted by biomethane (METH-1) and CeO 2 (CER-OXY). Particle size is defined from the particle distribution reported in Supplementary Table S1 [64]. At the outlet of this reactor, the flow PROD-1 is separated in cyclone CYC-1. Two different streams are produced. The first stream (Ce 2 O 3 (CER-RED)) is sent to the oxidation reactor, while the second stream is the syngas fraction (SYN-SOF1). This stream is compressed and cooled to be stored in the AISI316L tank [65]. This fraction is stored at 10 bar and 800 • C, following the mechanical and physical properties of AISI316L. However, in Aspen Plus, there is no model for the storage system simulation. Thus, downstream COMP-SYR and SRG1, there is a splitter (SPL-SYNS), which separates the flow entering the anode of the SOFC (SYN-SOF4) from the flow that should remain in the storage (SYN-S-ST), and a valve (VAL-SY) that reduces the pressure of the portion of syngas sent to the SOFC at 5 bar (SYN-SOF5). To sum up, for a continuous operation of the SOFC, the SYN-SOF1 mole flow is spread over the whole year and the molar flow stream SYN-SOF4 is obtained according to Equation (5).
kmol min operating miuntes of the CL operating minutes of the SOFC : 60 s min (5) Going back to the chemical looping unit, the oxidation reactor (OXY-CL) is fed by the reduced cerium oxide (CER-RED) and a mixture of H 2 O and CO 2 coming from the anodic exhausts of the SOFC. The product of the oxidation reactor (PROD-2) goes to a cyclone CYC-2 to separate the solid fraction (CeO 2 ) and the gaseous fraction (SYN-DME1). The solid flow is sent to the reduction reactor to close the cycle (CER-OX1) after a re-integration of the lost ceria (CER-REG2) for the nonunitary efficiency of the cyclones (Supplementary Table S2). The gaseous flow is sent to the DME production section of the plant.

Solid Oxide Fuel Cell (Secondary Device)
To produce electricity, a fuel cell is used to consider its higher efficiency and fuel flexibility compared to a thermal machine [66]. Particularly, an anode-supported planar Solid Oxide Fuel Cell is chosen for its advantages compared to the other fuel cells in the industry [67][68][69][70]. The SOFC is simulated in Aspen Plus for both situations when the CL is in the ON-state ( Figure 9) and OFF-state (Supplementary Figure S1). To produce electricity, a fuel cell is used to consider its higher efficiency and fuel flexibility compared to a thermal machine [66]. Particularly, an anode-supported planar Solid Oxide Fuel Cell is chosen for its advantages compared to the other fuel cells in the industry [67][68][69][70]. The SOFC is simulated in Aspen Plus for both situations when the CL is in the ON-state ( Figure 9) and OFF-state (Supplementary figure S1). The operating conditions of the SOFC system are: 1) Temperature = 850 °C, which is compatible with syngas storage; 2) Pressure = 5 bar, to improve the cell's performance [71]. The SOFC operating temperature could be lowered using new and advanced materials. It is selected as a condition suitable for the bloom energy generator [72].
As specified in Section 2.2.1, the syngas mole flow entering the anode of the SOFC is SYN-SOF5; another inlet stream in the anode is given by the oxygen ions (OXY-AN). The molar flow rate of this stream is equal to the stoichiometric molar flow rate of the oxygen O2 needed to generate in the cathode an amount of ions O 2-sufficient to oxidize the fuel. It is calculated according to the Faraday law (4): &#x ˙0078; The operating conditions of the SOFC system are: (1) Temperature = 850 • C, which is compatible with syngas storage; (2) Pressure = 5 bar, to improve the cell's performance [71].
The SOFC operating temperature could be lowered using new and advanced materials. It is selected as a condition suitable for the bloom energy generator [72].
As specified in Section 2.2.1, the syngas mole flow entering the anode of the SOFC is . n SYN-SOF5 ; another inlet stream in the anode is given by the oxygen ions (OXY-AN). The molar flow rate of this stream is equal to the stoichiometric molar flow rate of the oxygen O 2 needed to generate in the cathode an amount of ions O 2− sufficient to oxidize the fuel. It is calculated according to the Faraday law (4): .
where C tot = current produced by the cell and Z O 2 = 4, so . n oxy-an = 148.42 mol s . This molar flow rate of oxygen is set equal to the molar flow rate of the stream OXY-CATH exiting from the component SOF-CA-S to simulate the conduction of the ions from the anode to the cathode in the electrolytic layer of the SOFC. The cathode is supplied by air that enters the plant at ambient conditions to then be compressed and heated to 5 bar and 600 • C. The molar flow rate of the air stream is calculated iteratively to be heated from 600 • C to 850 • C, absorbing the entire heat rejected by the SOFC [73]: . n AIR-CA-3 = 3.044 kmol s . Once the air stream (AIR-CA-3) enters the cathode, it first encounters a heater that simulates the heating of the air and then a separator in which the streams OXY-CATH and CATH-EXH1 are produced. The latter represents the cathodic exhausts. They are expanded and then split into two streams. One stream (CAT-EXH4) feeds an HRSG (SRG4) to produce steam before the external release (CAT-EXH5); the other stream (CAT-EXH3) feeds the post-combustion unit. Considering a cell voltage equal to around 0.8 V, the Solid Oxide Fuel Cell can produce an electric power equal to 45.82 MW e . Using the Bloom Energy Servers [72] with an electric power production of 200 kW, the number of necessary modules is 230.

DME Production (Secondary Device) DME Synthesis
The single-step pathway with a dual catalyst and the co-production of methanol and DME is chosen for the DME synthesis. This process is characterized by a significant increase in the total methanol yield compared to the two-step pathway with methanol and DME produced in two different reactors [74]. Three main reactions occur: syngas conversion to methanol (reaction (7)), water gas shift (reaction (8)), and methanol dehydration to DME (reaction (9)).
CO + H 2 O → CO 2 + H 2 ∆H 0 = −41.2 kJ/mol (8) The global reaction is: According to the Le Châtelier principle [42], the global reaction is favored at high pressure and low temperature. The chosen catalytic reactor is a multi-tube fixed-bed reactor that is kept at a constant temperature of 250 • C and at 50 bar [75]. The bi-functional catalyst Cu/ZnO/Al 2 O 3 : γ-Al 2 O 3 is selected considering a loading ratio of 1:2 from the literature. The syngas produced in the oxidation reactor of the chemical looping (SYN-DME2), before being sent to the DME synthesis reactor (see Figure 10), is: • Cooled in the heat recovery steam generator (SRG2), producing steam; • Separated from water (in FLA-SYDM); • Integrated with additional syngas (SYN-RE3R) produced in the biomethane reforming reactor and compressed to the operating pressure.
literature. The syngas produced in the oxidation reactor of the chemical looping (SYN-DME2), before being sent to the DME synthesis reactor (see Figure 10), is:  Cooled in the heat recovery steam generator (SRG2), producing steam;  Separated from water (in FLA-SYDM);  Integrated with additional syngas (SYN-RE3R) produced in the biomethane reforming reactor and compressed to the operating pressure. Figure 10. Pre-treatment of the syngas sent to the DME reactor in Aspen Plus. Figure 10. Pre-treatment of the syngas sent to the DME reactor in Aspen Plus.
The DME reactor is simulated in Aspen plus ( Figure 11) and solved with the SRK-EOS property method. The parameters for the modeling of this component are listed in The DME reactor is simulated in Aspen plus ( Figure 11) and solved with the SRK-EOS property method. The parameters for the modeling of this component are listed in Table 2 [53]. Figure 11. Simulation of the DME reactor in Aspen Plus. A Langmuir-Hinshelwood-Hougen-Watson (LHHW) kinetic model is implemented considering the three simultaneous reactions defined above ((7), (8), and (9)). The expressions of the rates for CO2 hydrogenation, WGS, and methanol dehydration are calculated according to Equations (11)-(13), respectively [74,76]. These reaction rates are expressed in  Figure 11. Simulation of the DME reactor in Aspen Plus. A Langmuir-Hinshelwood-Hougen-Watson (LHHW) kinetic model is implemented considering the three simultaneous reactions defined above ((7), (8), and (9)). The expressions of the rates for CO 2 hydrogenation, WGS, and methanol dehydration are calculated according to Equations (11)-(13), respectively [74,76]. These reaction rates are expressed in kmol kg cat ·s : with: p, gas partial pressure (Pa); -C, concentration (kmol/m 3 ).

Distillation Unit
The DME synthesized in the reactor (PROD-3) is characterized by a high amount of impurities; thus, a separation process and a distillation process are necessary to obtain pure DME. Before the distillation columns, there is a vapor-liquid separation unit (V-L-SEP) at −45 • C and 10 bar, as shown in Figure 12. The output streams of this component are incondensable gases (H2COCO2) and a liquid stream (LIQ1). The first stream (H 2 COCO 2 ) is pre-heated to 400 • C and expanded. The resulting stream is sent to a post-combustion chamber (POST-COM) (see Figure 12) to produce H2OCO2-1, which is cooled in a heat recovery steam generator (SRG3) to be separated in a flash unit (SEP-RE) (see Figure 13). The obtained water at the bottom of the reactor is used to produce steam in counter-flow with the post-combustion exhausts in the previously mentioned SRG3, while the gaseous stream feeds the biomethane reforming reactor (METH-REF) after being mixed with other streams of CO 2 . These additional streams come from the portion of the anodic exhausts of the SOFC (CO2-RS-3) stored when the chemical looping is in OFF-state and not used in the oxidation reactor of the CL itself throughout the year. Another stream that enters the reforming reactor is biomethane (METH-2), taken from the pipeline. The aim of the biomethane reforming reactor (METH-REF) is to produce syngas following this reaction: This reaction is favored at low pressures and high temperatures. In the present study, it takes place at 800 • C and 1 bar. It can be fed by concentrated solar thermal energy in an additional receiver-reactor placed close to the reduction reactor of the CL. The syngas stream obtained from this reactor (SYN1) is cooled in the heat recovery steam generator (SRG5) and split (in SPL-SYN) to be partially recirculated to the DME synthesis reactor (SYN-RE3R) and partially sent to a syngas duct (SYN-DUCT). This last stream is an additional output of the overall system.
This reaction is favored at low pressures and high temperatures. In the present st it takes place at 800 °C and 1 bar. It can be fed by concentrated solar thermal energy additional receiver-reactor placed close to the reduction reactor of the CL. The sy stream obtained from this reactor (SYN1) is cooled in the heat recovery steam gene (SRG5) and split (in SPL-SYN) to be partially recirculated to the DME synthesis re (SYN-RE3R) and partially sent to a syngas duct (SYN-DUCT). This last stream is an a tional output of the overall system.  The second stream from V-L-SEP is a liquid stream mainly constituted of dissolved CO2, DME, and CH3OH (LIQ1). This stream is further treated in three different distillation columns ( Figure 14): 1. Column for CO2 separation (DIST-CO2); 2. Column for DME production (DIST-DME); 3. Column for methanol separation from water (DIST-MET). The water stream (H2O-DST1) is mixed with other water streams for steam production in SRG5. A valve and a heat exchanger are placed before each column to adjust the pressure to the optimal value and to have 50% of vapor in the column inlet stream [75]. The number of stages used in the distillation columns is estimated by increasing them until a certain change in composition is detected. In Table 3, the data inserted in the distillation columns are listed; while in Table 4 are listed the composition of thermodynamic properties of selected streams when the chemical looping is in ON-State.   The second stream from V-L-SEP is a liquid stream mainly constituted of dissolved CO 2 , DME, and CH 3 OH (LIQ1). This stream is further treated in three different distillation columns ( Figure 14):
Column for methanol separation from water (DIST-MET). The water stream (H2O-DST1) is mixed with other water streams for steam production in SRG5. The second stream from V-L-SEP is a liquid stream mainly constituted of dissolved CO2, DME, and CH3OH (LIQ1). This stream is further treated in three different distillation columns ( Figure 14): 1. Column for CO2 separation (DIST-CO2); 2. Column for DME production (DIST-DME); 3. Column for methanol separation from water (DIST-MET). The water stream (H2O-DST1) is mixed with other water streams for steam production in SRG5. A valve and a heat exchanger are placed before each column to adjust the pressure to the optimal value and to have 50% of vapor in the column inlet stream [75]. The number of stages used in the distillation columns is estimated by increasing them until a certain change in composition is detected. In Table 3, the data inserted in the distillation columns are listed; while in Table 4 are listed the composition of thermodynamic properties of selected streams when the chemical looping is in ON-State.   A valve and a heat exchanger are placed before each column to adjust the pressure to the optimal value and to have 50% of vapor in the column inlet stream [75]. The number of stages used in the distillation columns is estimated by increasing them until a certain change in composition is detected. In Table 3, the data inserted in the distillation columns are listed; while in Table 4 are listed the composition of thermodynamic properties of selected streams when the chemical looping is in ON-State.

Steam Production
In the analyzed plant, there are five heat recovery steam generators (SRG1, SRG2, SRG3, SRG4, and SRG5) with the primary purpose of cooling certain streams of the plant without wasting their high-temperature heat. In certain generators (SRG1 and SRG2), the water needs are completely satisfied through the recirculation of water produced in the system itself from different processes. In the others, the water is re-integrated from external sources. In Table 5, the steam streams produced in these HRSGs are listed.

Solar Tower and Heliostats
The reduction reactor of the chemical looping and the biomethane reforming reactor work simultaneously; they can be supplied by concentrated solar energy. As a result of the high temperatures and the overall thermal requirements (W th,tot~2 94 MW t ), it is considered to use a solar tower system [77,78]. This system could be installed in Turin, with two receiver-reactors where the reduction and the oxidation take place. These receiverreactors are similar to the receiver-reactor of the solar dish located on the Energy Center roof, both in terms of irradiance and temperature values. Furthermore, these data are available: The experimental daily values of the direct normal irradiance in the neighborhood of the Energy Center, measured by the Politecnico meteorological station;

2.
The experimental seasonal daily average temperatures of the receiver-reactor at the focus of the dish system installed on the Energy Center roof.
The yearly average direct normal irradiance (DNI) hitting the receiver-reactors and making the CL work properly is calculated to be around 828 W m 2 . Consequently, considering the optical efficiency of the solar field equal to 80%, its extension area should be around 44.27 hectares. To cover this area, 3170 heliostats with an aperture area of 140 m 2 are necessary. The dimensions of the solar field are in the midpoint between those of the Khi Solar One system by Abengoa in the Northern Cape (South Africa) [79] and the Spanish Gemasolar system [80]. This aspect is highlighted because the electricity production of the polygeneration plant analyzed in this work is equal to the electricity production of the Spanish system (24.19 MW e vs. 20 MW e ), underlining the possible competitiveness of this plant with other systems supplied by concentrated solar energy.

Plant Operating Components in the Absence of a Sufficient Irradiance
When the chemical looping cannot operate properly (T receiver-reactor < 900 • C), the plant operates as shown in Figure 4.

Thermal Balance of the Plant
The thermal balance is fundamental for assessing system costs. In the present study, the plant is analyzed at steady-state conditions. The thermal needs of the different units of the plant are evaluated.

Thermal Analysis of the Chemical Looping
Solar energy satisfies the heat duty of the reduction reaction and the streams entering the reactor. On the other hand, the oxidation reaction is exothermic, and a fraction of heat can be used in other plant sections. In Supplementary Table S3, the heat duties of both these reactors are shown. The heating of the inlet flows of the oxidation reaction is executed exploiting a portion of the steam produced in SRG2 and SRG1 ( Supplementary  Figures S2-S4).

Thermal Analysis of the SOFC
In the SOFC scheme, a large amount of heat is required to bring air to 600 • C. This thermal requirement (in HEAT-AIR) can be satisfied as follows: 1.
Air is heated in REG1 from 200 • C to 265 • C through the steam produced in SRG4; 2.
After this pre-heating, the air stream can be further heated to 600 • C in two different ways: (a) If the CL unit operates, air can be heated in REG2; If the CL does not operate, air can be heated through an electric heater. This drastically lowers the plant efficiency.

Thermal Analysis of the Reforming Unit
Different thermal flows are required for the biomethane reforming: 1.
The vapor-liquid separation in V-L-SEP requires heat, and a steam portion produced in SRG5 is exploited (see Supplementary Figure S5); 2.
The stream of CO 2 coming from the storage of the anodic exhausts of the SOFC (CO2-RS-1) is heated in HEAT-ST through a steam portion produced in SRG5 (see Supplementary Figure S6); 3.
Concentrated solar irradiance supplies the reforming reaction and heats the entering streams in the reactor. The thermal requirements of this section are listed in Supplementary Table S4.

Thermal Analysis of the Distillation Unit
Before the entry into the different distillation columns, the liquid is heated through a portion of steam produced in SRG2 (see Supplementary Figures S7 and S8). Additionally, the reboilers of the distillation columns [81] are supplied by the thermal load produced in the post-combustion reactor.

Plant Efficiencies
In this section, different plant efficiencies are evaluated considering the two operating conditions.

Electric Efficiency of the System
The net electric power produced from the plant is: W EL,NET,CL-ON = 6.17 MW e when the CL is in the ON-state and W EL,NET,CL-OFF = 28.96 MW e when the CL is in the OFFstate. The electric efficiency of the system is evaluated considering, respectively, the operation (15) or not (16) The large difference between these two electric efficiencies is mainly related to the following relation W EL,TOT,CL−OFF >> W EL,TOT,CL−ON . When the CL is in the OFF-state, steam streams can be exploited, and the following components are in the OFF-state (COMP-SYR, COMPR-1, and COMPR-2).

Thermal Efficiency of the System
The waste steam and hot water streams produced in the five HRSGs could be mixed (see Supplementary Figures S9 and S10) and cooled to their initial temperature (15 • C) to produce thermal power. When the chemical looping is in the ON-state, W th,tot,CL-ON = 111.97 MW t ; otherwise, W th,tot,CL-OFF = 35.82 MW t . The thermal efficiencies are, respectively, calculated as (17)  This large difference is mainly related to the different fuels required by the system whether the CL is in the ON-state or OFF-state.

Solar/Biomethane-to-Fuel Efficiency of the System
The fuel production only occurs when the CL operates. These results highlight the slight decrease in the global efficiency of the polygeneration plant at low solar irradiances.

Results Discussion and Comparison with Similar Plant
The results obtained from the polygenerative plant are examined and compared with a similar study, by Farooqui et al. [53]. The main differences are: 1.
The unit chosen for electric power production by Farooqui et al. [53] is a gas turbine. The thermal energy generated by the oxyfuel combustion chamber is used to operate continuously the chemical looping cycle; therefore, there is no exploitation of solar energy; 2.
The chemical looping executed by Farooqui et al. [53] is not isothermal, the two reactions take place at different temperatures (1312 • C and 900 • C), and it occurs at a higher operating pressure (2 bar).
A comparison between the average yearly outputs and inputs of the two plants is shown in Table 6. As can be seen in Table 6, the main difference between the present study plant remains in the lower electric power and DME production. This is strictly linked to the solar intermittence, which makes the chemical looping not continuous throughout the year. Consequently, although the size of the two plants is comparable, the chemical looping of this work produces a lower average yearly amount of syngas to be used for electricity and DME production. However, it is important to emphasize the high amount of thermal power produced, methanol and syngas, the lower amount of bio-methane needed, and the higher average global efficiency. In addition, among its outputs, there is no carbon dioxide to be sequestrated, because the plant completely reuses the produced CO 2 and executes the CCU.

Improvement in the Productivity of the Plant
In the attempt to increase the plant productivity, it is thought to decouple the chemical looping operation from the intermittence of solar energy. The reduction reactor of the chemical looping is supplied with the thermal load deriving from the SOFC system. It is necessary to lower the reduction reaction temperature at 800 • C (the SOFC operates at 850 • C) and, consequently, for the complete ceria reduction, it is necessary to increase the amount of biomethane, sent to the reduction reactor (as seen in Figure 6), to 0.95. In this plant, the syngas storage system is eliminated because the SOFC anode can be continuously supplied by the reduction reactor of the CL (continuous working). Thus, the SOFC stacks can produce W EL, SOFC = 243.48 MW e .
To fully eliminate the dependence of the polygeneration plant on the solar energy source, it is necessary to introduce a biomethane oxyfuel combustion unit.
The total biomethane requirement of this new unit is: . m BIO-METH-THERMAL-TOT = 5.57 kg/s. This stream is sent to the oxyfuel combustion, according to reaction (23 [82,83], leads to combustion exhausts mainly constituted of H 2 O and CO 2 . However, to obtain pure O 2 , an Air Separation Unit is necessary. The specific energy expenditure of the Air Separation Unit, to produce an O 2 stream with 98% purity, can be assumed equal to 925 kJ/kg O2 [84,85]. As a result, the total power expenditure of the ASU is 18.02 MW e . The new plant performances on a yearly basis are shown in the first column of Table 7. Comparing the new plant with the solar-aided plant and the plant of Farooqui et al. [53], the new system is the best from the point of view of electricity, DME, methanol, and syngas production. However, there is a much higher amount of biomethane consumption. Additionally, both Farooqui et al.'s plant [53] and the newly studied plant have the non-negligible issue of CO 2 sequestration. The biomethane requirement of the new plant could be reduced by eliminating the reforming unit, even if the CO 2 produced by the plant increases. This occurs because, beyond the CO 2 production from the biomethane oxyfuel combustion, the distillation and the SOFC units contribute to CO 2 production. The syngas lacks and CO 2 production decreases the plant efficiency (even if the produced electric power is the highest). The performance of this plant is shown in the second column of Table 7. From a pure planning point of view, there is no absolute optimum solution for the choice of the plant; this choice depends on the main goal: 1.
If there is a need for a medium-sized plant (W el,average requirement < 50 MW) to be installed in a location with a high availability of solar energy and space, the best choice would be the first polygeneration system analyzed with CS integration and no CO 2 emissions; 2.
If a large-sized plant is required (W el,average requirement > 100 MW), with a high availability of biomethane, the best system is the new plant analyzed with the integration of the reforming unit, proving that adequate sequestration of the CO 2 produced is achieved; 3.
If a large-sized plant is required (W el,average requirement > 100 MW), but no biomethane is available to feed the system, the optimal solution is the new system analyzed without the integration of the reforming unit. However, in the latter case, it is important to have an adequate and much larger storage system for the CO 2 sequestering.
Considering the fuel production, the preferred system is the new system without solar energy integration but with the reforming unit; the other two systems are comparable.

Evaluation of the Chemical Looping Model Performance
This section focuses on the evaluation of the chemical looping performance comparing the results obtained from the CL model of the present study with other experimental [54] and modeling [41,55] literature studies. The present study simulates the CeO 2 /Ce 2 O 3 chemical looping using two Gibbs reactors [41]. Thus, only the thermodynamics of the reactions is considered, and the kinetic limitations of the chemical reactions are neglected. The discrepancy between the thermodynamic simulation and the actual chemical reactions can be seen by comparing the model results of this work with the experimental study of the combined reduction of ceria and methane reforming in a particle transport reactor driven by solar energy [54]. However, these experiments [54] are executed at a temperature equal to 1302 • C (Table 8). Table 8. Comparison between the results of the reduction reactor model of this paper and the experimental data in the literature [54]. The main points of difference between thermodynamic simulation and real experiments are:

Reduction Reaction Set-Up
1.
Absence of carbon deposits in the model. Thermodynamically, the phenomenon is verifiable at methane-to-ceria ratios above 1 and a temperature above 900 • C [41]; 2.
High methane conversion compared to the ideal case.
It follows that the evaluation of the chemical looping performance is made by considering other similar thermodynamic models, such as Bose et al. [41] and Warren et al. [55]. An agreement between the results is obtained, as shown in Supplementary Tables S5-S7. The operating conditions between Bose and the present study are similar, with an operating temperature of 900 • C and a ratio of 0.8 (CH 4 /CeO 2 ). The molar concentration of carbon monoxide is similar, while the hydrogen fraction is slightly higher for Bose. This is due to a more uniform temperature distribution.

Conclusions
This study aims to evaluate the performance of a polygeneration plant to produce heat, electricity, DME, syngas, and methanol with the integration of a chemical looping CeO 2 /Ce 2 O 3 driven by solar energy and coupled to the biomethane reforming system. Some advantages of this system can be listed below: 1.
The plant using only renewable energy sources contributes to the industrial decarbonization; 2.
A carbon capture and utilization unit is implemented. The system is very versatile and environmentally sustainable while still relying on fossil fuel energy production to a small extent. The SOFC unit could be replaced by an existing fossil-fuel power plant, whose effluent can be sent to the CL oxidation reactor and the reforming unit; 3.
Production of innovative, green fuels that can help reduce the dependence on liquefied petroleum gases; 4.
The use of the SOFC stack for electricity and heat generation. The efficiency of these systems is high, especially considering the electric share; 5.
The production of syngas as an energy carrier to store and exploit the solar aleatoriness.
However, the coupling between the chemical looping and a solar energy system produces intermittence in the CL system. The absence of solar exploitation leads to a more complex system with greater production of electricity and green fuels. The following consequences arise:
CO 2 production in the oxyfuel unit. This leads to problems with the CO 2 produced.
Considering the decarbonization, the energy transition, and the reduction of greenhouse emissions as the main focus, the best system to install would be a polygeneration plant with the integration of a solar-powered chemical looping. Further studies could be conducted on this system, such as: a.
Evaluation of a kinetic model for chemical looping operation; b.
An economic analysis to assess the net present value and payback time of the investment; c.
An exergetic analysis to identify the components that can be improved in the complex process; d.
A more detailed assessment of the environmental impact and a discussion related to the material recovery of the various components.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/ma16010315/s1, Figure S1: SOFC scheme in Aspen Plus when the CL is in OFF-state; Figure S2: Thermal balance of HEA-CDS1 and HEA-CDS2 in Aspen Plus; Figure S3: Thermal balance of HEAT-CO2 in Aspen Plus; Figure S4: Thermal balance of HEAT-H2O in Aspen Plus; Figure S5: Thermal balance of HEAT-RE in Aspen Plus; Figure S6: Thermal balance of HEAT-ST in Aspen Plus; Figure S7: Thermal balance HEAT-D2 in Aspen Plus; Figure S8: Thermal balance of HEAT-D3 in Aspen Plus; Figure S9: Mixer of the steam and hot water streams of the plant and simulation of the cooling of the output stream (STE-TOT) to obtain useful thermal power (THE-REQ) when CL is in ON-state; Figure S10: Simulation of the thermal recovery from steam when CL is in OFF-state; Table S1: Particle size distribution of the ceria as a function of the frequency obtained from ELPI measurements [63]; Table S2: Cyclones' characteristics in Aspen Plus [65]; Table S3: Thermal balance of the chemical looping; Table S4: Thermal balance of the reforming unit; Table S5: Comparison of results of the reduction reactor model of this paper with the results reported by Bose et al. [36] at the chosen operating conditions of the solar aided CL;  Funding: This research received no external funding.

Data Availability Statement:
Original data can be requested to the corresponding author.

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