Robust superconvergence analysis of physics-preserving RMAC scheme for the Stokes and Navier--Stokes equations on non-uniform grids at high Reynolds numbers

The velocity errors of the classical marker and cell (MAC) scheme are dependent on the pressure approximation errors, which is non-pressure-robust and will cause the accuracy of the velocity approximation to deteriorate when the pressure approximation is poor. In this paper, we first propose the reconstructed MAC scheme (RMAC) based on the finite volume method to obtain the pressure-robustness for the time-dependent Stokes equations and then construct the $\mu$-robust and physics-preserving RMAC scheme on non-uniform grids for the Navier--Stokes equations, where $\mu$-robustness means that the velocity errors do not blow up for small viscosity $\mu$ when the true velocity is sufficiently smooth. Compared with the original MAC scheme, which was analyzed in [SIAM J. Numer. Anal. 55 (2017): 1135-1158], the RMAC scheme is different only on the right-hand side for Stokes equations. It can also be proved that the constructed scheme satisfies the local mass conservation law, the discrete unconditional energy dissipation law, the momentum conservation, and the angular momentum conservation for the Stokes and Navier--Stokes equations. Furthermore, by constructing the new auxiliary function depending on the velocity and using the high-order consistency analysis, we can obtain the pressure-robust and $\mu$-robust error estimates for the velocity and derive the second-order superconvergence for the velocity and pressure in the discrete $l^{\infty}(l^2)$ norm on non-uniform grids and the discrete $l^{\infty}(l^{\infty})$ norm on uniform grids. Finally, numerical experiments using the constructed schemes are demonstrated to show the robustness for our constructed schemes.