In this study, a numerical approach is adopted to investigate the effectiveness factors for distributed electrochemical reactions in thin active reaction layers of solid oxide fuel cells (SOFCs), taking into account the ButlerVolmer reaction kinetics. The mathematical equations for the electrochemical reaction and charge conduction process were formulated by assuming that the active reaction layer has a small thickness, homogeneous microstructure, and high effective electronic conductivity. The effectiveness factor is defined as the ratio of the actual reaction rate (or equivalently, current generation rate) in the active reaction layer to the nominal reaction rate. From extensive numerical calculations, the effectiveness factors were obtained for various charge transfer coefficients of 0.30.8. These effectiveness data were then fitted to simple correlation equations, and the resulting correlation coefficients are presented along with estimated magnitude of error.
1. Introduction
Solid oxide fuel cells (SOFCs) are a promising hightemperature fuel cell technology that can efficiently convert the chemical energy of fuel into electricity
[1
,
2]
. The high operating temperature (~8001000℃) of SOFCs reduces manufacturing costs by employing nonprecious metal catalyst, directly uses hydrocarbon fuels by utilizing internal reforming, and achieves higher overall efficiency through the construction of hybrid or combined heat and power (CHP) systems. However, the high operating temperature of SOFCs tends to accelerate the microstructural degradation of electrodes, thereby negatively affecting the longterm performance stability. Thus, intermediatetemperature SOFCs (ITSOFCs) that can operate at ~650800℃ have emerged as a solution to the longterm performance degradation problem
[3

5]
.
To achieve better performance of ITSOFCs, material scientists have focused on the development of new materials, such as electrolyte materials with higher ionic conductivities and electrode materials with higher electrochemical activities, at these intermediate temperature ranges
[6

8]
. Concurrently, other researchers have attempted to determine the optimal electrode structure that can maximize the electrode performance for a given material. For this purpose, multilayer
[9

12]
, and functionally or microstructurally graded electrode structures
[13

15]
have been proposed to enhance the electrochemical performance of ITSOFCs. In addition, recent experimental studies have directly reconstructed threedimensional electrode microstructures using advanced imaging techniques, through which the relationship between the electrode microstructural parameters and electrochemical efficiencies has been investigated
[16

18]
.
In multilayer electrodes of ITSOFCs, a thin layer dedicated to the electrochemical reactions is in direct contact with the electrolyte. This layer is called the active reaction layer and is made by mixing fine electronic and ionic conductor particles, to provide rich threephase boundaries (TPBs) for electrochemical reactions. In the exterior of the active reaction layer, a bulk transport layer is formed to facilitate fast mass transport and electron conduction. Almost all electrochemical reactions occur inside the active reaction layer, which highlights that microstructural optimization of the active reaction layer is essential for improved ITSOFC performance. In this regard, the electrochemical effectiveness factor has been investigated as an appropriate measure for evaluating the efficiency of electrodes.
Previously, Costamagna et al.
[19
,
20]
proposed the electrochemical effectiveness concept in search of the optimal microstructure for SOFC electrodes. However, their model is based on the linear transfer currentlocal overpotential (
i

η
) relationship, and thus is valid only for very low current density (or equally, for very low overpotential) conditions. Recently, Shin and Nam
[21]
resolved this limitation by considering a nonlinear
i

η
relationship, the symmetric ButlerVolmer reaction kinetics. They showed that the electrochemical effectiveness factor can be decomposed into the base effectiveness at zero overpotential and the relative effectiveness at finite overpotentials. In addition, a simple correlation equation was also proposed along with the correlation coefficients relevant to the symmetric ButlerVolmer reaction kinetics. Their model was successfully used for onedimensional simulation of a singlecell SOFC
[22]
and the theoretical prediction of electrode microstructural effects
[23
,
24]
.
In this study, the work of Shin and Nam
[21]
is further extended by considering both the symmetric and asymmetric ButlerVolmer reaction kinetics. As before, the active reaction layer is assumed to be homogeneous in microstructure and sufficiently thin to ensure uniform conditions (temperature, pressure, and species concentration). In addition, the electronic potential is assumed to be uniform inside the reaction layer, since the effective electronic conductivity is much higher that the effective ionic conductivity. The effectiveness factors for current generation performance (or equally, electrochemical reaction efficiency) in the active reaction layer were numerically obtained by varying the charge transfer coefficient from 0.3 to 0.8. Finally, the calculated effectiveness data were fitted to the correlation equation proposed by Shin and Nam
[21]
, from which the correlation coefficients and estimated errors were determined.
2. Theory and Calculations
 2.1 Physical model and assumptions
