▎ 摘　　要

In this paper, we develop a new variational form to simulate the propagation of surface plasmon polaritons on graphene sheets. Here the graphene is treated as a thin sheet of current with an effective conductivity, and modeled as a lower-dimensional interface. A novel time-domain finite element method is proposed for solving this graphene model, which coupled an ordinary differential equation on the interface with Maxwell's equations in the physical domain. Discrete stability and error estimate are proved for our proposed method. Numerical results are presented to demonstrate the effectiveness of this graphene model for simulating the surface plasmon polaritons propagating on graphene sheets.