Solution of time-convolutionary Maxwell's equations using parameter-dependent Krylov subspace reduction

We suggest a new algorithm for the solution of the time domain Maxwell equations in dispersive media. After spacial discretization we obtain a large system of time-convolution equations. Then this system is projected onto a small subspace consisting of the Laplace domain solutions for a preselected set of Laplace parameters. This approach is a generalization of the rational Krylov subspace approach for the solution of non-dispersive Maxwell's systems. We show that the projected system preserves such properties of the initial system as stability and passivity. As an example we consider the 3D quasistationary induced polarization problem with the Cole-Cole conductivity model important for geophysical oil exploration. Our numerical experiments show that the introduction of the induced polarization does not have significant effect on convergence. (C) 2010 Elsevier Inc. All rights reserved.