Dynamics and rheology of a dilute suspension of vesicles: higher-order theory

A theoretical study on the dynamics of a compound vesicle in shear flow

The dynamics of nucleate cells in shear flow is of great relevance in cancer cells and circulatory tumor cells where they determine the flow properties of blood. Buoyed by the success of giant unilamellar vesicles in explaining the dynamics of anucleate cells such as red blood cells, compound vesicles have been suggested as a simple model for nucleate cells. A compound vesicle consists of two concentric unilamellar vesicles with the inner, annular and outer regions filled with aqueous Newtonian solvents. In this work, a theoretical model is presented to study the deformation and dynamics of a compound vesicle in linear shear flow using small deformation theory and spherical harmonics with higher order approximation to the membrane forces. A coupling of viscous and membrane stresses at the membrane interface of the two vesicles results in highly nonlinear shape evolution equations for the inner and the outer vesicles which are solved numerically. The results indicate that the size of the inner vesicle (χ) does not affect the tank-treading dynamics of the outer vesicle. The inner vesicle admits a greater inclination angle than the outer vesicle. However, the transition to trembling/swinging and tumbling is significantly affected. The inner and outer vesicles exhibit identical dynamics in the parameter space defined by the nondimensional rotational (Λ_an) and extensional (S) strength of the general shear flow. At moderate χ, a swinging mode is observed for the inner vesicle while the outer vesicle exhibits tumbling. The inner vesicle also exhibits modification of the TU mode to IUS (intermediate tumbling swinging) mode. Moreover, synchronization of the two vesicles at higher χ and a Capillary number sensitive motion at lower χ is observed in the tumbling regime. These results are in accordance with the few experimental observations reported by Levant and Steinberg. A reduction in the inclination angle is observed with an increase in χ when the inner vesicle is replaced by a solid inclusion. Additionally, a very elaborate phase diagram is presented in the Λ_an-S parameter space, which could be tested in future experiments or numerical simulations.

Breakage of vesicles in a simple shear flow

The breakage of micron-size SOPC lipid vesicles in an aqueous suspension was studied in a simple shear flow over a range of shear rates (1000 s-1 to 4000 s-1). Evolution of the vesicle size distribution with time was determined using optical microscopy. The number average vesicle diameter was found to reduce continuously with time; at the highest shear rate (4000 s-1), the reduction was 38% after 6 h of shearing. The distributions indicated the existence of a critical diameter such that the number of vesicles larger than the critical diameter decreased and vesicles smaller than the critical diameter increased. The capillary number for the system (ratio of the characteristic viscous stress to the characteristic stress for stretching the lipid membrane) was two orders of magnitude lower than the values reported for breakage in an ultrasonically generated flow, indicating that vesicles do not rupture due to lysis of the membrane. Direct observation of the process in a shear cell fitted in a microscope stage revealed the mechanism of rupture. Measurements of the vesicle dimensions indicated an increase in aspect ratio with time as a result of leakage of fluid from inside vesicles. Once the aspect ratio increased above a threshold value, the vesicles elongated into long thread-like shapes, which broke into small daughter vesicles by pearling.

Effect of ac electric field on the dynamics of a vesicle under shear flow in the small deformation regime

Vesicles or biological cells under simultaneous shear and electric field can be encountered in dielectrophoretic devices or designs used for continuous flow electrofusion or electroporation. In this work, the dynamics of a vesicle subjected to simultaneous shear and uniform alternating current (ac) electric field is investigated in the small deformation limit. The coupled equations for vesicle orientation and shape evolution are derived theoretically, and the resulting nonlinear equations are handled numerically to generate relevant phase diagrams that demonstrate the effect of electrical parameters on the different dynamical regimes such as tank treading (TT), vacillating breathing (VB) [called trembling (TR) in this work], and tumbling (TU). It is found that while the electric Mason number (Mn), which represents the relative strength of the electrical forces to the shear forces, promotes the TT regime, the response itself is found to be sensitive to the applied frequency as well as the conductivity ratio. While higher outer conductivity promotes orientation along the flow axis, orientation along the electric field is favored when the inner conductivity is higher. Similarly a switch of orientation from the direction of the electric field to the direction of flow is possible by a mere change of frequency when the outer conductivity is higher. Interestingly, in some cases, a coupling between electric field-induced deformation and shear can result in the system admitting an intermediate TU regime while attaining the TT regime at high Mn. The results could enable designing better dielectrophoretic devices wherein the residence time as well as the dynamical states of the vesicular suspension can be controlled as per the application.

Viscoelastic transient of confined red blood cells

