Theory and experiments on the stagnant cap regime in the motion of spherical surfactant-laden bubbles

The buoyant motion of a bubble rising through a continuous liquid phase can be retarded by the adsorption onto the bubble surface of surfactant dissolved in the liquid phase. The reason for this retardation is that adsorbed surfactant is swept to the trailing pole of the bubble where it accumulates and lowers the surface tension relative to the front end. The difference in tension creates a Marangoni force which opposes the surface flow, rigidifies the interface and increases the drag coefficient. Surfactant molecules adsorb onto the bubble surface by diffusing from the bulk to the sublayer of liquid adjoining the surface, and kinetically adsorbing from the sublayer onto the surface. The surface surfactant distribution which defines the Marangoni force is determined by the rate of kinetic adsorption and bulk diffusion relative to the rate of surface convection. In the limit in which the rate of either kinetic or diffusive transport of surfactant to the bubble surface is slow relative to surface convection and surface diffusion is also slow, surfactant collects in a stagnant cap at the back end of the bubble while the front end is stress free and mobile. The size of the cap and correspondingly the drag coefficient increases with the bulk concentration of surfactant until the cap covers the entire surface and the drag coefficient is that of a bubble with a completely tangentially immobile surface. Previous theoretical research on the stagnant cap regime has not studied in detail the competing roles of bulk diffusion and kinetic adsorption in determining the size of the stagnant cap angle, and there have been only a few studies which have attempted to quantitatively correlate simulations with measurements.This paper provides a more complete theoretical study of and a validating set of experiments on the stagnant cap regime. We solve numerically for the cap angle and drag coefficient as a function of the bulk concentration of surfactant for a spherical bubble rising steadily with inertia in a Newtonian fluid, including both bulk diffusion and kinetic adsorption. For the case of diffusion-limited transport (infinite adsorption kinetics), we show clearly that very small bulk concentrations can immobilize the entire surface, and we calculate the critical concentrations which immobilize the surface as a function of the surfactant parameters. We demonstrate that the effect of kinetics is to reduce the cap angle (hence reduce the drag coefficient) for a given bulk concentration of surfactant. We also present experimental results on the drag of a bubble rising in a glycerol–water mixture, as a function of the dissolved concentration of a polyethoxylated non-ionic surfactant whose bulk diffusion coefficient and a lower bound on the kinetic rate constants have been obtained separately by measuring the reduction in dynamic tension as surfactant adsorbs onto a clean interface. For low concentrations of surfactant, the experiments measure drag coefficients which are intermediate between the drag coefficient of a bubble whose surface is tangentially mobile and one whose surface is completely immobilized. Using the separately obtained value for the diffusion coefficient of the polyethoxylate, we undertake simulations which provide, upon comparison with the measured drag coefficients, a tighter bound on the kinetic rate constants than were otherwise obtained using dynamic surface tension measurement, and this suggests a new method for the measurement of kinetic rate constants.