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.

Законы Стокса и Дарси имеют широкое применение в области гидродинамики. За последние десятилетия многие значительные исследовательские усилия были сосредоточены на изучении взаимосвязанной модели Стокса-Дарси, чтобы получить более глубокое понимание явлений гидродинамики. В этом контексте были предложены и тщательно изучены различные типы условий на границе между подобластями и улучшенные модели. Более того, стохастическая модель Стокса-Дарси стала ценным инструментом для учета неопределенностей и уточнения нашего понимания этих процессов. В этом обзоре мы проанализируем классическую и стохастические модели Стокса-Дарси, стремясь всесторонне изучить их преимущества и недостатки, а также некоторые численные методы. Мы фокусируемся на важности производных дробного порядка в моделях гидродинамики и анализируем преимущества новой обобщенной стохастической модели Стокса-Дарси дробного порядка.