We consider a finite element approximation of a phase field model for the evolution of voids by surface diffusion in an electrically conducting solid. The phase field equations are given by the nonlinear degenerate parabolic system, 1] on u and flux boundary conditions on all three equations. Here γ ∈ R >0 , α ∈ R ≥0 , is a non-smooth double well potential, and c(u) := 1 + u, b(u) := 1 − u 2 are degenerate coefficients. On extending existing results for the simplified two dimensional phase field model, we show stability bounds for our approximation and prove convergence, and hence existence of a solution to this nonlinear degenerate parabolic system in three space dimensions. Furthermore, a new iterative scheme for solving the resulting nonlinear discrete system is introduced and some numerical experiments are presented.Keywords Void electromigration · Surface diffusion · Phase field model · Degenerate Cahn-Hilliard equation · Fourth order degenerate parabolic system · Finite elements · Convergence analysis · Multigrid methods
IntroductionIn the recent paper [9], abbreviated to BNS throughout this paper, the authors proposed and analysed a fully practical finite element approximation for a phase field model describing