Section: New Results

Electromagnetic wave propagation

Numerical study of the non-linear Maxwell equations for Kerr media

Participants : Loula Fezoui, Stéphane Lanteri.

The system of Maxwell equations describes the evolution of the interaction of an electromagnetic field with a propagation medium. The different properties of the medium, such as isotropy, homogeneity, linearity, among others, are introduced through constitutive laws linking fields and inductions. In the present study, we focus on non-linear effects and address non-linear Kerr materials specifically. In this model, any dielectric may become non-linear provided the electric field in the material is strong enough. As a first setp, we considered the one-dimensional case and study the numerical solution of the non-linear Maxwell equations thanks to DG methods. In particular, we make use of an upwind scheme and limitation techniques because they have a proven ability to capture shocks and other kinds of singularities in the fluid dynamics framework. The numerical results obtained in this preliminary study gave us confidence towards extending them to higher spatial dimensions. This year, we have completed the development of a first version a parallel DGTD solver for the three-dimensional based on our past contributions on DGTD methods for the case of linear propagation media.

Numerical treatment of non-local dispersion for nanoplasmonics

Participants : Stéphane Lanteri, Claire Scheid, Nikolai Schmitt, Jonathan Viquerat.

When metallic nanostructures have sub-wavelength sizes and the illuminating frequencies are in the regime of metal's plasma frequency, electron interaction with the exciting fields have to be taken into account. Due to these interactions, plasmonic surface waves can be excited and cause extreme local field enhancements (surface plasmon polariton electromagnetic waves). Exploiting such field enhancements in applications of interest requires a detailed knowledge about the occurring fields which can generally not be obtained analytically. For the numerical modeling of light/matter interaction on the nanoscale, the choice of an appropriate model is a crucial point. Approaches that are adopted in a first instance are based on local (no interaction between electrons) dispersion models e.g. Drude or Drude-Lorentz. From the mathematical point of view, these models lead to an additional ordinary differential equation in time that is coupled to Maxwell's equations. When it comes to very small structures in a regime of 2 nm to 25 nm, non-local effects due to electron collisions have to be taken into account. Non-locality leads to additional, in general non-linear, partial differential equations and is significantly more difficult to treat, though. In this work, we study a DGTD method able to solve the system of Maxwell equations coupled to a linearized non-local dispersion model relevant to nanoplasmonics. This year, we have developed a parallel DGTD solver for the three-dimentional Maxwell equations coupeld to a non-local Drude model. Both centered flux-based and upwind flux-based DG schemes have been considered, in combination with with leap-frog and Runge-Kutta time stepping respectively.

Corner effects in nanoplasmonics

Participants : Camille Carvalho [ENSTA, POEMS project-team] , Patrick Ciarlet [ENSTA, POEMS project-team] , Claire Scheid.

In this work, we study nanoplasmonic structures with corners (typically a diedral/truangular structure). This is the central subject considered in the PhD thesis of Camille Carvalho. In the latter, the focus is made on a lossles Drude dispersion model with a frequency-domain approach. Several well posedness problems arise due to the presence of corners and are addressed in the PhD thesis. A time-domain approach in this context is also relevant and we propose to use the techniques developed in the team in this prospect. Even if both approaches (time-domain and frequency-domain) represent similar physical phenomena, problems that arise are different. These two approaches appear as complementary; it is thus worth bridging the gap between the two frameworks. We are currently performing a thorough comparison in the case of theses 2D structures with corners and we especially focus on the amplitude principle limit that raises a lot of questions.

Travelling waves for the non-linear Schrödinger equation in 2D

Participants : David Chiron [J.A. Dieudonné Laboratory, Université Nice Sophia Antipolis] , Claire Scheid.

Participants : Serge Nicaise [Université de Valenciennes et du Hainaut-Cambrésis] , Claire Scheid.

We are interested in the numerical study of the two-dimensional travelling waves of the non-linear Schrödinger equation for a general non-linearity and with nonzero condition at infinity. This equation is appearing in models of nonlinear optics. It has a variational structure that we propose to exploit to design a numerical method. We continue the sudy initiated in [1] and investigate excited states of the Kadomtsev-Petviashvili-I (KP-I) and Gross-Pitaevskii (GP) equations in dimension 2. We address numerically the question of the Morse index of some explicit solutions of KP-I. The results confirm that the lump solitary wave has Morse index one and that the other explicit solutions correspond to excited states. We then turn to the 2D GP equation which in some long wave regime converges to the KP-I equation. We finally perform numerical simulations showing that the other explicit solitary waves solutions to the KP-I equation give rise to new branches of travelling waves of GP corresponding to excited states.

In this ongoing work, we are interested in fundamental properties of the non local linearized hydrodynamic Drude model introduced in the context of nanoplasmonics. We propose an existence and detailed (polynomial/exponential) stability study for these models. We also investigate the discrete stability results. We propose to study the impact of the DG schemes developped in the team on these properties. This study complements the numerical approach that we already propose in the context of the PhD of Nikolai Schmitt for this model, towards a thorough understanding of its fundamentals properties.

A structure preserving numerical discretization framework for the Maxwell Klein Gordon equation in 2D.

Participants : Snorre Christiansen [Department of Mathematics, University of Oslo, Norway] , Claire Scheid.

