This paper presents a review of the available mathematical models and corresponding non-conforming numerical approximations which describe single-phase fluid flow in a fractured porous medium. One focus is on the geometrical difficulties that may arise in realistic simulations such as intersecting and immersed fractures. Another important aspect is the choice of the approximation spaces for the discrete problem: in mixed formulations, both the Darcy velocity and the pressure are considered as unknowns, while in classical primal formulations, a richer space for the pressure is considered and the Darcy velocity is computed a posteriori. In both cases, the extended finite element method is used, which allows for a complete geometrical decoupling among the fractures and rock matrix grids. The fracture geometries can thus be independent of the underlying grid thanks to suitable enrichments of the spaces that are able to represent possible jumps of the solution across the fractures. Finally, due to the dimensional reduction, a better approximation of the resulting boundary conditions for the fractures is addressed.
A Review of the XFEM-Based Approximation of Flow in Fractured Porous Media
FUMAGALLI, ALESSIO;SCOTTI, ANNA
2016-01-01
Abstract
This paper presents a review of the available mathematical models and corresponding non-conforming numerical approximations which describe single-phase fluid flow in a fractured porous medium. One focus is on the geometrical difficulties that may arise in realistic simulations such as intersecting and immersed fractures. Another important aspect is the choice of the approximation spaces for the discrete problem: in mixed formulations, both the Darcy velocity and the pressure are considered as unknowns, while in classical primal formulations, a richer space for the pressure is considered and the Darcy velocity is computed a posteriori. In both cases, the extended finite element method is used, which allows for a complete geometrical decoupling among the fractures and rock matrix grids. The fracture geometries can thus be independent of the underlying grid thanks to suitable enrichments of the spaces that are able to represent possible jumps of the solution across the fractures. Finally, due to the dimensional reduction, a better approximation of the resulting boundary conditions for the fractures is addressed.I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.