A Locally-Implicit Discontinuous Galerkin Time-Domain Method to Simulate Metasurfaces Using Generalized Sheet Transition Conditions

Creative Commons License

Chen L., ÖZAKIN M. B., Zhao R., Bagci H.

IEEE Transactions on Antennas and Propagation, vol.71, no.1, pp.869-881, 2023 (SCI-Expanded, Scopus) identifier identifier

  • Publication Type: Article / Article
  • Volume: 71 Issue: 1
  • Publication Date: 2023
  • Doi Number: 10.1109/tap.2022.3215079
  • Journal Name: IEEE Transactions on Antennas and Propagation
  • Journal Indexes: Science Citation Index Expanded (SCI-EXPANDED), Scopus, Academic Search Premier, PASCAL, Aerospace Database, Applied Science & Technology Source, Business Source Elite, Business Source Premier, Communication Abstracts, Compendex, Computer & Applied Sciences, INSPEC, Metadex, zbMATH, Civil Engineering Abstracts
  • Page Numbers: pp.869-881
  • Keywords: Discontinuous Galerkin time-domain method (DGTD), finite-element method (FEM), generalized sheet transition conditions (GSTCs), metasurface, numerical flux, time integration, time-domain analysis
  • Open Archive Collection: AVESIS Open Access Collection
  • Yıldız Technical University Affiliated: Yes

Abstract

The generalized sheet transition conditions (GSTCs) are incorporated into a discontinuous Galerkin time-domain (DGTD) method to efficiently simulate metasurfaces. The numerical flux for GSTCs is derived for the first time using the Rankine-Hugoniot jump conditions. This numerical flux includes the time derivatives of fields, and therefore, explicit time integration schemes, which are traditionally used with DGTD, do not yield a stable time marching method. To alleviate this bottleneck, a new time marching scheme, which solves a local matrix system for the unknowns of the elements touching the same GSTC face, is developed. This locally implicit method maintains its high-parallel efficiency just like the traditional explicit DGTD schemes. Numerical results, which validate the accuracy of the proposed method against analytical solutions and demonstrate its applicability to the simulation of curved and space/time-varying metasurfaces, are presented.