Inertial migration and axial control of deformable capsules

Computational Fluid-Structure Interaction in Microfluidics

Micro elastofluidics is a transformative branch of microfluidics, leveraging the fluid-structure interaction (FSI) at the microscale to enhance the functionality and efficiency of various microdevices. This review paper elucidates the critical role of advanced computational FSI methods in the field of micro elastofluidics. By focusing on the interplay between fluid mechanics and structural responses, these computational methods facilitate the intricate design and optimisation of microdevices such as microvalves, micropumps, and micromixers, which rely on the precise control of fluidic and structural dynamics. In addition, these computational tools extend to the development of biomedical devices, enabling precise particle manipulation and enhancing therapeutic outcomes in cardiovascular applications. Furthermore, this paper addresses the current challenges in computational FSI and highlights the necessity for further development of tools to tackle complex, time-dependent models under microfluidic environments and varying conditions. Our review highlights the expanding potential of FSI in micro elastofluidics, offering a roadmap for future research and development in this promising area.

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.

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.

A numerical lift force analysis on the inertial migration of a deformable droplet in steady and oscillatory microchannel flows

Inertial migration of deformable particles has become appealing in recent years due to its numerous applications in microfluidics and biomedicine. The physics underlying the motion of these particles is contingent upon the presence of lift forces in microchannels. Therefore, in this work, we present a lift force analysis for such migration of a deformable droplet in steady and oscillatory flow regimes and identify the effects of varying Capillary number and oscillation frequency on its dynamics. We then propose an expression that mimics the lift force behavior in oscillatory flows accurately. Finally, we introduce a procedure to derive and predict a simple expression for the steady and averaged oscillatory lift for any given combination of Capillary number and oscillation frequency within a continuous range.

Steady State of Motion of Two Particles in Poiseuille Flow of Power-Law Fluid

The steady state of motion of two particles in Poiseuille flow of power-law fluid is numerically studied using the lattice Boltzmann method in the range of Reynolds number 20 ≤ Re ≤ 60, diameter ratio of two particles 0.125 ≤ β ≤ 2.4, and power-law index of the fluid 0.4 ≤ n ≤ 1.2. Some results are validated by comparing with other available results. The effects of Re, β, and n on the steady state of motion of two particles are discussed. The results show that, for two particles of the same diameter, the particle spacing l in the steady state is independent of n. In shear-thinning fluid, l increases rapidly at first and then slowly, finally approaching a constant for different Re. In shear-thickening fluid, although l tends to be stable in the end, the values of l after stabilization are different. For two particles of different sizes, l does not always reach a stable state, and whether it reaches a stable state depends on n. When the small particle is downstream, l increases rapidly at first and then slowly in shear-thickening fluid, but increases rapidly at first and then decreases slowly, finally approaching a constant in a shear-thinning fluid. In shear-thinning fluid, the larger n is, the smaller l is. In shear-thickening fluid, β has no effect on l in steady-state. When the large particle is downstream, l increases rapidly at first and then slowly in shear-thinning fluid but increases rapidly at first and then decreases in a shear-thickening fluid. The effect of n on l in the steady state is obvious. In shear-thinning fluid, l increases rapidly at first and then slowly, the larger Re is, the smaller l is. In shear- thickening fluid, l will reach a stable state.

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

We investigate the formation and stability of a pair of identical soft capsules in channel flow under mild inertia. We employ a combination of the lattice Boltzmann, finite element and immersed boundary methods to simulate the elastic particles in flow. Validation tests show excellent agreement with numerical results obtained by other research groups. Our results reveal new trajectory types that have not been observed for pairs of rigid particles. While particle softness increases the likelihood of a stable pair forming, the pair stability is determined by the lateral position of the particles. A key finding is that stabilisation of the axial distance occurs after lateral migration of the particles. During the later phase of pair formation, particles undergo damped oscillations that are independent of initial conditions. These damped oscillations are driven by a strong hydrodynamic coupling of the particle dynamics, particle inertia and viscous dissipation. While the frequency and damping coefficient of the oscillations depend on particle softness, the pair formation time is largely determined by the initial particle positions: the time to form a stable pair grows exponentially with the initial axial distance. Our results demonstrate that particle softness has a strong impact on the behaviour of particle pairs. The findings could have significant ramifications for microfluidic applications where a constant and reliable axial distance between particles is required, such as flow cytometry.

