This paper is concerned with numerical algorithms for Biot model. By introducing an intermediate variable, the classical 2-field Biot model is written into a 3-field formulation. Based on such a 3-field formulation, we propose a coupled algorithm, some time-extrapolation based decoupled algorithms, and an iterative decoupled algorithm. Our focus is the analysis of the iterative decoupled algorithm. It is shown that the convergence of the iterative decoupled algorithm requires no extra assumptions on physical parameters or stabilization parameters. Numerical experiments are provided to demonstrate the accuracy and efficiency of the proposed method.

This paper concentrates on a priori error estimates of two monolithic schemes for Biot's consolidation model based on the three‐field formulation introduced by Oyarzúa et al. (SIAM J Numer Anal, 2016). The spatial discretizations are based on the Taylor–Hood finite elements combined with Lagrange elements for the three primary variables. We employ two different schemes to discretize the time domain. One uses the backward Euler method, and the other applies the combination of the backward Euler and Crank‐Nicolson methods. A priori error estimates show that both schemes are unconditionally convergent with optimal error orders. Detailed numerical experiments are presented to validate the theoretical analysis.

No abstract