The fuel cell processes inside the twolayer anode of ITSOFCs are illustrated in
Fig. 1
, where the species transport, charge conduction, and distributed electrochemical reactions are indicated. The active reaction layer in the anode is very thin (~20 μm) and is composed of fine electronic (Ni) and ionic (YSZ; yttriastabilized zirconia) conductor particles (whose mean particle diameter is ~0.5 μm) to provide rich TPB sites for electrochemical reactions. In contrast, the bulk transport layer in the anode is relatively thick and is made of coarse particles to enable fast mass transport (also to provide structural support in the case of anodesupported ITSOFCs). In
Fig. 1
, the electronic current,
I
_{el}
, ionic current,
I
_{io}
, and charge transfer current,
I
_{tr}
, are also illustrated to explain the distributed nature of electrochemical reactions in the active reaction layer. In the anode, the electronic current,
I
_{el}
, decreases in magnitude as it flows towards the electrolyte in direct proportion to the charge transfer current,
I
_{tr}
(or equally, the electrochemical reaction rate). Accordingly, the ionic current,
I
_{io}
, increases in magnitude by collecting
I
_{tr}
, and flows through the electrolyte and towards the cathode.
Electrochemical reactions and transport processes (mass diffusion and charge conduction) in the twolayer anode of ITSOFCs, along with the resultant electronic current, I_{el}, ionic current, I_{io}, and charge transfer current, I_{tr}.
The conservation of electron, oxygen ion, and gas species, along with the electrochemical reactions, should be considered in the simulation of fuel cell processes in the active reaction layer. In this study, it is assumed that the active reaction layer is homogeneous in microstructure, sufficiently thin to ensure uniform operating condition (temperature, pressure, and species concentration), and has a significantly higher effective electronic conductivity,
σ
_{el,eff}
, compared to its effective ionic conductivity,
σ
_{io,eff}
. These are believed to be practical assumptions that lead to negligible errors in the results
[21
,
22]
. The physical interpretation of each assumption can be summarized as follows.
• Homogeneous microstructure: The volumespecific TPB length (TPBL),
λ
_{tpb,V}
, and the effective ionic conductivity,
σ
_{io,eff}
, are uniform inside the active reaction layer.
• Uniform operating condition: The TPBLspecific exchange current density,
i
_{tpb}
, is uniform inside the active reaction layer. In addition, the Nernst potential,
φ
^{o}
, and the concentration overpotential,
η
_{conc}
, are also uniform.
• High effective electronic conductivity: The electronic potential,
φ
_{el}
, is relatively uniform inside the active reaction layer, compared with the variation of the ionic potential,
φ
_{io}
.
 2.2 Governing equations and reaction kinetics
