Efficient implementation of adaptive order reconstructions

by   Matteo Semplice, et al.

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 rely on a hierarchic approach and thus does not compute multiple sets of nonlinear weights in each cell.


page 1

page 2

page 3

page 4


WENO interpolations and reconstructions using data bounded polynomial approximation

This work characterizes the structure of third and forth order WENO weig...

Dissipation-based WENO stabilization of high-order finite element methods for scalar conservation laws

We present a new perspective on the use of weighted essentially nonoscil...

A Hermite WENO scheme with artificial linear weights for hyperbolic conservation laws

In this paper, a fifth-order Hermite weighted essentially non-oscillator...

A low-dissipation shock-capturing framework with flexible nonlinear dissipation control

In this work, a framework to construct arbitrarily high-order low-dissip...

Fast, reliable and unrestricted iterative computation of Gauss--Hermite and Gauss--Laguerre quadratures

Methods for the computation of classical Gaussian quadrature rules are d...

3D EIT Reconstructions from Electrode Data using Direct Inversion D-bar and Calderon Methods

The first numerical implementation of a D-bar method in 3D using electro...

Please sign up or login with your details

Forgot password? Click here to reset