The paper deals with the computational framework for the numerical simulation of the three dimensional fluid-filled fracture evolution in a poroelastic medium. The model consists of several groups of equations including the Biot poroelastic model to describe a bulk medium behavior, Reynold’s lubrication equations to describe a flow inside fracture and corresponding bulk/fracture interface conditions. The geometric model of the fracture assumes that it is described as an arbitrary sufficiently smooth surface with a boundary. Main attention is paid to describing numerical algorithms for particular problems (poroelasticity, fracture fluid flow, fracture evolution) as well as an algorithm for the coupled problem solution. An implicit fracture mid-surface representation approach based on the closest point projection operator is a particular feature of the proposed algorithms. Such a representation is used to describe the fracture mid-surface in the poroelastic solver, Reynold’s lubrication equation solver and for simulation of fracture evolutions. The poroelastic solver is based on a special variant of X-FEM algorithms, which uses the closest point representation of the fracture. To solve Reynold’s lubrication equations, which model the fluid flow in fracture, a finite element version of the closet point projection method for PDEs surface is used. As a result, the algorithm for the coupled problem is purely Eulerian and uses the same finite element mesh to solve equations defined in the bulk and on the fracture mid-surface. Finally, we present results of the numerical simulations which demonstrate possibilities of the proposed numerical techniques, in particular, a problem in a media with a heterogeneous distribution of transport, elastic and toughness properties.
