Numerical Investigation of the Effects of Red Blood Cell Cytoplasmic Viscosity Contrasts on Single Cell and Bulk Transport Behaviour

Numerical simulation of cellular blood flow in curved micro-vessels with saccular aneurysms: Effect of curvature degree and hematocrit level

The dynamics of cellular blood flow in curved vessels considerably differ from those in straight vessels. It is reported that clotting development is significantly affected by vessel shape irregularities. Thus, the current study aims to investigate the effect of curvature degree and hematocrit level on cellular blood flow in a curved micro-vessel with a saccular aneurysm. Accordingly, a three-dimensional numerical simulation is performed using a validated code developed for cellular blood flow problems. The obtained results show that the cell-free layer thickness is highly dependent on the curvature degree and hematocrit level, which may have a remarkable impact on the apparent viscosity of blood as well as the dynamics of other particles such as drug particulates. The near-wall region exhibits the highest degree of cell deformation, whereas the red blood cells within the aneurysm zone remain nearly undeformed. Meanwhile, the velocity of the red blood cells decreases with the increase in curvature degree, which can affect the quality of the oxygenation process. Because of the saccular aneurysm, a considerable decrease in plasma velocity is predicted. Moreover, no secondary flows are detected in the curved vessel except in the aneurysm zone. An increase in the curvature degree is expected to reduce the blood flow rate by about 10%. Furthermore, low wall shear stress values are predicted in the straight case compared to the values at the apex of the curved vessel, which may affect the structure and function of the endothelial cells of the vessel wall and, hence, increase the aneurysm rupture possibility.

Cellular Blood Flow Modeling with HemoCell

