We propose a degenerate Allen-Cahn/Cahn-Hilliard system coupled to a quasi-static diffusion equation to model the motion of intergranular voids. The system can be viewed as a phase field system with an interfacial parameter γ . In the limit γ → 0, the phase field system models the evolution of voids by surface diffusion and electromigration in an electrically conducting solid with a grain boundary. We introduce a finite element approximation for the proposed system, show stability bounds, prove convergence, and hence existence of a weak solution to this non-linear degenerate parabolic system in two space dimensions. An iterative scheme for solving the resulting non-linear discrete system at each time level is introduced and analysed, and some numerical experiments are presented. In the Appendix we discuss the sharp interface limit of the above degenerate system as the interfacial parameter γ tends to zero.