ON STABILITY AND CONVERGENCE OF FINITE-ELEMENT APPROXIMATIONS OF BIOTS CONSOLIDATION PROBLEM

Stability and convergence analysis of finite element approximations of Blot'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.