In this paper, a three-dimensional Petrov-Galerkin finite element interface (3D PGFEI) method with non-body-fitted grids is proposed for solving the inhomogeneous anisotropic time-harmonic Maxwell's equations in irregular regions. We derive the weak variational form of the anisotropic Maxwell interface problem. The irregular region is discretized by uniform tetrahedral grids while the media interface is not aligned with the grids. The level-set functions are applied to capture the discontinuity on the anisotropic media interface. Meanwhile, we construct the special basis functions according to the jump conditions across the inhomogeneous interface, where the permittivity and permeability are discontinuous three-dimensional tensors. The resulting linear system is a banded sparse matrix and demands fewer memory resources than dense matrix. Numerical experiments illustrate the efficiency and accuracy of the proposed 3D PGFEI method and the solutions can achieve second-order accuracy in terms of the relative L∞(Ω) and relative L2(Ω) error norms.
Read full abstract