The unique ability of a red blood cell to flow through extremely small microcapillaries depends on the viscoelastic properties of its membrane. Here, we study in vitro the response time upon flow startup exhibited by red blood cells confined into microchannels. We show that the characteristic transient time depends on the imposed flow strength, and that such a dependence gives access to both the effective viscosity and the elastic modulus controlling the temporal response of red cells. A simple theoretical analysis of our experimental data, validated by numerical simulations, further allows us to compute an estimate for the two-dimensional membrane viscosity of red blood cells, η(mem)(2D) ∼ 10(-7) N ⋅ s ⋅ m(-1). By comparing our results with those from previous studies, we discuss and clarify the origin of the discrepancies found in the literature regarding the determination of η(mem)(2D), and reconcile seemingly conflicting conclusions from previous works.

Diffuse interface models of locally inextensible vesicles in a viscous fluid

We present a new diffuse interface model for the dynamics of inextensible vesicles in a viscous fluid with inertial forces. A new feature of this work is the implementation of the local inextensibility condition in the diffuse interface context. Local inextensibility is enforced by using a local Lagrange multiplier, which provides the necessary tension force at the interface. We introduce a new equation for the local Lagrange multiplier whose solution essentially provides a harmonic extension of the multiplier off the interface while maintaining the local inextensibility constraint near the interface. We also develop a local relaxation scheme that dynamically corrects local stretching/compression errors thereby preventing their accumulation. Asymptotic analysis is presented that shows that our new system converges to a relaxed version of the inextensible sharp interface model. This is also verified numerically. To solve the equations, we use an adaptive finite element method with implicit coupling between the Navier-Stokes and the diffuse interface inextensibility equations. Numerical simulations of a single vesicle in a shear flow at different Reynolds numbers demonstrate that errors in enforcing local inextensibility may accumulate and lead to large differences in the dynamics in the tumbling regime and smaller differences in the inclination angle of vesicles in the tank-treading regime. The local relaxation algorithm is shown to prevent the accumulation of stretching and compression errors very effectively. Simulations of two vesicles in an extensional flow show that local inextensibility plays an important role when vesicles are in close proximity by inhibiting fluid drainage in the near contact region.

Rheological properties of a vesicle suspension

The rheological behavior of a dilute suspension of vesicles in linear shear flow at a finite concentration is analytically examined. In the quasispherical limit, two coupled nonlinear equations that describe the vesicle orientation in the flow and its shape evolution were derived [Phys. Rev. Lett. 96, 028104 (2006)PRLTAO0031-900710.1103/PhysRevLett.96.028104] and serve here as a starting point. Of special interest is to provide, for the first time, an exact analytical prediction of the time-dependent effective viscosity η_{eff} and normal stress differences N_{1} and N_{2}. Our results shed light on the effect of the viscosity ratio λ (defined as the inner over the outer fluid viscosities) as the main controlling parameter. It is shown that η_{eff},N_{1}, and N_{2} either tend to a steady state or describe a periodic time-dependent rheological response, previously reported numerically and experimentally. In particular, the shear viscosity minimum and the cusp singularities of η_{eff},N_{1}, and N_{2} at the tumbling threshold are brought to light. We also report on rheology properties for an arbitrary linear flow. We were able to obtain a constitutive law in a closed form relating the stress tensor to the strain rate tensor. It is found that the resulting constitutive markedly contrasts with classical laws known for other complex fluids, such as emulsions, capsule suspensions, and dilute polymer solutions (Oldroyd B model). We highlight the main differences between our law and classical laws.

Prediction of anomalous blood viscosity in confined shear flow

Red blood cells play a major role in body metabolism by supplying oxygen from the microvasculature to different organs and tissues. Understanding blood flow properties in microcirculation is an essential step towards elucidating fundamental and practical issues. Numerical simulations of a blood model under a confined linear shear flow reveal that confinement markedly modifies the properties of blood flow. A nontrivial spatiotemporal organization of blood elements is shown to trigger hitherto unrevealed flow properties regarding the viscosity η, namely ample oscillations of its normalized value [η] = (η-η(0))/(η(0)ϕ) as a function of hematocrit ϕ (η(0) = solvent viscosity). A scaling law for the viscosity as a function of hematocrit and confinement is proposed. This finding can contribute to the conception of new strategies to efficiently detect blood disorders, via in vitro diagnosis based on confined blood rheology. It also constitutes a contribution for a fundamental understanding of rheology of confined complex fluids.

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.

Microcapsule mechanics: from stability to function

Microcapsules are reviewed with special emphasis on the relevance of controlled mechanical properties for functional aspects. At first, assembly strategies are presented that allow control over the decisive geometrical parameters, diameter and wall thickness, which both influence the capsule's mechanical performance. As one of the most powerful approaches the layer-by-layer technique is identified. Subsequently, ensemble and, in particular, single-capsule deformation techniques are discussed. The latter generally provide more in-depth information and cover the complete range of applicable forces from smaller than pN to N. In a theory chapter, we illustrate the physics of capsule deformation. The main focus is on thin shell theory, which provides a useful approximation for many deformation scenarios. Finally, we give an overview of applications and future perspectives where the specific design of mechanical properties turns microcapsules into (multi-)functional devices, enriching especially life sciences and material sciences.

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.

