Breakup of finite-size liquid filaments: Transition from no-breakup to breakup including substrate effects⋆.

This work studies the breakup of finite-size liquid filaments, when also including substrate effects, using direct numerical simulations. The study focuses on the effects of three parameters: Ohnesorge number, the ratio of the viscous forces to inertial and surface tension forces, the liquid filament aspect ratio, and where there is a substrate, a measure of the fluid slip on the substrate, i.e. slip length. Through these parameters, it is determined whether a liquid filament breaks up during the evolution toward its final equilibrium state. Three scenarios are identified: a collapse into a single droplet, the breakup into one or multiple droplets, and recoalescence into a single droplet after the breakup (or even possibly another breakup after recoalescence). The results are compared with the ones available in the literature for free-standing liquid filaments. The findings show that the presence of the substrate promotes the breakup of the filament. The effect of the degree of slip on the breakup is also discussed. The parameter domain regions are comprehensively explored when including the slip effects. An experimental case is also carried out to illustrate the collapse and breakup of a finite-size silicon oil filament supported on a substrate, showcasing a critical length of the breakup in a physical configuration. Finally, direct numerical simulations reveal striking new details into the breakup pattern for low Ohnesorge numbers, where the dynamics are fast and the experimental imaging is not available; our results therefore significantly extend the range of Ohnesorge number over which filament breakup has been considered.

Continuum simulation of the discharge of the granular silo: a validation test for the µ(I) visco-plastic flow law.

Using a continuum Navier-Stokes solver with the µ(I) flow law implemented to model the viscous behavior, and the discrete Contact Dynamics algorithm, the discharge of granular silos is simulated in two dimensions from the early stages of the discharge until complete release of the material. In both cases, the Beverloo scaling is recovered. We first do not attempt a quantitative comparison, but focus on the qualitative behavior of velocity and pressure at different locations in the flow. A good agreement for the velocity is obtained in the regions of rapid flows, while areas of slow creep are not entirely captured by the continuum model. The pressure field shows a general good agreement, while bulk deformations are found to be similar in both approaches. The influence of the parameters of the µ(I) flow law is systematically investigated, showing the importance of the dependence on the inertial number I to achieve quantitative agreement between continuum and discrete discharge. However, potential problems involving the systems size, the configuration and "non-local" effects, are suggested. Yet the general ability of the continuum model to reproduce qualitatively the granular behavior is found to be very encouraging.

Numerical simulation of electrospray in the cone-jet mode.

We present a robust and computationally efficient numerical scheme for simulating steady electrohydrodynamic atomization processes (electrospray). The main simplification assumed in this scheme is that all the free electrical charges are distributed over the interface. A comparison of the results with those calculated with a volume-of-fluid method showed that the numerical scheme presented here accurately describes the flow pattern within the entire liquid domain. Experiments were performed to partially validate the numerical predictions. The simulations reproduced accurately the experimental shape of the liquid cone jet, providing correct values of the emitted electric current even for configurations very close to the cone-jet stability limit.

Coalescence of liquid drops by surface tension.

The merging of two mercury drops at very low kinetic energy is observed using fast, digital, and analog imaging techniques. Sequences showing the time evolution of the overall-surface shape as well as an amplified view of the contact region are shown. Qualitative and quantitative comparisons with computations of the Navier-Stokes equation with a free surface are made. In the model, the surface is tracked by a marker-chain method.