Efficient Solvers for a Stabilized Three-Field Mixed Formulation of Poroelasticity

We focus on a three-field (displacement-velocity-pressure) stabilized mixed method for poroelasticity based on piecewise trilinear (Q1), lowest order Raviart-Thomas (RT0), and piecewise constant (P0) approximations for displacement, Darcy’s velocity and fluid pore pressure, respectively. Since the selected discrete spaces do not intrinsically satisfy the inf-sup condition in the undrained/incompressible limit, we propose a stabilization strategy based on local pressure jumps. Then, we focus on the efficient solution of the stabilized formulation by a block preconditioned Krylov method. Robustness and efficiency of the proposed approach are demonstrated in two sets of numerical experiments.