×

Multi-level incomplete factorizations for the iterative solution of nonlinear finite element problems. (English) Zbl 1176.74184

Summary: Several engineering applications give rise quite naturally to linearized FE systems of equations possessing a multi-level structure. An example is provided by geomechanical models of layered and faulted geological formations. For such problems the use of a multi-level incomplete factorization (MIF) as a preconditioner for Krylov subspace methods can prove a robust and efficient solution accelerator, allowing for a fine tuning of the fill-in degree with a significant improvement in both the solver performance and the memory consumption. The present paper develops two novel MIF variants for the solution of multi-level symmetric positive definite systems. Two correction algorithms are proposed with the aim of preserving the positive definiteness of the preconditioner, thus avoiding possible breakdowns of the preconditioned conjugate gradient solver. The MIF variants are experimented with in the solution of both a single system and a long-term quasi-static simulation dealing with a multi-level geomechanical application. The numerical results show that MIF typically outperforms by up to a factor 3 a more traditional algebraic preconditioner such as an incomplete Cholesky factorization with partial fill-in. The advantage is emphasized in a long-term simulation where the fine fill-in tuning allowed for by MIF yields a significant improvement for the computer memory requirement as well.

MSC:

74S05 Finite element methods applied to problems in solid mechanics
74L05 Geophysical solid mechanics
74F10 Fluid-solid interactions (including aero- and hydro-elasticity, porosity, etc.)

Software:

