Abstract
We derive and analyze an efficient algorithm to incorporate the anomalously dispersive Havriliak-Negami dielectric model of induced polarization in the Finite-difference time-domain (FD-TD) method. Our algorithm implements this dielectric model, which in the time-domain involves fractional derivatives and fractional differential operators, with a preset error over the desired computational time interval [0,Tcomp] and correctly takes into account the singularity at t=0+ of the corresponding time-domain dielectric susceptibility. The overall algorithm is shown to be second-order accurate in space and time, and to obey the standard FD-TD stability condition. Numerical experiments confirm our analysis.
Original language | English (US) |
---|---|
Pages (from-to) | 3884-3899 |
Number of pages | 16 |
Journal | Journal of Computational Physics |
Volume | 230 |
Issue number | 10 |
DOIs | |
State | Published - May 10 2011 |
All Science Journal Classification (ASJC) codes
- Numerical Analysis
- Modeling and Simulation
- Physics and Astronomy (miscellaneous)
- General Physics and Astronomy
- Computer Science Applications
- Computational Mathematics
- Applied Mathematics
Keywords
- Computational electromagnetics
- Dispersive dielectrics
- Finite-difference time-domain
- Time-domain