Efficient High Order Semi-implicit Time Discretization and Local Discontinuous Galerkin Methods for Highly Nonlinear PDEs

In this paper, we develop a high order semi-implicit time discretization method for highly nonlinear PDEs, which consist of the surface diffusion and Willmore flow of graphs, the Cahn-Hilliard equation and the Allen-Cahn/Cahn-Hilliard system. These PDEs are high order in spatial derivatives, which motivates us to develop implicit or semi-implicit time marching methods to relax the severe time step restriction for stability of explicit methods. In addition, these PDEs are also highly nonlinear, fully implicit methods will incredibly increase the difficulty of implementation. In particular, we can not well separate the stiff and non-stiff components for these problems, which leads to traditional implicit-explicit methods nearly meaningless. In this paper, a high order semi-implicit time marching method and the local discontinuous Galerkin (LDG) spatial method are coupled together to achieve high order accuracy in both space and time, and to enhance the efficiency of the proposed approaches, the resulting linear or nonlinear algebraic systems are solved by multigrid solver. Specially, we develop a first order fully discrete LDG scheme for the Allen-Cahn/Cahn-Hilliard system and prove the unconditional energy stability. Numerical simulation results in one and two dimensions are presented to illustrate that the combination of the LDG method for Institut Camille Jordan, Université Claude Bernard Lyon I, 43 boulevard 11 novembre 1918, 69622 Villeurbanne cedex, France. Email: [email protected]. Institut de Mathématiques de Toulouse, Université Paul Sabatier Toulouse III, 118 route de Narbonne, 31062 Toulouse cedex, France. Email: [email protected]. School of Mathematical Sciences, University of Science and Technology of China, Hefei, Anhui 230026, P.R. China. Email: [email protected]. Research supported by NSFC grant No. 11371342, 11426236.