Abstract
AbstractThere exist several discretization techniques for the numerical solution of partial differential equations. In addition to classical finite difference, finite element and finite volume techniques, a more recent approach employs radial basis functions to generate differentiation stencils on unstructured point sets. This approach, abbreviated by RBF-FD (radial basis function-finite difference), has gained in popularity since it enjoys several advantages: It is (relatively) straightforward, does not require a mesh and generalizes easily to higher spatial dimensions. However, its application is not quite as blackbox as it may appear at first sight. The computed solution might suffer severely from various sources of errors if RBF-FD parameters are not selected carefully. Through comprehensive numerical experiments, we study the influence of several of these parameters on the condition numbers of intermediate (local) weight matrices, on the condition number of the resulting (global) stiffness matrix and ultimately on the approximation error of the computed discrete solution to the partial differential equation. The parameters of investigation include the type of RBF (and its shape or other parameters if applicable), the degree of polynomial augmentation, the discretization stencil size, the underlying type of point set (structured/unstructured), and the total number of (interior and boundary) points to discretize the PDE, here chosen as a three-dimensional Poisson’s problem with Dirichlet boundary conditions. Numerical tests on a sphere as well as tests for the convection-diffusion equation are included in a supplement and demonstrate that the results obtained for the Laplace problem on a cube generalize to wider problem classes. The purpose of this paper is to provide a comprehensive survey on the various components of the basic algorithms for RBF-FD discretization and steer away from potential pitfalls such as computationally more expensive setups which not always lead to more accurate numerical solutions. We guide toward a compatible selection of the multitude of RBF-FD parameters in the basic version of RBF-FD. For many of its components we refer to the literature for more advanced versions.
Funder
Technische Universität Hamburg
Publisher
Springer Science and Business Media LLC
Subject
Computational Theory and Mathematics,General Engineering,Theoretical Computer Science,Software,Applied Mathematics,Computational Mathematics,Numerical Analysis
Reference82 articles.
1. Baxter, B., Brummelhuis, R.: Convergence of stationary radial basis function-schemes for evolution equations, arXiv preprint arXiv:1905.01128 (2019)
2. Bayona, V.: An insight into RBF-FD approximations augmented with polynomials. Comput. Math. Appl. 77, 2337–2353 (2019)
3. Bayona, V., Flyer, N., Fornberg, B.: On the role of polynomials in RBF-FD approximations: III. Behavior near domain boundaries. J. Comput. Phys. 380, 378–399 (2019)
4. Bayona, V., Flyer, N., Fornberg, B., Barnett, G.A.: On the role of polynomials in RBF-FD approximations: II. Numerical solution of elliptic PDEs. J. Comput. Phys. 332, 257–273 (2017)
5. Bayona, V., Moscoso, M., Carretero, M., Kindelan, M.: RBF-FD formulas and convergence properties. J. Comput. Phys. 229, 8281–8295 (2010)
Cited by
4 articles.
订阅此论文施引文献
订阅此论文施引文献,注册后可以免费订阅5篇论文的施引文献,订阅后可以查看论文全部施引文献