skip to main content

Title: Interplay of deformability and adhesion on localization of elastic micro-particles in blood flow
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 more » 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. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
Journal of Fluid Mechanics
Page Range or eLocation-ID:
55 to 87
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 processmore »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.« less
  2. A flow vessel with an elastic wall can deform significantly due to viscous fluid flow within it, even at vanishing Reynolds number (no fluid inertia). Deformation leads to an enhancement of throughput due to the change in cross‐sectional area. The latter gives rise to a non‐constant pressure gradient in the flow‐wise direction and, hence, to a nonlinear flow rate–pressure drop relation (unlike the Hagen–Poiseuille law for a rigid tube). Many biofluids are non‐Newtonian, and are well approximated by generalized Newtonian (say, power‐law) rheological models. Consequently, we analyze the problem of steady low Reynolds number flow of a generalized Newtonian fluid through a slender elastic tube by coupling fluid lubrication theory to a structural problem posed in terms of Donnell shell theory. A perturbative approach (in the slenderness parameter) yields analytical solutions for both the flow and the deformation. Using matched asymptotics, we obtain a uniformly valid solution for the tube's radial displacement, which features both a boundary layer and a corner layer caused by localized bending near the clamped ends. In doing so, we obtain a “generalized Hagen–Poiseuille law” for soft microtubes. We benchmark the mathematical predictions against three‐dimensional two‐way coupled direct numerical simulations (DNS) of flow and deformation performedmore »using the commercial computational engineering platform by ANSYS. The simulations show good agreement and establish the range of validity of the theory. Finally, we discuss the implications of the theory on the problem of the flow‐induced deformation of a blood vessel, which is featured in some textbooks.« 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 marginationmore »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.« less
  4. The diverse chemical and physical reactions encountered during cooking connect us to science every day. Here, we theoretically and experimentally investigate the swelling and softening of pasta due to liquid imbibition as well as the elastic deformation and adhesion of pasta due to capillary force. As water diffuses into the pasta during cooking, it softens gradually from the outside inward as starch swells and relaxes. The softening follows three sequential regimes: regime I, the hard-glassy region, shows a slow decrease in modulus with cooking time; regime II, the glassy to rubbery transition region, or leathery region, is characterized by a very fast, several orders of magnitude drop in elastic modulus and regime III, the rubbery region, has an asymptotic modulus about four orders of magnitude lower than the raw pasta. We present experiments and theory to capture these regimes and relate them to the heterogeneous microstructure changes associated with swelling. Interestingly, we observe a modulus drop of two orders of magnitude within the range of “al dente” cooking duration, and we find the modulus to be extremely sensitive to the amount of salt added to the boiling water. While most chefs can gauge the pasta by tasting its texture, ourmore »proposed experiment, which only requires a measurement with a ruler, can precisely provide an optimal cooking time finely tuned for various kinds of pasta shapes.

    « less
  5. Introduction: Vaso-occlusive crises (VOCs) are a leading cause of morbidity and early mortality in individuals with sickle cell disease (SCD). These crises are triggered by sickle red blood cell (sRBC) aggregation in blood vessels and are influenced by factors such as enhanced sRBC and white blood cell (WBC) adhesion to inflamed endothelium. Advances in microfluidic biomarker assays (i.e., SCD Biochip systems) have led to clinical studies of blood cell adhesion onto endothelial proteins, including, fibronectin, laminin, P-selectin, ICAM-1, functionalized in microchannels. These microfluidic assays allow mimicking the physiological aspects of human microvasculature and help characterize biomechanical properties of adhered sRBCs under flow. However, analysis of the microfluidic biomarker assay data has so far relied on manual cell counting and exhaustive visual morphological characterization of cells by trained personnel. Integrating deep learning algorithms with microscopic imaging of adhesion protein functionalized microfluidic channels can accelerate and standardize accurate classification of blood cells in microfluidic biomarker assays. Here we present a deep learning approach into a general-purpose analytical tool covering a wide range of conditions: channels functionalized with different proteins (laminin or P-selectin), with varying degrees of adhesion by both sRBCs and WBCs, and in both normoxic and hypoxic environments. Methods: Our neuralmore »networks were trained on a repository of manually labeled SCD Biochip microfluidic biomarker assay whole channel images. Each channel contained adhered cells pertaining to clinical whole blood under constant shear stress of 0.1 Pa, mimicking physiological levels in post-capillary venules. The machine learning (ML) framework consists of two phases: Phase I segments pixels belonging to blood cells adhered to the microfluidic channel surface, while Phase II associates pixel clusters with specific cell types (sRBCs or WBCs). Phase I is implemented through an ensemble of seven generative fully convolutional neural networks, and Phase II is an ensemble of five neural networks based on a Resnet50 backbone. Each pixel cluster is given a probability of belonging to one of three classes: adhered sRBC, adhered WBC, or non-adhered / other. Results and Discussion: We applied our trained ML framework to 107 novel whole channel images not used during training and compared the results against counts from human experts. As seen in Fig. 1A, there was excellent agreement in counts across all protein and cell types investigated: sRBCs adhered to laminin, sRBCs adhered to P-selectin, and WBCs adhered to P-selectin. Not only was the approach able to handle surfaces functionalized with different proteins, but it also performed well for high cell density images (up to 5000 cells per image) in both normoxic and hypoxic conditions (Fig. 1B). The average uncertainty for the ML counts, obtained from accuracy metrics on the test dataset, was 3%. This uncertainty is a significant improvement on the 20% average uncertainty of the human counts, estimated from the variance in repeated manual analyses of the images. Moreover, manual classification of each image may take up to 2 hours, versus about 6 minutes per image for the ML analysis. Thus, ML provides greater consistency in the classification at a fraction of the processing time. To assess which features the network used to distinguish adhered cells, we generated class activation maps (Fig. 1C-E). These heat maps indicate the regions of focus for the algorithm in making each classification decision. Intriguingly, the highlighted features were similar to those used by human experts: the dimple in partially sickled RBCs, the sharp endpoints for highly sickled RBCs, and the uniform curvature of the WBCs. Overall the robust performance of the ML approach in our study sets the stage for generalizing it to other endothelial proteins and experimental conditions, a first step toward a universal microfluidic ML framework targeting blood disorders. Such a framework would not only be able to integrate advanced biophysical characterization into fast, point-of-care diagnostic devices, but also provide a standardized and reliable way of monitoring patients undergoing targeted therapies and curative interventions, including, stem cell and gene-based therapies for SCD. Disclosures Gurkan: Dx Now Inc.: Patents & Royalties; Xatek Inc.: Patents & Royalties; BioChip Labs: Patents & Royalties; Hemex Health, Inc.: Consultancy, Current Employment, Patents & Royalties, Research Funding.« less