Deformation and breakup of a liquid droplet past a solid circular cylinder: a lattice Boltzmann study

Impact Dynamics of a Single Droplet on Hydrophobic Cylinders: A Lattice Boltzmann Study

This study numerically investigates the effects of the Weber number (We) and cylinder-to-droplet radius ratio (R^*) on the impact dynamics of a low-viscosity droplet on a hydrophobic cylinder by the lattice Boltzmann method. The intrinsic contact angle of the surface is chosen as θ₀ = 122°± 2°, which ensures a representative hydrophobicity. The regime diagram of the impact dynamics in the parameter space of We versus R^* is established with categories of split and nonsplit regimes. The droplet would split during impact as α = We/R^* exceeds a critical value. In the nonsplit regime, the droplet bounces off the cylinder at most Weber numbers unless the impact velocity is minuscule (We < 2). The contact time of the droplet on the cylinder surface decreases with increasing R^* or decreasing We, indicating bouncing is facilitated under such conditions. This can be explained by the suppressed adhesion dissipation between the droplet and surface due to a reduction in the contact area. In the split regime, sufficient kinetic energy inside the impacting droplet determines whether the whole droplet could detach from the surface. With a small cylinder (R^* < 0.83) and large We (>25), the adhesion effect is weakened for the side fragments because of the small contact area, and it facilitates the dripping of fragments. For other conditions, the detachment, especially for the tiny droplet on the cylinder top, only occurs if the deformation is prominent at We > 35. Moreover, the spreading dynamics of the impacting droplet are also highlighted in this work.

On-Demand Droplet Merging with an AC Electric Field for Multiple-Volume Droplet Generation

We introduce a unique system to achieve on-demand droplet merging and splitting using a perpendicular AC electric field. The working mechanism involves a micropillar to split droplets, followed by electrocoalescence using an AC electric field. Adjusting the parameters of the AC signal and conductivity of the fluid result in different merging regimes. We observed a minimum threshold voltage and a strong influence of the surfactant. We hypothesize that the merging process is caused by dipole-dipole coalescence between the daughter droplets. At the same time, adjustment of the conductivity reveals a shift in the merging regimes and can be explained with an electric circuit diagram. Size-based sorting using this merging phenomenon is subsequently demonstrated, where alternate, single, double, and triple droplets sorting were achieved. The concept presented in this paper is potentially useful for drug dispensing or multivolume digital polymerase chain reaction, as droplets of multiple sizes can be generated simultaneously.

Lattice Boltzmann Simulation on Droplet Flow through 3D Metal Foam

The hydrodynamics of droplets passing through metal foam is investigated using the lattice Boltzmann method (LBM). The accurate 3D porous structure for the simulation is generated by X-ray micro-computed tomography. The simulated results are in good agreement with the experimental ones using high-speed video. The simulated results show that for droplets passing metal foam, there is a critical capillary number, Cac (around 0.061), above which the droplet continues to deform until it breaks up. The simulated results show that the capillary number, droplet size, pores diameter, and thickness of metal foam have the significant effect of droplets deforming and breaking up when the droplets pass through the metal foam. To avoid the calescence of two droplets at the inlet zone of the metal foam, the distance between droplets should be larger than three times the diameter of the droplet.

Microfluidic guillotine for single-cell wound repair studies

Wound repair is a key feature distinguishing living from nonliving matter. Single cells are increasingly recognized to be capable of healing wounds. The lack of reproducible, high-throughput wounding methods has hindered single-cell wound repair studies. This work describes a microfluidic guillotine for bisecting single Stentor coeruleus cells in a continuous-flow manner. Stentor is used as a model due to its robust repair capacity and the ability to perform gene knockdown in a high-throughput manner. Local cutting dynamics reveals two regimes under which cells are bisected, one at low viscous stress where cells are cut with small membrane ruptures and high viability and one at high viscous stress where cells are cut with extended membrane ruptures and decreased viability. A cutting throughput up to 64 cells per minute-more than 200 times faster than current methods-is achieved. The method allows the generation of more than 100 cells in a synchronized stage of their repair process. This capacity, combined with high-throughput gene knockdown in Stentor, enables time-course mechanistic studies impossible with current wounding methods.

