Geometry Effects on Interfacial Dynamics of Gas‐Driven Drainage in a Gradient Capillary

Abstract Unfavorable fluid‐fluid displacement, where a low‐viscosity fluid displaces a higher‐viscosity fluid in permeable media, is commonly encountered in various subsurface processes. Understanding the formation and evolution of the resulting interfacial instability can have practical benefits for engineering applications. Using gradient capillary tubes as surrogate models of permeable media, we numerically investigate interfacial dynamics during gas‐driven drainage. Our focus is on understanding the impact of tube geometry on interface stability. In a gradient tube, since the interface shape changes during the drainage process, we measure interfacial stability using the difference between the contact‐line velocity U cl and the meniscus tip velocity U tip . We define instability as a rapid reduction in the contact line velocity U cl compared to the tip velocity U tip . Beyond the onset of this instability, gas penetrates into the liquid, forming a finger, and entraining a liquid film on the tube wall. The observed stability transition can be rationalized to a large extent by adaptation of an existing theory for cylindrical tubes in terms of a critical capillary number Ca crit . For an expanding tube, simulations suggest that a stability transition from an initially unstable meniscus to a final stable one, with U cl catching up with U tip , can occur if the local capillary number is initially slightly larger than Ca crit and then drops below Ca crit . The insights gained from this study can be beneficial in estimating the mode and efficiency of subsurface fluid displacement.