Multiple linear system techniques for 3D finite element method modeling of direct current resistivity

Abstract

The strategies that minimize the overall solution time of multiple linear systems in 3D finite element method (FEM) modeling of direct current (DC) resistivity were discussed. A global stiff matrix is assembled and stored in two parts separately. One part is associated with the volume integral and the other is associated with the subsurface boundary integral. The equivalent multiple linear systems with closer right-hand sides than the original systems were constructed. A recycling Krylov subspace technique was employed to solve the multiple linear systems. The solution of the seed system was used as an initial guess for the subsequent systems. The results of two numerical experiments show that the improved algorithm reduces the iterations and CPU time by almost 50%, compared with the classical preconditioned conjugate gradient method.

Keywords

finite element method modeling / direct current resistivity / multiple linear systems / preconditioned conjugate gradient / recycling Krylov subspace

References

[1] DeyA., MorrisonH. F.. Resistivity modeling for arbitrarily shaped three-dimensional structures [J]. Geophysics, 1979, 44(4): 753-780

[2] LiY., SpitzerK.. Three-dimensional DC resistivity forward modeling using finite elements in comparison with finite-difference solutions [J]. Geophys J Int, 2002, 151(3): 924-934

[3] RuanB. Y., XiongB., XuS. Z.. Finite element method for modeling resistivity sounding 3-D geo-electric section [J]. Earth Science-Journal of China University of Geosciences, 2001, 26(1): 73-77

[4] TangJ. T., RenZ. R., HuaX. R.. Triangle and tetrahedral finite element meshing from arbitrary geophysical model data [J]. Progress in Geophysics, 2006, 21(4): 1272-1280

[5] WuX. P., WangT. T.. A 3D finite element resistivity forward modeling using conjugate gradient algorithm [J]. Chinese Journal of Geophysics, 2003, 46(3): 428-432

[6] XuS. Z.The finite element method in geophysics [M], 1994, Beijing, Science Press: 173-183

[7] RückerC., GüntherT., SpitzerK.. Three-dimensional modeling and inversion of dc resistivity data incorporating topography: I. Modeling [J]. Geophys J Int, 2006, 166(2): 495-505

[8] BlomeM., MaurerH. R., SchmidtK.. Advances in three-dimensional geoelectric forward solver techniques [J]. Geophys J Int, 2009, 176(3): 740-752

[9] SAAD Y. Iterative methods for sparse linear systems [M]. 2nd Edition. Society for Industrial and Applied Mathematics, 2003.

[10] ÓlearyD. P.. The block conjugate gradient algorithm and related methods [J]. Linear Algebra Appl, 1980, 29: 293-322

[11] FreundR. W., MalhotraM.. A block QMR algorithm for non-Hermitian linear systems with multiple right-hand sides [J]. Linear Algebra Appl, 1997, 254(1): 119-157

[12] MorganR. B.. Restarted block-GMRES with deflation of eigenvalues [J]. Appl Numer Math, 2005, 54(2): 222-236

[13] ParksM. L., de SturlerE., MackeyG., et al.. Recycling Krylov subspaces for sequences of linear systems [J]. SIAM Journal on Scientific Computing, 2006, 28(5): 1651-1674

[14] GiraudL., RuizD., TouhamiA.. A comparative study of iterative solvers exploiting spectral information for SPD systems [J]. SIAM Journal on Scientific Computing, 2006, 27(5): 1760-1786

[15] SalkuyehD. K.. CG-type algorithms to solve symmetric matrix equations [J]. Applied Mathematics and Computation, 2006, 172(2): 985-999

[16] GolubG. H., RuizD., TouhamiA.. A hybrid approach combining Chebyshev filter and conjugate gradient for solving linear systems with multiple right-hand sides [J]. SIAM J Matrix Anal Appl, 2007, 29(3): 774-795

[17] ChanT. F., NgM. K.. Galerkin projection methods solving multiple linear systems [J]. SIAM Journal on Scientific Computing, 1999, 21(3): 836-850

[18] ErhelJ., GuyomarćhF.. An augmented conjugate gradient method for solving consecutive symmetric positive definite linear systems [J]. SIAM J Matrix Anal Appl, 2000, 21(4): 1279-1299

[19] WangS., de SturlerE., PaulinoG. H.. Large-scale topology optimization using preconditioned Krylov subspace methods with recycling [J]. Int J Numer Meth Engng, 2007, 69(12): 2441-2468

[20] JinC., CaiX. C.. A preconditioned recycling GMRES solver for Stochastic Helmholtz problems [J]. Commun Comput Phys, 2009, 6(2): 342-353

[21] ClemensM., HeliasM., SteinmetzT., WimmerG.. Multiple right-hand side techniques for the numerical simulation of quasistatic electric and magnetic fields [J]. Journal of Computational and Applied Mathematics, 2008, 215(2): 328-338

[22] MelloL. A. M., de SturlerE., PaulinoG. H., SilvaE. C.. Recycling Krylov subspaces for efficient large-scale electrical impedance tomography [J]. Comput Methods Appl Mech Engrg, 2010, 199(49): 3101-3110

[23] KilmerM. E., de SturlerE.. Recycling subspace information for diffuse optical tomography [J]. SIAM Journal on Scientific Computing, 2006, 27(6): 2140-2166