Open AccessCombining standard with optimised split‐step finite‐difference time‐domain methods for the study of graphene configurations
Author(s) -
Amanatiadis Stamatios A.,
Zygiridis Theodoros T.,
Kantartzis Nikolaos V.
Publication year - 2019
Publication title -
iet science, measurement and technology
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.418
H-Index - 49
eISSN - 1751-8830
pISSN - 1751-8822
DOI - 10.1049/iet-smt.2018.5641
Subject(s) - graphene , discretization , robustness (evolution) , finite difference time domain method , finite difference method , time domain , computer science , finite difference , algorithm , partial differential equation , finite element method , ordinary differential equation , mathematics , mathematical optimization , differential equation , mathematical analysis , materials science , physics , computer vision , nanotechnology , biochemistry , chemistry , quantum mechanics , gene , thermodynamics
In this study, the authors propose a hybrid computational model for the reliable simulation of electromagnetic‐wave phenomena emerging in graphene structures, which incorporates two variations of an unconditionally‐stable finite‐difference time‐domain algorithm. The new approach features: (a) a dispersive model that relies on the auxiliary differential equation technique and takes into consideration the graphene's surface conductivity, and (b) a modified discretisation algorithm implementing error‐optimised spatial approximations, which is appropriate for graphene‐free mesh nodes. The resulting update equations are devoid of time‐step restrictions, thus formulating a suitable and useful computational framework for the investigation of contemporary graphene‐based problems, which commonly call for densely‐sampled grids and/or prolonged simulations. After theoretically assessing the fundamental features of the proposed algorithm, the propagation properties of surface waves on graphene are extracted numerically and compared with analytical estimations, so that computational validation is provided. In addition, a realistic surface‐wave coupler configuration is simulated and its main characteristics are extracted accurately, exemplifying the method's robustness on more complex setups.