Instability of a liquid sheet with viscosity contrast in inertial microfluidics

Flows at moderate Reynolds numbers in inertial microfluidics enable high throughput and inertial focusing of particles and cells with relevance in biomedical applications. In the present work, we consider a viscosity-stratified three-layer flow in the inertial regime. We investigate the interfacial instability of a liquid sheet surrounded by a density-matched but more viscous fluid in a channel flow. We use linear stability analysis based on the Orr-Sommerfeld equation and direct numerical simulations with the lattice Boltzmann method (LBM) to perform an extensive parameter study. Our aim is to contribute to a controlled droplet production in inertial microfluidics. In the first part, on the linear stability analysis we show that the growth rate of the fastest growing mode [Formula: see text] increases with the Reynolds number [Formula: see text] and that its wavelength [Formula: see text] is always smaller than the channel width w for sufficiently small interfacial tension [Formula: see text]. For thin sheets we find the scaling relation [Formula: see text], where m is viscosity ratio and [Formula: see text] the sheet thickness. In contrast, for thicker sheets [Formula: see text] decreases with increasing [Formula: see text] or m due to the nearby channel walls. Examining the eigenvalue spectra, we identify Yih modes at the interface. In the second part on the LBM simulations, the thin liquid sheet develops two distinct dynamic states: waves traveling along the interface and breakup into droplets with bullet shape. For smaller flow rates and larger sheet thicknesses, we also observe ligament formation and the sheet eventually evolves irregularly. Our work gives some indication how droplet formation can be controlled with a suitable parameter set [Formula: see text].

Flow-driven competition between two capsules passing through a narrow pore

By incorporating a distance function into the finite element simulation, we investigate the flow-driven competition between two soft capsules passing through a narrow pore, employing the arbitrary Lagrangian-Eulerian formulation to satisfy the boundary conditions for fluid flow and capsule deformation. In our simulations, the motion and deformation of the capsules can be described in an intuitive manner, and the order in which capsules of different sizes pass through a pore can be clearly determined. Meanwhile, when the capsules are near the narrow pore, the change of the flow field is also very interesting and can be expressed intuitively. It is shown that, driven by the Poiseuille flow, the larger capsule has a stronger tendency to pass through the pore than the small one, which can be attributed to the greater resistance and the volume advantage of the larger capsule. In addition, we demonstrate that this tendency can be reversed by changing the inlet velocity and setting the initial position of the smaller capsule closer to the axis of the pore. And as long as the large one passes through first, the small one will offset the axis to the same orientation as the initial, while the large one always moves along the axis.

EFFECT OF ARTERIOLAR DISTENSIBILITY ON THE LATERAL MIGRATION OF LIQUID-FILLED MICROPARTICLES FLOWING IN A HUMAN ARTERIOLE

A promising advance of bioengineering consists in the development of micro-nanoparticles as drug delivery vehicles injected intravenously or intraarterialy for targeted treatment. Proficient functioning of drug carries is conditioned by a reliable prediction of pharmacokinetics in human as well as their dynamical behavior once injected in blood stream. In this study, we aim to provide a reliable numerical prediction of dynamical behavior of microparticles in human arteriole focusing on the crucial mechanism of lateral migration. The dynamical response of the microparticle upon blood flow and arteriolar distensibility is investigated by varying main controlling parameters: viscosity ratio, confinement and capillary number. The influence of the hyperelastic arteriolar wall is highlighted through comparison with an infinitely rigid arteriolar wall. The hydrodynamic interaction in a microparticle train is examined. Fluid–structure interaction is solved by the Arbitrary Lagrangian–Eulerian method using the COMSOL Multiphysics software.

A pair of particles in inertial microfluidics: effect of shape, softness, and position

