Badr Kaoui

Noninertial lateral migration of vesicles in bounded Poiseuille flow

Cross-streamline noninertial migration of a vesicle in a bounded Poiseuille flow is investigated experimentally and numerically. The combined effects of the walls and of the curvature of the velocity profile induce a movement toward the center of the channel. A migration law (as a function of relevant structural and flow parameters) is proposed that is consistent with experimental and numerical results. This similarity law markedly differs from its analog in unbounded geometry. The dependency on the reduced volume ν and viscosity ratio λ is also discussed. In particular, the migration velocity becomes nonmonotonous as a function of ν beyond a certain λ.

Lateral migration of a two-dimensional vesicle in unbounded Poiseuille flow

The migration of a suspended vesicle in an unbounded Poiseuille flow is investigated numerically in the low Reynolds number limit. We consider the situation without viscosity contrast between the interior of the vesicle and the exterior. Using the boundary integral method we solve the corresponding hydrodynamic flow equations and track explicitly the vesicle dynamics in two dimensions. We find that the interplay between the nonlinear character of the Poiseuille flow and the vesicle deformation causes a cross-streamline migration of vesicles toward the center of the Poiseuille flow. This is in a marked contrast with a result [L. G. Leal, Annu. Rev. Fluid Mech. 12, 435 (1980)] according to which the droplet moves away from the center (provided there is no viscosity contrast between the internal and the external fluids). The migration velocity is found to increase with the local capillary number (defined by the time scale of the vesicle relaxation toward its equilibrium shape times the local shear rate), but reaches a plateau above a certain value of the capillary number. This plateau value increases with the curvature of the parabolic flow profile. We present scaling laws for the migration velocity.

Numerical simulations of complex fluid-fluid interface dynamics

Interfaces between two fluids are ubiquitous and of special importance for industrial applications, e.g., stabilisation of emulsions. The dynamics of fluid-fluid interfaces is difficult to study because these interfaces are usually deformable and their shapes are not known a priori. Since experiments do not provide access to all observables of interest, computer simulations pose attractive alternatives to gain insight into the physics of interfaces. In the present article, we restrict ourselves to systems with dimensions comparable to the lateral interface extensions. We provide a critical discussion of three numerical schemes coupled to the lattice Boltzmann method as a solver for the hydrodynamics of the problem: (a) the immersed boundary method for the simulation of vesicles and capsules, the Shan-Chen pseudopotential approach for multi-component fluids in combination with (b) an additional advection-diffusion component for surfactant modelling and (c) a molecular dynamics algorithm for the simulation of nanoparticles acting as emulsifiers.

Two-dimensional vesicle dynamics under shear flow: effect of confinement.

Dynamics of a single vesicle under shear flow between two parallel plates is studied in two-dimensions using lattice-Boltzmann simulations. We first present how we adapted the lattice-Boltzmann method to simulate vesicle dynamics, using an approach known from the immersed boundary method. The fluid flow is computed on an Eulerian regular fixed mesh while the location of the vesicle membrane is tracked by a Lagrangian moving mesh. As benchmarking tests, the known vesicle equilibrium shapes in a fluid at rest are found and the dynamical behavior of a vesicle under simple shear flow is being reproduced. Further, we focus on investigating the effect of the confinement on the dynamics, a question that has received little attention so far. In particular, we study how the vesicle steady inclination angle in the tank-treading regime depends on the degree of confinement. The influence of the confinement on the effective viscosity of the composite fluid is also analyzed. At a given reduced volume (the swelling degree) of a vesicle we find that both the inclination angle, and the membrane tank-treading velocity decrease with increasing confinement. At sufficiently large degree of confinement the tank-treading velocity exhibits a nonmonotonous dependence on the reduced volume and the effective viscosity shows a nonlinear behavior.

Interplay of inertia and deformability on rheological properties of a suspension of capsules

The interplay of inertia and deformability has a substantial impact on the transport of soft particles suspended in a fluid. However, to date a thorough understanding of these systems is still missing, and only a limited number of experimental and theoretical studies are available. We combine the finite-element, immersed-boundary and lattice-Boltzmann methods to simulate three-dimensional suspensions of soft particles subjected to planar Poiseuille flow at finite Reynolds numbers. Our findings confirm that the particle deformation and inclination increase when inertia is present. We observe that the Segre‐Silberberg effect is suppressed with respect to the particle deformability. Depending on the deformability and strength of inertial effects, inward or outward lateral migration of the particles takes place. In particular, for increasing Reynolds numbers and strongly deformable particles, a hitherto unreported distinct flow focusing effect emerges, which is accompanied by a non-monotonic behaviour of the apparent suspension viscosity and thickness of the particle-free layer close to the channel walls. This effect can be explained by the behaviour of a single particle and the change of the particle collision mechanism when both deformability and inertia effects are relevant.

