A new discontinuous Galerkin spectral element timedomain (DG-SETD) method for Maxwell's equations based on the field variables E and B is proposed to analyze three-dimensional (3-D) transient electromagnetic phenomena. Compared to the previous SETD method based on the field variables E and H (the EH scheme), in which different orders of interpolation polynomials for electric and magnetic field intensities are required, the newly proposed method can eliminate spurious modes using basis functions with the same order interpolation for electric field intensity and magnetic flux density (the EB scheme). Consequently, it can reduce the number of unknowns and computation load. Domain decomposition for the EB scheme SETD method is completed via the DG method. In addition, the EB scheme SETD method is extended to the well-posed time-domain perfectly matched layer (PML) to truncate the computation domain when solving open-region problems. The effectiveness and advantages of the new DG-SETD method are validated by eigenvalue analysis and numerical results. Index Terms-Discontinuous Galerkin (DG) method, EB scheme, Maxwell's equations, perfectly matched layer (PML), Riemann solver, spectral element time-domain (SETD) method. I. INTRODUCTION I N computational electromagnetics, finite element timedomain (FETD) method in full-wave simulation has received much attention due to the ability to handle transient and nonlinear phenomena in complex structures [1]-[8]. It has been widely used in IC packaging, smart antennas, electromagnetic compatibility (EMC), bioelectromagnetics, etc. The spectral element time-domain (SETD) method can be viewed as a special case of the FETD method [9]-[13]. While FETD usually employs low-order tetrahedrons to mesh, the SETD method uses hexahedron elements to discretize the computation region, and the interpolation polynomials in the hexahedrons can be either of low order or high order. The SETD method is Manuscript