Surface tension causes the edge of a fluid sheet to retract. If the sheet is also stretched along its edge then the flow and the rate of retraction are modified. A universal similarity solution for the Stokes flow in a stretched edge shows that the scaled shape of the edge is independent of the stretching rate, and that it decays exponentially to its far-field thickness. This solution justifies the use of a stress boundary condition in long-wavelength models of stretched viscous sheets, and gives the detailed shape of the edge of such a sheet, resolving the position of the sheet edge to the order of the thickness.

The coalescence of bubbles and drops plays a central role in nature and industry. During coalescence, two bubbles or drops touch and merge into one as the neck connecting them grows from microscopic to macroscopic scales. The hydrodynamic singularity that arises when two bubbles or drops have just touched and the flows that ensue have been studied thoroughly when two drops coalesce in a dynamically passive outer fluid. In this paper, the coalescence of two identical and initially spherical bubbles, which are idealized as voids that are surrounded by an incompressible Newtonian liquid, is analyzed by numerical simulation. This problem has recently been studied (a) experimentally using high-speed imaging and (b) by asymptotic analysis in which the dynamics is analyzed by determining the growth of a hole in the thin liquid sheet separating the two bubbles. In the latter, advantage is taken of the fact that the flow in the thin sheet of nonconstant thickness is governed by a set of one-dimensional, radial extensional flow equations. While these studies agree on the power law scaling of the variation of the minimum neck radius with time, they disagree with respect to the numerical value of the prefactors in the scaling laws. In order to reconcile these differences and also provide insights into the dynamics that are difficult to probe by either of the aforementioned approaches, simulations are used to access both earlier times than has been possible in the experiments and also later times when asymptotic analysis is no longer applicable. Early times and extremely small length scales are attained in the new simulations through the use of a truncated domain approach. Furthermore, it is shown by direct numerical simulations in which the flow within the bubbles is also determined along with the flow exterior to them that idealizing the bubbles as passive voids has virtually no effect on the scaling laws relating minimum neck radius and time.

When two spherical bubbles touch, a hole is formed in the fluid sheet between them, and capillary pressure acting on its tightly curved edge drives an outward radial flow which widens the hole joining the bubbles. Recent images of the early stages of this process (Paulsen et al., Nat. Commun., vol. 5, 2014) show that the radius of the hole $r_{\!E}$ at time $t$ grows proportional to $t^{1/2}$ , and that the rate is dependent on the fluid viscosity. Here, we explain this behaviour in terms of similarity solutions to a third-order system of radial extensional-flow equations for the thickness and velocity of the sheet of fluid between the bubbles, and determine the growth rate as a function of the Ohnesorge number $\mathit{Oh}$ . The initially quadratic sheet profile allows the ratio of viscous and inertial effects to be independent of time. We show that the sheet is slender for $r_{\!E}\ll a$ if $\mathit{Oh}\gg 1$ , where $a$ is the bubble radius, but only slender for $r_{\!E}\ll \mathit{Oh}^{2}a$ if $\mathit{Oh}\ll 1$ due to a compressional boundary layer of length $L\propto \mathit{Oh}\,r_{\!E}$ , after which there is a change in the structure but not the speed of the retracting sheet. For $\mathit{Oh}\ll 1$ , the detailed analysis justifies a simple momentum-balance argument, which gives the analytic prediction $r_{\!E}\sim (32a{\it\gamma}/3{\it\rho})^{1/4}t^{1/2}$ , where ${\it\gamma}$ is the surface tension and ${\it\rho}$ is the density.