Research Article  Open Access
Tong Wang, Yong Tao, Uwitije Rongin, Zhongwen Xing, "A TwoDimensional Numerical Investigation of Transport of MalariaInfected Red Blood Cells in Stenotic Microchannels", BioMed Research International, vol. 2016, Article ID 1801403, 16 pages, 2016. https://doi.org/10.1155/2016/1801403
A TwoDimensional Numerical Investigation of Transport of MalariaInfected Red Blood Cells in Stenotic Microchannels
Abstract
The malariainfected red blood cells experience a significant decrease in cell deformability and increase in cell membrane adhesion. Blood hemodynamics in microvessels is significantly affected by the alteration of the mechanical property as well as the aggregation of parasitized red blood cells. In this study, we aim to numerically study the connection between celllevel mechanobiological properties of human red blood cells and related malaria disease state by investigating the transport of multiple red blood cell aggregates passing through microchannels with symmetric stenosis. Effects of stenosis magnitude, aggregation strength, and cell deformability on cell rheology and flow characteristics were studied by a twodimensional model using the fictitious domainimmersed boundary method. The results indicated that the motion and dissociation of red blood cell aggregates were influenced by these factors and the flow resistance increases with the increase of aggregating strength and cell stiffness. Further, the roughness of the velocity profile was enhanced by cell aggregation, which considerably affected the blood flow characteristics. The study may assist us in understanding cellularlevel mechanisms in disease development.
1. Introduction
Malaria infects 350–500 million people and kills more than two millions every year, mostly in Africa and other underdeveloped countries. In most of the malaria cases, red blood cells (RBCs) are parasitized by Plasmodium falciparum (falciparum), a protozoan parasite transmitted by mosquitos. When the RBCs are parasitized by falciparum, two critical effects have been observed in the infected cells, namely, the reduced deformability of the cell and increased adhesion of the RBCs to vessel endothelium and other blood cells [1]. Thus, the parasitized RBCs become stiffer (could be tenfold harder than healthy RBCs [2]) and tend to form aggregates in blood. Blood flow may be substantially affected by the change of structural and mechanical properties of the malariainfected RBCs.
Under healthy physiological conditions, red blood cells move individually with blood plasma or form temporal aggregates and they can undergo severe, reversible, large elastic deformation in blood flow. Although the formation of red blood cell aggregates could be due to other factors such as the decrease of the shear rate, increase of the hematocrit, and variation of the viscosity of the suspending plasma in the blood vessel, the malariainfected RBCs tend to form firm and irreversible aggregates. Severe aggregation of RBCs may decrease the surface area of cells contacting with the blood plasma. It also causes pressure and shear stress change in the region of the aggregation. These aggregates get stuck in the small blood vessels, blocking individual RBC from passing, thus decreasing the amount of oxygen and the nutrition transferred in human body.
Experimental measurements have been done on dynamic or rheological properties and behaviors of parasitized RBCs [3–10]. In the experiments, microfluidic channels have been employed to probe the deformability of healthy and malariainfected red blood cells. The hardening of the cell and the blockage of the channel can be qualitatively studied. However, a major difficulty in carrying out experimental investigations on malariainfected RBCs is the complexity of the microvessels. Moreover, experimental techniques are in general not fully capable of illustrating cellularlevel rheological behaviors of RBC aggregates. Because of the relatively large number of red blood cells, celltocell interactions and local rheological properties were unable to identify. Therefore, numerical simulations have been utilized as an alternative tool to study the rheology of parasitized RBCs in blood flow [11–15]. Numerical models that have been used for interpretation and prediction of mechanical properties and dynamic behavior of RBCs in malaria include dissipative particle dynamics (DPD) [13, 15], smoothed particle hydrodynamics (SPH) [14], lattice Boltzmann method (LBM) [16], and others [17]. Two or multiple RBCs are considered by numerical simulations to probe their interactions, for example, their aggregation and disaggregation in blood flow [18–22].
However, very few investigations on malariainfected RBCs have been conducted in stenotic vessels and their impact on hemodynamics in stenotic vessels is not well understood. In particular, there is a lack of studies considering celltocell interactions at microscopic scale. This paper aims to provide a qualitative analysis on the dynamics of malariainfected RBCs in stenotic microvessels. Since the malariainfected RBC gradually loses its deformability and develops adhesiveness through the three stages of infection ring, trophozoite (early trophozoite, late trophozoite), and schizont, this study considers the ring and trophozoite stages in which the cells become stiffer and more adhesive while the shape of the RBCs is not significantly modified. The RBC was simulated by a spring model and the hardening of the parasitized cell was mimicked by varying the membrane constants. The malariainfected RBCs were assumed existing as aggregates with firm or loose adhesion. By using fictitious domainimmersed boundary formalism in two dimensions, we simulated the transit of RBC aggregates through a prototypical channel with a contraction. By studying the motion and dissociation of RBC aggregates and flow characteristics, we delineated the contribution to hemodynamics from each of the three factors, hardening of the membrane, aggregation strength, and stenosis magnitude.
2. Methods
We considered a twodimensional microvessel with a symmetric stenosis and employed numerical simulations to study the rheological behavior of RBC aggregates in the blood flow. In this study, blood was assumed to be a suspension of RBCs in an incompressible, Newtonian fluid with constant density and viscosity. In order to simulate the blood flow and fluidcell interactions in this irregularshaped domain, the fictitious domain method was combined with the immersed boundary method and the RBCs have been modeled by the spring model.
2.1. Red Blood Cell Model
Two types of RBC model are widely used currently, namely, elastic membrane models [19, 23] and spring models [20, 24–27]. In this paper, we adopted the spring model introduced in [24, 28] and modeled individual RBC as cytoplasm enclosed by a membrane represented by a finite number of membrane particles connected by springs. The springs of stretch/compression and bending modulus change its length and the bending angle between two neighboring springs under external force. The elastic moduli are affected in disease conditions [29, 30]. Evidence shows that malariainfected RBC has increased stiffness, which is closely related to the cell bending modulus. In this paper, we modeled malariainfected RBC with different deformability by varying the bending constant of the spring. The shape of the RBCs was chosen with the reduced area .
2.2. Aggregation Kinetics
The hardening of the RBCs is due to the parasite of the falciparum. However, the stickiness of the cell surface and the mechanism of the aggregation still remain disputable. Two different theoretical models, namely, the bridging model [18, 31] and the depletion model [32], coexist nowadays. Although they are distinct in the cause of the aggregating force, both of them predict the force being attractive at far and repulsive at near distance. Because the investigation of the mechanism of aggregation is not within the scope of this research, we chose a simple model [20] in which aggregation of RBCs was achieved by introducing an aggregating force existing in a neighborhood of membrane particles and the aggregating force was derived from with the Morsetype potential aswhere is the energy constant corresponding to the density of the aggregating particles; is the scaling factor which relates to the thickness of the depletion layer or interactive distance; is the reference distance or zeroforce distance.
The scaling aggregation force () as a function of separation distance for different parameters is shown in Figure 1. When the value of the scaling force is negative, the force is attractive; when the value of the force is positive, the force is repulsive. The force decreases almost to zero at far enough distance. Aggregation could be increased either by increasing or by decreasing . On the other hand, decreasing in and increasing in may raise the depletion thickness and the equilibrium distance of the RBCs in aggregates.
2.3. The Fictitious DomainImmersed Boundary Scheme
The flow region we studied was a stenotic microchannel for which the regular structured mesh was not applicable at the boundary of the region. Thus, we adopted the fictitious domain method because in this method the irregularshaped domain is extended to regular shape so that simple structured mesh instead of unstructured mesh can be used, which substantially reduces computational complexity of the algorithm. The fictitious domain method and its applications to fluid flow problems have been extensively described [33, 34]. To employ the fictitious domain method, the flow region was embedded in a rectangular domain denoted by . Then the fluid flow containing RBCs was solved in the bigger domain , and the noflow condition in the solid region was treated as constraints. Therefore, the governing equations for the modeled problem were the following extended NavierStokes equations:where and are the fluid velocity and pressure anywhere in the flow; is the fluid density; is the fluid viscosity. The boundary conditions were such that, on , a noslip condition was applied and, at the inlet and outlet of the channel, a periodic flow condition was enforced. A detailed description of the solution method of (2) can be found elsewhere [33, 34]. In this study, the fluidcell interaction was dealt with by the immersed boundary method developed by Peskin et al. [20, 35].
3. Results and Discussion
We studied the hydrodynamic behavior of multiple red blood cell aggregates in a horizontal channel filled with Newtonian fluid with a symmetric stenosis at the central part of the vessel. The blood plasma density and plasma viscosity were presumed fixed values. We performed a series of simulations to study RBC deformation, flow field, and cellcell interactions as the aggregates traversed the stenotic vessel. The fluid flow was generated from left to right by a constant pressure drop. A stream of flow approaching the stenosis contracted to a high speed flow at the throat of the stenosis. Afterwards, the flow was allowed to develop fully along the straight rectangular channel. At the meantime, prelocated multiple RBC aggregates flowed with the fluid in the channel. The parameters used are given in Table 1. For the parameter used in this study, larger spring constant for RBC membrane and higher aggregation energy constant correspond to malariainfected RBCs at more severe infection stages.

