The effect of the 2-D Laplacian operator approximation on the performance of finite-difference time-domain schemes for Maxwell's equations

Anastasios H. Panaretos, James Aberle, Rodolfo Diaz

Research output: Contribution to journalArticlepeer-review

13 Scopus citations

Abstract

The behavior of the finite-difference time-domain method (FDTD) is investigated with respect to the approximation of the two-dimensional Laplacian, associated with the curl-curl operator. Our analysis begins from the observation that in a two-dimensional space the Yee algorithm approximates the Laplacian operator via a strongly anisotropic 5-point approximation. It is demonstrated that with the aid of a transversely extended-curl operator any 9-point Laplacian can be mapped onto FDTD update equations. Our analysis shows that the mapping of an isotropic Laplacian approximation results in an isotropic and less dispersive FDTD scheme. The properties of the extended curl are further explored and it is proved that a unity Courant number can be achieved without the resulting scheme suffering from grid decoupling. Additionally, the case of a 25-point isotropic Laplacian is examined and it is shown that the corresponding scheme is fourth order accurate in space and exhibits isotropy up to sixth order. Representative numerical simulations are performed that validate the theoretically derived results.

Original languageEnglish (US)
Pages (from-to)513-536
Number of pages24
JournalJournal of Computational Physics
Volume227
Issue number1
DOIs
StatePublished - Nov 10 2007

Keywords

  • Curl operator
  • Electromagnetics
  • Finite-difference time-domain method
  • Laplacian operator

ASJC Scopus subject areas

  • Numerical Analysis
  • Modeling and Simulation
  • Physics and Astronomy (miscellaneous)
  • Physics and Astronomy(all)
  • Computer Science Applications
  • Computational Mathematics
  • Applied Mathematics

Fingerprint Dive into the research topics of 'The effect of the 2-D Laplacian operator approximation on the performance of finite-difference time-domain schemes for Maxwell's equations'. Together they form a unique fingerprint.

Cite this