Journal
IMA JOURNAL OF NUMERICAL ANALYSIS
Volume 32, Issue 3, Pages 926-956Publisher
OXFORD UNIV PRESS
DOI: 10.1093/imanum/drr001
Keywords
Maxwell's equations; Debye; Lorentz; higher-order FDTD; stability; dispersion
Categories
Funding
- National Science Foundation (NSF) [DMS-0811223]
- NSF's Collaborations in Mathematical Geosciences program [CMG EAR-0724865]
- Direct For Mathematical & Physical Scien
- Division Of Mathematical Sciences [0811223] Funding Source: National Science Foundation
Ask authors/readers for more resources
We study the stability properties of, and the phase error present in, several higher-order (in space) staggered finite difference schemes for Maxwell's equations coupled with a Debye or Lorentz polarization model. We present a novel expansion of the finite difference approximations, of arbitrary (even) order, of the first-order spatial derivative operator. This alternative representation allows the derivation of a concise formula for the numerical dispersion relation for all even order schemes applied to each model, including the limiting (infinite-order) case. We further derive a closed-form analytical stability condition for these schemes as a function of the order of the method. Using representative numerical values for the physical parameters, we validate the stability criterion while quantifying numerical dissipation. Lastly, we demonstrate the effect that the spatial discretization order, and the corresponding stability constraint, has on the dispersion error.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available