Amoeboid swimming: a generic self-propulsion of cells in fluids by means of membrane deformations

Microorganisms, such as bacteria, algae, or spermatozoa, are able to propel themselves forward thanks to flagella or cilia activity. By contrast, other organisms employ pronounced changes of the membrane shape to achieve propulsion, a prototypical example being the Eutreptiella gymnastica. Cells of the immune system as well as dictyostelium amoebas, traditionally believed to crawl on a substratum, can also swim in a similar way. We develop a model for these organisms: the swimmer is mimicked by a closed incompressible membrane with force density distribution (with zero total force and torque). It is shown that fast propulsion can be achieved with adequate shape adaptations. This swimming is found to consist of an entangled pusher-puller state. The autopropulsion distance over one cycle is a universal linear function of a simple geometrical dimensionless quantity A/V(2/3) (V and A are the cell volume and its membrane area). This study captures the peculiar motion of Eutreptiella gymnastica with simple force distribution.

Noisy nonlinear dynamics of vesicles in flow

We present a model for the dynamics of fluid vesicles in linear flow which consistently includes thermal fluctuations and nonlinear coupling between different modes. At the transition between tank treading and tumbling, we predict a trembling motion which is at odds with the known deterministic motions and for which thermal noise is strongly amplified. In particular, highly asymmetric shapes are observed even though the deterministic flow only allows for axisymmetric ones. Our results explain quantitatively recent experimental observations [Levant and Steinberg, Phys. Rev. Lett. 109, 268103 (2012)].

Amplification of thermal noise by vesicle dynamics

A novel noise amplification mechanism resulting from the interaction of thermal fluctuations and nonlinear vesicle dynamics is reported. It is observed in a time-dependent vesicle state called trembling (TR). High spatial resolution and very long time series of TR compared to the vesicle period allow us to quantitatively analyze the generation and amplification of spatial and temporal modes of the vesicle shape perturbations. During a compression part of each TR cycle, a vesicle finds itself on the edge of the wrinkling instability, where thermally excited spatial modes are amplified.

Squaring, parity breaking, and S tumbling of vesicles under shear flow

The numerical study of 3D vesicles with a reduced volume equal to that of human red blood cells leads to the discovery of three types of dynamics: (i) squaring motion, in which the angle between the direction of the longest distance and the flow velocity undergoes discontinuous jumps over time, (ii) spontaneous parity breaking of the shape leading to cross-streamline migration, and (iii) S tumbling where the vesicle tumbles, exhibiting a pronounced S-like shape with a waisted morphology in the center. We report on the phase diagram within a wide range of relevant parameters. Our estimates reveal that healthy and pathological red blood cells are also prone to these types of motion, which may affect blood microcirculation and impact oxygen transport.

Dynamic modes of quasispherical vesicles: exact analytical solutions

In this paper we introduce a simple mathematical analysis to reexamine vesicle dynamics in the quasispherical limit (small deformation) under a shear flow. In this context, a recent paper [Misbah, Phys. Rev. Lett. 96, 028104 (2006)] revealed a dynamic referred to as the vacillating-breathing (VB) mode where the vesicle main axis oscillates about the flow direction and the shape undergoes a breathinglike motion, as well as the tank-treading and tumbling (TB) regimes. Our goal here is to identify these three modes by obtaining explicit analytical expressions of the vesicle inclination angle and the shape deformation. In particular, the VB regime is put in evidence and the transition dynamics is discussed. Not surprisingly, our finding confirms the Keller-Skalak solutions (for rigid particles) and shows that the VB and TB modes coexist, and whether one prevails over the other depends on the initial conditions. An interesting additional element in the discussion is the prediction of the TB and VB modes as functions of a control parameter Γ, which can be identified as a TB-VB parameter.

Hydrodynamic interaction between two vesicles in a linear shear flow: asymptotic study

Interactions between two vesicles in an imposed linear shear flow are studied theoretically, in the limit of almost spherical vesicles, with a large intervesicle distance, in a strong flow, with a large inner to outer viscosity ratio. This allows to derive a system of ordinary equations describing the dynamics of the two vesicles. We provide an analytic expression for the interaction law. We find that when the vesicles are in the same shear plane, the hydrodynamic interaction leads to a repulsion. When they are not, the interaction may turn into attraction instead. The interaction law is discussed and analyzed as a function of relevant parameters.

