Next Article in Journal
Antibacterial Ferroelectric Hybrid Membranes Fabricated via Electrospinning for Wound Healing
Next Article in Special Issue
Cobalt-Doped Carbon Nitride Frameworks Obtained from Calcined Aromatic Polyimines as Cathode Catalyst of Anion Exchange Membrane Fuel Cells
Previous Article in Journal
The Mechanism of Metal Homeostasis in Plants: A New View on the Synergistic Regulation Pathway of Membrane Proteins, Lipids and Metal Ions
Previous Article in Special Issue
Design Strategies for Alkaline Exchange Membrane–Electrode Assemblies: Optimization for Fuel Cells and Electrolyzers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impedance Spectroscopy Measurements of Ionomer Film Oxygen Transport Resistivity in Operating Low-Pt PEM Fuel Cell

by
Tatyana V. Reshetenko
1,* and
Andrei Kulikovsky
2
1
Hawaii Natural Energy Institute, University of Hawaii, Honolulu, HI 96822, USA
2
Theory and Computation of Energy Materials (IEK-13), Institute of Energy and Climate Research Forschungszentrum Jülich GmbH, D-52425 Jülich, Germany
*
Author to whom correspondence should be addressed.
Membranes 2021, 11(12), 985; https://doi.org/10.3390/membranes11120985
Submission received: 22 November 2021 / Revised: 10 December 2021 / Accepted: 12 December 2021 / Published: 16 December 2021

Abstract

:
The work presents a model for local impedance of low-Pt proton exchange membrane fuel cells (PEMFCs), including cathode pore size distribution and O2 transport along pores and through a thin ionomer film covering Pt/C agglomerates. The model was applied to fit the local impedance spectra of low-Pt fuel cells operated at current densities from 100 to 800 mA cm−2 and recorded by a segmented cell system. Assuming an ionomer film thickness of 10 nm, the fitting returned the product of the dimensionless Henry’s constant of oxygen dissolution in ionomer KH by the oxygen diffusivity DN in the ionomer (KHDN). This parameter allowed us to determine the fundamental O2 transport resistivity R N through the ionomer film in the working electrode under conditions relevant to the realistic operation of PEMFCs. The results show that variation of the operating current density does not affect R N , which remains nearly constant at ≃0.4 s cm−1.

1. Introduction

At the end of 1838, Christian Friedrich Schoenbein immersed two platinum wires surrounded by glass vessels in liquid electrolyte and connected the electrodes to a Volta battery. Experiments with this simple system led him to the discovery of the fuel cell effect [1]. Ever since, due to its unique catalytic properties, Pt has been the material of choice for low temperature proton exchange membrane fuel cells (PEMFCs).
Unlike Schoenbein’s Pt wires, the modern PEMFC cathode is a porous composite material containing Pt nanoparticles on a carbon support, and a decrease in Pt loading in this system is a high-priority task for fuel cell technology. A decade ago it appeared that Pt reduction was a more complicated problem than anticipated, leading to a drop in fuel cell performance, which exceeded the expected loss due to a lowering of the Pt active area [2,3].
Further studies have revealed that the observed performance decline is caused by increased mass transport losses [3,4,5,6,7,8,9,10]. In the cathode catalyst layer (CCL), the oxygen pathway to the oxygen reduction reaction (ORR) site includes diffusive transport along the void pore, dissolution in the thin ionomer film and diffusion through the film to the Pt surface (Figure 1). Qualitatively, a lower number of ORR sites in the low-Pt electrode means higher oxygen and proton flux per site, making the ionomer film quite a significant barrier for oxygen transport to the site. Thus, in high-Pt cells, the effects of oxygen transport through the film are marginal, while in low-Pt cells, this transport gives quite a substantial contribution to the potential loss. Another explanation of the over-linear transport loss by the vulnerability of low-Pt electrodes to flooding has been suggested in [11].
The limiting current method is an established technique for the determination and analysis of O2 mass transport resistance in PEMFC electrodes [3,4,9,12,13,14,15,16]. A variation of oxygen concentration, total gas pressure or oxygen balance gas results in the separation of: (1) transport resistance due to molecular diffusion in gas channels and large pores; and (2) resistance originating from Knudsen diffusion in fine pores and diffusion from ionomer and liquid water films [3,4,9,14,15]. The effective thickness of the ionomer films estimated using O2 mass transport resistance was found to be significantly higher than that observed by transmission electron microscopy of 1–10 nm [6,7,9]. This finding suggested that the physical properties of the thin ionomer films in confinements such as water uptake and elastic modulus differed from the properties of bulk ionomer, and the additional interaction of ionomer and Pt negatively affected mass transport [6,7,9,17,18,19,20,21,22]. In spite of the fact that the limiting current approach is widely used for PEMFCs characterization, its implementation typically requires a cell with an active area of 1–5 cm2 and high stoichiometry of reagents to reach currents relevant for PEMFCs. The chosen conditions seem to be not quite representative for realistic fuel cell operation. The application of a segmented cell system and the distribution of limiting currents along a flow field can resolve this issue, since a cell with a high active area (>50 cm2) can be used and operated at realistic conditions [14,23,24]. However, it should be mentioned that a segmented cell system is not a widely available equipment, which makes its application limited. Therefore, developing methods and models for the determination of O2 mass transport parameters at conditions which are relevant for realistic PEMFCs operation and easily available for researchers is a very important task for the fuel cell community.
Currently, electrochemical impedance spectroscopy (EIS) seems to be the only easily accessible tool which allows us to obtain an in-situ diagnostic of operating PEMFCs. Moreover, the EIS method is very sensitive to oxygen transport processes in a fuel cell. Recently, we have reported a model for cell impedance which takes into account oxygen transport through the ionomer film [25]. The model links, in a 1d + 1d manner, the oxygen diffusion equations along the cylindrical void pore and through the ionomer film surrounding the pore. In this work, the model is extended for the processing of local impedance spectra, considering pore size distribution (PSD) in the electrode. To the best of our knowledge, this is the first impedance model of an operating low-Pt fuel cell that considers the PSD of the CCL.
This work is a continuation of our previous studies [25]. We determined the ionomer film oxygen transport resistivity R N based on a fitting of localized impedance spectra acquired from 10 segments of a low-Pt cell in a segmented cell setup. The experimental PSD of the cathode was approximated by the pores of nine characteristic radii. Assuming a film thickness of 10 nm, from each experimental run we obtained 10 local values of the product KHDN, where KH and DN are the Henry’s constant for oxygen solubility and oxygen diffusivity in the ionomer film, respectively. The mean over the cell surface value of KHDN and the standard deviation allowed us to calculate a statistically significant R N for every cell current density fixed in the experiments. In the range of 100 to 800 mA cm−2, it was determined that R N ≃ 0.4 s cm−1 is nearly independent of cell current density.