Study of flow behaviors of droplet merging and splitting in microchannels using Micro-PIV measurement

Droplet merging and splitting are important droplet manipulations in droplet-based microfluidics. However, the fundamental flow behaviors of droplets were not systematically studied. Hence, we designed two different microstructures to achieve droplet merging and splitting respectively, and quantitatively compared different flow dynamics in different microstructures for droplet merging and splitting via micro-particle image velocimetry (micro-PIV) experiments. Some flow phenomena of droplets different from previous studies were observed during merging and splitting using a high-speed microscope. It was also found the obtained instantaneous velocity vector fields of droplets have significant influence on the droplets merging and splitting. For droplet merging, the probability of droplets coalescence (η) in a microgroove is higher (50% < η < 92%) than that in a T-junction microchannel (15% < η < 50%), and the highest coalescence efficiency (η = 92%) comes at the two-phase flow ratio e of 0.42 in the microgroove. Moreover, compared with a cylinder obstacle, Y-junction bifurcation can split droplets more effectively and the droplet flow during splitting is steadier. The results can provide better understanding of droplet behaviors and are useful for the design and applications of droplet-based microfluidics.

Lattice Boltzmann modeling of three-phase incompressible flows

In this paper, based on multicomponent phase-field theory we intend to develop an efficient lattice Boltzmann (LB) model for simulating three-phase incompressible flows. In this model, two LB equations are used to capture the interfaces among three different fluids, and another LB equation is adopted to solve the flow field, where a new distribution function for the forcing term is delicately designed. Different from previous multiphase LB models, the interfacial force is not used in the computation of fluid velocity, which is more reasonable from the perspective of the multiscale analysis. As a result, the computation of fluid velocity can be much simpler. Through the Chapman-Enskog analysis, it is shown that the present model can recover exactly the physical formulations for the three-phase system. Numerical simulations of extensive examples including two circular interfaces, ternary spinodal decomposition, spreading of a liquid lens, and Kelvin-Helmholtz instability are conducted to test the model. It is found that the present model can capture accurate interfaces among three different fluids, which is attributed to its algebraical and dynamical consistency properties with the two-component model. Furthermore, the numerical results of three-phase flows agree well with the theoretical results or some available data, which demonstrates that the present LB model is a reliable and efficient method for simulating three-phase flow problems.

Analytical relations for long-droplet breakup in asymmetric T junctions

We develop accurate analytical relations for the droplet volume ratio, droplet length during breakup process, and pressure drop of asymmetric T junctions with a valve in each of the branches for producing unequal-sized droplets. An important advantage of this system is that after manufacturing the system, the size of the generated droplets can be changed simply by adjusting the valves. The results indicate that if the valve ratio is smaller than 0.65, the system enters a nonbreakup regime. Also the pressure drop does not depend on the time and decreases by increasing the valve ratio, namely, opening the degree of valve 1 to valve 2. In addition, the results reveal that by decreasing (increasing) the valve ratio, the droplet length of branch 1 decreases (increases) and the droplet length of branch 2 increases (decreases) linearly while the whole length of the droplet remains unchanged.

Lattice Boltzmann investigation of droplet inertial spreading on various porous surfaces

The spreading of liquid drops on solid surfaces is a wide-spread phenomenon of both fundamental and industrial interest. In many applications, surfaces are porous and spreading patterns are very complex with respect to the case on smooth surfaces. Focusing on the inertial spreading just before the Tanner-like viscous regime, this work investigates the spreading of a low-viscosity droplet on a porous surface using lattice Boltzmann numerical simulations. The case of a flat surface is first considered, and it reveals a dependence on the solid equilibrium contact angle θ(s)(eq), which is in good agreement with published experimental data. We conducted numerical experiments with various surfaces perforated by a regular pattern of holes of infinite length. The results show that the global spreading dynamics is independent of the porosity morphology. Through the assumption that, for wetting, the pores can be regarded as surface patches with a contact angle of θ(pore)(eq)=180°, we deduce an effective equilibrium contact angle θ(eff)(eq) on the porous surface from the Cassie-Baxter law. A spreading model is then proposed to describe both a prefactor and an exponent that are similar to a flat surface whose equilibrium contact angle is θ(eff)(eq). This model compares satisfactorily with a large number of numerical experiments under varying conditions.