Pure-positivity-preserving methods with an optimal sufficient CFL number for fifth-order MR-WENO schemes on structured meshes

Abstract

In this paper, one-dimensional and two-dimensional pure-positivity-preserving (PPP) methods are proposed for fifth-order finite volume multi-resolution WENO (MR-WENO) schemes to solve some extreme problems on structured meshes. The MR-WENO spatial reconstruction procedures only require one five-cell, one three-cell, and one one-cell stencils for achieving uniform fifth-order accuracy in smooth regions and keeping essentially non-oscillatory property in non-smooth regions in one dimension. One redefines five new cell averages vectors after performing such spatial reconstructions and design one quartic polynomials vector and three quadratic polynomials vectors based on them. After that, a new detective process is used to examine the positivity of density and pressure of three quadratic polynomials vectors inside the whole target cell. If the negativity happens, a new compression limiter is carried out to enable the positivity of density and pressure of three quadratic polynomials vectors over the whole target cell and the positivity of density and pressure of one quartic polynomials vector at the midpoint of the target cell. It is a new way to design the positivity-preserving methods to keep fifth-order accuracy and the positivity over the target cell instead of only at some discrete Gauss-Lobatto quadrature points, since the precise minimum values of the density and pressure are now available. Then a theoretically proof is given to increase the optimal sufficient CFL number from 1/12 to 1/6 for the fifth-order WENO schemes. This methodology can be expanded to multi-dimensions easily. Unlike some classical positivity-preserving methods, the PPP methods could apply a special four-point Gauss-Lobatto quadrature formula or any other quadrature formulas on condition that their numerical precision is no smaller than four. Since the optimal CFL number of 1/6 is a sufficient but not necessary condition, the novelty PPP methods for fifth-order finite volume MR-WENO schemes with a larger practical CFL number of 0.6 are also available and robust enough when simulating some extreme problems without timely halving its value.