Numerical simulations of proton exchange water electrolysis for hydrogen production were performed for the purpose of examining the phenomena occurring within the proton exchange membranes (PEM) water splitting cell. A two-dimensional steady-state isothermal model of the cell has been developed. Finite element method was used to solve the multicomponent transport model coupled with flow in porous medium, charge balance and electrochemical kinetics. The Maxwell-Stefan equation is applied for the multi-component diffusion and convection in water distribution electrodes. The Butler-Volmer kinetic equation is used to obtain the local current density distribution at the catalyst reactive boundaries. Darcy’s law was applied for the flow of species in the porous electrodes. Parametric studies are performed based on appropriate mass balances, transport, and electrochemical kinetics applied to the electrolysis cell. There are significant current spikes present at the corners of the current collector. The current density varies significantly in the cell, being highest at the corners of the current collector. As the water on the anode side flows from the inlet to the outlet, the mass fraction of oxygen increases. This is the effect of oxygen concentration due to the effect oxidation of water. On the cathode side, as the mass fraction of water decreases there is little variation in the hydrogen mass fraction content due to the effect of hydrogen reduction.

This content is only available via PDF.
You do not currently have access to this content.