scholarly journals Modeling Immiscible Fluid Displacement in a Porous Medium Using Lattice Boltzmann Method

Fluids ◽  
2021 ◽  
Vol 6 (2) ◽  
pp. 89
Author(s):  
Magzhan Atykhan ◽  
Bagdagul Kabdenova (Dauyeshova) ◽  
Ernesto Monaco ◽  
Luis R. Rojas-Solórzano

The numerical investigation of the interpenetrating flow dynamics of a gas injected into a homogeneous porous media saturated with liquid is presented. The analysis is undertaken as a function of the inlet velocity, liquid–gas viscosity ratio (D) and physical properties of the porous medium, such as porous geometry and surface wettability. The study aims to improve understanding of the interaction between the physical parameters involved in complex multiphase flow in porous media (e.g., CO2 sequestration in aquifers). The numerical simulation of a gaseous phase being introduced through a 2D porous medium constructed using seven staggered columns of either circular- or square-shaped micro-obstacles mimicking the solid walls of the pores is performed using the multiphase Lattice Boltzmann Method (LBM). The gas–liquid fingering phenomenon is triggered by a small geometrical asymmetry deliberately introduced in the first column of obstacles. Our study shows that the amount of gas penetration into the porous medium depends on surface wettability and on a set of parameters such as capillary number (Ca), liquid–gas viscosity ratio (D), pore geometry and surface wettability. The results demonstrate that increasing the capillary number and the surface wettability leads to an increase in the effective gas penetration rate, disregarding porous medium configuration, while increasing the viscosity ratio decreases the penetration rate, again disregarding porous medium configuration.

Author(s):  
Eslam Ezzatneshan ◽  
Reza Goharimehr

In the present study, a pore-scale multicomponent lattice Boltzmann method (LBM) is employed for the investigation of the immiscible-phase fluid displacement in a homogeneous porous medium. The viscous fingering and the stable displacement regimes of the invading fluid in the medium are quantified which is beneficial for predicting flow patterns in pore-scale structures, where an experimental study is extremely difficult. Herein, the Shan-Chen (S-C) model is incorporated with an appropriate collision model for computing the interparticle interaction between the immiscible fluids and the interfacial dynamics. Firstly, the computational technique is validated by a comparison of the present results obtained for different benchmark flow problems with those reported in the literature. Then, the penetration of an invading fluid into the porous medium is studied at different flow conditions. The effect of the capillary number (Ca), dynamic viscosity ratio (M), and the surface wettability defined by the contact angle (θ) are investigated on the flow regimes and characteristics. The obtained results show that for M<1, the viscous fingering regime appears by driving the invading fluid through the pore structures due to the viscous force and capillary force. However, by increasing the dynamic viscosity ratio and the capillary number, the invading fluid penetrates even in smaller pores and the stable displacement regime occurs. By the increment of the capillary number, the pressure difference between the two sides of the porous medium increases, so that the pressure drop Δp along with the domain at θ=40∘ is more than that of computed for θ=80∘. The present study shows that the value of wetting fluid saturation Sw at θ=40∘ is larger than its value computed with θ=80∘ that is due to the more tendency of the hydrophilic medium to absorb the wetting fluid at θ=40∘. Also, it is found that the magnitude of Sw computed for both the contact angles is decreased by the increment of the viscosity ratio from Log(M)=−1 to 1. The present study demonstrates that the S-C LBM is an efficient and accurate computational method to quantitatively estimate the flow characteristics and interfacial dynamics through the porous medium.


2020 ◽  
Vol 10 (1) ◽  
Author(s):  
Krzysztof M. Graczyk ◽  
Maciej Matyka

