Modeling the cation transport in an operating polymer electrolyte fuel cell (PEFC)

https://doi.org/10.1016/j.ijhydene.2010.02.137Get rights and content

Abstract

A computational fluid dynamics model is developed to investigate the multicomponent cation transport in polymer electrolyte membranes and to predict the performance degradation of the polymer electrolyte fuel cell (PEFC) due to the cationic contamination. A Maxwell–Stefan approach is implemented by modifying the Nernst–Planck equations to model the multicomponent cationic species transport in the membrane. Langmuir isotherms are used to model the non-ideal species adsorption in the membrane. Cation transport model shows good agreement with the experimental data found in the literature. Following the validation of the cation transport model, it is incorporated in a PEFC model framework, which solves for mass, momentum, species and charge conservations. Both fuel side and air side impurities are considered in analyses, which show that for air side impurities current density drops from 0.9 to 0.7 A/cm2 whereas for fuel side impurities current density drops to impractical values as low as 0.05 A/cm2, at the steady-state. Effect of cationic impurity on water transport in the membrane is also investigated and found that cathode dry-out occurs due to decreased water generation in case of fuel side contamination.

Introduction

Owing to their environment friendly operation with low pollutant levels, quiet operation and the potential to reduce fossil fuel dependency, PEFCs are soaring as a very promising candidate for an alternative power generation method amidst the discussions of renewable energy and global warming.

Although PEFCs are becoming commercially available, the improvement of the technology relies on overcoming certain issues. Other than relatively high cost, reliability and durability of PEFCs still remain among the biggest hurdles. Reliability requires sustaining a robust operation and having a stable electrochemical performance as well as a structural integrity without mechanical or chemical deterioration. However, a PEFC may be required to operate in hazardous environments that may cause chemical, electrochemical or mechanical degradation sometimes to the extent of the complete failure of the cell.

For example, PEFCs employed in transportation applications rely on ambient air, which can be rich of salt suspensions. Such an environment is expected to have an impact on the cell performance as the fuel cell does not use pure air but that is contaminated with sea salt containing NaCl. Other than salt particles suspended in the air, a broad range of impurities can be found in the PEFC air and fuel streams as a result of residuals of the fuel reforming process and metal particles due to the corrosion in the stack and balance of plant materials. These impurities are then adsorbed in the ionomer in cationic forms. Cationic contamination of PEFCs may cause serious degradation in cell performance, which sometimes can be recoverable depending on the operating conditions and the exposure time, but many times it is irreversible [1], [2].

There have been efforts in characterizing the effects of cationic contaminants in PEFC operation [3], [4], [5], [6], as well as characterizing the properties of the cations in the membrane [7], [8], [9], [10], [11], [12], [13]. However, there is a little work in the literature describing the performance degradation mechanisms of the cationic contaminants. The first attempts are by Okada who developed a simple model to outline the effects of contaminants on the water transport in the membrane for both anode side [14] and cathode side [15] impurities. These studies neglect the transport of the cations and assume an infected zone in the membrane where a prescribed profile (either linear or constant) for cations is considered. Sodaye et al. [16] used Maxwell–Stefan approach to model the cationic exchange in a polymer electrolyte membrane in the presence of two different foreign cations. They only focus on the competitive transport of different cations in the membrane and do not consider fuel cell operation. Kienitz et al. developed one dimensional models for a fuel cell membrane [17] and a half cell [18] and predicted the cation concentration profiles in the studied domains. Their model assumes constant current density and does not provide an expression to describe the reaction kinetics; hence it does not relate the effects of cations on the performance loss substantially. Greszler et al. [19] developed a model for a hydrogen pump to investigate the effects of local proton concentration on the reversible potential of the reaction. Their study is significant as they put forward the thermodynamic aspects of the effects of contamination the first time. The most rigorous study up to now is provided by Weber and Delacourt [20] who used concentrated solution theory to model the interactions between each cation and water while they neglect the interactions between two cation pairs in a polymer electrolyte membrane. The model is utilized to predict the operation of a hydrogen pump with cationic impurities.

Albeit their individual strengths in explaining the cation transport in the polymer electrolyte membrane, all of the models cited above have some simplifications and none of them considers a complete fuel cell with both anode and cathode reactions. None of these models incorporate the coupled transport phenomena related to oxygen, hydrogen and water either. Since the transport properties of the contaminants depend on fuel cell operational variables such as water content and temperature, a comprehensive model requires a good representation of the transport phenomena in the full cell operation. All of the above studies are one dimensional only. In order to have an accurate representation of the fuel cell operation two-dimensional effects need to be accounted for, since 1-D models are not capable of predicting the water transport inside a PEFC. Due to the water accumulation under the bipolar plate rib area, there is a significant water concentration gradient in the in-plane direction (i.e. channel to the current collector) in addition to the through-plane direction (i.e. from anode to cathode). The gradient along the flow direction can be neglected if the fuel flow rates are high enough which minimizes the variations in the flow direction.

In this study we present a modeling framework to elucidate the cation transport in the polymer electrolyte membrane as well as to determine the mechanisms leading to performance degradation of a PEFC. We present the equations to describe the cation transport and validate them with the available experimental data. We then incorporate these equations into our fuel cell model, which solves for the transport phenomena pertaining to the fuel cell operation that impact the cation transport significantly. With our fuel cell model, we explain how the foreign cations affect the performance degradation of an operating fuel cell.

