Search

The open circuit voltage (OCV) in a proton exchange membrane fuel cell (PEMFC) is typically recorded as being approximately 300 mV lower than the equilibrium voltage computed by the Nernst equation. While a number of causes have been proposed, the voltage drop is generally attributed to the oxidation of crossover hydrogen in the cathode. A single phase, through-the-channel model is presented that includes hydrogen transport across the membrane, an empirical model for the hydrogen oxidation reaction (HOR) fit to experimental data obtained at high potentials and a multi-step kinetic model to describe the oxygen reduction reaction (ORR). Model predictions were compared to experimentally obtained OCVs and the results show that the model is capable of capturing the experimentally observed changes in OCV with platinum loading, as well as fuel cell performance; and that, at low Pt loadings, small quantities of unreacted hydrogen leave the cathode because the HOR is kinetically limited by oxide blocking and anion adsorption. A parametric study is used to show that a minimum OCV is achieved at ultra-low loadings. Results also show that only a multi-step ORR model can simultaneously capture polarization data and the OCV.

The open-circuit voltage (OCV) is the measured potential when no current is drawn from a cell. Theoretically this potential is given by the equilibrium potential, which is determined from thermodynamics based on the operating conditions. In practice, the actual potential at zero current is measured to be approximately 300 mV lower than the theoretical value. There are a number of possible causes for the lower OCV [1,2], however the most commonly accepted is that hydrogen crossing from the anode to the cathode will react in the presence of the platinum catalyst, thereby developing a mixed potential [2,3,4]. The OCV will therefore depend on the relative kinetics of the hydrogen oxidation reaction (HOR) and oxygen reduction reaction (ORR) in the cathode such that when the ORR is equal to the HOR in the cathode, the OCV is reached, see Figure 1A. The overpotential required to achieve the required ORR will determine the potential at which OCV occurs and the overpotential required will be impacted by the platinum loading in the cathode. This effect is shown experimentally, where the OCV is seen to drop significantly at loadings less than 0.05mg/cm2 , as shown in Figure 1B. The two-dimensional, macro-homogeneous model in OpenFCST [5] is extended to allow multiple reactions in an electrode and used to study the impact of platinum loading on OCV. The model accounts for the transport of oxygen, hydrogen, water vapour and solved water, as well as the potentials of both the solid and electrolyte phases. Both the ORR and HOR are considered in the cathode, and as the hydrogen crossover rate is only of the order of 1 mA/cm2, the accuracy of the kinetic models used for the ORR and HOR are studied. Furthermore, the impact of the HOR kinetics is examined by considering recent experimental studies [7,8] that have shown the HOR is severely limited at high overpotentials by anion adsorption and oxide formation. This work shows that the double-trap kinetic model [6] for the ORR is able to match the polarization data in the kinetic region as well as correctly predict the OCVs for platinum loadings of 0.014-0.113 mg/cm2 . The commonly used Tafel equation, in contrast, cannot simultaneously fit both the kinetic region of the polarisation curve and the OCVs. Using a HOR model that reproduces the kinetic data showing limited activity at high overpotentials does not accurately capture the OCV and predicts that hydrogen would escape through the cathode. The impact of microscale transport is studied using the ionomer covered carbon particle (ICCP) model. A parametric study on the parameters controlling the dissolution of hydrogen into the ionomer thin film and its diffusion to the reaction site indicate that hydrogen transport is not limiting at these scales for the cell structural parameters used in this study, but may be important at higher ionomer film thicknesses. Changing the operating oxygen partial pressure shows that the double trap model in ref. [6] does not accurately capture the change in reaction rate with the change in oxygen concentration. This work demonstrates that modeling the OCV is helpful in fitting the kinetic region of the polarisation curve, as an accurate fit here and at OCV means that the kinetic model is accurate over several orders of magnitude in current and approximately 100mV of potential. References: [1] J. J. Giner-Sanz et al., Int. J. Hydrogen Energ., 39 (2014) 13206-13216 [2] S. A. Vilekar et al., J. Power Sources, 195 (2010) 2241–2247 [3] V. M. Ehlinger et al., J. Electrochem. Soc., 166, (7) F3255-F3267 (2019) [4] A. A. Kulikovsky, Electrochimica Acta, 247 (2017) 730–735 [5] M. Secanell et al., Electrochem. Soc. Trans., 64 (3) 655-680 (2014) [6] M. Moore et al., J. Electrochem. Soc., 160 (6) F670-F681 (2013) [7] J. Durst et al., J. Electrochem. Soc., 162 (1) F190-F203 (2015) [8] C. M. Zalitis et al., Phys. Chem. Chem. Phys., 2013, 15, 4329 Figure 1