Next Article in Journal
Homochiral Metal-Organic Framework Based Mixed Matrix Membrane for Chiral Resolution
Previous Article in Journal
Block Copolymer Membranes—Progress and Challenges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analytical Model for Concentration (Pressure) Impedance of a Low-Pt PEM Fuel Cell Oxygen Electrode

by
Andrei Kulikovsky
Theory and Computation of Energy Materials (IEK-13), Institute of Energy and Climate Research, Forschungszentrum Jülich GmbH, D-52425 Jülich, Germany
Membranes 2022, 12(4), 356; https://doi.org/10.3390/membranes12040356
Submission received: 28 February 2022 / Revised: 18 March 2022 / Accepted: 22 March 2022 / Published: 24 March 2022

Abstract

:
In this study, a model for concentration/pressure impedance ζ of the cathode catalyst layer of a low-Pt PEM fuel cell is developed. The model is based on transient oxygen mass transport equations through the cathode catalyst layer modeled as a single pore with a thin Nafion film covering the pore surface. This structure is used to simulate oxygen transport through the catalyst layer depth and through the ionomer film covering Pt/C agglomerates in low-Pt cells. Analytical solution for zeta-impedance at high cell current is derived; this solution can be used for fast fitting of experimental zeta-spectra. Optimal conditions for measuring the zeta-spectra of a low-Pt cell are discussed. Zeta impedance is not affected by faradaic processes in the cell, which makes this technique a useful alternative to standard EIS.

1. Introduction

Lowering of Pt loading is of large importance for the success of PEM fuel cells on the mass market. However, the first attempts to lower Pt loading on the cathode side have already revealed unexpected “overlinear” transport loss of the cell potential [1,2]. This effect has been attributed to oxygen transport through a thin ionomer film covering Pt/C agglomerates in the cathode catalyst layer (CCL). Over the past decade, a lot of research has been done to measure the film transport resistivity R N in PEMFCs [2,3,4,5,6,7,8,9]. Most of the works employed a limiting current method for measuring R N , moreover, in [9] hydrogen instead of oxygen was used in the limiting current experiments to avoid undesired effects due to oxygen reduction reaction (ORR) in the electrode.
An alternative method, requiring standard cell, equipment, and procedures, is electrochemical impedance spectroscopy (EIS). A numerical model for low-Pt PEMFC impedance has been developed and fitted to the experimental impedance spectra [10]. Fitting returned the Nafion film thickness and oxygen diffusivity. Nonetheless, complexity of the impedance model [10] stimulates searching for simpler ways for measuring R N . Generally, every transport process in a fuel cell is expected to form a separate peak in the distribution of relaxation times (DRT) spectrum [11]. The area under the DRT peak gives the contribution of the respective process to the total cell polarization resistivity. Thus, calculation of R N from experimental impedance spectra would be much simpler if the DRT peak due to oxygen transport in Nafion film were well separated from other peaks. However, a recent model [12] shows that at low currents, the characteristic frequency of the film-transport peak is only 1.73 times less than the frequency of faradaic peak and the two peaks tend to merge. An impedance spectroscopy method insensitive to faradaic processes in the cell would be of great help.
The idea to perturb cell potential by harmonic variation of the oxygen concentration or pressure (EPIS) has been suggested in [13] and developed further in [14,15,16,17,18]. A review of pressure impedance works can be found in [16]. The great advantage of this technique is that the zeta-spectrum of the cell is independent of the faradaic processes. So far, most of the reported EPIS models have been numerical. An analytical model for the PEMFC zeta-impedance taking into account oxygen transport in the gas diffusion layer (GDL) and CCL has been derived in [19]; however, this mean-field model ignores oxygen transport in ionomer film covering Pt/C agglomerates.
In this work, we report a model for the concentration impedance ζ of a low-Pt cell, with explicit account of oxygen transport through the ionomer film. We derive a formula for the ζ -impedance of a low-Pt cell suitable for fast fitting of experimental zeta-spectra. We show that measurements of the oxygen diffusion coefficient in the Nafion film using zeta-spectroscopy should be done at low oxygen concentration providing high rate of oxygen consumption in the CCL and at the cell current density close to the limiting current density due to Nafion film.

2. Model

2.1. Basic Transient Equations

