This paper is devoted to the numerical approximation of the solutions of a system of conservation laws arising in the modeling of two-phase flows in pipelines. The PDEs are closed by two highly nonlinear algebraic relations, namely, a pressure law and a hydrodynamic one. The severe nonlinearities encoded in these laws make the classical approximate Riemann solvers virtually intractable at a reasonable cost of evaluation. We propose a strategy for relaxing solely these two nonlinearities. The relaxation system we introduce is of course hyperbolic but all associated eigenfields are linearly degenerate. Such a property not only makes it trivial to solve the Riemann problem but also enables us to enforce some further stability requirements, in addition to those coming from a Chapman-Enskog analysis. The new method turns out to be fairly simple and robust while achieving desirable positivity properties on the density and the mass fractions. Extensive numerical evidences are provided.

This paper presents a finite volume discretization of two-phase Darcy flows in discrete fracture networks taking into account the mass exchange between the matrix and the fracture. We consider the asymptotic model for which the fractures are represented as interfaces of codimension one immersed in the matrix domain, leading to the so called hybrid dimensional Darcy flow model. The pressures at the interfaces between the matrix and the fracture network are continuous corresponding to a ratio between the normal permeability of the fracture and the width of the fracture assumed to be large compared with the ratio between the permeability of the matrix and the size of the domain. The discretization is an extension of the Vertex Approximate Gradient (VAG) scheme to the case of hybrid dimensional Darcy flow models. Compared with Control Volume Finite Element (CVFE) approaches, the VAG scheme has the advantage to avoid the mixing of the fracture and matrix rocktypes at the interfaces between the matrix and the fractures, while keeping the low cost of a nodal discretization on unstructured meshes. The convergence of the scheme is proved under the assumption that the relative permeabilities are bounded from below by a strictly positive constant. This assumption is needed in the convergence proof in order to take into account discontinuous capillary pressures in particular at the matrix fracture interfaces. The efficiency of our approach compared with CVFE discretizations is shown on two numerical examples of fracture networks in 2D and 3D.

Flow in fractured porous media occurs in the earth's subsurface, in biological tissues, and in man-made materials. Fractures have a dominating influence on flow processes, and the last decade has seen an extensive development of models and numerical methods that explicitly account for their presence. To support these developments, four benchmark cases for single-phase flow in three-dimensional fractured porous media are presented. The cases are specifically designed to test the methods' capabilities in handling various complexities common to the geometrical structures of fracture networks. Based on an open call for participation, results obtained with 17 numerical methods were collected. This paper presents the underlying mathematical model, an overview of the features of the participating numerical methods, and their performance in solving the benchmark cases.

We investigate the discretization of Darcy flow through fractured porous media on general meshes. We consider a hybrid dimensional model, invoking a complex network of planar fractures. The model accounts for matrix-fracture interactions and fractures acting either as drains or as barriers, i.e. we have to deal with pressure discontinuities at matrixfracture interfaces. The numerical analysis is performed in the general framework of gradient discretizations which is extended to the model under consideration. Two families of schemes namely the Vertex Approximate Gradient scheme (VAG) and the Hybrid Finite Volume scheme (HFV) are detailed and shown to satisfy the gradient scheme framework, which yields, in particular, convergence. Numerical tests confirm the theoretical results.

The gradient scheme family, which includes the conforming and mixed finite elements as well as the mimetic mixed hybrid family, is used for the approximation of Richards equation and the two-phase flow problem in heterogeneous porous media. We prove the convergence of the approximate saturation and of the approximate pressures and approximate pressure gradients thanks to monotony and compactness arguments under an assumption of non-degeneracy of the phase relative permeabilities. Strong convergence results stem from the convergence of the norms of the gradients of pressures, which demand handling the nonlinear time term. Numerical results show the efficiency on these problems of a particular gradient scheme, called the Vertex Approximate Gradient scheme.

No abstract

International audienceIn this work, we extend, to two-phase flow, the single-phase Darcy flow model proposed in [26], [12] in which the (d − 1)-dimensional flow in the fractures is coupled with the d-dimensional flow in the matrix. Three types of so called hybrid-dimensional two-phase Darcy flow models are proposed. They all account for fractures acting either as drains or as barriers, since they allow pressure jumps at the matrix-fracture interfaces. The models also permit to treat gravity dominated flow as well as discontinuous capillary pressure at the material interfaces. The three models differ by their transmission conditions at matrix fracture interfaces: while the first model accounts for the nonlinear two-phase Darcy flux conservations, the second and third ones are based on the linear single phase Darcy flux conservations combined with different approximations of the mobilities. We adapt the Vertex Approximate Gradient (VAG) scheme to this problem, in order to account for anisotropy and heterogeneity aspects as well as for applicability on general meshes. Several test cases are presented to compare our hybrid-dimensional models to the generic equi-dimensional model, in which fractures have the same dimension as the matrix, leading to deep insight about the quality of the proposed reduced models

International audienceThis paper concerns the discretization on general 3D meshes of multiphase compositional Darcy flows in heterogeneous anisotropic porous media. Extending Coats' formulation to an arbitrary number of phases, the model accounts for the coupling of the mass balance of each component with the pore volume conservation and the thermodynamical equilibrium, and dynamically manages phase appearance and disappearance. The spatial discretization of the multiphase compositional Darcy flows is based on a generalization of the Vertex Approximate Gradient scheme (VAG), already introduced for single phase diffusive problems. It leads to an unconditionally coercive scheme for arbitrary meshes and permeability tensors. The stencil of this vertex-centred scheme typically comprises 27 points on topologically Cartesian meshes, and the number of unknowns on tetrahedral meshes is considerably reduced, compared with usual cell-centred approaches. The efficiency of our approach is exhibited on the nearwell injection of miscible CO2 in a saline aquifer taking into account the vaporization of H2O in the gas phase as well as the precipitation of salt