Three-dimensional transient electromagnetic modeling using rational Krylov methods

A computational method is given for solving the forward modeling problem for transient electromagnetic exploration. Its key features are discretization of the quasi-static Maxwell's equations in space using the first-kind family of curl-conforming Nedelec elements combined with time integration using rational Krylov subspace methods. We show how rational Krylov subspace methods may be used to solve the same problem in the frequency domain followed by a synthesis of the transient solution using the fast Hankel transform, arguing that the pure time-domain is more efficient. We also propose a simple method for selecting the pole parameters of the rational Krylov subspace method which leads to convergence within an a priori determined number of iterations independent of mesh size and conductivity structure. These poles are repeated in a cyclic fashion, which, in combination with direct solvers for the discrete problem, results in significantly faster solution times than previously proposed schemes.