We perform an analytic and numerical study of an inviscid contracting bubble in a two-dimensional Hele-Shaw cell, where the effects of both surface tension and kinetic undercooling on the moving bubble boundary are not neglected. In contrast to expanding bubbles, in which both boundary effects regularize the illposedness arising from the viscous (Saffman-Taylor) instability, we show that in contracting bubbles the two boundary effects are in competition, with surface tension stabilizing the boundary, and kinetic undercooling destabilizing it. This competition leads to interesting bifurcation behaviour in the asymptotic shape of the bubble in the limit it approaches extinction. In this limit, the boundary may tend to become either circular, or approach a line or 'slit' of zero thickness, depending on the initial condition and the value of a nondimensional surface tension parameter. We show that over a critical range of surface tension values, both these asymptotic shapes are stable. In this regime there exists a third, unstable branch of limiting self-similar bubble shapes, with an asymptotic aspect ratio (dependent on the surface tension) between zero and one. We support our asymptotic analysis with a numerical scheme that utilizes the applicability of complex variable theory to Hele-Shaw flow.

Consider the dynamics of a healing film driven by surface tension, that is, the inward spreading process of a liquid film to fill a hole. The film is modelled using the lubrication (or thin-film) approximation, which results in a fourth-order nonlinear partial differential equation. We obtain a self-similar solution describing the early-time relaxation of an initial step-function condition and a family of self-similar solutions governing the finite-time healing. The similarity exponent of this family of solutions is not determined purely from scaling arguments; instead, the scaling exponent is a function of the finite thickness of the prewetting film, which we determine numerically. Thus, the solutions that govern the finite-time healing are self-similar solutions of the second kind. Laboratory experiments and time-dependent computations of the partial differential equation are also performed. We compare the self-similar profiles and exponents with both measurements in experiments and time-dependent computations near the healing time, and we observe good agreement in each case. . ‡ Z. Zheng and M. A. Fontelos contributed equally to this work.

We examine the effect of a kinetic undercooling condition on the evolution of a free boundary in Hele-Shaw flow, in both bubble and channel geometries. We present analytical and numerical evidence that the bubble boundary is unstable and may develop one or more corners in finite time, for both expansion and contraction cases. This loss of regularity is interesting because it occurs regardless of whether the less viscous fluid is displacing the more viscous fluid, or vice versa. We show that small contracting bubbles are described to leading order by a well-studied geometric flow rule. Exact solutions to this asymptotic problem continue past the corner formation until the bubble contracts to a point as a slit in the limit. Lastly, we consider the evolving boundary with kinetic undercooling in a Saffman-Taylor channel geometry. The boundary may either form corners in finite time, or evolve to a single long finger travelling at constant speed, depending on the strength of kinetic undercooling. We demonstrate these two different behaviours numerically. For the travelling finger, we present results of a numerical solution method similar to that used to demonstrate the selection of discrete fingers by surface tension. With kinetic undercooling, a continuum of corner-free travelling fingers exists for any finger width above a critical value, which goes to zero as the kinetic undercooling vanishes. We have not been able to compute the discrete family of analytic solutions, predicted by previous asymptotic analysis, because the numerical scheme cannot distinguish between solutions characterised by analytic fingers and those which are corner-free but non-analytic.

We study a simplified model of ice-ocean interaction beneath a floating ice shelf, and investigate the possibility for channels to form in the ice shelf base due to spatial variations in conditions at the grounding line. The model combines an extensional thin-film description of viscous ice flow in the shelf, with melting at its base driven by a turbulent ocean plume. Small transverse perturbations to the one-dimensional steady state are considered, driven either by ice thickness or subglacial discharge variations across the grounding line. Either forcing leads to the growth of channels downstream, with melting driven by locally enhanced ocean velocities, and thus heat transfer. Narrow channels are smoothed out due to turbulent mixing in the ocean plume, leading to a preferred wavelength for channel growth. In the absence of perturbations at the grounding line, linear stability analysis suggests that the one dimensional state is stable to initial perturbations, chiefly due to the background ice advection.