2. Materials and Methods

2.1. Electrochemical Evaluation

The impedance spectra were acquired using a segmented cell system designed at the Hawaii Natural Energy Institute (HNEI). The system includes a custom-built test station, segmented cell hardware, current/voltage sensors and a PXI data acquisition system [26]. The test station controls gas flow rates, humidification and back pressure, cell and gas line temperature and the overall cell current or voltage. The segments’ currents are measured but not controlled in order to mimic the experimental conditions of a non-segmented cell.
In this work, we used the Gore PRIMEA (A510.1/M715.18/C510.1) catalyst coated membrane (CCM) with an active area of 100 cm2. The Pt/C catalyst loading for anode and cathode was 0.1 mgPt cm−2. The electrode thickness varied from 3 to 4 μm. Sigracet 25BC, with a thickness of 220–230 μm, was employed as the gas diffusion layer (GDL). To keep an optimal membrane electrode assembly (MEA) compression ratio of 25–30%, we applied Teflon gaskets with a thickness of 125 μm on both sides of the cell.
Electrochemical measurements were performed using modified 100 cm2 hardware which had segmented and non-segmented flow fields. The segmented flow field and current collector were separated into 10 segments with an area of 7.6 cm2, connected in order from segment 1, located at the gas inlet, to segment 10, at the outlet. Due to the segmentation of the cathode flow field, the resulting active MEA area was 76 cm2. The segmented and unchanged flow fields had the same 10-channel parallel serpentine design and were arranged in co-flow configuration. Details of the segmented cell system operation are given in [26,27].
Impedance spectra were measured under galvanostatic control of the whole cell with 11 steps/decade; the cell operating conditions are listed in Table 1. The frequency range was from 10 kHz to 0.05 Hz. The amplitude of current perturbation was adjusted to obtain a cell voltage response of 10 mV or less. The sampling rate of the current/voltage response from 10 segments was 1 MHz.

2.2. Impedance Model

The cell used in the experiments has been modeled as a cell with the straight cathode channel separated into 10 segments (Figure 2a). The model of the local segment impedance Zseg is based on a transient CCL performance model, which includes a mass conservation equation for oxygen transport in the GDL. The CCL is modeled as a system of cylindrical pores; each pore is inserted into a coaxial Pt/C layer and separated from this layer by a thin ionomer film (Figure 2b). Oxygen is transported along the pore axis and in the radial direction through the ionomer film to the ionomer/Pt interface, where the oxygen reduction reaction occurs. Model equations have been described in details in [25]; here, for completeness, we list the basic equations.

2.3. Model Equations

