Chaouqi Misbah

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.

Rheology of a dilute two-dimensional suspension of vesicles

The rheology of a dilute two-dimensional suspension of vesicles (closed bags of a lipid bilayer membrane) is studied by numerical simulations. The numerical methods used are based on the boundary integral formulation (Greens function technique) and the phase field approach, which has become a quite popular and powerful tool for the numerical study of free-boundary problems. The imposed flow is an unbounded linear shear. The goal of the present study is to elucidate the link between the rheology of vesicle suspensions and the microscopic dynamics of the constituent particles (tank-treading and tumbling motions). A comparison with emulsion rheology reveals the central role played by the membrane. In particular, at low viscosity ratio λ (defined as the viscosity of the internal fluid over that of the ambient one), the effective viscosity decreases with λ, while the opposite trend is exhibited by emulsions, according to the classical Taylor result. This fact is explained by considering the velocity field of the ambient fluid. The area-incompressibility of the vesicle membrane modifies the surrounding velocity field in a quite different manner than what a drop does. The overall numerical results in two dimensions are in reasonable agreement with the three-dimensional analytical theory derived recently in the small deformation limit (quasi-spherical shapes). The finding that the simulations in two dimensions capture the essential features of the three-dimensional rheology opens the way for extensive and large-scale simulations for semi-dilute and concentrated vesicle suspensions. We discuss some peculiar effects exhibited by the instantaneous viscosity in the tumbling regime of vesicles. Finally, the rheology is found to be relatively insensitive to shear rate.

The plasma protein fibrinogen stabilizes clusters of red blood cells in microcapillary flows

The supply of oxygen and nutrients and the disposal of metabolic waste in the organs depend strongly on how blood, especially red blood cells, flow through the microvascular network. Macromolecular plasma proteins such as fibrinogen cause red blood cells to form large aggregates, called rouleaux, which are usually assumed to be disaggregated in the circulation due to the shear forces present in bulk flow. This leads to the assumption that rouleaux formation is only relevant in the venule network and in arterioles at low shear rates or stasis. Thanks to an excellent agreement between combined experimental and numerical approaches, we show that despite the large shear rates present in microcapillaries, the presence of either fibrinogen or the synthetic polymer dextran leads to an enhanced formation of robust clusters of red blood cells, even at haematocrits as low as 1%. Robust aggregates are shown to exist in microcapillaries even for fibrinogen concentrations within the healthy physiological range. These persistent aggregates should strongly affect cell distribution and blood perfusion in the microvasculature, with putative implications for blood disorders even within apparently asymptomatic subjects.

Dynamics of aeolian sand ripples

Abstract:We analyze theoretically the dynamics of aeolian sand ripples. In order to put the study in the context, we first review existing models. This paper is a continuation of two previous papers (Z. Csahók et al., Physica D 128, 87 (1999); A. Valance et al., Eur. Phys. J. B 10, 543 (1999)), the first one is based on symmetries and the second on a hydrodynamical model. We show how the hydrodynamical model may be modified to recover the missing terms that are dictated by symmetries. The symmetry and conservation arguments are powerful in that the form of the equation is model-independent. We then present an extensive numerical and analytical analysis of the generic sand ripple equation. We find that at the initial stage the wavelength of the ripple is that corresponding to the linearly most dangerous mode. At later stages the profile undergoes a coarsening process leading to a significant increase of the wavelength. We find that including the next higher-order nonlinear term in the equation leads naturally to a saturation of the local slope. We analyze both analytically and numerically the coarsening stage, in terms of a dynamical exponent for the mean wavelength increase. We discuss some future lines of investigations.

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.

When does coarsening occur in the dynamics of one-dimensional fronts?

Dynamics of a one-dimensional growing front with an unstable straight profile are analyzed. We argue that a coarsening process occurs if and only if the period lambda of the steady-state solution is an increasing function of its amplitude A. This statement is rigorously proved for two important classes of conserved and nonconserved models by investigating the phase diffusion equation of the steady pattern. We further provide clear numerical evidence for the growth equation of a stepped crystal surface.

Vesicle tumbling inhibited by inertia

Vesicles under flow constitute a model system for the study of red blood cells (RBCs) dynamics and blood rheology. In the blood circulatory system the Reynolds number (at the scale of the RBC) is not always small enough for the Stokes limit to be valid. We develop a numerical method in two dimensions based on the level set approach and solve the fluid/membrane coupling by using an adaptive finite element technique. We find that a Reynolds number of order one can destroy completely the vesicle tumbling motion obtained in the Stokes regime. We analyze in details this phenomenon and discuss some of the far reaching consequences. We suggest experimental tests on vesicles.

Red blood cell clustering in Poiseuille microcapillary flow

Red blood cells (RBC) flowing in microcapillaries tend to associate into clusters, i.e., small trains of cells separated from each other by a distance comparable to cell size. This process is usually attributed to slower RBCs acting to create a sequence of trailing cells. Here, based on the first systematic investigation of collective RBC flow behavior in microcapillaries in vitro by high-speed video microscopy and numerical simulations, we show that RBC size polydispersity within the physiological range does not affect cluster stability. Lower applied pressure drops and longer residence times favor larger RBC clusters. A limiting cluster length, depending on the number of cells in a cluster, is found by increasing the applied pressure drop. The insight on the mechanism of RBC clustering provided by this work can be applied to further our understanding of RBC aggregability, which is a key parameter implicated in clotting and thrombus formation.

Analytical analysis of a vesicle tumbling under a shear flow.

Vesicles under a shear flow exhibit a tank-treading motion of their membrane, while their long axis points with an angle <pi/4 with respect to the shear stress if the viscosity contrast between the interior and the exterior is not large enough. Above a certain viscosity contrast, the vesicle undergoes a tumbling bifurcation, a bifurcation which is known for red blood cells. We have recently presented the full numerical analysis of this transition. In this paper, we introduce an analytical model that has the advantage of being both simple enough and capturing the essential features found numerically. The model is based on general considerations and does not resort to the explicit computation of the full hydrodynamic field inside and outside the vesicle.

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.