Hydrodynamic interactions generate a diffusive motion in particulates in a shear flow, which plays seminal roles in overall particulate rheology and its microstructure. Here we investigate the shear induced diffusion in a red-blood cell (RBC) suspension using a numerical simulation resolving individual motion and deformation of RBCs. The non-spherical resting shape of RBCs gives rise to qualitatively different regimes of cell dynamics in a shear flow such as tank-treading, breathing, tumbling and swinging, depending on the cell flexibility determined by the elastic capillary number. We show that the transition from tumbling to tank-treading causes a reduction in the gradient diffusivity. The diffusivity is computed using a continuum approach from the evolution of a randomly packed cell-layer width with time as well as by the dynamic structure factor of the suspension. Both approaches, although operationally different, match and show that for intermediate capillary numbers RBCs cease tumbling accompanied by a drop in the coefficient of gradient diffusivity. A further increase of capillary number increases the diffusivity due to increased deformation. The effects of bending modulus and viscosity ratio variations are also briefly investigated. The computed shear induced diffusivity was compared with values in the literature. Apart from its effects in margination of cells in blood flow and use in medical diagnostics, the phenomenon broadly offers important insights into suspensions of deformable particles with non-spherical equilibrium shapes, which also could play a critical role in using particle flexibility for applications such as label free separation or material processing.
more »
« less
The effect of rigid cells on blood viscosity: linking rheology and sickle cell anemia
Sickle cell anemia (SCA) is a disease that affects red blood cells (RBCs). Healthy RBCs are highly deformable objects that under flow can penetrate blood capillaries smaller than their typical size. In SCA there is an impaired deformability of some cells, which are much stiffer and with a different shape than healthy cells, and thereby affect regular blood flow. It is known that blood from patients with SCA has a higher viscosity than normal blood. However, it is unclear how the rigidity of cells is related to the viscosity of blood, in part because SCA patients are often treated with transfusions of variable amounts of normal RBCs and only a fraction of cells will be stiff. Here, we report systematic experimental measurements of the viscosity of a suspension varying the fraction of rigid particles within a suspension of healthy cells. We also perform systematic numerical simulations of a similar mixed suspension of soft RBCs, rigid particles, and their hydrodynamic interactions. Our results show that there is a rheological signature within blood viscosity to clearly identify the fraction of rigidified cells among healthy deformable cells down to a 5% volume fraction of rigidified cells. Although aggregation of RBCs is known to affect blood rheology at low shear rates, and our simulations mimic this effect via an adhesion potential, we show that such adhesion, or aggregation, is unlikely to provide a physical rationalization for the viscosity increase observed in the experiments at moderate shear rates due to rigidified cells. Through numerical simulations, we also highlight that most of the viscosity increase of the suspension is due to the rigidity of the particles rather than their sickled or spherical shape. Our results are relevant to better characterize SCA, provide useful insights relevant to rheological consequences of blood transfusions, and, more generally, extend to the rheology of mixed suspensions having particles with different rigidities, as well as offering possibilities for developments in the field of soft material composites.
more »
« less
- Award ID(s):
- 1951526
- PAR ID:
- 10340743
- Date Published:
- Journal Name:
- Soft Matter
- Volume:
- 18
- Issue:
- 3
- ISSN:
- 1744-683X
- Page Range / eLocation ID:
- 554 to 565
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
AbstractSickle cell disease (SCD) is canonically characterized by reduced red blood cell (RBC) deformability, leading to microvascular obstruction and inflammation. Although the biophysical properties of sickle RBCs are known to influence SCD vasculopathy, the contribution of poor RBC deformability to endothelial dysfunction has yet to be fully explored. Leveraging interrelated in vitro and in silico approaches, we introduce a new paradigm of SCD vasculopathy in which poorly deformable sickle RBCs directly cause endothelial dysfunction via mechanotransduction, during which endothelial cells sense and pathophysiologically respond to aberrant physical forces independently of microvascular obstruction, adhesion, or hemolysis. We demonstrate that perfusion of sickle RBCs or pharmacologically-dehydrated healthy RBCs into small venule-sized “endothelialized” microfluidics leads to pathologic physical interactions with endothelial cells that directly induce inflammatory pathways. Using a combination of computational simulations and large venule-sized endothelialized microfluidics, we observed that perfusion of heterogeneous sickle RBC subpopulations with varying deformability, as well as suspensions of dehydrated normal RBCs admixed with normal RBCs, leads to aberrant margination of the less-deformable RBC subpopulations toward the vessel walls, causing localized, increased shear stress. Increased wall stress is dependent on the degree of subpopulation heterogeneity and oxygen tension and leads to inflammatory endothelial gene expression via mechanotransductive pathways. Our multifaceted approach demonstrates that the presence of sickle RBCs with reduced deformability leads directly to pathological physical (ie, direct collisions and/or compressive forces) and shear-mediated interactions with endothelial cells and induces an inflammatory response, thereby elucidating the ubiquity of vascular dysfunction in SCD.more » « less
-
Computational modeling and simulations can tackle a broad range of morphological, mechanical, and rheological problems relevant to blood and blood cells. Here, we review some continuum-based and particle-based computational approaches towards the modeling of healthy and diseased red blood cells (RBCs) with focus on the most recent contributions, including the three-level multiscale RBC model coupled with the boundary integral method of surrounding flows and two-component RBC models with explicit descriptions of lipid bilayer, cytoskeleton, and transmembrane proteins.more » « less
-
Networks and dense suspensions frequently reside near a boundary between soft (or fluidlike) and rigid (or solidlike) regimes. Transitions between these regimes can be driven by changes in structure, density, or applied stress or strain. In general, near the onset or loss of rigidity in these systems, dissipation-limiting heterogeneous nonaffine rearrangements dominate the macroscopic viscoelastic response, giving rise to diverging relaxation times and power-law rheology. Here, we describe a simple quantitative relationship between nonaffinity and the excess viscosity. We test this nonaffinity-viscosity relationship computationally and demonstrate its rheological consequences in simulations of strained filament networks and dense suspensions. We also predict critical signatures in the rheology of semiflexible and stiff biopolymer networks near the strain stiffening transition.more » « less
-
Discrete-particle simulations of bidisperse shear thickening suspensions are reported. The work considers two packing parameters, the large-to-small particle radius ratio ranging from [Formula: see text] (nearly monodisperse) to [Formula: see text], and the large particle fraction of the total solid loading with values [Formula: see text], 0.5, and 0.85. Particle-scale simulations are performed over a broad range of shear stresses using a simulation model for spherical particles accounting for short-range lubrication forces, frictional interaction, and repulsion between particles. The variation of rheological properties and the maximum packing fraction [Formula: see text] with shear stress [Formula: see text] are reported. At a fixed volume fraction [Formula: see text], bidispersity decreases the suspension relative viscosity [Formula: see text], where [Formula: see text] is the suspension viscosity and [Formula: see text] is the suspending fluid viscosity, over the entire range of shear stresses studied. However, under low shear stress conditions, the suspension exhibits an unusual rheological behavior: the minimum viscosity does not occur as expected at [Formula: see text], but instead decreases with further increase of [Formula: see text] to [Formula: see text]. The second normal stress difference [Formula: see text] acts similarly. This behavior is caused by particles ordering into a layered structure, as is also reflected by the zero slope with respect to time of the mean-square displacement in the velocity gradient direction. The relative viscosity [Formula: see text] of bidisperse rate-dependent suspensions can be predicted by a power law linking it to [Formula: see text], [Formula: see text] in both low and high shear stress regimes. The agreement between the power law and experimental data from literature demonstrates that the model captures well the effect of particle size distribution, showing that viscosity roughly collapses onto a single master curve when plotted against the reduced volume fraction [Formula: see text].more » « less