It is assumed that the following conditions hold
  • Oxygen transport loss in the gas-diffusion layer is small;
  • Proton transport in the CCL is fast;
  • Oxygen consumption in the CCL is large.
The limiting current density due to oxygen transport in the GDL is typically much larger, than the limiting current due to Nafion film [20] and to a first approximation oxygen transport loss in the GDL can be neglected. The second and third assumptions are discussed in Section 3.
The low–Pt CCL is modeled by a single cylindrical pore penetrating through the whole CCL depth. The pore volume is separated from the coaxial Pt/C tube by a thin Nafion film (Figure 1). Oxygen is transported along the pore and in the radial direction through the ionomer to Pt surface, where the ORR occurs. The static version of this model has been reported in [20]; here we briefly repeat the basic equations necessary for understanding the impedance model.
Oxygen transport along the pore is described by
c t D p 2 c x 2 = 2 N N , p R p , c x x = 0 = 0 , c ( l t ) = c 1
where R p is the pore radius, D p is the oxygen diffusion coefficient in the pore, c is the oxygen concentration in the pore, c 1 is the oxygen concentration at the pore/GDL interface and
N N , p = D N c N r r = R p
is the radial oxygen flux in the Nafion film at the pore/film interface.
Radial oxygen transport through the film is described by the diffusion equation
c N t D N r r r c N r = 0 , c N ( R p ) = K H c ( x ) , D N c N r r = R m = R p i * 2 ( 4 F ) c N , m c r e f exp η b
where c N is the oxygen concentration in the Nafion film, c N , m c N ( R m ) , D N is the oxygen diffusion coefficient in the film, i * is the ORR exchange current density, η is the positive by convention ORR overpotential, and b is the ORR Tafel slope. The left boundary condition for Equation (3) is Henry’s law for oxygen dissolution in Nafion. The right boundary condition to Equation (3) describes consumption of dissolved oxygen in the ORR; the factor R p / 2 provides correct transition to the standard oxygen mass conservation equation in the CCL for the limiting case of zero Nafion film thickness.
To simplify the calculations we introduce dimensionless variables
x ˜ = x l t , R ˜ = r l t , t ˜ = t i * 4 F c r e f , η ˜ = η b , j ˜ = j i * l t , D ˜ = 4 F D c r e f i * l t 2 , ω ˜ = ω 4 F c r e f i * , ζ ˜ = ζ c r e f b
where l t is the pore length (CCL thickness), ω is the angular frequency of AC signal, and ζ is the concentration impedance (see below).
With Equation (4), Equations (1) and (3) transform to
c ˜ t ˜ D ˜ p 2 c ˜ x ˜ 2 = 2 N ˜ N , p R ˜ p
c ˜ N t ˜ D ˜ N r ˜ r ˜ r ˜ c ˜ N r ˜ = 0 , c ˜ N ( R ˜ p ) = K H c ˜ ( x ˜ ) , D ˜ N c ˜ N r ˜ r ˜ = R ˜ m = R ˜ p 2 c ˜ N , m exp η ˜
Systems (5) and (6) form the basis for the concentration impedance model.

2.2. Equations for Perturbation Amplitudes

