Multiscale Modeling for Reversible Solid Oxide Cell Operation

Solid Oxide Cells (SOCs) can work efficiently in reversible operation, allowing the energy storage as hydrogen in power to gas application and providing requested electricity in gas to power application. They can easily switch from fuel cell to electrolyzer mode in order to guarantee the production of electricity, heat or directly hydrogen as fuel depending on energy demand and utilization. The proposed modeling is able to calculate effectively SOC performance in both operating modes, basing on the same electrochemical equations and system parameters, just setting the current density direction. The identified kinetic core is implemented in different simulation tools as a function of the scale under study. When the analysis mainly focuses on the kinetics affecting the global performance of small-sized single cells, a 0D code written in Fortran and then executed in Aspen Plus is used. When larger-scale single or stacked cells are considered and local maps of the main physicochemical properties on the cell plane are of interest, a detailed in-home 2D Fortran code is carried out. The presented modeling is validated on experimental data collected on laboratory SOCs of different scales and electrode materials, showing a good agreement between calculated and measured values and so confirming its applicability for multiscale approach studies.


Introduction
In view of the continuous increase in energy needs and derived environmental issues, a strict normative has been introduced to favor the development of more sustainable processes [1]. The high emissions of Greenhouse Gas (GHG) mainly from artificial sources such as power plants and transport sectors are one of the leading causes of pollution. Indeed, GHG atmospheric concentration is almost doubled from the preindustrial period to the current days [2]. The global goal is the achievement of thermoneutrality by 2050. This means that a balance between the emissions due to human activities and the needed GHG reduction must be achieved [3]. For this reason, renewable sources have been incentivized to decrease common fossil fuel extraction. However, their high uncertainty and irregular geographical distribution have limited their spread on the energy market. Several storage systems have been studied to provide the energy produced through renewable sources when it is actually required. Some possible storage examples are capacitors, batteries, superconducting magnetic energy systems, flywheels, aquifer energy storage systems, compressed air devices, pumped hydro and chemical storages [4,5]. Among these, only the last ones can guarantee long-term applications and an easy control discharge. Moreover, since the pumped hydro storage needs specific land features, the best solution results in the use of chemical energy carriers of which hydrogen is the most important [6]. H 2 can be used directly as fuel [7,8] or alternatively converted into other fuels, if a proper net distribution is not available. The addition of a downstream methanation unit permits to obtain synthetic natural gas that can be directly fed in existing pipelines [9]. Light hydrocarbons are derived from the combination of electrolyzer and Fischer-Tropsch reactor [10]. In addition, also methanol [11,12] and ammonia [13] are interesting chemical carriers due to their high energy density.
There are several possible ways to synthesize H 2 using thermochemical, biological or electrochemical processes. Natural gas steam reforming and coal gasification are the main routes in industrial applications, despite relevant environmental issues [6,14]. A more sustainable alternative consists of the use of biomass through gasification [15][16][17], pyrolysis [18] or biological treatments [19]. In all these cases, the product is a syngas with variable percentages of H 2 , H 2 O, CO, CO 2 and light hydrocarbons. Differently, water dissociation by means of electrolysis cells allows for the production of emission-free H 2 with a marked purity level [20]. Among different types of electrochemical cells, high-temperature molten carbon and solid oxide cells are the preferential designs as they guarantee optimized efficiency and low cost in both operation modes. Thus, fuel cell technologies might become the core of cogeneration plants producing electricity and heat at the same time [21,22]. On the other hand, a fraction of the requested power for the endothermicity of water dissociation is provided by the high working temperature as thermal energy [23,24]. In particular, reversible Solid Oxide Cells (rSOCs) have good electrical conductivity and fast kinetics, without requiring noble metals as catalysts. Since all layers are solid, their handling is easy and liquid leaks are nonexistent compared to other cell designs, such as Polymer Electrolyte Membrane Cells (PEMCs) [25]. Moreover, these low-temperature cells cannot reach high efficiency since the common catalysts used for oxygen reduction in fuel cell applications have low performance for oxygen evolution in electrolysis operation [26], therefore requiring the use of two distinct units [27]. Contrarily, ceramic materials are effective catalysts for both modes, making SOC a bifunctional unit.
rSOCs have been widely investigated from theoretical, environmental and economical points of view. The basic application consists of standalone systems that provide the requested energy as electricity and heat in isolated microgrids or self-sufficient buildings [28]. rSOCs can be integrated into existing electricity-gas networks where H 2 is also directly used for public transport. In comparison with current fossil fuel-based configurations, the GHG emissions can be decreased by 5-50% [29,30]. Moreover, these cells allow revamping industrial plants in order to optimize the internal energy cycle [31]. In all these fields, a proper simulation is fundamental to study the feasibility of proposed alternative energy production routes. SOC electrochemistry prediction is the main challenging aspect. Different specific kinetics are presented for Solid Oxide Fuel Cell (SOFC) [32,33] and Solid Oxide Electrolyzer Cell (SOEC) [24,34]. Nonetheless, only a minor number of studies have been validated simultaneously for two modes using either a macroscale formulation or detailed microkinetic approach [35][36][37][38].
In this context, this study aims at presenting a multiscale rSOC modeling that can describe single and stacked cell performance for both operation modes. Depending on the application goal, the most suitable simulation tool has been chosen to minimize computational efforts. When a small-sized single cell has been tested, a 0D stationary model, written in Fortran and then integrated into the process simulation software Aspen Plus, has been used. In this case, the analysis has been focused mainly on the electrochemical kinetics identification to evaluate the global system performance. Differently from the common literature approach, both activation and diffusion polarizations have been derived from a rigorous formulation of the Butler-Volmer equation where temperature and gas composition dependencies are underlined. According to this approach, the diffusion overpotential term is influenced by reactant kinetic orders, whereas thermodynamic formulation is commonly followed in literature [39]. A home-made 2D Fortran code has been proposed for stack simulation to guarantee the local control of system behavior on the cell plane and to reduce the evolution of possible stressed working conditions. At both scales, the work has studied an anionic conductive electrolyte planar SOC configuration with different perovskite air electrodes, using as reference a semi-empirical kinetics approach developed and preliminarily validated only for fuel cell operation [40,41]. Then, the model has been tuned to simulate also SOEC behavior, considering that the detected SOFC formulation is still efficient when the same cell materials and designs are used. Firstly, fuel composition has been changed in order to verify the need to introduce the reactant dependency on activation polarization terms in both modes. Subsequently, the stack thermal behavior has been also modeled to evaluate temperature gradients due to occurring reactions. A good match between simulated and experimental rSOC performance in a wide range of working conditions has confirmed the validity of the proposed general approach at both single and stacked cell level and in two operation modes.

