Adaptive WENO methods based on radial basis functions reconstruction

We explore the use of radial basis functions (RBF) in the weighted essentially non-oscillatory (WENO) reconstruction process used to solve hyperbolic conservation laws, resulting in a numerical method of arbitrarily high order to solve problems with discontinuous solutions. Thanks to the mesh-less property of the RBFs, the method is suitable for non uniform grids and mesh adaptation. We focus on multiquadric radial basis functions and propose a simple strategy to choose the inherent shape parameter to control the balance between theoretical achievable accuracy and the numerical stability. We also develop an original smoothness indicator independent of the chosen RBF for the WENO reconstruction step. Moreover, we introduce type I and type II RBF-WENO methods by computing specific linear weights. The RBF-WENO method is used to solve linear and nonlinear problems for both scalar conservation laws and systems of conservation laws, including Burgers equation, the Buckley-Leverett equation, and the Euler equations. Numerical results confirm the performance of the proposed method. We finally consider an effective conservative adaptive algorithm that captures moving shocks and rapidly varying solutions well. Numerical results on moving grids are presented for both Burgers equation and the more complex Euler equations.


Published in:
Journal of Scientific Computing, 72, 3, 986-1020
Year:
2017
Publisher:
New York, Springer Verlag
ISSN:
0885-7474
Keywords:
Laboratories:




 Record created 2016-05-16, last modified 2018-09-13

Preprint:
Download fulltext
PDF

Rate this document:

Rate this document:
1
2
3
 
(Not yet reviewed)