Deformation of liquid capsules enclosed by elastic membranes in simple shear flow: large deformations and the effect of fluid viscosities

Abstract

The deformation of a liquid capsule enclosed by an elastic membrane in an infinite simple shear flow is studied numerically at vanishing Reynolds numbers using a boundary-element method. The surface of the capsule is discretized into quadratic triangular elements that form an evolving unstructured grid. The elastic membrane tensions are expressed in terms of the surface deformation gradient, which is evaluated from the position of the grid points. Compared to an earlier formulation that uses global curvilinear coordinates, the triangular-element formulation suppresses numerical instabilities due to uneven discretization and thus enables the study of large deformations and the investigation of the effect of fluid viscosities. Computations are performed for capsules with spherical, spheroidal, and discoidal unstressed shapes over an extended range of the dimensionless shear rate and for a broad range of the ratio of the internal to surrounding fluid viscosities. Results for small deformations of spherical capsules are in quantitative agreement with the predictions of perturbation theories. Results for large deformations of spherical capsules and deformations of non-spherical capsules are in qualitative agreement with experimental observations of synthetic capsules and red blood cells. We find that initially spherical capsules deform into steady elongated shapes whose aspect ratios increase with the magnitude of the shear rate. A critical shear rate above which capsules exhibit continuous elongation is not observed for any value of the viscosity ratio. This behaviour contrasts with that of liquid drops with uniform surface tension and with that of axisymmetric capsules subject to a stagnation-point flow. When the shear rate is sufficiently high and the viscosity ratio is sufficiently low, liquid drops exhibit continuous elongation leading to breakup. Axisymmetric capsules deform into thinning needles at sufficiently high rates of elongation, independent of the fluid viscosities. In the case of capsules in shear flow, large elastic tensions develop at large deformations and prevent continued elongation, stressing the importance of the vorticity of the incident flow. The long-time behaviour of deformed capsules depends strongly on the unstressed shape. Oblate capsules exhibit unsteady motions including oscillation about a mean configuration at low viscosity ratios and continuous rotation accompanied by periodic deformation at high viscosity ratios. The viscosity ratio at which the transition from oscillations to tumbling occurs decreases with the sphericity of the unstressed shape. Results on the effective rheological properties of dilute suspensions confirm a non-Newtonian shear-thinning behaviour.

Wrinkling, tumbling and swinging microcapsules in simple shear flow

“…We observed the swinging mode in a shear rate range between 7.5 -19 s -1 until the oscillations disappeared. The tank-treading frequency was almost half of the oscillation frequency as predicted for purely elastic membranes [10,16,19,21].…”

Swinging of Red Blood Cells under Shear Flow

“…In particular, the model does not account for the possible elastic energy storage induced by the local deformations of the cytoskeleton during tanktreading. Approaches including membrane elasticity are either restricted to spherical resting shapes because of analytical complexities [10], or propose encouraging but still limited numerical analysis on tanktreading elastic biconcave capsules [11].Here, we reveal a new regime of motion for RBCs under small shear flow, characterized by an elastic capsule-like oscillation of the cell inclination superimposed to tanktreading that we name swinging. We develop a model, which predicts both swinging and the shear-stress dependency of the tumbling-tanktreading transition.…”

A fast algorithm for simulating vesicle flows in three dimensions

“…Once such triangulation is available, either collocation or Galerkin schemes can be used to discretize the integral equation [4]. Examples include [35,49,53,57]. High-order B-spline methods were used in [69].…”