A Memory-Efficient Implementation of Perfectly Matched Layer With Smoothly Varying Coefficients in Discontinuous Galerkin Time-Domain Method

Wrapping a computation domain with a perfectly matched layer (PML) is one of the most effective methods of imitating/approximating the radiation boundary condition in the Maxwell and wave equation solvers. Many PML implementations often use a smoothly increasing attenuation coefficient to increase the absorption for a given layer thickness, and, at the same time, to reduce the numerical reflection from the interface between the computation domain and the PML. In discontinuous Galerkin time-domain (DGTD) methods, using a PML coefficient that varies within a mesh element requires a different mass matrix to be stored for every element and therefore significantly increases the memory footprint. In this work, this bottleneck is addressed by applying a weight-adjusted approximation to these mass matrices. The resulting DGTD scheme has the same advantages as the scheme that stores individual mass matrices, namely, higher accuracy (due to reduced numerical reflection) and increased meshing flexibility (since the PML does not have to be defined layer by layer), but it requires significantly less memory.

[1]  Maarten V. de Hoop,et al.  A weight-adjusted discontinuous Galerkin method for the poroelastic wave equation: Penalty fluxes and micro-heterogeneities , 2020, J. Comput. Phys..

[2]  Kurt Busch,et al.  Higher-order time-domain methods for the analysis of nano-photonic systems , 2009 .

[3]  Li Jun Jiang,et al.  Analysis of Electromagnetic Scattering from Penetrable Conductive Objects with IBC , 2015 .

[4]  H. Bağcı,et al.  A Discontinuous Galerkin Framework for Multiphysics Simulation of Photoconductive Devices , 2019, 2019 International Applied Computational Electromagnetics Society Symposium (ACES).

[5]  Stephen D. Gedney,et al.  Introduction to the Finite-Difference Time-Domain (FDTD) Method for Electromagnetics , 2011, Synthesis Lectures on Computational Electromagnetics.

[6]  Allen Taflove,et al.  Computational Electrodynamics the Finite-Difference Time-Domain Method , 1995 .

[7]  H. Bağcı,et al.  Exact Absorbing Boundary Conditions for Periodic Three-Dimensional Structures: Derivation and Implementation in Discontinuous Galerkin Time-Domain Method , 2018, IEEE Journal on Multiscale and Multiphysics Computational Techniques.

[8]  Timothy C. Warburton,et al.  Weight-Adjusted Discontinuous Galerkin Methods: Wave Propagation in Heterogeneous Media , 2016, SIAM J. Sci. Comput..

[9]  D. Angulo,et al.  Time domain discontinuous Galerkin methods for Maxwell equations , 2014 .

[10]  Jian-Ming Jin,et al.  Perfectly Matched Layers in the Discretized Space: An Analysis and Optimization , 1996 .

[11]  Jean-Pierre Bérenger,et al.  Perfectly Matched Layer (PML) for Computational Electromagnetics , 2007, PML for Computational Electromagnetics.

[12]  Pingwen Zhang,et al.  Discontinuous Galerkin methods for dispersive and lossy Maxwell's equations and PML boundary conditions , 2004 .

[13]  J. Hesthaven,et al.  Nodal high-order methods on unstructured grids , 2002 .

[14]  Jesse Chan,et al.  Bernstein-Bézier weight-adjusted discontinuous Galerkin methods for wave propagation in heterogeneous media , 2018, J. Comput. Phys..

[15]  Jian-Ming Jin,et al.  A General Scheme for the Discontinuous Galerkin Time-Domain Modeling and S-Parameter Extraction of Inhomogeneous Waveports , 2018, IEEE Transactions on Microwave Theory and Techniques.

[16]  Weng Cho Chew,et al.  A 3D perfectly matched medium from modified maxwell's equations with stretched coordinates , 1994 .

[18]  Hakan Bagci,et al.  Multiphysics Modeling of Plasmonic Photoconductive Devices using Discontinuous Galerkin Methods , 2019, ArXiv.

[19]  Stephen D. Gedney,et al.  The Discontinuous Galerkin Finite-Element Time-Domain Method Solution of Maxwell ’ s Equations , 2009 .

[20]  Hakan Bagci,et al.  Modeling Floating Potential Conductors Using Discontinuous Galerkin Method , 2020, IEEE Access.

[21]  Ping Li,et al.  A hybridizable discontinuous Galerkin method for simulation of electrostatic problems with floating potential conductors , 2020, ArXiv.

[22]  M. F. Pantoja,et al.  Discontinuous Galerkin Time Domain Methods in Computational Electrodynamics : State of the Art , 2015 .

[23]  H. Bağcı,et al.  A Unit-Cell Discontinuous Galerkin Scheme for Analyzing Plasmonic Photomixers , 2019, 2019 IEEE International Symposium on Antennas and Propagation and USNC-URSI Radio Science Meeting.

[24]  Liang Chen,et al.  Steady-State Simulation of Semiconductor Devices Using Discontinuous Galerkin Methods , 2019, IEEE Access.

[25]  Gary Cohen,et al.  Finite Element and Discontinuous Galerkin Methods for Transient Wave Equations , 2016 .

[26]  Vipin Kumar,et al.  A Parallel Algorithm for Multilevel Graph Partitioning and Sparse Matrix Ordering , 1998, J. Parallel Distributed Comput..

[27]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[28]  L. J. Jiang,et al.  A Resistive Boundary Condition Enhanced DGTD Scheme for the Transient Analysis of Graphene , 2015, IEEE Transactions on Antennas and Propagation.

[29]  Optimization of the Periodic PML for SEM , 2019, IEEE Transactions on Electromagnetic Compatibility.

[30]  Meilin Liu,et al.  Incorporation of Exact Boundary Conditions into a Discontinuous Galerkin Finite Element Method for Accurately Solving 2D Time-Dependent Maxwell Equations , 2013, IEEE Transactions on Antennas and Propagation.

[31]  Qing Huo Liu,et al.  A Novel Coupling Algorithm for Perfectly Matched Layer With Wave Equation-Based Discontinuous Galerkin Time-Domain Method , 2018, IEEE Transactions on Antennas and Propagation.

[32]  R. Cools.  Monomial cubature rules since “Stroud”: a compilation—part 2 , 1999 .

[33]  Hakan Bagci,et al.  Multiphysics Simulation of Plasmonic Photoconductive Devices Using Discontinuous Galerkin Methods , 2020, IEEE Journal on Multiscale and Multiphysics Computational Techniques.

[34]  Carretera de Valencia,et al.  The finite element method in electromagnetics , 2000 .

[35]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[36]  Krishnaswamy Sankaran,et al.  Accurate domain truncation techniques for time-domain conformal methods , 2007 .

[37]  Chi-Wang Shu,et al.  Locally divergence-free discontinuous Galerkin methods for the Maxwell equations , 2004, Journal of Computational Physics.

[38]  Li Jun Jiang,et al.  Transient Analysis of Dispersive Power-Ground Plate Pairs With Arbitrarily Shaped Antipads by the DGTD Method With Wave Port Excitation , 2017, IEEE Transactions on Electromagnetic Compatibility.

[39]  Meilin Liu,et al.  An Efficient Discontinuous Galerkin Finite Element Method for Highly Accurate Solution of Maxwell Equations , 2012, IEEE Transactions on Antennas and Propagation.

[40]  Zydrunas Gimbutas,et al.  A numerical algorithm for the construction of efficient quadrature rules in two and higher dimensions , 2010, Comput. Math. Appl..