A RBF-WENO finite volume method for hyperbolic conservation laws with the monotone polynomial interpolation method

Abstract Essentially non-oscillatory (ENO) and weighted ENO (WENO) methods are efficient high order numerical methods for solving hyperbolic conservation laws designed to reduce the Gibbs oscillations. The original ENO and WENO methods are based on the polynomial interpolation and the overall convergence rate is uniquely determined by the total number of interpolation points involved for the approximation. In this paper, we propose non-polynomial ENO and WENO finite volume methods in order to enhance the local accuracy and convergence. The infinitely smooth radial basis functions (RBFs) are adopted as a non-polynomial interpolation basis. Particularly we use the multi-quadratic and Gaussian RBFs. The non-polynomial interpolation such as the RBF interpolation offers the flexibility to control the local error by optimizing the free parameter. Then we show that the non-polynomial interpolation can be represented as a perturbation of the polynomial interpolation. To guarantee the essentially non-oscillatory property, the monotone polynomial interpolation method is introduced as a switching method to the polynomial reconstruction adaptively near the non-smooth area. The numerical results show that the developed non-polynomial ENO and WENO methods with the monotone polynomial interpolation method enhance the local accuracy and give sharper solution profile than the ENO/WENO methods based on the polynomial interpolation.

[1]  Wai-Sun Don,et al.  An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws , 2008, J. Comput. Phys..

[2]  Antonio Marquina,et al.  Power ENO methods: a fifth-order accurate weighted power ENO method , 2004 .

[3]  Chi-Wang Shu,et al.  High Order Weighted Essentially Nonoscillatory Schemes for Convection Dominated Problems , 2009, SIAM Rev..

[4]  Armin Iske,et al.  Adaptive ADER Methods Using Kernel-Based Polyharmonic Spline WENO Reconstruction , 2010, SIAM J. Sci. Comput..

[5]  J. M. Powers,et al.  Mapped weighted essentially non-oscillatory schemes: Achieving optimal order near critical points , 2005 .

[6]  Jingyang Guo,et al.  Radial Basis Function ENO and WENO Finite Difference Methods Based on the Optimization of Shape Parameters , 2016, Journal of Scientific Computing.

[7]  B. Fornberg,et al.  Theoretical and computational aspects of multivariate interpolation with increasingly flat radial basis functions , 2003 .

[8]  Martin Berzins.  Nonlinear data-bounded polynomial approximations and their applications in ENO methods , 2010, Numerical Algorithms.

[9]  Tobin A. Driscoll,et al.  Eigenvalue stability of radial basis function discretizations for time-dependent problems , 2006, Comput. Math. Appl..

[10]  Peng Li,et al.  Hybrid Fourier-Continuation Method and Weighted Essentially Non-oscillatory Finite Difference Scheme for Hyperbolic Conservation Laws in a Single-Domain Framework , 2015, J. Sci. Comput..

[11]  Wai-Sun Don,et al.  Multi-domain hybrid spectral-WENO methods for hyperbolic conservation laws , 2007, J. Comput. Phys..

[12]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[13]  Jingyang Guo,et al.  A numerical study of the local monotone polynomial edge detection for the hybrid WENO method , 2017, J. Comput. Appl. Math..

[14]  Yeon Ju Lee,et al.  An improved weighted essentially non-oscillatory scheme with a new smoothness indicator , 2013, J. Comput. Phys..

[15]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[16]  A. Harten.  Multiresolution representation of data: a general framework , 1996 .

[17]  Nathan Albin,et al.  Multi-domain Fourier-continuation/WENO hybrid solver for conservation laws , 2011, J. Comput. Phys..

[18]  Chao Yang,et al.  A new smoothness indicator for improving the weighted essentially non-oscillatory scheme , 2014, J. Comput. Phys..

[19]  Charles A. Micchelli,et al.  On Convergence of Flat Multivariate Interpolation by Translation Kernels with Finite Smoothness , 2014 .

[20]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[21]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[22]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .