1. Introduction
Fuel cell (FC) converts storable fossil fuels into instantly usable electricity in a higher efficiency and with CO
2 capture ready [
1], and also is attractive due to less energy losses during its energy conversion process [
2]. The solid oxide fuel cell (SOFC) is one of the fuel cells family and converts chemical energy directly into electrical energy [
3], resulting in considerable attention in both its engineering fabrications and theoretic modeling that provides the in-depth mechanical and practical insights for its speed-up development [
4]. Among them, the planar-design SOFC is simple in its fabrication and more flexible in its cell geometry and gas manifolding [
5]. The SOFC generally operate under elevated temperatures (typically 800 °C), resulting in considerable concerns on the durability of its electrode and electrolyte materials and further its higher costs in manufacturing and operation. Higher operational temperatures also make start-up and shutdown longer, less applicable in portable power and transportation markets, and less flexible in practical situations. These key technical issues addressed the critical demands on scientific and engineering research breakthroughs in balance parameters and performance of SOFC and the development of advanced materials.
The operation temperature of SOFC has to mainly equilibrate two opposing factors. One is a higher temperature favoring for an increased rate in the oxygen incorporation, reduction, and ionic conductivity, and the other is a lower temperature favoring for the chemical and mechanical durability and the rapidness of the expansion-contraction cycle. Therefore, all aspects addressed on the highest efficiency and durability of the cathode-electrolyte-anode combination at the lowest possible temperatures. This inspired the current considerable progress in searching and making proper material candidates for a shift of the higher operational temperature (800 °C–850 °C) down to an intermediate temperature (such as 600 °C–650 °C) [
1]. Practical development involved the use of electrolyte materials of high ionic conductivity and the decrease of electrolyte thickness [
6,
7]. The perovskite-type metal oxides are traditional materials in the cathode-electrolyte-anode combination with the mobile oxygen vacancy [
8,
9], showing the promising ionic conductivity, structural rigidness, simple preparation, and excellent match with each other, such as LSMs (La
1−xSr
xMnO
3-delta) [
10], LSCFs(La
1−xSr
xCo
yFe
1−yO
3-delta) [
11,
12,
13] and many promising new candidates [
14]. The current perovskite SOFC materials considered a toxic cobalt-free formula, such as La
0.5Sr
0.5Fe
0.9Nb
0.1O
3-delta (LSFNb) by a simple sol-gel method [
15] and La
0.25Sr
0.75TiO
3 (LST) as the electrolyte in the super conduction [
16].
Alternatively, the development of SOFC took advantage of mathematical models, via better understandings of reaction mechanisms of SOFC for its optimization in terms of cost, efficiency, and durability [
17]. A few examples covered the current status of the development mathematical models and their focuses, such as the functionality of the direct internal reforming inside the SOFC system [
18], a two-dimensional model on mass transport in the porous electrodes and fuel and air channels based on the lattice Boltzmann method (LBM) [
19] and an anode-supported SOFC button cell [
20]. Further, the more accurate three-dimensional mathematical models were developed to predict the key operating parameters on the performance of the planar type of SOFC and parasitic loads [
3], optimized interconnect designs [
21], different flow stack designs of the planar SOFC cells [
22], the examination of the thermo-fluid reacting environment and local thermodynamic state in SOFC stacks [
23]. A mathematical model was presented within the transmission network for dynamic stability analysis and control [
24]. A few started to extract information on the restricted step in electrochemical properties of SOFC, such as effects of the electrolyte thickness and the operating temperature, but limited in investigations on the heat and mass transfer characteristics of SOFC [
25]. A novel application of interior search optimizer was presented to define the necessary parameters to model SOFCs [
26].
In summary, there were few explorations deep into the most critical parameters on electrode and electrolyte materials themselves, such as combining the ionic conductivity of electrode and electrolyte materials and the activation energy of electrochemical reactions. In searching for dominating factors restricting the overall performance and identification of critical parameters upon variable operational temperatures in a planar SOFC, a three-dimensional SOFC model was established based on COMSOL Multiphysics 5.4, COMSOL AB, Stockholm, Sweden. The numerical model couples mass transfer and electrochemical reactions. It served as a simulation platform to explore detailed reactant profiles on the multispecies mass transportations in gas channels and porous electrodes, ionic conductions, and electrochemical reactions, and further their impacts on the overall performance of SOFC. Based on this, the most significant and restricting parameters upon the development of the intermediate-temperature SOFC were critically investigated and identified. The study sketched views on the significance of switching between transportations and reactions of reactants in SOFCs, as well as assisting the material development in a precise way.
3. Results and Discussion
3.1. The Model Validation
The numerical simulation was conducted to evaluate the effects of parameters on the thermo-fluidic-reaction characteristics of the unit cell stack of the SOFC, using the commercial finite element software COMSOL Multiphysics. The comparison of the performance curve determined a properly applied grid to efficiently reach accurate results.
The numerical simulations were initially carried out following the same cell designs and operating parameters in the published literature [
25,
29]. The numerical results on the cell voltage, power density, and their associated current densities based on our current model were compared to those of published experimental and numerical results, in
Figure 2.
The errors between the present numerical results and the experimental results were less than 7%, implying that the current model agrees well with both experimental and previously-published numerical results [
25,
29]. The minimal difference in the current model and experimental results was mainly derived from the assumption and simplification of models which has primarily been introduced in the literature [
25], such as an ideal transport among multilayers of SOFC and diversified sources of physical parameters of the porous electrodes [
36,
37]. The accuracy of the simulation would be expected upon the future readiness of experimental material characterizations.
3.2. The Temperature-Dependence of the Regular YSZ Electrolyte (the Absolutely Electrolyte-Dominating Case)
To realize the intermediate temperature operation of the SOFC, the significance of temperature effect on the performance of SOFC components was initially investigated, based on the established model with capability in the quantitative evaluation of all concerned parameters of SOFC. Taking the regular and most popular YSZ electrolyte as an example, the sensitivity of the YSZ electrolyte with the operation temperature can be calculated. It’s assumed that major parameters on fuel and oxidize channels and electrodes, such as their velocity and constituents, were kept constant.
Figure 3a exhibited the numerical results on the temperature dependence of the cell voltage, the power density, and their associated current densities following a trend of temperatures of 800 °C, 700 °C, and as low as 650 °C. The current density of the electrolyte for different temperatures was depicted
Figure 3b at an operating voltage of 0.8 V.
As shown in
Figure 3a, the decrease of the operational temperatures from 800 °C to 700 °C and 650 °C, the current density significantly dropped by 63% and 80% at a constant operational voltage, implying the dramatic diminution of the cell performance as a consequence of lowing operating temperature by only 100 °C and 150 °C, respectively. Further investigation revealed the important temperature dependence of the ionic conductivity of the electrolyte, such as a 79% loss of the ionic conductivity in the applied YSZ electrolyte as a temperature dropping from 800 °C to 650 °C. Meanwhile, for a similar reason, the reference exchange current density
i0 of both anode and cathode dropped by 93.3% and 92.9%, respectively. In many anode-supported SOFC applications, a thinner electrolyte was usually applied at lower temperature operations to decrease its ohmic overpotential. But notably, the unexpected performance was still its major issues [
5]. Considering a minimum thickness of about 10 mm of YSZ in terms of its acceptable structural strength, the limit to operating temperature can only be as low as 700 °C [
1]. Therefore, further dropping to a desirable temperature as 650 °C unavoidably depended on the development of a new electrolyte with a higher ionic conductivity. It should be noticeable that the current designs of electrodes in the SOFC was exclusively the active-material-on-electrolyte. Therefore, the temperature-dependent properties of electrolytes had great impacts on the electrode performance in such structured SOFC.
3.3. The Temperature-Dependent Ionic Conductivity in Novel Electrolytes
The ion conductivity in the electrolyte was specifically noticed on its crucial influence on the overall conductivity, considering that it’s significantly lower than the electric conductivities by several orders of magnitude. This was especially true when they were applied in the intermediate-temperature SOFC. Under this inspiration, a new group of novel electrolyte materials were developed, leading both the increase of their ionic conductivity as well as further understandings of the fundamental oxide ion conduction mechanisms. Various electrolytes, alternative to YSZ, have been prepared and investigated [
38], mainly two groups involving aliovalent-doped ceria and isovalent-cation–stabilized bismuth oxides, which were particularly attractive in their superior ionic conductivity and thus applicable in the SOFC at the intermediate temperatures. Among them, dysprosium- and tungsten-stabilized bismuth oxide (DWSB) has a closer inherent structure to the pure d-Bi
2O
3. Both lattice parameter and conductivity linearly extrapolated with the total dopant concentration in the pure d-Bi
2O
3, resulting in its ability to stabilize d-phase at its lower dopant concentration, and thus achieving higher conductivity [
39]. The singly doped Bi
2O
3s, Er
2O
3-stabilized Bi
2O
3 (Er
0.2Bi
0.8O
1.5, ESB) was known to have the highest oxygen ion conductivity in its lowest dopant (Er
3+) concentration required to achieve the cubic δ-phase stabilization [
40]. Gadolinium-doped ceria (Ce
0.9Gd
0.1O
2−δ, GDC) [
41] has a cubic structure and a density of around 7.2 g cm
−3 in its oxidized form. It is one of a class of ceria-doped electrolytes with higher ionic conductivity and lower operating temperatures (<700 °C). Optimization of dopant concentration in the Sm
x/2Nd
x/2Ce
1−xO
2−δ system resulted in the development of Sm
0.075Nd
0.075Ce
0.85O
2−δ (SNDC), which exhibited a higher grain ionic conductivity than that of GDC in air [
42].
Table 2 listed temperature-dependent ionic conductivity properties of these 4 novel electrolytes in open literature [
1], and YSZ as a referenced electrolyte.
To evaluate the intermediate temperature operation of the SOFC with these new electrolytes, the overall performance and current densities for the different electrolytes were discussed. In the process, Except for the electrolyte, the geometric and other physical model parameters were the same as depicted in
Section 2.1.
Figure 4 revealed the expected increase in the overall current density of SOFC upon the improvement of the ionic conductivity in applied electrolytes at the same operational voltage, but without consideration of activation energy variation, there was no straightforward linearity dependence between the overall current density of SOFC and the ionic conductivity the electrolyte.
Figure 5 exhibited the current densities of the different electrolytes when different electrolytes were applied at an operating voltage of 0.8 V and 700 °C.
For example, in a case of 700 °C and 0.8 V operational voltage, the correspondent current density were only 1.67, 1.66, 1.52, 1.47 times of that of YSZ as an electrolyte, versus aforementioned ionic conductivities of 4 electrolytes (DWSB, ESB, SNDC, GDC) at their 71.88, 40.38, 6.14, and 4.55 times of that of YSZ. This implied the improvement of the current density of the SOFC was very limited upon the significant improvement in ionic conductivity of applied electrolytes. Notably, the ionic conductivity of DWSB was 14.8 times higher than that of GDC, versus the difference in the current density by only 13.6%. This was also true when the operational temperature further dropped to 650 °C. The most significant improvement in the SOFC performance, in terms of the current density, was when the switch from YSZ to GDC electrolytes, implying that the significant impact of the ionic conductivity of the electrolyte only when it was lower. And the further improvement of the ionic conductivity should not contribute any more benefits of the SOFC performance in the current specific situation. This was a sign of the transition of restriction steps of the overall electrochemical process inside SOFC.
3.4. The Transition of the Dominating Factor
Actually, both the ohmic overpotential and the activation overpotential contributed to the overall SOFC performance. The ohmic overpotential consists of all electric conductions including the aforementioned ionic conductivity and electronic conductivity. The activation overpotential was mainly relative to electrochemical reactions of reactants in electrodes of SOFC. At a constant temperature and same electrode materials, the activation overpotential can be regarded as equivalent, and thus the only ohmic overpotential can be focused. The Ohmic overpotential (
ηohm) is proportional to the current density
i as given by:
where
Rohm is the overall internal resistance of the cell that can be estimated from the effective distance between the cell components coupled with conductivity data, as given by [
43,
44]:
where
τ is the thickness of each cell component, and
σanode,
σcathode, and
σinterconnect are the electronic conductivity of the anode, cathode and interconnect, and
σelectrolyte is the ionic conductivity of the electrolyte.
In our current calculations, the electronic and ionic current densities were both governed by the charge continuity equations involving both the ionic conducting phase and the electronic conducting phase. Thus, the analysis on the Ohmic overpotential for the calculation of the internal resistance was based on Equation (14).
Table 3 lists the ratio for each selected electrolyte. This ratio was defined as a proportion between the electrical resistance of the applied electrolyte and the overall internal resistance of the cell unit (involving electrical resistance of the applied interconnect). The larger the ratio, the greater the contribution of electrical resistance of the applied electrolyte toward the overall internal resistance of the cell unit.
Table 3 indicated that the YSZ electrolyte has a large aforementioned ratio, and lower operational temperature lead to a larger ratio, such as 89.3% at 650 °C and 83.4% at 700 °C, which both were over 80%. Under this case, the improvement in the ionic conductivity of the electrolyte should be significant in the enhancement of the SOFC performance, just as that shown in
Section 3.3 when the regular YSZ electrolyte was switched with the GDC electrolyte. However, this trend was not continued as the ratio dropped to about 60%, as referring to the case of the switch from the GDC electrolyte to either SNDC, ESB, and DWSB, respectively, which corresponded to less effective improvement of the SOFC performance. This coincided with ratios of other three advanced electrolytes with excellent ionic conductivities, such as ratios of the electrical resistance in the internal resistance of the cell dropped to 64.7%, 57.6%, 17.1%, 10.4% at 650 °C for GDC, SNDC, ESB, and DWSB, respectively. These ratios were further dropped upon temperature toward 700 °C. The lower the operational temperature, the higher critical ratio of the electrical resistance of the electrolyte in the overall internal resistance of the cell. Alternatively, the cell voltage and the power density with their associated current densities can be calculated, if the increase of the ionic conductivity of the YSZ electrolyte can be assumed to be 2, 5, and 10 times of the ionic conductivity of the regular YSZ electrolyte (2.2669 S m
−1).
Figure 6 exhibits the comparison results on the cell voltage and the power density with their associated current densities in these cases. It was found the current density was enhanced by only 57.8%, 48.3%, and 25.6% at the same cell voltage, noticing that the ionic conductivity has been greatly increased by 900%, 400%, and 100%. A similar result was obtained, especially the move of 400% to 900% in the ionic conductivity to validate the critical value of the ratio of the electrical resistance.
There indeed existed a critical value of the addressed ratio, referring to the electrical resistance of the electrolyte in the overall internal resistance of the cell. The corresponding enhancement of the SOFC performance to the improvement only occurred when the ionic conductivity of electrolytes accounted for a substantial part of the electrical resistance. If the increase of ionic conductivity of the YSZ substrate in anode and cathode are considered, the current density can be increased by a further 15%, 43%, and 64%, respectively, corresponding to assuming ionic conductivity of YSZ 2, 5, and 10 times greater. This well explained the mechanism on the effect of the operational temperature on the overall performance of SOFC. This would be a moment on the switch of restriction steps of the electrochemical process in the SOFC.
3.5. The Switch of the Restriction Steps of Electrochemical Reactions in SOFC
Since the switch of the restriction steps of electrochemical reactions in the SOFC was caused by the temperature, other parameters of temperature dependence were further investigated. It was noticed that the activation overpotential became significantly lower to cause a cell voltage drop when the operating temperature decreased (such as 650 °C), although it was insignificantly small at higher temperatures, such as 800 °C [
5]. The activation overpotential (
ηact) refers to the restriction effect of chemisorption and electrochemical activation of gas reactants on electrode materials. The higher the activation overpotential, the lower reaction rates of electrochemical reactions in SOFC. This implied that the drastic degradation of the cell performance at the decrease of operating temperatures can also be caused by an increase in not only the ohmic overpotential, also the activation overpotential. The former seemed less significant to restrict SOFC performance than the latter.
The generally-accepted Butler–Volmer equation (Equation (6)) was applied to describe the activation overpotential in this study, and the relevant reference exchange current density
i0 was defined by Equation (7), in which the temperature-dependent
i0 was solely attributed to the activation energy. It’s noticeable that the activation energy of the regular 800 °C-adaptable electrode materials increases upon the operational temperature decreases. The activation energy was assumed to be both a 10% increase upon the operational temperatures dropping to 700 °C and 650 °C, respectively using YSZ as the electrolyte. The calculated
i0 were found to be 82.3% and 81.6% lower in the anode and cathode, respectively, and 83.9% and 83.2%, respectively at further lower to be 650 °C. The increase of the activation energy by only 10% significantly decreased
i0.
Table 4 lists the calculated overall cell current density when temperatures were at 700 °C and 650 °C at the cell voltage of 0.8V, covering five different electrolytes. For the regular YSZ, the drops of the overall cell current density were 42.8% and 49% at temperatures of 700 °C and 650 °C, respectively, almost 50% lower. Switching to GDC electrolyte, of which ionic conductivity was several times higher than that of YSZ at the same temperature, only 52.1% and 58.4% of the overall cell current density dropped. This trend was followed by the SNDC, ESB, and DWSB, where ionic conductivities were much higher than that of YSZ. It implied the electrochemical reactions dominated the overall SOFC performance, rather than the ionic conduction. The activation energy of electrode materials exhibited more sensitivity to the temperature than the ionic conductivity of electrolyte to the temperature.
It’s important to decrease the activation energy of the electrode material during the development of novel advanced materials. A similar calculation was conducted, but the activation energy was decreased by 10% upon the decrease of the temperature, to show its significance of the development of currently advanced electrode materials.
Table 5 exhibits the increased multiple of calculated reference exchange current density
i0 versus the change of the activation energy, referenced by
i0 at the true activation energy.
Table 6 exhibits the increased proportion of current density of SOFC cell unit when a 10% decrease of the activation energy at an operating voltage of 0.8 V.
Table 5 indicates the significant impact of the activation energy on the reference exchange current density
i0. The increase of only 10% activation energy can increase
i0 in 380%, 464%, and 520% at 800 °C, 700 °C and 650 °C in anode, and similarly 364%, 444%, and 496% in the cathode. Furthermore, it was interested to find the decrease of the activation energy by only 10% increased the overall cell current density by 27%, 43%, and 54% at temperatures of 800 °C, 700 °C, and 650 °C, respectively, in the case of the regular YSZ electrolyte material at the same cell voltage. Similar to the case of the ionic conductivity, the change of the overall current density was not linearly proportional to the temperature of the SOFC.
Table 6 indicated that switching to other electrolytes with increased ionic conductivities, the cell current density seemed not to change this trend, either. Moreover, the cell current density was independent from the variation of electrolytes, which was solely attributed to the activation energy of electrochemical reactions.
4. Conclusions
This study established a three-dimensional SOFC model based on the Comsol software package, focusing on the identification of critical parameters upon variable operational temperatures in a planar SOFC. The numerical results were in good agreement with the available data in the literature, and accuracy and reliability of the numerical model were validated. The current studies specially addressed major aspects of restricting parameters upon the development of the intermediate-temperature SOFC.
In searching for dominating factors restricting the overall SOFC performance, the ionic conductivity of the electrolyte was primarily focused, which was critical not only as the electrolyte itself but also the substrates applied for electrodes. The temperature-dependent nature of the ionic conductivity of electrolyte materials was identified to be unavoidably concerned in the development of the intermediate-temperature SOFC, because the current density can be 63% and 80% decreases when the operational temperature dropped to the ideal 700 °C and 650 °C from the regular 800 °C. However, the increase of the ionic conductivity of concerned materials was not linearly proportional to the current density of the SOFC. There identified a critical value on the ratio of electric resistance of the electrolyte in the overall internal resistance of the SOFC, above which the further increase in the ionic conductivity would not significantly improve the SOFC performance in terms of current density at a constant voltage.
The lower the operational temperature, the higher critical ratio of the electrical resistance in the overall internal resistance of the cell. The minimal decrease in the activation energy during the development of intermediate electrode materials can significantly enhance the overall SOFC cell performance, implied the restricting step in the whole electrochemical reactions in SOFC was the activation energy of electrode material. Considering the development trend toward the intermediate temperature SOFC, advanced electrode material with the decreased activation energy should be primarily focused. It’s all about the switch of restriction of steps in the whole electrochemical reactions in SOFC. A simulation model exhibited the potentials to assist the material development in a precise way.