Dynamics and breakup of a contracting liquid filament

Abstract

Contraction of a filament of an incompressible Newtonian liquid in a passive ambient fluid is studied computationally to provide insights into the dynamics of satellite drops created during drop formation. This free boundary problem, which is composed of the Navier–Stokes system and the associated initial and boundary conditions that govern the evolution in time of the filament shape and the velocity and pressure fields within it, is solved by the method of lines incorporating the finite element method for spatial discretization. The finite element algorithm developed here utilizes an adaptive elliptic mesh generation technique that is capable of tracking the dynamics of the filament up to the incipience of pinch-off without the use of remeshing. The correctness of the algorithm is verified by demonstrating that its predictions accord with (a) previously published results of Basaran (1992) on the analysis of finite-amplitude oscillations of viscous drops, (b) simulations of the dynamics of contracting filaments carried out with the well-benchmarked algorithm of Wilkes et al. (1999), and (c) scaling laws governing interface rupture and transitions that can occur from one scaling law to another as pinch-off is approached. In dimensionless form, just two parameters govern the problem: the dimensionless half-length $L_o$ and the Ohnesorge number $Oh$ which measures the relative importance of viscous force to capillary force. Regions of the parameter space are identified where filaments (a) contract to a sphere without breaking into multiple droplets, (b) break via the so-called endpinching mechanism where daughter drops pinch-off from the ends of the main filament, and (c) break after undergoing a series of complex oscillations. Predictions made with the new algorithm are also compared to those made with a model based on the slender-jet approximation. A region of the parameter space is found where the slender-jet approximation fares poorly, and its cause is elucidated by examination of the vorticity dynamics and flow fields within contracting filaments.