“…Computational analysis of fields in increasingly intricate periodic unit cells plays a significant role in their design and optimization. In the frequency domain, Integral Equation (IE) [5], [6], Finite Element (FE) [7], [8], and Discontinuous Galerkin (DG) [9] methods have been successfully applied to a variety of periodic electromagnetic systems. Timedomain (TD) methods for studying periodic systems include FE [10], [11], IE [12], and Finite Difference Time Domain (FDTD) [13], while DG methods remain relatively unexplored.…”