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

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.