Vesicle dynamics under weak flows: application to large excess area

Dynamics of a vesicle under simple shear flow is studied in the limit of small capillary number. A perturbative approach is used to derive the equation of vesicle dynamics. The expansions are shown to converge for significantly deflated vesicles (with excess area from the sphere as high as 2). In particular, we provide an explicit analytical expression for the tank-treading to tumbling bifurcation point. This expression is valid for excess areas up to 2.5. The results are compared with full 3D numerical simulations. The proposed method can be used for analytical or numerical solution of vesicle dynamics under weak flow of general form.

Three-dimensional numerical simulation of vesicle dynamics using a front-tracking method

Three-dimensional numerical simulation using the front-tracking method is presented on the dynamics of a vesicle in a linear shear flow. The focus here is to elucidate the parametric dependence and the self-similarity of the vesicle dynamics, quantification of vesicle deformation, and the analysis of shape dynamics. A detailed comparison of the numerical results is made with various theoretical models and experiments. It is found that the applicability of the theoretical models is limited despite some general agreement with the simulations and experiments. The deviations between the perturbative results and the simulation results occur even in the absence of thermal noise. Specifically, we find that the vesicle dynamics does not follow a self-similar behavior in a two-parameter phase space, as proposed in a theoretical model. Rather, the dynamics is governed by three controlling parameters, namely, the excess area, viscosity ratio, and dimensionless shear rate. Additionally, we find that a linear scaling of the tank-treading angle, as proposed in the theoretical model, is possible only for nearly spherical vesicles. The breakdown of the scaling occurs at higher values of the excess area even in the absence of thermal noise. We find that the vesicle deformation saturates at large shear rates, and the asymptotic deformation matches well with a theoretical prediction for nearly spherical vesicles. The dependence of the critical viscosity ratio associated with the onset of unsteady dynamics on the vesicle excess area is in excellent agreement with the experimental observation. We show that near the transition between the tank-treading and tumbling dynamics, both the vacillating-breathing-like motion characterized by a smooth ellipsoidal shape and the trembling-like motion characterized by a highly deformed shape are possible. For the trembling-like motion, the shape is highly three-dimensional with concavities and lobes, and the vesicle deforms more in the vorticity direction than in the shear plane. A Fourier spectral analysis of the vesicle shape shows the presence of the odd harmonics and higher order modes beyond fourth order.

Orientation and internal flow of a vesicle in tank-treading motion in shear flow

Deformation, orientation and internal flow of lipid bilayer vesicles in linear shear flows are investigated using phase contrast microscopy. We construct a rotating-cylinder apparatus, which can generate a linear shear flow with constant shear rates. Vesicles are prepared from 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) by the gentle hydration method. When visualizing internal flows, polystyrene tracer particles are mixed with the hydration water solution. In our observation, vesicles deform to steady ellipsoidal shapes and show constant orientations given by θ(i), which is the angle between the major axis and the flow direction. The tracer particles inside a vesicle rotate around the center of the vesicle along ellipsoidal orbits, which are homothetic to the shape of the vesicle. It is shown that the relationship between θ(i) and the swelling ratio (volume/surface ratio) S(w) agrees quantitatively with the experimental result of Abkarian et al. [Biophys. J.89, 1055 (2005)], which was obtained with vesicles in wall-bounded shear flows. It also agrees with a theoretical analysis of Keller and Skalak [J. Fluid Mech.120, 27 (1982)] and other numerical simulations. It is also shown that angular velocities of the particles near the membrane change periodically and agree quantitatively with the experimental result for the motion of a particle adhering to the membrane of a tank-treading vesicle [Kantsler and Steinberg, Phys. Rev. Lett.95, 258101 (2005)]. A statistical analysis indicates that the velocity of the internal fluid close to the membrane is not constant along the circumference, which implies the possibility of a three-dimensional flow field of the lipid molecules or an apparent stretching motion of the membrane by the effect of hidden surface area due to thermal fluctuation.

Symmetry breaking of vesicle shapes in Poiseuille flow

Vesicle behavior under unbounded axial Poiseuille flow is studied analytically. Our study reveals subtle features of the dynamics. It is established that there exists a stable off-centerline steady-state solution for low enough flow strength. This solution appears as a symmetry-breaking bifurcation upon lowering the flow strength and includes slipper shapes, which are characteristic of red blood cells in the microvasculature. A stable axisymmetric solution exists for any flow strength provided the excess area is small enough. It is shown that the mechanism of the symmetry breaking depends on the geometry of the flow: The bifurcation is subcritical in axial Poiseuille flow and supercritical in planar flow.

Tank-treading and tumbling frequencies of capsules and red blood cells

