Advanced
Electrochemical Effectiveness Factors for Butler-Volmer Reaction Kinetics in Active Electrode Layers of Solid Oxide Fuel Cells
Electrochemical Effectiveness Factors for Butler-Volmer Reaction Kinetics in Active Electrode Layers of Solid Oxide Fuel Cells
Journal of Electrochemical Science and Technology. 2017. Dec, 8(4): 344-355
Copyright © 2017, The Korean Electrochemical Society
  • Received : October 26, 2017
  • Accepted : December 06, 2017
  • Published : December 31, 2017
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
About the Authors
Jin Hyun Nam
jinhnam@gmail.com

Abstract
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 Butler-Volmer 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.3-0.8. These effectiveness data were then fitted to simple correlation equations, and the resulting correlation coefficients are presented along with estimated magnitude of error.
Keywords
1. Introduction
Solid oxide fuel cells (SOFCs) are a promising high-temperature fuel cell technology that can efficiently convert the chemical energy of fuel into electricity [1 , 2] . The high operating temperature (~800-1000℃) of SOFCs reduces manufacturing costs by employing non-precious 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 long-term performance stability. Thus, intermediate-temperature SOFCs (IT-SOFCs) that can operate at ~650-800℃ have emerged as a solution to the long-term performance degradation problem [3 - 5] .
To achieve better performance of IT-SOFCs, 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 IT-SOFCs. In addition, recent experimental studies have directly reconstructed three-dimensional 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 IT-SOFCs, 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 three-phase 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 IT-SOFC 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 current-local 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 Butler-Volmer 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 Butler-Volmer reaction kinetics. Their model was successfully used for one-dimensional simulation of a single-cell 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 Butler-Volmer 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 two-layer anode of IT-SOFCs 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 anode-supported IT-SOFCs). 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.
PPT Slide
Lager Image
Electrochemical reactions and transport processes (mass diffusion and charge conduction) in the two-layer anode of IT-SOFCs, along with the resultant electronic current, Iel, ionic current, Iio, and charge transfer current, Itr.
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 volume-specific TPB length (TPBL), λ tpb,V , and the effective ionic conductivity, σ io,eff , are uniform inside the active reaction layer.
• Uniform operating condition: The TPBL-specific 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 IT-SOFCs (shown in Fig. 1 ) as [20 , 21]
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
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
PPT Slide
Lager Image
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
PPT Slide
Lager Image
PPT Slide
Lager Image
Here, η conc is the total activation overpotential applied on the active reaction layer, determined as
PPT Slide
Lager Image
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 Butler-Volmer equation, the volumetric charge transfer current density, i tr,V ( η ), is expressed as
PPT Slide
Lager Image
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 TPBL-specific exchange current density, i tpb , and the volume-specific TPBL, λ tpb,V as i ex,V = i tpb × λ tpb,V . Note that Eq. (9) may also be expressed in terms of the TPBL-specific polarization resistance, r tpb (Ωm), as
PPT Slide
Lager Image
Fig. 2 shows the charge transfer current density vs. local activation overpotential ( i tr,V - η ) curves for the Butler-Volmer 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 Butler-Volmer 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
PPT Slide
Lager Image
Charge transfer current density vs. local activation overpotential (itr,V-η) curves for general Butler-Volmer reaction kinetics of Eq. (9) with various charge transfer coefficients. These curves are obtained for anodic reaction (ne = 2) at T = 1000 K.
PPT Slide
Lager Image
where r p,V is the volume-specific 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 Butler-Volmer 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.6-0.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 Butler-Volmer 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]
PPT Slide
Lager Image
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
PPT Slide
Lager Image
where ϕ T is the electrochemical Thiele modulus, defined as
PPT Slide
Lager Image
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 5-20 for the anode reaction layers and 0.5-2.0 for the cathode reaction layers, for ordinary operation of IT-SOFCs [21 , 22] .
Shin and Nam [21] studied the effectiveness factor for the symmetric Butler-Volmer 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,
PPT Slide
Lager Image
, as
PPT Slide
Lager Image
Here,
PPT Slide
Lager Image
is the dimensionless total activation over-potential applied on the active reaction layer
PPT Slide
Lager Image
. In Eq. (15),
PPT Slide
Lager Image
depends on
PPT Slide
Lager Image
, not separately on α , n e , T or η tot . It is also shown that
PPT Slide
Lager Image
has a nearly constant functional form with respect to
PPT Slide
Lager Image
for ϕ T ≥ 3.
It is not convenient to present all the effectiveness factors for various ϕ T and
PPT Slide
Lager Image
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
PPT Slide
Lager Image
where a , b , c , and d are the correlation coefficients dependent on ϕ T . Eq. (16) accurately describes the behavior of the relative effectiveness,
PPT Slide
Lager Image
, 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 TPBL-specific exchange current density, i tpb , or the TPBL-specific 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,
PPT Slide
Lager Image
, 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
PPT Slide
Lager Image
- 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 one-dimensional 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,
PPT Slide
Lager Image
. In this study, numerical calculations were performed for both the symmetric and asymmetric Butler-Volmer 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 Butler-Volmer reaction kinetics
Previously, Shin and Nam [21] investigated the electrochemical effectiveness factors for current generation in the active reaction layers, subject to the symmetric Butler-Volmer 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,
PPT Slide
Lager Image
, which was numerically determined by solving Eqs. (6), (7), and (10) with α = 0.5 (symmetric Butler-Volmer equation). In Fig. 3 ,
PPT Slide
Lager Image
generally decreases from 1.0 towards 0.0 as increases. For a given
PPT Slide
Lager Image
,
PPT Slide
Lager Image
becomes smaller as ϕ T increases. It is interesting to note that
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
Relative effectiveness factor, , for the current generation performance of the active reaction layer in IT-SOFCs: Symmetric Butler-Volmer 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
PPT Slide
Lager Image
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 Butler-Volmer reaction kinetics, described by Eqs. (6), (7), and (10).
Correlation coefficients to determinefor symmetric Butler-Volmer reaction kinetics withα= 0.5.
PPT Slide
Lager Image
aErr0.3V refers to the maximum correlation error in the estimated effectiveness factor for the cathode total overpotential range of 0.0-0.3 V at T=1000 K (or equally, for the anode total overpotential range of 0.0-0.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 current-overpotential performance curves obtained by the effectiveness model and those obtained by the detailed electrode microscale model. In addition, a one-dimensional simulation model was also developed to predict the current-voltage performance curves of a single-cell SOFC with two-layer 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 Butler-Volmer 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
PPT Slide
Lager Image
for α > 0.5 shown in Fig. 4 is similar to the behavior for α = 0.5 shown in Fig. 3 , except that the range of
PPT Slide
Lager Image
increases with α according to
PPT Slide
Lager Image
= α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 over-potential in the range of 0-0.3 V (equivalently, the anode total overpotential range of 0-0.6 V) at T = 1000 K.
PPT Slide
Lager Image
Relative effectiveness factor, , for the current generation performance of the active reaction layer in IT-SOFCs: Asymmetric Butler-Volmer reaction kinetics with (a) α = 0.6, (b) α = 0.7, and (c) α = 0.8.
Correlation coefficients to determinefor asymmetric Butler-Volmer reaction kinetics withα= 0.6.
PPT Slide
Lager Image
Correlation coefficients to determine for asymmetric Butler-Volmer reaction kinetics with α = 0.6.
Correlation coefficients to determinefor asymmetric Butler-Volmer reaction kinetics withα= 0.7
PPT Slide
Lager Image
Correlation coefficients to determine for asymmetric Butler-Volmer reaction kinetics with α = 0.7
Correlation coefficients to determinefor asymmetric Butler-Volmer reaction kinetics withα= 0.8.
PPT Slide
Lager Image
Correlation coefficients to determine for asymmetric Butler-Volmer 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 0-0.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
PPT Slide
Lager Image
relevant to asymmetric Butler-Volmer 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 Butler-Volmer 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 .
PPT Slide
Lager Image
Relative effectiveness factor, , for the current generation performance of the active reaction layer in IT-SOFCs: Asymmetric Butler-Volmer reaction kinetics with (a) α = 0.4 and (b) α = 0.3.
Correlation coefficients to determinefor asymmetric Butler-Volmer reaction kinetics withα= 0.4.
PPT Slide
Lager Image
Correlation coefficients to determine for asymmetric Butler-Volmer reaction kinetics with α = 0.4.
Correlation coefficients to determinefor asymmetric Butler-Volmer reaction kinetics withα= 0.3.
PPT Slide
Lager Image
Correlation coefficients to determine for asymmetric Butler-Volmer reaction kinetics with α = 0.3.
Fig. 5 shows that the relative effectiveness,
PPT Slide
Lager Image
, has values higher than 1.0 for small
PPT Slide
Lager Image
ranges, which is different from the behavior of
PPT Slide
Lager Image
for symmetric ( α = 0.5) and asymmetric ( α > 0.5) Butler-Volmer 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 Butler-Volmer 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,
PPT Slide
Lager Image
, 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
PPT Slide
Lager Image
which is smaller than 1.0 in Figs. 3 and 4 . In contrast, the Butler-Volmer 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 0-0.12 V with α = 0.4 and for η in the range of 0-0.26 V with α = 0.3 in Fig. 2 ). Thus, the relative effectiveness,
PPT Slide
Lager Image
, becomes larger than 1.0 in those low overpotential regions (for
PPT Slide
Lager Image
in the range of 0-1.1 with α = 0.4 and for
PPT Slide
Lager Image
in the range of 0-1.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
PPT Slide
Lager Image
greater than 1.0. The lines in Fig. 5 represent the fitting curves of Eq. (18).
PPT Slide
Lager Image
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,
PPT Slide
Lager Image
, correspond to the cathode total overpotential range of 0-0.05 V or 0-0.08 V at T = 1000 K. For
PPT Slide
Lager Image
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 Butler-Volmer 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 Butler-Volmer 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 ielElectronic current density (A m-2) iex,VVolume-specific exchange current density (A m-3) iioIonic current density (A m-2) imax,AMaximum current density in the active reaction layer (A m-2) ireal,AActual current density in the active reaction layer (A m-2) itpbTPBL-specific exchange current density (A m-1) itr,VVolumetric charge transfer current density (A m-3) LActive reaction layer thickness (m) neNumber of electrons for electrochemical reactions (2 for anode, 4 for cathode) RUniversal gas constant (8.314 J kg-1 K-1) rp,VVolume-specific linear polarization resistance (Ω m3) rtpbTPBL-specific linear polarization resistance (Ω m) TTemperature (K) zCoordinate (m) αCharge transfer coefficient ϕTThiele modulus ΓeffElectrochemical effectiveness factor Γeff,0VBase effectiveness factor at zero overpotential ηActivation overpotential (V) ηconcConcentration overpotential (V) ηtotTotal activation overpotential applied to the active reaction layer (V) Dimensionless total activation overpotential, αneFηtot ⁄ / (RT) φoNernst potential (V) φelElectronic potential (V) φioIonic potential (V) λtpb,VVolume-specific TPBL (m m-3) σeffEffective charge conductivity, σeff = (1/σio,eff + 1/σel,eff)-1 (S m-1), σel,effEffective electronic conductivity (S m-1) σio,effEffective ionic conductivity (S m-1)
Acknowledgements
This research was supported by the Daegu University Research Grant (Grant No. 20150130).
References
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
Huijsmans J. P. P. , Van Berkel F. P. F. , Christie G. M. 1998 J. Power Sources 71 (1) 107 - 110    DOI : 10.1016/S0378-7753(97)02789-4
Brett D. J. , Atkinson A. , Brandon N. P. , Skinner S. J. 2008 Chem. Soc. Rev. 37 (8) 1568 - 1578    DOI : 10.1039/b612060c
Wachsman E. D. , Lee K. T. 2011 Science 334 935 - 939    DOI : 10.1126/science.1204090
Maguire E. , Gharbage B. , Marques F. M. B. , Labrincha J. A. 2000 Solid State Ionics 127 (3) 329 - 335    DOI : 10.1016/S0167-2738(99)00286-6
Ishihara T. , Yan J. , Shinagawa M. , Matsumoto H. 2006 Electrochim. Acta 52 (4) 1645 - 1650    DOI : 10.1016/j.electacta.2006.03.103
Fu C. , Sun K. , Zhang N. , Chen X. , Zhou D. 2007 Electrochim. Acta 52 (13) 4589 - 4594    DOI : 10.1016/j.electacta.2007.01.001
Holtappels P. , Bagger C. 2002 J. Eur. Ceram. Soc. 22 (1) 41 - 48    DOI : 10.1016/S0955-2219(01)00238-2
Zhao F. , Virkar A. V. 2005 J. Power Sources 141 (1) 79 - 95    DOI : 10.1016/j.jpowsour.2004.08.057
Haanappel V. A. C. , Mertens J. , Rutenbeck D. , Tropartz C. , Herzhof W. , Sebold D. , Tietz F. 2005 J. Power Sources 141 (2) 216 - 226    DOI : 10.1016/j.jpowsour.2004.09.016
Jeon D. H. , Nam J. H. , Kim C. J. 2006 J. Electrochem. Soc. 153 (2) A406 - A417    DOI : 10.1149/1.2139954
Ni M. , Leung M. K. , Leung D. Y. 2007 J. Power Sources 168 (2) 369 - 378    DOI : 10.1016/j.jpowsour.2007.03.005
Wang Z. , Zhang N. , Qiao J. , Sun K. , Xu P. 2009 Electrochem. Commun. 11 (6) 1120 - 1123    DOI : 10.1016/j.elecom.2009.03.027
Chen Y. , Bunch J. , Li T. , Mao Z. , Chen F. 2012 J. Power Sources 213 93 - 99    DOI : 10.1016/j.jpowsour.2012.03.109
Iwai H. 2010 J. Power Sources 195 (4) 955 - 961    DOI : 10.1016/j.jpowsour.2009.09.005
Wilson J. R. , Cronin J. S. , Barnett S. A. 2011 Scripta Mater. 65 (2) 67 - 72    DOI : 10.1016/j.scriptamat.2010.09.025
Cronin J. S. , Chen-Wiegart Y. K. , Wang J. , Barnett S. A. 2013 J. Power Sources 233 174 - 179    DOI : 10.1016/j.jpowsour.2013.01.060
Costamagna P. , Costa P. , Arato E. 1998 Electrochim. Acta 43 (3) 375 - 394    DOI : 10.1016/S0013-4686(97)00063-7
Costamagna P. , Costa P. , Arato E. 1998 Electrochim. Acta 43 (8) 967 - 972    DOI : 10.1016/S0013-4686(97)00262-4
Shin D. , Nam J. H. 2015 Electrochim. Acta 171 1 - 6    DOI : 10.1016/j.electacta.2015.04.171
Shin D. , Baek S. M. , Nam J. H. , Kim C. J. 2016 Comp. Chem. Eng. 90 268 - 277    DOI : 10.1016/j.compchemeng.2016.04.032
Baek S. M. , Shin D. , Sohn S. , Nam J. H. 2016 Fuel Cells 16 (5) 591 - 599    DOI : 10.1002/fuce.201500204
Nam J. H. 2016 Electrochim. Acta 221 8 - 13    DOI : 10.1016/j.electacta.2016.10.145
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. , Ivers-Tiffée E. 2010 J. Electrochem. Soc. 157 (6) B920 - B930    DOI : 10.1149/1.3383041
Mizusaki J. 1994 Solid State Ionics 70-71 52 - 58    DOI : 10.1016/0167-2738(94)90286-0
de Boer B. , Gonzalez M. , Bouwmeester H. J. M. , Verweij H. 2000 Solid State Ionics 127 (3) 269 - 276    DOI : 10.1016/S0167-2738(99)00299-4
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 F-Chart Software Madison