Low-Reynolds-number gravity-driven migration and deformation of bubbles near a free surface

Abstract

We investigate numerically the axisymmetric migration of bubbles toward a free surface, using a boundary-integral technique. Our careful numerical implementation allows to study the bubble(s) deformation and film drainage; it is benchmarked against several tests. The rise of one bubble toward a free surface is studied and the computed bubble shape compared with the results of Princen [J. Colloid Interface Sci. 18, 178 (1963)]. The liquid film between the bubble and the free surface is found to drain exponentially in time in full agreement with the experimental work of Debrégeas et al. [Science 279, 1704 (1998)]. Our numerical results also cast some light on the role played by the deformation of the fluid interfaces and it turns out that for weakly deformed interfaces (high surface tension or a tiny bubble) the film drainage is faster than for a large fluid deformation. By introducing one or two additional bubble(s) below the first one, we examine to which extent the previous trends are affected by bubble-bubble interactions. For instance, for a 2-bubble chain, decreasing the bubble-bubble separation increases the deformation of the last bubble in the chain. Finally, the exponential drainage of the film between the free surface and the closest bubble is preserved, yet the drainage is enhanced.