A spectral numerical method for electromagnetic diffusion
- 1 January 2006
- journal article
- Published by Society of Exploration Geophysicists in Geophysics
- Vol. 71 (1) , I1-I9
Abstract
I present a pseudospectral explicit scheme that can simulate low-frequency electromagnetic (EM) propagation in the earth. This scheme solves linear periodic parabolic equations, having accuracy within machine precision, both temporally and spatially. The method is based on a Chebyshev expansion of the evolution operator, with the spatial derivatives computed via a staggered Fourier pseudospectral technique. The results match analytical solutions of the initial-value problem and the Green's function. An example of the EM field produced by a set of magnetic sources in a heterogeneous model illustrates the algorithm's performance.Keywords
This publication has 27 references indexed in Scilit:
- Solving the Maxwell equations by the Chebyshev method: a one-step finite-difference time-domain algorithmIEEE Transactions on Antennas and Propagation, 2003
- Electric drill-string telemetryJournal of Computational Physics, 2003
- Parallel finite element algorithm with domain decomposition for three-dimensional magnetotelluric modellingJournal of Applied Geophysics, 2000
- Staggered mesh for the anisotropic and viscoelastic wave equationGeophysics, 1999
- A numerical evaluation of electromagnetic methods in geothermal explorationGeophysics, 1996
- Modeling elastic fields across irregular boundariesGeophysics, 1992
- Efficient Solution of Parabolic Equations by Krylov Approximation MethodsSIAM Journal on Scientific and Statistical Computing, 1992
- Spectral Methods in Time for Parabolic ProblemsSIAM Journal on Numerical Analysis, 1989
- Implementation of a prime factor FFT algorithm on CRAY-1Parallel Computing, 1988
- Absorbing boundaries for wave propagation problemsJournal of Computational Physics, 1986