Lab-on-a-chip devices based on inertial microfluidics have emerged as a promising technique to manipulate particles in a precise way. Inertial microfluidics exploits internal hydrodynamic forces and the mechanical structure of particles to achieve separation and focusing. The article focuses on the hydrodynamic interaction of two particles. This will help to develop an understanding of the dynamics of particle trains in inertial microfluidics, which are typical structures in multi-particle systems. We perform three-dimensional lattice Boltzmann simulations combined with the immersed boundary method to unravel the dynamics of various mono- and bi-dispersed pairs in inertial microfluidics. We study the influence of different starting positions for mono- and bi-dispersed pairs. We also change their deformability from relatively soft to rigid and choose spherical and biconcave particle shapes. The observed two-particle motions in the present work can be categorized into four types: stable pair, stable pair with damped oscillations, stable pair with bounded oscillations, and unstable pair. We show that stable pairs become unstable when increasing the particle stiffness. Furthermore, a pair with both capsules in the same channel half is more prone to become unstable than a pair with capsules in opposite channel halves.

Public-Health-Driven Microfluidic Technologies: From Separation to Detection

Separation and detection are ubiquitous in our daily life and they are two of the most important steps toward practical biomedical diagnostics and industrial applications. A deep understanding of working principles and examples of separation and detection enables a plethora of applications from blood test and air/water quality monitoring to food safety and biosecurity; none of which are irrelevant to public health. Microfluidics can separate and detect various particles/aerosols as well as cells/viruses in a cost-effective and easy-to-operate manner. There are a number of papers reviewing microfluidic separation and detection, but to the best of our knowledge, the two topics are normally reviewed separately. In fact, these two themes are closely related with each other from the perspectives of public health: understanding separation or sorting technique will lead to the development of new detection methods, thereby providing new paths to guide the separation routes. Therefore, the purpose of this review paper is two-fold: reporting the latest developments in the application of microfluidics for separation and outlining the emerging research in microfluidic detection. The dominating microfluidics-based passive separation methods and detection methods are discussed, along with the future perspectives and challenges being discussed. Our work inspires novel development of separation and detection methods for the benefits of public health.

Inertial migration of a deformable capsule in an oscillatory flow in a microchannel

Dynamics of a deformable capsule in an oscillatory flow of a Newtonian fluid in a microchannel has been studied numerically. The effects of oscillation frequency, capsule deformability, and channel flow rate have been explored by simulating the capsule within a microchannel. In addition, the simulation captures the effect of the type of imposed pressure oscillations on the migration pattern of the capsule. An oscillatory channel flow enables the focusing of extremely small biological particles by eliminating the need to design impractically long channels. The presented results show that the equilibrium position of the capsule changes not only by the addition of an oscillatory component to the pressure gradient but it also is influenced by the capsule deformability and channel flow rate. Furthermore, it has been shown that the amplitude of oscillation of capsules decreases as the channel flow rate and the rigidity of the capsule increases.

Optimal Control of Colloidal Trajectories in Inertial Microfluidics Using the Saffman Effect

In inertial microfluidics colloidal particles in a Poiseuille flow experience the Segré-Silberberg lift force, which drives them to specific positions in the channel cross section. An external force applied along the microchannel induces a cross-streamline migration to a new equilibrium position because of the Saffman effect. We apply optimal control theory to design the time protocol of the axial control force in order to steer a single particle as precisely as possible from a channel inlet to an outlet at a chosen target position. We discuss the influence of particle radius and channel length and show that optimal steering is cheaper than using a constant control force. Using a single optimized control-force protocol, we demonstrate that even a pulse of particles spread along the channel axis can be steered to a target and that particles of different radii can be separarted most efficiently.

Computational inertial microfluidics: a review

