🔬 condensed matter

Enucleated incompressible red blood cells in shear flow: theoretical analysis of shape instabilities

This paper presents a perturbative theoretical framework for enucleated, incompressible red blood cells in shear flow that extends existing models to account for excess membrane area, revealing how shape instabilities and the resulting stomatocyte and trilobe morphologies depend critically on membrane tension, material moduli, and viscosity ratios.

Original authors: Avraham Moriel, Howard A. Stone, Simon Mendez

Published 2026-06-18

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Avraham Moriel, Howard A. Stone, Simon Mendez

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: The Shape-Shifting Red Blood Cell

Imagine a red blood cell (RBC) not as a rigid coin, but as a tiny, water-filled balloon made of a stretchy, elastic skin. Its job is to squeeze through narrow capillaries to deliver oxygen. Usually, this balloon is shaped like a flattened donut (a discocyte).

However, when these cells flow through blood vessels, the fluid pushes against them. At certain speeds, instead of just stretching or spinning, these cells suddenly snap into weird, complex shapes—like multi-lobed flowers or asymmetrical blobs. Scientists have seen this happen, but they didn't fully understand why the balloon decided to change its shape so drastically.

This paper is a theoretical investigation (a mathematical simulation) that tries to figure out the rules behind these sudden shape changes.

The "Before and After" of the Cell's Life

To understand the instability, the authors had to account for a specific biological event: enucleation.

The Analogy: Imagine a stress ball that is perfectly round and tight. Now, imagine you squeeze the air out of the inside of the ball, but you don't let the rubber skin shrink. The skin is now too big for the smaller ball inside. It has "extra skin" that has to fold or wrinkle.
The Science: When red blood cells mature in the bone marrow, they kick out their nucleus (the command center). This makes the cell smaller inside, but the outer skin stays the same size. This creates "excess area"—a bit of extra slack in the membrane. The paper argues that this "slack" is a crucial ingredient that makes the cell unstable and prone to changing shape when pushed by blood flow.

The Experiment: Pushing the Balloon

The researchers created a mathematical model to simulate what happens when these cells are subjected to shear flow (imagine the cell being squeezed between two sliding plates of glass).

They looked at three main forces fighting against each other:

The Flow: The water trying to stretch and spin the cell.
The Skin's Elasticity: The cell's desire to snap back to its original shape.
The Bending Resistance: The energy required to make the skin curve sharply.

The "Tension" Tug-of-War

The core discovery of the paper revolves around membrane tension. Think of this as the tightness of the balloon's skin.

The Setup: When the cell is just floating, the tension is balanced.
The Trigger: As the flow speeds up, the cell stretches. Because of the "excess area" (the extra skin mentioned earlier), the tension in the skin drops.
The Snap: If the tension drops too low, the skin can no longer hold its smooth, ellipsoidal shape. It becomes unstable. It's like a drumhead that is so loose it starts to wobble uncontrollably.

The authors found that this instability happens when the "slack" (excess area) is high, the fluid is moving fast, and the cell's skin is flexible enough to bend but not so stiff that it resists all movement.

The Result: From Smooth to Spiky

When the instability kicks in, the cell doesn't just stretch; it redistributes its "excess area."

The Analogy: Imagine a smooth, round balloon. If you push on it hard enough while it has too much rubber, it doesn't just flatten; it might suddenly pop out a few bumps or lobes.
The Finding: The math shows that the cell shifts its shape from a simple oval (like a rugby ball) to complex, multi-lobed shapes (like a trilobe or a stomatocyte, which looks like a mouth). The "excess area" moves from the main body of the cell into these new, sharp lobes.

Why Orientation Matters

The paper also notes that how the cell starts matters.

If the cell is aligned one way (perpendicular to the flow), it might just roll like a tire (tank-treading).
If it's aligned another way (parallel to the flow), the instability is more likely to trigger, causing it to tumble or flip over chaotically.

The Bottom Line

This paper provides a mathematical "recipe" for when a red blood cell will lose its smooth shape and turn into a complex, multi-lobed blob. The key ingredients are:

The "Slack": The extra skin left over after the cell matures (enucleation).
The Push: The speed of the blood flow.
The Tension: The moment the skin gets too loose to hold its shape.

The authors conclude that this "loose skin" is the secret sauce that allows red blood cells to undergo these dramatic shape transformations, which are essential for their function but also a sign of dynamic instability in the flow.

Technical Summary: Enucleated Incompressible Red Blood Cells in Shear Flow: Theoretical Analysis of Shape Instabilities

Problem Statement
Red blood cells (RBCs) are essential for oxygen transport, relying on their ability to undergo significant deformations in flow. While rigid-body models (e.g., Jeffery orbits) successfully describe RBC dynamics at low shear rates (tumbling, tank-treading, swinging), they fail to explain the emergence of complex, multi-lobed shapes and symmetry-broken morphologies (such as stomatocytes) observed at intermediate shear rates ( $100\text{--}1000 \text{ s}^{-1}$ ). The physical mechanisms driving these shape instabilities—specifically the interplay between membrane elasticity, bending resistance, membrane tension, excess area, and viscosity ratios—remain debated. Furthermore, existing theoretical frameworks often treat the RBC as a static ellipsoid or capsule, neglecting the specific geometric changes resulting from the maturation process of enucleation, where the cell nucleus is expelled, altering the relationship between the stress-free membrane area and the enclosed volume.

