Numerical analysis of coupled electron and mass transport in the gas diffusion layer of a PEMFC

Numerical investigation and parametric simulations of the coupled electrical conduction and mass diffusion in the GDL of a PEMFC are reported. The coupled transport is solved using a 2D, "under-the-rib" computational domain for the cathode GDL. The current density on the GDL/catalyst layer interface, which constitutes one of the boundary conditions for the GDL domain and reflects the activation overpotential in the catalyst layer and the ohmic loss in the membrane, is solved iteratively using a novel numerical algorithm. A parametric study is performed to investigate the effects on current density distribution of various operating conditions such as oxygen concentration and membrane resistance, and of design factors such as GDL geometry, anisotropic transport properties, and deformation under the land area. The results show that the current density distribution under the land area can be dominated by either electron transport or mass transport, depending on the operating regime.