Oxygen transport along the void pore (Figure 2b) is described by the transient diffusion equation:
c t D p 2 c t 2 = 2 N N , p R p , c x x = 0 = 0 , c l t = c 1
where c1 is the oxygen concentration at the CCL/GDL interface, lt is the CCL thickness (the pore length), Dp is oxygen diffusion coefficient in the pore,
N N , p = D N c N r r = R p +
is the oxygen flux in the ionomer film at the pore/film interface, and cN is the oxygen concentration in ionomer film.
Oxygen transport through the ionomer film along the radial direction is described by:
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 = i * R p 2 4 F c N , m c i n exp η b
where KH is the dimensionless Henry’s constant (mol/mol) for oxygen solubility in ionomer, DN is the oxygen diffusivity in the ionomer film, Rm is the Pt/C tube radius (Figure 2b), c N , m c N R m , c i n is the reference (inlet) oxygen concentration, b is the ORR Tafel slope, and i * is the ORR exchange current density.
Proton current conservation in the ionomer film and the Ohm’s law for this current lead to an equation for the positive by convention ORR overpotential η :
C d l η t σ p 2 η x 2 = i * c N , m c h i n exp η b
where Cdl is the double layer volumetric capacitance, and σp is the ionomer film proton conductivity.
Oxygen transport in the gas diffusion layer is described by:
c b t D b 2 c b x 2 = 0 , c b x x = l t + = c x x = l t , c b l t + l b = c h
where cb is the oxygen concentration in the GDL, Db and lb are the GDL oxygen diffusivity and thickness, respectively, and ch is the oxygen concentration in the air channel. The latter value obeys to plug-flow mass transport equation along the channel coordinate z:
c h t + v c h z = D b h c b x x = l b
where v is the air flow velocity and h is the channel depth.
As can be seen, transient diffusion equations for oxygen transport along the pore, Equation (1), and through the ionomer film, Equation (3), are linked in a 1d+1d manner. The ORR rate at the Pt/C interface is described by the Tafel law. The proton charge conservation equation along the ionomer film, Equation (4), takes into account finite proton conductivity of the film and charging/discharging of the double layer. Radial oxygen flux in the film at the pore/film interface represents a sink in the oxygen mass transport equation along the void pore [25].
The total oxygen flux along the void pore is proportional to the square of the pore radius, while this flux through the ionomer film covering the side surface of the pore is proportional to the first power of the radius. Evidently, the electrode O2 transport properties strongly depend on the PSD [28,29,30]. The PSD of a PRIMEA Gore CCM was presented in [31]. Since in the current work we applied the PRIMEA CCM as well, the reported PSD adequately described the pore size spectrum of the CCL under the study. To accelerate calculations, the PSD measured in [31] has been approximated by a stepwise function with the nine characteristic radii (Figure 3).
The pore radii and their relative fractions in the electrode volume calculated from the experimental PSD are depicted in Table 2. The CCL impedance Zccl,n in the n-th segment has been calculated as an impedance of these parallel pores:
1 Z c c l , n = k = 1 9 w k Z k , n
where Z k , n and w k are the pore impedance and volume fraction, respectively, of the k-th pore in the spectrum.
Linearization and Fourier-transform of the transport equations lead to a system of equations for the AC perturbation amplitudes in the ω-space [25]. This system is solved numerically for every individual segment, taking into account transport of oxygen concentration perturbations with the channel flow (Figure 4). Air flow in the cathode channel transports local oxygen concentration perturbations c h 1 to the segments located downstream. Thus, the local c h 1 sums up with the transported perturbations, which affect the local impedance.
A custom parallel Python code was developed using the message passing interface library. Calculated segment impedance Zseg was fitted to the local experimental spectrum using the SciPy least-squares procedure least_squares with the method = ‘trf’ option [32]. The nine-pore model described above for each of the 10 segments requires 90 cores on a parallel cluster. Detailed description of the parent whole-cell impedance model can be found in [25].

3. Results and Discussion

