Abstract:The present work deals with the numerical resolution of coupled 3D-2D problems arising from the simulation of fluid flow in fractured porous media modeled via the Discrete Fracture and Matrix (DFM) model. According to the DFM model, fractures are represented as planar interfaces immersed in a 3D porous matrix and can behave as preferential flow paths, in the case of conductive fractures, or can actually be a barrier for the flow, when, instead, the permeability in the normal-to-fracture direction is small compared to the permeability of the matrix. Consequently, the pressure solution in a DFM can be discontinuous across a barrier, as a result of the geometrical dimensional reduction operated on the fracture.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop a numerical scheme suitable for simulating fluid flow in porous media with fractures and barriers. Specifically, it aims to address the following challenges:
1. **Influence of fractures and barriers**: In the discrete fracture and matrix (DFM) model, fractures can serve as preferential flow paths (when they have higher permeability) or as barriers to flow (when their permeability in the direction perpendicular to the fracture is lower). Therefore, the pressure solution may be discontinuous at the barriers in DFM.
2. **Use of non - conforming grids**: To simplify calculations and improve efficiency, the paper proposes a method that can perform simulations when the grid of the 3D matrix does not conform to the fractures. This helps avoid the difficulties associated with generating complex grids that conform to the fractures.
3. **Domain decomposition strategy**: By introducing an optimization - based domain decomposition method, this method can effectively handle large - scale problems and can be solved on non - conforming grids.
4. **Application of the extended finite element method (XFEM)**: To reproduce the discontinuous pressure solution around the barriers, the paper adopts the extended finite element method, which is an effective method for describing irregular solutions on non - conforming grids.
### Specific objectives
- **Develop a five - field formulation**: By introducing three interface unknowns (a total of five independent fields), decouple the 3D problem from the 2D fracture problem, so that the global solution can be regarded as a problem of minimizing the matching error on the interface.
- **Prove the well - posedness of the discrete problem**: Ensure that the proposed discrete problem is well - posed mathematically, that is, there is a unique solution and the solution is stable with respect to changes in the data.
- **Verify the effectiveness of the method**: Demonstrate the feasibility and effectiveness of this method through four verification examples.
### Mathematical formulas
The key formulas involved in the paper include, but are not limited to:
- **Governing equations**:
\[
-\nabla\cdot\left(K_D\nabla H_D\right)=g\quad\text{in}\;\mathring{D}
\]
\[
-\nabla_F\cdot\left(K_F\nabla_F H_F\right)= - Q^+ + Q^-\quad\text{in}\;F
\]
\[
Q^\pm = -\eta(H_D^\pm - H_F)\quad\text{on}\;\omega^\pm
\]
- **Weak form**:
\[
a_D(H_D, v)+b(H_D^+, v^+)+b(H_D^-, v^-)= (g, v)_{\mathring{D}}+b(H_F, v^+)+b(H_F, v^-),\quad\forall v\in H_0^1(\mathring{D})
\]
\[
a_F(H_F, w)+2b(H_F, w)=b(H_D^+, w)+b(H_D^-, w),\quad\forall w\in V
\]
- **Optimization objective function**:
\[
J=\|\Psi^+ - H_D^+\|_Q^2+\|\Psi^- - H_D^-\|_Q^2+\|\Psi_F - H_F\|_V^2
\]
Through these formulas, the paper shows how to effectively simulate fluid flow in porous media under complex geometric structures and physical conditions.