A stable FE method for the space-time solution of the Cahn-Hilliard equation

Abstract

In its application to the modeling of a mineral separation process, we propose the numerical analysis of the Cahn-Hilliard equation by employing space-time discretizations of the automatic variationally stable finite element (AVS-FE) method. The AVS-FE method is a Petrov-Galerkin method which employs the concept of optimal discontinuous test functions of the discontinuous Petrov-Galerkin (DPG) method by Demkowicz and Gopalakrishnan. The trial space, however, consists of globally continuous Hilbert spaces such as H1(Ω) and H(div,Ω). Hence, the AVS-FE approximations employ classical C0 or Raviart-Thomas FE basis functions. The optimal test functions guarantee the numerical stability of the AVS-FE method and lead to discrete systems that are symmetric and positive definite. Hence, the AVS-FE method can solve the Cahn-Hilliard equation in both space and time without a restrictive CFL condition to dictate the space-time element size. We present multiple numerical verifications of both stationary and transient problems. The verifications show optimal rates of convergence in L2(Ω) and H1(Ω) norms. Results for mesh adaptive refinements in both space and time using a built-in error estimator of the AVS-FE method are also presented.