Investigation on spontaneous liquid–liquid imbibition in capillaries with varying axial geometries using lattice Boltzmann method

Abstract

Spontaneous liquid–liquid imbibition in capillaries with irregular axial geometries is common in the petroleum industry. Monitoring the real-time dynamic contact angle (DCA) of the meniscus is crucial during such processes. In this work, we extend the Bell–Cameron–Lucas–Washburn (BCLW) equation by considering the axial shape of the capillaries, inertial force, and non-wetting fluid viscosity. We also develop a cascaded multi-component Shan–Chen lattice Boltzmann method (CLBM) with a modified mass-conservative curved boundary scheme to accurately simulate imbibition processes in sinusoidal capillaries. The results indicate that the DCA is highly sensitive to variations in the axial geometry of the capillary during imbibition, displaying a periodic time evolution pattern. When the axial geometry diverges, the DCA increases, and when it converges, the DCA decreases. The viscosity ratio affects the imbibition velocity, controlling the evolution period and extreme values of the DCA. A critical contact angle exists for a fixed capillary axial geometry and viscosity ratio. Continuous spontaneous imbibition occurs if the static contact angle is smaller than this critical value. However, if it exceeds this threshold, imbibition ceases within regions where axial geometry divergence. Moreover, we noticed a discrepancy in imbibition lengths predicted by the extended BCLW equation that ignores the DCA compared to those computed through the CLBM. To address this issue, we employed CLBM to monitor the DCA in real time and used the gathered data to refine the extended BCLW equation. As a result, the prediction of imbibition lengths by the extended BCLW equation for coupling the DCA became more accurate.