Unconditional stability and convergence analysis of novel regularization schemes for the Navier-Stokes equations
This work addresses computational efficiency and stability issues in fluid dynamics simulations for researchers and engineers, though it appears incremental as it builds on existing regularization and SAV methods.
The paper tackles the problem of solving the Navier-Stokes equations by constructing novel first- and second-order decoupled schemes based on penalty and sequential regularization methods, achieving unconditional stability and convergence without time step constraints, with proven optimal error estimates in two dimensions and local estimates in three dimensions.
In this paper, we construct novel first- and second-order decoupled schemes for the Navier-Stokes equations based on the penalty method and the sequential regularization method (SRM), respectively. These schemes do not require the boundary condition on the pressure and thus preserve the original velocity boundary conditions. By using the idea of the scalar auxiliary variable (SAV), the nonlinear terms of these schemes are treated explicitly, which improves computational efficiency while maintaining stability. It is important to note that we carefully reformulated the Navier-Stokes system to ensure convergence of the proposed scheme without any restriction on the time step. For the Penalty-SAV (P-SAV) schemes, at each time step it is only necessary to solve elliptic equations with constant coefficients. We prove the high-order stability (high-order regularities of the solution) of the schemes, and establish an unconditional (without time step constraints) global optimal error estimate in two dimensions as well as a local error estimate in three dimensions for the first-order scheme. Furthermore, to more accurately approximate the incompressibility constraint without introducing extra stiffness into the system, the sequential regularization-SAV (SR-SAV) schemes are developed, and their error estimates are provided. In addition, we compare our proposed scheme with the classic linearized projection scheme to demonstrate its accuracy and efficiency.