A fast direct solver for two dimensional quasi-periodic multilayered medium scattering problems

07/14/2019
by   Yabin Zhang, et al.
0

This manuscript presents a fast direct solution technique for solving two dimensional wave scattering problems from quasi-periodic multilayered structures. The fast solver is built from the linear system that results from the discretization of a boundary integral formulation that is robust at Wood's anomalies. When the interface geometries are complex, the linear system is too large to be handled via dense linear algebra. The key building block of the proposed solver is a fast direct direct solver for the large sparse block system that corresponds to the discretization of boundary integral equations. The solver makes use of hierarchical matrix inversion techniques, has a cost that scales linearly with respect to the number of unknowns on the interfaces and the precomputation can be used for all choices of boundary data. By partitioning the remainder of the precomputation into parts based on their dependence on incident angle, the proposed direct solver is efficient for problems involving many incident angles like those that arise in applications. For example for a problem on an eleven layer geometry where the solution is desired for 287 incident angles, the proposed solution technique is 87 times faster than building a new fast direct solver for each new incident angle. An additional feature of the proposed solution technique is that solving a problem where an interface or layer property is changed require an update in the precomputation that cost linearly with respect to the number points on the affected interfaces with a small constant prefactor. The efficiency for modified geometries and multiple solves make the solution technique well suited for optimal design and inverse scattering applications.

READ FULL TEXT
research
04/13/2022

A fast direct solver for two dimensional quasi-periodic multilayered media scattering problems, Part II

This manuscript is the second in a series presenting fast direct solutio...
research
11/26/2022

Robust fast direct integral equation solver for three-dimensional quasi-periodic scattering problems with a large number of layers

A boundary integral equation method for the 3-D Helmholtz equation in mu...
research
12/04/2021

An iterative solver for the HPS discretization applied to three dimensional Helmholtz problems

This manuscript presents an efficient solver for the linear system that ...
research
01/07/2022

A Direct Parallel-in-Time Quasi-Boundary Value Method for Inverse Space-Dependent Source Problems

Inverse source problems arise often in real-world applications, such as ...
research
10/12/2021

Fast Block Linear System Solver Using Q-Learning Schduling for Unified Dynamic Power System Simulations

We present a fast block direct solver for the unified dynamic simulation...
research
11/05/2020

Fast Solver for Quasi-Periodic 2D-Helmholtz Scattering in Layered Media

We present a fast spectral Galerkin scheme for the discretization of bou...
research
04/23/2017

Boundary First Flattening

A conformal flattening maps a curved surface to the plane without distor...

Please sign up or login with your details

Forgot password? Click here to reset