Numerical Simulation of Ostwald Ripening of Bubbles in Porous Media Based on a Three-Dimensional Pore Network Model

Zhang Muan,Wang Jinqing,Wu Rui,Feng Zhi,Zhan Mingxiu,Xu Xu,Chi Zuohe,,

DOI: https://doi.org/10.7498/aps.72.20230695

IF: 0.906

2023-01-01

Acta Physica Sinica

Abstract:Ostwald ripening of bubbles in porous media is a common phenomenon observed in various fields, including CO 2 geological storage, preparation of porous materials, and fuel cells. A three-dimensional pore network model based on concentration coupling calculation has been developed to investigate the ripening characteristics of bubbles in porous media at the pore scale. This model takes into account the shape of bubbles, the structure of porous media, and mass transfer between gas and liquid. By solving the gas phase concentration of each pore body in the three-dimensional pore network, the model can track the evolution process of each bubble. A microfluidic chip with a four-pore structure was used to validate the reliability of the model through visual experiments. To analyze the effect of porous media heterogeneity on the bubble ripening process, the model constructed two different three-dimensional pore network structures with varying pore sizes and numerically simulated the ripening process of bubbles in two regions. The results show that the initial distribution of bubbles can affect the ripening process of porous media. When bubbles are uniformly distributed, during the ripening process, they exhibit regular and systematic changes in their spacing. However, in the case of uneven bubble distribution, as the bubbles transfer from smaller pore regions to larger pore regions, they also undergo individual mass transfer towards the larger bubble regions within their respective areas. Consequently, the remaining bubbles no longer maintain a spaced distribution pattern. Additionally, differences in the initial size of bubbles can accelerate the ripening process, resulting in a significantly shorter ripening time compared to a uniform distribution. The choice of pore number has a significant impact on continuous scale equivalent parameters, such as average capillary pressure and saturation. As the number of pores increases, the capillary pressure and saturation exhibit a more regular, nonlinear variation. A relationship between capillary pressure and saturation in the small and large pore regions has been established, which deviates from the assumptions made in existing literature. This result provides important guidance for the construction of continuous-scale ripening models. This model can be used to predict the evolution process of CO 2 during geological storage and provide guidance for studying the impact mechanism of heterogeneity in the long-term storage of CO 2 .

physics, multidisciplinary