The local impedance spectra were measured at the operating conditions listed in Table 1. Figure 5 demonstrates experimental and fitted Nyquist spectra for the current density of 400 and 800 mA cm−2. Typical impedance spectra for a PEM fuel cell are represented by a semicircle at high-frequency, reflecting a charge transfer, proton and oxygen transport within the cathode layer, and a low-frequency (LF) arc caused by oxygen transport in the channel [33,34]. The local EIS results show an increase in the LF arc size from the cell inlet (segment 1) to the outlet (segment 10) and clearly illustrate the effects of O2 depletion along the flow field (Figure 5a). At the outlet segments 8 to 10, the LF part of the spectra becomes noisy, and the quality of the fitting decays (Figure 5). However, this arc represents oxygen transport in the gas channel, while the characteristic frequency of ionomer film impedance is at least an order of magnitude higher [35]. Overall, the results show that the proposed model ensures a good fitting of the local EIS data at low and high current densities.
The model equations include pore radius Rp and “metal” radius Rm; the film thickness is l N = R m R p (Figure 2b). Code testing shows that the fitting is more sensitive to the ratio lN/DN, rather than to lN and DN separately. On the other hand, equations in Section 2.3 show that DN appears in the model only as a product KHDN, where KH is the dimensionless Henry’s constant for oxygen solubility in the ionomer film. To reliably obtain the parameter KHDN from the fitting, the ionomer film thickness was fixed at lN = 10 nm. This thickness value agrees well with our previous measurements [25] and the other literature sources [6,7,9,36,37,38]. The starting value for iterations of KHDN = 6.76 × 10−9 cm2 s−1 was taken, which is the dimensionless Henry’s constant for oxygen solubility in water at 80 °C, multiplied by the oxygen diffusivity in bulk Nafion of 10−6 cm2 s−1 [39].
The product KHDN that resulted from the fitting is nearly constant as the cell current density increases (Figure 6). The only imaginable difference between the low- and high-current cell operation is the amount of liquid water in the catalyst layer, which increased with the cell operating current. The authors [40,41,42,43,44] measured the ionomer film oxygen diffusivity DN using electrochemical methods and reported DN and O2 permeability growth with the relative humidity, which agreed with the classic works of T. Sakai and K. Broka [45,46,47].
Moreover, O2 mass transport resistance was shown to decrease with inlet gas humidification, and this effect was more pronounced during transition from dry conditions to 50% RH, while a further increase in reagent inlet humidification to 100% RH resulted in a slight reduction of this parameter [7,9,19,24,48,49,50]. In our case, no significant variation of KHDN with the current was detected, most likely due to the chosen operating conditions when fully humidified H2 was supplied to the anode and 50% RH air was fed to the cathode, ensuring proper hydration of both electrodes.
A single-pore oxygen transport model gives the following formula for the transport resistivity of ionomer film [35]:
R N = R p l N 2 K H D N l t
The film transport resistivity R N in the CCL was calculated using Equation (8) with the mean over the pore size distribution Rp:
R p = k = 1 9 w k R p , k
here, R p , k is the radius of the k-th pore in the spectrum. The determined film resistivity is close to the literature data of R N ≃ 0.3 s cm−1, corresponding to Pt loading of 0.1 mgPt cm−2 [3,4,9]. In the studied range of currents, R N is found to be nearly constant at ≃0.4 s cm−1 (Figure 7).
The effects of operating currents on oxygen transport resistance was firstly evaluated by D. Caulk and D. Baker using the limiting current method [51,52]. The authors reported that the oxygen transport resistance was constant for current density below 1 A cm−2, after which it increased rapidly with current density until 1.6–2.0 A cm−2, where it reached a second plateau. The observation was explained by the water production rate and its condensation to liquid at the cathode. The initial constant transport resistance corresponded to dry conditions when all produced water could be removed from the cathode in vapor form. As the current density increased, the water vapor pressure reached saturation, and water started to condense in the CCL porous structure leading to a growth of the transport resistance. At some point, the transport resistance reached the second plateau, indicating the maximum level of liquid saturation of the cathode electrode when the produced water could be removed mainly in liquid form. These results were confirmed by using in-situ neutron imaging [53] and synchrotron X-ray radiation [54,55] showing the formation of liquid water in the CCL and linking these data with results on oxygen mass transport resistance obtained by the limiting current method.
In our case, we operated at current densities below 1.0 A cm−2, and we could assume that the produced water was transferred from the catalyst surface to the gas channel in vapor form and did not noticeably flood the CCL [56]. Since initially the CCL was properly humidified, and there were no conditions for flooding the ionomer, film resistivity should not have been significantly affected by the operating current. We expect that a more accurate R N could be obtained from impedance measurements using a more detailed approximation of an experimental PSD.

4. Conclusions

A numerical model for local cell impedance, taking into account oxygen transport through the ionomer covering Pt/C agglomerates in the CCL and pore size distribution in the electrode, has been developed. The model has been fitted to experimental spectra acquired from10 segments in the segmented cell setup under conditions relevant to realistic PEMFC operation. Assuming an ionomer film thickness of 10 nm, from each experimental run, the fitting returns 10 values of the product of Henry’s constant for oxygen solubility in the ionomer film KH by oxygen diffusivity DN in the film (KHDN). Knowing KHDN allows us to calculate the statistical mean and confidence interval for the oxygen transport resistivity through the ionomer film in the range of the current densities from 100 to 800 mA cm−2. It was determined that the R N remains nearly constant at 0.4 s cm−1 for the chosen current density values.

Author Contributions

Conceptualization, A.K. and T.V.R.; Investigation, A.K. and T.V.R.; Methodology, A.K. and T.V.R.; Software, A.K.; Writing—Original Draft Preparation, A.K.; Writing—Review and Editing, A.K. and T.V.R. All authors have read and agreed to the published version of the manuscript.

Funding

T.V. Reshetenko gratefully acknowledges funding from the US Office of Naval Research (N00014-18-1-2127).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

The authors are thankful to Günter Randolf for valuable help in the system operation.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

