Numerical simulations of the energy-stable scheme for Swift-Hohenberg equation

A collocation Fourier scheme for Swift-Hohenberg equation based on the convex splitting idea is implemented. To ensure an efficient numerical computation, we propose a general framework with linear iteration algorithm to solve the non-linear coupled equations which arise with the semi-implicit scheme. Following the contraction mapping theorem, we present a detailed convergence analysis for the linear iteration algorithm. Various numerical simulations, including verification of accuracy, dissipative property of discrete energy and pattern formation, are presented to demonstrate the efficiency and the robustness of proposed method.