Error estimates in the fast multipole method for scattering problems. Part 2 : truncation of the Gegenbauer series

Carayol, Quentin  ; Collino, Francis

ESAIM: Mathematical Modelling and Numerical Analysis , Tome 39 (2005) no. 1, pp. 183-221

Résumé

We perform a complete study of the truncation error of the Gegenbauer series. This series yields an expansion of the Green kernel of the Helmholtz equation, $\frac{{\mathrm{e}}^{i|\overrightarrow{u}-\overrightarrow{v}|}}{4\pi i|\overrightarrow{u}-\overrightarrow{v}|}$, which is the core of the Fast Multipole Method for the integral equations. We consider the truncated series where the summation is performed over the indices $\ell \le L$. We prove that if $v=|\overrightarrow{v}|$ is large enough, the truncated series gives rise to an error lower than $ϵ$ as soon as $L$ satisfies $L+\frac{1}{2}\simeq v+C{W}^{\frac{2}{3}}\left(K\left(\alpha \right){ϵ}^{-\delta }{v}^{\gamma }\right)v^{\frac{1}{3}}$ where $W$ is the Lambert function, $K\left(\alpha \right)$ depends only on $\alpha =\frac{|\overrightarrow{u}|}{|\overrightarrow{v}|}$ and $C,\delta ,\gamma$ are pure positive constants. Numerical experiments show that this asymptotic is optimal. Those results are useful to provide sharp estimates of the error in the fast multipole method for scattering computation.

Détail

DOI : 10.1051/m2an:2005008

Classification : 33C10, 33C55, 41A80
Keywords: Gegenbauer, fast multipole method, truncation error