Singular boundary method

Fig. 1. Problem sketch and nodes distribution using the MFS: (a) interior problems, (b) exterior problems (please click to see big pictures)
Fig. 2. Problem sketch and nodes distribution using the SBM: (c) interior problems, (d) exterior problems (please click to see big pictures)

In numerical analysis, the singular boundary method (SBM) belongs to a family of meshless boundary collocation techniques which include the method of fundamental solutions (MFS),[1][2][3] boundary knot method (BKM),[4] regularized meshless method (RMM),[5] boundary particle method (BPM),[6] modified MFS,[7] and so on. This family of strong-form collocation methods is designed to avoid singular numerical integration and mesh generation in the traditional boundary element method (BEM) in the numerical solution of boundary value problems with boundary nodes, in which a fundamental solution of the governing equation is explicitly known.

The salient feature of the SBM is to overcome the fictitious boundary in the method of fundamental solution, while keeping all merits of the latter. The method offers several advantages over the classical domain or boundary discretization methods, among which are:

  • meshless. The method requires neither domain nor boundary meshing but boundary-only discretization points;
  • integration-free. The numerical integration of singular or nearly singular kernels could be otherwise troublesome, expensive, and complicated, as in the case, for example, the boundary element method;
  • boundary-only discretization for homogeneous problems. The SBM shares all the advantages of the BEM over domain discretization methods such as the finite element or finite difference methods;
  • to overcome the perplexing fictitious boundary in the method of fundamental solutions (see Figs. 1 and 2), thanks to the introduction of the concept of the origin intensity factor, which isolates the singularity of the fundamental solutions.

The SBM provides a significant and promising alternative to popular boundary-type methods such as the BEM and MFS, in particular, for infinite domain, wave, thin-walled structures, and inverse problems.

History of the singular boundary method

The methodology of the SBM was firstly proposed by Chen and his collaborators in 2009.[8][9] The basic idea is to introduce a concept of the origin intensity factor to isolate the singularity of the fundamental solutions so that the source points can be placed directly on the real boundary. In comparison, the method of fundamental solutions requires a fictitious boundary for placing the source points to avoid the singularity of fundamental solution. The SBM has since been successfully applied to a variety of physical problems, such as potential problems,[10][11] infinite domain problem,[12] Helmholtz problem,[13] and plane elasticity problem.[14]

There are the two techniques to evaluate the origin intensity factor. The first approach is to place a cluster of sample nodes inside the problem domain and to calculate the algebraic equations. The strategy leads to extra computational costs and makes the method is not as efficient as expected compared to the MFS. The second approach[15][16] is to employ a regularization technique to cancel the singularities of the fundamental solution and its derivatives. Consequently, the origin intensity factors can be determined directly without using any sample nodes. This scheme makes the method more stable, accurate, efficient, and extends its applicability.

Recent developments

Boundary layer effect problems

Like all the other boundary-type numerical methods, also it is observed that the SBM encounters a dramatic drop of solution accuracy at the region nearby boundary. Unlike singularity at origin, the fundamental solution at near-boundary regions remains finite. However, instead of being a flat function, the interpolation function develops a sharp peak as the field point approaches the boundary. Consequently, the kernels become “nearly singular” and can not accurately be calculated. This is similar to the so-called boundary layer effect encountered in the BEM-based methods.

A nonlinear transformation, based on the sinh function, can be employed to remove or damp out the rapid variations of the nearly singular kernels.[17] As a result, the troublesome boundary layer effect in the SBM has been successfully remedied. The implementation of this transformation is straightforward and can easily be embedded in existing SBM programs. For the test problems studied, very promising results are obtained even when the distance between the field point and the boundary is as small as 1×10−10.

Large-scale problems

Like the MFS and BEM, the SBM will produce dense coefficient matrices, whose operation count and the memory requirements for matrix equation buildup are of the order of O(N2) which is computationally too expensive to simulate large-scale problems.

The fast multipole method (FMM) can reduce both CPU time and memory requirement from O(N2) to O(N) or O(NlogN). With the help of FMM, the SBM can be fully capable of solving a large scale problem of several million unknowns on a desktop. This fast algorithm dramatically expands the applicable territory of the SBM to far greater problems than were previously possible.

See also

References

  1. ^ method of fundamental solutions (MFS)
  2. ^ Golberg MA, Chen CS, Ganesh M, "Particular solutions of 3D Helmholtz-type equations using compactly supported radial basis functions", Eng Anal Bound Elem 2000;24(7–8): 539–47.
  3. ^ Fairweather G, Karageorghis A, "The method of fundamental solutions for elliptic boundary value problems", Adv Comput Math 1998;9(1): 69–95.
  4. ^ Chen W, Tanaka M, "A meshless, integration-free, and boundary-only RBF technique Archived 2016-03-04 at the Wayback Machine", Comput Math Appl 2002;43(3–5): 379–91.
  5. ^ D.L. Young, K.H. Chen, C.W. Lee, "Novel meshless method for solving the potential problems with arbitrary domain", J Comput Phys 2005;209(1): 290–321.
  6. ^ boundary particle method (BPM)
  7. ^ Sarler B, "Solution of potential flow problems by the modified method of fundamental solutions: Formulations with the single layer and the double layer fundamental solutions", Eng Anal Bound Elem 2009;33(12): 1374–82.
  8. ^ Chen W, "Singular boundary method: A novel, simple, meshfree, boundary collocation numerical method", Chin J Solid Mech 2009;30(6): 592–9.
  9. ^ Chen W, Wang FZ, "A method of fundamental solutions without fictitious boundary Archived 2015-06-06 at the Wayback Machine", Eng Anal Bound Elem 2010;34(5): 530–32.
  10. ^ Wei X, Chen W, Fu ZJ, "Solving inhomogeneous problems by singular boundary method", J Mar SCI Tech 2012; 20(5).
  11. ^ Chen W, Fu ZJ, Wei X, "Potential Problems by Singular Boundary Method Satisfying Moment Condition", Comput Model Eng Sci 2009;54(1): 65–85.
  12. ^ Chen W, Fu Z, "A novel numerical method for infinite domain potential problems", Chin Sci Bull 2010;55(16): 1598–603.
  13. ^ Fu ZJ, Chen W, "A novel boundary meshless method for radiation and scattering problems", Advances in Boundary Element Techniques XI, Proceedings of the 11th international Conference, 12–14 July 2010, 83–90, Published by EC Ltd, United Kingdom (ISBN 978-0-9547783-7-8)
  14. ^ Gu Y, Chen W, Zhang CZ., "Singular boundary method for solving plane strain elastostatic problems", Int J Solids Struct 2011;48(18): 2549–56.
  15. ^ Chen W, Gu Y, "Recent advances on singular boundary method", Joint International Workshop on Trefftz Method VI and Method of Fundamental Solution II, Taiwan 2011.
  16. ^ Gu Y, Chen, W, "Improved singular boundary method for three dimensional potential problems", Chinese Journal of Theoretical and Applied Mechanics, 2012, 44(2): 351-360 (in Chinese)
  17. ^ Gu Y, Chen W, Zhang J, "Investigation on near-boundary solutions by singular boundary method", Eng Anal Bound Elem 2012;36(8): 117–82.