Abstract:We consider a poromechanics model including frictionless contact mechanics. The resulting model consists of the Biot equations with contact boundary conditions leading to a variational inequality modelling mechanical deformations coupled to a linear parabolic flow equation. We propose a fully discrete iterative scheme for solving this model. This scheme decoupled the flow and mechanics equations and extends the fixed-stress splitting scheme for the Biot equations. We use finite elements in space and a backward Euler discretization in time. We show that the fixed stress split scheme is a contraction.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the modeling and numerical solution of contact problems in porous media, especially when fluid action is involved. Specifically, the author focuses on the frictionless contact mechanics problem in porous media and combines it with the Biot equation to form a variational inequality model with contact boundary - value conditions.
### Problem Background
1. **Contact Mechanics Problems in Porous Media**: In porous media, especially in the presence of fractures or faults, contact mechanics problems become particularly important. These areas are usually the places where mechanical failure is most likely to occur and also have a significant impact on fluid flow.
2. **Coupling of Fluid Flow and Mechanical Deformation**: Underground engineering processes (such as oil and gas extraction, geothermal energy production, geological hazard assessment) involve the injection or extraction of fluids, which will lead to pressure changes and stress imbalances, and may further trigger seismic events. Therefore, understanding the mechanical behavior of fractures and faults and their interaction with fluids is an important topic in earth science.
### Core Problems of the Paper
The paper aims to propose a fixed - stress splitting method for solving contact problems in porous media and prove its convergence. Specifically:
- **Model Equations**: The author considers the quasi - static Biot model and imposes contact conditions on part of its boundaries, including the Signorini condition in the case of frictionless contact.
- **Numerical Method**: A fully - discrete iterative scheme is proposed to solve this model, which extends the existing fixed - stress splitting scheme. The finite element method is used for spatial discretization, and the backward Euler method is used for time discretization.
- **Convergence Analysis**: It is proved that the proposed fixed - stress splitting scheme is a contraction mapping, thus ensuring the convergence of the algorithm.
### Mathematical Model
The contact problem in porous media can be expressed in the following form:
\[
\text{Problem}: \text{(Signorini condition for frictionless contact with gap function in porous media)}
\]
Seek the displacement field \( \mathbf{u}: \Omega \to \mathbb{R}^d \) and the stress field \( \sigma: \Omega \to S^d \), satisfying:
\[
\sigma_{\text{por}}(\mathbf{u}, p)=\sigma(\mathbf{u})-\alpha pI, \quad \sigma(\mathbf{u})=\lambda (\nabla \cdot \mathbf{u})I + 2\mu \varepsilon(\mathbf{u}), \quad \varepsilon(\mathbf{u})=\frac{1}{2}(\nabla \mathbf{u}+\nabla \mathbf{u}^T), \quad \nabla \cdot \sigma_{\text{por}}+\mathbf{f} = 0.
\]
### Boundary Conditions
\[
\mathbf{u} = 0 \text{ on } \Gamma_1, \quad \sigma \mathbf{n}=\mathbf{t}_2 \text{ on } \Gamma_2,
\]
\[
u_n \leq g, \quad t_n \leq 0, \quad t_n(u_n - g)=0 \text{ on } \Gamma_3, \quad \sigma \mathbf{t} = 0 \text{ on } \Gamma_3.
\]
### Flow Equations
The flow equations include the mass conservation equation:
\[
\frac{\partial}{\partial t}\left(\frac{1}{M}p+\alpha \nabla \cdot \mathbf{u}\right)+\nabla \cdot \mathbf{q}=f,
\]
where \( M \) is the total compressibility coefficient, and \( \mathbf{q} \) is