---
title: Blood Symmetries in Flow Dynamics
url: https://www.emergentmind.com/topics/blood-symmetries
type: topic
---

# Blood Symmetries in Flow Dynamics

Blood symmetries comprise the geometric invariances, orientational order, spatial regularities, and symmetry-breaking transitions that arise when red blood cells (RBCs), plasma, and interfaces evolve under flow, confinement, and evaporation. In this broad sense, symmetry in blood is not limited to static morphology. It includes axisymmetric membrane equilibria, preferred orientations in shear, crystal-like ordering under confinement, radial segregation in porous media, reproducible axisymmetric drying patterns, and a constitutive universality that appears when stress and strain rate are conditioned on local RBC concentration rather than on global hematocrit or channel width [2507.21010], [1010.6196], [1711.08772], [2408.13824].

## 1. Elastic-surface geometry and the symmetry of the single RBC

A central geometric formulation treats the RBC membrane as an equilibrium elastic surface governed by the Helfrich–Canham functional,
\[
\Theta(M)=\iint_M \beta(2H - c_0)^2 \, dS + \lambda \int dA + \Delta P \int dV,
\]
where \(\beta>0\) is the bending rigidity, \(c_0\) is the spontaneous curvature, \(H = -\frac{k_1+k_2}{2}\) is the mean curvature, and \(\lambda\) and \(\Delta P\) are Lagrange multipliers enforcing fixed area and fixed enclosed volume. Critical points satisfy the Helfrich shape equation
\[
2 \Delta H + (2H - c_0)\left[2H^2 - 2K + c_0 H\right] + \bar{P} - 2\bar{\lambda} H = 0,
\]
with \(\bar{P} := \Delta P/\beta\) and \(\bar{\lambda} := \lambda/\beta\). In this framework, the RBC is a closed axisymmetric elastic surface of spherical topology rather than a rigid particulate object [2507.21010].

For a surface of revolution generated by a profile curve \((r(s),z(s))\), the parameterization
\[
X(s,v)=\big(r(s)\cos v,\ r(s)\sin v,\ z(s)\big)
\]
reduces the problem to an axisymmetric ODE in the tangent angle \(\psi(r)\). This axisymmetric reduction is technically important because the full Helfrich equation is a nonlinear fourth-order elliptic PDE. The paper emphasizes that axisymmetry is a useful structural restriction, but not every visually plausible axisymmetric profile is an exact equilibrium.

That distinction is made explicit for Cassinian ovals. Using
\[
z(r)=\pm\sqrt{\sqrt{4\epsilon^2 r^2+1}-\epsilon^2-r^2},
\]
the paper proves that Cassinian ovals do not satisfy the Helfrich shape equation for \(\epsilon>0\). The only exception is the limiting case \(\epsilon=0\), where the profile degenerates to a circle and the surface of revolution is a round sphere. The geometric consequence is precise: Cassinian ovals may serve as effective approximations, but they are not exact Helfrich equilibria except in the spherical limit [2507.21010].

This establishes an important boundary condition for the subject. “Blood symmetry” at the single-cell level is not merely the visual symmetry of a biconcave outline. It is the curvature-constrained symmetry of an elastic membrane subject to area and volume constraints. A common overidentification of RBC symmetry with a convenient planar curve is therefore not supported by the elastic-surface model.

## 2. Shear alignment, nematic order, and rotational symmetry under flow

In suspension, RBC symmetry becomes dynamical. A mesoscale lattice Boltzmann study models plasma with a D3Q19/BGK method and each RBC as a rigid ellipsoid coupled hydrodynamically by bounce-back with correction for moving boundaries and interacting through a soft, anisotropic repulsive cell–cell potential. This coarse-grained framework preserves the particulate character of blood while allowing simulations of suspensions of order \(10^6\) cells; the orientation study itself is reported for about \(3\times10^4\) cells at \(\Phi=45\%\) [1010.6196].

In Couette flow, with shear gradient along \(x\), flow along \(z\), and vorticity axis along \(y\), RBCs do not remain perfectly aligned with the velocity direction. They develop a preferred nonzero inclination angle \(\theta\) relative to the flow, with the sign convention chosen so that \((\hat{\mathbf{o}}_i)_x>0\). At physiological hematocrit, the most probable inclination decreases from \(\theta^* \approx 14.5^\circ\) to \(\theta^* \approx 3.4^\circ\) as the shear rate increases from about \(300\,\mathrm{s}^{-1}\) to \(6000\,\mathrm{s}^{-1}\). Stronger shear therefore produces closer alignment with the flow but not exact parallelism [1010.6196].

