z-logo
Premium
Finite volume discretization for poroelastic media with fractures modeled by contact mechanics
Author(s) -
Berge Runar L.,
Berre Inga,
Keilegavlen Eirik,
Nordbotten Jan M.,
Wohlmuth Barbara
Publication year - 2019
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.6238
Subject(s) - poromechanics , biot number , discretization , lagrange multiplier , mathematics , finite element method , unilateral contact , mathematical analysis , variational inequality , finite volume method , piecewise , mathematical optimization , mechanics , physics , porous medium , materials science , porosity , composite material , thermodynamics
Summary A fractured poroelastic body is considered where the opening of the fractures is governed by a nonpenetration law, whereas slip is described by a Coulomb‐type friction law. This physical model results in a nonlinear variational inequality problem. The variational inequality is rewritten as a complementary function, and a semismooth Newton method is used to solve the system of equations. For the discretization, we use a hybrid scheme where the displacements are given in terms of degrees of freedom per element, and an additional Lagrange multiplier representing the traction is added on the fracture faces. The novelty of our method comes from combining the Lagrange multiplier from the hybrid scheme with a finite volume discretization of the poroelastic Biot equation, which allows us to directly impose the inequality constraints on each subface. The convergence of the method is studied for several challenging geometries in 2D and 3D, showing that the convergence rates of the finite volume scheme do not deteriorate when it is coupled to the Lagrange multipliers. Our method is especially attractive for the poroelastic problem because it allows for a straightforward coupling between the matrix deformation, contact conditions, and fluid pressure.

This content is not available in your region!

Continue researching here.

Having issues? You can contact us here