A PWDG method for the Maxwell system in anisotropic media with piecewise constant coefficient matrix

¹ College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, P.R. China
² LSEC, Institute of Computational Mathematics and Scientic/Engineering Computing, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, P.R. China
³ School of Mathematical Sciences, University of Chinese Academy of Sciences, Beijing 100049, P.R. China

^* Corresponding author: sdbjbjsd@163.com; hqy@lsec.cc.ac.cn

Received: 21 August 2023
Accepted: 27 November 2023

Abstract

In this paper we are concerned with plane wave discontinuous Galerkin (PWDG) methods for time-harmonic Maxwell equations in three-dimensional anisotropic media, for which the coefficients of the equations are piecewise constant symmetric matrices, where each constant symmetric matrix is defined on a medium (subdomain). By using suitable scaling transformations and coordinate (complex) transformations on every subdomain, the original Maxwell equation in anisotropic media is transformed into a Maxwell equation in isotropic media occupying a union domain of specific subdomains of complex Euclidean space. Based on these transformations, we define anisotropic plane wave basis functions and discretize the considered Maxwell equations by PWDG method with the proposed plane wave basis functions. We derive error estimates of the resulting approximate solutions, and further introduce a practically feasible local hp-refinement algorithm, which substantially improves accuracies of the approximate solutions. Numerical results indicate that the approximate solutions generated by the proposed PWDG methods possess high accuracy for the case of strong discontinuity media.