A second-order, weakly energy-stable pseudo-spectral scheme for the Cahn-Hilliard equation and its solution by the homogeneous linear iteration method (Q2399231)

The paper presents an unconditionally energy stable second-order numerical scheme for the Cahn-Hilliard equation with Fourier pseudo-spectral approximation in space. The temporal discretization follows the second-order convex splitting. In order to carry out Fourier pseudo-spectral analysis, a second-order artificial diffusion term, and a Douglas-Dupont-type regularization are introduced. A homogeneous linear iteration algorithm is developed to solve the nonlinear equation. The scheme is numerically stable when it is implemented with the help of the Fft Fourier transform (FFT) in a pseudo-spectral setting. The numerical simulation experiments show that the second-order scheme, combined with the linear iteration algorithm is able to produce accurate long time numerical results. Several theorems and lemmas are stated and proved.