Parallel numerical modeling of hybrid-dimensional compositional non-isothermal Darcy flows in fractured porous media
For researchers in subsurface flow simulation, this work provides a rigorous derivation and efficient parallel implementation of a hybrid-dimensional model for complex fracture networks, though it is an incremental extension of existing methods.
This paper introduces a new discrete fracture model for non-isothermal compositional multiphase Darcy flows in fractured porous media, using a hybrid-dimensional approach (2D in fractures, 3D in matrix) discretized with a fully implicit VAG finite volume scheme and parallel SPMD implementation. Numerical tests demonstrate parallel scalability, nonlinear and linear convergence on various meshes and fracture networks.
This paper introduces a new discrete fracture model accounting for non-isothermal compositional multiphase Darcy flows and complex networks of fractures with intersecting, immersed and non immersed fractures. The so called hybrid-dimensional model using a 2D model in the fractures coupled with a 3D model in the matrix is first derived rigorously starting from the equi-dimensional matrix fracture model. Then, it is dis-cretized using a fully implicit time integration combined with the Vertex Approximate Gradient (VAG) finite volume scheme which is adapted to polyhedral meshes and anisotropic heterogeneous media. The fully coupled systems are assembled and solved in parallel using the Single Program Multiple Data (SPMD) paradigm with one layer of ghost cells. This strategy allows for a local assembly of the discrete systems. An efficient preconditioner is implemented to solve the linear systems at each time step and each Newton type iteration of the simulation. The numerical efficiency of our approach is assessed on different meshes, fracture networks, and physical settings in terms of parallel scalability, nonlinear convergence and linear convergence.