Energy stable and $L^2$ norm convergent BDF3 scheme for the Swift-Hohenberg equation

Abstract: A fully discrete implicit scheme is proposed for the Swift-Hohenberg model, combining the third-order backward differentiation formula (BDF3) for the time discretization and the second-order finite difference scheme for the space discretization. Applying the Brouwer fixed-point theorem and the positive definiteness of the convolution coefficients of BDF3, the presented numerical algorithm is proved to be uniquely solvable and unconditionally energy stable, further, the numerical solution is shown to be bounded in the maximum norm. The proposed scheme is rigorously proved to be convergent in $L^2$ norm by the discrete orthogonal convolution (DOC) kernel, which transfer the four-level-solution form into the three-level-gradient form for the approximation of the temporal derivative. Consequently, the error estimate for the numerical solution is established by utilization of the discrete Gronwall inequality. Numerical examples in 2D and 3D cases are provided to support the theoretical results.