Advances in Numerical Analysis
Volume 2010 (2010), Article ID 923832, 28 pages
doi:10.1155/2010/923832
Research Article

Finite Element Approximation of Maxwell's Equations with Debye Memory

BICOM, Institute of Computational Mathematics, Brunel University, Uxbridge, Middlesex UB8 3PH, UK

Received 18 May 2010; Accepted 31 August 2010

Academic Editor: Dongwoo Sheen

Copyright © 2010 Simon Shaw. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Maxwell's equations in a bounded Debye medium are formulated in terms of the standard partial differential equations of electromagnetism with a Volterra-type history dependence of the polarization on the electric field intensity. This leads to Maxwell's equations with memory. We make a correspondence between this type of constitutive law and the hereditary integral constitutive laws from linear viscoelasticity, and we are then able to apply known results from viscoelasticity theory to this Maxwell system. In particular, we can show long-time stability by shunning Gronwall's lemma and estimating the history kernels more carefully by appeal to the underlying physical fading memory. We also give a fully discrete scheme for the electric field wave equation and derive stability bounds which are exactly analogous to those for the continuous problem, thus providing a foundation for long-time numerical integration. We finish by also providing error bounds for which the constant grows, at worst, linearly in time (excluding the time dependence in the norms of the exact solution). Although the first (mixed) finite element error analysis for the Debye problem was given by Li (2007), this seems to be the first time sharp constants have been given for this problem.