In this paper, the electromagnetic (EM) features of graphene are characterized by a discontinuous Galerkin time-domain (DGTD) algorithm with a resistive boundary condition (RBC). The atomically thick graphene is equivalently modeled using an RBC by regarding the graphene as an infinitesimally thin conductive sheet. To incorporate RBC into the DGTD analysis, the surface conductivity of the graphene composed of contributions from both intraband and interband terms is first approximated by rational basis functions using the fast-relaxation vector-fitting (FRVF) method in the Laplace domain. Next, through the inverse Laplace transform, the corresponding time-domain matrix equations in integral can be obtained. Finally, these matrix equations are solved by time-domain finite integral technique (FIT). For elements not touching the graphene sheet, however, the well-known Runge-Kutta (RK) method is employed to solve the two first-order time-derivative Maxwell's equations. The application of the surface boundary condition significantly alleviates the memory consuming and the limitation of time step size required by Courant-Friedrichs-Lewy (CFL) condition. To validate the proposed algorithm, various numerical examples are presented and compared with available references.


Electrical and Computer Engineering

Keywords and Phrases

Discontinuous Galerkin time-domain (DGTD) method; fast-relaxation vector-fitting (FRVF); finite integral technique (FIT); grapheme; Laplace transform; resistive boundary condition (RBC); surface conductivity

International Standard Serial Number (ISSN)


Document Type

Article - Journal

Document Version


File Type





© 2024 Institute of Electrical and Electronics Engineers, All rights reserved.

Publication Date

01 Jul 2015