The governing equations and boundary conditions for electronic and ionic charge conservation can be expressed for the active reaction layer in the anode of ITSOFCs (shown in
Fig. 1
) as
[20
,
21]
where
i
_{el}
is the electronic current density,
i
_{io}
is the ionic current density, and
i
_{tr,V}
(
η
) is the volumetric charge transfer current density at the local activation overpotential of
η
. In addition,
φ
_{el,0}
denotes the electronic potential at
z
= 0 and
φ
_{io,L}
denotes the ionic potential at
z
= 0.
In this study,
φ
_{el}
is assumed to be uniform inside the active reaction layer, such that
φ
_{el}
(
z
) =
φ
_{el,0}
; thus, it is not necessary to solve Eqs. (1) and (2). It should be noted that the local activation overpotential,
η
, is defined as
where
φ
^{o}
is the Nernst potential and
η
_{conc}
is the concentration overpotential. In Eq. (5),
φ
_{el}
,
φ
^{o}
, and
η
_{conc}
are all constant inside the active reaction layer according to the assumptions of this study. Then, the governing equations and boundary conditions for ionic charge conservation, Eqs. (3) and (4), can be expressed in terms of
η
as
Here,
η
_{conc}
is the total activation overpotential applied on the active reaction layer, determined as
Eq. (6) can be alternatively obtained by subtracting Eq. (3) from Eq. (1), based on the relation ▽
^{2}
η
= ▽
^{2}
φ
_{el}
– ▽
^{2}
φ
_{io}
[20]
. In this case,
σ
_{io,eff}
in Eq. (6) should be replaced with the effective charge conductivity,
σ
_{eff}
, defined as
σ
_{eff}
= (1/
σ
_{io,eff}
+ 1/
σ
_{el,eff}
)
^{1}
. When
σ
_{el,eff}
is much larger than
σ
_{io,eff}
, as considered in this study,
σ
_{eff}
is approximately equal to
σ
_{io,eff}
. It should be noted that Eqs. (6) and (7) are valid for both the anode and cathode reaction layers.
According to the ButlerVolmer equation, the volumetric charge transfer current density,
i
_{tr,V}
(
η
), is expressed as
where
α
is the charge transfer coefficient,
n
_{e}
is the number of electrons (2 for the anodic reaction and 4 for the cathodic reaction),
F
is the Faraday constant (96,485 C mol
^{1}
),
R
is the universal gas constant (8.314 J mol
^{1}
K
^{1}
), and
T
is the temperature in Kelvin. In Eq. (9), the volumetric exchange current density,
i
_{ex,V}
, can be obtained by multiplying the TPBLspecific exchange current density,
i
_{tpb}
, and the volumespecific TPBL,
λ
_{tpb,V}
as
i
_{ex,V}
=
i
_{tpb}
×
λ
_{tpb,V}
. Note that Eq. (9) may also be expressed in terms of the TPBLspecific polarization resistance,
r
_{tpb}
(Ωm), as
Fig. 2
shows the charge transfer current density vs. local activation overpotential (
i
_{tr,V}

η
) curves for the ButlerVolmer reaction kinetics (
n
_{e}
= 2 and
T
= 1000 K are assumed) determined using Eqs. (9) and (10). A larger charge transfer coefficient,
α
, leads to a smaller activation overpotential for a given transfer current density. As shown in
Fig. 2
, the ButlerVolmer equation is reduced to a single linear polarization equation at very small overpotentials, irrespective of the charge transfer coefficient. The asymptotic expansion of Eqs. (9) and (10) results in
Charge transfer current density vs. local activation overpotential (i_{tr,V}η) curves for general ButlerVolmer reaction kinetics of Eq. (9) with various charge transfer coefficients. These curves are obtained for anodic reaction (n_{e} = 2) at T = 1000 K.
where
r
_{p,V}
is the volumespecific linear polarization resistance (
r
_{p,V}
=
r
_{tpb}
/
λ
_{tpb,V}
).
The charge transfer coefficient is an indicator of the symmetry of the activation energy barrier when a positive or negative overpotential is applied
[25]
. In SOFC modeling studies, the coefficient is generally assumed to be 0.5 (symmetric ButlerVolmer equation), primarily due to the lack of experimental data
[26]
. However, detailed electrochemical kinetic studies
[26

29]
on Ni pattern anodes or Ni cermet anodes has shown that the charge transfer coefficient can have values different from 0.5. For example, the charge transfer coefficient,
α
, for hydrogen oxidation on Ni/YSZ anodes was estimated to be 0.60.7 by Utz et al.
[26]
and approximately 0.7 by Holtappels et al.
[29]
. Thus, the development of the general electrochemical effectiveness model for both the symmetric and asymmetric ButlerVolmer reaction kinetics is important.
 2.3 Effectiveness model and correlations
