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

Creative Commons License

Chen L., ÖZAKIN M. B., Ahmed S., Bagci H.

IEEE Transactions on Antennas and Propagation, cilt.69, sa.6, ss.3605-3610, 2021 (SCI-Expanded) identifier identifier

  • Yayın Türü: Makale / Tam Makale
  • Cilt numarası: 69 Sayı: 6
  • Basım Tarihi: 2021
  • Doi Numarası: 10.1109/tap.2020.3037651
  • Dergi Adı: IEEE Transactions on Antennas and Propagation
  • Derginin Tarandığı İndeksler: 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
  • Sayfa Sayıları: ss.3605-3610
  • Anahtar Kelimeler: Absorbing boundary conditions, discontinuous Galerkin (DG) method, perfectly matched layer (PML), time-domain analysis, weight-adjusted approximation (WAA)
  • Yıldız Teknik Üniversitesi Adresli: Evet

Özet

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.