Instabilities and singularities in Hele–Shaw flow

Abstract

A mechanism by which smooth initial conditions evolve towards a topological reconfiguration of fluid interfaces is studied in the context of Darcy’s law. In the case of thin fluid layers, nonlinear PDEs for the local thickness are derived from an asymptotic limit of the vortex sheet representation. A particular example considered is the Rayleigh–Taylor instability of stratified fluid layers, where the instability of the system is controlled by a Bond number B. It is proved that, for a range of B and initial data “subharmonic” to it, interface pinching must occur in at least infinite time. Numerical simulations suggest that “pinching” singularities occur generically when the system is unstable, and in particular immediately above a bifurcation point to instability. Near this bifurcation point an approximate analytical method describing the approach to a finite-time singularity is developed. The method exploits the separation of time scales that exists close to the first instability in a system of finite extent, with a discrete spectrum of modes. In this limit, slowly growing long-wavelength modes entrain faster short-wavelength modes, and thereby, allow the derivation of a nonlinear evolution equation for the amplitudes of the slow modes. The initial-value problem is solved in this slaved dynamics, yielding the time and analytical structure of a singularity that is associated with the motion of zeros in the complex plane, suggesting a general mechanism of singularity formation in this system. The discussion emphasizes the significance of several variational principles, and comparisons are made between the numerical simulations and the approximate theory.