The orientational statistics are quantified by the nematic order tensor
\[
S_{kl}=\frac{1}{2}\left\langle 3\hat{o}_k\hat{o}_l-\delta_{kl}\right\rangle_{i,t},
\]
whose largest eigenvalue \(\lambda_+\) is the nematic order parameter. In the reported regime, increasing shear rate decreases \(\lambda_+\). At fixed shear rate, \(\lambda_+\) varies with hematocrit; at higher hematocrit, order can remain relatively high because tumbling is hindered, while at lower hematocrit, order is reduced because there is less coherent alignment in the vorticity plane. The suspension is therefore partially ordered, not perfectly ordered, and its order cannot be reduced to a single-particle Jeffery orbit.

The rotational dynamics confirm this point. RBCs undergo intermittent flips by an angle of \(\pi\) around the vorticity direction, and the rotation period \(T\) is defined as the time between successive flips. The measured distribution of \(T\) is broad rather than sharp. At lower hematocrit it becomes narrower and shifts toward shorter times, yet its width remains comparable to the mean. The peak lies near the theoretical rotation period for a freely tumbling ellipsoidal particle, which indicates that the dominant rotational timescale remains close to single-particle shear physics even though interactions substantially broaden the dynamics [1010.6196].

A separate theoretical analysis of enucleated incompressible RBCs in shear extends this symmetry discussion from orientation to shape. In a quasi-spherical harmonic expansion, the current and stress-free membrane shapes are decomposed into spherical harmonics, and instability appears when higher modes satisfy
\[
\Pi^2+\sigma(t)\Pi+4\alpha<0,\qquad \Pi\equiv l(l+1),
\]
with critical tension
\[
\sigma_*=-4\sqrt{\alpha},\qquad \alpha=\frac{\mu R^2}{\kappa}.
\]
When the instability is triggered, excess area is redistributed from the \(l=2\) ellipsoidal sector into higher modes. From an out-of-plane initial condition, the long-time morphology can become a stomatocyte; from an in-plane initial condition, it can become a trilobe. The paper interprets these morphologies as changes in symmetry class driven by shear, bending, elasticity, viscosity ratio, and excess area [2606.19072].

Taken together, these results show that shear does not impose a single symmetry state on blood. It produces small-angle alignment, finite nematic order, broad symmetry-related tumbling statistics, and, under suitable mechanical conditions, higher-mode symmetry breaking of the cell shape itself.

## 3. Hydrodynamic self-organization and crystal-like order in confined shear

Under wall-confined shear flow, deformable RBCs can assemble into regular spatial patterns. Experiments and simulations show that healthy RBCs form flow-aligned chains at low concentration and, at higher concentration, two-dimensional lattices and triangular arrangements. The final configurations are symmetric about the flow direction, and the order is described as purely hydrodynamic and inertialess in origin [1711.08772].

The organizing mechanisms are twofold. First, deformable RBCs near confining walls experience wall-induced hydrodynamic repulsion or migration toward the channel mid-plane. Second, intercellular hydrodynamic interactions can be attractive or repulsive depending on cell–cell separation. A single tank-treading RBC in the midplane generates a quadrupolar flow field and behaves at leading order like a force dipole or stresslet with strength scaling as
\[
\Sigma \sim -\eta_0 \dot{\gamma} R^3.
\]
At larger separations this far-field structure favors alignment along the flow direction, but attraction alone is not sufficient for a stable lattice. Stability requires a balance between quadrupolar attraction, wall-induced lift, and a short-range hydrodynamic interaction that is oscillatory and exponentially decaying.

Within this analytical picture, a preferred spacing is selected. In the weak-confinement limit, the stable pair distance obeys the universal law
\[
\frac{\Delta Z_0}{W} = \frac{\pi/2+\phi}{a} \approx 1.805.
\]
The existence of a selected spacing, largely independent of detailed cell mechanics, is one of the clearest examples of hydrodynamically selected blood symmetry in the literature [1711.08772].

Deformability is decisive. Healthy deformable RBCs tank-tread, migrate to the mid-plane, and reach stable separations. Hardened RBCs in experiments and rigid particles in simulations remain disordered under the same conditions. Even when rigid particles are initially placed near the mid-plane, their mutual distance keeps increasing without saturation. The order is therefore not attributable to adhesion, electrostatics, Brownian motion, or inertia. It is a deformability-mediated Stokes-regime self-organization process [1711.08772].