The effectiveness factor,
Γ
_{eff}
, has been defined as the ratio of the actual current generation rate (or equally, the electrochemical reaction rate) in the active reaction layer to the maximum current generation rate as
[20
,
21]
In Eq. (12), the actual current density,
i
_{real,A}
, is calculated by integrating
i
_{tr,V}
(
η
) in the active reaction layer (0 ≤
z
≤
L
), while the maximum current density,
i
_{max,A}
, is directly obtained by multiplying
i
_{tr,V}
(
η
_{tot}
) and the layer thickness,
L
. Note that the maximum current density is obtained when all the TPB sites inside the active reaction layer are uniformly subject to
η
_{tot}
. In real situations, the local overpotential,
η
, is usually smaller than
η
_{tot}
(
η
≤
η
_{tot}
) inside the reaction layer and thus the effectiveness factor,
Γ
_{eff}
, is always equal to or smaller than 1 (0 ≤
Γ
_{eff}
≤ 1).
It was shown by Costamagna et al.
[20]
that the effectiveness factor for a linear charge transfer reaction, such as Eq. (11), can be expressed as
where
ϕ
_{T}
is the electrochemical Thiele modulus, defined as
Note that
Γ
_{eff}
in Eq. (13) has the same functional form as the effectiveness factor for chemical reaction/mass transfer in heterogeneous catalysis, or the fin efficiency for conduction/convection heat transfer in extended surfaces (fins). The Thiele modulus,
ϕ
_{T}
, usually has a value in the range of 520 for the anode reaction layers and 0.52.0 for the cathode reaction layers, for ordinary operation of ITSOFCs
[21
,
22]
.
Shin and Nam
[21]
studied the effectiveness factor for the symmetric ButlerVolmer reaction kinetics, such as Eqs. (9) and (10) with
α
= 0.5. From extensive numerical calculations, they showed that the effectiveness factor,
Γ
_{eff}
, for nonlinear reaction kinetics can be decomposed into two parts, namely, the base effectiveness at zero activation overpotential,
Γ
_{eff,0V}
, and the relative effectiveness at finite activation overpotential,
, as
Here,
is the dimensionless total activation overpotential applied on the active reaction layer
. In Eq. (15),
depends on
, not separately on
α
,
n
_{e}
,
T
or
η
_{tot}
. It is also shown that
has a nearly constant functional form with respect to
for
ϕ
_{T}
≥ 3.
It is not convenient to present all the effectiveness factors for various
ϕ
_{T}
and
in a tabulated form because of their large data size. Thus, Shin and Nam
[21]
proposed a simple correlation equation for easy and accurate determination of the electrochemical effectiveness data, which is written as
where
a
,
b
,
c
, and
d
are the correlation coefficients dependent on
ϕ
_{T}
. Eq. (16) accurately describes the behavior of the relative effectiveness,
, starting from 1.0 at a very small
ϕ
_{T}
, and decreasing towards 0.0 as
ϕ
_{T}
increases.
Using the electrochemical effectiveness factor, the current generation in the active reaction layer can be determined according to the following steps. First, the TPBLspecific exchange current density,
i
_{tpb}
, or the TPBLspecific linear polarization resistance,
r
_{tpb}
, is calculated for given operating conditions. Second, the electrochemical Thiele modulus,
ϕ
_{T}
, is calculated using Eq. (14). Third, the base effectiveness,
Γ
_{eff,0V}
, and the relative effectiveness,
, are calculated using Eqs. (15) and (16), from which the electrochemical effectiveness factor,
Γ
_{eff}
, is determined. Once
Γ
_{eff}
is known, the actual current density,
i
_{real,A}
, generated in the active reaction layer is obtained as
 2.4 Numerical calculations
The governing equations and boundary conditions provided in Eqs. (6), (7), and (10) were discretized using the finite difference method (FDM). Uniform 2000 grid points were placed in the onedimensional calculation domain of the active reaction layer (0 ≤
z
≤
L
). The discretized nonlinear algebraic equations were solved using the Engineering Equation Solver (EES) software
[30]
. Once the converged solution is obtained, the effectiveness factor was decomposed into the base effectiveness,
Γ
_{eff,0V}
, and the relative effectiveness,
. In this study, numerical calculations were performed for both the symmetric and asymmetric ButlerVolmer reaction kinetics, by varying the charge transfer coefficient,
α
, as 0.3, 0.4, 0.5, 0.6, 0.7, and 0.8.
3. Results and Discussion
 3.1 Symmetric ButlerVolmer reaction kinetics
