Unified convergence analysis of numerical schemes for a miscible displacement problem

This article performs a unified convergence analysis of a variety of numerical methods for a model of the miscible displacement of one incompressible fluid by another through a porous medium. The unified analysis is enabled through the framework of the gradient discretisation method for diffusion operators on generic grids. We use it to establish a novel convergence result in $L^\infty(0,T; L^2(\Omega))$ of the approximate concentration using minimal regularity assumptions on the solution to the continuous problem. The convection term in the concentration equation is discretised using a centred scheme. We present a variety of numerical tests from the literature, as well as a novel analytical test case. The performance of two schemes are compared on these tests; both are poor in the case of variable viscosity, small diffusion and medium to small time steps. We show that upstreaming is not a good option to recover stable and accurate solutions, and we propose a correction to recover stable and accurate schemes for all time steps and all ranges of diffusion

An HMM--ELLAM scheme on generic polygonal meshes for miscible incompressible flows in porous media

We design a numerical approximation of a system of partial differential equations modelling the miscible displacement of a fluid by another in a porous medium. The advective part of the system is discretised using a characteristic method, and the diffusive parts by a finite volume method. The scheme is applicable on generic (possibly non-conforming) meshes as encountered in applications. The main features of our work are the reconstruction of a Darcy velocity, from the discrete pressure fluxes, that enjoys a local consistency property, an analysis of implementation issues faced when tracking, via the characteristic method, distorted cells, and a new treatment of cells near the injection well that accounts better for the conservativity of the injected fluid

Discontinuous Galerkin Finite Element Convergence for Incompressible Miscible Displacement Problems of Low Regularity

In this article we analyse the numerical approximation of incompressible miscible displacement problems with a combined mixed finite element and discontinuous Galerkin method under minimal regularity assumptions. The main result is that sequences of discrete solutions weakly accumulate at weak solutions of the continuous problem. In order to deal with the non-conformity of the method and to avoid overpenalisation of jumps across interelement boundaries, the careful construction of a reflexive subspace of the space of bounded variation, which compactly embeds into $L^2(\Omega)$, and of a lifting operator, which is compatible with the nonlinear diffusion coefficient, are required. An equivalent skew-symmetric formulation of the convection and reaction terms of the nonlinear partial differential equation allows to avoid flux limitation and nonetheless leads to an unconditionally stable and convergent numerical method. Numerical experiments underline the robustness of the proposed algorithm