Local Discontinuous Galerkin Methods for the Functionalized Cahn-Hilliard Equation

In this paper, we develop a local discontinuous Galerkin (LDG) method for the sixth order nonlinear functionalized Cahn-Hilliard (FCH) equation. We address the accuracy and stability issues from simulating high order stiff equations in phase-field modeling. Within the LDG framework, various boundary conditions associated with the background physics can be naturally implemented. We prove the energy stability of the LDG method for the general nonlinear case. A semi-implicit time marching method is applied to remove the severe time step restriction () for explicit methods. The adaptive capability of the LDG method allows for capturing the interfacial layers and the complicated geometric structures of the solution with high resolution. To enhance the efficiency of the proposed approach, the multigrid (MG) method is used to solve the system of linear equations resulting from the semi-implicit temporal integration at each time step. We show numerically that the MG solver has mesh-independent convergence rates. Numerical simulation results for the FCH equation in two and three dimensions are provided to illustrate that the combination of the LDG method for spatial approximation, semi-implicit temporal integration with the MG solver provides a practical and efficient approach when solving this family of problems.