Comparison between iterative solution of symmetric and non‐symmetric forms of Biot's FEM equations using the generalized Jacobi preconditioner

Finite element discretization of Biot's consolidation equations can produce a symmetric indefinite system (commonly used in geomechanics) or a non‐symmetric system. While this difference appears to be minor, however, it will require the selection of entirely different Krylov subspace solvers with potentially significant impact on solution efficiency. The former is solved using the symmetric quasi‐minimal residual whereas the latter is solved using the popular bi‐conjugate gradient stabilized. This paper presents an extensive comparison of the symmetric and non‐symmetric forms by varying the time step, size of the spatial domain, choice of physical units, and left versus left–right preconditioning. The generalized Jacobi (GJ) preconditioner is able to handle the non‐symmetric version of Biot's finite element method equation, although there are no practical incentives to do so. The convergence behaviour of GJ‐preconditioned systems and its relation to the spectral condition number or the complete spectrum are studied to clarify the concept of ill‐conditioning within the context of iteration solvers. Copyright © 2007 John Wiley & Sons, Ltd.