Stable Computing with an Enhanced Physics Based Scheme for the 3D Navier-Stokes Equations

Abstract

We study extensions of the energy and helicity preserving scheme for the 3D Navier-Stokes equations, developed in [23], to a more general class of problems. The scheme is studied together with stabilizations of grad-div type in order to mitigate the effect of the Bernoulli pressure error on the velocity error. We prove stability, convergence, discuss conservation properties, and present numerical experiments that demonstrate the advantages of the scheme.