REVIEW 3 major objections 6 minor 13 references
Red blood cell partitioning and segregation through vascular bifurcations in a model of sickle cell disease
T0 review · 3 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read Simulations show stiff sickle cells prefer the low-flow branch at vascular bifurcations, reversing the Zweifach-Fung effect for normal red cells.
desk verdict Stiff sickle cells preferring the low-flow branch is a plausible and nicely explained qualitative result, but the extreme 20%-volume cell model and lack of error bars cap how far the quantitative claims should be trusted. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the separatrix, the dividing surface in the parent vessel that marks whether a cell's trajectory ends in the left or right daughter branch, together with the cell-free layer (CFL), the near-wall region depleted of healthy red cells. The argument works because healthy cells are concentrated at the channel center, sickle cells are marginated in the CFL, and shifting the separatrix under unequal flow sends the CFL preferentially into the low-flow branch. The sickle cell itself is modeled as a stiff capsule with a curved prolate spheroidal rest shape, roughly 20% of normal cell volume and five times the membrane shear modulus.
What would settle it
Measure, in a microfluidic Y-junction with physiological flow splits and a suspension of stiffened or sickled red cells, the cell partition ratio into the high-flow branch: if the ratio equals or exceeds the fluid partition ratio for the high-flow branch at flow splits of 0.7 and 0.9, the central claim is contradicted. A complementary observation would be imaging the parent-vessel cell distribution and showing that stiff cells do not concentrate in the cell-free layer upstream of the bifurcation.
Extended reading notes
Core claim
The central claim is that margination inverts the Zweifach-Fung effect for aberrant red blood cells. In a simulated binary suspension of 90% normal and 10% sickle cells flowing through a symmetric bifurcation, normal cells follow the classical rule (cell partition ratio exceeds fluid partition ratio when a branch receives more flow), while sickle cells show an anti-Zweifach-Fung effect: they preferentially enter the low-flow branch. The mechanism is spatial: upstream of the fork, normal cells occupy the channel center while sickle cells reside in the cell-free layer, and the separatrix—the surface dividing fluid that enters each daughter branch—shifts toward the high-flow side when flow is asymmetric, sweeping the near-wall cell-free layer into the low-flow branch. The segregation persists downstream, with sickle cells accumulating on the outer walls of daughter branches, and this accumulation increases the probability of large wall-shear-stress events, particularly on the outer side of the high-velocity branch. In geometrically asymmetric bifurcations, cells of both types favor the larger-radius branch at equal flow, and sickle cells still exhibit the anti-Zweifach-Fung effect.
Load-bearing premise
Everything hinges on the idealized sickle cell: stiff, small, and curved enough to marginate strongly into the cell-free layer; if real sickle cells in a patient are softer or heterogeneous enough not to marginate, the anti-Zweifach-Fung effect and the wall-stress consequences would weaken or disappear.
Editorial extensions
If this is right
- In mixed sickle-cell blood, downstream branches with lower flow will receive a disproportionate share of sickle cells, concentrating stiff cells where flow is already sluggish.
- The outer walls of daughter branches, especially the high-velocity branch, see more high wall-shear-stress events in sickle-cell suspensions than in healthy suspensions.
- Geometric asymmetry alone biases cell entry toward the larger-radius branch even when flow rates are equal, and this bias compounds the flow-driven partitioning.
- The anti-Zweifach-Fung effect provides a purely physical, non-adhesive mechanism connecting margination to endothelial stress, complementing experiments showing that sickle-cell margination upregulates an endothelial dysfunction marker.
Reading between the lines
- If real sickle-cell populations contain many cells that only partially stiffen, the magnitude of the anti-Zweifach-Fung effect may be weaker or distributed across a spectrum of partition ratios, since the current model uses a single idealized sickle cell.
- The same separatrix-shift logic should apply to other marginated blood elements such as platelets, white cells, or rigid inclusions, so bifurcations may generally sort stiff or small particles into low-flow branches—a possible design principle for microfluidic separators.
- A testable extension is to measure the cell partition ratio at varying flow splits in a microfluidic Y-junction with stiffened or sickled cells; the crossover where the cell partition ratio drops below the fluid partition ratio for the high-flow branch would confirm the mechanism.
- The wall-shear-stress increase is reported as event frequency at fixed points; mapping these events to endothelial calcium signals or adhesion-molecule expression would connect the physical mechanism to tissue-level pathology.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents immersed-boundary simulations of binary red blood cell (RBC) suspensions in model vascular bifurcations, with one population of normal biconcave capsules and one population of stiff, smaller, curved prolate spheroidal capsules intended to represent sickle RBCs. The authors report that normal RBCs reproduce the classical Zweifach-Fung effect, preferentially entering the higher-flow branch, while the model sickle cells instead preferentially enter the lower-flow branch—an 'anti-Zweifach-Fung effect'—which they attribute to the marginated position of the stiff cells in the cell-free layer coupled with the flow-dependent shift of the separatrix. Downstream of the bifurcation, sickle cells accumulate along the outer branch walls, and the simulations show that the SCD suspension produces a higher probability of elevated wall shear stress events than a healthy suspension, especially on the outer side of high-velocity branches. An additional set of simulations in a geometrically asymmetric bifurcation shows that cells preferentially enter the larger-radius branch and that the anti-Zweifach-Fung effect persists.
Significance. If the central finding is robust, it is novel and potentially clinically relevant: it provides a mechanistic, purely hydrodynamic explanation for how a stiff, marginating RBC subpopulation could be routed to low-flow branches in the microcirculation and how this could locally amplify wall shear stress fluctuations, contributing to endothelial dysfunction in sickle cell disease. The paper's strengths include validation of the healthy-cell partitioning behavior against published experimental microcirculation data (Fig. 1B, Ref. 61), a clear mechanistic account of the partitioning via the separatrix and the cell-free layer, and the use of an established simulation method (the authors' prior Science Advances work, Ref. 18). The anti-Zweifach-Fung effect emerges from the simulations rather than being imposed as an input, which is a genuine model output. However, the sickle cell is represented by a single, strongly marginating parameter set, and the quantitative claims are not supported by sensitivity analysis or uncertainty quantification, so the breadth of the physiological conclusion currently outruns the evidence.
major comments (3)
- [Formulation] The sickle cell model uses a rest volume approximately 20% that of the normal RBC and a membrane shear modulus five times larger. This parameter choice is physiologically extreme: reported sickle RBC volumes are typically 60–80% of normal, not 20%. Because the entire anti-Zweifach-Fung mechanism rests on the near-wall (cell-free layer) location of sickle cells, and because smaller, stiffer particles marginate more strongly, this extreme parameter set likely inflates the reported effect. The authors themselves call the single-class representation 'a substantial simplification,' yet they provide no sensitivity analysis in cell volume, stiffness ratio, or aberrant-cell fraction. Without such a test (e.g., repeating the key partitioning simulations at 50% and 75% of normal volume, or at stiffness ratios of 2 and 3), the reader cannot tell whether the anti-Zweifach-Fung effect is a robust feature of diseased RBCs or an artifact of the chosen idealized parameter point. This is the load-bearing assumption for the paper's headline result, so the omission is a major gap.
- [Results, Partitioning and Wall Shear Stress] No uncertainty quantification or replicate statistics are provided for any of the central quantitative claims. Figures 2(B), 4(D–F), and 5(B) show smooth curves and probability densities from what appear to be single simulations per condition, with no error bars, confidence bands, or indication of run-to-run variability. Given that cell trajectories through a bifurcation are stochastic, the statements that ηN for sickle cells lies below ηQ for ηQ > 0.5 (and above for ηQ < 0.5), and that the SCD suspension raises the probability of high WSS events, cannot be assessed for statistical significance as presented. At a minimum, the authors should report the number of independent simulations or the length of the stationary time window, and provide bootstrap or replicate-based error estimates on ηN and on the WSS probability distributions. This is necessary to support the quantitative, and especially the inverted-branch-preference, claim.
- [Results, Geometrically Asymmetric Bifurcation] The asymmetric-bifurcation result that cells preferentially enter the larger-radius branch even at equal volumetric flow (ηQ = 0.5) is presented without any sensitivity to the degree of geometric asymmetry. Only one pair of daughter radii (14 µm and 11 µm) is reported, and the claim that 'the geometric asymmetry of the bifurcation increases the curvature of the separatrix, further amplifying the uneven distribution' is not backed by a systematic variation of the radius ratio or bifurcation angle. Because the separatrix curvature is the proposed explanation, the lack of a parameter sweep leaves the generality of this secondary conclusion unsupported, even if the symmetric-bifurcation mechanism is accepted.
minor comments (6)
- [Discussion] The word 'modling' in the first paragraph of the Discussion should be 'modeling'; the same paragraph also contains a missing article ('this is an anti-Zweifach-Fung effect' reads awkwardly in context).
- [Formulation] The sentence 'The Chorin projection method is utilized to advance the velocity field u. This method involves solving an advection-diffusion equation...' repeats the earlier statement about the projection method; the redundancy should be removed.
- [Results, Validation] In Fig. 1(B), the experimental data from Ref. 61 are plotted without error bars; reproducing the experimental uncertainty would help the reader judge the agreement quantitatively.
- [Results, Wall Shear Stress] The claim that 'the average WSS at point a remains close to zero' should be supported by a numerical value, since the distribution in Fig. 4(D) is broad and symmetric but the mean is not reported.
- [Results, Cell Distribution and Segregation] The 'regions of interest' (ROI) in Fig. 3(A) are only defined in the caption; the main text should state the axial extent of the ROI and how the cross-sectional number density is normalized.
- [Formulation] The verification at Re = 0.05 is mentioned only as a sentence without any results; a supplementary figure or a quantitative statement of the change in ηN would make this verification credible.
Circularity Check
No significant circularity; the anti-Zweifach-Fung effect is an emergent simulation result, and the model is validated against independent experimental data.
full rationale
The paper is a computational experiment, not a derivation. The sickle-cell model inputs (rest volume ~20% of normal, fivefold membrane stiffness) are physiological assumptions, not the target outcome. The headline anti-Zweifach-Fung effect is measured from simulation outputs (ηN vs ηQ) and is not encoded in the inputs; it emerges from the simulated marginated distribution interacting with the flow separatrix. The healthy-cell case is validated against independent experimental results (Pries et al., ref. 61; Fig. 1B), breaking any self-referential loop. The separatrix analysis is a post-hoc logistic-regression diagnostic fitted to trajectory data; it is used to explain the measured partitioning, not to generate the predicted partition ratios, so it does not reduce to a fitted-parameter-as-prediction step. Self-citations to the authors' prior work (ref. 18) provide the numerical method and previous validation, but the current paper contains its own external benchmark, so those citations are not load-bearing. No uniqueness theorem or ansatz is smuggled in via self-citation, and the 'anti-Zweifach-Fung effect' naming is contextualized by earlier observations of low-flow-branch preference for stiff cells (refs. 48, 62). No step reduces by construction to its own inputs.
Assumptions & free parameters
free parameters (5)
- Sickle cell membrane stiffness ratio =
5x normal RBC
- Sickle cell volume ratio =
0.2x normal RBC
- Aberrant cell number fraction =
0.1
- Viscosity ratio (internal to plasma) =
1
- Particle Reynolds number Rep =
0.1
assumptions (4)
- standard math Plasma is an incompressible Newtonian fluid obeying the Navier-Stokes equations.
- domain assumption Normal RBCs can be represented as deformable biconcave capsules with shear and bending elasticity.
- ad hoc to paper Sickle RBCs can be represented as a single class of stiff, smaller prolate spheroidal capsules.
- domain assumption Inflow/outflow single-phase zones with Dirichlet velocity boundary conditions produce fully developed particulate flow at the bifurcation.
Cite this review
Pith. "Pith review of Red blood cell partitioning and segregation through vascular bifurcations in a model of sickle cell disease." pith.science (2026). https://pith.science/paper/LBSZDHOW
@misc{pith2026250104017,
author = {Pith},
title = {Pith review of: Red blood cell partitioning and segregation through vascular bifurcations in a model of sickle cell disease},
year = {2026},
howpublished = {\url{https://pith.science/paper/LBSZDHOW}},
note = {Machine review of arXiv:2501.04017}
}
read the original abstract
The impact of cell segregation and margination in blood disorders on microcirculatory hemodynamics within bifurcated vessels are physiologically significant, yet poorly understood. This study presents a comprehensive computational investigation of red blood cell (RBC) suspension dynamics, with a focus on a model of sickle cell disease (SCD) as an example of a disorder associated with subpopulations of aberrant RBCs. The findings reveal how cell margination influences cellular partitioning and distributions as well as vessel wall shear stress (WSS) at vascular bifurcations. Normal RBCs, which migrate toward the channel center, exhibit the Zweifach-Fung effect, preferentially entering high-flow-rate branches. In contrast, sickle cells, which marginate near the vessel wall, demonstrate an anti-Zweifach-Fung effect, favoring lower-flow-rate branches due to their position within the cell-free layer (CFL). The upstream segregation of cells remains downstream through the bifurcation, where sickle cells accumulate along the outer branch walls. This accumulation of sickle cells increases the frequency of high WSS events via direct physical interactions, particularly on the outer side of high-velocity branches, potentially contributing to the vascular damage and endothelial disruption observed in many disorders that affect RBCs. In geometrically asymmetric bifurcations, cells preferentially enter branches with larger radii, underscoring the influence of geometric complexity on microcirculatory blood flow. These findings provide insights into microvascular hemodynamics in SCD and other blood disorders.
Reference graph
Works this paper leans on
-
[1]
Red blood cell partitioning and segregation through vas- cular bifurcations in a model of sickle cell disease Xiaopo Cheng, a Christina Caruso , b Wilbur A. Lam, bc and Michael D. Graham ∗a The impact of cell segregation and margination in blood disorders on microcirculatory hemodynamics within bifurcated vessels are physiologically significant, yet poorl...
work page Pith review arXiv 2024
-
[10]
10 A. Kumar and M. D. Graham, Physical Review Letters , 2012, 109, 108102. 11 V. Narsimhan, H. Zhao and E. S. G. Shaqfeh, Physics of Fluids, 8 | 1–10 + P V S O B M / B N F < Z F B S > < W P M > 2013, 25, 061901. 12 R. G. Henríquez Rivera, K. Sinha and M. D. Graham, Physical Review Letters, 2015, 114, 188101. 13 R. G. Henríquez Rivera, X. Zhang and M. D. G...
work page 2012
-
[18]
To gain a comprehensive under- standing of the hydrodynamic ramifications associated with cell margination and segregation at the vascular surface, we conduct an evaluation of the WSS near the bifurcation. As depicted in Fig. 4 (A), illustrating the transient distribution of WSS on the bifurcated vessel surface, distinct positive local WSS fluctuations ar...
work page 1971
- [36]
-
[249]
57 K. Sinha and M. D. Graham, Physical Review E , 2015, 92, 042710. 58 E. Evans and Y.-C. Fung, Microvascular Research , 1972, 4, 335–347. 59 P. Balogh and P. Bagchi, Journal of Computational Physics , 2017, 334, 280–307. 60 R. Mittal, H. Dong, M. Bozkurttas, F. Najjar, A. Vargas and A. Von Loebbecke, Journal of Computational Physics , 2008, 227, 4825–485...
work page 2015
-
[425]
43 C. Sun and L. L. Munn, Computers & Mathematics with Appli- cations, 2008, 55, 1594–1600. 44 S. S. Shevkoplyas, T. Yoshida, L. L. Munn and M. W. Bitensky , Analytical chemistry, 2005, 77, 933–937. 45 H. Wei Hou, H. Y. Gan, A. A. S. Bhagat, L. D. Li, C. T. Lim and J. Han, Biomicrofluidics, 2012, 6, 024115. 46 R. Enjalbert, D. Hardman, T. Krüger and M. O....
work page 2008
- [462]
-
[712]
16 X. Zhang, C. Caruso, W. A. Lam and M. D. Graham, Physical Review Fluids, 2020, 5, 053101. 17 C. Caruso, M. E. Fay , X. Cheng, A. Y. Liu, S. I. Park, T. A. Sulchek, M. D. Graham and W. A. Lam, iScience, 2022, 25, 104606. 18 X. Cheng, C. Caruso, W. A. Lam and M. D. Graham, Science Advances, 2023, 9, eadj6423. 19 A. Pries, T. W. Secomb and P. Gaehtgens, C...
work page 2020
Show all 13 references
-
[845]
Klitzman and B
51 B. Klitzman and B. R. Duling, American Journal of Physiology- Heart and Circulatory Physiology, 1979, 237, H481–H490. 52 I. H. Sarelius and B. R. Duling, American Journal of Physiology-Heart and Circulatory Physiology , 1982, 243, H1018–H1026. 53 Z. L. Liu, J. R. Clausen, J...
1979
-
[1700]
8 Q. M. Qi and E. S. G. Shaqfeh, Physical Review Fluids, 2017, 2, 093102. 9 K. Mueller, D. A. Fedosov and G. Gompper, Medical Engineer- ing & Physics, 2016-01, 38, 2
2017
-
[3560]
Caruso, X
4 C. Caruso, X. Cheng, M. E. Michaud, H. M. Szafraniec, B. E. Thomas, M. E. Fay , R. G. Mannino, X. Zhang, Y. Sakurai, W. Li et al., Blood, 2024, 144, 2050–2062. 5 A. Kumar and M. D. Graham, Physical Review E , 2011, 84, 066316. 6 A. Kumar, R. G. Henríquez Rivera and M. D. Gra...
2024
-
[4304]
Rashidi, O
38 Y. Rashidi, O. Aouane, A. Darras, T. John, J. Harting, C. Wag- ner and S. M. Recktenwald, Soft Matter , 2023, 19, 6255–
2023
-
[6266]
Ye and L
39 T. Ye and L. Peng, Physics of Fluids, 2019, 31, 021903. 40 T. Ye, L. Peng and G. Li, Biomechanics and modeling in mechanobiology, 2019, 18, 1821–1835. 41 M. Sugihara-Seki, T. Onozawa, N. Takinouchi, T. Itano and J. Seki, Biorheology, 2021, 57, 101–116. 42 C. Bächer, A. Kihm...
2019
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.