This study is motivated in part by the discrepancy that exists in the literature with regard to the dependence of the tank-treading frequency of red blood cells on the shear rate and suspending medium viscosity. Here we consider three-dimensional numerical simulations of deformable capsules of initially spherical and oblate spheroidal shapes and biconcave discoid representing the red blood cell resting shape. By considering a much broader range of the viscosity ratio (ratio of capsule or cell interior to suspending fluid viscosity), shear rate, and aspect ratio (ratio of minor to major axes) than that considered in the previous experiments, we find several new characteristics of the tank-treading and tumbling frequencies that have not been reported earlier. These new characteristics are the result of the large shape deformation and the coupling between shape and angular oscillations of the capsules or cells. For the spherical and oblate spheroidal capsules, the tank-treading frequency shows a nonmonotonic trend that is characterized by an initial decrease leading to a minimum followed by an increase with increasing viscosity ratio. For red blood cells, we find two regimes of the viscosity dependence of the tank-treading frequency: an exponential regime in which the tank-treading frequency decreases at a slower rate with increasing viscosity ratio, and a logarithmic range in which it decreases at a much faster rate. While this trend agrees well with different theoretical models of shape-preserving capsules, it was not evident in previous experimental results. When the shear rate dependence is considered, the tank-treading frequency of red blood cells and capsules of highly elongated initial shapes exhibits a nonmonotonic trend that is characterized by an initial increase leading to a maximum followed by a sharp decrease with decreasing shear rate. This anomalous behavior of the tank-treading frequency is shown to be due to a breathing-like dynamics of the capsule or cell that is characterized by a repeated emergence and absence of deep, crater-like dimples, and a large swinging motion. We further observe that the tumbling frequency exhibits a decreasing trend with increasing viscosity ratio that is in contrast to the theoretical result for the shape-preserving capsules and is due to the periodic deformation and preferential alignment of the capsules in the extensional quadrant of the flow.

Vesicle electrohydrodynamics

A small amplitude perturbation analysis is developed to describe the effect of a uniform electric field on the dynamics of a lipid bilayer vesicle in a simple shear flow. All media are treated as leaky dielectrics and fluid motion is described by the Stokes equations. The instantaneous vesicle shape is obtained by balancing electric, hydrodynamic, bending, and tension stresses exerted on the membrane. We find that in the absence of ambient shear flow, it is possible that an applied stepwise uniform dc electric field could cause the vesicle shape to evolve from oblate to prolate over time if the encapsulated fluid is less conducting than the suspending fluid. For a vesicle in ambient shear flow, the electric field damps the tumbling motion, leading to a stable tank-treading state.

Three-dimensional vesicles under shear flow: numerical study of dynamics and phase diagram

The study of vesicles under flow, a model system for red blood cells (RBCs), is an essential step in understanding various intricate dynamics exhibited by RBCs in vivo and in vitro. Quantitative three-dimensional analyses of vesicles under flow are presented. The regions of parameters to produce tumbling (TB), tank-treating, vacillating-breathing (VB), and even kayaking (or spinning) modes are determined. New qualitative features are found: (i) a significant widening of the VB mode region in parameter space upon increasing shear rate γ and (ii) a robustness of normalized period of TB and VB with γ. Analytical support is also provided. We make a comparison with existing experimental results. In particular, we find that the phase diagram of the various dynamics depends on three dimensionless control parameters, while a recent experimental work reported that only two are sufficient.

Analytical progress in the theory of vesicles under linear flow

Vesicles are becoming a quite popular model for the study of red blood cells. This is a free boundary problem which is rather difficult to handle theoretically. Quantitative computational approaches constitute also a challenge. In addition, with numerical studies, it is not easy to scan within a reasonable time the whole parameter space. Therefore, having quantitative analytical results is an essential advance that provides deeper understanding of observed features and can be used to accompany and possibly guide further numerical development. In this paper, shape evolution equations for a vesicle in a shear flow are derived analytically with precision being cubic (which is quadratic in previous theories) with regard to the deformation of the vesicle relative to a spherical shape. The phase diagram distinguishing regions of parameters where different types of motion (tank treading, tumbling, and vacillating breathing) are manifested is presented. This theory reveals unsuspected features: including higher order terms and harmonics (even if they are not directly excited by the shear flow) is necessary, whatever the shape is close to a sphere. Not only does this theory cure a quite large quantitative discrepancy between previous theories and recent experiments and numerical studies, but also it reveals a phenomenon: the VB mode band in parameter space, which is believed to saturate after a moderate shear rate, exhibits a striking widening beyond a critical shear rate. The widening results from excitation of fourth-order harmonic. The obtained phase diagram is in a remarkably good agreement with recent three-dimensional numerical simulations based on the boundary integral formulation. Comparison of our results with experiments is systematically made.

Dynamic modes of red blood cells in oscillatory shear flow

