1 Introduction
We use the framework of a contraction map to arrive at the convergence criterion for a staggered solution algorithm coupling small strain anisotropic poroelastoplasticity with single phase flow. As shown in Figure 1
, the flow subproblem is solved with stress tensor fixed followed by the poromechanics subproblem in every coupling iteration at each time step. The coupling iterations are repeated until convergence and Backward Euler is employed for time marching. The analysis is motivated by the results in our previous works as follows

In [1], the contraction map for twogrid staggered algorithm lent us closed form expressions for coarse scale moduli in terms of fine scale data. The flow equations were solved on a fine grid and the isotropic poroelasticity equations were solved on a coarse grid.

In [2], we used contraction map to demonstrate convergence of staggered solution algorithm for anisotropic poroelasticity coupled with single phase flow. The speciality of this algorithm was that the stress tensor was fixed during the flow solve as an extension to the case with isotropic poroelasticity in which the mean stress was fixed during the flow solve.
2 Model equations
2.1 Flow model
Let the boundary where is the Dirichlet boundary and is the Neumann boundary. The equations are
(2.1) 
where is the fluid pressure, is the fluid flux, is the increment in fluid content^{1}^{1}1[3] defines the increment in fluid content as the measure of the amount of fluid which has flowed in and out of a given element attached to the solid frame, is the source or sink term, is the uniformly symmetric positive definite absolute permeability tensor, is the fluid viscosity, is a reference density, is the porosity, is a measure of the hydraulic conductivity of the pore fluid, is the fluid compressibility and is the time interval.
2.2 Poromechanics model
The important phenomenological aspects of small strain elastoplasticity are

The existence of an elastic domain, i.e. a range of stresses within which the behaviour of the material can be considered as purely elastic, without evolution of permanent (plastic) strains. The elastic domain is delimited by the socalled yield stress. A scalar yield function is introduced. The yield locus is the boundary of the elastic domain where and the corresponding yield surface is defined as .

If the material is further loaded at the yield stress, then plastic yielding (or plastic flow), i.e. evolution of plastic strains, takes place.
Let the boundary where is the Dirichlet boundary and is the Neumann boundary. The equations are
(2.2) 
where is the solid displacement, is the rock density, is the body force per unit volume, is the traction specified on , is the strain tensor, and are the elastic and plastic parts of strain tensor respctively, is the Cauchy stress tensor, is the fourth order symmetric positive definite anisotropic elasticity tensor, is the Biot tensor and is the plastic multiplier satisfying the complementarity condition
(2.3) 
The inverse of the constitutive law is
(2.4) 
where is a generalized Hooke’s law constant (see [18]) and is the Skempton pore pressure coefficient (see [5]).
2.3 Increment in fluid content
3 Statement of convergence
As elucidated in Figure 2, we use the notations and for the change in the quantity during the flow and poromechanics solves respectively over the coupling iteration and for the change in the quantity over the coupling iteration at any time level such that
The problem statement is: find , and such that
(3.7)  
(3.8)  
(3.9) 
where the finite dimensional spaces , and are
where represents the space of constants, represents the space of trilinears and the details of are given in [2]. The details of (3.7), (3.8) and (3.9) are given in Appendices A, B and C respectively. The fixed stress split iterative scheme is a contraction map given by
where the term is driven to a small value by the convergence criterion.
Step 1: Flow equations
Testing (3.7) with and (3.8) with , we get
(3.10)  
(3.11) 
From (3.10) and (3.11), we get
(3.12) 
Step 2: Poromechanics equations
Testing (3.9) with , we get
(3.13) 
We now invoke (2.4) to arrive at . Substituting in (3.13), we get
(3.14) 
Step 3: Combining flow and poromechanics equations
Adding (3.12) and (3.14), we get
(3.15) 
Step 4: Variation in fluid content
In lieu of (2.5), we write
(3.16) 
which can be written as
Dividing throughout by , we get
(3.17) 
From (3.15) and (3.17), we get
(3.18) 
Adding and subtracting to the LHS of (3.18) results in
Multiplying throughout by results in
(3.19) 
Step 5: Invoking the fixed stress constraint
In lieu of (2.5) and the fixed stress constraint during the flow solve, we get
Further, since the pore pressure is frozen during the poromechanical solve, we have . As a result, we can write
(3.20) 
Subtracting (3.20) from (3.16), we can write
(3.21) 
which implies that
(3.22) 
In lieu of (3.22), we can write (3.19) as
(3.23) 
Step 6: Invoking the Young’s inequality
We invoke the Young’s inequality (see [8])
for the RHS of (3.23) as follows
(3.24) 
In lieu of (3.24), we write (3.23) as
4 Convergence criterion
We desire to drive the following quantity to zero
(4.25) 
In lieu of (3.21), we can write
(4.26) 
In lieu of (4.26), we can write (4.25) as
which can also be written as
As a result, we pose the convergence criterion as
(4.27) 
where is a prespecified tolerance and represents a small value.
4.1 Computation of quantities of interest
In lieu of (2.5) and (2.6), we can write
(4.28)  
(4.29)  
(4.30) 
where we keep in mind that the pore pressure is frozen during the poromechanics solve and stress tensor is fixed during the flow solve. The quantity is obtained after the flow solve while the quantities and are obtained after the poromechanical solve.
To understand why , we present the basic algorithmic framework for the solution of elastoplastic equations: The system of equations (2.2) is first solved with for a trial stress state .

If , then we proceed with .
In summary, the solution is such that . During the subsequent flow solve, since the stress tensor is fixed, the value of the yield function does not change i.e. still. This implies that during the flow solve and the porous solid does not accumulate any plastic strain during the flow solve.
Appendix A Discrete variational statement of mass conservation
In lieu of (2.5), we write mass conservation equation as
(1.32) 
The discrete in time form of (1.32) in the time step is written as
where is the time step. The fixed stress split constraint implies that gets replaced by as is fixed during the flow solve. The modified equation is written as
As a result, the discrete weak form of mass conservation is given by
Replacing by and subtracting the two equations, we get
Comments
There are no comments yet.