Now we apply small-amplitude perturbations of the form
c ˜ ( x ˜ , t ˜ ) = c ˜ 0 ( x ˜ ) + c ˜ 1 ( x ˜ , ω ˜ ) exp ( i ω ˜ t ˜ ) c ˜ N ( x ˜ , t ˜ ) = c ˜ N 0 ( x ˜ ) + c ˜ N 1 ( x ˜ , ω ˜ ) exp ( i ω ˜ t ˜ ) η ˜ ( t ˜ ) = η ˜ 0 + η ˜ 1 ( ω ˜ ) exp ( i ω ˜ t ˜ ) .
where the superscripts 0 and 1 mark the static functions and the perturbation amplitudes, respectively. Note that fast proton transport means that the static and perturbed ORR overpotentials are independent of x ˜ .
Substituting Equation (7) into Equations (5) and (6) and performing standard procedure of linearization, we come to the system of linear equations for c ˜ 1 and c ˜ N 1 :
D ˜ p 2 c ˜ 1 x ˜ 2 = 2 N ˜ N , p 1 R ˜ p + i ω ˜ c ˜ 1 , c ˜ 1 x ˜ x ˜ = 0 = 0 , c ˜ 1 ( 1 ) = c ˜ 1 1
D ˜ N r ˜ r ˜ r ˜ c ˜ N 1 r ˜ = i ω ˜ c ˜ N 1 , c ˜ N 1 ( R ˜ p ) = K H c ˜ 1 ( x ˜ ) , D ˜ N c ˜ N 1 r ˜ r ˜ = R ˜ m = R ˜ p 2 e η ˜ 0 c ˜ N 1 ( R ˜ m ) + c ˜ N , m 0 η ˜ 1 ,
where c ˜ N , m 0 = c ˜ N 0 ( R ˜ m ) , c ˜ 1 1 is the oxygen perturbation at the pore/GDL interface, and
N ˜ N , p 1 = D ˜ N c ˜ N 1 r ˜ r ˜ = R ˜ p
is the perturbed oxygen flux in the Nafion film at the pore interface. Note that due to assumption of fast oxygen transport in the GDL we have
c ˜ 1 1 = c ˜ h 1 ,
i.e., c ˜ 1 1 is equal to the applied oxygen concentration perturbation in the channel c ˜ h 1 .
The system of Equations (8) and (9) with the coefficient functions given by Equations (13) and (14) in the next section determine the CCL concentration impedance ζ ˜
ζ ˜ = η ˜ 1 c ˜ h 1 .
Here, c ˜ h 1 is the oxygen concentration perturbation applied in the gas channel, and η ˜ 1 is the measured perturbation of the cell potential.

2.3. Static Equations and Solutions

Static equations for c ˜ 0 and c ˜ N 0 are obtained from Equations (5) and (6) simply by chalking out the time derivatives:
D ˜ p 2 c ˜ 0 x ˜ 2 = 2 R ˜ p D ˜ N c ˜ N 0 r ˜ r ˜ = R ˜ p , c ˜ 0 x ˜ x ˜ = 0 = 0 , c ˜ 0 ( 1 ) = c ˜ 1
D ˜ N r ˜ r ˜ r ˜ c ˜ N 0 r ˜ = 0 , c ˜ N 0 ( R ˜ p ) = K H c ˜ 0 ( x ˜ ) , D ˜ N c ˜ N 0 r ˜ r ˜ = R ˜ m = R ˜ p 2 c ˜ N , m 0 exp η ˜ 0
The solution to Equation (14) is
c ˜ N 0 ( r ˜ ) = R ˜ p R ˜ m ln R ˜ m / r ˜ e η ˜ 0 + 2 D ˜ N R ˜ p R ˜ m ln R ˜ m / R ˜ p e η ˜ 0 + 2 D ˜ N K H c ˜ 0 ( x ˜ )
Equation (15) allows us to calculate the flux D ˜ N c ˜ N 0 / r ˜ | r ˜ = R ˜ p :
D ˜ N c ˜ N 0 r ˜ r ˜ = R ˜ p = R ˜ m D ˜ N e η ˜ 0 K H c ˜ 0 ( x ˜ ) R ˜ p R ˜ m ln R ˜ m / R ˜ p e η ˜ 0 + 2 D ˜ N
With Equation (16), Equation (13) takes the form
D ˜ p 2 c ˜ 0 x ˜ 2 = ρ c ˜ 0 , c ˜ 0 x ˜ x ˜ = 0 = 0 , c ˜ 0 ( 1 ) = c ˜ 1 0
where
ρ = 2 R ˜ m D ˜ N e η ˜ 0 K H R ˜ p R ˜ p R ˜ m ln R ˜ m / R ˜ p e η ˜ 0 + 2 D ˜ N
The solution of Equation (17) reads
c ˜ 0 ( x ˜ ) = c ˜ 1 0 cosh x ˜ ρ / D ˜ p cosh ρ / D ˜ p .
Using Equation (19) in Equation (15) and setting r ˜ = R ˜ m , we get c ˜ N , m 0 , which appears in Equation (9):
c ˜ N , m 0 = ρ c ˜ 1 0 cosh x ˜ ρ / D ˜ p e η ˜ 0 cosh ρ / D ˜ p .
For further calculations we need the polarization curve of the CCL. The steady-state proton current conservation equation in the CCL is
j ˜ 0 x ˜ = c ˜ N , m 0 exp η ˜ 0
where j ˜ 0 is the local proton current density in the film. Integrating this equation over x ˜ from 0 to 1 with Equation (20) and the boundary conditions j ˜ 0 ( 0 ) = j ˜ 0 , j ˜ 0 ( 1 ) = 0 , we find
j ˜ 0 = c ˜ 1 0 ρ tanh ρ / D ˜ p ρ / D ˜ p .
With ρ given by Equation (18), this equation provides the explicit dependence of j ˜ 0 vs. the ORR overpotential η ˜ 0 , i.e., the static polarization curve of the system.
In two limiting cases, Equation (22) can be simplified. If the argument of the tanh-function is small, we may approximate tanh x x and Equation (22) simplifies to
j ˜ 0 = c ˜ 1 0 ρ .
With ρ , Equation (18), we get expression of η ˜ 0 through j ˜ 0 :
e η ˜ 0 = 2 D ˜ N R ˜ p j ˜ 0 R ˜ m 2 D ˜ N K H c ˜ 1 0 R ˜ p 2 ln R ˜ m / R ˜ p j ˜ 0 , when ρ / D ˜ p 1
Equating the zero denominator of Equation (24) we get a limiting current density due to the oxygen transport in the Nafion film [20]
j ˜ N lim = 2 D ˜ N K H c ˜ 1 0 R ˜ p 2 ln R ˜ m / R ˜ p .
Equations (24) and (25) do not contain D ˜ p , meaning that this case corresponds to fast oxygen transport in the void pore.
In the limit of ρ / D ˜ p 2 , we may replace the tanh-function in Equation (22) by unity, which leads to
j ˜ 0 = c ˜ 1 0 ρ ρ / D ˜ p .
With Equation (18) we find
e η ˜ 0 = 2 D ˜ N R ˜ p j ˜ 0 2 R ˜ m 2 D ˜ p D ˜ N K H ( c ˜ 1 0 ) 2 R ˜ p 2 ln R ˜ m / R ˜ p j ˜ 0 2 , when ρ / D ˜ p 2
Equation (27) determines the limiting current density
j ˜ N lim = 2 D ˜ p D ˜ N K H ( c ˜ 1 0 ) 2 R ˜ p 2 ln R ˜ m / R ˜ p
Equation (27) correlates with the assumption of large oxygen transport loss in the CCL and below, this equation will be used in numerical calculations.