How does confinement affect the dynamics of viscous vesicles and red blood cells

Despite its significance in microfluidics, the effect of confinement on the transition from the tank-treading (steady motion) to the tumbling (unsteady motion) dynamical state of deformable micro-particles has not been studied in detail. In this paper, we investigate the dynamics of a single viscous vesicle under confining shear as a general model system for red blood cells, capsules, or viscous droplets. The transition from tank-treading to tumbling motion can be triggered by the ratio between internal and external fluid viscosities. Here, we show that the transition can be induced solely by reducing the confinement, keeping the viscosity contrast constant. The observed dynamics results from the variation of the relative importance of viscous-, pressure-, and lubrication-induced torques exerted upon the vesicle. Our findings are of interest for designing future experiments or microfluidic devices: the possibility to trigger the tumbling-to-tank-treading transition either by geometry or viscosity contrast alone opens attractive possibilities for microrheological measurements as well as the detection and diagnosis of diseased red blood cells in confined flow.

Inversion of hematocrit partition at microfluidic bifurcations

Partitioning of red blood cells (RBCs) at the level of bifurcations in the microcirculatory system affects many physiological functions yet it remains poorly understood. We address this problem by using T-shaped microfluidic bifurcations as a model. Our computer simulations and in vitro experiments reveal that the hematocrit (ϕ0) partition depends strongly on RBC deformability, as long as ϕ0<20% (within the normal range in microcirculation), and can even lead to complete deprivation of RBCs in a child branch. Furthermore, we discover a deviation from the Zweifach-Fung effect which states that the child branch with lower flow rate recruits less RBCs than the higher flow rate child branch. At small enough ϕ0, we get the inverse scenario, and the hematocrit in the lower flow rate child branch is even higher than in the parent vessel. We explain this result by an intricate up-stream RBC organization and we highlight the extreme dependence of RBC transport on geometrical and cell mechanical properties. These parameters can lead to unexpected behaviors with consequences on the microcirculatory function and oxygen delivery in healthy and pathological conditions.

Forced transport of deformable containers through narrow constrictions

We study, numerically and analytically, the forced transport of deformable containers through a narrow constriction. Our central aim is to quantify the competition between the constriction geometry and the active forcing, regulating whether and at which speed a container may pass through the constriction and under what conditions it gets stuck. We focus, in particular, on the interrelation between the force that propels the container and the radius of the channel, as these are the external variables that may be directly controlled in both artificial and physiological settings. We present lattice Boltzmann simulations that elucidate in detail the various phases of translocation and present simplified analytical models that treat two limiting types of these membrane containers: deformational energy dominated by the bending or stretching contribution. In either case we find excellent agreement with the full simulations, and our results reveal that not only the radius but also the length of the constriction determines whether or not the container will pass.

Complex dynamics of a bilamellar vesicle as a simple model for leukocytes

The influence of the internal structure of a biological cell (e.g., a leukocyte) on its dynamics and rheology is not yet fully understood. By using 2D numerical simulations of a bilamellar vesicle (BLV) consisting of two vesicles as a cell model, we find that increasing the size of the inner vesicle (mimicking the nucleus) triggers a tank-treading-to-tumbling transition. A new dynamical state is observed, the undulating motion: the BLV inclination with respect to the imposed flow oscillates while the outer vesicle develops rotating lobes. The BLV exhibits a non-Newtonian behavior with a time-dependent apparent viscosity during its unsteady motion. Depending on its inclination and on its inner vesicle dynamical state, the BLV behaves like a solid or a liquid.

Interplay between microdynamics and macrorheology in vesicle suspensions

The microscopic dynamics of objects suspended in a fluid determines the macroscopic rheology of a suspension. For example, as shown by Danker and Misbah [Phys. Rev. Lett., 2007, 98, 088104], the viscosity of a dilute suspension of fluid-filled vesicles is a non-monotonic function of the viscosity contrast (the ratio between the viscosities of the internal encapsulated and the external suspending fluids) and exhibits a minimum at the critical point of the tank-treading-to-tumbling transition. By performing numerical simulations, we recover this effect and demonstrate that it persists for a wide range of vesicle parameters such as the concentration, membrane deformability, or swelling degree. We also explain why other numerical and experimental studies lead to contradicting results. Furthermore, our simulations show that this effect even persists in non-dilute and confined suspensions, but that it becomes less pronounced at higher concentrations and for more swollen vesicles. For dense suspensions and for spherical (circular in 2D) vesicles, the intrinsic viscosity tends to depend weakly on the viscosity contrast.