Asymptotic behavior of two-phase flows in heterogeneous porous media for capillarity depending only on space. I. Convergence to the optimal entropy solution

We consider an immiscible two-phase flow in a heterogeneous one-dimensional porous medium. We suppose particularly that the capillary pressure field is discontinuous with respect to the space variable. The dependence of the capillary pressure with respect to the oil saturation is supposed to be weak, at least for saturations which are not too close to 0 or 1. We study the asymptotic behavior when the capillary pressure tends to a function which does not depend on the saturation. In this paper, we show that if the capillary forces at the spacial discontinuities are oriented in the same direction that the gravity forces, or if the two phases move in the same direction, then the saturation profile with capillary diffusion converges toward the unique optimal entropy solution to the hyperbolic scalar conservation law with discontinuous flux functions

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 phase-by-phase upstream scheme that converges to the vanishing capillarity solution for countercurrent two-phase flow in two-rocks media

International audienceWe discuss the convergence of the upstream phase-by-phase scheme (or upstream mobility scheme) towards the vanishing capillarity solution for immiscible incompressible two-phase flows in porous media made of several rock types. Troubles in the convergence where recently pointed out in [S. Mishra & J. Jaffré, Comput. Geosci., 2010] and [S. Tveit & I. Aavatsmark, Comput. Geosci, 2012]. In this paper, we clarify the notion of vanishing capillarity solution, stressing the fact that the physically relevant notion of solution differs from the one inferred from the results of [E. F. Kaasschieter, Comput. Geosci., 1999]. In particular, we point out that the vanishing capillarity solution de- pends on the formally neglected capillary pressure curves, as it was recently proven in by the authors [B. Andreianov & C. Canc'es, Comput. Geosci., 2013]. Then, we propose a numerical procedure based on the hybridization of the interfaces that converges towards the vanishing capillarity solution. Numerical illustrations are provided

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