Please use this identifier to cite or link to this item: http://bura.brunel.ac.uk/handle/2438/28755
Title: An efficient frequency-independent numerical method for computing the far-field pattern induced by polygonal obstacles
Authors: Gibbs, A
Langdon, S
Keywords: embedding formula;far-field pattern;scattering;Cauchy integral;oversampling
Issue Date: 5-Mar-2024
Publisher: Society for Industrial and Applied Mathematics (SIAM)
Citation: Gibbs, A. and Langdon, S. (2024) 'An efficient frequency-independent numerical method for computing the far-field pattern induced by polygonal obstacles', SIAM Journal on Scientific Computing, 0 (accepted, in press), pp. 1 - 26.
Abstract: For problems of time-harmonic scattering by rational polygonal obstacles, embedding formulae express the far-field pattern induced by any incident plane wave in terms of the far-field patterns for a relatively small (frequency-independent) set of canonical incident angles. Although these remarkable formulae are exact in theory, here we demonstrate that: (i) they are highly sensitive to numerical errors in practice, and (ii) direct calculation of the coefficients in these formulae may be impossible for particular sets of canonical incident angles, even in exact arithmetic. Only by overcoming these practical issues can embedding formulae provide a highly efficient approach to computing the far-field pattern induced by a large number of incident angles. Here we address challenges (i) and (ii), supporting our theory with numerical experiments. Challenge (i) is solved using techniques from computational complex analysis: we reformulate the embedding formula as a complex contour integral and prove that this is much less sensitive to numerical errors. In practice, this contour integral can be efficiently evaluated by residue calculus. Challenge (ii) is addressed using techniques from numerical linear algebra: we oversample, considering more canonical incident angles than are necessary, thus expanding the set of valid coefficient vectors. The coefficient vector can then be selected using either a least squares approach or column subset selection.
Description: MSC codes. 35J05, 78A45, 30E20, 65F20
The preprint archived on this institutional repository is available on arXiv at https://doi.org/10.48550/arXiv.2310.17603. It is licensed under a CC-BY licence (https://creativecommons.org/licenses/by/4.0/) to meet UKRI terms and conditions. It has not been certified by peer review. Cite as arXiv:2310.17603v2 [math.NA] for this version).
URI: https://bura.brunel.ac.uk/handle/2438/28755
DOI: https://doi.org/10.48550/arXiv.2310.17603
ISSN: 1064-8275
Appears in Collections:Dept of Mathematics Embargoed Research Papers

Files in This Item:
File Description SizeFormat 
Preprint.pdfEmbargoed until publication3.84 MBAdobe PDFView/Open


This item is licensed under a Creative Commons License Creative Commons