Abstract
We investigate the error between any discrete solution of the implicit marker-and-cell (MAC) numerical scheme for compressible Navier–Stokes equations in the low Mach number regime and an exact strong solution of the incompressible Navier–Stokes equations. The main tool is the relative energy method suggested on the continuous level in Feireisl et al. (2012, Relative entropies, suitable weak solutions, and weak–strong uniqueness for the compressible Navier–Stokes system. J. Math. Fluid Mech., 14, 717–730). Our approach highlights the fact that numerical and mathematical analyses are not two separate fields of mathematics. The result is achieved essentially by exploiting in detail the synergy of analytical and numerical methods. We get an unconditional error estimate in terms of explicitly determined positive powers of the space–time discretization parameters and Mach number in the case of well-prepared initial data and in terms of the boundedness of the error if the initial data are ill prepared. The multiplicative constant in the error estimate depends on a suitable norm of the strong solution but it is independent of the numerical solution itself (and of course, on the discretization parameters and the Mach number). This is the first proof that the MAC scheme is unconditionally and uniformly asymptotically stable in the low Mach number regime.
H2-stability of the first order Galerkin method for the Boussinesq equations with smooth and non-smooth initial data