Boundary particle method


In applied mathematics, the boundary particle method is a boundary-only meshless collocation technique, in the sense that none of inner nodes are required in the numerical solution of nonhomogeneous partial differential equations. Numerical experiments show that the BPM has spectral convergence. Its interpolation matrix can be symmetric.

History and recent developments

In recent decades, the dual reciprocity method and multiple reciprocity method have been emerging as promising techniques to evaluate the particular solution of nonhomogeneous partial differential equations in conjunction with the boundary discretization techniques, such as boundary element method. For instance, the so-called DR-BEM and MR-BEM are popular BEM techniques in the numerical solution of nonhomogeneous problems.
The DRM has become a common method to evaluate the particular solution. However, the DRM requires inner nodes to guarantee the convergence and stability. The MRM has an advantage over the DRM in that it does not require using inner nodes for nonhomogeneous problems. Compared with the DRM, the MRM is computationally more expensive in the construction of the interpolation matrices and has limited applicability to general nonhomogeneous problems due to its conventional use of high-order Laplacian operators in the annihilation process.
The recursive composite multiple reciprocity method, was proposed to overcome the above-mentioned problems. The key idea of the RC-MRM is to employ high-order composite differential operators instead of high-order Laplacian operators to eliminate a number of nonhomogeneous terms in the governing equation. The RC-MRM uses the recursive structures of the MRM interpolation matrix to reduce computational costs.
The boundary particle method is a boundary-only discretization of an inhomogeneous partial differential equation by combining the RC-MRM with strong-form meshless boundary collocation discretization schemes, such as the method of fundamental solution, boundary knot method, regularized meshless method, singular boundary method, and Trefftz method. The BPM has been applied to problems such as nonhomogeneous Helmholtz equation and convection–diffusion equation. The BPM interpolation representation is of a wavelet series.
For the application of the BPM to Helmholtz, Poisson and plate bending problems, the high-order fundamental solution or general solution, harmonic function or Trefftz function are often used, for instance, those of Berger, Winkler, and vibrational thin plate equations. The method has been applied to inverse Cauchy problem associated with Poisson and nonhomogeneous Helmholtz equations.

Further comments

The BPM may encounter difficulty in the solution of problems having complex source functions, such as non-smooth, large-gradient functions, or a set of discrete measured data. The solution of such problems involves:
The complex functions or a set of discrete measured data can be interpolated by a sum of polynomial or trigonometric function series. Then, the RC-MRM can reduce the nonhomogeneous equation to a high-order homogeneous equation, and the BPM can be implemented to solve these problems with boundary-only discretization.
The domain decomposition may be used to in the BPM boundary-only solution of large-gradient source functions problems.

Free software and Matlab codes