BILUTM; ILUT; ILUM; ARMS
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Saad, Iterative solution of linear systems in the 20th century, Journal of Computational and Applied Mathematics 23 pp 1– (2000) · Zbl 0965.65051
[2] van der Vorst, Iterative Krylov Methods for Large Linear Systems (2003) · doi:10.1017/CBO9780511615115
[3] Janna, A comparison of projective and direct solvers for Finite Elements in elastostatics, Advances in Engineering Software 40 pp 675– (2009) · Zbl 1272.74613 · doi:10.1016/j.advengsoft.2008.11.010
[4] Elman, Performance and analysis of saddle point preconditioners for the discrete steady-state Navier-Stokes equations, Numerische Mathematik 90 pp 665– (2002) · Zbl 1143.76531
[5] Bergamaschi, Novel preconditioners for the iterative solution to FE-discretized consolidation equations, Computer Methods in Applied Mechanics and Engineering 196 pp 2647– (2007) · Zbl 1173.76330
[6] Ferronato, Mixed constraint precoditioning in computational contact mechanics, Computer Methods in Applied Mechanics and Engineering 197 pp 3922– (2008) · Zbl 1194.74522
[7] Benzi, Preconditioning techniques for large linear systems: a survey, Journal of Computational Physics 182 pp 418– (2002) · Zbl 1015.65018
[8] Trottenberg, Multigrid (2001)
[9] Bank, Multilevel ILU decomposition, Numerische Mathematik 82 pp 543– (1999) · Zbl 0938.65063
[10] Saad, BILUM: block version of multi-elimination and multi-level ILU preconditioner for general sparse linear systems, SIAM Journal on Scientific Computing 20 pp 2103– (1999)
[11] Saad, BILUTM: A domain-based multilevel block ILUT preconditioner for general sparse matrices, SIAM Journal on Matrix Analysis and Applications 21 pp 279– (1999) · Zbl 0942.65045
[12] Meurant, A multilevel AINV preconditioner, Numerical Algorithms 29 pp 107– (2002)
[13] Saad, ARMS: an algebraic recursive multilevel solver for general sparse linear systems, Numerical Linear Algebra with Applications 9 pp 359– (2002) · Zbl 1071.65001
[14] Kraus, Algebraic multilevel preconditioning of finite element matrices using local Schur complements, Numerical Linear Algebra with Applications 13 pp 49– (2006) · Zbl 1174.65398
[15] Axelsson, Algebraic multilevel preconditioning methods, I, Numerische Mathematik 56 pp 157– (1989) · Zbl 0661.65110
[16] Axelsson, Algebraic multilevel preconditioning methods, II, SIAM Journal on Numerical Analysis 27 pp 1569– (1990) · Zbl 0715.65091
[17] Vassilevski, Multilevel Block Factorization Preconditioners: Matrix-based Analysis and Algorithms for Solving Finite Element Equations (2008) · Zbl 1170.65001
[18] Ferronato, Ill-conditioning of finite element poroelasicity equations, International Journal of Solids and Structures 38 pp 5995– (2001) · Zbl 1075.74643
[19] Pao, A fully coupled hydro-thermo-poro-mechanical model for black oil reservoir simulation, International Journal for Numerical and Analytical Methods in Geomechanics 25 pp 1229– (2001) · Zbl 1016.74023
[20] Schrefler, A fully coupled dynamic model for two phase flow in deformable porous media, Computer Methods in Applied Mechanics and Engineering 190 pp 3223– (2001) · Zbl 0977.74019
[21] Mazzia, On the reliability of numerical solutions for density dependent flow and transport in groundwater, Transport in Porous Media 43 pp 65– (2001)
[22] Zohdi, Modeling and simulation of a class of coupled thermo-chemo-mechanical processes in multiphase solids, Computer Methods in Applied Mechanics and Engineering 193 pp 679– (2004) · Zbl 1060.74528
[23] Saad, ILUT: a dual threshold incomplete LU factorization, Numerical Linear Algebra with Applications 1 pp 387– (1994) · Zbl 0838.65026
[24] Lin, Incomplete Cholesky factorizations with limited memory, SIAM Journal on Scientific Computing 21 pp 24– (1999) · Zbl 0941.65033
[25] Saad, ILUM: a multi-elimination ILU preconditioner for general sparse matrices, SIAM Journal on Scientific Computing 17 pp 830– (1996) · Zbl 0858.65029
[26] Axelsson, Diagonally compensated reduction and related preconditioning methods, Numerical Linear Algebra with Applications 1 pp 155– (1994) · Zbl 0837.65023
[27] Benzi, A robust incomplete factorization preconditioner for positive definite matrices, Numerical Linear Algebra with Applications 10 pp 385– (2003)
[28] Ajiz, A robust incomplete Cholesky-conjugate gradient algorithm, International Journal for Numerical Methods in Engineering 20 pp 949– (1984) · Zbl 0541.65019
[29] Saint-George, Problem-dependent preconditioners for iterative solvers in FE elastostatics, Computers and Structures 73 pp 33– (1999) · Zbl 1048.74591
[30] Tismenetsky, A new preconditioning technique for solving large sparse linear systems, Linear Algebra and its Applications 154-156 pp 331– (1991)
[31] Kaporin, High quality preconditioning of a general symmetric positive definite matrix based on its UTU+ UTR+ RTU-decomposition, Numerical Linear Algebra with Applications 5 pp 483– (1998)
[32] Biot, General theory of three-dimensional consolidation, Journal of Applied Physics 12 pp 155– (1941) · JFM 67.0837.01
[33] Geertsma, Land subsidence above compacting gas and oil reservoirs, Journal of Petroleum Technology 25 pp 734– (1973)
[34] Boade, Forecasting of Ekofisk reservoir compaction and subsidence by numerical simulation, Journal of Petroleum Technology 41 pp 723– (1989)
[35] Gambolati, Importance of poro-elastic coupling in dynamically active aquifers of the Po river basin, Italy, Water Resources Research 36 pp 2443– (2000)
[36] Ferronato, Numerical modelling of regional faults in land subsidence prediction above gas/oil reservoirs, International Journal for Numerical and Analytical Methods in Geomechanics 32 pp 633– (2008) · Zbl 1273.74532
[37] Baù, Basin-scale compressibility of the northern Adriatic by the radioactive marker technique, Gèotechnique 52 pp 605– (2002)
[38] Ferronato, Geomechanical issues of anthropogenic CO2 sequestration in exploited gas fields, Energy Conversion and Management (2009)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.