This distinction also clarifies a frequent misunderstanding. Regularity in confined blood flow does not imply crystalline order in the solid-state sense, nor does it require an attractive biochemical potential. The relevant symmetry is a dynamic lattice-like organization selected by viscous hydrodynamics and confinement.

## 4. Symmetry breaking in bifurcations and porous microstructured media

Microvascular branching introduces another class of blood symmetries: the relation between geometric symmetry and partition symmetry. In two-dimensional simulations of diverging and converging bifurcated microvessels, blood plasma is treated as an incompressible Newtonian fluid and RBC membranes are modeled as spring networks coupled to the flow by the immersed boundary method. Both symmetric and asymmetric bifurcations are considered, with RBC counts corresponding to hematocrits of approximately \(3.2\%\), \(6.4\%\), \(17.6\%\), and \(32\%\) [1602.07044].

In a symmetric bifurcation, the imposed geometry and initial flow field are symmetric, yet RBC partitioning can still be unequal because cell deformation and cell–cell interactions break the symmetry dynamically. For 8 cells at Hct \(3.2\%\), the paper reports that 5 cells entered the upper branch and 3 the lower branch. This asymmetry is associated with asymmetric deformation below a critical reduced area \(s^* = 0.736\). The simulations also show lateral migration, blunted velocity profiles in straight and branched sections, and a cell-free layer near the vessel wall whose thickness decreases as hematocrit increases [1602.07044].

In an asymmetric bifurcation, the branch with the higher flow rate receives disproportionately more cells. For 8 cells at Hct \(3.2\%\), essentially all RBCs enter the larger branch. As source hematocrit rises, the smaller branch receives more cells than in the dilute case, but the larger branch remains dominant. Near the diverging apex, the velocity profile develops two asymmetric peaks; after reconvergence it becomes single-peaked again. A cell-free region is also observed around the tip of the confluence, and this depleted zone shrinks as hematocrit increases [1602.07044].

A different but related symmetry-breaking phenomenon appears when RBC suspensions wick through fibrous porous media. In filter paper, the stain can separate into a dense RBC-rich core surrounded by a PBS annulus that is almost free of RBCs. This core–ring structure is described as phase separation between RBCs and suspending fluid and represents a spontaneous radial organization emerging in a statistically disordered porous matrix [2509.01408].

The reported controls are hematocrit, aggregation strength, membrane flexibility, and confinement. For healthy cells on Grade 1 paper, the phase-separated area decreases with HCT and disappears above roughly \(40\%\) HCT. In simulations, a cross-over in relative mean speed occurs for healthy RBCs at about \(29\%\) HCT, beyond which the suspending fluid no longer outruns the cells. Larger pores reduce the range over which phase separation is seen. Increasing dextran concentration suppresses separation by producing larger and more complex aggregates, whereas rigidified RBCs show a much stronger phase-separated pattern that persists even at very high hematocrits [2509.01408].

These studies show that blood-flow symmetry is not simply a property of external geometry. A symmetric vessel can yield asymmetric partitioning, while a disordered porous substrate can yield a robust radial core–annulus symmetry. This suggests that deformability, crowding, and transport resistance can select macroscopic organization even when the underlying geometry does not prescribe it directly.

## 5. Axisymmetry and impact morphology in blood droplets

At free surfaces, blood exhibits both preserved and altered symmetries. During evaporation of a drop of human blood, the final stage reveals a regular pattern with good reproducibility for a healthy person. The same axisymmetric pattern formation is observed, and can be forecast for different blood drop diameters. The paper attributes the evaporation process to Marangoni flow only, with evaporation, wettability, and RBC motion at the drop periphery shaping the final deposit. The evaporation mass flux can be predicted with good agreement assuming only the knowledge of the colloids mass concentration [1010.2510].

The axisymmetry is therefore not treated as incidental. It is the organized outcome of Marangoni-driven internal circulation acting in a radially distributed geometry. Anaemic and hyperlipidemic samples produce different final structures, but these are presented as modifications of the same evaporation-driven process rather than wholly separate regimes. The dried pattern can thus be read as a frozen record of fluid mechanics and hematological composition [1010.2510].

