Lattice Boltzmann simulation of immiscible displacement in the cavity with different channel configurations

Abstract

In this work, the immiscible displacement in a cavity with different channel configurations is studied using an improved pseudo-potential lattice Boltzmann equation (LBE) model. This model overcomes the drawback of the dependence of the fluid properties on the grid size, which exists in the original pseudo-potential LBE model. The approach is first validated by the Laplace law. Then, it is employed to study the immiscible displacement process. The influences of different factors, such as the surface wettability, the distance between the gas cavity and liquid cavity and the surface roughness of the channel are investigated. Numerical results show that the displacement efficiency increases and the displacement time decreases with the increase of the surface contact angle. On the other hand, the displacement efficiency increases with increasing distance between the gas cavity and the liquid cavity at first and finally reaches a constant value. As for the surface roughness, two structures (a semicircular cavity and a semicircular bulge) are studied. The comprehensive results show that although the displacement processes for both the structures depend on the surface wettability, they present quite different behaviors. Specially, for the roughness structure constituted by the semicircular cavity, the displacement efficiency decreases and displacement time increases evidently with the size of the semicircular cavity for the small contact angle. The trend slows down as the increase of the contact angle. Once the contact angle exceeds a certain value, the size of the semicircular cavity almost has no influence on the displacement process. While for the roughness structure of a semicircular bulge, the displacement efficiency increases with the size of bulge first and then it decreases for the small contact angle. The displacement efficiency increases first and finally reaches a constant for the large contact angle. The results also show that the displacement time has an extreme value in these cases for the small contact angles.