Stationary Stokes solver for single-phase flow in porous media: A blastingly fast solution based on Algebraic Multigrid Method using GPU

The paper is focused on high efficiency Stokes solver that is applied to the incompressible flow in porous media. Computational domains are represented by binarized 3D computed tomography voxel models. The solution procedure is constructed around a fast algebraic multigrid solver that utilizes the p...

Full description

Saved in:
Bibliographic Details
Published in:Advances in water resources Vol. 171; p. 104340
Main Authors: Evstigneev, Nickolay M., Ryabkov, Oleg I., Gerke, Kirill M.
Format: Journal Article
Language:English
Published: Elsevier Ltd 01.01.2023
Subjects:
Algebraic Multigrid Method
computed tomography
domain
General purpose GPU
Krylov subspace solvers
memory
Permeability
Permeability tensor
porosity
Porous media
Stokes flow
topology
water
Porous media
Permeability tensor
Krylov subspace solvers
Permeability
General purpose GPU
Stokes flow
Algebraic Multigrid Method
ISSN:0309-1708
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:The paper is focused on high efficiency Stokes solver that is applied to the incompressible flow in porous media. Computational domains are represented by binarized 3D computed tomography voxel models. The solution procedure is constructed around a fast algebraic multigrid solver that utilizes the power of graphics processing units (GPUs). In order to minimize memory footprint and accelerate the solver, a simple MAC-type staggered finite difference discretization is used and a coupled Stokes saddle-point type system is solved directly on a GPU. The MAC discretization is discussed, taking particular care of internal solid boundaries around pores. The method includes topological domain analysis on a GPU, which removes isolated volumes with no flow in the domain (based on the connected component labeling), depending on the boundary conditions. We consider various types of boundary conditions and efficient parallel strategies for the GPUs, including fast matrix assembly and residual regularization. Our method is extensively benchmarked both against analytical solutions and applied problems from digital rock physics in terms of computational wall time, precision, approximation order, and convergence. We demonstrate that it takes up to 5–23 s on a modern GPU card to obtain a solution with 1⋅10−6 error residuals for 3D geometries with 300–4503 voxels and a porosity range of 5–37%. [Display omitted] •Coupled second order Stokes solver for pore-scale modeling is presented.•The solution for 4503 voxels is obtained within 20 s on a modern GPU card.•The solver uses approximate commutator, multigrid preconditioners, Krylov methods.•The topological domain analysis removes null subspaces in linear operators.•The software is implemented on modern GPU architecture; it is extensively tested.
Bibliography:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 23
ISSN:0309-1708
DOI:10.1016/j.advwatres.2022.104340