A mixed discontinuous Galerkin, convex splitting scheme for a modified Cahn-Hilliard equation and an efficient nonlinear multigrid solver

Abstract

In this paper we devise and analyze a mixed discontinuous Galerkin finite element method for a modified Cahn-Hilliard equation that models phase separation in diblock copolymer melts. The time discretization is based on a convex splitting of the energy of the equation. We prove that our scheme is unconditionally energy stable with respect to a spatially discrete analogue of the continuous free energy of the system, unconditionally uniquely solvable, and convergent in the natural energy norm with optimal rates. We describe an efficient nonlinear multigrid solver for advancing our semi-implicit scheme in time and conclude the paper with numerical tests confirming the predicted convergence rates and suggesting the near-optimal complexity of the solver.