Previously, Shin and Nam
[21]
investigated the electrochemical effectiveness factors for current generation in the active reaction layers, subject to the symmetric ButlerVolmer reaction kinetics (
α
= 0.5). The results and the accuracy of the electrochemical effectiveness model are briefly reviewed in this section.
Fig. 3
presents the relative effectiveness factor,
, which was numerically determined by solving Eqs. (6), (7), and (10) with
α
= 0.5 (symmetric ButlerVolmer equation). In
Fig. 3
,
generally decreases from 1.0 towards 0.0 as increases. For a given
,
becomes smaller as
ϕ
_{T}
increases. It is interesting to note that
seems to have a constant functional shape when
ϕ
_{T}
is higher than 3 (
ϕ
_{T}
≥ 3). Shin and Nam
[21]
showed that this behavior occurs when the active reaction layer thickness,
L
, becomes sufficiently thick, such that further increase of
L
does not contribute to the enhancement of current generation.
Relative effectiveness factor, , for the current generation performance of the active reaction layer in ITSOFCs: Symmetric ButlerVolmer reaction kinetics with α = 0.5.
In
Fig. 3
, “symbols” denote the relative effectiveness data obtained by numerical calculation, while “lines” denote the curves that are fitted to the correlation equation of Eq. (16). The correlation coefficients,
a
,
b
,
c
, and
d
, are listed in
Table 1
, along with the estimated errors. The results clearly indicate that using Eq. (16) for estimating
results in a small error (less than 1%). Thus, Eq. (16) and
Table 1
can be viewed as a complete solution for the electrochemical effectiveness for symmetric ButlerVolmer reaction kinetics, described by Eqs. (6), (7), and (10).
Correlation coefficients to determinefor symmetric ButlerVolmer reaction kinetics withα= 0.5.
^{a}Err_{0.3V} refers to the maximum correlation error in the estimated effectiveness factor for the cathode total overpotential range of 0.00.3 V at T=1000 K (or equally, for the anode total overpotential range of 0.00.6 V).
The accuracy of the electrochemical effectiveness factors provided in
Fig. 3
and
Table 1
was fully validated by Shin and Nam
[21]
. Excellent agreement was observed between the currentoverpotential performance curves obtained by the effectiveness model and those obtained by the detailed electrode microscale model. In addition, a onedimensional simulation model was also developed to predict the currentvoltage performance curves of a singlecell SOFC with twolayer electrodes (anode and cathode)
[22]
. The results obtained were essentially the same as the more detailed comprehensive microscale model results
[12]
. Thus, the electrochemical effectiveness factors presented in this study enable efficient calculation of the current density in the anode and cathode functional layers of SOFCs, without addressing the detailed electrochemical reaction/charge transport processes therein.
 3.2 Asymmetric ButlerVolmer reaction kinetics withα> 0.5
To extend the work of Shin and Nam
[21]
, extensive numerical calculations were performed to determine the effectiveness factors for current generation subject to asymmetric ButlerVolmer reaction kinetics with
α
> 0.5 (see
Fig. 2
).
Fig. 4
show the relative effectiveness data (symbols) and correlation equations (lines) for
α
= 0.6, 0.7, and 0.8. The general behavior of
for
α
> 0.5 shown in
Fig. 4
is similar to the behavior for
α
= 0.5 shown in
Fig. 3
, except that the range of
increases with
α
according to
=
αn
_{e}
F
η
_{tot}
/
RT
. The correlation coefficients are provided in
Tables 2
,
3
, and
4
. These tables summarize the maximum correlation errors encountered for the cathode total overpotential in the range of 00.3 V (equivalently, the anode total overpotential range of 00.6 V) at
T
= 1000 K.
Relative effectiveness factor, , for the current generation performance of the active reaction layer in ITSOFCs: Asymmetric ButlerVolmer reaction kinetics with (a) α = 0.6, (b) α = 0.7, and (c) α = 0.8.
Correlation coefficients to determinefor asymmetric ButlerVolmer reaction kinetics withα= 0.6.
Correlation coefficients to determine for asymmetric ButlerVolmer reaction kinetics with α = 0.6.
Correlation coefficients to determinefor asymmetric ButlerVolmer reaction kinetics withα= 0.7
Correlation coefficients to determine for asymmetric ButlerVolmer reaction kinetics with α = 0.7
Correlation coefficients to determinefor asymmetric ButlerVolmer reaction kinetics withα= 0.8.
Correlation coefficients to determine for asymmetric ButlerVolmer reaction kinetics with α = 0.8.
Relatively good agreement between the numerical effectiveness data denoted by “symbols” and the correlation equations denoted by “lines” is observed in
Fig. 4
. In
Tables 2

