Characteristic spatial scale of vesicle pair interactions in a plane linear flow

Intermediate regime and a phase diagram of red blood cell dynamics in a linear flow

In this paper we investigate the in vitro dynamics of a single rabbit red blood cell (RBC) in a planar linear flow as a function of a shear stress σ and the dynamic viscosity of outer fluid η_{o}. A linear flow is a generalization of previous studies dynamics of soft objects including RBC in shear flow and is realized in the experiment in a microfluidic four-roll mill device. We verify that the RBC stable orientation dynamics is found in the experiment being the in-shear-plane orientation and the RBC dynamics is characterized by observed three RBC dynamical states, namely tumbling (TU), intermediate (INT), and swinging (SW) [or tank-treading (TT)] on a single RBC. The main results of these studies are the following. (i) We completely characterize the RBC dynamical states and reconstruct their phase diagram in the case of the RBC in-shear-plane orientation in a planar linear flow and find it in a good agreement with that obtained in early experiments in a shear flow for human RBCs. (ii) The value of the critical shear stress σ_{c} of the TU-TT(SW) transition surprisingly coincides with that found in early experiments in spite of a significant difference in the degree of RBC shape deformations in both the SW and INT states. (iii) We describe the INT regime, which is stationary, characterized by strong RBC shape deformations and observed in a wide range of the shear stresses. We argue that our observations cast doubts on the main claim of the recent numerical simulations that the only RBC spheroidal stress-free shape is capable to explain the early experimental data. Finally, we suggest that the amplitude dependence of both θ and the shape deformation parameter D on σ can be used as the quantitative criterion to determine the RBC stress-free shape.

Stagnation point flows in analytical chemistry and life sciences

Isolated microfluidic stagnation points – formed within microfluidic interfaces – have come a long way as a tool for characterizing materials, manipulating micro particles, and generating confined flows and localized chemistries.

Experimental observation of the asymmetric instability of intermediate-reduced-volume vesicles in extensional flow

Vesicles provide an attractive model system to understand the deformation of living cells in response to mechanical forces. These simple, enclosed lipid bilayer membranes are suitable for complementary theoretical, numerical, and experimental analysis. A recent study [Narsimhan, Spann, Shaqfeh, J. Fluid Mech., 2014, 750, 144] predicted that intermediate-aspect-ratio vesicles extend asymmetrically in extensional flow. Upon infinitesimal perturbation to the vesicle shape, the vesicle stretches into an asymmetric dumbbell with a cylindrical thread separating the two ends. While the symmetric stretching of high-aspect-ratio vesicles in extensional flow has been observed and characterized [Kantsler, Segre, Steinberg, Phys. Rev. Lett., 2008, 101, 048101] as well as recapitulated in numerical simulations by Narsimhan et al., experimental observation of the asymmetric stretching has not been reported. In this work, we present results from microfluidic cross-slot experiments observing this instability, along with careful characterization of the flow field, vesicle shape, and vesicle bending modulus. The onset of this shape transition depends on two non-dimensional parameters: reduced volume (a measure of vesicle asphericity) and capillary number (ratio of viscous to bending forces). We observed that every intermediate-reduced-volume vesicle that extends forms a dumbbell shape that is indeed asymmetric. For the subset of the intermediate-reduced-volume regime we could capture experimentally, we present an experimental phase diagram for asymmetric vesicle stretching that is consistent with the predictions of Narsimhan et al.

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.

Rheology of dense suspensions of elastic capsules: normal stresses, yield stress, jamming and confinement effects

We study the shearing rheology of dense suspensions of elastic capsules, taking aggregation-free red blood cells as a physiologically relevant example. Particles are non-Brownian and interact only via hydrodynamics and short-range repulsive forces. An analysis of the different stress mechanisms in the suspension shows that the viscosity is governed by the shear elasticity of the capsules, whereas the repulsive forces are subdominant. Evidence for a dynamic yield stress above a critical volume fraction is provided and related to the elastic properties of the capsules. The shear stress is found to follow a critical jamming scenario and is rather insensitive to the tumbling-to-tank-treading transition. The particle pressure and normal stress differences display some sensitivity to the dynamical state of the cells and exhibit a characteristic scaling, following the behavior of a single particle, in the tank-treading regime. The behavior of the viscosity in the fluid phase is rationalized in terms of effective medium models. Furthermore, the role of confinement effects, which increase the overall magnitude and enhance the shear-thinning of the viscosity, is discussed.

