Numerical simulation of bubble rising in porous media using lattice Boltzmann method

Rising bubble systems in porous media exist in a variety of industrial processes. However, the flow characteristics of the issue are not well understood. In this work, the rising of bubble/bubbles through two types of porous structures, namely, in-line structured pore and staggered structured pore, are studied using a large density ratio lattice Boltzmann model. The effects of Eötvös number, pore shape, viscosity ratio, initial bubble number, and arrangement manner of the initial bubbles on the bubble deformation, bubble rising velocity, residual bubble mass, bubble perimeter, and the number of bubble breakups are investigated. It is found that as the Eötvös number increases, the bubbles are more easily broken during the process of passing through the porous media, the shapes of the sub-bubbles deviate from the original ones more and more, the bubble perimeter increases, and the difference between the bubble dynamics obtained by the in-line and staggered porous media decreases. Compared to the results of circular and rectangular pores, the bubble rising through the diamondoid pore has a more considerable deformation, which causes a slower rising speed. Furthermore, in the case that two bubbles are originally placed under the porous medium, the bubble deformation is greater and the bubble fracture times increase if the initial bubbles are aligned vertically. The findings of this work can contribute to the understanding of gas–liquid two-phase flow in porous media.