Summary

In this paper we present an error analysis for a second order accurate numerical scheme for the 2-D and 3-D Cahn-Hilliard (CH) equation, with an improved convergence constant. The unique solvability, unconditional energy stability, and a uniform-in-time H-2 stability of this numerical scheme have already been established. However, a standard error estimate gives a convergence constant in an order of exp(CT epsilon(-m0)), with m(0) a positive integer and the interface width parameter epsilon being small, which comes from the application of discrete Gronwall inequality. To overcome this well-known difficulty, we apply a spectrum estimate for the linearized Cahn-Hilliard operator (Alikakos and Fusco, 1993; Chen, 1994; Feng and Prohl, 2004), perform a detailed numerical analysis, and get an improved estimate, in which the convergence constant depends on 1/epsilon only in a polynomial order, instead of the exponential order.

  • Institution
    苏州大学

Full-Text