A Dirichlet-to-Neumann approach for the exact computation of guided modes in photonic crystal waveguides

This works deals with one dimensional infinite perturbation - namely line defects - in periodic media. In optics, such defects are created to construct an (open) waveguide that concentrates light. The existence and the computation of the eigenmodes is a crucial issue. This is related to a self-adjoint eigenvalue problem associated to a PDE in an unbounded domain (in the directions orthogonal to the line defect), which makes both the analysis and the computations more complex. Using a Dirichlet-to-Neumann (DtN) approach, we show that this problem is equivalent to one set on a small neighborhood of the defect. On contrary to existing methods, this one is exact but there is a price to be paid : the reduction of the problem leads to a nonlinear eigenvalue problem of a fixed point nature

Guided modes in a hexagonal periodic graph like domain

This paper deals with the existence of guided waves and edge states in particular two-dimensional media obtained by perturbing a reference periodic medium with honeycomb symmetry. This reference medium is a thin periodic domain (the thickness is denoted $\delta$ > 0) with an hexagonal structure, which is close to an honeycomb quantum graph. In a first step, we show the existence of Dirac points (conical crossings) at arbitrarily large frequencies if $\delta$ is chosen small enough. We then perturbe the domain by cutting the perfectly periodic medium along the so-called zigzag direction, and we consider either Dirichlet or Neumann boundary conditions on the cut edge. In the two cases, we prove the existence of edges modes as well as their robustness with respect to some perturbations, namely the location of the cut and the thickness of the perturbed edge. In particular, we show that different locations of the cut lead to almost-non dispersive edge states, the number of locations increasing with the frequency. All the results are obtained via asymptotic analysis and semi-explicit computations done on the limit quantum graph. Numerical simulations illustrate the theoretical results

Time-harmonic wave propagation in junctions of two periodic half-spaces

We are interested in the Helmholtz equation in a junction of two periodic half-spaces. When the overall medium is periodic in the direction of the interface, Fliss and Joly (2019) proposed a method which consists in applying a partial Floquet-Bloch transform along the interface, to obtain a family of waveguide problems parameterized by the Floquet variable. In this paper, we consider two model configurations where the medium is no longer periodic in the direction of the interface. Inspired by the works of G\'erard-Varet and Masmoudi (2011, 2012), and Blanc, Le Bris, and Lions (2015), we use the fact that the overall medium has a so-called quasiperiodic structure, in the sense that it is the restriction of a higher dimensional periodic medium. Accordingly, the Helmholtz equation is lifted onto a higher dimensional problem with coefficients that are periodic along the interface. This periodicity property allows us to adapt the tools previously developed for periodic media. However, the augmented PDE is elliptically degenerate (in the sense of the principal part of its differential operator) and thus more delicate to analyse.Comment: 60 pages, 29 figure

Time harmonic wave propagation in one dimensional weakly randomly perturbed periodic media

International audienceIn this work we consider the solution of the time harmonic wave equation in a one dimensional periodic medium with weak random perturbations. More precisely, we study two types of weak perturbations: (1) the case of stationary, ergodic and oscillating coefficients, the typical size of the oscillations being small compared to the wavelength and (2) the case of rare random perturbations of the medium, where each period has a small probability to have its coefficients modified, independently of the other periods. Our goal is to derive an asymptotic approximation of the solution with respect to the small parameter. This can be used in order to construct absorbing boundary conditions for such media

Trapped modes in thin and infinite ladder like domains. Part 1 : existence results

