To facilitate the numerical solution of this problem, Bourdin et al 26 introduced a phase-field approximation of Eq.(1) that takes the form: