Multiple linear system techniques for 3D finite element method modeling of direct current resistivity |
| |
Authors: | Chang-wei Li Bin Xiong Jian-ke Qiang Yu-zeng Lü |
| |
Affiliation: | [1]School of Geoscience and Info-physics, Central South University, Changsha 410083, China [2]College of Earth Sciences, Guilin University of Technology, Guilin 541004, China |
| |
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 |
本文献已被 维普 SpringerLink 等数据库收录! |