Abstract:Conventionally, piecewise polynomials have been used in the boundary elements method (BEM) to approximate unknown boundary values. Since infinitely smooth radial basis functions (RBFs) are more stable and accurate than the polynomials for high dimensional domains, the unknown values are approximated by the RBFs in this paper. Therefore, a new formulation of BEM, called radial BEM, is obtained. To calculate singular boundary integrals of the new method, we propose a new distribution for boundary source points that removes singularity from the integrals. Therefore, the boundary integrals are calculated precisely by the standard Gaussian quadrature rule (GQR) with n = 16 quadrature nodes. Several numerical examples are presented to check the efficiency of the radial BEM versus standard BEM and RBF collocation method for solving partial differential equations (PDEs). Analytical and numerical studies presented in this paper admit the radial BEM as a perfect version of BEM which will enrich the contribution of BEM and RBFs in solving PDEs, impressively.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the poor stability and insufficient accuracy in high - dimensional domains of the polynomial basis functions used in the traditional boundary element method (BEM) when solving partial differential equations (PDEs). Specifically, the paper proposes a new method - the radial boundary element method (Radial BEM), which uses radial basis functions (RBFs) to approximate the unknown boundary values in order to improve the stability and accuracy of the calculation.
### Main contributions of the paper
1. **Introduction of radial basis functions (RBFs)**: The paper proposes to use infinitely smooth radial basis functions (RBFs) to replace the traditional polynomial basis functions in order to improve the stability and accuracy of solving partial differential equations in high - dimensional domains.
2. **Treatment of singular integrals**: In order to accurately calculate the singular boundary integrals in BEM, the paper proposes a new method for distributing boundary source points. By optimizing the positions of the source points, the singular integrals can be accurately calculated by the standard Gaussian quadrature rule (GQR).
3. **Verification by numerical experiments**: Through multiple numerical experiments, the paper verifies the high efficiency and accuracy of the radial BEM in solving partial differential equations, especially in dealing with singular integrals.
### Specific content of the paper
- **Background introduction**: The traditional BEM uses polynomial basis functions to approximate the unknown boundary values, but there are problems of poor stability and insufficient accuracy in high - dimensional domains. RBFs become an ideal choice to replace polynomial basis functions because of their high accuracy and stability in data approximation.
- **Method description**:
- **Derivation of the formula for radial BEM**: The paper derives in detail the BEM formula using RBFs and obtains a new form of radial BEM.
- **Treatment of singular integrals**: A new method for distributing boundary source points is proposed. By optimizing the positions of the source points, the singular integrals can be accurately calculated by the standard Gaussian quadrature rule.
- **Numerical experiments**:
- **Comparison of the accuracy of different RBFs**: Through numerical experiments, the accuracy of different types of RBFs (such as Gaussian RBF, inverse multi - quadratic RBF, thin - plate spline RBF, etc.) in solving the Laplace equation, Helmholtz equation and convection equation is compared.
- **Comparison with the traditional BEM**: Through numerical experiments, it is verified that the radial BEM is more efficient and accurate in solving partial differential equations, especially in dealing with singular integrals, which is superior to the traditional BEM.
- **Comparison with the RBF collocation method**: Through numerical experiments, it is verified that the radial BEM is superior to the RBF collocation method in terms of stability.
### Conclusion
The radial BEM method proposed in the paper shows higher stability and accuracy in solving partial differential equations, especially in dealing with singular integrals. By optimizing the positions of the boundary source points, the radial BEM can effectively overcome the problems of poor stability and insufficient accuracy of the traditional BEM in high - dimensional domains. The results of numerical experiments further verify the effectiveness and superiority of this method.