Premium
Scaling improves stability of preconditioned CG‐like solvers for FE consolidation equations
Author(s) -
Gambolati Giuseppe,
Pini Giorgio,
Ferronato Massimiliano
Publication year - 2003
Publication title -
international journal for numerical and analytical methods in geomechanics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.419
H-Index - 91
eISSN - 1096-9853
pISSN - 0363-9061
DOI - 10.1002/nag.311
Subject(s) - preconditioner , solver , scaling , consolidation (business) , mathematics , logarithm , rate of convergence , conjugate gradient method , factorization , matrix (chemical analysis) , finite element method , incomplete lu factorization , mathematical optimization , computer science , mathematical analysis , algorithm , iterative method , matrix decomposition , geometry , physics , eigenvalues and eigenvectors , materials science , channel (broadcasting) , computer network , accounting , business , composite material , thermodynamics , quantum mechanics
Preconditioned projection (or conjugate gradient like) methods are increasingly used for the accurate and efficient solution to finite element (FE) coupled consolidation equations. Theory indicates that preliminary row/column scaling does not affect the eigenspectrum of the iteration matrix controlling convergence as long as the preconditioner relies on the incomplete factorization of the FE coefficient matrix. However, computational experience with mid‐large size problems shows that the above inexpensive operation can significantly accelerate the solver convergence, and to a minor extent also improve the final accuracy, as a result of a better solver stability to the accumulation and propagation of floating point round‐off errors. This is demonstrated with the aid of the least square logarithm (LSL) scaling algorithm on FE consolidation problems of increasing size up to more than 100 000. It is shown that a major source of numerical instability rests with the sub‐matrix which couples the structural to the fluid part of the underlying mathematical model. It is concluded that for mid‐large size, possibly difficult, FE consolidation problems left/right LSL scaling is to be always recommended when the incomplete factorization is used as a preconditioning technique. Copyright © 2003 John Wiley & Sons, Ltd.