Open Access

ARTICLE

Influence of Fractal Dimension on Gas-Driven Two-Phase Flow in Fractal Porous Media: A VOF Model-Based Simulation

Xiaolin Wang, Richeng Liu^*, Kai Qiu, Zhongzhong Liu, Shisen Zhao, Shuchen Li

State Key Laboratory of Intelligent Construction and Healthy Operation and Maintenance of Deep Underground Engineering, China University of Mining and Technology, Xuzhou, 221116, China

* Corresponding Author: Richeng Liu. Email:

Computer Modeling in Engineering & Sciences 2025, 144(1), 289-307. https://doi.org/10.32604/cmes.2025.066716

Received 15 April 2025; Accepted 30 June 2025; Issue published 31 July 2025

Abstract

Gas-liquid two-phase flow in fractal porous media is pivotal for engineering applications, yet it remains challenging to be accurately characterized due to complex microstructure-flow interactions. This study establishes a pore-scale numerical framework integrating Monte Carlo-generated fractal porous media with Volume of Fluid (VOF) simulations to unravel the coupling among pore distribution characterized by fractal dimension (D_f), flow dynamics, and displacement efficiency. A pore-scale model based on the computed tomography (CT) microstructure of Berea sandstone is established, and the simulation results are compared with experimental data. Good agreement is found in phase distribution, breakthrough behavior, and flow path morphology, confirming the reliability of the numerical simulation method. Ten fractal porous media models with D_f ranging from 1.25~1.7 were constructed using a Monte-Carlo approach. The gas-liquid two-phase flow dynamics was characterized using the VOF solver across gas injection rates of 0.05–5 m/s, in which the time-resolved two-phase distribution patterns were systematically recorded. The results reveal that smaller fractal dimensions (D_f = 1.25~1.45) accelerate fingering breakthrough (peak velocity is 1.73 m/s at D_f = 1.45) due to a bimodal pore size distribution dominated by narrow channels. Increasing D_f amplifies vorticity generation by about 3 times (eddy viscosity is 0.033 Pa·s at D_f = 1.7) through reduced interfacial curvature, while tortuosity-driven pressure differentials transition from sharp increases (0.4~6.3 Pa at D_f = 1.25~1.3) to inertial plateaus (4.8 Pa at D_f = 1.7). A nonlinear increase in equilibrium gas volume fraction (f_av = 0.692 at D_f = 1.7) emerges from residual gas saturation and turbulence-enhanced dispersion. This behavior is further modulated by flow velocity, with f_av peaking at 0.72 under capillary-dominated conditions (0.05 m/s), but decreasing to 0.65 in the inertial regime (0.5 m/s). The work quantitatively links fractal topology to multiphase flow regimes, demonstrating the critical role of D_f in governing preferential pathways, energy dissipation, and phase distribution.

---

Keywords

---