Optimization Design of Rib Width and Performance Analysis of Solid Oxide Electrolysis Cell

: Structure design is of great value for the performance improvement of solid oxide electrolysis cells (SOECs) to diminish the gap between scientiﬁc research and industrial application. A comprehensive multi-physics coupled model is constructed to conduct parameter sensitivity analysis to reveal the primary and secondary factors on the SOEC performance and optimal rib width. It is found that the parameters of the O 2 electrode have almost no inﬂuence on the optimal rib width at the H 2 electrode side and vice versa. The optimized rib width is not sensitive to the electrode porosity, thickness, electrical conductivity and gas composition. The optimal rib width at the H 2 electrode side is sensitive to the contact resistance at the interface between the electrode and interconnect rib, while the extremely small concentration loss at the O 2 electrode leads to the insensitivity of optimal rib width to the parameters inﬂuencing the O 2 di ﬀ usion. In addition to the contact resistance, the applied cell voltage and pitch width also has a dramatic inﬂuence on the optimal rib width of the fuel electrode. An analytical expression considering the inﬂuence of total cell polarization loss, the pitch width and the contact resistance is further developed for the beneﬁt of the engineering society. The maximum error in the cell performance between the numerically obtained and analytically acquired optimal rib width is only 0.14% and the predictive power of the analytical formula is fully veriﬁed.


