Abstract:<p>A novel fast multipole boundary element method (FM-BEM) is proposed to analyze 2-D potential problems by using linear and three-node quadratic elements. In FM-BEM, fast multipole expansions are used for the integrals on elements which are far away from the source point, whereas direct evaluations are used for the integrals on elements which are close to the source point. The use of higher-order elements results in more complex forms of the integrands, which increases the burden in direct evaluations, especially for singular and nearly singular integrals. Herein, the complex notation is introduced to simplify the computational formulations in boundary integral equations for 2-D potential problems. The singular and nearly singular integrals on linear elements are calculated by the analytic formulas, and those on three-node quadratic elements are evaluated by a robust semi-analytical algorithm. Numerical examples show that the proposed FM-BEM possesses higher accuracy than the conventional FM-BEM. Besides, the present method can analyze thin structures and evaluate accurately the physical quantities at interior points much closer to the boundary.</p>
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to improve the computational efficiency and accuracy of the fast multipole boundary element method (FM - BEM) in the analysis of two - dimensional potential problems, especially the accurate evaluation of singular and near - singular integrals when using high - order elements. Specifically:
1. **Limitations of the traditional BEM**: The traditional boundary element method (BEM) faces difficulties in dealing with large - scale problems because the matrices it generates are dense and asymmetric, resulting in a time complexity of O(𝑁²) for solving linear equations, with huge storage requirements and CPU time consumption.
2. **Introduction of the fast multipole method (FMM)**: By introducing the fast multipole method (FMM), the computational complexity can be reduced from O(𝑁²) to O(N), enabling the BEM to handle large - scale problems. However, when using high - order elements, the calculation of singular and near - singular integrals becomes complicated, affecting the computational efficiency and accuracy.
3. **Challenges brought by high - order elements**: When using linear and three - node quadratic elements, the integral form becomes more complex, especially for singular and near - singular integrals. The direct Gaussian quadrature method cannot effectively calculate these integrals, and special strategies are required to improve the computational accuracy.
4. **The proposed new method**: The paper proposes a new fast multipole boundary element method (FM - BEM), which combines linear and three - node quadratic elements. By introducing complex number representation, the calculation formula in the boundary integral equation is simplified, and a robust semi - analytical algorithm is developed to evaluate the singular and near - singular integrals on high - order elements.
5. **Specific contributions**:
- A calculation method based on complex number representation is proposed, which simplifies the calculation of singular and near - singular integrals on high - order elements.
- A semi - analytical algorithm is developed to evaluate the singular and near - singular integrals on three - node quadratic elements.
- The high efficiency and accuracy of the new method are verified through numerical experiments, especially in dealing with slender structures and evaluating physical quantities at internal points close to the boundary.
In conclusion, this paper aims to solve the computational efficiency and accuracy problems of high - order elements in the analysis of two - dimensional potential problems by improving the fast multipole boundary element method (FM - BEM), enabling it to handle large - scale and complex geometric shape problems more effectively.