Abstract
Polymer electrolyte membrane (PEM) electrolysis has a promising future for large-scale hydrogen production. As PEM electrolysis technology develops, larger operating current densities are required. In order to increase current density, more water should be available at the reaction sites. Moreover, the removal rate of oxygen and hydrogen needs to be effectively improved. This, in turn, necessitates a better understanding of the main mass transport and electrochemical processes. On the anode side, mass transport is particularly crucial because water must be supplied to the catalyst layer (CL) while, at the same time, oxygen bubbles must be eliminated in a parallel flow from the reaction sites into the flow channels. Hence, simulating the two-phase bubbly flow across the cell thickness is necessary to predict PEM electrolysis performance more accurately as a function of the operating current density. This study provides a systematic understanding of how morphological and geometrical features contribute to the polarization curve and performance characteristics of a PEM electrolysis cell. Hence, a multi-phase PEM electrolysis model has been implemented using MATLAB R2022a. Polarization curves have been calibrated against experimental data and then assessed to provide a fundamental understanding of the relationship between the two-phase flow and cell performance.