Numerical algorithms based on Galerkin methods for the modeling of reactive interfaces in photoelectrochemical solar cells