Local RBF Algorithms for Elliptic Boundary Value Problems in Annular Domains

Document Type

Article

Publication Date

1-1-2019

Department

Mathematics

School

Mathematics and Natural Sciences

Abstract

A local radial basis function method (LRBF) is applied for the solution of boundary value problems in annular domains governed by the Poisson equation, the inhomogeneous biharmonic equation and the inhomogeneous Cauchy-Navier equations of elasticity. By appropriately choosing the collocation points we obtain linear systems in which the coefficient matrices possess block sparse circulant structures and which can be solved efficiently using matrix decomposition algorithms (MDAs) and fast Fourier transforms (FFTs). The MDAs used are appropriately modified to take into account the sparsity of the arrays involved in the discretization. The leave-oneout cross validation (LOOCV) algorithm is employed to obtain a suitable value for the shape parameter in the radial basis functions (RBFs) used. The selection of the nearest centres for each local influence domain is carried out using a modification of the kdtree algorithm. In several numerical experiments, it is demonstrated that the proposed algorithm is both accurate and capable of solving large scale problems.

Publication Title

Communications In Computational Physics

Volume

25

Issue

1

First Page

41

Last Page

67

Link to Full Text

Find in your library

Share

COinS