A coupling of mixed and continuous Galerkin finite element methods for poroelasticity I: the continuous in time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot's consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a continuous in-time setting for a strictly positive constrained specific storage coefficient. Of particular interest is the case when the lowest-order Raviart-Thomas approximating space or cell-centered finite differences are used in the mixed formulation, and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.