Fluid vesicles in flow

We review the dynamical behavior of giant fluid vesicles in various types of external hydrodynamic flow. The interplay between stresses arising from membrane elasticity, hydrodynamic flows, and the ever present thermal fluctuations leads to a rich phenomenology. In linear flows with both rotational and elongational components, the properties of the tank-treading and tumbling motions are now well described by theoretical and numerical models. At the transition between these two regimes, strong shape deformations and amplification of thermal fluctuations generate a new regime called trembling. In this regime, the vesicle orientation oscillates quasi-periodically around the flow direction while asymmetric deformations occur. For strong enough flows, small-wavelength deformations like wrinkles are observed, similar to what happens in a suddenly reversed elongational flow. In steady elongational flow, vesicles with large excess areas deform into dumbbells at large flow rates and pearling occurs for even stronger flows. In capillary flows with parabolic flow profile, single vesicles migrate towards the center of the channel, where they adopt symmetric shapes, for two reasons. First, walls exert a hydrodynamic lift force which pushes them away. Second, shear stresses are minimal at the tip of the flow. However, symmetry is broken for vesicles with large excess areas, which flow off-center and deform asymmetrically. In suspensions, hydrodynamic interactions between vesicles add up to these two effects, making it challenging to deduce rheological properties from the dynamics of individual vesicles. Further investigations of vesicles and similar objects and their suspensions in steady or time-dependent flow will shed light on phenomena such as blood flow.

Fully automated digital holographic processing for monitoring the dynamics of a vesicle suspension under shear flow

We investigate the dynamics of a vesicle suspension under shear flow between plates using DHM with a spatially reduced coherent source. Holograms are grabbed at a frequency of 24 frames/sec. The distribution of the vesicle suspension is obtained after numerical processing of the digital holograms sequence resulting in a 4D distribution. Obtaining this distribution is not straightforward and requires special processing to automate the analysis. We present an original method that fully automates the analysis and provides distributions that are further analyzed to extract physical properties of the fluid. Details of the numerical implementation, as well as sample experimental results are presented.

Complex dynamics of compound vesicles in linear flow

We report first experimental observations of dynamics of compound vesicles in linear flow realized in a microfluidic four-roll mill. We show that while a compound vesicle undergoes the same main tank-treading, trembling (TR), and tumbling regimes, its dynamics are far richer and more complex than that of unilamellar vesicles. A new swinging motion of the inner vesicle is found in TR in accord with simulations. The inner and outer vesicles can exist simultaneously in different dynamical regimes and can undergo either synchronized or unsynchronized motions depending on the filling factor. A compound vesicle can be used as a physical model to study white blood cell dynamics in flow similar to a unilamellar vesicle used successfully to model anucleate cells.

Rheology of a vesicle suspension with finite concentration: a numerical study

Vesicles, closed membranes made of a bilayer of phospholipids, are considered as a biomimetic system for the mechanics of red blood cells. The understanding of their dynamics under flow and their rheology is expected to help the understanding of the behavior of blood flow. We conduct numerical simulations of a suspension of vesicles in two dimensions at a finite concentration in a shear flow imposed by countertranslating rigid bounding walls by using an appropriate Green's function. We study the dynamics of vesicles, their spatial configurations, and their rheology, namely, the effective viscosity η(eff). A key parameter is the viscosity contrast λ (the ratio between the viscosity of the encapsulated fluid over that of the suspending fluid). For small enough λ, vesicles are known to exhibit tank treading (TT), while at higher λ they exhibit tumbling (TB). We find that η(eff) decreases in the TT regime, passes a minimum at a critical λ=λ(c), and increases in the TB regime. This result confirms previous theoretical and numerical works performed in the extremely dilute regime, pointing to the robustness of the picture even in the presence of hydrodynamic interactions. Our results agree also with very recent numerical simulations performed in three dimensions both in the dilute and more concentrated regime. This points to the fact that dimensionality does not alter the qualitative features of η(eff). However, they disagree with recent simulations in two dimensions. We provide arguments about the possible sources of this disagreement.

Continuum- and particle-based modeling of shapes and dynamics of red blood cells in health and disease

We review recent advances in multiscale modeling of the mechanics of healthy and diseased red blood cells (RBCs), and blood flow in the microcirculation. We cover the traditional continuum-based methods but also particle-based methods used to model both the RBCs and the blood plasma. We highlight examples of successful simulations of blood flow including malaria and sickle cell anemia.