Local multiple traces formulation for high-frequency scattering problems

Abstract

a b s t r a c tWe present an efficient method to solve high-frequency scattering problems by heterogeneous penetrable objects in two dimensions. This is achieved by extending the so-called Local Multiple Traces Formulation, introduced recently by Hiptmair and Jerez-Hanckes, to purely spectral discretizations employing weighted Chebyshev polynomials. Together with regularization strategies to handle boundary integral operators singularities, matrix entries are quickly computed via the Fast Fourier Transform. The resulting Fredholm first-kind formulation is free from spurious resonances, and though ill-conditioned, it possesses built-in Calderón-type preconditioners. Numerical results obtained for different settings validate our claims and greatly motivate future research in this direction.

Fast Spectral Galerkin Method for Logarithmic Singular Equations on a Segment

“…First, define the following first-order differential operator: 18) and denote by D n the composition D • . .…”

“…Numerical quadratures are computed as described in [18,Section 3.3]. Table 6.5: Helmholtz convergence when using an overkill ϕ ok with N ok dofs taking as right-hand side g = exp(ıx).…”

Local Multiple Traces Formulation for Electromagnetics: Stability and Preconditioning for Smooth Geometries

“…With this definition, we have (A j κ,µ ) 2 = Id. Now, let us rewrite (12) in a matrix form. We first introduce the continuous map…”

“…Though theoretical aspects for Maxwell scattering are now fully available for global MTFs [9,10] their implementation is extremely cumbersome. Opposingly, local versions of the MTF [1,11,12,13,14] are easily implemented and parallelized and though theoretical results for acoustic and static versions are available, in the electromagnetic case similar results remain elusive. Regarding Maxwell's equations, preconditioning is crucial as STF and MTFs incorporate electric field integral operators.…”

High-order Galerkin method for Helmholtz and Laplace problems on multiple open arcs

“…For implementation purposes, we follow the scheme introduced in [16] wherein all integral kernel singularities are subtracted. This gives rise to smooth functions and singular functions whose integrals are respectively computed via the Fast Fourier Transform (FFT) [21] or analytically using a Chebyshev polynomial expansion of the fundamental solution [11]. Recently, Slevinsky and Olver [33] devised a similar construction based on Chebyshev polynomials for more general integral equations, but limited to line segments and focused exclusively on the spectral properties of collocation method.…”