The radial basis integral equation method for 2D Helmholtz problems

摘要

A meshless method for the solution of 2D Helmholtz equation has been developed by using the Boundary Integral Equation (BIE) combined with Radial Basis Function (RBF) interpolations. BIE is applied by using the fundamental solution of the Helmholtz equation, therefore domain integrals are not encountered in the method. The method exploits the advantage of placing the source point always in the centre of circular sub-domains in order to avoid singular or near-singular integrals. Three equations for two-dimensional (2D) or four for three-dimensional (3D) potential problems are required at each node. The first equation is the integral equation arising from the application of the Green's identities and the remaining equations are the derivatives of the first equation in respect to space coordinates. RBF interpolation is applied in order to obtain the values of the field variable and partial derivatives at the boundary of the circular sub-domains, providing in this way the boundary conditions for solution of the integral equations at the nodes (centres of circles). The accuracy and robustness of the method has been tested on some analytical solutions of the problem. Two different RBFs are used, namely f_1 (R) = R~2 ln(R) + x + y + 1 and f_2 (R) = R~4 ln(R) + x~2 + y~2 + xy + x + 1. The latter has been found to produce more accurate results.