4
, the maximum correlation error, Err
_{0.3V}
, for the cathode total overpotential in the range of 00.3 V at
T
= 1000 K, is smaller than 0.6% for
α
= 0.6, smaller than 0.5% for
α
= 0.7, and smaller than 1.7% for
α
= 0.8. This result indicates that the proposed correlation equation of Eq. (16) accurately describes the behavior of
relevant to asymmetric ButlerVolmer reaction kinetics with 0.5 <
α
≤ 0.8.
 3.3 Asymmetric ButlerVolmer reaction kinetics withα< 0.5
Numerical calculations were also conducted to determine the electrochemical effectiveness for electrochemical reactions in the thin active reaction layers when subject to asymmetric ButlerVolmer reaction kinetics with
α
< 0.5 (see
Fig. 2
).
Figs. 5(a)
and
5(b)
show the relative effectiveness data (symbols) and correlation equations (lines) for
α
= 0.4 and 0.3, respectively, while the correlation coefficients are provided in
Tables 5
and
6
.
Relative effectiveness factor, , for the current generation performance of the active reaction layer in ITSOFCs: Asymmetric ButlerVolmer reaction kinetics with (a) α = 0.4 and (b) α = 0.3.
Correlation coefficients to determinefor asymmetric ButlerVolmer reaction kinetics withα= 0.4.
Correlation coefficients to determine for asymmetric ButlerVolmer reaction kinetics with α = 0.4.
Correlation coefficients to determinefor asymmetric ButlerVolmer reaction kinetics withα= 0.3.
Correlation coefficients to determine for asymmetric ButlerVolmer reaction kinetics with α = 0.3.
Fig. 5
shows that the relative effectiveness,
, has values higher than 1.0 for small
ranges, which is different from the behavior of
for symmetric (
α
= 0.5) and asymmetric (
α
> 0.5) ButlerVolmer reaction kinetics shown in
Figs. 3
and
4
, respectively. This trend can be explained by the ButlerVolmer reaction kinetics curves shown in
Fig. 2
. The ButlerVolmer equation with
α
≥ 0.5 always results in a higher transfer current density at a given overpotential, compared with the linear polarization equation. The relative effectiveness,
, is equal to 1.0 when the transfer current density vs. local activation overpotential (
i
_{tr,V}

η
) curve strictly follows the linear polarization relationship. The deviation of
i
_{tr,V}

