Tomáš Dohnal

Tomáš Dohnal, Ben Schweizer

We present a new numerical scheme to solve the Helmholtz equation in a wave-guide. We consider a medium that is bounded in the $x_2$-direction, unbounded in the $x_1$-direction and $\varepsilon$-periodic for large $|x_1|$, allowing different media on the left and on the right. We suggest a new numerical method that is based on a truncation of the domain and the use of Bloch wave ansatz functions in radiation boxes. We prove the existence and a stability estimate for the infinite dimensional version of the proposed problem. The scheme is tested on several interfaces of homogeneous and periodic media and it is used to investigate the effect of negative refraction at the interface of a photonic crystal with a positive effective refractive index.

Tomáš Dohnal

This paper constructs perfectly matched layers (PML) for a system of 2D Coupled Nonlinear Schrödinger equations with mixed derivatives which arises in the modeling of gap solitons in nonlinear periodic structures with a non-separable linear part. The PML construction is performed in Laplace Fourier space via a modal analysis and can be viewed as a complex change of variables. The mixed derivatives cause the presence of waves with opposite phase and group velocities, which has previously been shown to cause instability of layer equations in certain types of hyperbolic problems. Nevertheless, here the PML is stable if the absorption function $σ$ lies below a specified threshold. The PML construction and analysis are carried out for the linear part of the system. Numerical tests are then performed in both the linear and nonlinear regimes checking convergence of the error with respect to the layer width and showing that the PML performs well even in many nonlinear simulations.