Numerical investigation of a viscous regularization of the Euler equations by entropy viscosity

Abstract

The Navier–Stokes viscous fluxes are a well-known viscous regularization of the Euler equations. However, since these fluxes do not add any viscosity to the mass equation, the positivity of density is violated. This paper investigates a new class of viscous regularization of the Euler equations, which was recently proposed by Guermond & Popov (2014). In contrast to the Navier–Stokes fluxes, the new regularization adds a viscous term to the mass equation. Since non-physical viscous terms are used, it is important to show that the exact solution’s properties, such as the location of shocks, contact and rarefaction waves are not violated. The present study concerns a careful numerical investigation of the new viscous regularization in a number of well-known 1D and 2D benchmark problems. Also, a direct numerical comparison with respect to the physical Navier–Stokes regularization is shown. The numerical tests show that the entropy viscosity method can achieve high order accuracy for any polynomial degrees. Detailed algorithms for the implementation of a slip wall boundary condition are presented in a weak and a strong form.