Vortical flow structures induced by red blood cells in capillaries.

OBJECTIVE: Knowledge about the flow field of the plasma around the red blood cells in capillary flow is important for a physical understanding of blood flow and the transport of micro- and nanoparticles and molecules in the flowing plasma. We conducted an experimental study on the flow field around red blood cells in capillary flow that is complemented by simulations of vortical flow between red blood cells. METHODS: Red blood cells were injected in a 10 × 12 µm rectangular microchannel at a low hematocrit, and the flow field around one or two cells was captured by a high-speed camera that tracked 250 nm nanoparticles in the flow field, acting as tracers. RESULTS: While the flow field around a steady "croissant" shape is found to be similar to that of a rigid sphere, the flow field around a "slipper" shape exhibits a small vortex at the rear of the red blood cell. Even more pronounced are vortex-like structures observed in the central region between two neighboring croissants. CONCLUSIONS: The rotation frequency of the vortices is to a good approximation, inversely proportional to the distance between the cells. Our experimental data are complemented by numerical simulations.

Numerical-experimental observation of shape bistability of red blood cells flowing in a microchannel.

Red blood cells flowing through capillaries assume a wide variety of different shapes owing to their high deformability. Predicting the realized shapes is a complex field as they are determined by the intricate interplay between the flow conditions and the membrane mechanics. In this work we construct the shape phase diagram of a single red blood cell with a physiological viscosity ratio flowing in a microchannel. We use both experimental in vitro measurements as well as 3D numerical simulations to complement the respective other one. Numerically, we have easy control over the initial starting configuration and natural access to the full 3D shape. With this information we obtain the phase diagram as a function of initial position, starting shape and cell velocity. Experimentally, we measure the occurrence frequency of the different shapes as a function of the cell velocity to construct the experimental diagram which is in good agreement with the numerical observations. Two different major shapes are found, namely croissants and slippers. Notably, both shapes show coexistence at low (<1 mm s-1) and high velocities (>3 mm s-1) while in-between only croissants are stable. This pronounced bistability indicates that RBC shapes are not only determined by system parameters such as flow velocity or channel size, but also strongly depend on the initial conditions.

Theory and algorithms to compute Helfrich bending forces: a review.

Cell membranes are vital to shield a cell's interior from the environment. At the same time they determine to a large extent the cell's mechanical resistance to external forces. In recent years there has been considerable interest in the accurate computational modeling of such membranes, driven mainly by the amazing variety of shapes that red blood cells and model systems such as vesicles can assume in external flows. Given that the typical height of a membrane is only a few nanometers while the surface of the cell extends over many micrometers, physical modeling approaches mostly consider the interface as a two-dimensional elastic continuum. Here we review recent modeling efforts focusing on one of the computationally most intricate components, namely the membrane's bending resistance. We start with a short background on the most widely used bending model due to Helfrich. While the Helfrich bending energy by itself is an extremely simple model equation, the computation of the resulting forces is far from trivial. At the heart of these difficulties lies the fact that the forces involve second order derivatives of the local surface curvature which by itself is the second derivative of the membrane geometry. We systematically derive and compare the different routes to obtain bending forces from the Helfrich energy, namely the variational approach and the thin-shell theory. While both routes lead to mathematically identical expressions, so-called linear bending models are shown to reproduce only the leading order term while higher orders differ. The main part of the review contains a description of various computational strategies which we classify into three categories: the force, the strong and the weak formulation. We finally give some examples for the application of these strategies in actual simulations.

Long-lived anomalous thermal diffusion induced by elastic cell membranes on nearby particles.

The physical approach of a small particle (virus, medical drug) to the cell membrane represents the crucial first step before active internalization and is governed by thermal diffusion. Using a fully analytical theory we show that the stretching and bending of the elastic membrane by the approaching particle induces a memory in the system, which leads to anomalous diffusion, even though the particle is immersed in a purely Newtonian liquid. For typical cell membranes the transient subdiffusive regime extends beyond 10 ms and can enhance residence times and possibly binding rates up to 50%. Our analytical predictions are validated by numerical simulations.