In this paper, an approach combining the DG method in space with CG method in time (CG-DG method) is developed to solve time-dependent Maxwell's equations when meta-materials are involved. Both the unconditional L2-stability and error estimate of order are obtained when polynomials of degree at most r is used for the temporal discretization and at most k for the spatial discretization. Numerical results in 3D are given to validate the theoretical results.

A fully discrete discontinuous Galerkin method is introduced for solving time-dependent Maxwell’s equations. Distinguished from the Runge-Kutta discontinuous Galerkin method (RKDG) and the finite element time domain method (FETD), in our scheme, discontinuous Galerkin methods are used to discretize not only the spatial domain but also the temporal domain. The proposed numerical scheme is proved to be unconditionally stable, and a convergent rate is established under the L2-norm when polynomials of degree atmost r and k are used for temporal and spatial approximation, respectively. Numerical results in both 2-D and 3-D are provided to validate the theoretical prediction. An ultra-convergence of order in time step is observed numerically for the numerical fluxes w.r.t. temporal variable at the grid points.