3. Results and Discussion

Static solutions of the previous section allow us to calculate the concentration impedance. The solution to Equation (9) is a rather cumbersome expression containing Bessel functions. This solution is only needed for calculation of the perturbed oxygen flux N ˜ N , p 1 appearing in Equation (8). Rather tedious algebra leads to
N ˜ N , p 1 = D ˜ N Q P c K H c ˜ 1 ( x ˜ ) + P η c ˜ N , m 0 ( x ˜ ) e η ˜ 0 η ˜ 1
where the independent of x ˜ coefficients P c , P η and Q are given in Appendix A. Equation (29) with (20) allow us to solve the key Equation (8):
c ˜ 1 ( x ˜ ) = 2 ρ P η D ˜ N c ˜ 1 0 η ˜ 1 R ˜ p Q cos ( ϕ ) ( D ˜ p ϕ 2 + ρ ) cosh ρ / D ˜ p × cos ( ϕ x ˜ ) cosh ρ D ˜ p cos ( ϕ ) cosh x ˜ ρ D ˜ p + c ˜ h 1 cos ( ϕ x ˜ ) cos ϕ
where
ϕ = 2 P c D ˜ N K H R ˜ p D ˜ p Q i ω ˜ D ˜ p .
As discussed above, the model is valid in the limit of strong oxygen transport loss through the CCL depth. In this regime, the perturbation of oxygen concentration at the membrane surface is nearly zero: c ˜ 1 ( 0 ) 0 . Setting in Equation (30) x ˜ = 0 , c ˜ 1 ( 0 ) = 0 and dividing the resulting equation by c ˜ h 1 , we get an algebraic equation for zeta-impedance
0 = 2 ρ P η D ˜ N c ˜ 1 0 ζ ˜ cosh ρ / D ˜ p cos ϕ R ˜ p Q D ˜ p ϕ 2 + ρ cosh ρ / D ˜ p + 1
Solving Equation (32) for ζ ˜ , we finally find
ζ ˜ = R ˜ p Q D ˜ p ϕ 2 + ρ cosh ρ / D ˜ p 2 ρ P η D ˜ N c ˜ 1 0 cos ϕ cosh ρ / D ˜ p .
Equation (33) is the main result of this work.
The spectra of Equation (33) in the dimension form ζ = ζ ˜ b / c r e f for the two oxygen diffusion coefficients D N in the Nafion film are shown in Figure 2. The base-case set of parameters used in the calculations is collected in Table 1. As can be seen, for these parameters, the zeta-spectrum is quite sensitive to the value of D N (Figure 2a), which makes zeta-spectroscopy a good candidate for measuring film transport properties. Note that the variation of D N does not change the characteristic frequency of the curve in Figure 2b.
Variation of the zeta-spectrum with pore diffusivity D p is illustrated in Figure 3: the growth of D p increases the static “resistivity” of the spectrum and shifts the peak of − Im ζ to higher frequency (Figure 3b). The dependence of − Im ζ peak frequency on D p is close to linear: twice higher D p shifts the peak to twice larger frequency (Figure 3b). The characteristic frequency of peaks in Figure 3b is between 1 and 10 kHz, far above the characteristic frequency of faradaic processes in the cell, which typically does not exceed 100 Hz. It is interesting to note that variation of D N and D p lead to the opposite trends in the spectrum diameter: the curl increases with the decrease in D N and with the growth of D p (cf. Figure 2a and Figure 3a).
Of particular interest is the static value ζ ˜ 0 , corresponding to the rightmost point of the spectra in Figure 2 and Figure 3. Unfortunately, Maple ® fails to calculate ζ ˜ 0 = lim ω ˜ 0 ζ ˜ . Numerically calculated ζ 0 as a function of the Nafion film diffusivity is shown in Figure 4. For D N below 1.0 × 10 6 cm 2 s 1 the curve is very steep due to effect of limiting current: with the decrease in D N , the right side of the dimension version of Equation (28) tends to the current of 1.5 A cm 2 fixed in the calculations and the zeta-impedance rapidly increases. From Figure 4 it follows that measurements of D N should be best done close to the limiting current j N lim due to oxygen transport in the film. In this range of currents, the zeta-impedance is most sensitive to Nafion film transport parameters.
The model above is developed assuming fast proton transport in the CCL. This means that the cell current density must be much less than the characteristic current for proton transport in the CCL:
j 0 j * = σ N b l t
where σ N is the CCL proton conductivity. In a working PEMFC, σ N 0.02 S cm 1 , hence with b = 0.03 V and l t = 3 × 10 4 cm we get j * = 2 A cm 2 .
On the other hand, the condition c ˜ 1 1 0 holds for the cell current density satisfying to
j 0 4 F D p c 1 0 l t
while the largest sensitivity of the zeta-spectrum to the Nafion film oxygen diffusivity is achieved for the cell currents satisfying to
j 0 4 F c 1 0 2 D p D N K H R p 2 ln ( R m / R p )
which is the dimension version of Equation (28). Equations (35) and (36) allow one to select optimal c 1 0 and j 0 for experiments.