bORR Tafel slope, V
CdlDouble layer volumetric capacitance, F cm−3
cOxygen molar concentration in the pore, mol cm−3
cbOxygen molar concentration in the GDL, mol cm−3
chOxygen molar concentration in the channel, mol cm−3
cNOxygen molar concentration in the ionomer (Nafion) film, mol cm−3
cinReference (inlet) oxygen concentration, mol cm−3
DbOxygen diffusion coefficient in the GDL, cm2 s−1
DpOxygen diffusion coefficient in the pore, cm2 s−1
DNOxygen diffusion coefficient in the Nafion film, cm2 s−1
FFaraday constant, C mol−1
hChannel depth, cm
i * ORR volumetric exchange current density, A cm−3
jLocal proton current density along the pore, A cm−2
j0Cell current density, A cm−2
KHHenry’s constant for oxygen solubility in the film, mol/mol
ltPore length (CCL thickness), cm
lNIonomer film thickness, cm
lbGDL thickness, cm
NNOxygen flux in the ionomer film, mol cm−2s−1
RmRadius of a Pt/C tube, cm
RpPore radius, cm
Rp,kRadius of k-th pore in the spectrum, cm
R N Ionomer film O2 transport resistivity, s cm−1
rRadial coordinate, cm
vFlow velocity in the cathode channel, cm s−1
w k Volume fraction of pores with the k-th radius
xCoordinate along the pore, cm
zCoordinate along the channel, cm
Zccl, nCCL impedance in the n-th segment, Ohm cm2
Zk, nPore impedance with k-th radius, Ohm cm2
ZsegSegment impedance, Ohm cm2
Subscripts:
0Membrane/CCL interface
1CCL/GDL interface
tCatalyst layer
Superscripts:
0Steady-state value
1Small-amplitude perturbation
Greek:
ηORR overpotential, positive by convention, V
σpIonomer film proton conductivity, S cm−1
ωAngular frequency of the AC signal, s−1
Abbreviations:
PEMFCProton exchange membrane fuel cell
CCLCathode catalyst layer
ORROxygen reduction reaction
EISElectrochemical impedance spectroscopy
PSDPore size distribution
CCMCatalyst coated membrane
GDLGas diffusion layer
MEAMembrane electrode assembly
LFLow frequency

