In this paper, we study some discretization schemes for diffusive flows in heterogeneous anisotropic porous media. We first introduce the notion of gradient scheme, and show that several existing schemes fall into this framework. Then, we construct two new gradient schemes which have the advantage of a small stencil. Numerical results obtained for real reservoir meshes show the efficiency of the new schemes, compared to existing ones.
The topic of this work is to obtain discrete Sobolev inequalities for piecewise constant functions, and to deduce error estimates on the approximate solutions of convection diffusion equations by finite volume schemes.
Models of two phase flows in porous media, used in petroleum engineering, lead to a system of two coupled equations with elliptic and parabolic degenerate terms, and two unknowns, the saturation and the pressure. For the purpose of their approximation, a coupled scheme, consisting in a finite volume method together with a phase-by-phase upstream weighting scheme, is used in the industrial setting. This paper presents a mathematical analysis of this coupled scheme, first showing that it satisfies...
The topic of this work is to obtain discrete Sobolev inequalities for piecewise constant functions, and to deduce error estimates on the approximate solutions of convection diffusion equations by finite volume schemes.
Models of two phase flows in porous media, used in petroleum engineering, lead to a system of two coupled equations with elliptic and parabolic degenerate terms, and two unknowns, the saturation and the pressure. For the purpose of their approximation, a coupled scheme, consisting in a finite volume method together with a phase-by-phase upstream weighting scheme, is used in the industrial setting. This paper presents a mathematical analysis of this coupled scheme, first showing that it satisfies...
In this paper, we study some discretization schemes for diffusive flows in heterogeneous anisotropic porous media. We first introduce the notion of gradient scheme, and show that several existing schemes fall into this framework. Then, we construct two new gradient schemes which have the advantage of a small stencil. Numerical results obtained for real reservoir meshes show the efficiency of the new schemes, compared to existing ones.
We develop gradient schemes for the approximation of the Perona-Malik equations and nonlinear tensor-diffusion equations. We prove the convergence of these methods to the weak solutions of the corresponding nonlinear PDEs. A particular gradient scheme on rectangular meshes is then studied numerically with respect to experimental order of convergence which shows its second order accuracy. We present also numerical experiments related to image filtering by time-delayed Perona-Malik and tensor diffusion...
We design an abstract setting for the approximation in Banach spaces of operators acting in duality. A typical example are the gradient and divergence operators in Lebesgue-Sobolev spaces on a bounded domain. We apply this abstract setting to the numerical approximation of Leray-Lions type problems, which include in particular linear diffusion. The main interest of the abstract setting is to provide a unified convergence analysis that simultaneously covers (i) all usual boundary conditions, (ii)...
We present and analyse in this paper a novel colocated Finite Volume scheme for the solution of the Stokes problem. It has been developed following two main ideas. On one hand, the discretization of the pressure gradient term is built as the discrete transposed of the velocity divergence term, the latter being evaluated using a natural finite volume approximation; this leads to a non-standard interpolation formula for the expression of the pressure on the edges of the control volumes. On the other...
We present and analyse in this paper a novel cell-centered collocated finite volume scheme for incompressible flows. Its definition involves a partition of the set of control volumes; each element of this partition is called a cluster and consists in a few neighbouring control volumes. Under a simple geometrical assumption for the clusters, we obtain that the pair of discrete spaces associating the classical cell-centered approximation for the velocities and cluster-wide constant pressures is ...
We present in this paper a pressure correction scheme for the barotropic compressible Navier-Stokes equations, which enjoys an unconditional stability property, in the sense that the energy and maximum-principle-based estimates of the continuous problem also hold for the discrete solution. The stability proof is based on two independent results for general finite volume discretizations, both interesting for their own sake: the -stability of the discrete advection operator provided...
We present in this paper a pressure correction scheme for the drift-flux model combining finite element and finite volume discretizations, which is shown to enjoy essential stability features of the continuous problem: the scheme is conservative, the unknowns are kept within their physical bounds and, in the homogeneous case ( when the drift velocity vanishes), the discrete entropy of the system decreases; in addition, when using for the drift velocity a closure law which takes the form of a Darcy-like...
Modeling the kinetics of a precipitation dissolution reaction occurring in a porous medium where diffusion also takes place leads to a system of two parabolic equations and one ordinary differential equation coupled with a stiff reaction term. This system is discretized by a finite volume scheme which is suitable for the approximation of the discontinuous reaction term of unknown sign. Discrete solutions are shown to exist and converge towards a weak solution of the continuous problem. Uniqueness...
This paper concerns the discretization of multiphase Darcy flows, in the case of heterogeneous anisotropic porous media and general 3D meshes used in practice to represent reservoir and basin geometries. An unconditionally coercive and symmetric vertex centred approach is introduced in this paper. This scheme extends the Vertex Approximate Gradient scheme (VAG), already introduced for single phase diffusive problems in [9], to multiphase Darcy flows....
We review in this paper a class of schemes for the numerical simulation of compressible flows. In order to ensure the stability of the discretizations in a wide range of Mach numbers and introduce sufficient decoupling for the numerical resolution, we choose to implement and study pressure correction schemes on staggered meshes. The implicit version of the schemes is also considered for the theoretical study. We give both algorithms for the barotropic Navier-Stokes equations, for the full Navier-Stokes...
Page 1