A hp-adaptive discontinuous Galerkin method for plasmonic waveguides

In this paper we propose and analyse a hphp-adaptive discontinuous finite element method for computing electromagnetic modes of propagation supported by waveguide structures comprised of a thin lossy metal film of finite width embedded in an infinite homogeneous dielectric. We propose a goal-oriented or dual weighted residual error estimator based on the solution of a dual problem that we use to drive the adaptive refinement with the aim to compute accurate approximation of the modes. We illustrate in the last section the benefits of the resulting hphp-adaptive method in practice, which consist in fast convergence and accurate estimation of the error. We tested the method computing the vanishing modes for a metallic waveguide of square section

Aplicação do método de Galerkin descontínuo para a análise de guias fotônicos

Orientador: Hugo Enrique Hernández FigueroaDissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Elétrica e de ComputaçãoResumo: Um novo método de onda completo para realizar a análise modal em guias de onda é introduzido nesta dissertação. A ideia central por trás do método é baseada na discretização da equação de onda vetorial com o Método de Galerkin Descontínuo com Penalidade Interior (IPDG, do inglês Interior Penalty Discontinuous Galerkin). Com uma função de penalidade apropriada, um método de alta precisão e sem modos espúrios é obtido. A eficiência do método proposto é provada em vários guias de onda, incluindo complicados guias de ondas ópticos com modos vazantes e também em guias de onda plasmônicos. Os resultados foram comparados com os métodos do estado-da-arte descritos na literatura. Também é discutida a importância dessa nova abordagem. Além disso, os resultados indicam que o método é mais preciso do que abordagens anteriores baseadas em Elementos Finitos. As principais contribuições deste trabalho são: foi desenvolvido um novo método robusto e de alta precisão para a análise de guias de ondas arbitrários, uma nova função de penalidade para o IPDG foi proposta e aplicações práticas do método proposto são apresentadas. Adicionalmente, no apêndice é apresentado uma aplicação da análise modal em simulação eletromagnética 3D com um método de Galerkin DescontínuoAbstract: A novel full-wave method to perform mode analysis on waveguides is introduced in this dissertation. The core of the method is based on an Interior Penalty Discontinuous Galerkin (IPDG) discretization of the vector wave equation. With an appropriate penalty function a spurious-free and high accuracy method is achieved. The efficiency of the proposed method was proved in several waveguides, including intricate optical waveguides with leaky modes and also on plasmonic waveguides. The obtained results were compared with the state-of-the-art mode solvers described in the literature. Also, a discussion on the importance of this new approach is presented. Moreover, the results indicate that the proposed method is more accurate than the previous approaches based on Finite Elements Methods. The main contributions of this work are: the development of a novel robust and accurate method for the analysis of arbitrary waveguides, a new penalty function for the IPDG was proposed and practical applications of the methods are discussed. In addition, in the appendix an application of modal analysis on 3D electromagnetic simulations with a Discontinuous Galerkin method is detailedMestradoTelecomunicações e TelemáticaMestre em Engenharia ElétricaCAPE

Different Approaches of Numerical Analysis of Electromagnetic Phenomena in Shaded Pole Motor with Application of Finite Elements Method

In this paper is used Finite Element Method-FEM for analysis of electromagnetic quantities of small micro motor – single phase shaded pole motor-SPSPM. FEM is widely used numerical method for solving nonlinear partial differential equations with variable coefficients. For that purpose motor model is developed with exact geometry and material’s characteristics. Two different approaches are applied in FEM analysis of electromagnetic phenomena inside the motor: magneto-static where all electromagnetic quantities are analysed in exact moment of time meaning frequency f=0 Hz and timeharmonic magnetic approach where the magnetic field inside the machine is time varying, meaning frequency f=50 Hz. Obtained results are presented and compared with available analytical result

A DGTD method for the numerical modeling of the interaction of light with nanometer scale metallic structures taking into account non-local dispersion effects