References

  1. Schoenbein, C.F. On the voltaic polarization of certain solid and fluid substances. Phil. Mag. 1839, 14, 43–45. [Google Scholar] [CrossRef]
  2. 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]
  3. 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]
  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 low platinum 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. Comm. 2017, 79, 14–17. [Google Scholar] [CrossRef] [Green Version]
  9. 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]
  10. Choo, M.-J.; Oh, K.-H.; Park, J.-K.; Kim, H.-T. Analysis of oxygen transport in cathode catalyst layer of low-Pt-loaded fuel cells. ChemElectroChem 2015, 2, 382–388. [Google Scholar] [CrossRef]
  11. Muzaffar, T.; Kadyk, T.; Eikerling, M. Tipping water balance and the Pt loading effect in polymer electrolyte fuel cells: A model-based analysis. Sustain. Energy Fuels 2018, 2, 1189–1196. [Google Scholar] [CrossRef] [Green Version]
  12. Beuscher, U. Experimental method to determine the mass transport resistance of a polymer electrolyte fuel cell. J. Electrochem. Soc. 2006, 153, A1788–A1793. [Google Scholar] [CrossRef]
  13. Mashio, T.; Ohma, A.; Yamamoto, S.; Shinohara, K. Analysis of reactant gas transport in a catalyst layer. ECS Trans. 2007, 11, 529–540. [Google Scholar] [CrossRef]
  14. Baker, D.R.; Caulk, D.; Neyerlin, K.C.; Murphy, M.W. Measurement of oxygen transport resistance in PEM fuel cells by limiting current methods. J. Electrochem. Soc. 2009, 156, B991–B1003. [Google Scholar] [CrossRef]
  15. Reshetenko, T.; St-Pierre, J. Separation method for oxygen mass transport coefficient in gas and ionomer phases in PEMFC GDE. J. Electrochem. Soc. 2014, 161, F1089–F1100. [Google Scholar] [CrossRef]
  16. Tabe, Y.; Sakai, T.; Iiri, T.; Hayashi, T.; Chikahisa, T. Experimental evaluation of dominant transport resistances of oxygen in catalyst layers of PEFC. ECS Trans. 2017, 80, 205–214. [Google Scholar] [CrossRef]
  17. Kusoglu, A.; Weber, A.Z. New insights into perfluorinated sulfonic-acid ionomers. Chem. Rev. 2017, 117, 987–1104. [Google Scholar] [CrossRef]
  18. Page, K.A.; Kusoglu, A.; Stafford, C.M.; Kim, S.R.J.; Kline, R.J.; Weber, A.Z. Confinement-driven increase in ionomer thin-film modulus. Nano Lett. 2014, 14, 2299–2304. [Google Scholar] [CrossRef]
  19. Jinnouchi, R.; Kudo, K.; Kitano, N.; Morimoto, Y. Molecular dynamics simulations on O2 permeation through Nafion ionomer on Platinum surface. Electrochim. Acta 2016, 188, 767–776. [Google Scholar] [CrossRef]
  20. Iiri, T.; Tabe, Y.; Ckikahisa, T. Experimental analysis of oxygen transport resistance for different types of ionomer in PEFC catalyst. ECS Trans. 2018, 86, 141–150. [Google Scholar] [CrossRef]
  21. Jomori, S.; Komatsubara, K.; Nonoyama, N.; Kato, M.; Yoshida, T. An experimental study of the effects of ioperational history on activity changes in a PEMFC. J. Electrochem. Soc. 2013, 160, F1067–F1073. [Google Scholar] [CrossRef]
  22. Liu, H.; Epting, W.K.; Lister, S. Gas transport resistance in polymer electrolyte thin films on oxygen reduction reaction catalysts. Langmuir 2015, 31, 9853–9858. [Google Scholar] [CrossRef]
  23. Reshetenko, T.; Ben, B.L. Impact of a gas diffusion layer’s structural and textural properties on oxygen mass transport resistance in the cathode and performance of proton exchange membrane fuel cells. Electrochim. Acta 2021, 371, 137752. [Google Scholar] [CrossRef]
  24. Reshetenko, T.; Ben, B.L. Exploration of operating conditions on oxygen mass transport resistance and performance of PEM fuel cells: Effects of inlet gas humidification. Electrochem. Sci. Adv. 2021, e2100134. [Google Scholar] [CrossRef]
  25. Kulikovsky, A.; Reshetenko, T. Nafion film transport properties in a low-Pt PEM fuel cell: Impedance spectroscopy study. RSC Adv. 2019, 9, 38797–38806, Erratum in RSC Adv. 2021, 11, 6764–6765. [Google Scholar] [CrossRef] [Green Version]
  26. Reshetenko, T.V.; Bender, G.; Bethune, K.; Rocheleau, R. Systematic study of back pressure and anode stoichiometry effects on spatial PEMFC performance distribution. Electrochim. Acta 2011, 56, 8700–8710. [Google Scholar] [CrossRef]
  27. Reshetenko, T.V.; Bender, G.; Bethune, K.; Rocheleau, R. A segmented cell approach for studying the effects of serpentine ow field parameters on PEMFC current distribution. Electrochim. Acta 2012, 88, 571–579. [Google Scholar] [CrossRef]
  28. Liu, J.; Eikerling, M. Model of cathode catalyst layers for polymer electrolyte fuel cells: The role of porous structure and water accumulation. Electrochim. Acta 2008, 53, 4435–4446. [Google Scholar] [CrossRef]
  29. Yu, Z.; Carter, R.N.; Zhang, J. Measurements of pore size distribution, porosity, effective oxygen diffusivity, and tortuosity of PEM fuel cell electrodes. Fuel Cells 2012, 12, 557–565. [Google Scholar] [CrossRef]
  30. Zhou, J.; Putz, A.; Secanell, M. A mixed wettability pore size distribution based mathematical model for analyzing two-phase flow in porous electrodes. J. Electrochem. Soc. 2017, 164, F530–F539. [Google Scholar] [CrossRef] [Green Version]
  31. Yu, H.; Bonville, L.; Maric, R. Analysis of H2/Air polarization curves: The influence of low Pt loading and fabrication process. J. Electrochem Soc. 2018, 165, F272–F284. [Google Scholar] [CrossRef]
  32. Scipy Optimize Least Squares. 2018. Available online: https://docs.scipy.org/doc/scipy/reference/generated/scipy.optimize.least_squares.html (accessed on 1 January 2020).
  33. Schneider, I.A.; Freunberger, S.A.; Kramer, D.; Wokaun, A.; Scherer, G.G. Oscillations in gas channels. Part, I. The forgotten player in impedance spectroscopy in PEFCs. J. Electrochem. Soc. 2007, 154, B383–B388. [Google Scholar] [CrossRef]
  34. Springer, T.E.; Zawodzinski, T.A.; Wilson, M.S.; Gottesfeld, S. Characterization of polymer electrolyte fuel cells using AC impedance spectroscopy. J. Electrochem. Soc. 1996, 143, 587–599. [Google Scholar] [CrossRef]
  35. Kulikovsky, A. Impedance and resistivity of low-Pt cathode in a PEM fuel cell. J. Electrochem. Soc. 2021, 168, 044512. [Google Scholar] [CrossRef]
  36. Lopez-Haro, M.; Guétaz, L.; Printemps, T.; Morin, A.; Escribano, S.; Jouneau, P.-H.; Bayle-Guillemaud, P.; Chandezon, F.; Gebel, G. Three-dimensional analysis of Nafion layers in fuel cell electrodes. Nat. Comm. 2013, 5, 5229. [Google Scholar] [CrossRef]
  37. Cetinbas, F.C.; Wang, X.; Ahluwalia, R.K.; Kariuki, N.N.; Winarski, R.; Yang, Z.; Sharman, J.; Myers, D.J. Analysis and transport resistances of low-Platinum-loaded PEFC electrodes. J. Electrochem. Soc. 2017, 164, F1595–F1607. [Google Scholar] [CrossRef]
  38. Suzuki, T.; Okada, S.; Tsushima, S. Analysis of ionomer distribution and Pt/C agglomerate size in catalyst layers by two-stage ion-beam processing. J. Electrochem. Soc. 2020, 167, 124513. [Google Scholar] [CrossRef]
  39. Sethuraman, V.A.; Khan, S.; Jur, J.S.; Haug, A.T.; Weidner, J. Measuring oxygen, carbon monoxide and hydrogen sulfide diffusion coefficient and solubility in Nafion membranes. Electrochim. Acta 2009, 54, 6850–6860. [Google Scholar] [CrossRef] [Green Version]
  40. Novitski, D.; Holdcroft, S. Determination of O2 mass transport at the Pt/PFSA ionomer interface under reduced relative humidity. ACS Appl. Mater. Interfaces 2015, 7, 27314–27323. [Google Scholar] [CrossRef]
  41. Sambandam, S.; Parrondo, J.; Ramani, V. Estimation of electrode ionomer oxygen permeability and ionomer-phase oxygen transport resistance in polymer electrolyte fuel cells. Phys. Chem. Chem. Phys. 2013, 15, 14994–15002. [Google Scholar] [CrossRef]
  42. Takamura, Y.; Nakashima, E.; Yamada, H.; Tasaka, A.; Inaba, M. Effects of temperature and relative humidity on oxygen permeation in Nafion® and sulfonated poly(arylene ether sulfone). ECS Trans. 2008, 16, 881–889. [Google Scholar] [CrossRef]
  43. Ogumi, Z.; Takehara, Z.; Yoshizawa, S. Gas permeation in SPE method: I. Oxygen permeation through Nafion and Neosepta. J. Electrochem. Soc. 1984, 131, 769–773. [Google Scholar] [CrossRef]
  44. Ogumi, Z.; Kuroe, T.; Takehara, Z. Gas permeation in SPE method: II. Oxygen and hydrogen permeation through Nafion. J. Electrochem. Soc. 1985, 132, 2601–2605. [Google Scholar] [CrossRef]
  45. Sakai, T.; Takenaka, H.; Wakabayashi, N.; Kawami, Y.; Torikai, E. Gas permeation properties of solid polymer electrolyte (SPE) membranes. J. Electrochem. Soc. 1985, 132, 1328–1332. [Google Scholar] [CrossRef]
  46. Sakai, T.; Takenaka, H.; Torikai, E. Gas diffusion in the dried and hydrated Nafions. J. Electochem. Soc. 1986, 133, 88–92. [Google Scholar] [CrossRef]
  47. Broka, K.; Ekdunge, P. Oxygen and hydrogen permeation properties and water uptake of Nafion® 117 membrane and recast film for PEM fuel cell. J. Appl. Electrochem. 1997, 27, 117–123. [Google Scholar] [CrossRef]
  48. Oh, H.; Lee, Y.-I.; Lee, G.; Min, K.; Yi, J.-S. Experimental dissection of oxygen transport resistance in the components of a polymer electrolyte membrane fuel cell. J. Power Sources 2017, 345, 67–77. [Google Scholar] [CrossRef]
  49. Shen, S.; Cheng, X.; Wang, C.; Yan, X.; Ke, C.; Yin, J.; Zhang, J. Exploaration of significant influence of the operating conditions on the local O2 transport in proton exchange membrane fuel cells (PEMFCs). Phys. Chem. Chem. Phys. 2017, 19, 26221–26229. [Google Scholar] [CrossRef]
  50. Du, F.; Dao, T.A.; Peitl, P.V.J.; Bauer, A.; Preuss, K.; Bonastre, A.M.; Sharman, J.; Spikes, G.; Perchthaler, M.; Schmidt, T.J.; et al. effects of PEMFC operational history under dry/wet conditions on additional voltage losses due to ionomer migration. J. Electrochem. Soc. 2020, 167, 144513. [Google Scholar] [CrossRef]
  51. Caulk, D.A.; Baker, D.R. Heat and water transport in hydrophobic diffusion media of PEM fuel cells. J. Electrochem. Soc. 2010, 157, B1237–B1244. [Google Scholar] [CrossRef]
  52. Caulk, D.A.; Baker, D.R. Modeling two-phase water transport in hydrophobic diffusion media for PEM fuel cells. J. Electrochem. Soc. 2011, 158, B384–B393. [Google Scholar] [CrossRef]
  53. Owejan, J.P.; Trabold, T.A.; Mench, M.M. Oxygen transport resistance correlated to liquid water saturation in the gas diffusion layer of PEM fuel cells. Int. J. Heat Mass Transfer 2014, 71, 585–592. [Google Scholar] [CrossRef]
  54. Muirhead, D.; Banerjee, R.; George, M.G.; Ge, N.; Shrestha, P.; Liu, H.; Lee, J.; Bazylak, A. Liquid water saturation and oxygen transport resistance in polymer electrolyte membrane fuel cell gas diffusion layers. Electrochim. Acta 2018, 274, 250–263. [Google Scholar] [CrossRef] [Green Version]
  55. Muirhead, D.; Banerjee, R.; Lee, J.; George, M.G.; Ge, N.; Liu, H.; Chevalier, S.; Hinebaugh, J.; Han, K.; Bazylak, A. Simultaneous characterization of oxygen transport resistance and spatially resolved liquid water saturation at high-current density of polymer membrane fuel cells with varied cathode humidity. Int. J. Hydrogen Energy 2017, 42, 29472–29483. [Google Scholar] [CrossRef]
  56. Owejan, J.P.; Owejan, J.E.; Gu, W.; Trabold, T.A.; Tighe, T.W.; Mathias, M.F. Water transport mechanisms in PEMFC gas diffusion layer. J. Electrochem. Soc. 2010, 157, B1456–B1464. [Google Scholar] [CrossRef]
