1887
Volume 31, Issue 1-2
  • ISSN: 0812-3985
  • E-ISSN: 1834-7533

Abstract

Electromagnetic modelling of 3D conductive targets using staggered grids results in matrix systems that are severely ill conditioned and unsuitable for many iterative solvers. Either a minimal residual solver, MINRES, or a restarted biconjugate gradient stabilised solver, BiCGSTAB, can be used to bring the solution to convergence. The BiCGSTAB, method, combining the bi-conjugate gradient method with a restarted minimal residual correction, is very efficient in solving the matrix system arising in our applications. To speed up the solver further we employ various preconditioning techniques. The use of a symmetric Jacobian preconditioner greatly improves the convergence of the MINRES method. In homogeneous regions where the divergence of the electric field is zero the Maxwell equations reduce to vector Helmholtz equations with decoupled components. These Helmholtz equations can be readily approximated by central differences, resulting in diagonally dominant matrices that can be used as preconditioners for the original curl curl equation of the staggered grid. The diagonals of the Helmholtz equation can also be used as a Jacobian preconditioner, which is superior to Jacobian preconditioning using the diagonal elements of the original curl curl equation. Numerical results have demonstrated the efficiency of these approaches.

© 2000 ASEG
Loading

Article metrics loading...

/content/journals/10.1071/EG00158
2000-03-01
2026-01-19

  • From This Site

    /content/journals/10.1071/EG00158
    dcterms_title,dcterms_subject,pub_keyword
    -contentType:Journal -contentType:Contributor -contentType:Concept -contentType:Institution
    10
    5
Loading full text...

Full text loading...

References

  1. Agarwal, A. K. and J. T. Weaver, 1998, Comparison of electric fields obtained from 3D magnetic field calculations on fixed staggered grids, The 14th Workshop on Electromagnetic Induction in the Earth, Sinaia, Romania, Aug. 16-22, 1998.
  2. Arnason, K., 1998, A consistent discretization of the electromagnetic field in conducting media and application to the TEM problem, in Oristaglio, M. and Spies, B., Eds., 3D Electromagnetics, Society of Exploration Geophysicists, 1998.
  3. Druskin,V. and Knizhnerman, L., 1988, A spectral semi-discrete method for numerical solution of 3-D non-stationary problems in electrical prospecting: Phys. Sol. Earth, 24, 641-648.
  4. Druskin, V. and Knizhnerman, L., 1994, Spectral approach to solving threedimensional Maxwell’s diffusion equations in the time and frequency domain: Radio Science, 29, 937-953.
  5. Gupta, P.K., Raiche, A.P. and Sugeng, F., 1989, Three-dimensional time-domain electromagnetic modeling using a compact finite-element frequency-stepping method: Geophys. J., 96, 457-468.
  6. Hohmann, G. W., 1975, Three-dimensional induced polarization and electromagnetic modeling: Geophysics, 40, 309-324.
  7. Mackie, R. L., Smith, J.T., and Madden, T.R., 1994, Three-dimensional electromagnetic modeling using finite difference equations, The magnetotelluric example: Radio Science, 29, 923-935.
  8. Paige, C.C. and Saunders, M.A., 1975, Solution of sparse indefinite systems of linear equations: SIAM J. Numer. Anal., 12, 617-629.
  9. Raiche, A. P., 1974, An integral equation approach to three-dimensional modeling: Geophys. J. R. Astr. Soc., 36, 363-376.
  10. Raiche, A.P. and Sugeng, F., 1989, Predicting the transient response of complex structures using the compact finite element method: Expl. Geophys., 20, 51-56.
  11. Sleijpen, G.L.G. and Fokkema, D.R., 1993, BiCGSTAB(l) for linear matrices involving nonsymmetric matrices with complex spectrum: ETNA, 1, 11-32.
  12. Sleijpen, G.L.G. and van der Horst, H.A., 1993, Krylov subspace methods for large linear systems of equations, Mathematical Institute, Univ. of Utrecht.
  13. Smith, J.T., 1996, Conservative modeling of 3-D electromagnetic fields; Part I: Properties and error analysis: Geophysics, 61, 1308-1318.
  14. Sugeng, F., 1998, Modelling the 3D TDEM response using the 3D full-domain finite element method based on the hexahedral edge-element technique: Expl. Geophys., 29, 615-619.
  15. Sugeng, F., Raiche, A. and Rijo, L., 1993, Comparing the time-domain response of 2- D and elongated 3-D conductors excited by a rectangular loop source: J. Geomag. Geoelectr, 45, 873-875.
  16. Wang, T., and Hohmann, G.W., 1993, A finite-difference, time-domain solution for three-dimensional electromagnetic modeling: Geophysics, 58, 797-809.
  17. Weaver, J.T., Agarwal, A. K. and Pu, X.H., 1993, Recent developments in three dimensional finite difference modelling of the magnetic field in geoelectromagnetic induction, International Symposium on Three-dimensional Electromagnetics, Schlumberger-Doll Research, Ridgefield, Connecticut, USA, Oct. 4-6, 1993.
  18. Weidelt, P., 1975, Electromagnetic induction in three-dimensional structures: J. Geophys., 41, 85-109.
  19. Weidelt, P., 1998, Three-dimensional conductivity models: implications of electrical anisotropy, in Oristaglio, M. and Spies, B., Eds., 3D Electromagnetics, Society of Exploration Geophysicists, 1998.
  20. Wesseling, P., 1992, An introduction to multigrid methods: John Wiley & Sons Ltd.
  21. Xiong, Z., 1989a, Electromagnetic fields of electrical dipoles embedded in a stratified anisotropic earth, Geophysics, 54, 1643-1646.
  22. Xiong, Z., 1989b, Electromagnetic fields of electrical dipoles embedded in a stratified anisotropic earth - errata, Geophysics, 56, 414.
  23. Xiong, Z., 1999, Domain decomposition for 3D electromagnetic modeling: Earth, Planets and Space, 51, 1013-1018.
  24. Yee, K. S., 1966, Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media: IEEE Trans. Antennas Propagat., AP-14, 302-307.
/content/journals/10.1071/EG00158
Loading
  • Article Type: Research Article
Keyword(s): Full-domain 3D electromagnetic modelling; preconditioning

Most Cited This Month Most Cited RSS feed

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error