High order parallel discontinuous Galerkin time domain method with curvilinear element and continuously varying material properties for Maxwell's equations

Discontinuous Galerkin (DG) finite element method is well-suited on unstructured meshes for high order approximation with the freedom of choosing the order of basis functions in each element locally. Moreover, DG can handle complicated geometries with curved boundary easily. Further, since information exchange in DG only involves neighboring elements, high efficiency in parallelization can be achieved. The conformal perfect matched layer (PML) with continuously changing material property tensors in space is adopted to truncate the computational domain. And the model presented in this paper treats the material property tensor in PML region with polynomial representation using a universal array approach.