Three‐dimensional electromagnetic modeling using finite difference equations: The magnetotelluric example

DOI Web Site Web Site 28 Citations

Abstract

<jats:p>We have developed a robust and efficient finite difference algorithm for computing the magnetotelluric response of general three‐dimensional (3‐D) models using the minimum residual relaxation method. The difference equations that we solve are second order in <jats:bold>H</jats:bold> and are derived from the integral forms of Maxwell's equations on a staggered grid. The boundary <jats:bold>H</jats:bold> field values are obtained from two‐dimensional transverse magnetic mode calculations for the vertical planes in the 3‐D model. An incomplete Cholesky decomposition of the diagonal subblocks of the coefficient matrix is used as a preconditioner, and corrections are made to the <jats:bold>H</jats:bold> fields every few iterations to ensure there are no <jats:bold>H</jats:bold> divergences in the solution. For a plane wave source field, this algorithm reduces the errors in the <jats:bold>H</jats:bold> field for simple 3‐D models to around the 0.01% level compared to their fully converged values in a modest number of iterations, taking only a few minutes of computation time on our desktop workstation. The <jats:bold>E</jats:bold> fields can then be determined from discretized versions of the curl of <jats:bold>H</jats:bold> equations.</jats:p>

Journal

  • Radio Science

    Radio Science 29 (4), 923-935, 1994-07

    American Geophysical Union (AGU)

Citations (28)*help

See more

Keywords

Details 詳細情報について

Report a problem

Back to top