3.1. Aggregation of Red Blood Cells
Aggregates of four cells used in the simulations were formed in static plasma before putting them in the flow channels. The red blood cells were placed facetoface in static plasma initially. The centertocenter distance was close enough in order for the aggregation force to take effect. When the simulation starts, red blood cells approached each other and reached an equilibrium configuration. The red blood cells were modeled by two membrane constants for two different stiffness and three aggregation strengths for different adhesiveness. Therefore, six different configurations were obtained and they are shown in Figure 2. The configuration in Figure 2(a) corresponds to healthy RBC. The configurations in Figures 2(b)–2(d) may correspond to the ring stage and early trophozoite stage infection, while Figures 2(e) and 2(f) are assumed to correspond to the late stage of trophozoite in which both the stiffness and stickiness of the cell membrane increase considerably.
(a)
(b)
(c)
(d)
(e)
(f)
3.2. Motion and Dissociation of Aggregates
The microvessel in this study is a 140 μm long, 30 μm wide twodimensional channel with symmetric stenosis formed at the central location. The geometry of the fluid domain and the initial placement of the RBC aggregates are illustrated in Figure 3. Blood flow in the stenotic vessel was driven by a constant pressure gradient imposed at the inlet and the outlet. The pressure gradient was determined such that the maximum flow velocity was about 20 cm/s (a typical value in arterioles [36]).
3.2.1. In the Channel with 40% Stenosis
The motion and dissociation of 6 RBC aggregates in a 40% stenosis microvessel have been simulated. The results for two membrane bending constants and three aggregating strengths are shown in Figures 4–9 for the initial configurations in Figure 2. The snapshots are for four time instants with the second subfigure showing the instant when the mean location of all cells was at the throat of the stenosis. When the RBCs are softer (Figures 4, 6, and 8), they experienced more deformation in the flow even after the dissociation. The aggregates with weak aggregating strength (Figures 4 and 5) were dissociated easily by the hydrodynamic force and the RBCs existed as dispersed individual cells after long enough elapsed time. The aggregates with mild aggregating strength (Figures 6 and 7) were partially disaggregated from their initial configurations. Smaller aggregates of two to three cells were observed in the simulation. However, when the aggregating strength was strong (Figures 8 and 9), the hydrodynamic force was unable to disaggregate the RBC aggregates. Aggregates underwent some deformation but however kept their initial configuration even at the stenosis. Because the velocity of the flow at the throat of the stenosis increased significantly due to the contraction, the aggregates lag behind which caught up the upstream ones. Eventually, agglomeration and rouleaux of RBCs of bigger size were formed by the aggregating force. In addition, more disturbance of the flow field around the aggregates or the cells has been noted.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
3.2.2. In the Channel with 50% Stenosis
Figures 10–15 show the motion of 6 RBC aggregates in the microvessel with 50% stenosis. The parameters and the initial configuration of the aggregates were the same as in the last section. It has been shown in the simulations that the velocity of the blood flow at the throat of the stenosis increased comparing to the 40% stenosis vessel. The disassociation of the aggregates and the deformation of the RBCs were similar to the results for the 40% stenosis when the aggregating strength was weak or mild. However, when the aggregating strength became strong, the larger agglomeration or rouleaux of RBCs were formed. The reason for this phenomenon is because that the increase of stenosis severity increased blood velocity at the stenosis. Thus, the downstream cells easily caught up the upstream cells which were slowed down by the friction of the vessel wall.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
3.2.3. In the Channel with 60% Stenosis
Simulation also has been done in the microchannel with 60% stenosis and the results are shown in Figures 16–21. As the stenosis severity increased to 60%, the microvessel was blocked more than the last two cases. Unlike the cases of 40% and 50% stenosis, the aggregates with mild aggregating strength did not dissociate completely even when the elapsed time was long enough. On the other hand, they almost kept their original configuration for a long time (Figures 18(c) and 19(c)) before some of them disassociating into smaller aggregates or individual cells. Furthermore, aggregates formed by RBCs with strong aggregating strength attracted together and more compact aggregates were observed than in 50% stenosis case. It is reasonable to conclude that the increase of stenosis severity facilitated aggregation of malariainfected RBCs. More disturbance of the flow at the throat of the stenosis was also observed.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
In general, for all three stenosis levels, healthy RBCs and aggregates with ring stage or early trophozoitestage infection (configurations in Figures 2(b)–2(d)) passed through the stenosis easily and dissociate into individual cells or smaller aggregates. They also recovered their normal biconcave shape quickly after passage through the constriction. However, configurations in Figures 2(e) and 2(f) maintained and even formed larger aggregates in the stenotic vessel. They tended to block the flow at the mouth of the stenosis and this can be seen from the decrease of the flow velocities in the straight section of the channel. Thus, the delivery of nutrients and removal of toxins by RBCs will be severely decreased and the phenomena show the behavior of latestage trophozoite. These results are qualitatively consistent with the in vivo experimental findings [3].
3.3. Flow Characteristics
3.3.1. Effect of Stenosis Magnitude
In Figures 22(a)–22(f), velocity at the stenosis verses radial location for three aggregation forces and two cell membrane constants are shown in the same panel for three stenosis magnitudes. The comparison shows that the velocity increased as the stenosis increased. This trend was more profound for the soft cells with mild and strong aggregating strength. It can also be seen from Figure 22 that velocity profiles were more distorted when the aggregating strength was strong. The asymmetry of the velocity profiles was mainly due to the initial asymmetric location of the aggregates.
(a)
(b)
(c)
(d)
(e)
(f)
3.3.2. Effect of Aggregation Force
Next we investigated the effect of aggregation force on the velocity profile at the throat of stenosis (Figure 23). It is noted that the flow velocity at the throat of the stenosis decreased as the aggregation force increased. This trend was found for all the stenosis magnitudes and both the soft and the rigid red blood cells. However, it is also observed that the decrease of the velocity was not linear. When the aggregation force increased from weak (blue lines) to mild (red lines), the results showed little to no effect on the velocity. While when the aggregation force became strong (black lines), the effect was significant. The results indicated that increasing of aggregation strength up to a certain level leaded to increase of flow resistance in the stenotic vessels.
(a)
(b)
(c)
(d)
(e)
(f)
3.3.3. Effect of Cell Deformability
The cell membrane deformability also had an effect on the velocity of blood flow at the stenosis and the results are revealed in Figures 24(a)–24(i). Two membrane constant values have been used in the simulations and two types of RBCs were modeled with these constants, namely, soft cells and rigid cells. In Figure 24, black lines represented the velocity profile at the throat of the stenosis of the blood flow with soft RBCs, and red lines are for the blood flow with rigid cells. Overall, the flow velocity decreased with the increase in stiffness of the cell. In particular, the effect was more profound when the stenosis was more severe or when the aggregation force was stronger. Figure 24(a) shows the velocity profile for the case of 40% stenosis and the weak aggregation force. The velocity profile was relatively smooth and less disturbance was observed on the curve. On the contrary, the velocity profiles were rougher in narrower stenotic vessels or when the aggregation of the RBCs was more severe.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)
4. Conclusions
The main objective of this work is to probe how changes in the cell membrane rigidity, the aggregating strength, and the magnitude of the stenosis affect the transit of the malariainfected RBC aggregates through stenotic microchannels. Two types of RBCs, namely, soft and rigid, have been studied under three different levels of aggregating strength. The simulations were performed over a range of stenosis magnitudes: 40%, 50%, and 60%.
It has been found in this study that soft RBCs with strong aggregating force formed the most compact aggregates. In the aggregates, the RBCs experienced large deformation. Rigid ones with weak aggregating force formed the loosest aggregates. In these aggregates, the RBCs maintained their normal biconcave shape. In blood flow, the aggregates of weak aggregating strength were more likely to dissociate to individual cells, while the aggregates with strong aggregating strength would form even larger aggregates or rouleaux. The rigidity of the cell membrane hindered passage of the cell through microchannels and this result is qualitatively in agreement with experimental findings on the rheological behaviors of malariainfected RBCs through a narrow constriction in a microchannel. The flow velocity at the throat of the stenosis decreased with increasing membrane modulus and more deviation from the parabolic profile was observed, especially for the narrower stenosis and stronger aggregating force.
The present study used a twodimensional model to qualitatively simulate malariainfected RBCs in stenotic microvessels. It is important to extend this algorithm to the physiologically relevant hematocrit contents in threedimensional situation and compare the results with experimental observations quantitatively. It is also quite interesting to adapt the model for the more thorough investigation of RBC rheology at various stages of malaria infection as well as its impact on leukocyte migration. In addition, it provides a potential approach to investigate drug delivery at cellular level in microvessels involving infected RBCs.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
The authors acknowledge the support of the National Natural Science Foundation of China (11671203) and the Foundation for Graduate Student Innovation at Nanjing University of Aeronautics and Astronautics, China.
References
 L. H. Miller, D. I. Baruch, K. Marsh, and O. K. Doumbo, “The pathogenic basis of malaria,” Nature, vol. 415, no. 6872, pp. 673–679, 2002. View at: Publisher Site  Google Scholar
 C. T. Lim, “Single cell mechanics study of the human disease malaria,” Journal of Biomechanical Science and Engineering, vol. 1, no. 1, pp. 82–92, 2006. View at: Publisher Site  Google Scholar
 J. P. Shelby, J. White, K. Ganesan, P. K. Rathod, and D. T. Chiu, “A microfluidic model for singlecell capillary obstruction by Plasmodium falciparuminfected erythrocytes,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 25, pp. 14618–14622, 2003. View at: Publisher Site  Google Scholar
 S. Kim, A. S. Popel, M. Intaglietta, and P. C. Johnson, “Aggregate formation of erythrocytes in postcapillary venules,” American Journal of Physiology  Heart and Circulatory Physiology, vol. 288, no. 2, pp. H584–H590, 2005. View at: Publisher Site  Google Scholar
 M. Abkarian, M. Faivre, and H. Stone, “Red blood cell dynamics, deformation and separation in microfluidic devices,” Journal of Biomechanics, vol. 39, no. 1, p. S332, 2006. View at: Publisher Site  Google Scholar
 S. S. Lee, Y. Yim, K. H. Ahn, and S. J. Lee, “Extensional flowbased assessment of red blood cell deformability using hyperbolic converging microchannel,” Biomedical Microdevices, vol. 11, no. 5, pp. 1021–1027, 2009. View at: Publisher Site  Google Scholar
 H. Bow, I. V. Pivkin, M. DiezSilva et al., “A microfabricated deformabilitybased flow cytometer with application to malaria,” Lab on a Chip, vol. 11, no. 6, pp. 1065–1073, 2011. View at: Publisher Site  Google Scholar
 V. Leble, R. Lima, R. Dias et al., “Asymmetry of red blood cell motions in a microchannel with a diverging and converging bifurcation,” Biomicrofluidics, vol. 5, no. 4, Article ID 044120, 2011. View at: Publisher Site  Google Scholar
 G. Tomaiuolo, M. Barra, V. Preziosi, A. Cassinese, B. Rotoli, and S. Guido, “Microfluidics analysis of red blood cell membrane viscoelasticity,” Lab on a Chip—Miniaturisation for Chemistry and Biology, vol. 11, no. 3, pp. 449–454, 2011. View at: Publisher Site  Google Scholar
 Q. Guo, S. J. Reiling, P. Rohrbach, and H. Ma, “Microfluidic biomechanical assay for red blood cells parasitized by Plasmodium falciparum,” Lab on a Chip—Miniaturisation for Chemistry and Biology, vol. 12, no. 6, pp. 1143–1150, 2012. View at: Publisher Site  Google Scholar
 Y. Imai, H. Kondo, T. Ishikawa, C. Teck Lim, and T. Yamaguchi, “Modeling of hemodynamics arising from malaria infection,” Journal of Biomechanics, vol. 43, no. 7, pp. 1386–1393, 2010. View at: Publisher Site  Google Scholar
 Y. Imai, K. Nakaaki, H. Kondo, T. Ishikawa, C. Teck Lim, and T. Yamaguchi, “Margination of red blood cells infected by Plasmodium falciparum in a microvessel,” Journal of Biomechanics, vol. 44, no. 8, pp. 1553–1558, 2011. View at: Publisher Site  Google Scholar
 D. A. Fedosov, B. Caswell, S. Suresh, and G. E. Karniadakis, “Quantifying the biophysical characteristics of Plasmodiumfalciparumparasitized red blood cells in microcirculation,” Proceedings of the National Academy of Sciences of the United States of America, vol. 108, no. 1, pp. 35–39, 2011. View at: Publisher Site  Google Scholar
 T. Wu and J. J. Feng, “Simulation of malariainfected red blood cells in microfluidic channels: passage and blockage,” Biomicrofluidics, vol. 7, no. 4, Article ID 044115, 2013. View at: Publisher Site  Google Scholar
 T. Ye, N. PhanThien, B. Cheong Khoo, and C. Teck Lim, “Numerical modelling of a healthy/malariainfected erythrocyte in shear flow using dissipative particle dynamics method,” Journal of Applied Physics, vol. 115, no. 22, Article ID 224701, 2014. View at: Publisher Site  Google Scholar
 M. Navidbakhsh and M. Rezazadeh, “An immersed boundarylattice Boltzmann model for simulation of malariainfected red blood cell in microchannel,” Scientia Iranica, vol. 19, no. 5, pp. 1329–1336, 2012. View at: Publisher Site  Google Scholar
 M. Aingaran, R. Zhang, S. K. Law et al., “Host cell deformability is linked to transmission in the human malaria parasite Plasmodium falciparum,” Cellular Microbiology, vol. 14, no. 7, pp. 983–993, 2012. View at: Publisher Site  Google Scholar
 P. Bagchi, P. C. Johnson, and A. S. Popel, “Computational fluid dynamic simulation of aggregation of deformable cells in a shear flow,” Journal of Biomechanical Engineering, vol. 127, no. 7, pp. 1070–1080, 2005. View at: Publisher Site  Google Scholar
 Y. Liu and W. K. Liu, “Rheology of red blood cell aggregation by computer simulation,” Journal of Computational Physics, vol. 220, no. 1, pp. 139–154, 2006. View at: Publisher Site  Google Scholar
 T. Wang, T.W. Pan, Z. W. Xing, and R. Glowinski, “Numerical simulation of rheology of red blood cell rouleaux in microchannels,” Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, vol. 79, no. 4, Article ID 041916, 2009. View at: Publisher Site  Google Scholar
 T. Wang, Z. Xing, and D. Xing, “Structureinduced dynamics of erythrocyte aggregates by microscale simulation,” Journal of Applied Mathematics, vol. 2013, Article ID 409387, 13 pages, 2013. View at: Publisher Site  Google Scholar
 H. Li, T. Ye, and K. Y. Lam, “Computational analysis of dynamic interaction of two red blood cells in a capillary,” Cell Biochemistry and Biophysics, vol. 69, no. 3, pp. 673–680, 2014. View at: Publisher Site  Google Scholar
 P. Bagchi, “Mesoscale simulation of blood flow in small vessels,” Biophysical Journal, vol. 92, no. 6, pp. 1858–1877, 2007. View at: Publisher Site  Google Scholar
 K. I. Tsubota, S. Wada, and T. Yamaguchi, “Simulation study on effects of hematocrit on blood flow properties using particle method,” Journal of Biomechanical Science and Engineering, vol. 1, no. 1, pp. 159–170, 2006. View at: Publisher Site  Google Scholar
 H. Noguchi and G. Gompper, “Shape transitions of fluid vesicles and red blood cells in capillary flows,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 40, pp. 14159–14164, 2005. View at: Publisher Site  Google Scholar
 D. A. Fedosov, B. Caswell, and G. E. Karniadakis, “A multiscale red blood cell model with accurate mechanics, rheology, and dynamics,” Biophysical Journal, vol. 98, no. 10, pp. 2215–2225, 2010. View at: Publisher Site  Google Scholar
 D. A. Fedosov, B. Caswell, and G. E. Karniadakis, “Systematic coarsegraining of spectrinlevel red blood cell models,” Computer Methods in Applied Mechanics and Engineering, vol. 199, no. 2932, pp. 1937–1948, 2010. View at: Publisher Site  Google Scholar
 T. Wang, U. Rongin, and Z. Xing, “A microscale simulation of red blood cell passage through symmetric and asymmetric bifurcated vessels,” Scientific Reports, vol. 6, Article ID 20262, 2016. View at: Publisher Site  Google Scholar
 K. Tsukada, E. Sekizuka, C. Oshio, and H. Minamitani, “Direct measurement of erythrocyte deformability in diabetes mellitus with a transparent microchannel capillary model and highspeed video camera system,” Microvascular Research, vol. 61, no. 3, pp. 231–239, 2001. View at: Publisher Site  Google Scholar
 J. M. Higgins, D. T. Eddington, S. N. Bhatia, and L. Mahadevan, “Sickle cell vasoocclusion and rescue in a microfluidic device,” Proceedings of the National Academy of Sciences of the United States of America, vol. 104, no. 51, pp. 20496–20500, 2007. View at: Publisher Site  Google Scholar
 S. Chien and K.M. Jan, “Ultrastructural basis of the mechanism of rouleaux formation,” Microvascular Research, vol. 5, no. 2, pp. 155–166, 1973. View at: Publisher Site  Google Scholar
 B. Neu and H. J. Meiselman, “Depletionmediated red blood cell aggregation in polymer solutions,” Biophysical Journal, vol. 83, no. 5, pp. 2482–2490, 2002. View at: Publisher Site  Google Scholar
 R. Glowinski, T.W. Pan, and J. Periaux, “A fictitious domain method for Dirichlet problem and applications,” Computer Methods in Applied Mechanics and Engineering, vol. 111, no. 34, pp. 283–303, 1994. View at: Publisher Site  Google Scholar
 R. Glowinski, T.W. Pan, and J. Periaux, “A fictitious domain method for external incompressible viscous flow modeled by NavierStokes equations,” Computer Methods in Applied Mechanics and Engineering, vol. 112, no. 14, pp. 133–148, 1994. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 C. S. Peskin, “Numerical analysis of blood flow in the heart,” Journal of Computational Physics, vol. 25, no. 3, pp. 220–252, 1977. View at: Publisher Site  Google Scholar  MathSciNet
 A. Quarteroni, “Cardiovascular mathematics,” in Proceedings of the International Congress of Mathematicians, Madrid, Spain, August 2006. View at: Google Scholar
Copyright
Copyright © 2016 Tong Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.