Abstract
Recently, there has been an increased interest in the problem of numerical implementation of multiphase filtration models due to its enormous economic importance in the oil industry, hydrology, and nuclear waste management. In contrast to the classical models of filtration, filtration models in highly porous fractured formations with the fractal geometry of wells are not fully understood. The solution to this problem reduces to solving a system of differential equations with fractional derivatives. In the paper, a finite-difference scheme is constructed for solving the initial-boundary value problem for the convection-diffusion equation with a fractional time derivative in the sense of Caputo-Fabrizio. A priori estimates are obtained for solving a difference problem under the assumption that there is a solution to the problem in the class of sufficiently smooth functions that prove the uniqueness of the solution and the stability of the difference scheme. The convergence of the solution of the difference problem to the solution of the original differential problem with the second order in time and space variables is shown. The results of computational experiments confirming the reliability of theoretical analysis are presented.