The dynamics of red blood cells (RBCs) in oscillatory shear flow was studied using differential equations of three variables: a shape parameter, the inclination angle θ, and phase angle ϕ of the membrane rotation. In steady shear flow, three types of dynamics occur depending on the shear rate and viscosity ratio. (i) tank-treading (TT): ϕ rotates while the shape and θ oscillate. (ii) tumbling (TB): θ rotates while the shape and ϕ oscillate. (iii) intermediate motion: both ϕ and θ rotate synchronously or intermittently. In oscillatory shear flow, RBCs show various dynamics based on these three motions. For a low shear frequency with zero mean shear rate, a limit-cycle oscillation occurs, based on the TT or TB rotation at a high or low shear amplitude, respectively. This TT-based oscillation well explains recent experiments. In the middle shear amplitude, RBCs show an intermittent or synchronized oscillation. As shear frequency increases, the vesicle oscillation becomes delayed with respect to the shear oscillation. At a high frequency, multiple limit-cycle oscillations coexist. The thermal fluctuations can induce transitions between two orbits at very low shear amplitudes. For a high mean shear rate with small shear oscillation, the shape and θ oscillate in the TT motion but only one attractor exists even at high shear frequencies. The measurement of these oscillatory modes is a promising tool for quantifying the viscoelasticity of RBCs, synthetic capsules, and lipid vesicles.

Dynamic modes of microcapsules in steady shear flow: effects of bending and shear elasticities

The dynamics of microcapsules in steady shear flow were studied using a theoretical approach based on three variables: the Taylor deformation parameter αD , the inclination angle θ , and the phase angle ϕ of the membrane rotation. It is found that the dynamic phase diagram shows a remarkable change with an increase in the ratio of the membrane shear and bending elasticities. A fluid vesicle (no shear elasticity) exhibits three dynamic modes: (i) tank treading at low viscosity ηin of internal fluid (αD and θ relaxes to constant values), (ii) tumbling (TB) at high ηin (θ rotates), and (iii) swinging (SW) at middle ηin and high shear rates γ (θ oscillates). All of three modes are accompanied by a membrane (ϕ) rotation. For microcapsules with low shear elasticity, the TB phase with no ϕ rotation and the coexistence phase of SW and TB motions are induced by the energy barrier of ϕ rotation. Synchronization of ϕ rotation with TB rotation or SW oscillation occurs with integer ratios of rotational frequencies. At high shear elasticity, where a saddle point in the energy potential disappears, intermediate phases vanish and either ϕ or θ rotation occurs. This phase behavior agrees with recent simulation results of microcapsules with low bending elasticity.

Rheology of a dilute suspension of liquid-filled elastic capsules

Rheology of a dilute suspension of liquid-filled elastic capsules in linear shear flow is studied by three-dimensional numerical simulations using a front-tracking method. This study is motivated by a recent discovery that a suspension of viscous vesicles exhibits a shear viscosity minimum when the vesicles undergo an unsteady vacillating-breathing dynamics at the threshold of a transition between the tank-treading and tumbling motions. Here we consider capsules of spherical resting shape for which only a steady tank-treading motion is observed. A comprehensive analysis of the suspension rheology is presented over a broad range of viscosity ratio (ratio of internal-to-external fluid viscosity), shear rate (or, capillary number), and capsule surface-area dilatation. We find a result that the capsule suspension exhibits a shear viscosity minimum at moderate values of the viscosity ratio, and high capillary numbers, even when the capsules are in a steady tank-treading motion. It is further observed that the shear viscosity minimum exists for capsules with area-dilating membranes but not for those with nearly incompressible membranes. Nontrivial results are also observed for the normal stress differences which are shown to decrease with increasing capillary number at high viscosity ratios. Such nontrivial results neither can be predicted by the small-deformation theory nor can be explained by the capsule geometry alone. Physical mechanisms underlying these results are studied by decomposing the particle stress tensor into a contribution due to the elastic stresses in the capsule membrane and a contribution due to the viscosity differences between the internal and suspending fluids. It is shown that the elastic contribution is shear-thinning, but the viscous contribution is shear thickening. The coupling between the capsule geometry and the elastic and viscous contributions is analyzed to explain the observed trends in the bulk rheology.

Vesicles under simple shear flow: elucidating the role of relevant control parameters