4. Conclusions

In this study, a model for concentration (zeta-) impedance of the low-Pt cathode catalyst layer in a PEM fuel cell is developed. The model is based on the transient oxygen mass transport equations in the cylindrical pore surrounded by a thin Nafion film separating the pore volume from Pt/C surface. An analytical solution for zeta-impedance is obtained, which can be used for fast fitting experimental zeta-spectra. It is shown that the zeta-spectrum is very sensitive to the value of oxygen diffusion coefficient in the Nafion film, provided that the cell operates at a high rate of oxygen consumption in the CCL and close to the limiting current density due to oxygen transport in the film.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

˜ Marks dimensionless variables
bORR Tafel slope, V
C d l Double layer volumetric capacitance, F cm 3
cOxygen molar concentration in the pore, mol cm 3
c 1 Static oxygen concentration at
the CCL/GDL interface, mol cm 3
c h Oxygen concentration in the channel, mol cm 3
c N Oxygen concentration in the Nafion film, mol cm 3
c N , m = c N ( R m )
c r e f Reference (inlet) oxygen concentration, mol cm 3
D b Oxygen diffusion coefficient in the GDL, cm 2  s 1
D p Oxygen diffusion coefficient in the pore, cm 2  s 1
D N Oxygen diffusion coefficient in the Nafion film, cm 2  s 1
FFaraday constant, C mol 1
i * ORR volumetric exchange current density, A cm 3
i Imaginary unit
jLocal proton current density along the pore, A cm 2
j N lim Limiting current density due to oxygen transport in Nafion film, A cm 2
j * Characteristic current density of proton transport, A cm 2 , Equation (34)
j 0 Cell current density, A cm 2
K H Dimensionless Henry’s constant for oxygen solubility in Nafion, mol/mol
l b GDL thickness, cm
l t Pore length (CCL thickness), cm
l N Nafion film thickness, cm
N N Radial oxygen flux in the Nafion film, mol cm 2 s 1
qAuxiliary parameter, Equation (A4)
R m Radius of a Pt/C tube, cm
R p Pore radius, cm
rRadial coordinate, cm
xCoordinate along the pore, cm
Subscripts:
0Membrane/CCL interface
1CCL/GDL interface
hChannel
mPt/C (metal) surface
NNafion film
pPore/Nafion film interface
Superscripts:
0Steady-state value
1Small-amplitude perturbation
limLimiting
Greek:
γ                             = 2 / R ˜ p
η ORR overpotential, positive by convention, V
ζ Concentration impedance, V cm 3  mol 1
ρ Auxiliary dimensionless parameter, Equation (18)
σ N Nafion film proton conductivity, S cm 1
ϕ Auxiliary dimensionless parameter, Equation (31)
ω ˜ = ω 4 F c r e f / i * , dimensionless frequency
ω Angular frequency of the AC signal, s 1