Methodology
The authors employ a perturbative theoretical framework modeling the RBC as a closed, fluid-filled, incompressible membrane suspended in a Newtonian shear flow. The approach extends the linearized non-linear framework previously established by Vlahovska et al. [26] by explicitly incorporating two distinct excess areas:

$\Delta_0$ (Stress-free excess area): Corresponding to the ellipsoidal shape of the RBC prior to enucleation.
$\Delta$ (Current excess area): Corresponding to the discocyte shape after enucleation, where the volume decreases ( $V < V_0$ ) while the surface area remains conserved, leading to locked internal stresses.

The RBC shape is expanded in spherical harmonics ( $Y_{l,m}$ ), separating the current configuration $\rho$ and the stress-free configuration $\rho_0$ . The governing equations of motion are derived by balancing fluid stresses (Stokes flow) with membrane forces, which are decomposed into elastic (Skalak model), bending (Helfrich energy), and membrane tension contributions. Membrane tension ( $\Sigma$ ) acts as a Lagrange multiplier to enforce constant local area and couples the different spatial modes.

The analysis proceeds in two stages:

Reduced System ( $l=2$ ): The system is first reduced to ellipsoidal modes to revisit the dynamics of tumbling, tank-treading, and swinging, explicitly accounting for the mismatch between $\Delta$ and $\Delta_0$ .
Full Perturbative System ( $l \geq 2$ ): Additional spatial modes are introduced to derive a linear instability criterion. The stability of the system is analyzed by examining the growth rate of higher-order modes ( $l > 2$ ) as a function of dimensionless parameters: the bending number ( $B$ ), capillary number ( $C$ ), viscosity ratio ( $\lambda$ ), and the excess areas.

Key Contributions and Results

Extension of the Theoretical Framework: The study generalizes the existing perturbative approach to explicitly account for the geometric and stress state changes induced by enucleation. By distinguishing between the stress-free shape ( $\Delta_0$ ) and the current shape ( $\Delta$ ), the model captures the "memory" of the cell's maturation process, a feature absent in fluid vesicle models.
Role of Excess Area in Reduced Dynamics: In the $l=2$ subspace, the authors demonstrate that the mismatch between $\Delta$ and $\Delta_0$ renormalizes the effective elastic constant ( $S$ ). This allows for dynamical transitions (e.g., from swinging to tumbling) driven by changes in the reference shape's excess area, even at fixed shear rates. The initial orientation (in-plane vs. out-of-plane) significantly influences the resulting trajectories.
Instability Criterion and Membrane Tension: The primary contribution is the derivation of a linear instability criterion for shape transitions. The analysis reveals that shape instabilities (the emergence of multi-lobed shapes) are critically dependent on the membrane tension $\sigma(t)$ $σ (t)$ . An instability occurs when the membrane tension drops below a critical threshold ( $\sigma < \sigma^*$ $σ < σ^{*}$ ), determined by the ratio of shear elasticity to bending stiffness ( $\alpha = B/C$ $α = B / C$ ).
- The instability is initiated when the bending and elastic responses are insufficient to counteract external flow and internal stresses.
- The most unstable mode ( $l^*$ ) and the critical tension depend on $\alpha$ , the viscosity ratio $\lambda$ , and the excess areas.
Saturation and Redistribution: The study analyzes the long-time behavior following an instability. The growth of unstable modes causes an increase in membrane tension, which eventually stabilizes the system. This process leads to a redistribution of excess area from the dominant $l=2$ mode to higher-order modes (e.g., $l=3, 4$ ), resulting in a new, stable, non-ellipsoidal shape configuration.
Phase Diagrams: Numerical integration of the coupled equations validates the analytical predictions. The authors construct phase diagrams in the $(B, C)$ plane, identifying regions where specific modes ( $l=3, 4, \dots$ ) become unstable. These diagrams show that increasing the viscosity ratio $\lambda$ and the excess area $\Delta$ destabilizes RBCs at high capillary numbers (weak elasticity), promoting the emergence of complex shapes.

Significance
The paper provides a theoretical explanation for the emergence of complex RBC morphologies, such as stomatocytes and trilobe shapes, under external shear flow. By linking these instabilities to the specific mechanical consequences of enucleation (the mismatch between stress-free and current excess areas), the work offers a mechanistic understanding of how RBCs transition from simple ellipsoidal dynamics to complex, multi-lobed states. The framework identifies membrane tension as the key driver coupling spatial modes and highlights the critical role of the bending-to-shear modulus ratio and viscosity in determining the onset of shape instabilities. This theoretical picture bridges the gap between rigid-body models and full numerical simulations, offering analytical insight into the conditions under which RBCs lose their ellipsoidal symmetry.