Collocation BEM for 3D Laplace Equation

Download CBEM_LAP, a package for solving the 3D Laplace equation based on a piecewise constant Collocation Boundary Element Method.

 

Piecewise constant collocation method

The potential $u$ and flux $t$ are assumed to be constant over each boundary element. It can be shown (see cbem_lapGuide.pdf provided in the package CBEM_LAP) that the potentials

$\displaystyle g_j(\boldsymbol{x}) = \int_{\Gamma _{\!j}}G(\boldsymbol{x},\boldsymbol{y}) {\rm d}\Gamma _{\!\boldsymbol{y}},$   and$\displaystyle \quad h_j(\boldsymbol{x}) = \int_{\Gamma _{\!j}}\boldsymbol{H}(\b... ...{y}) {\rm d}\Gamma _{\!\boldsymbol{y}},\quad\boldsymbol{x}\!\in\!\mathbb{R}^3,$   

due to a uniform source distribution over a flat triangle $\Gamma _j,$ can be employed to successfully compute $u$ and $t$ on $\Gamma =\bigcup \overline{\Gamma }_{\!j}$. In addition, these potentials can be utilized to effectively calculate $u$ at interior points $\boldsymbol{x}\!\in\!\Omega .$ The analytic expressions for$g_j(\boldsymbol{x})$ and$h_j(\boldsymbol{x})$ over a flat triangle are given in [1].

 

References

[1] S. Nintcheu Fata.
     Explicit expressions for 3D boundary integrals in potential theory.
     Int. J. Num. Meth. Eng., 78(1):32-47, 2009.