Comprehensive analytical model for polarization curve of a PEM fuel cell and experimental validation

• Published in: Chemical engineering science 2019-10, Vol.206, p.96-117
• Main Authors: Thosar, Aniket U., Agarwal, Harshal, Govarthan, S., Lele, Ashish K.
• Format: Article
• Language: English
• Subjects: Fuel cell equation; Modelling; PEMFC; Polarization curve
• ISSN: 0009-2509; 1873-4405
• DOI: 10.1016/j.ces.2019.05.022

•Generalized solution to FA-corrected MH model of the CCL is presented.•All the possible regimes of CCL operation are captured in the analytical model.•Derivation and experimental validation of ‘The PEM fuel cell equation’.•Experimental methodology to decipher regime of CCL operation is elucidated.•Procedure to determine the O2 diffusion coefficient in the CCL and GDL is provided. The kinetics of cathodic oxygen reduction reaction (ORR) in a proton exchange membrane fuel cell (PEMFC) is significantly modulated by the resistances for transport of reactants to the catalytic sites offered by different components of the fuel cell. This modulation governs the polarization curve of the PEMFC. Consequently, the various operating, geometric and material parameters of the fuel cell dictate the polarization curve. The effects of these parameters on the polarization curve over the entire range of current density, from zero to limiting current, can be predicted using detailed numerical simulations, which are however expensive. Analytical models, although simple can capture the essential details of physico-chemical processes occurring inside a PEMFC and are significantly inexpensive. In this article, we derive an analytical equation of the polarization curve which is valid over the entire range of current density. Specifically, the representative situation of a humidified low temperature PEMFC is considered wherein oxygen transport resistance in the cathode catalyst layer (CCL) is encountered at lower current density than proton transport resistance in the CCL. A novel experimental methodology is illustrated to confirm that this is indeed the case. Next, we elucidate a procedure to determine in-situ oxygen diffusion coefficients in the various domains of an operational PEMFC. Finally, it is shown that the analytical polarization curve predicted using these parameters is in excellent agreement with the experimental and numerically simulated polarization curves over the entire range of current density. The significance of this work is that the analytical model relates the performance of a PEMFC to all operating and geometric parameters as well as the average transport and kinetic properties of the materials used in its different components, without the need for computationally expensive numerical simulations. The model can therefore provide useful insights for enhancing the performance of PEMFC in different regimes of current density as well as for diagnostic purposes.