Abstract
An implicit-explicit (IMEX) method is combined with some so-called Richardson extrapolation (RiEx) methods for the numerical solution of reaction-diffusion equations with pure Neumann boundary conditions. The results are applied to a model for determining the overpotential in a Proton Exchange Membrane (PEM) fuel cell.