Published in

Elsevier, International Journal of Solids and Structures, 34-35(38), p. 5995-6014, 2001

DOI: 10.1016/s0020-7683(00)00352-8

Links

Tools

Export citation

Search in Google Scholar

Ill-conditioning of finite element poroelasicity equations

Journal article published in 2001 by Massimiliano Ferronato, Giuseppe Gambolati, Pietro Teatini
This paper is available in a repository.
This paper is available in a repository.

Full text: Download

Green circle
Preprint: archiving allowed
Orange circle
Postprint: archiving restricted
Red circle
Published version: archiving forbidden
Data provided by SHERPA/RoMEO

Abstract

The solution to Biot's coupled consolidation theory is usually addressed by the finite element (FE) method thus obtaining a system of first-order differential equations which is integrated by the use of an appropriate time marching scheme. For small values of the time step the resulting linear system may be severely ill-conditioned and hence the solution can prove quite difficult to achieve. Under such conditions efficient and robust projection solvers based on Krylov's subspaces which are usually recommended for non-symmetric large size problems can exhibit a very slow convergence rate or even fail. The present paper investigates the correlation between the ill-conditioning of FE poroelasticity equations and the time integration step Δt. An empirical relation is provided for a lower bound Δtcrit of Δt below which ill-conditioning may suddenly occur. The critical time step is larger for soft and low permeable porous media discretized on coarser grids. A limiting value for the rock stiffness is found such that for stiffer systems there is no ill-conditioning irrespective of Δt however small, as is also shown by several numerical examples. Finally, the definition of a different Δtcrit as suggested by other authors is reviewed and discussed.