rSOC Modeling
A multiscale study is a common approach in chemical engineering, since it allows the use of the most effective simulation tool, in view of the available and requested level of detail [42]. When a small-sized cell is considered, high-level modeling requires relevant simulation time without providing more useful data on the system operation. At this scale, the cell surface has a quite homogeneous behavior, and it can be considered as a perfect stirred reactor. On the contrary, larger-scale units need a more specific investigation to highlight the influence of physicochemical property changes on global cell performance. In this case, the model can also be refined to obtain a better match with system behavior, introducing different kinetics equations and parameters in function of the considered point on the cell plane [43].
In view of the specific scale of analyzed case studies, two simulation tools have been used. A 0D code written in Fortran and implemented in Aspen Plus has allowed evaluating global cell performance by solving material and energy balances at the anodic and cathodic sides. This modeling approach has been proposed to study small-scale single cells focusing on the electrokinetic formulation, thanks to a simulation tool that provides results in a few seconds due to minimum computational efforts. This simplified approach has allowed for a preliminary evaluation of the reversible cell operation, easily switching from fuel cell to electrolysis mode. In addition, the implementation in Aspen Plus has favored a quick and easy execution of sensitivity analysis and process investigation. The model has assumed that the system reaches the stationary working condition and all balances have been solved in the gas domain, considering the exchanges with cell solid elements from a macroscopic point of view.
When larger-scale single or stacked cells are under study, the local performance acquires relevance and the cell dimension can no longer be overlooked. In this case, a 2D simulation code, previously validated for SOFC [41] and here improved to model also the electrolysis operation, has been used. Assuming uniform operating conditions for stacked cells to reduce computational efforts, the simulation has converged in just one minute. Considering an isopotential behavior, anodic and cathodic planes have been divided into an optimized number of subelements, where steady-state material, energy, momentum and charge balances have been formulated (for a more detailed description refer to [44]). In particular, three domains have been considered: fuel electrode gas, air electrode gas and the single solid cell (i.e., anode, cathode, electrolyte and current collectors). The related differential equations have been solved introducing the finite difference approach to simplify the initial problem in a system of ordinary equations. Since each variable has been calculated at specific local conditions, a detailed description of the transport and reaction mechanisms occurring inside the cell has been obtained. Additional information has been provided, which cannot be easily measured at the experimental level such as gas composition, temperature and polarization loss distribution on the cell plane.