Toward a better understanding of non-linear optical phenomena, we focus on the case of the Maxwell Klein Gordon (MKG) equation in dimension 2. This equation appears in the context of quantum electrodynamics but also in relativity. We propose to develop a numerical discretization framework that takes advantage of the Hamiltonian structure of the equation. The gauge invariance is recovered at the discrete level with the help of the Lattice Gauge theory. We then propose a fully discrete scheme and prove its convergence. The strategy of proof, based on discrete energy principle, is developed in a more general context and next applied in the particular case of MKG equation. This work has been conducted and finalized during a of five month's stay of C. Scheid at the University of Oslo through an invitation in the contexty of the ERC Starting Grant project STUCCOFIELD of S. Christiansen.

Multiscale DG methods for the time-domain Maxwell equations

Participants : Stéphane Lanteri, Raphaël Léger, Diego Paredes Concha [Instituto de Matemáticas, Universidad Católica de Valparaiso, Chile] , Claire Scheid, Frédéric Valentin [LNCC, Petropolis, Brazil] .

Although the DGTD method has already been successfully applied to complex electromagnetic wave propagation problems, its accuracy may seriously deteriorate on coarse meshes when the solution presents multiscale or high contrast features. In other physical contexts, such an issue has led to the concept of multiscale basis functions as a way to overcome such a drawback and allow numerical methods to be accurate on coarse meshes. The present work, which is conducted in the context of the HOMAR Associate Team, is concerned with the study of a particular family of multiscale methods, named Multiscale Hybrid-Mixed (MHM) methods. Initially proposed for fluid flow problems, MHM methods are a consequence of a hybridization procedure which caracterize the unknowns as a direct sum of a coarse (global) solution and the solutions to (local) problems with Neumann boundary conditions driven by the purposely introduced hybrid (dual) variable. As a result, the MHM method becomes a strategy that naturally incorporates multiple scales while providing solutions with high order accuracy for the primal and dual variables. The completely independent local problems are embedded in the upscaling procedure, and computational approximations may be naturally obtained in a parallel computing environment. In this study, a family of MHM methods is proposed for the solution of the time-domain Maxwell equations where the local problems are discretized either with a continuous FE method or a DG method (that can be viewed as a multiscale DGTD method). Preliminary results have been obtained in the two-dimensional case.

Figure 4. Light propagation in a photonic crystal structure using a MHM-DGTD method for solving the 2D Maxwell's equations. Left: quadrangular mesh. Right: contour lines of the amplitude of the electric field.
IMG/mhm_mesh-quad.png IMG/mhm_sol-quad.png

HDG methods for the time-domain Maxwell equations

Participants : Alexandra Christophe-Argenvillier, Stéphane Descombes, Stéphane Lanteri.

This study is concerned with the development of accurate and efficient solution strategies for the system of 3D time-domain Maxwell equations coupled to local dispersion models (e.g. Debye, Drude or Drude-Lorentz models) in the presence of locally refined meshes. Such meshes impose a constraint on the allowable time step for explicit time integration schemes that can be very restrictive for the simulation of 3D problems. We consider here the possibility of using an unconditionally stable implicit time or a locally implicit time integration scheme combined to a HDG discretization method.

HDG methods for the frequency-domain Maxwell equations

Participants : Alexis Gobé, Stéphane Lanteri, Ludovic Moya.

In the context of the ANR TECSER project, we continue our efforts towards the development of scalable high order HDG methods for the solution of the system of 3D frequency-domain Maxwell equations. We aim at fully exploiting the flexibiity of the HDG discretization framework with regards to the adaptation of the interpolation order (p-adaptivity) and the mesh (h-adaptivity). In particular, we study the formulation of HDG methods on a locally refined non-conforming tetrahedral mesh and on a non-confoming hybrid cubic/tetrahedral mesh. We also investigate the coupling between the HDG formulation and a BEM (Boundary Element Method) discretization of an integral representation of the electromagnetic field in the case of propagation problems theoretically defined in unbounded domains. The associated methodological contributions are implemented in the HORSE simulation software.

HDG methods for the frequency-domain plasmonics

Participants : Stéphane Lanteri, Liang Li [UESTC, Chengdu, China] , Asger Mortensen [DTU Fotonik, Technical University of Denmark] , Martijn Wubs [DTU Fotonik, Technical University of Denmark] .

In this colleaboration with physicists at DTU Fotonik, we study HDG methods for solving the frequency-domain Maxwell's equations coupled to the Nonlocal Hydrodynamic Drude (NHD) and Generalized Nonlocal Optical Response (GNOR) models, which are employed to describe the optical properties of nanoplasmonic scatterers and waveguides. The formulations of the HDG method for these two models are exetnsion of our previous works for classical microwave applications. In teh present case, two conservativity conditions are globally enforced to make the problem solvable and to guarantee the continuity of the tangential component of the electric field and the normal component of the current density. Numerical results show that the proposed HDG methods converge at optimal rate. These new HDG formulations hace been implemented and numerically assessed for two-dimensional problems.

Exponential time integrators for a DGTD method

Participants : Stéphane Descombes, Stéphane Lanteri, Bin Li [UESTC, Chengdu, China] , Hao Wang [UESTC, Chengdu, China] , Li Xu [UESTC, Chengdu, China] .

The objective of this study is to design efficient and (high order) accurate time integration strategies for the system of time-domain Maxwell equations discretized in space by a high order discontinuous Galerkin scheme formulated on locally refined unstructured meshes. A new family of implicit-explicit (IMEX) schemes using exponential time integration is developed. The Lawson procedure is applied based on a partitioning of the underlying tetrahedral mesh in coarse and fine parts, allowing the contruction of a time advancing strategy that combines an exact integration of the semi-discrete system for the problem unknowns associated to the elements of the fine part, with an arbitrary high order explicit time integration scheme for the Lawson-transformed system.