Analysis of interior penalty discontinuous Galerkin methods for the Allen-Cahn equation and the mean curvature flow

This paper develops and analyzes two fully discrete interior penalty discontinuous Galerkin (IP-DG) methods for the Allen-Cahn equation, which is a nonlinear singular perturbation of the heat equation and originally arises from phase transition of binary alloys in materials science, and its sharp interface limit (the mean curvature flow) as the perturbation parameter tends to zero. Both fully implicit and energy-splitting time-stepping schemes are proposed. The primary goal of the paper is to derive sharp error bounds which depend on the reciprocal of the perturbation parameter $\epsilon$ (also called "interaction length") only in some lower polynomial order, instead of exponential order, for the proposed IP-DG methods. The derivation is based on a refinement of the nonstandard error analysis technique first introduced in [12]. The centerpiece of this new technique is to establish a spectrum estimate result in totally discontinuous DG finite element spaces with a help of a similar spectrum estimate result in the conforming finite element spaces which was established in [12]. As a nontrivial application of the sharp error estimates, they are used to establish convergence and the rates of convergence of the zero level sets of the fully discrete IP-DG solutions to the classical and generalized mean curvature flow. Numerical experiment results are also presented to gauge the theoretical results and the performance of the proposed fully discrete IP-DG methods.

[1]  Xinfu Chen,et al.  Spectrum for the allen-chan, chan-hillard, and phase-field equations for generic interfaces , 1994 .

[2]  J. Cahn,et al.  A microscopic theory for antiphase boundary motion and its application to antiphase domain coasening , 1979 .

[3]  Yun-Gang Chen,et al.  Uniqueness and existence of viscosity solutions of generalized mean curvature flow equations , 1989 .

[4]  Haijun Wu,et al.  A Posteriori Error Estimates and an Adaptive Finite Element Method for the Allen–Cahn Equation and the Mean Curvature Flow , 2005, J. Sci. Comput..

[5]  Béatrice Rivière,et al.  Discontinuous Galerkin methods for solving elliptic and parabolic equations - theory and implementation , 2008, Frontiers in applied mathematics.

[6]  G. McFadden Phase-Field Models of Solidification , 2002 .

[7]  Andreas Prohl,et al.  Numerical analysis of the Allen-Cahn equation and approximation for mean curvature flows , 2003, Numerische Mathematik.

[8]  Zhangxin Chen,et al.  Pointwise Error Estimates of Discontinuous Galerkin Methods with Penalty for Second-Order Elliptic Problems , 2004, SIAM J. Numer. Anal..

[9]  Sören Bartels,et al.  Robust A Priori and A Posteriori Error Analysis for the Approximation of Allen-Cahn and Ginzburg-Landau Equations Past Topological Changes , 2011, SIAM J. Numer. Anal..

[10]  D. J. Eyre,et al.  An Unconditionally Stable One-Step Scheme for Gradient Systems , 1997 .

[11]  G. Burton Sobolev Spaces , 2013 .

[12]  Ricardo H. Nochetto,et al.  Convergence Past Singularities for a Fully Discrete Approximation of Curvature-Drive Interfaces , 1997 .

[13]  P. Souganidis,et al.  Phase Transitions and Generalized Motion by Mean Curvature , 1992 .

[14]  Maurizio Paolini,et al.  Quasi-optimal error estimates for the mean curvature flow with a forcing term , 1995, Differential and Integral Equations.

[15]  ADAPTIVE DISCONTINUOUS GALERKIN APPROXIMATIONS OF SECOND-ORDER ELLIPTIC PROBLEMS , 2004 .

[16]  Yinnian He,et al.  Analysis of finite element approximations of a phase field model for two-phase fluids , 2006, Math. Comput..

[17]  B. G. Pachpatte Inequalities for finite difference equations , 2001 .

[18]  Tom Ilmanen,et al.  Convergence of the Allen-Cahn equation to Brakke's motion by mean curvature , 1993 .

[19]  Ricardo H. Nochetto,et al.  A posteriori error control for the Allen–Cahn problem: circumventing Gronwall's inequality , 2004 .