Authors : Hélène Barucq, M’Barek Fares (Cerfacs), Carola Kruse (Cerfacs), Sébastien Tordeux
Research article : Hélène Barucq, M’Barek Fares, Carola Kruse, Sébastien Tordeux, Sparsified discrete wave problem involving a radiation condition on a prolate spheroidal surface, IMA Journal of Numerical Analysis, Volume 41, Issue 1, January 2021, Pages 315–343 (link to the journal, link to the preprint)
We develop and analyse a high-order outgoing radiation boundary condition for solving three-dimensional scattering problems by elongated obstacles. This Dirichlet-to-Neumann condition is constructed using the classical method of separation of variables that allows one to define the scattered field in a truncated domain. It reads as an infinite series that is truncated for numerical purposes. The radiation condition is implemented in a finite element framework represented by a large dense matrix. Fortunately, the dense matrix can be decomposed into a full block matrix that involves the degrees of freedom on the exterior boundary and a sparse finite element matrix. The inversion of the full block is avoided by using a Sherman–Morrison algorithm that reduces the memory usage drastically. Despite being of high order, this method has only a low memory cost.
