This paper introduces a C0 weak Galerkin finite element method for a linear Cahn–Hilliard–Cook equation. The highlights of the proposed method are that the complexity of constructing the C1 finite element space for fourth order problem is avoided and the number of degree of freedom is apparently reduced compared to the fully discontinuous weak Galerkin finite element method. With the redefined discrete weak Laplace operator and the classical C0 Lagrange elements, the L2 optimal error estimates in spatial variable are obtained. In time, the classical Euler scheme is then used to do the numerical simulation. Finally, numerical experiments are presented to demonstrate the efficiency of the proposed numerical method.
Read full abstract