Steady displacement of long gas bubbles in channels and tubes filled by a Bingham fluid

Bingham fluids behave like solids below a von Mises stress threshold, the yield stress, while above it they behave like Newtonian fluids. They are characterized by a dimensionless parameter, Bingham number (Bn), which is the ratio of the yield stress to a characteristic viscous stress. In this study, the noninertial steady motion of a finite-size gas bubble in both a plane two-dimensional (2D) channel and an axisymmetric tube filled by a Bingham fluid has been studied numerically. The Bingham number, Bn, is in the range 0 <= Bn <= 3, where Bn = 0 is the Newtonian case, while the capillary number, which is the ratio of a characteristic viscous force to the surface tension, has values Ca = 0.05, 0.10, and 0.25. The volume of all axisymmetric and 2D bubbles has been chosen to be identical for all parameter choices and large enough for the bubbles to be long compared to the channel, tube height, and diameter. The Bingham fluid constitutive equation is approximated by a regularized equation. During the motion, the bubble interface is separated from the wall by a static liquid film. The film thickness scaled by the tube radius (axisymmetric) and half of the channel height (2D) is the dimensionless film thickness, h. The results show that increasing Bn initially leads to an increase in h; however, the profile h versus Bn can be monotonic or nonmonotonic depending on Ca values and 2D or axisymmetric configurations. The yield stress also alters the shape of the front and rear of the bubble and suppresses the capillary waves at the rear of the bubble. The yield stress increases the magnitude of the wall shear stress and its gradient and therefore increases the potential for epithelial cell injuries in applications to lung airway mucus plugs. The topology of the yield surfaces as well as the flow pattern in the bubble frame of reference varies significantly by Ca and Bn.