Stability and error analysis of a new class of higher-order consistent splitting schemes for the Navier-Stokes equations

Abstract

A new class of fully decoupled consistent splitting schemes for the Navier-Stokes equations are constructed and analyzed in this paper. The schemes are based on the Taylor expansion at $t^{n+\beta }$ with $\beta \ge 1$ being a free parameter. It is shown that by choosing {\color{black} $\beta =3,6,9$} respectively for the second-, third- and fourth-order schemes, their numerical solutions are uniformed bounded in a strong norm, and admit optimal global-in-time convergence rates in both 2D and 3D. {\color{black}These } results are the first stability and convergence results for any fully decoupled, higher than second-order schemes for the Navier-Stokes equations. Numerical results are provided to show that the third- and fourth-order schemes based on the usual BDF (i.e. $\beta =1$) are not unconditionally stable while the new third- and fourth-order schemes with suitable $\beta$ are unconditionally stable and lead to expected convergence rates.

arXiv commentThis article was accepted for publication in Mathematics of Computation on June 21, 2025