Inertial migration of deformable droplets in a microchannel

The microfluidic inertial effect is an effective way of focusing and sorting droplets suspended in a carrier fluid in microchannels. To understand the flow dynamics of microscale droplet migration, we conduct numerical simulations on the droplet motion and deformation in a straight microchannel. The results are compared with preliminary experiments and theoretical analysis. In contrast to most existing literature, the present simulations are three-dimensional and full length in the streamwise direction and consider the confinement effects for a rectangular cross section. To thoroughly examine the effect of the velocity distribution, the release positions of single droplets are varied in a quarter of the channel cross section based on the geometrical symmetries. The migration dynamics and equilibrium positions of the droplets are obtained for different fluid velocities and droplet sizes. Droplets with diameters larger than half of the channel height migrate to the centerline in the height direction and two equilibrium positions are observed between the centerline and the wall in the width direction. In addition to the well-known Segré-Silberberg equilibrium positions, new equilibrium positions closer to the centerline are observed. This finding is validated by preliminary experiments that are designed to introduce droplets at different initial lateral positions. Small droplets also migrate to two equilibrium positions in the quarter of the channel cross section, but the coordinates in the width direction are between the centerline and the wall. The equilibrium positions move toward the centerlines with increasing Reynolds number due to increasing deformations of the droplets. The distributions of the lift forces, angular velocities, and the deformation parameters of droplets along the two confinement direction are investigated in detail. Comparisons are made with theoretical predictions to determine the fundamentals of droplet migration in microchannels. In addition, existence of the inner equilibrium position is linked to the quartic velocity distribution in the width direction through a simple model for the slip angular velocities of droplets.