Hierarchical reconstruction with up to second degree remainder for solving nonlinear conservation laws^*

The hierarchical reconstruction (HR) (Liu et al 2007 SIAM J. Numer. Anal. 45 2442–67) can effectively reduce spurious oscillations without local characteristic decomposition for numerical capturing of discontinuous solutions. However, there are still small remaining overshoots/undershoots in the vicinity of discontinuities. HR with partial neighbouring cells (Xu et al 2009 J. Comput. Phys. 228 2194–212) essentially overcomes this drawback for the third order case, and in the mean time further improves the resolution of the numerical solution. Extending the technique to higher order cases we observe the returning of overshoots/undershoots. In this paper, we introduce a new technique to work with HR on partial neighbouring cells, which lowers the order of the remainder of the polynomial in the current cell while maintaining the theoretical order of accuracy, essentially eliminates overshoots/undershoots for the fourth and fifth order cases and reduces the numerical cost.