Recursive moment computation in Filon methods and application to high-frequency wave scattering in two dimensions

DOI10.1093/IMANUM/DRAC067MaRDI QIDQ6190827zbMATH OpenOpenAlexFDO

Authors Georg Maierhofer, Arieh Iserles, Nigel Peake

Publication date 6 February 2024

Published in IMA Journal of Numerical Analysis (Search for Journal in Brave)

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

zbMATH Keywords

numerical integrationwave scatteringhighly oscillatory integralsFilon quadrature

Mathematics Subject Classification ID

Numerical quadrature and cubature formulas (65D32) Numerical methods for trigonometric approximation and interpolation (65T40) Approximate quadratures (41A55)

Abstract: We study the efficient approximation of highly oscillatory integrals using Filon methods. A crucial step in the implementation of these methods is the accurate and fast computation of the Filon quadrature moments. In this work we demonstrate how recurrences can be constructed for a wide class of oscillatory kernel functions, based on the observation that many physically relevant kernel functions are in the null space of a linear differential operator whose action on the Filon interpolation basis is represented by a banded (infinite) matrix. We discuss in further detail the application to two classes of particular interest, integrals with algebraic singularities and stationary points and integrals involving a Hankel function. We provide rigorous stability results for the moment computation for the first of these classes and demonstrate how the corresponding Filon method results in an accurate approximation at truly frequency-independent cost. For the Hankel kernel, we derive error estimates which describe the convergence behaviour of the method in terms of frequency and number of Filon quadrature points. Finally, we show how Filon methods with recursive moment computation can be applied to compute efficiently integrals arising in hybrid numerical-asymptotic collocation methods for high-frequency wave scattering on a screen.

Recommendations

• Quadrature methods for highly oscillatory singular integrals
• An adaptive Filon algorithm for highly oscillatory integrals
• A generalization of Filon-Clenshaw-Curtis quadrature for highly oscillatory integrals
• Efficient quadrature of highly oscillatory integrals using derivatives
• A combined Filon/asymptotic quadrature method for highly oscillatory problems

This page was built for publication: Recursive moment computation in Filon methods and application to high-frequency wave scattering in two dimensions

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