η
curves with
α
≥ 0.5 from the linear polarization relationship leads to
which is smaller than 1.0 in
Figs. 3
and
4
. In contrast, the ButlerVolmer equation with
α
< 0.5 has a low overpotential region where the transfer current density is smaller than the linear polarization relationship (for
η
in the range of 00.12 V with
α
= 0.4 and for
η
in the range of 00.26 V with
α
= 0.3 in
Fig. 2
). Thus, the relative effectiveness,
, becomes larger than 1.0 in those low overpotential regions (for
in the range of 01.1 with
α
= 0.4 and for
in the range of 01.8 with
α
= 0.3 in
Fig. 5
).
In this study, the correlation equation of Eq. (16) was replaced with the following equation for
α
< 0.5 to describe the observed behavior of
greater than 1.0. The lines in
Fig. 5
represent the fitting curves of Eq. (18).
In
Tables 5
and
6
, the maximum correlation error, Err
_{0.3V}
, is smaller than 3.5% for
α
= 0.4 but as large as 13.9% for
α
= 0.3. As observed in
Fig. 5
, most correlation errors occur at low dimensionless total overpotential,
, correspond to the cathode total overpotential range of 00.05 V or 00.08 V at
T
= 1000 K. For
higher than these ranges, the correlation errors are relatively small as shown in
Fig. 5
. Thus, it may be necessary in future studies to develop a new correlation equation, other than Eq. (18), which is more appropriate for describing the electrochemical effectiveness for current generation in the active reaction layer subject to asymmetric ButlerVolmer reaction kinetics with
α
< 0.5.
4. Conclusions
In this study, extensive numerical calculations were performed to obtain the electrochemical effectiveness factors for current generation in the active reaction layer of SOFCs. Both the symmetric and asymmetric ButlerVolmer equations were considered to study the charge transfer process by varying the charge transfer coefficient from 0.3 to 0.8. Simple correlation equations were proposed for easy and accurate determination of the numerically determined effectiveness data, and the corresponding correlation coefficients and their estimated errors were summarized. It is anticipated that these results will prove to be useful for estimating current generation in the anodes and cathodes of SOFCs.
Nomenclature
a, b, c, dCorrelation coefficients FFaraday constant (96,485 C mol^{1}) Relative effectiveness factor at finite overpotential of i_{el}Electronic current density (A m^{2}) i_{ex,V}Volumespecific exchange current density (A m^{3}) i_{io}Ionic current density (A m^{2}) i_{max,A}Maximum current density in the active reaction layer (A m^{2}) i_{real,A}Actual current density in the active reaction layer (A m^{2}) i_{tpb}TPBLspecific exchange current density (A m^{1}) i_{tr,V}Volumetric charge transfer current density (A m^{3}) LActive reaction layer thickness (m) n_{e}Number of electrons for electrochemical reactions (2 for anode, 4 for cathode) RUniversal gas constant (8.314 J kg^{1} K^{1}) r_{p,V}Volumespecific linear polarization resistance (Ω m^{3}) r_{tpb}TPBLspecific linear polarization resistance (Ω m) TTemperature (K) zCoordinate (m) αCharge transfer coefficient ϕ_{T}Thiele modulus Γ_{eff}Electrochemical effectiveness factor Γ_{eff,0V}Base effectiveness factor at zero overpotential ηActivation overpotential (V) η_{conc}Concentration overpotential (V) η_{tot}Total activation overpotential applied to the active reaction layer (V) Dimensionless total activation overpotential, ≡ αn_{e}Fη_{tot} ⁄ / (RT) φ^{o}Nernst potential (V) φ_{el}Electronic potential (V) φ_{io}Ionic potential (V) λ_{tpb,V}Volumespecific TPBL (m m^{3}) σ_{eff}Effective charge conductivity, σ_{eff} = (1/σ_{io,eff} + 1/σ_{el,eff})^{1} (S m^{1}), σ_{el,eff}Effective electronic conductivity (S m^{1}) σ_{io,eff}Effective ionic conductivity (S m^{1})
Acknowledgements
This research was supported by the Daegu University Research Grant (Grant No. 20150130).
Larminie J.
,
Dicks A.
,
McDonald M. S.
2003
Fuel Cell Systems Explained
2nd Ed.
John Wiley & Sons
Chichester
O’Hayre R.
,
Cha S. W.
,
Colella W.
,
Prinz F. B.
2009
Fuel Cell Fundamentals
2nd Ed.
John Wiley & Sons
New York
Brett D. J.
,
Atkinson A.
,
Brandon N. P.
,
Skinner S. J.
2008
Chem. Soc. Rev.
37
(8)
1568 
1578
DOI : 10.1039/b612060c
Bard A. J.
,
Faulkner L. R.
2001
Electrochemical Methods: Fundamentals and Applications
John Wiley & Sons
New York
Utz A.
,
Störmer H.
,
Leonide A.
,
Weber A.
,
IversTiffée E.
2010
J. Electrochem. Soc.
157
(6)
B920 
B930
DOI : 10.1149/1.3383041
Holtappels P.
,
de Haart L. G. J.
,
Stimming U.
1999
J. Electrochem. Soc.
146
(5)
1620 
1625
DOI : 10.1149/1.1391816
Klein S. A.
,
Alvarado F. L.
2002
Engineering Equation Solver
FChart Software
Madison