Electrochemical Kinetics for Fuel Cell and Electrolyzer Modes
The electrochemical kinetics focusing on the specific description of cell operations is the core of both SOFC and SOEC models. Taking as reference a previous work of the authors [40], a more detailed Energies 2020, 13, 5058 4 of 16 formulation has been developed considering also the possibility of SOEC mode without modifying the basic assumptions validated for SOFC.
An anionic conductive electrolyte cell type is considered. In these cells, O = species are the charge carriers between two electrodes thanks to the migration of oxygen vacancies present in electrolyte lattice that move in the opposite direction of the ions [45]. At each side, semiredox reactions occur as a function of the operation mode: H 2 oxidation or steam reduction at the fuel electrode, O 2 production or consumption at the air electrode (refer to Table 1). Table 1. Occurring electrochemical reactions.

Reaction SOFC SOEC
Anodic The cell effective performance is evaluated starting from the equilibrium potential E eq penalized by several contributions η that cause a reduction of the measured voltage in fuel cell operation and an increase during electrolysis (Equation (1)).
The equilibrium voltage is formulated according to the Nernst equation, considering operating temperature T and reactant composition in terms of partial pressure p (Equation (2)). Since it represents the potential at equilibrium condition, it is valid for both directions of the redox reaction.
where the reversible potential E 0 is a function of temperature according to Equation (3) [39]: However, the effective value of OCV (Open Circuit Voltage) is lower due to the possible gas leakage phenomena or electron migrations inside the electrolyte. To consider this aspect, a term η l is introduced as a further polarization loss (Equation (4)) [46,47].
The overpotential terms consider different resistances that penalize the electrochemical phenomena. A relevant penalization is due to the electron conductivity and the migration of O = ions from the fuel to the air electrode and vice versa, which are assumed as activated thermal processes following the Ohm law [44]. Consequently, the ohmic overpotential η ohm is expressed as Equation (5), showing a linear dependence on the applied current density J.
where P 1 and P 2 are two empirical parameters specific for every case study.
The overpotential due to electrochemical reactions η el is derived from the Butler-Volmer formulation that represents the rate-limiting step in the function of the applied current density [48]. It depends on the direct (i.e., oxidation) and indirect (i.e., reduction) reactions occurring at each electrode in the function of bulk partial pressures p bulk and Three Phase Boundary (TPB) partial pressures p TPB of the reactants (Equation (6)). Thus, the same kinetic approach has to be followed in both operation modes.
When a low electric load is applied, the actual redox exchange is the slowest step and so reactant pressure gradients are a minimum value, i.e., the ratio between bulk and TPB values (p TPB /p bulk ) equal to one (Equation (7)).
Equation (7) can be rewritten according to a hyperbolic sine formulation for the activation polarization η act , depending on temperature and reactant compositions [49]. Only bulk values are considered since diffusion issues along electrode thickness are negligible at low current density [40]. Equations (8) and (9) are formulated assuming atmospheric working pressure and the charge transfer coefficient α equal to 0.5. It is worth to underline that the exchange current density J 0,el has to show the same dependencies on the working conditions (i.e., temperature and reactant composition) in both modes since it represents the forward and reverse electrode reaction rate at the equilibrium state [24,34].
where the pre-exponential coefficients γ, the activation energies E act and the reactant kinetics orders A, B and C are specific for every occurring reaction path. Differently, when a high current density is applied, the gas diffusion inside the microstructure of the electrode becomes the rate-limiting step. This effect is represented by the diffusion overpotential η diff , which considers the difference between bulk and TPB partial pressures as driving force [50]. It is again derived from Butler-Volmer formulation (Equation (6)), but in this case, the electrochemical reaction is faster than transport mechanisms. Thus, the exchange current density J 0,el representing the kinetic rate at equilibrium state tends to infinity at this condition (Equation (10)).
The concentration gradient at the oxygen electrode side is relevant only at very low O 2 partial pressure (p O2 < 0.05 atm), so it is usually neglected in common working conditions with air feed [51]. Differently, the contribution at the fuel electrode for both operations is represented by using an equal formulation through the introduction of the absolute value in order to guarantee the right ratio among reactant compositions, independently from the current density direction (Equation (11)). By solving Equation (10) for the fuel electrode side and fixing α equal to 0.5, the diffusion overpotential is derived from Butler-Volmer equation in contrast with the common literature formulation. As Equation (11) Energies 2020, 13, 5058 6 of 16 shows, this results in the direct influence of the kinetic orders, previously evaluated in the exchange current density expression (Equation (8)), on gas partial pressures.
To evaluate the reactant pressures at TBP active sites, material balances are solved along the fuel electrode thickness δ fuel assuming a Fickian mechanism as the main transport phenomenon. Regarding the diffusion coefficients D, a combination of molecular and Knudsen terms is taken into account. H 2 and H 2 O average values in the detected gas profile are considered as TPB partial pressures, as Equations (12) and (13) show for SOFC and SOEC modes, respectively (for a detailed description refer to [34,41] Substituting Equations (12) and (13) into Equation (11), the diffusion overpotential at fuel electrode η diff,fuel is estimated (Equation (14)). In this approach, the reaction is assumed to occur along electrode length and not only at the electrode-electrolyte interface [33], thanks to the mixed electron and ion conductivity of used ceramic materials.
The proposed approach has been followed in both 0D and 2D simulation. In the first case, an average between input and output gas partial pressures has been considered as bulk composition. In the second case, specific values have been used in every defined subelement on the cell plane thanks to the local modeling.
In light of this final observation, a further term is introduced. The Nernst overpotential η Nernst considers the reference Open Circuit Voltage owing to the inhomogeneous concentration distribution that occurs on the cell plane under electric load. In the 2D approach, this term has not been considered, since it is possible to evaluate directly the cell behavior knowing actual local features in each subdomain [44]. Thus, it has been only added to the 0D modeling code. For both operations, Equation (15) defines this further overpotential as the difference of Open Circuit Voltage calculated in initial p in and average p av (i.e., between input and output values) conditions [52].
where the first term is referred to the air electrode contribution and the second to the fuel electrode one. All these formulations are valid for both SOFC and SOEC operations, so that the model can easily switch between two operating modes by only changing the current direction without the need to modify kinetic parameters [37].

Zero-Dimensional (0D) Simulation of rSOC Small-Scale Single Cell
The 0D approach has been used to evaluate the global performance of a rSOC small-sized single cell. The analysis has been focused on electrochemical kinetics to verify if the proposed approach, previously validated only for fuel cell mode [40], can represent both operations.
The study has considered literature tests of a planar anionic conductive electrolyte cell with an active area equal to 8 cm 2 . The cell is structured by a porous Ni/YSZ (Ni/Yttria-stabilized Zirconia) layer fuel electrode that is the support for the YSZ electrolyte and LSM (Lanthanum Strontium Manganite) air electrode [53]. Dry air is fed to the air electrode, whereas an H 2 /H 2 O mixture is used at the fuel electrode side working at constant flow rates (0.14 and 0.025 m 3 h −1 , respectively).
The set of parameters necessary for the application of the above described semi-empirical kinetics model has been identified through the fitting of experimental polarization curves ( Table 2). The different kinetics parameters derived from the fitting of experimental results are comparable with previous reference studies. Activation energies are usually assumed equal to 100 and 120 kJ mol −1 for fuel and air electrode, respectively [38]. Pre-exponential coefficients are set between 10 4 and 10 6 A cm −2 [51], whereas reaction orders vary from −0.5 to 2 in SOFC operation [54]. These values should be the same for SOEC mode, yet the common literature approach [11,34,38] neglects reactant dependencies in activation overpotentials, differently from the proposed formulation (Equations (8) and (9)). All few studies follow this method: for instance, the fuel electrode reaction orders are equal to 0.11 and 0.67 for A and B, respectively, while γ fuel results in 2.5 × 10 5 A cm −2 in [36]. Such values are quite similar to the identified ones.
The model has been validated for both operating modes, testing cell behavior at 50/50 H 2 /H 2 O mixture as the fuel electrode feed and two different temperatures (i.e., 1023 and 1123 K). As Figure 1 shows, a good match between experimental and simulated results has been obtained confirming the possibility to use the same set of parameters. Only the current direction has to be updated switching from electrolysis to fuel cell mode. Detected trends confirm the relevant influence of the temperature on cell performance, due to a thermally occurring activate transport and reaction phenomena. As expected, the temperature increase favors the cell voltage reducing overpotential contributions.
Energies 2020, 13, 5058 8 of 16 mixture as the fuel electrode feed and two different temperatures (i.e., 1023 and 1123 K). As Figure 1 shows, a good match between experimental and simulated results has been obtained confirming the possibility to use the same set of parameters. Only the current direction has to be updated switching from electrolysis to fuel cell mode. Detected trends confirm the relevant influence of the temperature on cell performance, due to a thermally occurring activate transport and reaction phenomena. As expected, the temperature increase favors the cell voltage reducing overpotential contributions.

Two-Dimensional (2D) Simulation of rSOC Larger-Scale Stacked Cells
For 2D modeling validation, a rSOC short stack case study has been considered, consisting of six planar anode supported cells of 80 cm 2 active area supplied by SOLIDPower. The materials used are Ni/8YSZ for the fuel electrode, 8YSZ for the electrolyte and GDC + LSCF (Gadolinium-doped Ceria + Lanthanum Strontium Cobalt Ferrite) for the bilateral layers. Thus, in this case, a different type of perovskite is present, also adding a protective layer between the electrolyte and the air electrode to avoid possible poisoning effects. Tests have been performed at equal initial temperature (i.e., 1023 K), feeding dry air and varying H 2 /H 2 O composition as the fuel mixture. Two sets of data have been analyzed. Initially, system performance has been estimated through ten polarization curves at fixed flow rates and furnace temperature in order to underline fed reactant dependencies in both SOFC and SOEC modes. Later different working points have been tested setting fuel reactant utilization and reducing furnace temperature to underline the actual cell thermal behavior. The complete description of the test rig and the experimental campaign is explained in [55].
As for the 0D approach, a careful model tuning to evaluate cell kinetics through the fitting of experimental polarization curves is necessary. In this first step, the analysis has focused on highlighting the stack performance dependence on fed hydrogen in order to verify if a more general approach, where reactant influences are considered also for the activation overpotential (Equations (8) and (9)) differently from previous literature, could provide a more detailed simulation or only an increase of computational efforts. Focusing on this point, an isothermal model has been assumed referring to experimental tests where the furnace temperature has been set in order to balance temperature changes due to redox reactions. Table 3 shows the identified values.
Comparing the obtained set with the previous one related to the single cell, it is worth to underline that only the parameters concerning specific cell features (i.e., materials and design) have been changed. Kinetic order and activation energy values have been maintained equally as they refer to the same anodic and cathodic reaction paths despite the presence of a different air electrode. The resulting model has been again validated for both operation modes. As Figures 2 and 3 show, a good agreement has been obtained comparing experimental and simulated results with an error lower than 4%. The most relevant mismatch has been measured at high current density in SOEC operation where a relevant electric noise has been detected during lab tests. The polarization curves underline that the increase of reactant compositions (i.e., H 2 for fuel cell and H 2 O for electrolysis) improves stack performance. Under a current of 0.3 A cm −2 , a difference of 0.3 V has been measured by comparing tests in rich and poor H 2 percentages. This discrepancy reaches almost 0.4 V in SOEC operation at −0.3 A cm −2 . According to the proposed simulation, the polarization curve improvement is mainly due to activation Energies 2020, 13, 5058 9 of 16 overpotential reduction (Equation (8)), confirming the need for the proposed formulation for effective cell modeling. On the contrary, diffusion overpotentials have a lower influence as the cell works far from the limit current density.     After successful validation of the electrochemical kinetics at different feed compositions, the study has focused on the system thermal behavior to highlight the temperature increase or decrease After successful validation of the electrochemical kinetics at different feed compositions, the study has focused on the system thermal behavior to highlight the temperature increase or decrease due to redox reactions. In this case, also energy balances have been solved without varying the previously evaluated kinetic parameters except for the ohmic coefficients (Table 4). These are the parameters that show the main temperature dependence, and consequently, a minimum adjustment was required in changing from an isothermal to a temperature-dependent simulation. In this case, also the system parameters referred to heat exchange phenomena have to be taken into account. The main used parameters are summarized in Table 5, showing similar values to literature ones [56]. Different working conditions have been considered varying gas compositions and flow rates to guarantee a fixed Reactant Utilization (RU) for both operating modes (i.e., hydrogen RU equal to 0.7 and oxygen RU equal to 0.25 in SOFC, water RU equal to 0.7 in SOEC). Figure 4 shows the measured data expressing the rSOC performance as an average among different voltage values of stacked cells. The good comparison confirms again the validity of model assumptions.
Heat capacity [cal mol −1 K −1 ] 50 34 29 Global cell conductivity [W m −1 K −1 ] 2.0 Different working conditions have been considered varying gas compositions and flow rates to guarantee a fixed Reactant Utilization (RU) for both operating modes (i.e., hydrogen RU equal to 0.7 and oxygen RU equal to 0.25 in SOFC, water RU equal to 0.7 in SOEC). Figure 4 shows the measured data expressing the rSOC performance as an average among different voltage values of stacked cells. The good comparison confirms again the validity of model assumptions.  To also validate the thermal behavior of the simulated stack, the temperature evolution has been studied starting from the open circuit condition at an initial system temperature equal to 1021 K. It has been focused on the air electrode outlet temperature, since it is a representative value of the global system temperature due to its high flow rate and heat capacity [55]. Quite similar profiles have been detected for all the tested fuel compositions. The model validation of both SOFC and SOEC modes puts in evidence a good match with a maximum discrepancy of 3 K compared to experimental measurements. In fuel cell operation, the process exothermicity causes a temperature increase from 1020 to 1033 K, whereas a variety of low degrees characterizes tests in electrolysis.
In addition to lab test results, 2D simulation favors a detailed analysis of cell temperature, assuming an equal value for all layers of each single unit (i.e., electrodes, electrolyte and current collectors). Through this local approach, it has been possible to know the actual temperature on the system plane, detecting the maximum and average values reached. For instance, Table 6 shows the obtained results under common working loads at different fuel compositions. A quite uniform profile has been underlined comparing cell temperature gradients (first two columns), so safety working conditions are guaranteed without issues due to creep or delamination that represent relevant causes of degradation in high working temperature cells. Moreover, these data are also useful to validate the following experimental approach. A difference of only some degrees has been detected between the maximum simulated solid temperature and the tested air electrode outlet. Therefore, this last one can effectively represent the stack temperature with a good approximation permitting a simpler measurement during rSOC operation.

Conclusions
A multiscale approach has been proposed to simulate rSOC performance, starting from a 0D model executed in Aspen Plus for a small-scale single cell to a 2D Fortran code for stacked cells. Indeed, depending on system size, the most suitable tool has to be used guaranteeing minimum computational efforts without loss of information. The calculation of the inhomogeneous distribution of temperature and gas composition on the cell plane due to occurring local phenomena cannot be neglected for large-sized cells since they have a direct influence on global cell performance, so that this results in a more precise simulation tool. Both lumped and local approaches have been based on conservation equations to describe material and energy transport phenomena, while a semi-empirical formulation has been followed to represent the system electrochemistry. Firstly, common electrochemical kinetics has been derived from Ohm and Butler-Volmer equations for SOFC and SOEC operations and successfully validated through 0D model for a single cell. To simulate the stack configuration, a 2D local analysis has been introduced. The previous detected activation energies (E act,fuel = 100 kJ mol −1 and E act,air = 120 kJ mol −1 ) and kinetic orders (A = 0.5, B = 0.5 and C = 0.25) have been maintained unchanged, since they are linked to intrinsic redox reactions properties and they are validated for both air electrode materials. On the contrary, the parameters depending on the specific used design have been newly tuned. In this second configuration, the model has been tested at different fuel electrode feeds and a fixed operating temperature highlighting the importance to consider the same electrochemical kinetics dependencies for both operation modes (resulting in an error between experimental and simulated data lower than 4%). Later, the system thermal behavior has been evaluated: the model follows adequately the obtained experimental profiles at different operating conditions, with only a few mismatched degrees. SOFC operation results are more sensitive to temperature variation due to an increase of 10 K, whereas the thermal trend is more homogeneous in SOEC mode. Available measured temperatures of the air electrode outlet have been assumed as a representative value of the cell working one, after a preliminary assessment of effective occurring thermal gradients on the cell plane. The detected differences among maximum, average and outlet temperatures are lower than 2 K, confirming the feasibility of this simplified assumption.