The present paper is devoted to the construction and study of numerical methods for solving an initial boundary value problem for a differential equation containing several terms with fractional time derivatives in the sense of Caputo. This equation is suitable for describing the process of fluid flow in fractured porous media under some physical assumptions, and has an important applied significance in petroleum engineering. Two different approaches to constructing numerical schemes depending on orders of the fractional derivatives are proposed. The semi-discrete and fully discrete numerical schemes for solving the problem are analyzed. The construction of a fully discrete scheme is based on applying the finite difference approximation to time derivatives and the finite element method in the spatial direction. The approximation of the fractional derivatives in the sense of Caputo is carried out using the L1-method. The convergence of both numerical schemes is rigorously proved. The results of numerical tests conducted for model problems are provided to confirm the theoretical analysis. In addition, the proposed computational method is applied to study the flow of oil in a fractured porous medium within the framework of the considered model. Based on the results of the numerical tests, it was concluded that the model reproduces the characteristic features of the fluid flow process in the medium under consideration.

We study an elliptic equation in a regular domain with a condition on the boundary involving a generalized Riemann–Liouville derivative of fractional order. The Bitsadze-Samarskii type problem is formulated for that equation. The uniqueness is proved by the maximum principle for harmonic functions. The Poisson kernel of the Dirichlet problem for the Laplace equation is used for proving the existence of a solution for the formulated problem.