A level set method for simulating capillary‐controlled displacements at the pore scale with nonzero contact angles

Jettestuen, Espen; Helland, Johan Olav; Prodanović, Maša

doi:10.1002/wrcr.20334

Cited by 104 publications

(83 citation statements)

References 43 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…The fundamental idea of the LBM is to construct simplified kinetic models that incorporate the essential physics of microscopic or mesoscopic processes to ensure that the macroscopic averaged properties obey the desired macroscopic equations. LBM has several advantages over the conventional grid-based computational fluid dynamics (CFD) methods, such as volume-of-fluid (VOF) 17,18 and level-set (LS) methods, 19,20 especially in dealing with complex boundaries, incorporation of microscopic interactions, flexible reproduction of interface between multiple phases, and parallelisation of the algorithm. In the LBM community, a number of multiphase models have been proposed.…”

Section: Introductionmentioning

confidence: 99%

Lattice Boltzmann simulation of immiscible fluid displacement in porous media: Homogeneous versus heterogeneous pore network

2015

View full text Add to dashboard Cite

Injection of anthropogenic carbon dioxide (CO 2 ) into geological formations is a promising approach to reduce greenhouse gas emissions into the atmosphere. Predicting the amount of CO 2 that can be captured and its long-term storage stability in subsurface requires a fundamental understanding of multiphase displacement phenomena at the pore scale. In this paper, the lattice Boltzmann method is employed to simulate the immiscible displacement of a wetting fluid by a non-wetting one in two microfluidic flow cells, one with a homogeneous pore network and the other with a randomly heterogeneous pore network. We have identified three different displacement patterns, namely, stable displacement, capillary fingering, and viscous fingering, all of which are strongly dependent upon the capillary number (Ca), viscosity ratio (M), and the media heterogeneity. The non-wetting fluid saturation (S nw ) is found to increase nearly linearly with log Ca for each constant M. Increasing M (viscosity ratio of non-wetting fluid to wetting fluid) or decreasing the media heterogeneity can enhance the stability of the displacement process, resulting in an increase in S nw . In either pore networks, the specific interfacial length is linearly proportional to S nw during drainage with equal proportionality constant for all cases excluding those revealing considerable viscous fingering. Our numerical results confirm the previous experimental finding that the steady state specific interfacial length exhibits a linear dependence on S nw for either favorable (M ≥ 1) or unfavorable (M < 1) displacement, and the slope is slightly higher for the unfavorable displacement. C 2015 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License. [http://dx

show abstract

Section: Introductionmentioning

confidence: 99%

Lattice Boltzmann simulation of immiscible fluid displacement in porous media: Homogeneous versus heterogeneous pore network

2015

View full text Add to dashboard Cite

show abstract

“…To compute multiphase flow properties at the pore scale several direct numerical simulation methods have been developed using the volume-of-fluid method ( Raeini et al, 2015;Shams et al, 2018 ), the level set method ( Jettestuen et al, 2013;Prodanovi ć and Bryant, 2006 ) and the lattice Boltzmann method (LBM) ( Boek et al, 2017;Leclaire et al, 2017;Ramstad et al, 2012;Tölke et al, 2006 ). Several studies have used the LBM to compute flow through porespace images ( Ahrenholz et al, 2008;Boek et al, 2017;Leclaire et al, 2017;Ramstad et al, 2012;.…”

Section: Introductionmentioning

confidence: 99%

Wetting boundary condition for the color-gradient lattice Boltzmann method: Validation with analytical and experimental data

Akai

Bijeljic

Blunt

2018

Advances in Water Resources

102

View full text Add to dashboard Cite

a b s t r a c tIn the color gradient lattice Boltzmann model (CG-LBM), a fictitious-density wetting boundary condition has been widely used because of its ease of implementation. However, as we show, this may lead to inaccurate results in some cases. In this paper, a new scheme for the wetting boundary condition is proposed which can handle complicated 3D geometries. The validity of our method for static problems is demonstrated by comparing the simulated results to analytical solutions in 2D and 3D geometries with curved boundaries. Then, capillary rise simulations are performed to study dynamic problems where the three-phase contact line moves. The results are compared to experimental results in the literature (Heshmati and Piri, 2014). If a constant contact angle is assumed, the simulations agree with the analytical solution based on the Lucas-Washburn equation. However, to match the experiments, we need to implement a dynamic contact angle that varies with the flow rate.

show abstract

“…The progressive quasistatic algorithm based on the level set method was applied in an analysis of the matrix-fracture transfer and the shape of the two-phase interface during drainage and imbibition in low permeability rock [19,20]. The level set method was also employed to simulate quasistatic drainage and imbibition for different contact angles in porous media [21]. The level set method and phase field method were used to simulate twophase flow with viscosity contrast through complex porous media using COMSOL software [22].…”

Section: Introductionmentioning

confidence: 99%

Modeling of Two-Phase Flow in Rough-Walled Fracture Using Level Set Method

et al. 2017

View full text Add to dashboard Cite

To describe accurately the flow characteristic of fracture scale displacements of immiscible fluids, an incompressible two-phase (crude oil and water) flow model incorporating interfacial forces and nonzero contact angles is developed. The roughness of the two-dimensional synthetic rough-walled fractures is controlled with different fractal dimension parameters. Described by the Navier-Stokes equations, the moving interface between crude oil and water is tracked using level set method. The method accounts for differences in densities and viscosities of crude oil and water and includes the effect of interfacial force. The wettability of the rough fracture wall is taken into account by defining the contact angle and slip length. The curve of the invasion pressure-water volume fraction is generated by modeling two-phase flow during a sudden drainage. The volume fraction of water restricted in the rough-walled fracture is calculated by integrating the water volume and dividing by the total cavity volume of the fracture while the two-phase flow is quasistatic. The effect of invasion pressure of crude oil, roughness of fracture wall, and wettability of the wall on two-phase flow in rough-walled fracture is evaluated.

show abstract

A level set method for simulating capillary‐controlled displacements at the pore scale with nonzero contact angles

Cited by 104 publications

References 43 publications

Lattice Boltzmann simulation of immiscible fluid displacement in porous media: Homogeneous versus heterogeneous pore network

Lattice Boltzmann simulation of immiscible fluid displacement in porous media: Homogeneous versus heterogeneous pore network

Wetting boundary condition for the color-gradient lattice Boltzmann method: Validation with analytical and experimental data

Modeling of Two-Phase Flow in Rough-Walled Fracture Using Level Set Method

Contact Info

Product

Resources

About