Abstract:This article presents a three-dimensional numerical framework for the simulation of fluid-fluid immiscible compounds in complex geometries, based on the multiple-relaxation-time lattice Boltzmann method to model the fluid dynamics and the color-gradient approach to model multicomponent flow interaction. New lattice weights for the lattices D3Q15, D3Q19, and D3Q27 that improve the Galilean invariance of the color-gradient model as well as for modeling the interfacial tension are derived and provided in the Appendix. The presented method proposes in particular an approach to model the interaction between the fluid compound and the solid, and to maintain a precise contact angle between the two-component interface and the wall. Contrarily to previous approaches proposed in the literature, this method yields accurate solutions even in complex geometries and does not suffer from numerical artifacts like nonphysical mass transfer along the solid wall, which is crucial for modeling imbibition-type problems. The article also proposes an approach to model inflow and outflow boundaries with the color-gradient method by generalizing the regularized boundary conditions. The numerical framework is first validated for three-dimensional (3D) stationary state (Jurin's law) and time-dependent (Washburn's law and capillary waves) problems. Then, the usefulness of the method for practical problems of pore-scale flow imbibition and drainage in porous media is demonstrated. Through the simulation of nonwetting displacement in two-dimensional random porous media networks, we show that the model properly reproduces three main invasion regimes (stable displacement, capillary fingering, and viscous fingering) as well as the saturating zone transition between these regimes. Finally, the ability to simulate immiscible two-component flow imbibition and drainage is validated, with excellent results, by numerical simulations in a Berea sandstone, a frequently used benchmark case used in this field, using a complex geometry that originates from a 3D scan of a porous sandstone. The methods presented in this article were implemented in the open-source PALABOS library, a general C++ matrix-based library well adapted for massive fluid flow parallel computation.

APPLICATION OF THE PRECONDITIONED CONJUGATE GRADIENTS METHOD IN THE SIMULATION OF RELATIVE PERMEABILITY PROPERTIES OF POROUS MEDIA

Inversion Calculation of Permeability Coefficient with the Multi-Scale Asymptotic Expansion Method

Application of percolation, critical-path, and effective-medium theories for calculation of two-phase relative permeability

An Adaptive Preconditioner for Three-Dimensional Single-Phase Compressible Flow in Highly Heterogeneous Porous Media

Generalized network modeling of capillary-dominated two-phase flow

Modeling flow in porous media with double porosity/permeability: Mathematical model, properties, and analytical solutions

A Novel Finite Element-Based Method for Predicting the Permeability of Heterogeneous and Anisotropic Porous Microstructures

Numerical simulation of transport in porous media: some problems from micro to macro scale

Simulation of the impregnation in the porous media by the Self- organized Gradient Percolation method

A NEW RELATIVE PERMEABILITY MODEL OF UNSATURATED POROUS MEDIA BASED ON FRACTAL THEORY

Comparison between Van Genuchten and Brooks-Corey Parameterizations in the Solution of Multiphase Problems in Rigid One-Dimensional Porous Media

Pore Network Modeling of Fibrous Porous Media of Uniform and Gradient Porosity

Relative Permeability of Porous Media with Nonuniform Pores

Finite Analytic Numerical Method for Two-Dimensional Fluid Flow in Heterogeneous Porous Media

Prediction Of Two-Phase Relative Permeability In Porqus Media Based On Network Modeling Of Lattice Gas Automata

Three-dimensional mixed-wet random pore-scale network modeling of two- and three-phase flow in porous media. I. Model description

Network-inspired versus Kozeny-Carman based permeability-porosity relations applied to Biot's poroelasticity model

High Viscosity Ratio Multicomponent Flow Simulations in Porous Media using a Pseudo-Potential Central Moment Lattice Boltzmann Method

Scale-dependent permeability and formation factor in porous media: Applications from percolation theory

Generalized three-dimensional lattice Boltzmann color-gradient method for immiscible two-phase pore-scale imbibition and drainage in porous media

Computing Relative Permeability and Capillary Pressure of Heterogeneous Rocks Using Realistic Boundary Conditions