Energy stable interior penalty discontinuous Galerkin finite element method for Cahn-Hilliard equation

Abstract

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 CH equation. Numerical results confirm the theoretical convergence rates and the performance of the proposed approach.