Rheology of a dilute two-dimensional suspension of vesicles

Abstract

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 (Green's 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.