Impact on a solid surface probes a different interfacial symmetry problem: the radial spreading and splash morphology of a rapidly deforming blood drop. Comparative experiments on glass use whole blood (WB), plasma with hard particles (PwP), glycerol-water with hard particles (GWwP), and a commercial blood simulant (BS). The ranges are \(550<Re<1700\) and \(120<We<860\), with impact velocities \(V=2.0\)–\(4.5\ \mathrm{m/s}\). Despite similar effective high-shear viscosities of about \(5\ \mathrm{mPa\cdot s}\), the liquids do not behave equivalently [2201.06673].

| Liquid | Spreading and splashing signature | Threshold note |
|---|---|---|
| WB | Largest spreading; only finger-splashing | Finger-splashing \(We^\ast=350\pm25\); predicted \(381\pm9\) |
| PwP | Reduced spreading; finger-splashing and tiny-splashing | Finger \(300\pm25\); tiny \(450\pm25\) |
| GWwP | Similar to PwP; finger-splashing and tiny-splashing | Finger \(300\pm25\); tiny \(400\pm25\) |
| BS | Smallest spreading factor; only tiny-splashing | Tiny \(250\pm25\); no finger-splashing |

For WB, the maximum spreading radius follows
\[
\frac{R_{max}}{R_0} = C Re^{1/4},
\]
with fitted constant \(C=0.70\). WB therefore behaves like a Newtonian fluid in spreading under these impact conditions. However, the hard-particle suspensions and the commercial simulant spread less and exhibit particle-ejection-type tiny-splashing. The paper interprets this difference as a consequence of RBC deformability: deformable RBCs suppress the suspension-like frictional dissipation and momentum-transfer mechanisms that promote particle ejection in hard-particle systems [2201.06673].

This is a precise counterexample to the idea that high-shear viscosity alone determines impact symmetry and splash morphology. Matching viscosity is not sufficient. The data indicate that deformable particles are essential if a blood simulant is meant to reproduce the spreading and splashing of whole blood.

## 6. Constitutive universality as a non-geometric blood symmetry

Not all blood symmetries are geometric or spatial. A recent straight-channel theory identifies a constitutive universality that emerges when blood is described by the correct local state variable. The paper argues that blood flow appears non-universal when one uses position, channel width, or global hematocrit alone, but becomes universal when local stress and strain rate are conditioned on the local RBC concentration \(\bar\phi(y)\) [2408.13824].

The local fields are the concentration \(\bar\phi(y)\), shear stress \(\bar\sigma_{xy}(y)\), and strain rate \(\bar{\dot\gamma}(y)=\partial_y \bar u_x(y)\). The constitutive statement is
\[
\bar{\sigma}_{xy}(y)=\bar{\eta}(y)\,\bar{\dot{\gamma}}(y),
\qquad
\bar{\dot{\gamma}}(y)=\bar f(y)\,\bar{\sigma}_{xy}(y),
\]
with local fluidity empirically fitted by
\[
\bar f = \frac{(1-\bar\phi)\left[1+(1-\bar\phi)^2\right]}{2}.
\]
At fixed \(\bar\phi\), data from shear-driven and pressure-driven flows, and from different channel widths and global concentrations, collapse onto the same stress–strain-rate relation. In this sense, the symmetry is an invariance of constitutive behavior under changes in global setup once the local concentration is fixed [2408.13824].

The law is not closed until \(\bar\phi(y)\) is determined. The paper therefore introduces a concentration transport equation,
\[
\partial_t \bar\phi + \partial_y \bar J_y = 0,
\qquad
\bar J_y = \bar J_y^{\,l} + \bar J_y^{\,r},
\]
where the flux combines wall-induced lift migration and hydrodynamic diffusion from cell–cell interactions. A simple Fick law is found to be insufficient. The proposed non-local form
\[
\bar J_{cy}^{\,r} = -D(y)\,\partial_y \bar\phi_c(y) + \bar\phi_c(y)\int \partial_y U(y-y')\,\bar\phi_c(y')\,dy'
\]
reproduces the qualitative and often quantitative concentration profiles. The theory is restricted to straight channels, simplified cell models, and no explicit aggregation, but it gives a closed mesoscopic scheme in which local constitutive symmetry is paired with a non-local transport law for the symmetry-setting variable itself [2408.13824].

This non-geometric notion of blood symmetry unifies several of the preceding phenomena at the level of interpretation. A plausible implication is that many apparently distinct blood patterns—alignment in shear, ordering under confinement, branch bias, and phase separation—may be viewed as mechanisms that generate local concentration and microstructural fields, while the constitutive response remains universal once those fields are specified. In that formulation, blood symmetries are not only shapes and patterns; they are also invariances of description.

Source: https://www.emergentmind.com/topics/blood-symmetries