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.