In this study we present the following for the first time:

  • i)

    Cation transport model is incorporated in a complete fuel cell model with catalyst layers (CLs) and gas diffusion layers (GDLs); a complete representation of the fuel cell transport phenomena is aimed so that the effects of fuel cell variables on the cation transfer is accounted more accurately.

  • ii)

    Cation uptake in the catalyst layers; competitive site occupation between the protons and the foreign cations in the catalyst layers is represented in the model.

  • iii)

    Two-dimensional effects in description of the cation transfer; we investigate the effects of water accumulation under the current collector rib.

Section snippets

Mathematical model

Fig. 1 shows a schematic of a 3-D PEFC geometry and our reduction and the details of the subdomains. In anode and cathode GDLs and CLs, where fluid is assumed in gaseous form; continuity, momentum and species equations are solved. For water in the membrane a separate species balance is considered. Electronic charge balance are solved in the GDLs and the CLs whereas ionic charge balance is coupled with the species balance to model the cation transfer in the CLs and the membrane.

Table 1 shows the

Numerical method

The models developed in this study are implemented in the commercial multiphysics software, COMSOL 3.4, which uses finite element method to discretize the partial differential equations. Analyses are carried out with the built-in non-iterative, direct linear system solver PARDISO, which is developed for the multiprocessing architectures, with nested dissection method as the preordering algorithm. 1800 quadrilateral Lagrangian elements are used to discretize the computational domain, which

Results

To validate our contamination model, we compare our results with the experimental data of Sodaye et al. [16] who studied multicomponent adsorption of Na+ and Cs+ in a Nafion 117 membrane. They used radiotracers to track the temporal concentration changes of two cations in the membrane. Using equilibrating solutions of NaCl and CsCl with different molarities they investigated the ion exchange process of the initially protonated Nafion membrane. With the solution containing NaCl/CsCl ratio of 1:1

Conclusion

A model to describe the multicomponent cation transport in the polymer electrolyte membrane is developed. Model results show good agreement with the experimental data. The model is capable of representing the competitive adsorption and diffusion of cationic species in the ion exchange process. After validating the cation transport model, the equations are implemented in a fuel cell model to predict the distributions of the cations in an operating fuel cell and predict the performance drop

Acknowledgment

We gratefully acknowledge the financial support from US Department of Energy through a cooperative research agreement DE-FE-07GO17020.

List of

Symbols

a
water activity
cw
water concentration in the membrane
Djk
binary diffusivities
Dw
water diffusion coefficient in the membrane
F
Faraday constant
i
transfer current density
i0
exchange current density
J
ionic current density vector
K
permeability
Mi
molecular weight
n
surface normal
N
mass flux
nd
electro-osmotic drag coefficient
p
pressure
Qe
volumetric electronic current source
R
universal gas constant
Rj
volumetric consumption of the jth species
T
temperature
u
velocity vector
um
mobility
U0
equilibrium voltage
V
inlet velocity
w

References (36)

  • G. Suresh et al.

    Study on physical and electrostatic interactions of counterions in poly(perfluorosulfonic) acid matrix: Characterization of diffusion properties of membrane using radiotracers

    Electrochimica Acta

    (2007)
  • H. Ju et al.

    A single-phase, non-isothermal model for PEM fuel cells

    Int J Heat Mass Transfer

    (2005)
  • S.M. Park et al.

    Effects of several trace contaminants on fuel cell performance

    (1980)
  • P.R. Hayter et al.

    The effect of contaminants in the fuel and air streams on the performance of a solid polymer fuel cell

    (1997)
  • F.A. Uribe et al.

    Effect of ammonia as potential fuel impurity on proton exchange membrane fuel cell performance

    J Electrochem Soc

    (2002)
  • H.J. Soto et al.

    Effect of transient ammonia concentrations on PEMFC performance

    Electrochem Solid-State Lett

    (2003)
  • M.S. Mikkola et al.

    The effect of NaCl in the cathode air stream on PEMFC performance

    Fuel Cells

    (2007)
  • H.L. Yeager et al.

    Cation and water diffusion in Nafion ion exchange membranes: Influence of polymer structure

    J Electrochem Soc

    (1981)
  • Cited by (33)

    • Evaluation of cathode contamination with Ca<sup>2+</sup> in proton exchange membrane fuel cells

      2018, Electrochimica Acta
      Citation Excerpt :

      The Nafion® perfluorosulfonic acid (PFSA) ionomer, composed of a carbon-fluorine backbone and side chains ending with sulfonic acid groups, is conventionally used to prepare catalyst layers by exploiting its binding, proton-conduction and O2 transport properties. Foreign cations can replace the ionomer H+ and occupy sulfonic acid sites by ion exchange [7–10], and this replacement of H+ by foreign cations increases the mass transport resistances of H+ and O2 in the ionomer. The model results obtained by Uddin et al. [10] showed that the presence of Na+ in the ionomer of the cathode catalyst layer, which decreases the O2 permeability, caused a decrease in the dissolved O2 concentration.

    View all citing articles on Scopus
    View full text