On stability and convergence of finite element approximations of biot's consolidation problem

Stability and convergence analysis of finite element approximations of Biot's equations governing quasistatic consolidation of saturated porous media are, discussed. A family of decay functions, parametrized by the number of time steps, is derived for the fully discrete backward Euler–Galerkin formulation, showing that the pore-pressure oscillations, arising from an unstable approximation of the incompressibility constraint on the initial condition, decay in time. Error estimates holding over the unbounded time domain for both semidiscrete and fully discrete formulations are presented, and a post-processing technique is employed to improve the pore-pressure accuracy.