Numerical investigation of the formation and stability of homogeneous pairs of soft particles in inertial microfluidics

Numerical investigation of heterogeneous soft particle pairs in inertial microfluidics

The formation of pairs of particles or cells of different types in microfluidic channels can be desired or detrimental in healthcare applications. It is still unclear what role softness heterogeneity plays in the formation of these particle pairs. We use an in-house lattice-Boltzmann-immersed-boundary-finite-element solver to simulate a pair of particles with different softness flowing through a straight channel with a rectangular cross-section under initial conditions representative of a dilute suspension. We find that softness heterogeneity significantly affects the pair dynamics, determining whether a pair will form or not, and determining the lateral and inter-particle equilibrium behaviour in the pair. We also observe close matches between the transient deformation of particles in a linear pair and single particles in isolation. These results further our understanding of pair behaviour, providing a foundation for understanding particle train formation, and open up the potential to develop reduced-order models for particle pair formation based upon the behaviour of single particles.

Numerical investigation of the dynamics of a rigid spherical particle in a vortical cross-slot flow at moderate inertia

The study of flow and particle dynamics in microfluidic cross-slot channels is of high relevance for lab-on-a-chip applications. In this work, we investigate the dynamics of a rigid spherical particle in a cross-slot junction for a channel height-to-width ratio of 0.6 and at a Reynolds number of 120 for which a steady vortex exists in the junction area. Using an in-house immersed-boundary-lattice-Boltzmann code, we analyse the effect of the entry position of the particle in the junction and the particle size on the dynamics and trajectory shape of the particle. We find that the dynamics of the particle depend strongly on its lateral entry position in the junction and weakly on its vertical entry position; particles that enter close to the centre show trajectory oscillations. Larger particles have longer residence times in the junction and tend to oscillate less due to their confinement. Our work contributes to the understanding of particle dynamics in intersecting flows and enables the design of optimised geometries for cytometry and particle manipulation.

Fluid interfaces laden by force dipoles: towards active matter-driven microfluidic flows

In recent years, nonlinear microfluidics in combination with lab-on-a-chip devices has opened a new avenue for chemical and biomedical applications such as droplet formation and cell sorting. In this article, we integrate ideas from active matter into a microfluidic setting, where two fluid layers with identical densities but different viscosities flow through a microfluidic channel. Most importantly, the fluid interface is laden with active particles that act with dipolar forces on the adjacent fluids and thereby generate flows. We perform lattice-Boltzmann simulations and combine them with phase field dynamics of the interface and an advection-diffusion equation for the density of active particles. We show that only contractile force dipoles can destabilize the flat fluid interface. It develops a viscous finger from which droplets break up. For interfaces with non-zero surface tension, a critical value of activity equal to the surface tension is necessary to trigger the instability. Since activity depends on the density of force dipoles, the interface can develop steady deformation. Lastly, we demonstrate how to control droplet formation using switchable activity.

Rheology and structure of elastic capsule suspensions within rectangular channels

Three-dimensional simulations of the pressure-driven flow dynamics of elastic capsule suspensions within both slit and rectangular cross-section channels are presented. The simulations utilize the Immersed Boundary Method and the Lattice-Boltzmann Method models. The capsule volume fraction is fixed at 0.1 (i.e., a semi-dilute suspension), while the channel Reynolds number (Re), the capillary number (Ca), and the cross-sectional channel dimensions are systematically varied. Comparing results for slit and rectangular channels, it is found that multi-directional confinement hinders inertial focusing due to the capsule-free layers that develop in the two transverse directions. Furthermore, the thicknesses of the capsule-free layers in the two transverse directions differ when the height and width of the channel are not equal. Both the size and aspect ratio of the channel impact the apparent viscosity. It is found that square channels exhibit maximal viscosity and that holding one dimension fixed while increasing or decreasing the other results in a decrease in viscosity. The results therefore represent an expansion of the Fahraeus-Lindqvist effect from 1D cylindrical channels to 2D rectangular channels.

Microparticle separation using dielectrophoresis-assisted inertial microfluidics: A GPU-accelerated immersed boundary-lattice Boltzmann simulation

In this study, the migration of microparticles towards the inertial equilibrium positions in a straight microchannel with a square cross section in the presence of an inhomogeneous oscillating electric field was examined. The dynamics of microparticles were simulated using the immersed boundary-lattice Boltzmann method of fluid-structure interaction simulation. Moreover, the lattice Boltzmann Poisson solver was applied to calculate the electric field required for calculation of the dielectrophoretic force using the equivalent dipole moment approximation. These numerical methods were implemented on a single GPU coupled with the AA pattern of storing distribution functions in memory to speed up the computationally demanding simulation of microparticles dynamics. In the absence of an electric field, spherical polystyrene microparticles migrate to four symmetric stable equilibrium positions corresponding to the sidewalls of the square cross-sectional microchannel. The equilibrium distance from the sidewall was increased by increasing the particle size. The equilibrium positions near electrodes disappeared and particles migrated to the other equilibrium positions far from the electrodes by the application of the high-frequency oscillatory electric field at voltages beyond a threshold value. Finally, a two-step dielectrophoresis-assisted inertial microfluidics methodology was introduced for particle separation based on the crossover frequencies and the observed threshold voltages of different particles. The proposed method exploited the synergistic effect of dielectrophoresis and inertial microfluidics methods to remove their limitations, allowing the separation of a broad range of polydisperse particle mixtures with a single device in a short time.