The dynamics of vesicles under shear flow are carefully analyzed in the regime of a small vesicle excess area relative to a sphere. This regime corresponds to the quasispherical limit, for which several groups have analytically extracted simple nonlinear differential equations. Under shear flow, vesicles are known to exhibit three types of motion: (i) tank-treading (TT): the vesicle assumes a steady inclination angle with respect to the flow direction, while its membrane undergoes a tank-treading motion, (ii) tumbling (TB), and (iii) vacillating-breathing (VB): the vesicle main axis oscillates about the flow direction, whereas the overall shape undergoes a breathinglike motion. The region of existence for each regime depends on material and control parameters. The whole set of parameters can be cast into three dimensionless control parameters: (i) the viscosity ratio between the internal and external fluid, lambda , (ii) the excess area relative to a sphere (this parameter measures the degree of the vesicle deflation), Delta , and (iii) the capillary number (the ratio between the vesicle relaxation time toward its equilibrium shape after cessation of the flow and the flow time scale, which is the inverse shear rate), Ca. Recent studies [Danker, Phys. Rev. E 76, 041905 (2007)] have focused on the shape of the phase diagram (representing the TT, TB, and VB regimes in the Ca-lambda plane). In this paper, the physical quantities are analyzed in detail and attention is brought to features that are essential for future experimental studies. It is shown that the boundaries delimiting different dynamical regimes (TT, TB, and VB) in parameter space depend on the three dimensionless control parameters, in contrast with a recent study [V. V. Lebedev, Phys. Rev. Lett. 99, 218101 (2007)] where it is claimed that only two parameters are relevant. Consideration of the amplitude of oscillation (of the vesicle orientation angle and its shape deformation) in the VB mode reveals an even more significant dependence on the three parameters. It is also shown that the inclination angle in the TT regime significantly depends on the shear rate (Ca), which runs contrary to common belief. Finally, we show that the TB and VB periods are quite insensitive to Ca, in marked contrast with a recent study [H. Noguchi and G. Gompper, Phys. Rev. Lett. 98, 128103 (2007)].

Dynamics of nonspherical capsules in shear flow

Three-dimensional numerical simulations using a front-tracking method are presented on the dynamics of oblate shape capsules in linear shear flow by considering a broad range of viscosity contrast (ratio of internal-to-external fluid viscosity), shear rate (or capillary number), and aspect ratio. We focus specifically on the coupling between the shape deformation and orientation dynamics of capsules, and show how this coupling influences the transition from the tank-treading to tumbling motion. At low capillary numbers, three distinct modes of motion are identified: a swinging or oscillatory (OS) mode at a low viscosity contrast in which the inclination angle theta(t) oscillates but always remains positive; a vacillating-breathing (VB) mode at a moderate viscosity contrast in which theta(t) periodically becomes positive and negative, but a full tumbling does not occur; and a pure tumbling mode (TU) at a higher viscosity contrast. At higher capillary numbers, three types of transient motions occur, in addition to the OS and TU modes, during which the capsule switches from one mode to the other as (i) VB to OS, (ii) TU to VB to OS, and (iii) TU to VB. Phase diagrams showing various regimes of capsule dynamics are presented. For all modes of motion (OS, VB, and TU), a large-amplitude oscillation in capsule shape and a strong coupling between the shape deformation and orientation dynamics are observed. It is shown that the coupling between the shape deformation and orientation is the strongest in the VB mode, and hence at a moderate viscosity contrast, for which the amplitude of shape deformation reaches its maximum. The numerical results are compared with the theories of Keller and Skalak, and Skotheim and Secomb. Significant departures from the two theories are discussed and related to the strong coupling between the shape deformation, inclination, and transition dynamics.

Dynamical regimes and hydrodynamic lift of viscous vesicles under shear

The dynamics of two-dimensional viscous vesicles in shear flow, with different fluid viscosities etain and etaout inside and outside, respectively, is studied using mesoscale simulation techniques. Besides the well-known tank-treading and tumbling motions, an oscillatory swinging motion is observed in the simulations for large shear rate. The existence of this swinging motion requires the excitation of higher-order undulation modes (beyond elliptical deformations) in two dimensions. Keller-Skalak theory is extended to deformable two-dimensional vesicles, such that a dynamical phase diagram can be predicted for the reduced shear rate and the viscosity contrast etain/etaout. The simulation results are found to be in good agreement with the theoretical predictions, when thermal fluctuations are incorporated in the theory. Moreover, the hydrodynamic lift force, acting on vesicles under shear close to a wall, is determined from simulations for various viscosity contrasts. For comparison, the lift force is calculated numerically in the absence of thermal fluctuations using the boundary-integral method for equal inside and outside viscosities. Both methods show that the dependence of the lift force on the distance ycm of the vesicle center of mass from the wall is well described by an effective power law ycm(-2) for intermediate distances 0.8Rp< approximately ycm< approximately 3Rp with vesicle radius Rp. The boundary-integral calculation indicates that the lift force decays asymptotically as 1/[ycm ln(ycm)] far from the wall.

Dynamics of a vesicle in general flow

