Dynamics of deformable straight and curved prolate capsules in simple shear flow

Abstract

This work investigates the motion of neutrally-buoyant, slightly deformable straight and curved prolate fluid-filled capsules in unbounded simple shear flow at zero Reynolds number using direct simulations. The curved capsules serve as a model for the typical crescent-shaped sickle red blood cells in sickle cell disease (SCD). The effects of deformability and curvature on the dynamics are revealed. We show that with low deformability, straight prolate spheroidal capsules exhibit tumbling in the shear plane as their unique asymptotically stable orbit. This result contrasts with that for rigid spheroids, where infinitely many neutrally stable Jeffery orbits exist. The dynamics of curved prolate capsules are more complicated due to a combined effect of deformability and curvature. At short times, depending on the initial orientation, slightly deformable curved prolate capsules exhibit either a Jeffery-like motion such as tumbling or kayaking, or a non-Jeffery-like behavior in which the director (end-to-end vector) of the capsule crosses the shear-gradient plane back and forth. At long times, however, a Jeffery-like quasiperiodic orbit is taken regardless of the initial orientation. We further show that the average of the long-time trajectory can be well approximated using the analytical solution for Jeffery orbits with an effective orbit constant C _eff and aspect ratio ℓ _eff. These parameters are useful for characterizing the dynamics of curved capsules as a function of given deformability and curvature. As the capsule becomes more deformable or curved, C _eff decreases, indicating a shift of the orbit towards log-rolling motion, while ℓ _eff increases weakly as the degree of curvature increases but shows negligible dependency on deformability. These features are not changed substantially as the viscosity ratio between the inner and outer fluids is changed from 1 to 5. As cell deformability, cell shape, and cell-cell interactions are all pathologically altered in blood disorders such as SCD, these results will have clear implications on improving our understanding of the pathophysiology of hematologic disease.

FIG. 1:

Schematic of the initial orientation of a curved prolate capsule in unbounded simple shear flow.

FIG. 2:

Evolution of inclination angle ϕ for a stiff prolate spheroid (Ca = 0.06) with (a) r_p = 2.0 and (b) r_p = 4.8, respectively. Solid blue lines are predictions by Jeffery’s theory and red dashed lines are simulation results.

FIG. 3:

Trajectories on the unit sphere of (a) p and (b) n for a prolate spheroid (r_p = 4.8, Ca = 0.12) with initial orientation α = 0 (up to $\stackrel{.}{\gamma }t=2000$).

FIG. 4:

Trajectories on the unit sphere of p ((a) and (b) (y-z view)) and n ((c) and (d) (x-z view)) for a prolate spheroid (r_p = 4.8, Ca = 0.12) with initial orientation α = π/36.

FIG. 5:

The evolution of the p_z and n_z for a prolate spheroid (r_p = 4.8, Ca = 0.12) with initial orientation α = π/36.

FIG. 6:

Trajectories of p (a) and n (b) for a curved prolate capsule with initial orientation [α, β] = [0,0] at the early stage of simulation ($0<\stackrel{.}{\gamma }t<500$).

FIG. 7:

Early-stage trajectories on the unit sphere of the p (blue) and n (red) vectors of a curved prolate capsule with initial orientation [α, β] = [0,π/6] ($0<\stackrel{.}{\gamma }t<500$) (a,d), [0,π/2] ($0<\stackrel{.}{\gamma }t<100$) (b,e), and [π/6,π/3] ($0<\stackrel{.}{\gamma }t<400$) (c,f), respectively.

FIG. 8:

Time sequence images (front view) of the early-stage (a) and long-time (b) motions of a curved prolate capsule with initial orientation [α, β] = [0,π/6]. The dashed arrow in each image represents the p vector at the corresponding time spot.

FIG. 9:

Long-time trajectories on the unit sphere of the p (blue) and n (red) vectors of a curved prolate capsule with initial orientation [α, β] = [0,π/6] (a,d), [0,π/2] (b,e), and [π/6,π/3] (c,f), respectively.

FIG. 10:

Evolution of y- and z-components y_cm and z_cm of the center-of-mass position of curved prolate capsules with initial orientations [α, β] = [0,π/6] (a), [0,π/2] (b), and [π/6,π/3] ^(c).

FIG. 11:

The long-time trajectory (blue curve) of a curved prolate capsule in terms of ϕ and θ and the average trajectory (black circles). The red curve represents the fit of the average trajectory using Eq. 11. The black dashed line corresponds to ϕ = π. Δθ⁺ and Δθ⁻ denote the maximum positive and negative deviations, respectively, of the instantaneous trajectory from the mean when ϕ = π. (a) Ca = 0.12, K = 0.36, λ = 1; (b) Ca = 0.06, K = 0.39, λ = 5.

FIG. 12:

The dependencies of (a) the effective orbit constant C_eff of the long-time orbit and (b) the effective aspect ratio ℓ_eff on the deformability Ca and the degree of curvature K of the curved prolate capsules.