Abstract. The method of mapping function was first proposed by Henrick et al. [J.
Comput. Phys. 207:542-547 (2005)] to adjust nonlinear weights in [0,1] for the
fifth order WENO scheme, and through which the requirement of convergence order
is satisfied and the performance of the scheme is improved. Different from
Henrick’s method, a concept of piecewise polynomial function is
proposed in this study and corresponding WENO schemes are obtained. The
advantage of the new method is that the function can have a gentle profile at
the location of the linear weight (or the mapped nonlinear weight can be close
to its linear counterpart), and therefore is favorable for the resolution
enhancement. Besides, the function also has the flexibility of quick convergence
to identity mapping near two endpoints of [0,1], which is favorable for improved
numerical stability. The fourth-, fifth- and sixth-order polynomial functions
are constructed correspondingly with different emphasis on aforementioned
flatness and convergence. Among them, the fifth-order version has the flattest
profile. To check the performance of the methods, the 1-D Shu-Osher problem, the
2-D Riemann problem and the double Mach reflection are tested with the
comparison of WENO-M, WENO-Z and WENO-NS. The proposed new methods show the best
resolution for describing shear-layer instability of the Riemann problem, and
they also indicate high resolution in computations of double Mach reflection,
where only these proposed schemes successfully resolved the vortex-pairing
phenomenon. Other investigations have shown that the single polynomial mapping
function has no advantage over the proposed piecewise one, and it is of no
evident benefit to use the proposed method for the symmetric fifth-order WENO.
Overall, the fifth-order piecewise polynomial and corresponding WENO scheme are
suggested for resolution improvement.