Non-Fickian convection–diffusion models in porous media

Abstract

In this paper we propose a numerical scheme to approximate the solution of a non-Fickian coupled model that describes, e.g., miscible transport in porous media. The model is defined by a system of a quasilinear elliptic equation, which governs the fluid pressure, and a quasilinear integro-differential equation, which models the convection–diffusion transport process. The numerical scheme is based on a conforming piecewise linear finite element method for the discretization in space. The fully discrete approximations is obtained with an implicit–explicit method. Estimates for the continuous in time and the fully discrete methods are derived, showing that the numerical approximation for the concentrations and the pressure are second order convergent in a discrete L^2-norm and in a discrete H^1-norm, respectively.