Iteratively Coupled Flow and Geomechanics in Fractured Poroelastic Reservoirs: A Phase Field Fracture Model

Algorithm 1

Extended fixed stress iterative coupling method.

1 Set phase field Dirichlet boundary conditions for pre-existing fracture networks;

2 Set initial conditions;

3 Solve the phase field formulation Equation (33) for ;

4 At the time ;

5 repeat

6  Update the fracture outward normal using Equation (5);

7  Update the fracture width using Equation (13);

8  Solve the pressure formulation Equation (32) for using , , , , and ;

9  Solve the displacement formulation Equation (31) for using and ;

10 until;

11 increment the time .