Search

This paper investigates a model of fluid flow in a fractured porous medium under the assumption of a uniform distribution of fractures throughout the volume. This model is based on the use of a fractional differential analogue of Darcy's law, as well as on the assumption that the properties of rock and fluid depend on pressure and its fractional derivative. Unlike previous studies, this study uses a fractional derivative in the Caputo-Fabrizio sense with a non-singular kernel. In this paper, we propose a numerical method for solving this initial boundary value problem and theoretically investigate the order of its convergence. The formulation of a fully discrete scheme is based on application of the finite difference approximation for integer and fractional timederivatives, and the Galerkin method in the spatial variable. A second-order formula is used to approximate both integer derivative and the fractional derivative in the sense of Caputo-Fabrizio. A priori estimates are obtained for both semi-discrete and fully discrete schemes, which imply their second-order convergence in time and space variables. A number of computational experiments were carried out on the example of a model problem to validate the accuracy of the scheme. The results of the numerical tests fully confirm the outcome of the theoretical analysis.

This paper is aimed at efficient numerical implementation of the fractional-order generalization of the stochastic Stokes–Darcy model, which has important scientific, applied, and economic significance in hydrology, the oil industry, and biomedicine. The essence of this generalization of the stochastic model is the introduction of fractional time derivatives in the sense of Caputo’s definition to take into account long-term changes in the properties of media. An efficient numerical method for the implementation of the fractional-order Stokes–Darcy model is proposed, which is based on the use of a higher-order approximation formula for the fractional derivative, higher-order finite difference relations, and a finite element approximation of the problem in the spatial direction. In the paper, a rigorous theoretical analysis of the stability and convergence of the proposed numerical method is carried out, which is confirmed by numerous computational experiments. Further, the proposed method is applied to the implementation of the fractional-order stochastic Stokes–Darcy model using an ensemble technique, in which the approximation is carried out in such a way that the resulting systems of linear equations have the same coefficient matrix for all realizations. Furthermore, evaluation of the discrete fractional derivatives is carried out with the use of parallel threads. The efficiency of applying both approaches has been demonstrated in numerical tests.

The paper studies a model of fluid flow in a fractured porous medium in which fractures are distributed uniformly over the volume. This model includes a nonlinear equation containing several terms with fractional derivatives in the sense of Caputo of order belonging to the interval 1,2. The relevance of studying this problem is determined by its practical significance in the oil industry, since most of the world’s oil reserves are in these types of reservoirs. The uniqueness of the solution to the problem in a differential form and its dependence on the initial data and the right-hand side of the equation is proved. A numerical method is proposed based on the use of the finite difference approximation for integer and fractional time derivatives and the finite element method in the spatial direction. A change of variables is introduced to reduce the order of the fractional derivatives. Furthermore, the fractional derivative is approximated by using the L1-method. The stability and convergence of the proposed numerical method are rigorously proved. The theoretical order of convergence is confirmed by the results of numerical tests for a problem of fluid flow in fractured porous media with a known exact solution.

In this paper, we propose an efficient numerical method for solving an initial boundary value problem for a coupled system of equations consisting of a nonlinear parabolic partial integro-differential equation and an elliptic equation with a nonlinear term. This problem has an important applied significance in petroleum engineering and finds application in modeling two-phase nonequilibrium fluid flows in a porous medium with a generalized nonequilibrium law. The construction of the numerical method is based on employing the finite element method in the spatial direction and the finite difference approximation to the time derivative. Newton’s method and the second-order approximation formula are applied for the treatment of nonlinear terms. The stability and convergence of the discrete scheme as well as the convergence of the iterative process is rigorously proven. Numerical tests are conducted to confirm the theoretical analysis. The constructed method is applied to study the two-phase nonequilibrium flow of an incompressible fluid in a porous medium. In addition, we present two examples of models allowing for prediction of the behavior of a fluid flow in a porous medium that are reduced to solving the nonlinear integro-differential equations studied in the paper.

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.