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.

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.

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.

This paper presents an efficient numerical method for the fractional-order generalization of the stochastic Stokes–Darcy model, which finds application in various engineering, biomedical and environmental problems involving interaction between free fluid flow and flows in porous media. Unlike the classical model, this model allows taking into account the hereditary properties of the process under uncertainty conditions. The proposed numerical method is based on the combined use of the sparse grid stochastic collocation method, finite element/finite difference discretization, a fast numerical algorithm for computing the Caputo fractional derivative, and a cost-effective ensemble strategy. The hydraulic conductivity tensor is assumed to be uncertain in this problem, which is modeled by the reduced Karhunen–Loève expansion. The stability and convergence of the deterministic numerical method have been rigorously proved and validated by numerical tests. Utilizing the ensemble strategy allowed us to solve the deterministic problem once for all samples of the hydraulic conductivity tensor, rather than solving it separately for each sample. The use of the algorithm for computing the fractional derivatives significantly reduced both computational cost and memory usage. This study also analyzes the influence of fractional derivatives on the fluid flow process within the fractional-order Stokes–Darcy model under uncertainty conditions.

The paper focuses on constructing and investigation of cost-effective difference schemes for the numerical solution of the two-dimensional three-phase non-isothermal flow problem without capillary and gravitational forces. In this paper, the finite difference method is used to solve the problem numerically. The method of energy inequalities is applied to examine the stability of the finite difference scheme with respect to the initial data and the right-hand sides of the equations. Three cost-effective difference schemes are constructed on the base of the studied scheme. The efficiency of the proposed algorithms is analyzed on the basis of comparing the average time spent on the numerical implementation of one time layer.

This article presents an algorithm for the numerical solution of an initial-boundary value problem for a symmetric t-hyperbolic system of partial differential equations. This problem is based on continual filtration model, which describes the propagation of seismic waves in a poroelastic medium saturated with a fluid characterized by such physical parameters as the propagation velocities of longitudinal P- (fast and slow) and transverse S-waves, the density of the medium materials, and porosity. The system of linearized equations of saturated porous media is formulated in terms of physical variables of the velocity–stress tensor of the porous matrix and the velocity–pressure of the saturating fluid in the absence of energy dissipation. The solution is implemented numerically using an explicit finite difference upwind scheme built on a staggered grid to avoid the appearance of oscillations in the solution functions. The program code implementing parallel computing is developed in the high-performance Julia programming language. The possibility of using the approach is demonstrated by the example of solving the problem of propagation of seismic waves from a source located in the formation. Computational experiments based on real data from oil reservoirs have been implemented, and dynamic visualization of solutions consistent with the first waves arrival times has been obtained.