▎ 摘　　要

This letter extends wave equation-based discontinuous Galerkin time-domain (DGTD-WE) method to analyze the graphene-based devices. The numerical fluxes related to a general impedance transmission boundary condition in the DGTD-WE method have been proposed for the first time to model the graphene sheet as an infinite thin impedance surface. A vector-fitting technique is used to approximate surface conductivity of the graphene in terms of a series of partial fractions. Several graphene-based devices have been given to demonstrate the correctness and effectiveness of the proposed method.