Algebraic multigrid preconditioning within parallel finite-element solvers for 3-D electromagnetic modelling problems in geophysics
MetadataShow full item record
We present an elaborate preconditioning scheme for Krylov subspace methods which has been developed to improve the performance and reduce the execution time of parallel node-based finite-element (FE) solvers for 3-D electromagnetic (EM) numerical modelling in exploration geophysics. This new preconditioner is based on algebraic multigrid (AMG) that uses different basic relaxation methods, such as Jacobi, symmetric successive over-relaxation (SSOR) and Gauss-Seidel, as smoothers and the wave front algorithm to create groups, which are used for a coarse-level generation. We have implemented and tested this new preconditioner within our parallel nodal FE solver for 3-D forward problems in EM induction geophysics. We have performed series of experiments for several models with different conductivity structures and characteristics to test the performance of our AMG preconditioning technique when combined with biconjugate gradient stabilized method. The results have shown that, the more challenging the problem is in terms of conductivity contrasts, ratio between the sizes of grid elements and/or frequency, the more benefit is obtained by using this preconditioner.Compared to other preconditioning schemes, such as diagonal, SSOR and truncated approximate inverse, the AMG preconditioner greatly improves the convergence of the iterative solver for all tested models. Also, when it comes to cases in which other preconditioners succeed to converge to a desired precision, AMG is able to considerably reduce the total execution time of the forward-problem code-up to an order of magnitude. Furthermore, the tests have confirmed that ourAMGscheme ensures grid-independent rate of convergence, aswell as improvement in convergence regardless of howbig local mesh refinements are. In addition, AMGis designed to be a black-box preconditioner, which makes it easy to use and combine with different iterative methods. Finally, it has proved to be very practical and efficient in the parallel context.
Showing items related by title, author, creator and subject.
Algebraic multigrid preconditioning for finite-element methods for 3-D electromagnetic modelling in geophysicsKoldan, J.; Puzyrev, Volodymyr; Cela, J. (2013)Efficiency of a 3-D electromagnetic numerical modelling scheme is critical for its future use within a 3-D electromagnetic inversion algorithm. Therefore, we have developed and implemented a more elaborate preconditioning ...
Puzyrev, Volodymyr; Cela, J. (2015)Practical applications of controlled-source electromagnetic (EM) modelling require solutions for multiple sources at several frequencies, thus leading to a dramatic increase of the computational cost. In this paper, we ...
Puzyrev, Vladimir; De La Puente, J.; Cela, J. (2013)We present a massively parallel finite-element method for 3D electromagnetic forward modelling using vector and scalar potentials. The numerical scheme employs parallel Krylov subspace solvers with efficient preconditioning ...