Efficient implementation of adaptive order reconstructions

28 Jan 2020  ·  Semplice Matteo, Visconti Giuseppe ·

Including polynomials with small degree and stencil when designing very high order reconstructions is surely beneficial for their non oscillatory properties, but may bring loss of accuracy on smooth data unless special care is exerted. In this paper we address this issue with a new Central WENOZ (CWENOZ) approach, in which the reconstruction polynomial is computed from a single set of non linear weights, but the linear weights of the polynomials with very low degree (compared to the final desired accuracy) are infinitesimal with respect to the grid size. After proving general results that guide the choice of the CWENOZ parameters, we study a concrete example of a reconstruction that blends polynomials of degree six, four and two, mimicking already published Adaptive Order WENO reconstructions [Balsara, Garain and Shu, J.Comput.Phys., 2016], [Arbogast, Huang and Zhao, SIAM J.Numer.Anal., 2018]. The novel reconstruction yields similar accuracy and oscillations with respect to the previous ones, but saves up to 20% computational time since it does not rely on a hierarchic approach and thus does not compute multiple sets of nonlinear weights in each cell.

PDF Abstract
No code implementations yet. Submit your code now

Categories


Numerical Analysis Numerical Analysis