Appendix A. Coefficients in Equations (29)–(33)

P c = γ q 2 D ˜ N ( K 1 ( q R ˜ m ) I 1 ( q R ˜ p ) K 1 ( q R ˜ p ) I 1 ( q R ˜ m ) ) + q e η ˜ 0 ( K 0 ( q R ˜ m ) I 1 ( q R ˜ p ) K 1 ( q R ˜ p ) I 0 ( q R ˜ m ) )
P η = q ( K 0 ( q R ˜ p ) I 1 ( q R ˜ p ) K 1 ( q R ˜ p ) I 0 ( q R ˜ p ) )
Q = γ q D ˜ N ( K 1 ( q R ˜ m ) I 0 ( q R ˜ p ) K 0 ( q R ˜ p ) I 1 ( q R ˜ m ) ) + e η ˜ 0 ( K 0 ( q R ˜ m ) I 0 ( q R ˜ p ) K 0 ( q R ˜ p ) I 0 ( q R ˜ m ) )
q = i ω ˜ D ˜ N , γ = 2 R ˜ p
Here, I 0 , I 1 are the Bessel functions of the second kind, K 0 , K 1 are the modified Bessel functions of the second kind.

References

  1. Ono, Y.; Mashio, T.; Takaichi, S.; Ohma, A.; Kanesaka, H.; Shinohara, K. The Analysis of Performance Loss with Low Platinum Loaded Cathode Catalyst Layers. ECS Trans. 2010, 28, 69–78. [Google Scholar] [CrossRef]
  2. Greszler, T.A.; Caulk, D.; Sinha, P. The Impact of Platinum Loading on Oxygen Transport Resistance. J. Electrochem. Soc. 2012, 159, F831–F840. [Google Scholar] [CrossRef]
  3. Nonoyama, N.; Okazaki, S.; Weber, A.Z.; Ikogi, Y.; Yoshida, T. Analysis of Oxygen-Transport Diffusion Resistance in Proton-Exchange-Membrane Fuel Cells. J. Electrochem. Soc. 2011, 158, B416–B423. [Google Scholar] [CrossRef]
  4. Owejan, J.P.; Owejan, J.E.; Gu, W. Impact of Platinum Loading and Catalyst Layer Structure on PEMFC Performance. J. Electrochem. Soc. 2013, 160, F824–F833. [Google Scholar] [CrossRef]
  5. Weber, A.Z.; Kusoglu, A. Unexplained Transport Resistances for Low–Loaded Fuel–Cell Catalyst Layers. J. Mater. Chem. A 2014, 2, 17207–17211. [Google Scholar] [CrossRef]
  6. Kongkanand, A.; Mathias, M.F. The Priority and Challenge of High-Power Performance of LowPlatinum Proton-Exchange Membrane Fuel Cells. Phys. Chem. Lett. 2016, 7, 1127–1137. [Google Scholar] [CrossRef]
  7. Kudo, K.; Jinnouchi, R.; Morimoto, Y. Humidity and Temperature Dependences of Oxygen Transport Resistance of Nafion Thin Film on Platinum Electrode. Electrochim. Acta 2016, 209, 682–690. [Google Scholar] [CrossRef]
  8. Freiberg, A.T.S.; Tucker, M.C.; Weber, A.Z. Polarization Loss Correction Derived from Hydrogen Local–Resistance Measurement in Low Pt–Loaded Polymer–Electrolyte Fuel Cells. Electrochem. Commun. 2017, 79, 14–17. [Google Scholar] [CrossRef] [Green Version]
  9. Schuler, T.; Chowdhury, A.; Freiberg, A.T.; Sneed, B.; Spingler, F.B.; Tucker, M.C.; More, K.L.; Radke, C.J.; Weber, A.Z. Fuel-Cell Catalyst-Layer Resistance via Hydrogen Limiting-Current Measurementss. J. Electrochem. Soc. 2019, 166, F3020–F3031. [Google Scholar] [CrossRef] [Green Version]
  10. Reshetenko, T.; Kulikovsky, A. A Single–Pore Model for Cathode Catalyst Layer Impedance: The Effect of Nafion Film on PEM Fuel Cell Performance. RSC Adv. 2019, 9, 38797–38806. [Google Scholar] [CrossRef] [Green Version]
  11. Ivers-Tiffée, E.; Weber, A. Evaluation of Electrochemical Impedance Spectra by the Distribution of Relaxation Times. J. Ceram. Soc. Jpn. 2017, 125, 193–201. [Google Scholar] [CrossRef] [Green Version]
  12. Kulikovsky, A. Impedance and resistivity of low–Pt cathode in a PEM fuel cell. J. Electrochem. Soc. 2021, 168, 044512. [Google Scholar] [CrossRef]
  13. Niroumand, A.M.; Merida, W.; Eikerling, M.; Saif, M. Pressure–Voltage Oscillations as a Diagnostic Tool for PEFC Cathodes. Electrochem. Commun. 2010, 12, 122–124. [Google Scholar] [CrossRef]
  14. Engebretsen, E.; Mason, T.J.; Shearing, P.R.; Hinds, G.; Brett, D.J.L. Electrochemical Pressure Impedance Spectroscopy Applied to the Study of Polymer Electrolyte Fuel Cells. Electrochem. Commun. 2017, 75, 60–63. [Google Scholar] [CrossRef]
  15. Sorrentino, A.; Vidakovic-Koch, T.; Hanke-Rauschenbach, R.; Sundmacher, K. Concentration–Alternating Frequency Response: A New Method for Studying Polymer Electrolyte Membrane Fuel Cell Dynamics. Electrochim. Acta 2017, 243, 53–64. [Google Scholar] [CrossRef] [Green Version]
  16. Shirsath, A.V.; Rael, S.; Bonnet, C.; Schiffer, L.; Bessler, W.; Lapicque, F. Electrochemical Pressure Impedance Spectroscopy for Investigation of Mass Transfer in Polymer Electrolyte Membrane Fuel Cells. Curr. Opin. Electrochem. 2020, 20, 82–87. [Google Scholar] [CrossRef]
  17. Kubannek, F.; Krewer, U. Studying the Interaction of Mass Transport and Electrochemical Reaction Kinetics by Species Frequency Response Analysis. J. Electrochem. Soc. 2020, 167, 144510. [Google Scholar] [CrossRef]
  18. Schiffer, L.; Shirsath, A.V.; Raël, S.; Lapicque, F.; Bessler, W.G. Electrochemical Pressure 190 Impedance Spectroscopy for Polymer Electrolyte Membrane Fuel Cells: A Combined Modeling 191 and Experimental Analysis. J. Electrochem. Soc. 2021, 169, 034503. [Google Scholar] [CrossRef]
  19. Kulikovsky, A. Analytical Model for PEM Fuel Cell Concentration Impedance. J. Electroanal. Chem. 2021, 899, 115672. [Google Scholar] [CrossRef]
  20. Kulikovsky, A. The Effect of Nafion Film on the Cathode Catalyst Layer Performance in a Low–Pt PEM Fuel Cell. Electrochem. Commun. 2019, 103, 61–65. [Google Scholar] [CrossRef]
