A Discontinuous Control Volume Finite Element Method for Multi-Phase Flow in Heterogeneous Porous Media

We present a new, high-order, control-volume-finite-element (CVFE) method for multiphase porous media flow with discontinuous 1st-order representation for pressure and discontinuous 2nd-order representation for velocity. The method has been implemented using unstructured tetrahedral meshes to discretize space. The method locally and globally conserves mass. However, unlike conventional CVFE formulations, the method presented here does not require the use of control volumes (CVs) that span the boundaries between domains with differing material properties. We demonstrate that the approach accurately preserves discontinuous saturation changes caused by permeability variations across such boundaries, allowing efficient simulation of flow in highly heterogeneous models. Moreover, accurate solutions are obtained at significantly lower computational cost than using conventional CVFE methods. We resolve a long-standing problem associated with the use of classical CVFE methods to model flow in highly heterogeneous porous media

A Force-Balanced Control Volume Finite Element Method for Multi-Phase Porous Media Flow Modelling

Dr D. Pavlidis would like to acknowledge the support from the following research grants: Innovate UK ‘Octopus’, EPSRC ‘Reactor Core-Structure Re-location Modelling for Severe Nuclear Accidents’) and Horizon 2020 ‘In-Vessel Melt Retention’. Funding for Dr P. Salinas from ExxonMobil is gratefully acknowledged. Dr Z. Xie is supported by EPSRC ‘Multi-Scale Exploration of Multi-phase Physics in Flows’. Part funding for Prof Jackson under the TOTAL Chairs programme at Imperial College is also acknowledged. The authors would also like to acknowledge Mr Y. Debbabi for supplying analytic solutions.Peer reviewedPublisher PD

A discontinuous control volume finite element method for multi-phase flow in heterogeneous porous media

We present a new, high-order, control-volume-finite-element (CVFE) method for multiphase porous media flow with discontinuous 1st-order representation for pressure and discontinuous 2nd-order representation for velocity. The method has been implemented using unstructured tetrahedral meshes to discretize space. The method locally and globally conserves mass. However, unlike conventional CVFE formulations, the method presented here does not require the use of control volumes (CVs) that span the boundaries between domains with differing material properties. We demonstrate that the approach accurately preserves discontinuous saturation changes caused by permeability variations across such boundaries, allowing efficient simulation of flow in highly heterogeneous models. Moreover, accurate solutions are obtained at significantly lower computational cost than using conventional CVFE methods. We resolve a long-standing problem associated with the use of classical CVFE methods to model flow in highly heterogeneous porous media

A fully-coupled discontinuous Galerkin method for two-phase flow in porous media with discontinuous capillary pressure

In this paper we formulate and test numerically a fully-coupled discontinuous Galerkin (DG) method for incompressible two-phase flow with discontinuous capillary pressure. The spatial discretization uses the symmetric interior penalty DG formulation with weighted averages and is based on a wetting-phase potential / capillary potential formulation of the two-phase flow system. After discretizing in time with diagonally implicit Runge-Kutta schemes the resulting systems of nonlinear algebraic equations are solved with Newton's method and the arising systems of linear equations are solved efficiently and in parallel with an algebraic multigrid method. The new scheme is investigated for various test problems from the literature and is also compared to a cell-centered finite volume scheme in terms of accuracy and time to solution. We find that the method is accurate, robust and efficient. In particular no post-processing of the DG velocity field is necessary in contrast to results reported by several authors for decoupled schemes. Moreover, the solver scales well in parallel and three-dimensional problems with up to nearly 100 million degrees of freedom per time step have been computed on 1000 processors

A force-balanced control volume finite element method for multi-phase porous media flow modelling

A novel method for simulating multi-phase flow in porous media is presented. The approach is based on a control volume finite element mixed formulation and new force-balanced finite element pairs. The novelty of the method lies in: (a) permitting both continuous and discontinuous description of pressure and saturation between elements; (b) the use of arbitrarily high-order polynomial representation for pressure and velocity and (c) the use of high-order flux-limited methods in space and to time avoid introducing non-physical oscillations while achieving high-order accuracy where and when possible. The model is initially validated for two-phase flow. Results are in good agreement with analytically obtained solutions and experimental results. The potential of this method is demonstrated by simulating flow in a realistic geometry composed of highly permeable meandering channels

Parallel numerical modeling of hybrid-dimensional compositional non-isothermal Darcy flows in fractured porous media

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

A Cell-Centred CVD-MPFA Finite Volume Method for Two-Phase Fluid Flow Problems with Capillary Heterogeneity and Discontinuity

A novel finite-volume method is presented for porous media flow simulation that is applicable to discontinuous capillary pressure fields. The method crucially retains the optimal single of freedom per control-volume being developed within the flux-continuous control-volume distributed multi-point flux approximation (CVD-MPFA) framework (Edwards and Rogers in Comput Geosci 02(04):259–290, 1998; Friis et al. in SIAM J Sci Comput 31(02):1192–1220, 2008) . The new methods enable critical subsurface flow processes involving oil and gas trapping to be correctly resolved on structured and unstructured grids. The results demonstrate the ability of the method to resolve flow with oil/gas trapping in the presence of a discontinuous capillary pressure field for diagonal and full-tensor permeability fields. In addition to an upwind approximation for the saturation equation flux, the importance of upwinding on capillary pressure flux via a novel hybrid formulation is demonstrated