Many of the intriguing properties of blood originate from its cellular nature. Bulk effects, such as viscosity, depend on the local shear rates and on the size of the vessels. While empirical descriptions of bulk rheology are available for decades, their validity is limited to the experimental conditions they were observed under. These are typically artificial scenarios (e.g., perfectly straight glass tube or in pure shear with no gradients). Such conditions make experimental measurements simpler; however, they do not exist in real systems (i.e., in a real human circulatory system). Therefore, as we strive to increase our understanding on the cardiovascular system and improve the accuracy of our computational predictions, we need to incorporate a more comprehensive description of the cellular nature of blood. This, however, presents several computational challenges that can only be addressed by high performance computing. In this chapter, we describe HemoCell ( https://www.hemocell.eu ), an open-source high-performance cellular blood flow simulation, which implements validated mechanical models for red blood cells and is capable of reproducing the emergent transport characteristics of such a complex cellular system. We discuss the accuracy and the range of validity, and demonstrate applications on a series of human diseases.

Lingering Dynamics of Type 2 Diabetes Mellitus Red Blood Cells in Retinal Arteriolar Bifurcations

It has been proven that the deformability of red blood cells (RBC) is reduced owing to changes in mechanical properties, such as diabetes mellitus and hypertension. To probe the effects of RBC morphological and physical parameters on the flow field in bifurcated arterioles, three types of RBC models with various degrees of biconcave shapes were built based on the in vitro experimental data. The dynamic behaviors of the RBCs in shear flow were simulated to validate the feasibility of the finite element-Arbitrary Lagrangian-Eulerian method with a moving mesh. The influences of the shear rate and viscosity ratios on RBC motions were investigated. The motion of RBCs in arteriolar bifurcations was further simulated. Abnormal variations in the morphological and physical parameters of RBCs may lead to diminished tank-tread motion and enhanced tumbling motion in shear flow. Moreover, abnormal RBC variations can result in slower RBC motion at the bifurcation with a longer transmit time and greater flow resistance, which may further cause inadequate local oxygen supply. These findings would provide useful insights into the microvascular complications in diabetes mellitus.

The effect of stiffened diabetic red blood cells on wall shear stress in a reconstructed 3D microaneurysm

Blood flow within the vasculature of the retina has been found to influence the progression of diabetic retinopathy. In this research cell resolved blood flow simulations are used to study the pulsatile flow of whole blood through a segmented retinal microaneurysm. Images were collected using adaptive optics optical coherence tomography of the retina of a patient with diabetic retinopathy, and a sidewall (sacciform) microaneurysm was segmented from the volumetric data. The original microaneurysm neck width was varied to produce two additional aneurysm geometries in order to probe the influence of neck width on the transport of red blood cells and platelets into the aneurysm. Red blood cell membrane stiffness was also increased to resolve the impact of rigid red blood cells, as a result of diabetes, in blood flow. Wall shear stress and wall shear stress gradients were calculated throughout the aneurysm domains, and the quantification of the influence of the red blood cells is presented. Average wall shear stress and wall shear stress gradients increased due to the increase of red blood cell membrane stiffness. Stiffened red blood cells were also found to induce higher local wall shear stress and wall shear stress gradients as they passed through the leading and draining parental vessels. Stiffened red blood cells were found to penetrate the aneurysm sac more than healthy red blood cells, as well as decreasing the margination of platelets to the vessel walls of the parental vessel, which caused a decrease in platelet penetration into the aneurysm sac.

Modeling Red Blood Cell Viscosity Contrast Using Inner Soft Particle Suspension

The inner viscosity of a biological red blood cell is about five times larger than the viscosity of the blood plasma. In this work, we use dissipative particles to enable the proper viscosity contrast in a mesh-based red blood cell model. Each soft particle represents a coarse-grained virtual cluster of hemoglobin proteins contained in the cytosol of the red blood cell. The particle interactions are governed by conservative and dissipative forces. The conservative forces have purely repulsive character, whereas the dissipative forces depend on the relative velocity between the particles. We design two computational experiments that mimic the classical viscometers. With these experiments we study the effects of particle suspension parameters on the inner cell viscosity and provide parameter sets that result in the correct viscosity contrast. The results are validated with both static and dynamic biological experiment, showing an improvement in the accuracy of the original model without major increase in computational complexity.

Effect of cytosol viscosity on the flow behavior of red blood cell suspensions in microvessels

OBJECTIVE

The flow behavior of blood is strongly affected by red blood cell (RBC) properties, such as the viscosity ratio C between cytosol and suspending medium, which can significantly be altered in several pathologies (e.g. sickle-cell disease, malaria). The main objective of this study is to understand the effect of C on macroscopic blood flow properties such as flow resistance in microvessels, and to link it to the deformation and dynamics of single RBCs.

METHODS

We employ mesoscopic hydrodynamic simulations to investigate flow properties of RBC suspensions with different cytosol viscosities for various flow conditions in cylindrical microchannels.

RESULTS

Starting from a dispersed cell configuration which approximates RBC dispersion at vessel bifurcations in the microvasculature, we find that the flow convergence and development of RBC-free layer (RBC-FL) depend only weakly on C, and require a convergence length in the range of 25D-50D, where D is channel diameter. In vessels with D ≤ 20 μ m , the final resistance of developed flow is nearly the same for C = 5 and C = 1, while for D = 40 μ m , the flow resistance for C = 5 is about 10% larger than for C = 1. The similarities and differences in flow resistance can be explained by viscosity-dependent RBC-FL thicknesses, which are associated with the viscosity-dependent dynamics of single RBCs.

CONCLUSIONS

The weak effect on the flow resistance and RBC-FL explains why RBCs can contain a high concentration of hemoglobin for efficient oxygen delivery, without a pronounced increase in the flow resistance. Furthermore, our results suggest that significant alterations in microvascular flow in various pathologies are likely not due to mere changes in cytosolic viscosity.

The influence of red blood cell deformability on hematocrit profiles and platelet margination

The influence of red blood cell (RBC) deformability in whole blood on platelet margination is investigated using confocal microscopy measurements of flowing human blood and cell resolved blood flow simulations. Fluorescent platelet concentrations at the wall of a glass chamber are measured using confocal microscopy with flowing human blood containing varying healthy-to-stiff RBC fractions. A decrease is observed in the fluorescent platelet signal at the wall due to the increase of stiffened RBCs in flow, suggesting a decrease of platelet margination due to an increased fraction of stiffened RBCs present in the flow. In order to resolve the influence of stiffened RBCs on platelet concentration at the channel wall, cell-pair and bulk flow simulations are performed. For homogeneous collisions between RBC pairs, a decrease in final displacement after a collision with increasing membrane stiffness is observed. In heterogeneous collisions between healthy and stiff RBC pairs, it is found that the stiffened RBC is displaced most. The influence of RBC deformability on collisions between RBCs and platelets was found to be negligible due to their size and mass difference. For a straight vessel geometry with varying healthy-to-stiff RBC ratios, a decrease was observed in the red blood cell-free layer and platelet margination due to an increase in stiffened RBCs present in flow.

Simulation of circulating tumor cell transport and adhesion in cell suspensions in microfluidic devices

Understanding cell transport and adhesion dynamics under flow is important for many biotransport problems. We investigated the influence of cell size, ligand coating density, micropost size, and intercellular collisions on circulating tumor cell adhesion and transport in microfluidic devices. The cells were modeled as coarse-grained cell membranes and the adhesion was modeled as pairwise interacting potentials, while the fluid was solved using the lattice Boltzmann method. The coupling between the cell and the fluid was achieved through the immersed boundary method. The cell showed transient rolling adhesion in high shear regions and firm adhesion in low shear regions. The adhesive force for rolling cells on a micropost was increasing before the cell reached the crest of the post and then decreasing afterward. The adhesive strength for cells increases with ligand coating density. Cell trajectories in a microfluidic device with a shifted post design were studied as well. At low concentrations, the majority of the cells follow streamlines closely. However, the intercellular collision and collision from red blood cells impacted the cell trajectories. An L 2 norm of | e | was defined to characterize the difference between the cell trajectories and the associated streamlines. It was shown that | e | L 2 increases with micropost sizes and cell concentrations.