In addition to theoretical challenges, a more practical challenge arises from solving the (GM)PNP system of partial difference equations numerically at operating potentials. The Galerkin finite element method is found to be unstable for typical values of applied potentials unless a very fine time and space discretization is used to resolve the EDL region. The SUPG method was used in this study to stabilize the PNP system as described in section 2.3. However, there is still a need to develop stabilization techniques suited for the modified version of the Nernst-Planck equation with a non-linear volume exclusion term added. The development of better stabilization techniques will go a long way in making the application of the GMPNP system more feasible to a wider array of operating conditions for electrocataytic systems of relevance.