skip to main content


Title: Shear rate dependent margination of sphere-like, oblate-like and prolate-like micro-particles within blood flow
This study investigates the shear rate dependent margination of micro-particles (MPs) with different shapes in blood flow through numerical simulations. We develop a multiscale computational model to handle the fluid–structure interactions involved in the blood flow simulations. The lattice Boltzmann method (LBM) is used to solve the plasma dynamics and a coarse-grained model is employed to capture the dynamics of red blood cells (RBCs) and MPs. These two solvers are coupled together by the immersed boundary method (IBM). The shear rate dependent margination of sphere MPs is firstly investigated. We find that margination of sphere MPs dramatically increases with the increment of wall shear rate  ω under 800 s −1 , induced by the breaking of rouleaux in blood flow. However, the margination probability only slowly grows when  ω > 800 s −1 . Furthermore, the shape effect of MPs is examined by comparing the margination behaviors of sphere-like, oblate-like and prolate-like MPs under different wall shear rates. We find that the margination of MPs is governed by the interplay of two factors: hydrodynamic collisions with RBCs including the collision frequency and collision displacement of MPs, and near wall dynamics. MPs that demonstrate poor performance in one process such as collision frequency may stand out in the other process like near wall dynamics. Specifically, the ellipsoidal MPs (oblate and prolate) with small aspect ratio (AR) outperform those with large AR regardless of the wall shear rate, due to their better performance in both the collision with RBCs and near wall dynamics. Additionally, we find there exists a transition shear rate region 700 s −1 <  ω < 900 s −1 for all of these MPs: the margination probability dramatically increases with the shear rate below this region and slowly grows above this region, similar to sphere MPs. We further use the surface area to volume ratio (SVR) to distinguish different shaped MPs and illustrate their shear rate dependent margination in a contour in the shear rate–SVR plane. It is of significance that we can approximately predict the margination of MPs with a specific SVR. All these simulation results can be potentially applied to guide the design of micro-drug carriers for biomedical applications.  more » « less
Award ID(s):
1755779
PAR ID:
10085628
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Soft Matter
Volume:
14
Issue:
36
ISSN:
1744-683X
Page Range / eLocation ID:
7401 to 7419
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The margination and adhesion of micro-particles (MPs) have been extensively investigated separately, due to their important applications in the biomedical field. However, the cascade process from margination to adhesion should play an important role in the transport of MPs in blood flow. To the best of our knowledge, this has not been explored in the past. Here we numerically study the margination behaviour of elastic MPs to blood vessel walls under the interplay of their deformability and adhesion to the vessel wall. We use the lattice Boltzmann method and molecular dynamics to solve the fluid dynamics and particle dynamics (including red blood cells (RBCs) and elastic MPs) in blood flow, respectively. Additionally, a stochastic ligand–receptor binding model is employed to capture the adhesion behaviours of elastic MPs on the vessel wall. Margination probability is used to quantify the localization of elastic MPs at the wall. Two dimensionless numbers are considered to govern the whole process: the capillary number $Ca$ , denoting the ratio of viscous force of fluid flow to elastic interfacial force of MP, and the adhesion number $Ad$ , representing the ratio of adhesion strength to viscous force of fluid flow. We systematically vary them numerically and a margination probability contour is obtained. We find that there exist two optimal regimes favouring high margination probability on the plane $Ca{-}Ad$ . The first regime, namely region I, is that with high adhesion strength and moderate particle stiffness; the other one, region II, has moderate adhesion strength and large particle stiffness. We conclude that the existence of optimal regimes is governed by the interplay of particle deformability and adhesion strength. The corresponding underlying mechanism is also discussed in detail. There are three major factors that contribute to the localization of MPs: (i) near-wall hydrodynamic collision between RBCs and MPs; (ii) deformation-induced migration due to the presence of the wall; and (iii) adhesive interaction between MPs and the wall. Mechanisms (i) and (iii) promote margination, while (ii) hampers margination. These three factors perform different roles and compete against each other when MPs are located in different regions of the flow channel, i.e. near-wall region. In optimal region I, adhesion outperforms deformation-induced migration; and in region II, the deformation-induced migration is small compared to the coupling of near-wall hydrodynamic collision and adhesion. The finding of optimal regimes can help the understanding of localization of elastic MPs at the wall under the adhesion effect in blood flow. More importantly, our results suggest that softer MP or stronger adhesion is not always the best choice for the localization of MPs. 
    more » « less
  2. The motion of cells orthogonal to the direction of main flow is of importance in natural and engineered systems. The lateral movement of red blood cells (RBCs) distal to sudden expansion is considered to influence the formation and progression of thrombosis in venous valves, aortic aneurysms, and blood-circulating devices and is also a determining parameter for cell separation applications in flow-focusing microfluidic devices. Although it is known that the unique geometry of venous valves alters the blood flow patterns and cell distribution in venous valve sinuses, the interactions between fluid flow and RBCs have not been elucidated. Here, using a dilute cell suspension in an in vitro microfluidic model of a venous valve, we quantified the spatial distribution of RBCs by microscopy and image analysis, and using micro-particle image velocimetry and 3D computational fluid dynamics simulations, we analyzed the complex flow patterns. The results show that the local hematocrit in the valve pockets is spatially heterogeneous and is significantly different from the feed hematocrit. Above a threshold shear rate, the inertial separation of streamlines and lift forces contribute to an uneven distribution of RBCs in the vortices, the entrapment of RBCs in the vortices, and non-monotonic wall shear stresses in the valve pockets. Our experimental and computational characterization provides insights into the complex interactions between fluid flow, RBC distribution, and wall shear rates in venous valve mimics, which is of relevance to understanding the pathophysiology of thrombosis and improving cell separation efficiency.

     
    more » « less
  3. 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
  4. Through high-fidelity numerical simulation based on the lattice Boltzmann method, we have conducted an in-depth study on the heat and mass transport from an oblate spheroid neutrally suspended in a simple shear flow. In the simulation, the temperature and mass concentration are modeled as a passive scalar released at the surface of the spheroid. The fluid dynamics induced by the interaction between the carrier fluid and the suspended spheroid, as well as the resultant scalar transport process, have been extensively investigated. A coupled transport mechanism comprising several components of the flow around the oblate spheroid has been identified. The effects of the Reynolds number and the aspect ratio of the spheroid on the flow characteristics and scalar transport rate are examined. The variation of the nondimensional scalar transport rate suggests that the effect of spheroid shape on scalar transfer rate can be decoupled from the effects of Peclet and Reynolds numbers, which facilitates the development of a correlation of scalar transfer rate for oblate spheroids based on the well-developed correlations for a sphere.

     
    more » « less
  5. Abstract

    The wall shear stress (WSS) exerted by blood flowing through microvascular capillaries is an established driver of new blood vessel growth, or angiogenesis. Such adaptations are central to many physiological processes in both health and disease, yet three-dimensional (3D) WSS characteristics in real angiogenic microvascular networks are largely unknown. This marks a major knowledge gap because angiogenesis, naturally, is a 3D process. To advance current understanding, we model 3D red blood cells (RBCs) flowing through rat angiogenic microvascular networks using state-of-the-art simulation. The high-resolution fluid dynamics reveal 3D WSS patterns occurring at sub-endothelial cell (EC) scales that derive from distinct angiogenic morphologies, including microvascular loops and vessel tortuosity. We identify the existence of WSS hot and cold spots caused by angiogenic surface shapes and RBCs, and notably enhancement of low WSS regions by RBCs. Spatiotemporal characteristics further reveal how fluctuations follow timescales of RBC “footprints.” Altogether, this work provides a new conceptual framework for understanding how shear stress might regulate EC dynamics in vivo.

     
    more » « less