Algebraic construction and numerical behaviour of a new s-consistent difference scheme for the 2D Navier-Stokes equations

In this paper we consider a regular grid with equal spatial spacings and construct a new finite difference approximation (difference scheme) for the system of two-dimensional Navier-Stokes equations describing the unsteady motion of an incompressible viscous liquid of constant viscosity. In so doing, we use earlier constructed discretization of the system of three equations: the continuity equation and the proper Navier-Stokes equations. Then, we compute the canonical Gröbner basis form for the obtained discrete system. It gives one more difference equation which is equivalent to the pressure Poisson equation modulo difference ideal generated by the Navier-Stokes equations, and thereby comprises a new finite difference approximation (scheme). We show that the new scheme is strongly consistent. Besides, our computational experiments demonstrate much better numerical behaviour of the new scheme in comparison with the other strongly consistent schemes we constructed earlier and also with the scheme which is not strongly consistent.