Fast hybrid numerical-asymptotic boundary element methods for high frequency screen and aperture problems based on least-squares collocation

A. Gibbs (), D. P. Hewett, D. Huybrechs and E. Parolin
Additional contact information
A. Gibbs: University College London
D. P. Hewett: University College London
D. Huybrechs: KU Leuven
E. Parolin: Institut Polytechnique de Paris

Partial Differential Equations and Applications, 2020, vol. 1, issue 4, 1-26

Abstract: Abstract We present a hybrid numerical-asymptotic (HNA) boundary element method (BEM) for high frequency scattering by two-dimensional screens and apertures, whose computational cost to achieve any prescribed accuracy remains bounded with increasing frequency. Our method is a collocation implementation of the high order hp HNA approximation space of Hewett et al. (IMA J Numer Anal 35:1698–1728, 2015), where a Galerkin implementation was studied. An advantage of the current collocation scheme is that the one-dimensional highly oscillatory singular integrals appearing in the BEM matrix entries are significantly easier to evaluate than the two-dimensional integrals appearing in the Galerkin case, which leads to much faster computation times. Here we compute the required integrals at frequency-independent cost using the numerical method of steepest descent, which involves complex contour deformation. The change from Galerkin to collocation is nontrivial because naive collocation implementations based on square linear systems suffer from severe numerical instabilities associated with the numerical redundancy of the HNA basis, which produces highly ill-conditioned BEM matrices. In this paper we show how these instabilities can be removed by oversampling, and solving the resulting overdetermined collocation system in a weighted least-squares sense using a truncated singular value decomposition. On the basis of our numerical experiments, the amount of oversampling required to stabilise the method is modest (around 25% typically suffices), and independent of frequency. As an application of our method we present numerical results for high frequency scattering by prefractal approximations to the middle-third Cantor set.

Keywords: High frequency scattering; Hybrid numerical-asymptotic boundary element method; Diffractal; Numerical steepest descent; Oscillatory quadrature; 65N38; 65R20; 78A45; 78M15; 78M35 (search for similar items in EconPapers)
Date: 2020
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
http://link.springer.com/10.1007/s42985-020-00013-3 Abstract (text/html)
Access to the full text of the articles in this series is restricted.

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:spr:pardea:v:1:y:2020:i:4:d:10.1007_s42985-020-00013-3

Ordering information: This journal article can be ordered from
https://www.springer.com/journal/42985/

DOI: 10.1007/s42985-020-00013-3

Access Statistics for this article

Partial Differential Equations and Applications is currently edited by Zhitao Zhang

More articles in Partial Differential Equations and Applications from Springer
Bibliographic data for series maintained by Sonal Shukla () and Springer Nature Abstracting and Indexing ().