Premium
Direct, partitioned and projected solution to finite element consolidation models
Author(s) -
Gambolati Giuseppe,
Pini Giorgio,
Ferronato Massimiliano
Publication year - 2002
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.249
Subject(s) - conjugate gradient method , implementation , finite element method , consolidation (business) , robustness (evolution) , computer science , mathematics , software , mathematical optimization , projection method , algorithm , engineering , structural engineering , dykstra's projection algorithm , biochemistry , chemistry , accounting , business , gene , programming language
Direct, partitioned, and projected (conjugate gradient‐like) solution approaches are compared on unsymmetric indefinite systems arising from the finite element integration of coupled consolidation equations. The direct method is used in its most recent and computationally efficient implementations of the Harwell Software Library. The partitioned approach designed for coupled problems is especially attractive as it addresses two separate positive definite problems of a smaller size that can be solved by symmetric conjugate gradients. However, it may stagnate and when converging it does not prove competitive with a global projection method such as Bi‐CGSTAB, which may take full advantage of its flexibility in working on scaled and reordered equations, and thus may greatly improve its computational performance in terms of both robustness and convergence rate. The Bi‐CGSTAB superiority to the other approaches is discussed and demonstrated with a few representative examples in two‐dimensional (2‐D) and three‐dimensional (3‐D) coupled consolidation problems. Copyright © 2002 John Wiley & Sons, Ltd.