Mixed Finite-Element Time-Domain Method for Transient Maxwell Equations in Doubly Dispersive Media

We describe a mixed finite-element time-domain algorithm to solve transient Maxwell equations in inhomogeneous and doubly dispersive linear media where both the permittivity and permeability are functions of frequency. The mixed finite-element time-domain algorithm is based on the simultaneous use of both electric and magnetic field as state variables with a mix of edge (Whitney 1-form) and face (Whitney 2-form) elements for discretization of the coupled first-order Maxwell curl equations. The constitutive relations are decoupled from the curl equations and cast in terms of (auxiliary) ordinary differential equations involving time derivatives. Permittivity and permeability dispersion models considered here are quite general and recover Lorentz, Debye, and Drude models as special cases. The present finite-element time-domain algorithm also incorporates the perfectly matched layer absorbing boundary conditions in a natural way.