A finite volume/element approximation of a mean field model of superconducting vortices in one and two space dimensions is presented. The model incorporates positive and negative vorticity with time dependent applied magnetic field, flux pinning and boundary nucleation. Stability bounds on the solutions of these approximations are proved. Some computed results are presented.