Search

If the Stokes equations are properly discretized, it is well-known that the Schur complement matrix is spectrally equivalent to the identity matrix. Moreover, in the case of simple geometries, it is often observed that most of its eigenvalues are equal to one. These facts form the basis for the famous Uzawa and Krylov-Uzawa algorithms. However, in the case of complex geometries, the Schur complement matrix can become arbitrarily ill-conditioned having a significant portion of non-unit eigenvalues, which makes the established Uzawa preconditioner inefficient. In this article, we study the Schur complement formulation for the staggered finite-difference discretization of the Stokes problem in 3D CT images and synthetic 2D geometries. We numerically investigate the performance of the CG iterative method with the Uzawa and SIMPLE preconditioners and draw several conclusions. First, we show that in the case of low porosity, CG with the SIMPLE preconditioner converges faster to the discrete pressure and provides a more accurate calculation of sample permeability. Second, we show that an increase in the surface-to-volume ratio leads to an increase in the condition number of the Schur complement matrix, while the dependence is inverse for the Schur complement matrix preconditioned with the SIMPLE. As an explanation, we conjecture that the no-slip boundary conditions are the reason for non-unit eigenvalues of the Schur complement.

The paper presents a workflow for fast pore-scale simulation of single-phase flow in tight reservoirs typically characterized by low, multiscale porosity. Multiscale porosity implies that the computational domain contains porous voxels (unresolved porosity) in addition to pure fluid voxels. In this case, the Stokes-Brinkman equations govern the flow, with the Darcy term needed to account for the flow in the porous voxels. As the central part of our workflow, robust and efficient solvers for Stokes and Stokes-Brinkman equations are presented. The solvers are customized for low-porosity binary and multiclass images, respectively. Another essential component of the workflow is a preprocessing module for classifying images with respect to the connectivity of the multiscale pore space. Particularly, an approximation of the Stokes-Brinkman problem, namely, the Darcy problem, is investigated for the images that do not have pure fluid percolation paths. Thorough computational experiments demonstrate efficiency and robustness of the workflow for simulations on images from tight reservoirs. Raw files describing the used CT images are provided as supplementary materials to enable other researchers to use them.

In our work, we consider the classical density-based approach to the topology optimization. We propose to modify the discretized cost functional using a posteriori error estimator for the finite element method. It can be regarded as a new technique to prevent checkerboards. It also provides higher regularity of solutions and robustness of results.