Quotient-space boundary element methods for scattering at complex screens

Journal Article (2021)
Author(s)

Xavier Claeys (Université Paris Diderot-Sorbonne Paris Cité)

Lorenzo Giacometti (CERN)

Ralf Hiptmair (ETH Zürich)

Carolina Urzúa-Torres (TU Delft - Numerical Analysis)

Research Group
Numerical Analysis
Copyright
© 2021 Xavier Claeys, Lorenzo Giacometti, Ralf Hiptmair, Carolina Urzúa-Torres
DOI related publication
https://doi.org/10.1007/s10543-021-00859-y
More Info
expand_more
Publication Year
2021
Language
English
Copyright
© 2021 Xavier Claeys, Lorenzo Giacometti, Ralf Hiptmair, Carolina Urzúa-Torres
Research Group
Numerical Analysis
Issue number
4
Volume number
61
Pages (from-to)
1193-1221
Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Abstract

A complex screen is an arrangement of panels that may not be even locally orientable because of junction lines. A comprehensive trace space framework for first-kind variational boundary integral equations on complex screens has been established in Claeys and Hiptmair (Integr Equ Oper Theory 77:167–197, 2013. https://doi.org/10.1007/s00020-013-2085-x) for the Helmholtz equation, and in Claeys and Hiptmair (Integr Equ Oper Theory 84:33–68, 2016. https://doi.org/10.1007/s00020-015-2242-5) for Maxwell’s equations in frequency domain. The gist is a quotient space perspective that allows to make sense of jumps of traces as factor spaces of multi-trace spaces modulo single-trace spaces without relying on orientation. This paves the way for formulating first-kind boundary integral equations in weak form posed on energy trace spaces. In this article we extend that idea to the Galerkin boundary element (BE) discretization of first-kind boundary integral equations. Instead of trying to approximate jumps directly, the new quotient space boundary element method employs a Galerkin BE approach in multi-trace boundary element spaces. This spawns discrete boundary integral equations with large null spaces comprised of single-trace functions. Yet, since the right-hand-sides of the linear systems of equations are consistent, Krylov subspace iterative solvers like GMRES are not affected by the presence of a kernel and still converge to a solution. This is strikingly confirmed by numerical tests.