Modeling of Metasurfaces Using Discontinuous Galerkin Time-Domain Method Based on Generalized Sheet Transition Conditions

A new discontinuous Galerkin time-domain (DGTD) method based on generalized sheet transition conditions (GSTCs) for analyzing the transient response of electromagnetic metasurfaces is presented. The discontinuities in electromagnetic fields around the vicinity of the metasurface are modeled by placing virtual edges on both sides of the metasurface. The GSTCs are incorporated in the updating equations of the standard DGTD algorithm by integrating the time-domain GSTCs formula with Galerkin' s-weighted residual method. The DGTD-GSTC formulation can also simulate a bf piecewise linear approximation of a curved metasurfaces due to the unstructured mesh. This is the first time that DGTD is applied to the simulation of GSTCs. The implementation of this new method is verified by numerical cases.