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

DOI10.1007/S42985-020-00013-3MaRDI QIDQ2218211zbMATH OpenOpenAlexFDO

Authors A. Gibbs, David P. Hewett, Daan Huybrechs, Emile Parolin

Publication date 15 January 2021

Published in SN Partial Differential Equations and Applications (Search for Journal in Brave)

Copyright license Creative Commons Attribution 4.0 International

Full work available at URL https://arxiv.org/abs/1912.09916

zbMATH Keywords

high frequency scatteringoscillatory quadraturenumerical steepest descentdiffractalhybrid numerical-asymptotic boundary element method

Mathematics Subject Classification ID

Diffraction, scattering (78A45) Numerical methods for integral equations (65R20) Boundary element methods for boundary value problems involving PDEs (65N38) Asymptotic analysis in optics and electromagnetic theory (78M35) Boundary element methods applied to problems in optics and electromagnetic theory (78M15)

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 (2015), pp.1698- 1728, 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.

Recommendations

• A frequency-independent boundary element method for scattering by two-dimensional screens and apertures
• Spectral Galerkin boundary element methods for high-frequency sound-hard scattering problems
• A hybrid numerical-asymptotic boundary element method for high frequency scattering by penetrable convex polygons
• A fast hybrid Galerkin method for high-frequency acoustic scattering
• A high frequency \(hp\) boundary element method for scattering by convex polygons

Cites work

• scientific article; zbMATH DE number 647478 (Why is no real title available?)
• scientific article; zbMATH DE number 936298 (Why is no real title available?)
• A Galerkin BEM for high-frequency scattering problems based on frequency-dependent changes of variables
• A Galerkin Boundary Element Method for High Frequency Scattering by Convex Polygons
• A collocation method for high-frequency scattering by convex polygons
• A frequency-independent boundary element method for scattering by two-dimensional screens and apertures
• A fully discrete Galerkin method for high frequency exterior acoustic scattering in three dimensions
• A high frequency \(hp\) boundary element method for scattering by convex polygons
• A high frequency boundary element method for scattering by a class of nonconvex obstacles
• A high frequency boundary element method for scattering by convex polygons with impedance boundary conditions
• A hybrid numerical-asymptotic boundary element method for high frequency scattering by penetrable convex polygons
• A hybrid numerical-asymptotic boundary integral method for high-frequency acoustic scattering
• An augmented galerkin procedure for the boundary integral method applied to two-dimensional screen and crack problems
• Applying an iterative method numerically to solve \(n \times n\) matrix Wiener-Hopf equations with exponential factors
• Approximation theory and approximation practice
• Asymptotically derived boundary elements for the Helmholtz equation in high frequencies
• Complex Gaussian quadrature of oscillatory integrals
• Computing Highly Oscillatory Integrals
• Diffractals
• Extraction of uniformly accurate phase functions across smooth shadow boundaries in high frequency scattering problems
• Filon-Clenshaw-Curtis rules for highly oscillatory integrals with algebraic singularities and stationary points
• Frames and numerical approximation
• Frequency-adapted Galerkin boundary element methods for convex scattering problems
• High frequency scattering by convex curvilinear polygons
• Highly oscillatory quadrature
• Hybrid numerical-asymptotic approximation for high-frequency scattering by penetrable convex polygons
• Improved Integral Formulation for Acoustic Radiation Problems
• Least squares data fitting with applications
• Mathematical and numerical modeling of wave propagation in fractal trees
• Numerical-asymptotic boundary integral methods in high-frequency acoustic scattering
• On generalized Gaussian quadrature rules for singular and nearly singular integrals
• On quadrature methods for highly oscillatory integrals and their implementation
• On the Evaluation of Highly Oscillatory Integrals by Analytic Continuation
• Piezoelectric ultrasonic transducers with fractal geometry
• Prescribed error tolerances within fixed computational times for scattering problems of arbitrarily high frequency: the convex case
• Shadow boundary effects in hybrid numerical-asymptotic methods for high-frequency scattering
• Solving Ill-Conditioned and Singular Linear Systems: A Tutorial on Regularization
• The AZ algorithm for least squares systems with a known incomplete generalized inverse
• Wavenumber-explicit continuity and coercivity estimates in acoustic scattering by planar screens
• Well-posed PDE and integral equation formulations for scattering by fractal screens

Cited in

(8)

• Numerical evaluation of oscillatory integrals via automated steepest descent contour deformation
• An efficient frequency-independent numerical method for computing the far-field pattern induced by polygonal obstacles
• Convergence analysis of oversampled collocation boundary element methods in 2D
• An analysis of least-squares oversampled collocation methods for compactly perturbed boundary integral equations in two dimensions
• Oversampling collocation method for the Volterra integral equation with contaminated data
• Fast and spectrally accurate numerical methods for perforated screens (with applications to Robin boundary conditions)
• A frequency-independent boundary element method for scattering by two-dimensional screens and apertures
• High-frequency asymptotic expansions for multiple scattering problems with Neumann boundary conditions

Describes a project that uses

Uses Software

• DLMF

This page was built for publication: Fast hybrid numerical-asymptotic boundary element methods for high frequency screen and aperture problems based on least-squares collocation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2218211)