Energy Stable Interior Penalty Discontinuous Galerkin Finite Element Method for Cahn–Hilliard Equation

An energy stable conservative method is developed for the Cahn–Hilliard (CH) equation with the degenerate mobility. The CH equation is discretized in space with the mass conserving symmetric interior penalty discontinuous Galerkin (SIPG) method. The resulting semi-discrete nonlinear system of ordinary differential equations are solved in time by the unconditionally energy stable average vector field (AVF) method. We prove that the AVF method preserves the energy decreasing property of the fully discretized CH equation. Numerical results for the quartic double-well and the logarithmic potential functions with constant and degenerate mobility confirm the theoretical convergence rates, accuracy and the performance of the proposed approach.