Figure 1. Schematic of a single-pore model for impedance of the low-Pt cathode catalyst layer.
Figure 1. Schematic of a single-pore model for impedance of the low-Pt cathode catalyst layer.
Membranes 12 00356 g001
Figure 2. (a) The Nyquist spectra and (b) the frequency dependence of imaginary part of ζ -impedance, Equation (33), for the indicated values of oxygen diffusion coefficient in the Nafion film D N , cm 2 s 1 . The other parameters for calculations are listed in Table 1.
Figure 2. (a) The Nyquist spectra and (b) the frequency dependence of imaginary part of ζ -impedance, Equation (33), for the indicated values of oxygen diffusion coefficient in the Nafion film D N , cm 2 s 1 . The other parameters for calculations are listed in Table 1.
Membranes 12 00356 g002
Figure 3. (a) The Nyquist spectra and (b) the frequency dependence of imaginary part of ζ -impedance, Equation (33), for the indicated values of oxygen diffusion coefficient in the pore D p , cm 2 s 1 . The other parameters for calculations are listed in Table 1.
Figure 3. (a) The Nyquist spectra and (b) the frequency dependence of imaginary part of ζ -impedance, Equation (33), for the indicated values of oxygen diffusion coefficient in the pore D p , cm 2 s 1 . The other parameters for calculations are listed in Table 1.
Membranes 12 00356 g003
Figure 4. The dependence of static concentration impedance ζ 0 on the oxygen diffusion coefficient in the Nafion film D N . The spectra in Figure 2 are plotted at the points indicated by red dots.
Figure 4. The dependence of static concentration impedance ζ 0 on the oxygen diffusion coefficient in the Nafion film D N . The spectra in Figure 2 are plotted at the points indicated by red dots.
Membranes 12 00356 g004
Table 1. The base-case cell parameters used in calculations.
Table 1. The base-case cell parameters used in calculations.
Tafel slope b, V0.03
Exchange current density i * , A cm 3 10 3
Double layer capacitance C d l , F cm 3 20
Oxygen diffusion coefficient in
the Nafion film, D N , cm 2  s 1 1 × 10 6
Dimensionless Henry’s constant for
O 2 solubility in water at 80 C, K H 6.76 × 10 3
Catalyst layer thickness l t , cm 3 × 10 4 (3 μ m)
Nafion film thickness l N , cm 10 × 10 7 (10 nm)
Pore radius R p 30 × 10 7 (30 nm)
Cell current density j 0 , A cm 2 1.5
PressureStandard
Cell temperature T, K273 + 80
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kulikovsky, A. Analytical Model for Concentration (Pressure) Impedance of a Low-Pt PEM Fuel Cell Oxygen Electrode. Membranes 2022, 12, 356. https://doi.org/10.3390/membranes12040356

AMA Style

Kulikovsky A. Analytical Model for Concentration (Pressure) Impedance of a Low-Pt PEM Fuel Cell Oxygen Electrode. Membranes. 2022; 12(4):356. https://doi.org/10.3390/membranes12040356

Chicago/Turabian Style

Kulikovsky, Andrei. 2022. "Analytical Model for Concentration (Pressure) Impedance of a Low-Pt PEM Fuel Cell Oxygen Electrode" Membranes 12, no. 4: 356. https://doi.org/10.3390/membranes12040356

APA Style

Kulikovsky, A. (2022). Analytical Model for Concentration (Pressure) Impedance of a Low-Pt PEM Fuel Cell Oxygen Electrode. Membranes, 12(4), 356. https://doi.org/10.3390/membranes12040356

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop