Efficient and stable generalized auxiliary differential equation FDTD implementation of graphene dispersion. Issue 6 (27th September 2019)

• Record Type: Journal Article
• Authors: Ramadan, Omar
• Is Part Of: Compel. Volume 38:Issue 6(2019)
• Page Start: 2070
• Page End: 2083
• Publication Date: 2019-09-27
• Subjects: Finite difference time-domain analysis -- Computational electromagnetics -- Graphene -- Auxiliary differential equation (ADE) -- Generalized dispersive model (GDM) -- Stability analysis -- Von Neumann -- Routh-Hurwitz; Electrical engineering -- Data Processing -- Periodicals; Electrical engineering -- Mathematics -- Periodicals; Electrical engineering -- Periodicals; Electronics -- Data Processing -- Periodicals; Electronics -- Mathematics -- Periodicals; 621.3
• DOI: 10.1108/COMPEL-03-2019-0115
• Languages: English
• ISSNs: 0332-1649

Abstract

Abstract : Purpose: The purpose of this paper is to present efficient and stable generalized auxiliary differential equation finite difference time domain (G-ADE-FDTD) implementation of graphene dispersion. Design/methodology/approach: A generalized dispersive model is used for describing the graphene's intraband and interband conductivities in the terahertz and infrared frequencies. In addition, the von Neumann method combined with the Routh-Hurwitz criterion are used for studying the stability of the given implementation. Findings: The presented G-ADE-FDTD implementation allows modeling graphene's dispersion using the minimal number of additional auxiliary variables, which will reduce both the CPU time and memory storage requirements. In addition, the stability of the implementation retains the standard non-dispersive Courant–Friedrichs–Lewy (CFL) constraint. Practical implications: The given implementation is conveniently applicable for most commonly used dispersive models, such as Debye, Lorentz, complex-conjugate pole residue, etc. Originality/value: The presented G-ADE-FDTD implementation not only unifies the implementation of both graphene's intraband and interband conductivities, with the minimal computational requirements but also retains the standard non-dispersive CFL time step stability constraint.