Author(s)

E. H. Ooi & V. Popov

Abstract

In this paper, we propose an efficient implementation of the radial basis integral equation method (RBIEM) that does not involve discretization of the circular subdomains. By avoiding discretization on the boundaries of the subdomains, a major source of error in the numerical scheme can be eliminated. The proposed implementation is tested on the Helmholtz equation with higher gradients in the exact solution. Three different radial basis functions are investigated, namely the augmented thin plate spline, r3 and r4log(r). The latter two functions are augmented with the second order global polynomial. Numerical results show that the new implementation of the RBIEM produces more accurate results and is more robust in handling problems with highly variable solutions. By avoiding the boundary discretization, the tasks of keeping track of the boundary elements and the boundary nodes are not needed, which can be a daunting task especially in three-dimensional problems with complicated geometries. The proposed implementation of the RBIEM is promising and the feasibility of the approach in three-dimensional problems is currently being investigated. Keywords: meshless, radial basis function, discretization, Helmholtz, efficiency. 1 Introduction In the past decade, research in numerical methods has been moving towards meshless approaches for solving partial differential equations. In particular, meshless methods based on integral equations have been gaining wide attention due to the accuracy of the integral equation based methods and the reduced requirements for meshing. Meshless methods based on integral equations such as the local boundary integral equation (LBIE) method developed by Zhu et

Keywords

meshless, radial basis function, discretization, Helmholtz, efficiency