▎ 摘　　要

In this paper, a simple finite difference time domain (FDTD) approach is presented for the analysis and design of graphene based optical devices. Here, Maxwell's curl equations are expressed in normalized flux density. Then, the FDTD method is incorporated with that expression. As the flux density is considered, graphene conductivity which depends on frequency is easily incorporated with the method. The proposed method is used to investigate properties of single and multilayer graphene sheets. To probe the transmission, reflection and surface plasmon polariton (SPP) mode characteristics of graphene in terahertz and far infared frequency regime, numerical examples are demonstrated. Validation of the method has been demonostrated comparing with existing methods and results available in the literature. The accuracy, stability and benefits of proposed method are also presented. (C) 2017 Elsevier GmbH. All rights reserved.