Three-dimensional computational modeling of multiple deformable cells flowing in microvessels

Three-dimensional (3D) computational modeling and simulation are presented on the motion of a large number of deformable cells in microchannels. The methodology is based on an immersed boundary method, and the cells are modeled as liquid-filled elastic capsules. The model retains two important features of the blood flow in the microcirculation, that is, the particulate nature of blood and deformation of the erythrocytes. The tank-treading and tumbling motion and the lateral migration, as observed for erythrocytes in dilute suspension, are briefly discussed. We then present results on the motion of multiple cells in semidense suspension and study how their collective dynamics leads to various physiologically relevant processes such as the development of the cell-free layer and the Fahraeus-Lindqvist effect. We analyze the 3D trajectory and velocity fluctuations of individual cell in the suspension and the plug-flow velocity profile as functions of the cell deformability, hematocrit, and vessel size. The numerical results allow us to directly obtain various microrheological data, such as the width of the cell-free layer, and the variation in the apparent blood viscosity and hematocrit over the vessel cross section. We then use these results to calculate the core and plasma-layer viscosity and show that the two-phase (or core-annular) model of blood flow in microvessels underpredicts the blood velocity obtained in the simulations by as much as 40%. Based on a posteriori analysis of the simulation data, we develop a three-layer model of blood flow by taking into consideration the smooth variation in viscosity and hematocrit across the interface of the cell-free layer and the core. We then show that the blood velocity predicted by the three-layer model agrees very well with that obtained from the simulations.