Simulated two-dimensional red blood cell motion, deformation, and partitioning in microvessel bifurcations

Movement, deformation, and partitioning of mammalian red blood cells (RBCs) in diverging microvessel bifurcations are simulated using a two-dimensional, flexible-particle model. A set of viscoelastic elements represents the RBC membrane and the cytoplasm. Motion of isolated cells is considered, neglecting cell-to-cell interactions. Center-of-mass trajectories deviate from background flow streamlines due to migration of flexible cells towards the mother vessel centerline upstream of the bifurcation and due to flow perturbations caused by cell obstruction in the neighborhood of the bifurcation. RBC partitioning in the bifurcation is predicted by determining the RBC fraction entering each branch, for a given partition of total flow and for a given upstream distribution of RBCs. Typically, RBCs preferentially enter the higher-flow branch, leading to unequal discharge hematocrits in the downstream branches. This effect is increased by migration toward the centerline but decreased by the effects of obstruction. It is stronger for flexible cells than for rigid circular particles of corresponding size, and decreases with increasing parent vessel diameter. For unequally sized daughter vessels, partitioning is asymmetric, with RBCs tending to enter the smaller vessel. Partitioning is not significantly affected by branching angles. Model predictions are consistent with previous experimental results.