Introduction
Renewable energy resources, including solar, wind, tidal, and biomass, are of great significance as the fossil energy crisis is becoming increasingly serious. However, their intermittence leads to an undesirable imbalance between demand and supply [1]. An energy storage device is required so that energy can be stored and released as needed. Among various candidates for energy storage, solid oxide electrolysis cells (SOECs) have gained popularity for higher efficiency and lower pollution and even no pollution. SOEC is an energy conversion device that can convert electrical energy and heat to chemical energy, by splitting H 2 O/CO 2 to produce H 2 and CO. The products H 2 and CO can be used as fuel in a solid oxide fuel cell (SOFC) to produce electricity or be stored as raw materials for the synthesis of hydrocarbons via the Fischer-Tropsch reaction. SOEC technology is of great superiority and prospect. Among three main electrolysis configurations, it has been reported that the efficiency of hydrogen production by high temperature SOEC is more than twice of that by an Alkaline electrolysis cell, and is 1.5 times of that by proton exchange membrane electrolyzer [2]. There are mainly three kinds of SOEC according to reaction gas species: high temperature H 2 O electrolysis, CO 2 electrolysis, H 2 O and CO 2 co-electrolysis cell. Steam electrolysis can produce H 2 , which is a completely environmentally friendly fuel. In addition, H 2 O is rich in nature. CO 2 electrolysis is advantageous in that it can consume CO 2 and relieve the greenhouse effect, and the product CO is easier to store and transport than H 2 . However, it has a potential carbon deposit risk. H 2 O and CO 2 co-electrolysis can produce H 2 and CO mixtures, and by adjusting the inlet H 2 O/CO 2 ratio, it can produce applicable hydrocarbon synthesis. The steam electrolysis has the highest electrolysis efficiency while CO 2 electrolysis has the lowest efficiency, and the efficiency of H 2 O/CO 2 co-electrolysis is between them. Due to the environmental friendliness and higher efficiency of steam electrolysis, high temperature H 2 O electrolysis attracts increased attention. With the increase of temperature, the electricity needed to electrolyze H 2 O decreases, while the low quality heat needed increases. Moreover, high temperature SOEC is not only thermodynamically beneficial but also kinetically favorable. Hence, high temperature SOEC has a more promising application perspective.
Materials, performance and degradation issues are still three challenges of SOEC technology to be settled. Extensive research about SOEC concentrates on the optimization of material microstructure [3][4][5], geometrical [6] and operating parameters [7] and the analysis to improve SOEC performance. Notice, however, the structural size choice is also of great significance. For example, the authors of [8] studied the effect of cathode thickness on CO 2 /H 2 O co-electrolysis performance under various operating conditions by experiment, which reveals that SOEC performance can be substantially improved by decreasing cathode thickness. For the time being, most of the SOEC researches were conducted experimentally [2,[9][10][11][12]. Unfortunately, the experiment is expensive and time-consuming, so research about structure optimization by experiment is rare. Simulation is an efficient alternative to help the design of SOEC/SOFC to improve performance, especially when exploring a large combination of operating and structural parameters space. For example, Reference [1,[11][12][13] studied the influence of pressure on SOEC. Ni et al. [14] researched electrode thickness, support type, electrode porosity and pore size and operating pressure on SOEC performance. It is concluded that anode-supported SOEC has the best output performance. Kong et al. [15] examined the impacts of different electrode-rib contact resistances, fuel compositions, electrode porosity, electrode thicknesses and electrode conductivity on the optimal anode and cathode ribs of SOFC independently.
The development of SOEC is later than that of SOFC. As the reverse process of SOFC, SOEC has basically the same materials system as SOFC. The research of SOFC is enlightening to the development of SOEC. In the last 20 years, during the simulation of SOFC, it found out that rib width design is of great significance to improve the cell performance [15][16][17][18], which has also been experimentally confirmed [19]. Actually, as early as 2003, Lin [16] had provided a phenomenological model and analytical expressions to estimate the rib effects on the concentration and ohmic polarization of anode-supported SOFC stacks. Jeon et al. [20] described a microstructure model and examined the influence of the rib and pitch widths and the electrode-interconnect contact area specific resistance (ASR) on the stack-cell performance. The authors of [17] investigated the effect of ASR between the electrode and the rib on the performance of SOFC, and conducted the rib width optimization by 2D SOFC multi-physical modes. In [18], the authors primarily compared the optimal rib width result attained by 2D and 3D multi-physical models, and revealed that the optimization result of 2D and 3D models are in good agreement. The authors of [15,21] optimized the anode and cathode rib width for anode-supported and cathode-supported SOFCs, respectively. An analytical expression of optimal rib width is deduced to help the engineering design of SOFC. In [22], the authors conducted the optimization of the cylindrical interconnect rib width of SOFC, and it found out that anode and cathode rib width should be optimized separately, which is in consistent with the conclusion in Reference [21].
Reference [23] studied the influence of rib size on the performance of a reversible solid oxide cell. It takes into account the efficiency of a solid oxide fuel cell and electrolysis cell comprehensively.
However, the polarization process of SOEC is vastly different from SOFC. During past years, to our knowledge, it is still not clear how to choose the rib width when fabricating the SOEC interconnector characterized by the rib-channel structure, i.e., the optimization of rib width has not been clearly addressed. This paper conducts a comprehensive parameter sensitivity analysis on the rib width optimization of SOEC by 2D multi-physics simulations. Analytical expressions for the optimal cathode rib width design are obtained to provide an easy-to-use guide for designing the rib-channel layout of SOEC.

Model Description
A typical unit of SOECs contains a dense electrolyte sandwiched by a porous anode and cathode, channels and inter-connector. To obtain better electrochemical performance, a relatively denser and thinner porous layer are generally added between the electrode and electrolyte, which is called the function layer. The working principle of SOEC is illustrated in Figure 1a. Figure 1b shows a 2D cross section of the SOEC repeating unit, displaying the rib-channel design parameters: the pitch width, d pitch , and the rib width, d rib . The multi-physics coupled modeling is applied to a pitch unit of SOEC and considers the mass, species, momentum transfer processes, the current conduction and the electrochemical reaction.

Species Transfer Process
The species conservation equation can be expressed as: Here, N i is the total molar flux vector, which includes the convective flow and the diffusion flow. D i eff , C i , and S i are effective diffusion coefficient, molar concentration and source or sink of species i, respectively. S i is source or sink of species i by electrochemical consumption and production in the electrode. The species diffusion is modeled by a multi-component Dusty-Gas model [24], which is proved to be the most accurate model to simulate gas diffusion in a porous electrode of SOFC [25]. A binary molecular diffusion coefficient is generally used directly in channels, where gas convective transfer rather than the diffusion transfer is dominant, and the binary molecular mutual diffusion coefficient is expressed as [18]: where T, v i and M i is the temperature, diffusion volume and mole mass of species i. P is the total gas pressure.
In the electrode, the influence of porosity and tortuosity should be considered to describe the effective binary molecular diffusion coefficient: ε, τ are porosity and tortuosity of the electrode, respectively. In the electrode, molecular mutual diffusion dominates when the electrode pore is larger than the molecular mean free path, while Knudsen diffusion becomes primary when the pore size is small and the collision between the species and pore wall increases. The effective Knudsen diffusion coefficient is expressed as: Here r g is pore radius, and r g is expressed as: Considering both molecular diffusion and Knudsen diffusion simultaneously, the effective diffusion coefficient in the electrode is formally expressed as [18]:

Electrochemical Reaction Model
The electrochemical reactions include the H 2 O reduction reaction at the fuel electrode and the oxidation reaction at the O 2 electrode, which are expressed as: The electronic current and ion current are governed by charge continuity equations: where n is the number of electrons transferred in the electrochemical reaction, i el and i io are the electronic and ionic current density vector, φ el , φ io are the electronic and ionic potential. λ TPB,eff is the effective three phase boundary (TPB) density per unit. j TPB is the current density at TPB. Nernst potential is the minimum cell voltage needed to drive an electrochemical reaction in SOEC: Here, ∆G is the Gibbs free energy change of electrochemical reaction. P H2 , P O2 , P H2O are partial pressures of H 2 , O 2 and H 2 O. P 0 is the standard atmospheric pressure.
The applied voltage in an operating SOEC is: where η ohmic , η act , η con are ohmic, activation and concentration loss, respectively. Ohmic loss is induced by the conductivity resistance of electron and ion transfer in solid components and is calculated according to Ohm's law: Here, ASR ohm is total area specific resistance. The occurrence of the electrochemical reaction needs to overcome the reaction activation energy barrier, leading to irreversible activation loss. The relationship between current density and the activation loss is described by the Butler-Volmer equations: The concentration polarization is induced by a change of species concentration when SOEC is in operation: E O2 and E H2 are, respectively, the activation energies for the O 2 and H 2 electrode electrochemical reactions.

Effective Material Property Model
The intrinsic material conductivity is temperature dependent [1,6,15,26]: where R is the universal gas constant.
For the porous electrode, the relationship between the material macro property and microstructure can be expressed by coordination number theory and percolation theory [27]: where σ k 0 is the intrinsic electric conductivity of the k th phase material in the dense solid. φ k is the volume fraction of the k th phase particles in the composite material, φ k t is the percolation threshold volume fraction of the k phase particles, which is determined by where Z is the average coordination number for each particle and set as six for a random packing of spheres [15]. r el and r io are, respectively, the electronic conductivity particle radius and the ionic conductivity particle radius. LSCF is a material that can conduct electrons and ions simultaneously. For the LSCF-GDC composite electrode, the effective electronic conductivity is expressed as: The effective ionic conductivity can be considered as the parallel current conduction of LSCF and GDC. The ionic conductivity of the composite electrode is expressed as: where γ is a Bragg factor and is usually set as 1.5. P k is the percolation probability of phase k. Both the percolation probabilities of LSCF and GDC are assumed to be 1. The TPB density per unit volume of a composite electrode with a binary mixture is expressed as [27]: The relevant physical quantities are expresses as follows, More details about the specific symbols are referred to in [27].
The multi-physics simulations are conducted with the commercial finite element software COMSOL 4.3b. Simulations were first carried out to reproduce the performance of an experimental cell that is designated as the standard cell. The parameters used to reproduce the experiment are listed in Table 1, and are referred to as the standard parameter set. Unless explicitly stated otherwise, the parametric analysis was conducted by varying one parameter as the variable while all the other parameters are kept as the standard case.
As can be seen from Figure 2, the experimental and simulation results agree well, demonstrating the validity of the multi-physics model. Naturally, some discrepancies between the experimental and theoretical I-V curves are observed. Nevertheless, the theoretical and experimental difference is less than 0.02 V for any given current densities.

Parameters Sensitivity Analysis for Rib Width Optimization at the H 2 Electrode Side
As seen in Figure 3, the current density first increases to a peak value then decreases with the increase of rib width. This is quite understandable as the concentration loss increases and the ohmic loss decreases when the rib width increases. Therefore, there is an optimized rib width d rib-m where the current density reaches its maximum value. When d rib < d rib-m , with the increase of rib width, the decrease of ohmic loss is higher than the increase of concentration loss, the current density tends to increase. When d rib > d rib-m , with the increase of rib width, the decrease of ohmic loss is smaller than the increase of concentration loss, and the current density decreases. The competition between the concentration polarization and the ohmic polarization determines the final optimized rib width.   As shown in Figure 3a, when the rib/pitch width ratio at the O 2 electrode side increases from 0.25 to 0.5 then to 0.75, the optimized H 2 electrode rib width is 0.752, 0.720, 0.710 mm, respectively. Notice, however, although the optimized rib widths are different, the width of 0.72 mm is sufficiently optimal for all the O 2 electrode rib/pitch ratio conditions. The maximum current density difference between optimal rib width and 0.72 mm are only 0.05%, 0.007% for the rib/pitch ratio of 0.25 and 0.75, respectively. Therefore, it is concluded that the optimized rib width at the H 2 electrode side is not sensitive to the rib/pitch ratio at the O 2 electrode side.
The experimentally measured ASR of SOFC is in the range 0.01-0.05 Ωcm 2 [33]. SOEC has the same electrode and inter-connector materials. However, the oxidation atmosphere of both the electrodes of SOEC means potentially larger oxidation risk of inter-connector than SOFC, so the maximum ASR studied here is 0.08 Ωcm 2 . ASR of both the electrodes in the range of 0.04-0.08 Ωcm 2 is used to conduct parameter sensitivity analysis to optimize rib width. To study the parameters' influence on the rib width optimization at one electrode side, ASR of another electrode side is set as 0.03 Ωcm 2 . As displayed in Figure 3b, when ASR between the O 2 electrode and rib increases from 0.01 to 0.05 Ωcm 2 , the current density decreases gradually because of the increased ohmic loss. The optimized H 2 electrode rib width is 0.689, 0.705, 0.720, 0.734, 0.746 mm when ASR at the O 2 electrode side are 0.01, 0.02, 0.03, 0.04, 0.05 Ωcm 2 , respectively. However, when choosing 0.72 mm as the optimal rib width for the other four cases, the differences between the maximum current density at the optimal rib width and the current density with a rib width of 0.72 mm are 0.0756%, 0.0165%, 0.0128%, 0.04% for ASR 0.01, 0.02, 0.04, 0.05 Ωcm 2 , respectively. Hence, the optimized rib width on the H 2 electrode side is not sensitive to the ASR on the O 2 electrode side.
Along the fuel flow direction, H 2 O is consumed gradually, the molar fraction of H 2 O changes correspondingly. Furthermore, for practically operating SOEC, high H 2 O conversion rate is favorable, so the research about the rib width optimization should take into account different H 2 O molar fractions. Figure 3c shows, for a fixed rib width, the current density increases with the increase of H 2 O molar fraction. This is because the Nernst potential decreases and the electrochemical reaction rate increases with the increased H 2 O molar fraction, and the current density increases when the applied voltage is fixed. The optimal rib width is 0.720, 0.687, 0.648, 0.602 mm when H 2 O molar fractions are 90%, 80%, 70%, 60%, respectively. Because with the increase of H 2 O molar fraction, current density increases, concentration polarization and ohmic polarization both increase simultaneously, which leads to a slight change in the optimized rib width. This result is also in accordance with that of SOFC [15]. Nevertheless, 0.65 mm width is sufficiently optimal for SOEC with different H 2 O molar fraction, because the difference between the maximum current density and the current density at the rib width of 0.65 mm is only 0.0002%, 0.0004%, 0.0001%, 0.0028% for SOEC with the inlet H 2 O molar fraction of 90%, 80%, 70%, 60%, respectively. The optimized rib width can be considered independent of the H 2 O fraction.
The SOEC H 2 electrode porosity is in the range of 0.3-0.5 [34]. As shown in Figure 3d, the optimal rib width increases with the increase of electrode porosity due to the reduced gas diffusion resistance. For the porosity of 0.3, 0.4, 0.5, the optimized rib width is 0.616, 0.720, 0.800 mm, respectively. Even the optimized rib width for the three porosities seems to be quite different, the maximum current densities for porosity 0.3 and 0.5 are only 0.98% and 0.35% higher than the current density with a rib width of 0.720 mm. Therefore, the optimized rib width is not sensitive to the H 2 electrode porosity.
As shown in Figure 3e, when the H 2 electrode electrical conductivity increases from 717.28 to 7172.8 S/m, the current density increases relatively drastically. However, when H 2 electrode electrical conductivity increases from 7172.8 to 71,728 S/m, the current density has almost no change. It means that the electrical conductivity of 7172.8 S/m is large enough, electrical conductivity larger than 7172.8 S/m has no evident improvement for SOEC performance. The optima rib width for the electrode electric conductivity of 717.28, 7172.8 and 71,728 S/m is 0.741, 0.722 and 0.719 mm, respectively. However, using a rib width of 0.720 mm, the current density differs from the maximum only by 0.03% for the conductivity of 717.28 S/m, and virtually zero for the conductivity of 71,728 and 7172.8 S/m. Hence, the optimal rib width is basically independent of the electrode conductivity.
As seen in Figure 3f, the performance of SOEC decreases with the increase of H 2 electrode thickness. The increase of electrode thickness means longer gas diffuse path from the channel-electrode interface to the active three phase boundary. The diffusion process becomes more difficult, and the concentration loss increases accordingly. Moreover, a thicker electrode inevitably leads to larger ohmic loss for current in the thickness direction. Hence, the overall SOEC performance decreases when electrode thickness increases. However, as the electrode thickness increases, the optimized rib width changes only slightly. The optimized rib width for the electrode thickness of 0.2, 0.4 and 0.6 are 0.719, 0.720 and 0.646 mm, respectively. However, the width of 0.720 mm is sufficiently optimal. The current density for the rib width of 0.72 mm differs only 0.006% and 0.41% from the maximum for the electrode thickness of 0.2 and 0.6 mm, respectively. Therefore, the optimal rib width is insensitive to the electrode thickness. The reason behind this is that both the concentration and ohmic polarizations increase comparatively at the same time when the electrode thickness increases.
As displayed in Figure 4a, fixing the rib width, the current density decreases significantly with the increase of rib-electrode ASR due to directly the increased ohmic polarization loss. Moreover, the optimal H 2 electrode rib width increases, because the increased ohmic loss due to the increase of ASR can be partially offset by increasing the rib width. As seen in Figure 4b, the optimized rib width increases almost linearly with the increase of ASR. When ASR increases from 0.01 to 0.05 Ωcm 2 , the optimized rib width increases from 0.511 to 0.840 mm, or an increase of 64.4%. This ASR effect is quite dramatic. With the increase of ASR, the current density decreases, the concentration loss decreases passively. The increased ohmic loss and the decreased concentration loss collectively leads to the dramatic increase of the optimal rib width.

Parameters Sensitivity Analysis for the Rib Width Optimization on the O 2 Electrode Side
As can be seen in Figure 5a, with the increase of rib/pitch ratio on the H 2 electrode side, the current density decreases, which can be attributed to the increased concentration polarization in the cathode due to the increased rib width. When the rib/pitch ratio increases from 0.25 to 0.5, the SOEC performance decreases slightly. However, the when rib/pitch ratio increases from 0.5 to 0.75, the SOEC performance decreases drastically. It can be inferred that when the rib/pitch ratio exceeds 0.5, the negative effect of the concentration loss increase is far greater than the positive effect of the ohmic loss decrease as induced by the increase of the H 2 electrode rib width. However, when the rib width on the O 2 electrode side increases, the current density increases all the way, independent of the rib/pitch ratio on the H 2 electrode side, indicating the optimized rib width on the O 2 electrode side is not sensitive to the rib/pitch ratio on the H 2 electrode side.  As shown in Figure 5b, when ASR on the H 2 electrode side increases from 0.01 Ωcm 2 to 0.05 Ωcm 2 , the current density decreases gradually. For a fixed ASR, the current density increases with the increased rib width on the O 2 electrode side, independent of the ASR change on the H 2 electrode side. Hence, the optimized O 2 -electrode rib width is not sensitive to the ASR of the H 2 electrode. Figure 5c shows the current density increases extremely slightly with the increase of O 2 electrode porosity from 0.3 to 0.4 then to 0.5, indicating the polarization loss induced by O 2 diffusion inside the anode is rather small. This is also ascribed to the thinness of the O 2 electrode. Furthermore, the current density increases with the increase of rib width, for all the anode porosities examined. It can be concluded that the optimal O 2 -electrode rib width is independent of the O 2 electrode porosity. Figure 5d shows the current density increases with the increase of the O 2 electrode electrical conductivity, due to the decreased ohmic polarization. The current density increases more dramatically when the conductivity increases from 280.6 to 2806 S/m, but increases slightly when the conductivity increases from 2806 to 28,060 S/m, indicating the ohmic polarization is no longer a major factor limiting the cell performance when the anode conductivity is above 2806 S/m. For a fixed conductivity, the current density increases with the increase of rib width. The increase of current density with the rib width is quite large for the low conductivity of 280.6 S/m, but only moderate for the conductivity of 2806 and 28,060 S/m, indicating a conductivity of 2806 S/m is adequately high for the anode. Figure 5e shows the current density increases with the increase of O 2 electrode thickness, which is opposite to the trend of the cathode shown in Figure 4f. The opposite trends imply that the major polarization factors in the anode and cathode are different. The concentration loss is more influential than the ohmic polarization in the cathode, while the ohmic polarization is much larger than the concentration polarization in the anode. The thicker anode increases the difficulty of both the O 2 diffusion and the current conduction along the electrode thickness direction. However, the current passage through the narrow cross section of the electrode to the rib is the main ohmic polarization loss. The increase of electrode thickness can reduce the major ohmic loss and improves the SOEC performance, as observed experimentally [19]. For all the anode thicknesses considered, the current density increases with the increased O 2 electrode rib width, due to the reduced ohmic polarization by the shorter conduction path.
As shown in Figure 5f, fixing the rib width, the SOEC performance decreases with the increased anode ASR due to the increased ohmic loss. Meanwhile, the current density increases continuously with the increase of rib width, independent of the ASR value, confirming the concentration loss is far lower than the ohmic loss in the O 2 electrode. The optimal rib width can be quite close to the whole pitch size.

Analytical Expression of the Optimal Rib/Pitch Ratio on the H 2 Electrode Side
From Figure 6, it can be seen that with the increase of pitch width, the optimal rib width increases, while the optimal rib/pitch ratio (R) decreases. Furthermore, for a SOEC with pitch width 2 mm and ASR 0.01 Ωcm 2 , the optimal rib width is 0.387 and 0.603 mm when loaded with applied voltage 1.6 and 1.2 V, increasing by 56%, so the influence of the applied voltage on optimal rib width cannot be ignored. The specific voltage is meaningless, the possible influence factor is the total polarization loss η = Vcell−Enerst. However, here the explored total polarization change is induced only by applied voltage change. Combining the parameters sensitivity analysis above, it can be distinguished that ASR, polarization loss and pitch width from all the parameters studied as the major factors affecting optimal rib width. Here we take the rib/pitch ratio into consideration, and analytically express the relationship between optimal rib/pitch ratio with pitch width, rib-electrode ASR, total polarization loss. Finally, we get the expression for the applicable voltage range 1.2-1.6 V for SOEC. The optimal rib/pitch ratio is denoted as R 1.2V , R 1.4V , R 1.6V .  R can be expressed generally as R = (a 1 *ASR + a 2 ) + (a 3 *ASR + a 4 ). The total polarization for SOEC with applied voltage 1.2, 1.4, 1.6 is 0.34, 0.54, 0.74, respectively. The four sets of coefficients are listed in Table 2 and the relationship between the coefficients and η is displayed in Figure 7.  By simplifying the relationship between a 1 , a 2 , a 3 , a 4 and η, a 1 = 0.325η − 0.54523, a 2 = 0.0395η − 0.04719, a 3 = −2.06325η + 5.89582, a 4 = −0.33428η + 0.4488 are obtained to describe the relationship between η and a 1 , a 2 , a 3 , a 4 linearly. The synthetic expression of R considering the influence of pitch width, rib-electrode ASR, and the total polarization is formulated as: It should be pointed out that in Equation (41), the units of ASR, η, d pitch are Ωcm 2 , V, mm, respectively. Subsequently, the optimal rib can be formulated as: Here, d rib,a (mm) is the optimal rib width obtained analytically. It should be noted that even the d rib,a (mm) is the quadratic function of d pitch , the coefficient for the first order term is an order of magnitude larger than the quadratic term. Therefore, the relationship between optimal rib width and pitch width is mainly linear.
To verify the reliability of Equation (42), two accuracy parameters are defined as where i max-n is the maximum current density corresponding to the optimal rib width d rib-n obtained numerically. i max-a is the maximum current density for the analytically predicted optimal rib width d rib-a . λ 1 , λ 2 are defined to evaluate the effectiveness of the analytical expression for predicting optimal rib width. Even the maximum margin of error for optimal rib width (λ 2 ) is about 8.575%, the current density for the analytically predicted optimal rib is only 0.14% different from the maximum current density for the numerically attained optimal rib width. Hence, the cathode rib width predicted by Equation (42) is sufficiently optimal.

Conclusions
We have developed a comprehensive mathematical model for the performance research of SOEC. The impacts of the electrode rib widths on SOEC performance are systematically examined by varying the contact resistance, fuel composition, electrode porosity, electrode thickness and electrode electric conductivity. Different from rib width at the H 2 electrode side, the optimal rib width at the O 2 electrode side is not so sensitive to ASR and other parameters. With the increase of rib width, current density almost increases all the way, because the concentration polarization at the O 2 electrode side is too small, the ohmic loss induced by ASR is overwhelmingly greater than concentration loss. Nevertheless, the current density increase rate decreases when rib width is large enough. More importantly, it finds out three main factors affecting optimal rib width at the H 2 electrode side: ASR at the rib-electrode interface, the pitch width, and the applied cell voltage. Giving the pitch width and contact resistance, the optimal rib width decreases with the increase of cell voltage. Giving the cell voltage and contact resistance, the optimal rib width increases and the optimal rib/pitch ratio decreases with the increase of pitch width. When ASR at H 2 electrode side is in the range of (0.01, 0.05) Ωcm 2 , the optimal rib width at the cell voltage of 1.4 V is in the range (0.511, 0.840), (0.637, 1.065), (0.736, 1.248) and (0.818, 1.410) mm for the corresponding pitch width of 2, 3, 4 and 5 mm, respectively. Finally, an analytical expression is proposed to formulate the relationship between the optimal rib width and total polarization loss, ASR, and pitch width. The prediction error between the maximum current density obtained by numerically and analytically optimal rib width is within 0.14%, proving the predictive ability of the analytical expression.