Figure 1. Schematic of the cathode catalyst layer structure in PEM fuel cells. Pt/C agglomerates are covered by a thin ionomer film. To reach Pt sites, oxygen transported along the void pore must penetrate through the ionomer.
Figure 1. Schematic of the cathode catalyst layer structure in PEM fuel cells. Pt/C agglomerates are covered by a thin ionomer film. To reach Pt sites, oxygen transported along the void pore must penetrate through the ionomer.
Membranes 11 00985 g001
Figure 2. Schematic of the cell with segmented electrodes and gas diffusion layer (a). Sketch of a single pore used to model cathode catalyst layer impedance (b).
Figure 2. Schematic of the cell with segmented electrodes and gas diffusion layer (a). Sketch of a single pore used to model cathode catalyst layer impedance (b).
Membranes 11 00985 g002
Figure 3. Experimental [31] (dotted line) and approximate model (dashed line) pore size distributions in the cathode catalyst layer.
Figure 3. Experimental [31] (dotted line) and approximate model (dashed line) pore size distributions in the cathode catalyst layer.
Membranes 11 00985 g003
Figure 4. Schematic of mixing of the oxygen concentration AC perturbations c h 1 due to the air flow in the channel. The perturbations c h 1 produced upstream from the given segment are transported down the channel and sum up with the local perturbation. This interference affects local impedance.
Figure 4. Schematic of mixing of the oxygen concentration AC perturbations c h 1 due to the air flow in the channel. The perturbations c h 1 produced upstream from the given segment are transported down the channel and sum up with the local perturbation. This interference affects local impedance.
Membranes 11 00985 g004
Figure 5. Experimental (points) and fitted model (lines) local Nyquist spectra from the individual segments at 400 (a) and 800 mA cm−2 (b). Air inlet is at segment 1, outlet is at segment 10.
Figure 5. Experimental (points) and fitted model (lines) local Nyquist spectra from the individual segments at 400 (a) and 800 mA cm−2 (b). Air inlet is at segment 1, outlet is at segment 10.
Membranes 11 00985 g005
Figure 6. The product KHDN of Henry’s constant for oxygen solubility in ionomer film by oxygen diffusivity in the film vs. cell current density.
Figure 6. The product KHDN of Henry’s constant for oxygen solubility in ionomer film by oxygen diffusivity in the film vs. cell current density.
Membranes 11 00985 g006
Figure 7. Ionomer film oxygen transport resistivity.
Figure 7. Ionomer film oxygen transport resistivity.
Membranes 11 00985 g007
Table 1. Cell and operating parameters. A/C stands for anode/cathode.
Table 1. Cell and operating parameters. A/C stands for anode/cathode.
ParameterValue
Catalyst loading A/C, mgPt cm−20.1/0.1
CCL thickness lt, cm3.0 × 10−4
GDL thickness lb, cm230 × 10−4
Ionomer film thickness lN, cm10.0 × 10−7
Flow stoichiometry A/C2/4
Relative humidity A/C100%/50%
Back pressure, kPa150
Cell temperature, K353
Table 2. Parameters of the stepwise PSD in Figure 3.
Table 2. Parameters of the stepwise PSD in Figure 3.
Parameter123456789
Pore radius, nm16.535.445.054.365.175.188.8130152
Pore volume fraction0.1990.1000.1270.1370.1120.06850.07250.1110.0727
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Reshetenko, T.V.; Kulikovsky, A. Impedance Spectroscopy Measurements of Ionomer Film Oxygen Transport Resistivity in Operating Low-Pt PEM Fuel Cell. Membranes 2021, 11, 985. https://doi.org/10.3390/membranes11120985

AMA Style

Reshetenko TV, Kulikovsky A. Impedance Spectroscopy Measurements of Ionomer Film Oxygen Transport Resistivity in Operating Low-Pt PEM Fuel Cell. Membranes. 2021; 11(12):985. https://doi.org/10.3390/membranes11120985

Chicago/Turabian Style

Reshetenko, Tatyana V., and Andrei Kulikovsky. 2021. "Impedance Spectroscopy Measurements of Ionomer Film Oxygen Transport Resistivity in Operating Low-Pt PEM Fuel Cell" Membranes 11, no. 12: 985. https://doi.org/10.3390/membranes11120985

APA Style

Reshetenko, T. V., & Kulikovsky, A. (2021). Impedance Spectroscopy Measurements of Ionomer Film Oxygen Transport Resistivity in Operating Low-Pt PEM Fuel Cell. Membranes, 11(12), 985. https://doi.org/10.3390/membranes11120985

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