Finite-difference analysis of fully dynamic problems for saturated porous media

Finite-difference methods, using staggered grids in space, are considered for the numerical approximation of fully dynamic poroelasticity problems. First, a family of second-order schemes in time is analyzed. A priori estimates for displacements in discrete energy norms are obtained and the corresponding convergence results are proved. Numerical examples are given to illustrate the convergence properties of these methods. As in the case of an incompressible fluid and small permeability, these schemes suffer from spurious oscillations in time, a first order scheme is proposed and analyzed. For this new scheme a priori estimates and convergence results are also given. Finally, numerical examples in one and two dimensions are presented to show the good monotonicity properties of this method.