International audienceThe present paper deals with the wave propagation in a particular two dimensional structure, obtained from a localized perturbation of a reference periodic medium. This reference medium is a ladder like domain, namely a thin periodic structure (the thickness being characterized by a small parameter $\epsilon > 0$) whose limit (as $\epsilon$ tends to 0) is a periodic graph. The localized perturbation consists in changing the geometry of the reference medium by modifying the thickness of one rung of the ladder. Considering the scalar Helmholtz equation with Neumann boundary conditions in this domain, we wonder whether such a geometrical perturbation is able to produce localized eigenmodes. To address this question, we use a standard approach of asymptotic analysis that consists of three main steps. We first find the formal limit of the eigenvalue problem as the $\epsilon$ tends to 0. In the present case, it corresponds to an eigenvalue problem for a second order differential operator defined along the periodic graph. Then, we proceed to an explicit calculation of the spectrum of the limit operator. Finally, we prove that the spectrum of the initial operator is close to the spectrum of the limit operator. In particular, we prove the existence of localized modes provided that the geometrical perturbation consists in diminishing the width of one rung of the periodic thin structure. Moreover, in that case, it is possible to create as many eigenvalues as one wants, provided that ε is small enough. Numerical experiments illustrate the theoretical results

On the identification of defects in a periodic waveguide from far field data

International audienceThe aim of this paper is to apply the Linear Sampling Method and the Factorization Method to retrieve some defects in a known periodic 2D waveguide from scattering data. More precisely, some far field approximations of these two sampling methods are derived. They amount to consider the so-called propagating Floquet modes as incident waves. The efficiency of the far field formulation of the LSM is shown with the help of some numerical experiments

Solutions of the time-harmonic wave equation in periodic waveguides : asymptotic behaviour and radiation condition

International audienceIn this paper, we give the expression and the asymptotic behaviour of the physical solution of a time harmonic wave equation set in a periodic waveguide. This enables us to define a radiation condition and show well-posedness of the Helmholtz equation set in a periodic waveguide

Iterative methods for scattering problems in isotropic or anisotropic elastic waveguides

International audienceWe consider the time-harmonic problem of the diffraction of an incident propagative mode by a localized defect, in an infinite straight isotropic elastic waveguide. We propose several iterative algorithms to compute an approximate solution of the problem, using a classical finite element discretization in a small area around the perturbation, and a modal expansion in unbounded straight parts of the guide. Each algorithm can be related to a so-called domain decomposition method, with or without an overlap between the domains. Specific transmission conditions are used, so that only the sparse finite element matrix has to be inverted, the modal expansion being obtained by a simple projection, using the Fraser bi-orthogonality relation. The benefit of using an overlap between the finite element domain and the modal domain is emphasized, in particular for the extension to the anisotropic case. The transparency of these new boundary conditions is checked for two- and three-dimensional anisotropic waveguides. Finally, in the isotropic case, numerical validation for two- and three-dimensional waveguides illustrates the efficiency of the new approach, compared to other existing methods, in terms of number of iterations and CPU time

On the approximation of electromagnetic fields by edge finite elements. Part 2: A heterogeneous ultiscale method for Maxwell’s equations

In the second part of this series of papers we consider highly oscillatory media. In this situation, the need for a triangulation that resolves all microscopic details of the medium makes standard edge finite elements impractical because of the resulting tremendous computational load. On the other hand, undersampling by using a coarse mesh might lead to inaccurate results. To overcome these diffculties and to improve the ratio between accuracy and computational costs, homogenization techniques can be used. In this paper we recall analytical homogenization results and propose a novel numerical homogenization scheme for Maxwell\u27s equations in frequency domain. This scheme follows the design principles of heterogeneous multiscale methods. We prove convergence to the effective solution of the multiscale Maxwell\u27s equations in a periodic setting and give numerical experiments in accordance to the stated results

On the Approximation of Electromagnetic Fields by Edge Finite Elements. Part 2: A Heterogeneous Multiscale Method for Maxwell's equations

International audienceIn the second part of this series of papers we consider highly oscillatory media. In this situation, the need for a triangulation that resolves all microscopic details of the medium makes standard edge finite elements impractical because of the resulting tremendous computational load. On the other hand, undersampling by using a coarse mesh might lead to inaccurate results. To overcome these difficulties and to improve the ratio between accuracy and computational costs, homogenization techniques can be used. In this paper we recall analytical homogenization results and propose a novel numerical homogenization scheme for Maxwell’s equations in frequency domain. This scheme follows the design principles of heterogeneous multiscale methods. We prove convergence to the effective solution of the multiscale Maxwell’s equations in a periodic setting and give numerical experiments in accordance to the stated results