An approach to quantitatively study vesicle dynamics as well as biologically-related micro-objects in a fluid flow, which is based on the combination of a dynamical trap and a control parameter, the ratio of the vorticity to the strain rate, is suggested. The flow is continuously varied between rotational, shearing, and elongational in a microfluidic 4-roll mill device, the dynamical trap, that allows scanning of the entire phase diagram of motions, i.e., tank-treading (TT), tumbling (TU), and trembling (TR), using a single vesicle even at lambda = eta(in)/eta(out) = 1, where eta(in) and eta(out) are the viscosities of the inner and outer fluids. This cannot be achieved in pure shear flow, where the transition between TT and either TU or TR is attained only at lambda>1. As a result, it is found that the vesicle dynamical states in a general are presented by the phase diagram in a space of only 2 dimensionless control parameters. The findings are in semiquantitative accord with the recent theory made for a quasi-spherical vesicle, although vesicles with large deviations from spherical shape were studied experimentally. The physics of TR is also uncovered.

Electrohydrodynamic model of vesicle deformation in alternating electric fields

We develop an analytical theory to explain the experimentally observed morphological transitions of quasispherical giant vesicles induced by alternating electric fields. The model treats the inner and suspending media as lossy dielectrics, and the membrane as an impermeable flexible incompressible-fluid sheet. The vesicle shape is obtained by balancing electric, hydrodynamic, bending, and tension stresses exerted on the membrane. Our approach, which is based on force balance, also allows us to describe the time evolution of the vesicle deformation, in contrast to earlier works based on energy minimization, which are able to predict only stationary shapes. Our theoretical predictions for vesicle deformation are consistent with experiment. If the inner fluid is more conducting than the suspending medium, the vesicle always adopts a prolate shape. In the opposite case, the vesicle undergoes a transition from a prolate to oblate ellipsoid at a critical frequency, which the theory identifies with the inverse membrane charging time. At frequencies higher than the inverse Maxwell-Wagner polarization time, the electrohydrodynamic stresses become too small to alter the vesicle's quasispherical rest shape. The model can be used to rationalize the transient and steady deformation of biological cells in electric fields.

Vesicles in Poiseuille flow

Blood microcirculation critically depends on the migration of red cells towards the flow centerline. We identify theoretically the ratio of the inner over the outer fluid viscosities lambda as a key parameter. At low lambda, the vesicle deforms into a tank-treading ellipsoid shape far away from the flow centerline. The migration is always towards the flow centerline, unlike drops. Above a critical lambda, the vesicle tumbles or breaths and migration is suppressed. A surprising coexistence of two types of shapes at the centerline, a bulletlike and a parachutelike shape, is predicted.

Phase diagram of single vesicle dynamical states in shear flow

We report the first experimental phase diagram of vesicle dynamical states in a shear flow presented in a space of two dimensionless parameters suggested recently by V. Lebedev et al. To reduce errors in the control parameters, 3D geometrical reconstruction and determination of the viscosity contrast of a vesicle in situ in a plane Couette flow device prior to the experiment are developed. Our results are in accord with the theory predicting three distinctly separating regions of vesicle dynamical states in the plane of just two self-similar parameters.

Micro-macro link in rheology of erythrocyte and vesicle suspensions

We report on the rheology of dilute suspensions of red blood cells (RBC) and vesicles. The viscosity of RBC suspensions reveals a previously unknown signature: it exhibits a pronounced minimum when the viscosity of the ambient medium is close to the value at which the transition from tank-treading to tumbling occurs. This bifurcation is triggered by varying the viscosity of the ambient fluid. It is found that the intrinsic viscosity of the suspension varies by about a factor of 4 in the explored parameter range. Surprisingly, this significant change of the intrinsic viscosity is revealed even at low hematocrit (5%). We suggest that this finding may be used to detect blood flow disorders linked to pathologies that affect RBC shape and mechanical properties. This opens future perspectives on setting up new diagnostic tools, with great efficiency even at very low hematocrit. Investigations are also performed on giant vesicle suspensions, and compared to RBCs.

Two-dimensional fluctuating vesicles in linear shear flow

The stochastic motion of a two-dimensional vesicle in linear shear flow is studied at finite temperature. In the limit of small deformations from a circle, Langevin-type equations of motion are derived, which are highly nonlinear due to the constraint of constant perimeter length. These equations are solved in the low-temperature limit and using a mean-field approach, in which the length constraint is satisfied only on average. The constraint imposes non-trivial correlations between the lowest deformation modes at low temperature. We also simulate a vesicle in a hydrodynamic solvent by using the multi-particle collision dynamics technique, both in the quasi-circular regime and for larger deformations, and compare the stationary deformation correlation functions and the time autocorrelation functions with theoretical predictions. Good agreement between theory and simulations is obtained.

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.