1887

Abstract

Summary

We introduce two novel preconditioning methods for 3D finite difference (FD) frequency domain electromagnetic modeling. The first preconditioner is based on the layered earth Green’s function, the second one is based on a special transformation of the original system of FD equations into a system with a contraction operator. This transformation extends to the FD modeling the approach originally developed for the integral equation modeling method. We also examine the spectral properties of the two preconditioned methods.

For numerical study of the developed methods, we have designed a complex marine 3D geoelectrical model. The numerical experiments confirm the results of theoretical analysis: an iterative solver with the contraction preconditioner converges faster or at near the same speed as with the Green’s function preconditioner. We apply these two solvers to controlled source modeling and demonstrate that both approaches are very effective, and memory saving. We have also developed a parallel version of the algorithm and studied the scalability of shared and distributed memory parallelization.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.201701354
2017-06-12
2024-03-28
Loading full text...

Full text loading...

References

  1. Martikainen, J., Rossi, T., Rogovin, K. & Toivanen, J.
    [2003] Fast direct solver for a time-harmonic electromagnetic problem with an application, in Mathematical and Numerical Aspects of Wave Propagation, WAVES 2003, pp. 675–680, Springer.
    [Google Scholar]
  2. Weiss, C. J. and Newman, G.A.
    [2002] Electromagnetic induction in a fully 3-D anisotropic earth. Geophysics, 67, 1104–1114.
    [Google Scholar]
  3. Yavich, N. and Zhdanov, M.S.
    [2016], Contraction pre-conditioner in finite-difference electromagnetic modeling, Geophys. J. Int, 206 (3): 1718–1729.
    [Google Scholar]
  4. Yee, K.S.
    [1966] Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media. IEEE Trans. Antennas and Propagation, 16, 302–307.
    [Google Scholar]
  5. Zaslavsky, M., Druskin, V., Davydycheva, S., Knizhnerman, L., Abubakar, A. & Habashy, T.
    [2011], Hybrid finite-difference integral equation solver for 3D frequency domain anisotropic electromagnetic problems, Geophysics, 76(2), F123–F137.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.201701354
Loading
/content/papers/10.3997/2214-4609.201701354
Loading

Data & Media loading...

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