Mixed and Nitsche’s Discretizations of Frictional Contact-Mechanics in Fractured Porous Media

Abstract

This work deals with the discretization of single-phase Darcy flows in fractured and deformable porous media, including frictional contact at the matrix-fracture interfaces. Fractures are described as a network of planar surfaces leading to so-called mixed-dimensional models. Small displacements and a linear poro-elastic behavior are considered in the matrix. One key difficulty to simulate such coupled poro-mechanical models is related to the formulation and discretization of the contact mechanical sub-problem. Our starting point is based on the mixed formulation using facewise constant Lagrange multipliers along the fractures representing normal and tangential stresses. This is a natural choice for the discretization of the contact dual cone in order to account for complex fracture networks with corners and intersections. It leads to local expressions of the contact conditions and to efficient semi-smooth nonlinear solvers. On the other hand, such a mixed formulation requires to satisfy a compatibility condition between the discrete spaces restricting the choice of the displacement space and potentially leading to sub-optimal accuracy. This motivates the investigation of two alternative formulations based either on a stabilized mixed formulation or on the Nitsche’s method. These three types of formulations are first investigated theoretically in order to reveal the connections between them. Then, they are compared numerically in terms of accuracy and nonlinear convergence on a coupled poromechanical 2D model.