Since the discovery of inertial focusing in 1961, numerous theories have been put forward to explain the migration of particles in inertial flows, but a complete understanding is still lacking. Recently, computational approaches have been utilized to obtain better insights into the underlying physics. In particular, fundamental aspects of particle focusing inside straight and curved microchannels have been explored in detail to determine the dependence of focusing behavior on particle size, channel shape, and flow Reynolds number. In this review, we differentiate between the models developed for inertial particle motion on the basis of whether they are semi-analytical, Navier-Stokes-based, or built on the lattice Boltzmann method. This review provides a blueprint for the consideration of numerical solutions for modeling of inertial particle motion, whether deformable or rigid, spherical or non-spherical, and whether suspended in Newtonian or non-Newtonian fluids. In each section, we provide the general equations used to solve particle motion, followed by a tutorial appendix and specified sections to engage the reader with details of the numerical studies. Finally, we address the challenges ahead in the modeling of inertial particle microfluidics for future investigators.

A flowing pair of particles in inertial microfluidics

A flowing pair of particles in inertial microfluidics gives important insights into understanding and controlling the collective dynamics of particles like cells or droplets in microfluidic devices. They are applied in medical cell analysis and engineering. We study the dynamics of a pair of solid particles flowing through a rectangular microchannel using lattice Boltzmann simulations. We determine the inertial lift force profiles as a function of the two particle positions, their axial distance, and the Reynolds number. Generally, the profiles strongly differ between particles leading and lagging in flow and the lift forces are enhanced due to the presence of a second particle. At small axial distances, they are determined by viscous forces, while inertial forces dominate at large separations. We identify cross-streamline pairs as stable fixed points in the lift force profiles and argue that same-streamline configurations are only one-sided stable. Depending on the initial conditions, the two-particle lift forces in combination with the Poiseuille flow give rise to three types of unbound particle trajectories, called moving-apart, passing, and swapping, and one type of bound trajectory, where the particles perform damped oscillations towards the cross-stream line configuration. The damping rate scales with Reynolds number squared, since inertial forces are responsible for driving the particles to their steady-state positions.

Dynamics of a capsule flowing in a tube under pulsatile flow

We analyze numerically the behavior of a deformable micro-capsule confined in a pipe under a pulsatile flow. The capsule moves and is deformed by the action of a pulsatile flow inside the tube with a non-null mean velocity. This configuration can be found in the nature and in many bioengineering systems where artificial capsules are driven by micro-pumps through micro-channels. The capsule is considered as a thin hyperelastic membrane, which encloses an internal fluid. As it has been demonstrated in the literature, this model represents a wide range of artificial capsules, for example, the alginate-based capsules, typically used in bioengineering applications. A hybrid isogeometric finite element method and boundary element method based on a T-spline discretization and formulated in the time domain is used to solve the mechanical and hydrodynamical equations. The influence of the relative rigidity of the membrane, frequency and amplitude of the pulsatile flow is studied. Results show that the behavior of the capsule differs from steady flows and it depends strongly on the frequency of the flow and mechanical characteristic of the capsule.

Flow Behavior of Chain and Star Polymers and Their Mixtures

Star-shaped polymers show a continuous change of properties from flexible linear chains to soft colloids, as the number of arms is increased. To investigate the effect of macromolecular architecture on the flow properties, we employ computer simulations of single chain and star polymers as well as of their mixtures under Poiseuille flow. Hydrodynamic interactions are incorporated through the multi-particle collision dynamics (MPCD) technique, while a bead-spring model is used to describe the polymers. For the ultradilute systems at rest, the polymers are distributed uniformly in the slit channel, with a weak dependence on their number of arms. Once flow is applied, however, we find that the stars migrate much more strongly towards the channel center as the number of arms is increased. In the star-chain mixtures, we find a flow-induced separation between stars and chains, with the stars located in the channel center and the chains closer to the walls. In order to identify the origin of this flow-induced partitioning, we conduct additional simulations without hydrodynamic interactions, and find that the observed cross-stream migration originates from a combination of wall-induced hydrodynamic lift forces and viscoelastic effects. The results from our study give valuable insights for designing microfluidic devices for separating polymers based on their architecture.

Inertial manipulation of bubbles in rectangular microfluidic channels

