The PML-Method for a Scattering Problem for a Local Perturbation of an Open Periodic Waveguide

Andreas Kirsch,Ruming Zhang

2024-01-01

Abstract:The perfectly matched layers method is a well known truncation technique for its efficiency and convenience in numerical implementations of wave scattering problems in unbounded domains. In this paper, we study the convergence of the perfectly matched layers (PML) for wave scattering from a local perturbation of an open waveguide in the half space above the real line, where the refractive index is a function which is periodic along the axis of the waveguide and equals to one above a finite height. The problem is challenging due to the existence of guided waves, and a typical way to deal with the difficulty is to apply the limiting absorption principle. Based on the Floquet-Bloch transform and a curve deformation theory, the solution from the limiting absorption principle is rewritten as the integral of a coupled family of quasi-periodic problems with respect to the quasi-periodicity parameter on a particularly designed curve. By comparing the Dirichlet-to-Neumann maps on a straight line above the locally perturbed periodic layer, we finally show that the PML method converges exponentially with respect to the PML parameter. Finally, the numerical examples are shown to illustrate the theoretical results.

Numerical Analysis,Analysis of PDEs

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the influence of local perturbations on wave scattering problems in open - periodic waveguides. Specifically, the author studied the convergence of the Perfectly Matched Layer (PML) method for wave scattering problems in the two - dimensional half - space $\mathbb{R}^2_+$, when the refractive index function $n$ varies periodically along the waveguide axis direction and is equal to 1 above a certain finite height.

### Problem Background
In open - periodic waveguides, due to the existence of guided - wave modes, wave scattering problems become very complex. The traditional treatment method is to apply the Limiting Absorption Principle (LAP), but this method is more difficult in numerical implementation. Therefore, the author proposes to use the Perfectly Matched Layer (PML) method to truncate the infinite domain, thereby simplifying the solution of the problem.

### Research Objectives
1. **Convergence of the PML method**: Prove the exponential convergence of the PML method when dealing with wave scattering problems.
2. **Numerical verification**: Verify the correctness of the theoretical results through numerical examples.

### Main Methods
1. **Floquet - Bloch transformation**: Transform the original problem into an integral form of a series of quasi - periodic problems.
2. **Curve deformation technique**: By designing specific curves, make each quasi - periodic problem have a unique solution.
3. **Dirichlet - to - Neumann mapping comparison**: Prove the exponential convergence of the PML method by comparing the Dirichlet - to - Neumann mappings along the deformed curves.

### Key Formulas
- **Wave equation**:
\[
\Delta u + k^2 (n + q) u=-f \quad \text{in} \quad \mathbb{R}^2_+, \quad u = 0 \quad \text{for} \quad x_2 = 0
\]
- **Variational form**:
\[
\int_{\mathbb{R}^2_+} (\nabla u\cdot\nabla\psi - k^2 (n + q) u\psi) \, dx=\int_Q f\psi \, dx
\]
- **Equation after Floquet - Bloch transformation**:
\[
\Delta\tilde{u}(x,\alpha)+2i\alpha\partial_1\tilde{u}(x,\alpha)+[(k + i\epsilon)^2 n-\alpha^2]\tilde{u}(x,\alpha)=-f e^{-i\alpha x_1}-(k + i\epsilon)^2 q u e^{-i\alpha x_1}
\]

### Conclusions
By combining the Floquet - Bloch transformation and the curve deformation technique, the author successfully transformed the wave scattering problem into an integral form of a series of quasi - periodic problems. By comparing the Dirichlet - to - Neumann mappings along the deformed curves, the exponential convergence of the PML method in dealing with wave scattering problems was finally proved. Numerical experiments further verified the correctness of this theoretical result.

### Application Prospects
This research not only provides an effective numerical method for wave scattering problems in open - periodic waveguides, but also provides a theoretical basis and technical support for research in other related fields.