AbstractConvolutional neural networks (CNN) are utilized to encode the relation between initial configurations of obstacles and three fundamental quantities in porous media: porosity ($$\varphi$$ φ ), permeability (k), and tortuosity (T). The two-dimensional systems with obstacles are considered. The fluid flow through a porous medium is simulated with the lattice Boltzmann method. The analysis has been performed for the systems with $$\varphi \in (0.37,0.99)$$ φ ∈ ( 0.37 , 0.99 ) which covers five orders of magnitude a span for permeability $$k \in (0.78, 2.1\times 10^5)$$ k ∈ ( 0.78 , 2.1 × 10 5 ) and tortuosity $$T \in (1.03,2.74)$$ T ∈ ( 1.03 , 2.74 ) . It is shown that the CNNs can be used to predict the porosity, permeability, and tortuosity with good accuracy. With the usage of the CNN models, the relation between T and $$\varphi$$ φ has been obtained and compared with the empirical estimate.


1998 ◽  
Vol 09 (08) ◽  
pp. 1491-1503 ◽  
Author(s):  
David M. Freed

An extension to the basic lattice-BGK algorithm is presented for modeling a simulation region as a porous medium. The method recovers flow through a resistance field with arbitrary values of the resistance tensor components. Corrections to a previous algorithm are identified. Simple validation tests are performed which verify the accuracy of the method, and demonstrate that inertial effects give a deviation from Darcy's law for nominal simulation velocities.


Author(s):  
Haijing Li ◽  
Herman J. H. Clercx ◽  
Federico Toschi

A model based on the Lattice Boltzmann method is developed to study the flow of reactive electro-kinetic fluids in porous media. The momentum, concentration and electric/potential fields are simulated via the Navier–Stokes, advection–diffusion/Nernst–Planck and Poisson equations, respectively. With this model, the total density and velocity fields, the concentration of reactants and reaction products, including neutral and ionized species, the electric potential and the interaction forces between the fields can be studied, and thus we provide an insight into the interplay between chemistry, flow and the geometry of the porous medium. The results show that the conversion efficiency of the reaction can be strongly influenced by the fluid velocity, reactant concentration and by porosity of the porous medium. The fluid velocity determines how long the reactants stay in the reaction areas, the reactant concentration controls the amount of the reaction material and with different dielectric constant, the porous medium can distort the electric field differently. All these factors make the reaction conversion efficiency display a non-trivial and non-monotonic behaviour as a function of the flow and reaction parameters. To better illustrate the dependence of the reaction conversion efficiency on the control parameters, based on the input from a number of numerical investigations, we developed a phenomenological model of the reactor. This model is capable of capturing the main features of the causal relationship between the performance of the reactor and the main test parameters. Using this model, one could optimize the choice of reaction and flow parameters in order to improve the performance of the reactor and achieve higher production rates. This article is part of the theme issue ‘Progress in mesoscale methods for fluid dynamics simulation’.


2017 ◽  
Vol 20 (10) ◽  
pp. 899-919 ◽  
Author(s):  
Sajjad Foroughi ◽  
Mohsen Masihi ◽  
Saeid Jamshidi ◽  
Mahmoud Reza Pishvaie

Author(s):  
Wojciech Sobieski

AbstractThe paper describes the so-called Waterfall Algorithm, which may be used to calculate a set of parameters characterising the spatial structure of granular porous media, such as shift ratio, collision density ratio, consolidation ratio, path length and minimum tortuosity. The study is performed for 1800 different two-dimensional random pore structures. In each geometry, 100 individual paths are calculated. The impact of porosity and the particle size on the above-mentioned parameters is investigated. It was stated in the paper, that the minimum tortuosity calculated by the Waterfall Algorithm cannot be used directly as a representative tortuosity of pore channels in the Kozeny or the Carman meaning. However, it may be used indirect by making the assumption that a unambiguous relationship between the representative tortuosity and the minimum tortuosity exists. It was also stated, that the new parameters defined in the present study are sensitive on the porosity and the particle size and may be therefore applied as indicators of the geometry structure of granular media. The Waterfall Algorithm is compared with other methods of determining the tortuosity: A-Star Algorithm, Path Searching Algorithm, Random Walk technique, Path Tracking Method and the methodology of calculating the hydraulic tortuosity based on the Lattice Boltzmann Method. A very short calculation time is the main advantage of the Waterfall Algorithm, what meant, that it may be applied in a very large granular porous media.


Sign in / Sign up

Export Citation Format

Share Document