Inertial microfluidics is an active field of research that deals with crossflow positioning of the suspended entities in microflows. Until now, the majority of the studies have focused on the behavior of rigid particles in order to provide guidelines for microfluidic applications such as sorting and filtering. Deformable entities such as bubbles and droplets are considered in fewer studies despite their importance in multiphase microflows. In this paper, we show that the trajectory of bubbles flowing in rectangular and square microchannels can be controlled by tuning the balance of forces acting on them. A T-junction geometry is employed to introduce bubbles into a microchannel and analyze their lateral equilibrium position in a range of Reynolds (1 < Re < 40) and capillary numbers (0.1 < Ca < 1). We find that the Reynolds number (Re), the capillary number (Ca), the diameter of the bubble (D[combining macron]), and the aspect ratio of the channel are the influential parameters in this phenomenon. For instance, at high Re, the flow pushes the bubble towards the wall while large Ca or D[combining macron] moves the bubble towards the center. Moreover, in the shallow channels, having aspect ratios higher than one, the bubble moves towards the narrower sidewalls. One important outcome of this study is that the equilibrium position of bubbles in rectangular channels is different from that of solid particles. The experimental observations are in good agreement with the performed numerical simulations and provide insights into the dynamics of bubbles in laminar flows which can be utilized in the design of flow based multiphase flow reactors.

Numerical-experimental observation of shape bistability of red blood cells flowing in a microchannel

Red blood cells flowing through capillaries assume a wide variety of different shapes owing to their high deformability. Predicting the realized shapes is a complex field as they are determined by the intricate interplay between the flow conditions and the membrane mechanics. In this work we construct the shape phase diagram of a single red blood cell with a physiological viscosity ratio flowing in a microchannel. We use both experimental in vitro measurements as well as 3D numerical simulations to complement the respective other one. Numerically, we have easy control over the initial starting configuration and natural access to the full 3D shape. With this information we obtain the phase diagram as a function of initial position, starting shape and cell velocity. Experimentally, we measure the occurrence frequency of the different shapes as a function of the cell velocity to construct the experimental diagram which is in good agreement with the numerical observations. Two different major shapes are found, namely croissants and slippers. Notably, both shapes show coexistence at low (<1 mm s^-1) and high velocities (>3 mm s^-1) while in-between only croissants are stable. This pronounced bistability indicates that RBC shapes are not only determined by system parameters such as flow velocity or channel size, but also strongly depend on the initial conditions.

Feedback control of photoresponsive fluid interfaces

Photoresponsive surfactants provide a unique microfluidic driving mechanism. Since they switch between two molecular shapes under illumination and thereby affect surface tension of fluid interfaces, Marangoni flow along the interface occurs. To describe the dynamics of the surfactant mixture at a planar interface, we formulate diffusion-advection-reaction equations for both surfactant densities. They also include adsorption from and desorption into the neighboring fluids and photoisomerization by light. We then study how the interface responds when illuminated by spots of light. Switching on a single light spot, the density of the switched surfactant spreads in time and assumes an exponentially decaying profile in steady state. Simultaneously, the induced radial Marangoni flow reverses its flow direction from inward to outward. We use this feature to set up specific feedback rules, which couple the advection velocities sensed at the light spots to their intensities. As a result two neighboring spots switch on and off alternately. Extending the feedback rule to light spots arranged on the vertices of regular polygons, we observe periodic switching patterns for even-sided polygons, where two sets of next-nearest neighbors alternate with each other. A triangle and pentagon also show regular oscillations, while heptagon and nonagon exhibit irregular oscillations due to frustration. While our findings are specific to the chosen set of parameters, they show how complex patterns at photoresponsive fluid interfaces emerge from simple feedback coupling.

Brownian motion near an elastic cell membrane: A theoretical study

Elastic confinements are an important component of many biological systems and dictate the transport properties of suspended particles under flow. In this paper, we review the Brownian motion of a particle moving in the vicinity of a living cell whose membrane is endowed with a resistance towards shear and bending. The analytical calculations proceed through the computation of the frequency-dependent mobility functions and the application of the fluctuation-dissipation theorem. Elastic interfaces endow the system with memory effects that lead to a long-lived anomalous subdiffusive regime of nearby particles. In the steady limit, the diffusional behavior approaches that near a no-slip hard wall. The analytical predictions are validated and supplemented with boundary-integral simulations.