A fast direct solver for quasi-periodic scattering problems
From MaRDI portal
Publication:347892
DOI10.1016/J.JCP.2013.04.015zbMATH Open1349.78044arXiv1301.1247OpenAlexW2085884278MaRDI QIDQ347892
Author name not available (Why is that?)
Publication date: 5 December 2016
Published in: (Search for Journal in Brave)
Abstract: We consider the numerical solution of the scattering of time-harmonic plane waves from an infinite periodic array of reflection or transmission obstacles in a homogeneous background medium, in two dimensions. Boundary integral formulations are ideal since they reduce the problem to unknowns on the obstacle boundary. However, for complex geometries and/or higher frequencies the resulting dense linear system becomes large, ruling out dense direct methods, and often ill-conditioned (despite being 2nd-kind), rendering fast multipole-based iterative schemes also inefficient. We present an integral equation based solver with O(N) complexity, which handles such ill-conditioning, using recent advances in "fast" direct linear algebra to invert hierarchically the isolated obstacle matrix. This is combined with a recent periodizing scheme that is robust for all incident angles, including Wood's anomalies, based upon the free space Green's function kernel. The resulting solver is extremely efficient when multiple incident angles are needed, as occurs in many applications. Our numerical tests include a complicated obstacle several wavelengths in size, with and solution error of , where the solver is 66 times faster per incident angle than a fast multipole based iterative solution, and 600 times faster when incident angles are chosen to share Bloch phases.
Full work available at URL: https://arxiv.org/abs/1301.1247
No records found.
No records found.
This page was built for publication: A fast direct solver for quasi-periodic scattering problems
Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q347892)
