A 2D through-the-Membrane Transient Model for Polymer Electrolyte Membrane Fuel Cells
Authors: Goshtasbi Alireza , Pence Benjamin , Ersal Tulga 
Published: February 2020

Effective and efficient water and thermal management remain an integral part of the research and development efforts aimed at improving performance and durability of polymer electrolyte membrane (PEM) fuel cells. Understanding the transport phenomena occurring in the cell is key to developing successful water and thermal management strategies. To this end, recent experimental observations have shed light on some aspects of the transport in the cell, such as the phase change induced flow [1], and have significantly contributed to the current understanding in this area. Furthermore, mathematical modeling has been an indispensable tool in this endeavor and has provided insight into the physical phenomena in the cell. However, most of the models developed to date for this purpose have focused on steady state performance [2, 3]. While successful transient modeling efforts have been reported in the literature [4], a complete understanding of the dynamic behavior of the cell requires further research.
Motivated by the automotive industry’s requirements to better understand this dynamic behavior, in this talk we present our latest modeling efforts emphasizing the transient response of the cell. In particular, we present a fully coupled two-dimensional finite element model with two-phase flow and non-isothermal effects. The model resolves transport across various length scales, including the diffusion media and the catalyst layers, where a recently developed homogenous model is employed to account for the local oxygen transport resistance [5]. Moreover, anisotropy in material properties is considered to model in-plane variations. Oxygen reduction reaction (ORR) is described using the modified double-trap model that captures doubling of the Tafel slope [6]. ORR dependency on relative humidity is also taken into account [7, 8]. Values reported in the literature for time constants of membrane water sorption and desorption are discussed. Based on this discussion, a factor scaling the time constant of the original diffusion equation used for membrane water transport is introduced and its importance is discussed through transient modeling results.
In general, the model can be used to investigate the transient response of the fuel cell under various operating conditions. A steady state version of the model can help in making design decisions to improve cell performance. Finally, the model can serve as a benchmark to parameterize and validate more computationally efficient models that can be used for real-time applications.

Acknowledgement:

Financial support for this work was provided by Ford Motor Company through the Ford University Alliance Program.

References:

[1] I. V. Zenyuk, A. Lamibrac, J. Eller, D. Y. Parkinson, F. Marone, F. N. Büchi, and A. Z. Weber, "Investigating Evaporation in Gas Diffusion Layers for Fuel Cells with X-ray Computed Tomography," The Journal of Physical Chemistry C, vol. 120, no. 50, pp. 28701-28711, 2016.
[2] M. Moore, P. Wardlaw, P. Dobson, J. Boisvert, A. Putz, R. Spiteri, and M. Secanell, "Understanding the effect of kinetic and mass transport processes in cathode agglomerates," Journal of the Electrochemical Society, vol. 161, no. 8, pp. E3125-E3137, 2014.
[3] M. Bhaiya, A. Putz, and M. Secanell, "Analysis of non-isothermal effects on polymer electrolyte fuel cell electrode assemblies," Electrochimica Acta, vol. 147, pp. 294-309, 2014.
[4] I. V. Zenyuk, P. K. Das, and A. Z. Weber, "Understanding Impacts of Catalyst-Layer Thickness on Fuel-Cell Performance via Mathematical Modeling," Journal of The Electrochemical Society, vol. 163, no. 7, pp. F691-F703, 2016.
[5] L. Hao, K. Moriyama, W. Gu, and C.-Y. Wang, "Modeling and experimental validation of pt loading and electrode composition effects in pem fuel cells," Journal of The Electrochemical Society, vol. 162, no. 8, pp. F854-F867, 2015.
[6] M. Moore, A. Putz, and M. Secanell, "Investigation of the ORR using the double-trap intrinsic kinetic model," Journal of the Electrochemical Society, vol. 160, no. 6, pp. F670-F681, 2013.
[7] K. Neyerlin, H. A. Gasteiger, C. K. Mittelsteadt, J. Jorne, and W. Gu, "Effect of relative humidity on oxygen reduction kinetics in a PEMFC," Journal of The Electrochemical Society, vol. 152, no. 6, pp. A1073-A1080, 2005.
[8] Y. Liu, M. W. Murphy, D. R. Baker, W. Gu, C. Ji, J. Jorne, and H. A. Gasteiger, "Proton conduction and oxygen reduction kinetics in PEM fuel cell cathodes: effects of ionomer-to-carbon ratio and relative humidity," Journal of The Electrochemical Society, vol. 156, no. 8, pp. B970-B980, 2009.