Premium
Multi‐level incomplete factorizations for the iterative solution of non‐linear finite element problems
Author(s) -
Janna Carlo,
Ferronato Massimiliano,
Gambolati Giuseppe
Publication year - 2009
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.2664
Subject(s) - preconditioner , krylov subspace , cholesky decomposition , solver , incomplete cholesky factorization , conjugate gradient method , factorization , linear system , finite element method , computer science , positive definite matrix , mathematics , subspace topology , mathematical optimization , incomplete lu factorization , iterative method , algorithm , eigenvalues and eigenvectors , mathematical analysis , matrix decomposition , structural engineering , physics , quantum mechanics , artificial intelligence , engineering
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. Copyright © 2009 John Wiley & Sons, Ltd.