The interaction of light with metallic nanostructures is of increasing interest for various fields of research. When metallic structures 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 (e.g. 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 latter mentioned reason, numerical tools as well as a deeper understanding of the underlying physics, are absolutely necessary. For the numerical modeling of light/structure interaction on the nanoscale, the choice of an appropriate material model is a crucial point. Approaches that are adopted in a first instance are based on local (i.e. with no interaction between electrons) dispersive models e.g. Drude or Drude-Lorentz models. From the mathematical point of view, when a time-domain modeling is considered, these models lead to an additional system of ordinary differential equation which 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, system of partial differential equations and is significantly more difficult to treat, though. Nevertheless, dealing with a linear non-local dispersion model is already a setting that opens the route to numerous practical applications of plasmonics. In this work, we present a Discontinuous Galerkin Time-Domain (DGTD) method able to solve the system of Maxwell equations coupled to a linearized non-local dispersion model relevant to plasmonics. While the method is presented in the general 3d case, numerical results are given for 2d simulation settings only

Simulation of the interaction of light with 3‐D metallic nanostructures using a proper orthogonal decomposition‐Galerkin reduced‐order discontinuous Galerkin time‐domain method

International audienc

Discontinuous Galerkin discretised level set methods with applications to topology optimisation

This thesis presents research concerning level set methods discretised using discontinuous Galerkin (DG) methods. Whilst the context of this work is level set based topology optimisation, the main outcomes of the research concern advancements which are agnostic of application. The first of these outcomes are the development of two novel DG discretised PDE based level set reinitialisation techniques, the so called Elliptic and Parabolic reinitialisation methods, which are shown through experiment to be robust and satisfy theoretical optimal rates of convergence. A novel Runge-Kutta DG discretisation of a simplified level set evolution equation is presented which is shown through experiment to be high-order accurate for smooth problems (optimal error estimates do not yet exist in the literature based on the knowledge of the author). Narrow band level set methods are investigated, and a novel method for extending the level set function outside of the narrow band, based on the proposed Elliptic Reinitialisation method, is presented. Finally, a novel hp-adaptive scheme is developed for the DG discretised level set method driven by the degree with which the level set function can locally satisfy the Eikonal equation defining the level set reinitialisation problem. These component parts are thus combined to form a proposed DG discretised level set methodology, the efficacy of which is evaluated through the solution of numerous example problems. The thesis is concluded with a brief exploration of the proposed method for a minimum compliance design problem

Computational modeling and simulation of nonlinear electromagnetic and multiphysics problems

In this dissertation, nonlinear electromagnetic and multiphysics problems are modeled and simulated using various three-dimensional full-wave methods in the time domain. The problems under consideration fall into two categories. One is nonlinear electromagnetic problems with the nonlinearity embedded in either the permeability or the conductivity of the material's constitutive properties. The other is multiphysics problems that involve interactions between electromagnetic and other physical phenomena. A numerical solution of nonlinear magnetic problems is formulated using the three-dimensional time-domain finite element method (TDFEM) combined with the inverse Jiles-Atherton vector hysteresis model. A second-order nonlinear partial differential equation (PDE) that governs the nonlinear magnetic problem is constructed through the magnetic vector potential in the time domain, which is solved by applying the Newton-Raphson method. To solve the ordinary differential equation (ODE) representing the magnetic hysteresis accurately and efficiently, several ODE solvers are specifically designed and investigated. To improve the computational efficiency of the Newton-Raphson method, the multi-dimensional secant methods are incorporated in the nonlinear TDFEM solver. A nonuniform time-stepping scheme is also developed using the weighted residual approach to remove the requirement of a uniform time-step size during the simulation. Breakdown phenomena during high-power microwave (HPM) operation are investigated using different physical and mathematical models. During the breakdown process, the bound charges in solid dielectrics and air molecules break free and are pushed to move by the Lorentz force produced by the electromagnetic fields. The motion of free electrons produces plasma currents, which generate secondary electromagnetic fields that couple back to the externally applied fields and interact with the free electrons. When the incident field intensity is high enough, this will lead to an exponential increase of the charged particles known as breakdown. Such a process is first described by a nonlinear conductivity of the solid dielectric as a function of the electric field to model the dielectric breakdown phenomenon. The air breakdown problem encountered with HPM operation is then simulated with the plasma current modeled by a simplified plasma fluid equation. Both the dielectric and air breakdown problems are solved with the TDFEM together with a Newton's method, where the dielectric breakdown is treated as a pure nonlinear electromagnetic problem, while the air breakdown is treated as a multiphysics problem. To describe the plasma behavior more accurately, the plasma density and velocity are modeled by the equations of diffusion and motion, respectively. This results in a multiphysics and multiscale system depicted by the nonlinearly coupled full-wave Maxwell and plasma fluid equations, which are solved by a nodal discontinuous Galerkin time-domain (DGTD) method in three dimensions. The air breakdown during the HPM operation and the resulting plasma formation and shielding are modeled and simulated. Several important numerical issues in the simulation of nonlinear electromagnetic and multiphysics problems have been investigated and discussed. A continuity-preserving and divergence-cleaning scheme for electromagnetic problems involving inhomogeneous materials has been proposed based on the purely and damped hyperbolic Maxwell equations. A divergence-cleaning method is presented to enforce Gauss's laws and normal flux continuity by introducing auxiliary variables and damping terms into the original Maxwell's equations, which result in artificial propagation and dissipation of the numerical errors. Based on the DGTD method, dynamic h- and p-adaptation algorithms are developed for a full-wave analysis of electromagnetic and multiphysics problems. The dynamic h-adaptation algorithm can dynamically refine the mesh to resolve the local variation of the fields during the wave propagation, while the dynamic p-adaptation algorithm can determine and adjust the basis order in real time during the simulation. Both algorithms developed and investigated in this dissertation are highly flexible and efficient, and are powerful simulation tools in the solution of nonlinear electromagnetic and multiphysics problems

Discontinuous Galerkin Methods in Nanophotonics

In this thesis I present discontinuous Galerkin methods for Maxwell\u27s equations in both time- and frequency-domain. The method\u27s computational capabilities are extended by perfectly matched layers, dispersive and anisotropic materials, and sources. These techniques are applied to a study of coupling effects in split-ring resonator arrays