Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy

In this paper we design a class of numerical schemes that are higher-order extensions of the weighted essentially non-oscillatory (WENO) schemes of G.-S. Jiang and C.-W. Shu (1996) and X.-D. Liu, S. Osher, and T. Chan (1994). Used by themselves, the schemes may not always be monotonicity preserving but coupled with the monotonicity preserving bounds of A. Suresh and H. T. Huynh (1997) they perform very well. The resulting monotonicity preserving weighted essentially non-oscillatory (MPWENO) schemes have high phase accuracy and high order of accuracy. The higher-order members of this family are almost spectrally accurate for smooth problems. Nevertheless, they, have robust shock capturing ability. The schemes are stable under normal CFL numbers. They are also efficient and do not have a computational complexity that is substantially greater than that of the lower-order members of this same family of schemes. The higher accuracy that these schemes offer coupled with their relatively low computational complexity makes them viable competitors to lower-order schemes, such as the older total variation diminishing schemes, for problems containing both discontinuities and rich smooth region structure. We describe the MPWENO schemes here as well as show their ability to reach their designed accuracies for smooth flow. We also examine the role of steepening algorithms such as the artificial compression method in the design of very high order schemes. Several test problems in one and two dimensions are presented. For multidimensional problems where the flow is not aligned with any of the grid directions it is shown that the present schemes have a substantial advantage over lower-order schemes. It is argued that the methods designed here have great utility for direct numerical simulations and large eddy simulations of compressible turbulence. The methodology developed here is applicable to other hyperbolic systems, which is demonstrated by showing that the MPWENO schemes also work very well on magnetohydrodynamical test problems.