Finite Element Approximation of a Three Dimensional Phase Field Model for Void Electromigration

  • Authors:

  • L'ubomír Baňas;Robert Nürnberg

  • Affiliations:

  • Department of Mathematics, Imperial College London, London, UK SW7 2AZ;Department of Mathematics, Imperial College London, London, UK SW7 2AZ

  • Venue:
  • Journal of Scientific Computing
  • Year:
  • 2008

Quantified Score

Hi-index 0.01

Visualization

Abstract

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 $$\gamma\frac{\partial u}{\partial t}-\nabla.(b(u)\nabla[w+\alpha\phi])=0,\qquad w=-\gamma\Delta u+\gamma^{-1}\Psi'(u),\qquad\nabla.(c(u)\nabla\phi)=0$$ subject to an initial condition u 0(驴)驴[驴1,1] on u and flux boundary conditions on all three equations. Here 驴驴驴0, 驴驴驴驴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.