Positivity-preserving high order discontinuous Galerkin schemes for compressible Navier-Stokes equations

Xiangxiong Zhang, Purdue University
April 27th, 2016 at 3:30PM–4:30PM in 939 Evans Hall [Map]

For gas dynamics equations such as compressible Euler and Navier-Stokes equations, preserving the positivity of density and pressure without losing conservation is crucial to stabilize the numerical computation. The L1-stability of mass and energy can be achieved by enforcing the positivity of density and pressure during the time evolution. However, high order schemes such as DG methods do not preserve the positivity. It is difficult to enforce the positivity without destroying the high order accuracy and the local conservation in an efficient manner for time-dependent gas dynamics equations. For compressible Euler equations, a weak positivity property holds for any high order finite volume type schemes including DG methods, which was used to design a simple positivity-preserving limiter for high order DG schemes in Zhang and Shu, JCP 2010. Generalizations to compressible Navier-Stokes equations are however nontrivial. We show that weak positivity property still holds for DG method solving compressible Navier-Stokes equations if a proper penalty term is added in the scheme. This allows us to obtain the first high order positivity-preserving schemes for compressible Navier-Stokes equations.