Improvement of the WENO-NIP Scheme for Hyperbolic Conservation Laws

Abstract

The WENO-NIP scheme was obtained by developing a class of L1-norm smoothness indicators based on Newton interpolation polynomial. It recovers the optimal convergence order in smooth regions regardless of critical points and achieves better resolution than the classical WENO-JS scheme. However, the WENO-NIP scheme produces severe spurious oscillations when solving 1D linear advection problems with discontinuities at long output times, and it is also very oscillatory near discontinuities for 1D Riemann problems. In this paper, we find that the spectral property of WENO-NIP exhibits the negative dissipation characteristic, and this is the reason why WENO-NIP is unstable near discontinuities. Using this knowledge, we develop a way of improving the WENO-NIP scheme by introducing an additional term to eliminate the negative dissipation interval. The proposed scheme, denoted as WENO-NIP+, maintains the same convergence property, as well as the same low-dissipation property, as the corresponding WENO-NIP scheme. Numerical examples confirm that the proposed scheme is much more stable near discontinuities for 1D linear advection problems with large output times and 1D Riemann problems than the WENO-NIP scheme. Furthermore, the new scheme is far less dissipative in the region with high-frequency waves. In addition, the improved WENO-NIP+ scheme can remove or at least greatly decrease the post-shock oscillations that are commonly produced by the WENO-NIP scheme when simulating 2D Euler equations with strong shocks.