Developing and analyzing a DG method for modeling of metasurfaces

Metasurfaces, with their unique ability to manipulate electromagnetic waves at subwavelength scales, have garnered significant attention for their applications in various areas such as wavefront shaping, beam steering, and cloaking. Accurately modeling these structures presents a complex challenge due to the intricate interactions between electromagnetic waves and the metasurface's subwavelength geometry. In this work, we present a discontinuous Galerkin time-domain (DGTD) method for the efficient and accurate modeling of metasurfaces by coupling with the generalized sheet transition conditions (GSTCs). Numerical stability is proved for two types of metasurfaces. Through a series of benchmark numerical simulations, we demonstrate the accuracy and computational efficiency of the DGTD method in modeling metasurfaces, especially those with irregular geometries. This work opens new avenues for the accurate simulation and design of metasurfaces in a wide range of electromagnetic applications.