A positivity preserving strategy for entropy stable discontinuous Galerkin discretizations of the compressible Euler and Navier-Stokes equations

Yimin Lin,Jesse Chan,Ignacio Tomas
DOI: https://doi.org/10.1016/j.jcp.2022.111850
2022-09-16
Abstract:High-order entropy-stable discontinuous Galerkin methods for the compressible Euler and Navier-Stokes equations require the positivity of thermodynamic quantities in order to guarantee their well-posedness. In this work, we introduce a positivity limiting strategy for entropy-stable discontinuous Galerkin discretizations constructed by blending high order solutions with a low order positivity-preserving discretization. The proposed low order discretization is semi-discretely entropy stable, and the proposed limiting strategy is positivity preserving for the compressible Euler and Navier-Stokes equations. Numerical experiments confirm the high order accuracy and robustness of the proposed strategy.
Numerical Analysis
What problem does this paper attempt to address?
This paper attempts to solve the problem of maintaining the positivity of physical quantities (such as density and pressure) in the high - order entropy - stable discontinuous Galerkin (DG) method. Specifically: 1. **Stability issues of high - order DG methods**: Although high - order DG methods are efficient and accurate in dealing with convection - dominated problems, they usually encounter stability problems when the solution characteristics are not fully resolved. Entropy - stable DG discretization can improve stability by satisfying semi - discrete entropy balance while maintaining high - order accuracy. However, such schemes do not guarantee the positivity of thermodynamic quantities (such as density, pressure) during the simulation process. 2. **Importance of positivity preservation**: To ensure the well - posedness of the numerical method, thermodynamic quantities such as density and pressure must be maintained as positive values throughout the simulation process. If these quantities become negative, it may lead to unphysical results or numerical instability. 3. **Proposed method**: The author introduces a new positivity - limiting strategy, which is constructed by mixing the high - order solution with a low - order positivity - preserving discretization scheme. The proposed low - order discretization scheme is semi - discrete entropy - stable, and the proposed limiting strategy can ensure the positivity of the compressible Euler equations and Navier - Stokes equations. 4. **Technical details**: - The low - order discretization uses sparse low - order SBP operators, combined with a graph viscosity term to ensure positivity. - Through appropriate CFL conditions and forward Euler time - stepping, it can be proved that the low - order scheme can maintain the positivity of density and pressure. - For high - order DG methods, by introducing an element - level limiter based on the Zhang - Shu type, the high - order and low - order solutions can be mixed on each element, thereby achieving positivity preservation while maintaining high - order accuracy. 5. **Verification and application**: Numerical experiments show that this strategy can significantly improve the robustness of the numerical solution while maintaining high - order accuracy and performs well in practical applications. In summary, this paper aims to solve the problem of maintaining the positivity of thermodynamic quantities in high - order entropy - stable DG methods and proposes a limiting strategy that combines high - order and low - order discretization schemes to ensure the physical rationality and stability of the numerical solution.