The discontinuous Galerkin finite-element time-domain (DG-FETD) method is useful in transient simulations of multiscale electromagnetic systems. Its special capability in geometric modeling, by dividing the computational domain in several domains, transforms a large system into several moderate-sized matrix equations. The numerical fluxes, which communicate fields between domains, are defined by tangential components of E and H on the interfaces; for this reason, the conventional DG-FETD is traditionally based on these two variables.