High order complex contour discretization methods to simulate scattering problems in locally perturbed periodic waveguides

TL;DR

This paper introduces two high-order complex contour discretization methods for simulating wave scattering in locally perturbed periodic waveguides, addressing non-uniqueness issues and demonstrating convergence through numerical validation.

Contribution

The paper develops two novel high-order complex contour discretization techniques based on Floquet-Bloch transform for waveguide scattering problems with local perturbations.

Findings

Both methods effectively handle non-uniqueness due to guided modes.

High-order discretization improves computational performance.

Numerical examples confirm convergence and accuracy.

Abstract

In this paper, two high order complex contour discretization methods are proposed to simulate wave propagation in locally perturbed periodic closed waveguides. As is well known the problem is not always uniquely solvable due to the existence of guided modes. The limiting absorption principle is a standard way to get the unique physical solution. Both methods are based on the Floquet-Bloch transform which transforms the original problem to an equivalent family of cell problems. The first method, which is designed based on a complex contour integral of the inverse Floquet-Bloch transform, is called the CCI method. The second method, which comes from an explicit definition of the radiation condition, is called the decomposition method. Due to the local perturbation, the family of cell problems are coupled with respect to the Floquet parameter and the computational complexity becomes much larger. To this end, high order methods to discretize the complex contours are developed to have better performances. Finally we give the convergence results which we confirm with numerical examples.