ANALYSIS OF MIXED INTERIOR PENALTY DISCONTINUOUS GALERKIN METHODS FOR THE CAHN-HILLIARD EQUATION AND THE HELE-SHAW FLOW

This paper proposes and analyzes two fully discrete mixed interior penalty discontinuous Galerkin (DG) methods for the fourth order nonlinear Cahn-Hilliard equation. Both methods use the backward Euler method for time discretization and interior penalty DG methods for spatial discretization. They differ from each other on how the nonlinear term is treated; one of them is based on fully implicit time-stepping and the other uses energy-splitting time-stepping. The primary goal of the paper is to prove the convergence of the numerical interfaces of the DG methods to the interface of the Hele-Shaw flow. This is achieved by establishing error estimates that depend on epsilon(-1) only in some low polynomial orders, instead of exponential orders. Similar to [X. Feng and A. Prohl, Numer. Math., 74 (2004), pp. 47-84], the crux is to prove a discrete spectrum estimate in the discontinuous Galerkin finite element space. However, the validity of such a result is not obvious because the DG space is not a subspace of the (energy) space H-1(Omega) and it is larger than the finite element space. This difficulty is overcome by a delicate perturbation argument which relies on the discrete spectrum estimate in the finite element space proved by Feng and Prohl. Numerical experiment results are also presented to gauge the theoretical results and the performance of the proposed fully discrete mixed DG methods.