Gibbs, A;
Hewett, DP;
Huybrechs, D;
Parolin, E;
(2020)
Fast hybrid numerical-asymptotic boundary element methods for high frequency screen and aperture problems based on least-squares collocation.
SN Partial Differential Equations and Applications
, 1
, Article 21. 10.1007/s42985-020-00013-3.
Preview |
Text
Gibbs2020_Article_FastHybridNumerical-asymptotic.pdf - Published Version Download (4MB) | Preview |
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.
| Type: | Article |
|---|---|
| Title: | Fast hybrid numerical-asymptotic boundary element methods for high frequency screen and aperture problems based on least-squares collocation |
| Open access status: | An open access version is available from UCL Discovery |
| DOI: | 10.1007/s42985-020-00013-3 |
| Publisher version: | https://doi.org/10.1007/s42985-020-00013-3 |
| Language: | English |
| Additional information: | This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. |
| Keywords: | High frequency scattering, Hybrid numerical-asymptotic boundary element method, Diffractal, Numerical steepest descent, Oscillatory quadrature |
| UCL classification: | UCL UCL > Provost and Vice Provost Offices UCL > Provost and Vice Provost Offices > UCL BEAMS UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Maths and Physical Sciences UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Maths and Physical Sciences > Dept of Mathematics |
| URI: | https://discovery.ucl.ac.uk/id/eprint/10114664 |
Archive Staff Only
 |
View Item |
