Blood flow and cell-free layer in microvessels

Comment on "Dynamics and rheology of vesicles under confined Poiseuille flow" by Z. Gou, H. Zhang, A. Nait-Ouhra, M. Abbasi, A. Farutin and C. Misbah, Soft Matter, 2023, , 9101

In a recent paper, [Gou et al., Soft Matter, 2023, 19, 9101-9114] studied numerically the viscosity of a confined suspension of vesicles flowing in a channel as a function of vesicle concentration. In order to discuss the genericity of the observed behaviour, namely a nearly constant effective viscosity at low concentrations, we complement their study by a comparison with the few existing ones in the literature. In particular, we highlight that they fail to reproduce well established results for blood viscosity in microcirculation, thereby suggesting that the conclusions regarding the optimization of cell transport and oxygenation may not apply. We conclude with a quick discussion on potential improvements regarding numerical modeling, as long as physiological relevance is sought.

Hemodynamic Characteristics of a Tortuous Microvessel Using High-Fidelity Red Blood Cell Resolved Simulations

OBJECTIVE

Tortuous microvessels are characteristic of microvascular remodeling associated with numerous physiological and pathological scenarios. Three-dimensional (3D) hemodynamics in tortuous microvessels influenced by red blood cells (RBCs), however, are largely unknown, and important questions remain. Is blood viscosity influenced by vessel tortuosity? How do RBC dynamics affect wall shear stress (WSS) patterns and the near-wall cell-free layer (CFL) over a range of conditions? The objective of this work was to parameterize hemodynamic characteristics unique to a tortuous microvessel.

METHODS

RBC-resolved simulations were performed using an immersed boundary method-based 3D fluid dynamics solver. A representative tortuous microvessel was selected from a stimulated angiogenic network obtained from imaging of the rat mesentery and digitally reconstructed for the simulations. The representative microvessel was a venule with a diameter of approximately 20 μm. The model assumes a constant diameter along the vessel length and does not consider variations due to endothelial cell shapes or the endothelial surface layer.

RESULTS

Microvessel tortuosity was observed to increase blood apparent viscosity compared to a straight tube by up to 26%. WSS spatial variations in high curvature regions reached 23.6 dyne/cm2 over the vessel cross-section. The magnitudes of WSS and CFL thickness variations due to tortuosity were strongly influenced by shear rate and negligibly influenced by tube hematocrit levels.

CONCLUSIONS

New findings from this work reveal unique tortuosity-dependent hemodynamic characteristics over a range of conditions. The results provide new thought-provoking information to better understand the contribution of tortuous vessels in physiological and pathological processes and help improve reduced-order models.

Morphology, repulsion, and ordering of red blood cells in viscoelastic flows under confinement

Red blood cells (RBC), the primary carriers of oxygen in the body, play a crucial role across several biomedical applications, while also being an essential model system of a deformable object in the microfluidics and soft matter fields. However, RBC behavior in viscoelastic liquids, which holds promise in enhancing microfluidic diagnostic applications, remains poorly studied. We here show that using viscoelastic polymer solutions as a suspending carrier causes changes in the clustering and shape of flowing RBC in microfluidic flows when compared to a standard Newtonian suspending liquid. Additionally, when the local RBC concentration increases to a point where hydrodynamic interactions take place, we observe the formation of equally-spaced RBC structures, resembling the viscoelasticity-driven ordered particles observed previously in the literature, thus providing the first experimental evidence of viscoelasticity-driven cell ordering. The observed RBC ordering, unaffected by polymer molecular architecture, persists as long as the surrounding medium exhibits shear-thinning, viscoelastic properties. Complementary numerical simulations reveal that viscoelasticity-induced repulsion between RBCs leads to equidistant structures, with shear-thinning modulating this effect. Our results open the way for the development of new biomedical technologies based on the use of viscoelastic liquids while also clarifying fundamental aspects related to multibody hydrodynamic interactions in viscoelastic microfluidic flows.

Sequential Dual Alignments Introduce Synergistic Effect on Hexagonal Boron Nitride Platelets for Superior Thermal Performance

Planarly aligning 2D platelets is challenging due to their additional orientational freedom compared to 1D materials. This study reports a sequential dual-alignment approach, employing an extrusion-printing-induced shear force and rotating-magnetic-field-induced force couple for platelet planarly alignment in a yield-stress support bath. It is hypothesized that the partial alignment induced by a directional shear force facilitates subsequent axial rotation of the platelets for planar alignment under an external force couple, resulting in a synergistic alignment effect. This sequential dual-alignment approach achieves better planar alignment of 2D modified hexagonal boron nitride (mhBN). Specifically, the thermal conductivity of the 40 wt% mhBN/epoxy composite is significantly higher (692%) than that of unaligned composites, surpassing the cumulative effect of individual methods (only 133%) with a 5 times more synergistic effect. For 30, 40, and 50 wt% mhBN composites, the thermal conductivity values (5.9, 9.5, and 13.8 W m-1 K-1) show considerable improvement compared to the previously reported highest values (5.3, 6.6, and 8.6 W m-1 K-1). Additionally, a 3D mhBN/epoxy heat sink is printed and evaluated to demonstrate the feasibility of device fabrication. The approach enables the planar alignment of electrically or thermally conducting 2D fillers during 3D fabrication.

Numerical analysis of blood flow through stenosed microvessels using a multi-phase model

Blood flow in arterioles have attracted considerable research attention due to their clinical implications. However, the fluid structure interaction between red blood cells and plasma in the blood poses formidable difficulty to the computational efforts. In this contribution, we seek to represent the red blood cells in the blood as a continuous non-Newtonian phase and construct a multi-phase model for the blood flow in microvessels. The methods are presented and validated using a channel with sudden expansion. And the resulting blood flow inside a stenosed microvessel is investigated at different inlet velocity amplitudes and hematocrits. It is show that the increase of both inlet velocity amplitude and inlet hematocrit leads to longer and thicker cell-rich layer downstream the stenosis. Besides, it is found that the maximum values of wall shear stress scales up with inlet velocity amplitudes and hematocrits. These results show the validity of the proposed computational model and provide helpful insights into blood flow behaviors inside stenosed vessels.

Optical Coherence Tomography Angiography as a Diagnostic Tool for Diabetic Retinopathy

Diabetic retinopathy (DR) is the most common microvascular complication of diabetes mellitus, leading to visual impairment if left untreated. This review discusses the use of optical coherence tomography angiography (OCTA) as a diagnostic tool for the early detection and management of DR. OCTA is a fast, non-invasive, non-contact test that enables the detailed visualisation of the macular microvasculature in different plexuses. OCTA offers several advantages over fundus fluorescein angiography (FFA), notably offering quantitative data. OCTA is not without limitations, including the requirement for careful interpretation of artefacts and the limited region of interest that can be captured currently. We explore how OCTA has been instrumental in detecting early microvascular changes that precede clinical signs of DR. We also discuss the application of OCTA in the diagnosis and management of various stages of DR, including non-proliferative diabetic retinopathy (NPDR), proliferative diabetic retinopathy (PDR), diabetic macular oedema (DMO), diabetic macular ischaemia (DMI), and pre-diabetes. Finally, we discuss the future role of OCTA and how it may be used to enhance the clinical outcomes of DR.

Predicting capillary vessel network hemodynamics in silico by machine learning

Blood velocity and red blood cell (RBC) distribution profiles in a capillary vessel cross-section in the microcirculation are generally complex and do not follow Poiseuille's parabolic or uniform pattern. Existing imaging techniques used to map large microvascular networks in vivo do not allow a direct measurement of full 3D velocity and RBC concentration profiles, although such information is needed for accurate evaluation of the physiological variables, such as the wall shear stress (WSS) and near-wall cell-free layer (CFL), that play critical roles in blood flow regulation, disease progression, angiogenesis, and hemostasis. Theoretical network flow models, often used for hemodynamic predictions in experimentally acquired images of the microvascular network, cannot provide the full 3D profiles either. In contrast, such information can be readily obtained from high-fidelity computational models that treat blood as a suspension of deformable RBCs. These models, however, are computationally expensive and not feasible for extension to the microvascular network at large spatial scales up to an organ level. To overcome such limitations, here we present machine learning (ML) models that bypass such expensive computations but provide highly accurate and full 3D profiles of the blood velocity, RBC concentration, WSS, and CFL in every vessel in the microvascular network. The ML models, which are based on artificial neural networks and convolution-based U-net models, predict hemodynamic quantities that compare very well against the true data but reduce the prediction time by several orders. This study therefore paves the way for ML to make detailed and accurate hemodynamic predictions in spatially large microvascular networks at an organ-scale.

Partitioning of dense RBC suspensions in single microfluidic bifurcations: role of cell deformability and bifurcation angle

Red blood cells (RBCs) are a key determinant of human physiology and their behaviour becomes extremely heterogeneous as they navigate in narrow, bifurcating vessels in the microvasculature, affecting local haemodynamics. This is due to partitioning in bifurcations which is dependent on the biomechanical properties of RBCs, especially deformability. We examine the effect of deformability on the haematocrit distributions of dense RBC suspensions flowing in a single, asymmetric Y-shaped bifurcation, experimentally. Human RBC suspensions (healthy and artificially hardened) at 20% haematocrit (Ht) were perfused through the microchannels at different flow ratios between the outlet branches, and negligible inertia, and imaged to infer cell distributions. Notable differences in the shape of the haematocrit distributions were observed between healthy and hardened RBCs near the bifurcation apex. These lead to more asymmetric distributions for healthy RBCs in the daughter and outlet branches with cells accumulating near the inner channel walls, exhibiting distinct hematocrit peaks which are sharper for healthy RBCs. Although the hematocrit distributions differed locally, similar partitioning characteristics were observed for both suspensions. Comparisons with RBC distributions measured in a T-shaped bifurcation showed that the bifurcation angle affects the haematocrit characteristics of the healthy RBCs and not the hardened ones. The extent of RBC partitioning was found similar in both geometries and suspensions. The study highlights the differences between local and global characteristics which impact RBC distribution in more complex, multi-bifurcation networks.

Lift at low Reynolds number

Lift forces are widespread in hydrodynamics. These are typically observed for big and fast objects and are often associated with a combination of fluid inertia (i.e. large Reynolds numbers) and specific symmetry-breaking mechanisms. In contrast, the properties of viscosity-dominated (i.e. low Reynolds numbers) flows make it more difficult for such lift forces to emerge. However, the inclusion of boundary effects qualitatively changes this picture. Indeed, in the context of soft and biological matter, recent studies have revealed the emergence of novel lift forces generated by boundary softness, flow gradients and/or surface charges. The aim of the present review is to gather and analyse this corpus of literature, in order to identify and unify the questioning within the associated communities, and pave the way towards future research.

Effect of collision, size, and oscillation of RBCs on blood heat transfer in a bifurcated vessel

This study attempts to analyze the effect of red blood cells (RBCs) on blood heat transfer in a three-dimensional bifurcated vessel when treated with hyperthermia procedure. A two-phase granular model is used in this paper to study the various underlying factors that affect the flow dynamics of RBCs in a blood vessel. Separate cases are analyzed to study the effect of RBC size, RBC-RBC, and RBC-wall collision and oscillation on heat transfer in a three-dimensional bifurcated vessel under pulsatile flow condition. Blood temperature and transient Nusselt number are used as heat transfer representative parameters. A good agreement with the experimental results from the existing literature is observed when the numerical model used in this study is compared for accuracy. From this study, it has been found that an increase in the size of RBCs of a blood disorder patient can decrease the temperature of blood compared to a normal patient when subjected to hyperthermia treatment. A change in the nature of collision between RBCs does not affect the heat transfer of blood under pulsatile flow condition in a bifurcated vessel. Also, an increase in granular temperature or oscillation of RBCs slightly increases the blood temperature when exposed to thermal treatment.

Cell-free layer development and spatial organization of healthy and rigid red blood cells in a microfluidic bifurcation

Bifurcations and branches in the microcirculation dramatically affect blood flow as they determine the spatiotemporal organization of red blood cells (RBCs). Such changes in vessel geometries can further influence the formation of a cell-free layer (CFL) close to the vessel walls. Biophysical cell properties, such as their deformability, which is impaired in various diseases, are often thought to impact blood flow and affect the distribution of flowing RBCs. This study investigates the flow behavior of healthy and artificially hardened RBCs in a bifurcating microfluidic T-junction. We determine the RBC distribution across the channel width at multiple positions before and after the bifurcation. Thus, we reveal distinct focusing profiles in the feeding mother channel for rigid and healthy RBCs that dramatically impact the cell organization in the successive daughter channels. Moreover, we experimentally show how the characteristic asymmetric CFLs in the daughter vessels develop along their flow direction. Complimentary numerical simulations indicate that the buildup of the CFL is faster for healthy than for rigid RBCs. Our results provide fundamental knowledge to understand the partitioning of rigid RBC as a model of cells with pathologically impaired deformability in complex in vitro networks.

Image-based flow simulation of platelet aggregates under different shear rates

Hemodynamics is crucial for the activation and aggregation of platelets in response to flow-induced shear. In this paper, a novel image-based computational model simulating blood flow through and around platelet aggregates is presented. The microstructure of aggregates was captured by two different modalities of microscopy images of in vitro whole blood perfusion experiments in microfluidic chambers coated with collagen. One set of images captured the geometry of the aggregate outline, while the other employed platelet labelling to infer the internal density. The platelet aggregates were modelled as a porous medium, the permeability of which was calculated with the Kozeny-Carman equation. The computational model was subsequently applied to study hemodynamics inside and around the platelet aggregates. The blood flow velocity, shear stress and kinetic force exerted on the aggregates were investigated and compared under 800 s-1, 1600 s-1 and 4000 s-1 wall shear rates. The advection-diffusion balance of agonist transport inside the platelet aggregates was also evaluated by local Péclet number. The findings show that the transport of agonists is not only affected by the shear rate but also significantly influenced by the microstructure of the aggregates. Moreover, large kinetic forces were found at the transition zone from shell to core of the aggregates, which could contribute to identifying the boundary between the shell and the core. The shear rate and the rate of elongation flow were investigated as well. The results imply that the emerging shapes of aggregates are highly correlated to the shear rate and the rate of elongation. The framework provides a way to incorporate the internal microstructure of the aggregates into the computational model and yields a better understanding of the hemodynamics and physiology of platelet aggregates, hence laying the foundation for predicting aggregation and deformation under different flow conditions.

Influence of the vessel wall geometry on the wall-induced migration of red blood cells

The geometry of the blood vessel wall plays a regulatory role on the motion of red blood cells (RBCs). The overall topography of the vessel wall depends on many features, among which the endothelial lining of the endothelial surface layer (ESL) is an important one. The endothelial lining of vessel walls presents a large surface area for exchanging materials between blood and tissues. The ESL plays a critical role in regulating vascular permeability, hindering leukocyte adhesion as well as inhibiting coagulation during inflammation. Changes in the ESL structure are believed to cause vascular hyperpermeability and entrap immune cells during sepsis, which could significantly alter the vessel wall geometry and disturb interactions between RBCs and the vessel wall, including the wall-induced migration of RBCs and the thickening of a cell-free layer. To investigate the influence of the vessel wall geometry particularly changed by the ESL under various pathological conditions, such as sepsis, on the motion of RBCs, we developed two models to represent the ESL using the immersed boundary method in two dimensions. In particular, we used simulations to study how the lift force and drag force on a RBC near the vessel wall vary with different wall thickness, spatial variation, and permeability associated with changes in the vessel wall geometry. We find that the spatial variation of the wall has a significant effect on the wall-induced migration of the RBC for a high permeability, and that the wall-induced migration is significantly inhibited as the vessel diameter is increased.

Semi-automated red blood cell core detection in blood micro-flow

In microcirculation, red blood cells (RBCs) tend to migrate toward the centre of the vessel leaving a region of a cell depleted layer or cell-free layer (CFL) at the vessel wall and a core of RBCs at the centre. This heterogenous distribution of cells has an effect on the blood apparent viscosity and the exchanges of gases and nutrients between the RBCs and the vessel. Understanding the formation of the CFL and obtaining accurate measurement of it is paramount for furthering development of devices such as drug administration. This paper presents a general semi-automatic method to quantify the thickness of the CFL for different channel geometries and image quality. It enables the use of a method based on intensity, a method using the gradient of the intensity, or a method based on spatiotemporal variation. The main features are reported, the performance is demonstrated on experimentally obtained image sets and accuracy is validated using synthetic images with known CFL thickness. A pure automatic detection is limited by the most visually correct using the spatiotemporal method, however proposed thresholding through automatic detection allows for quality controls through manual adjustments. With a semi-automatic approach RBC core variability between 3 % to 8 % was found when the test user was tasked with repeating the analysis of the same set. The presented method provides, users without programming ability with a user-friendly interface that can extract CFL automatically with quality control through manual adjustments.

Modeling Reactive Hyperemia to Better Understand and Assess Microvascular Function: A Review of Techniques

Reactive hyperemia is a well-established technique for the non-invasive evaluation of the peripheral microcirculatory function, measured as the magnitude of limb re-perfusion after a brief period of ischemia. Despite widespread adoption by researchers and clinicians alike, many uncertainties remain surrounding interpretation, compounded by patient-specific confounding factors (such as blood pressure or the metabolic rate of the ischemic limb). Mathematical modeling can accelerate our understanding of the physiology underlying the reactive hyperemia response and guide in the estimation of quantities which are difficult to measure experimentally. In this work, we aim to provide a comprehensive guide for mathematical modeling techniques that can be used for describing the key phenomena involved in the reactive hyperemia response, alongside their limitations and advantages. The reported methodologies can be used for investigating specific reactive hyperemia aspects alone, or can be combined into a computational framework to be used in (pre-)clinical settings.

Influence of storage and buffer composition on the mechanical behavior of flowing red blood cells

On-chip study of blood flow has emerged as a powerful tool to assess the contribution of each component of blood to its overall function. Blood has indeed many functions, from gas and nutrient transport to immune response and thermal regulation. Red blood cells play a central role therein, in particular through their specific mechanical properties, which directly influence pressure regulation, oxygen perfusion, or platelet and white cell segregation toward endothelial walls. As the bloom of in-vitro studies has led to the apparition of various storage and sample preparation protocols, we address the question of the robustness of the results involving cell mechanical behavior against this diversity. The effects of three conservation media (EDTA, citrate, and glucose-albumin-sodium-phosphate) and storage time on the red blood cell mechanical behavior are assessed under different flow conditions: cell deformability by ektacytometry, shape recovery of cells flowing out of a microfluidic constriction, and cell-flipping dynamics under shear flow. The impact of buffer solutions (phosphate-buffered saline and density-matched suspension using iodixanol/Optiprep) are also studied by investigating individual cell-flipping dynamics, relative viscosity of cell suspensions, and cell structuration under Poiseuille flow. Our results reveal that storing blood samples up to 7 days after withdrawal and suspending them in adequate density-matched buffer solutions has, in most experiments, a moderate effect on the overall mechanical response, with a possible rapid evolution in the first 3 days after sample collection.

A microfluidic method to investigate platelet mechanotransduction under extensional strain

Background

Blood platelets have evolved a complex mechanotransduction machinery to rapidly respond to hemodynamic conditions. A variety of microfluidic flow-based approaches have been developed to explore platelet mechanotransduction; however, these experimental models primarily focus on the effects of increased wall shear stress on platelet adhesion events and do not consider the critical effects of extensional strain on platelet activation in free flow.

Objectives

We report the development and application of a hyperbolic microfluidic assay that allows for investigation of platelet mechanotransduction under quasi-homogenous extensional strain rates in the absence of surface adhesions.

Methods

Using a combined computational fluid dynamic and experimental microfluidic approach, we explore 5 extensional strain regimes (geometries) and their effect on platelet calcium signal transduction.

Results

We demonstrate that in the absence of canonical adhesion, receptor engagement platelets are highly sensitive to both initial increase and subsequent decrease in extensional strain rates within the range of 747 to 3319/s. Furthermore, we demonstrate that platelets rapidly respond to the rate of change in extensional strain and define a threshold of ≥7.33 × 106/s/m, with an optimal range of 9.21 × 107 to 1.32 × 108/s/m. In addition, we demonstrate a key role of both the actin-based cytoskeleton and annular microtubules in the modulation of extensional strain-mediated platelet mechanotransduction.

Conclusion

This method opens a window onto a novel platelet signal transduction mechanism and may have potential diagnostic utility in the identification of patients who are prone to thromboembolic complications associated with high-grade arterial stenosis or are on mechanical circulatory support systems, for which the extensional strain rate is a predominant hemodynamic driver.

Flow-induced "waltzing" red blood cells: Microstructural reorganization and the corresponding rheological response

We investigate flow-induced structural organization in a dilute suspension of tumbling red blood cells (RBCs) under confined shear flow. For small Reynolds (Re = 0.1) and capillary numbers (Ca), with fully coupled hydrodynamic interaction (HI) and without interparticle adhesion, we find that HI between the biconcave discoid particles prompts the formation of layered RBC chains and synchronized rotating RBC pairs, referred here as "waltzing doublets." As the volume fraction ϕ increases, more waltzing doublets appear in RBC files. Stronger shear stress disrupts structural arrangements at higher Ca. We find that the flow-induced organization of waltzing doublets changes how the suspension viscosity varies with ϕ qualitatively. The intrinsic viscosity is particularly sensitive to microstructural rearrangement, increasing (decreasing) with ϕ at low (high) Ca that correlates with the change in the fraction of doublets. We verified flow-induced collective motion with comparison to two-cell simulations in which the cell volume fraction is controlled by varying the domain volume.

Shear induced diffusion of platelets revisited

The transport of platelets in blood is commonly assumed to obey an advection-diffusion equation with a diffusion constant given by the so-called Zydney-Colton theory. Here we reconsider this hypothesis based on experimental observations and numerical simulations including a fully resolved suspension of red blood cells and platelets subject to a shear. We observe that the transport of platelets perpendicular to the flow can be characterized by a non-trivial distribution of velocities with and exponential decreasing bulk, followed by a power law tail. We conclude that such distribution of velocities leads to diffusion of platelets about two orders of magnitude higher than predicted by Zydney-Colton theory. We tested this distribution with a minimal stochastic model of platelets deposition to cover space and time scales similar to our experimental results, and confirm that the exponential-powerlaw distribution of velocities results in a coefficient of diffusion significantly larger than predicted by the Zydney-Colton theory.

Shear stress and intravascular pressure effects on vascular dynamics: two-phase blood flow in elastic microvessels accounting for the passive stresses

We study the steady hemodynamics in physiological elastic microvessels proposing an advanced fluid-structure interaction model. The arteriolar tissue is modeled as a two-layer fiber-reinforced hyperelastic material representing its Media and Adventitia layers. The constitutive model employed (Holzapfel et al. in J Elast 61:1-48, 2000) is parametrized via available data on stress-strain experiments for arterioles. The model is completed by simulating the blood/plasma flow in the lumen, using the thixotropic elasto-viscoplastic model in its core, and the linear Phan-Thien and Tanner viscoelastic model in its annular part. The Cell-Free Layer (CFL) and the Fåhraeus and Fåhraeus-Lindqvist effects are considered via analytical expressions based on experimental data (Giannokostas et al. in Materials (Basel) 14:367, 2021b). The coupling between tissue deformation and blood flow is achieved through the experimentally verified pressure-shear hypothesis (Pries et al. Circ Res 77:1017-1023, 1995). Our calculations confirm that the increase in the reference inner radius produces larger expansion. Also, by increasing the intraluminal pressure, the thinning of the walls is more pronounced and it may reach 40% of the initial thickness. Comparing our predictions with those in rigid-wall microtubes, we conclude that apart from the vital importance of vasodilation, there is an up to 25% reduction in wall shear stress. The passive vasodilation contributes to the decrease in the tissue stress fields and affects the hemodynamic features such as the CFL thickness, reducing the plasma layer when blood flows in vessels with elastic walls, in quantitative agreement with previous experiments. Our calculations verify the correctness of the pressure-shear hypothesis but not that of the Laplace law.

High-Throughput Imaging of Blood Flow Reveals Developmental Changes in Distribution Patterns of Hemodynamic Quantities in Developing Zebrafish

Mechanical forces from blood flow and pressure (hemodynamic forces) contribute to the formation and shaping of the blood vascular network during embryonic development. Previous studies have demonstrated that hemodynamic forces regulate signaling and gene expression in endothelial cells that line the inner surface of vascular tubes, thereby modifying their cellular state and behavior. Given its important role in vascular development, we still know very little about the quantitative aspects of hemodynamics that endothelial cells experience due to the difficulty in measuring forces in vivo. In this study, we sought to determine the magnitude of wall shear stress (WSS) exerted on ECs by blood flow in different vessel types and how it evolves during development. Utilizing the zebrafish as a vertebrate model system, we have established a semi-automated high-throughput fluorescent imaging system to capture the flow of red blood cells in an entire zebrafish between 2- and 6-day post-fertilization (dpf). This system is capable of imaging up to 50 zebrafish at a time. A semi-automated analysis method was developed to calculate WSS in zebrafish trunk vessels. This was achieved by measuring red blood cell flow using particle tracking velocimetry analysis, generating a custom-made script to measure lumen diameter, and measuring local tube hematocrit levels to calculate the effective blood viscosity at each developmental stage. With this methodology, we were able to determine WSS magnitude in different vessels at different stages of embryonic and larvae growth and identified developmental changes in WSS, with absolute levels of peak WSS in all vessel types falling to levels below 0.3 Pa at 6 dpf. Additionally, we discovered that zebrafish display an anterior-to-posterior trend in WSS at each developmental stage.

Ultrafast two-photon fluorescence imaging of cerebral blood circulation in the mouse brain in vivo

SignificanceCharacterizing blood flow by tracking individual red blood cells as they move through vessels is essential for understanding vascular function. With high spatial resolution, two-photon fluorescence microscopy is the method of choice for imaging blood flow at the cellular level. However, its application is limited to a low flow speed regimen in anesthetized animals by its slow focus scanning mechanism. Using an ultrafast scanning module, we demonstrated two-photon fluorescence imaging of blood flow at 1,000 two-dimensional frames and 1,000,000 one-dimensional line scans per second in the brains of awake mice. These ultrafast measurements enabled us to study hemodynamic and fluid mechanical regimens beyond the reach of conventional methods.

Assessment of Blood Biophysical Properties Using Pressure Sensing with Micropump and Microfluidic Comparator

To identify the biophysical properties of blood samples consistently, macroscopic pumps have been used to maintain constant flow rates in a microfluidic comparator. In this study, the bulk-sized and expensive pump is replaced with a cheap and portable micropump. A specific reference fluid (i.e., glycerin solution [40%]) with a small volume of red blood cell (RBC) (i.e., 1% volume fraction) as fluid tracers is supplied into the microfluidic comparator. An averaged velocity () obtained with micro-particle image velocimetry is converted into the flow rate of reference fluid (Qr) (i.e., Qr = CQ × Ac × , Ac: cross-sectional area, CQ = 1.156). Two control variables of the micropump (i.e., frequency: 400 Hz and volt: 150 au) are selected to guarantee a consistent flow rate (i.e., COV < 1%). Simultaneously, the blood sample is supplied into the microfluidic channel under specific flow patterns (i.e., constant, sinusoidal, and periodic on-off fashion). By monitoring the interface in the comparator as well as Qr, three biophysical properties (i.e., viscosity, junction pressure, and pressure-induced work) are obtained using analytical expressions derived with a discrete fluidic circuit model. According to the quantitative comparison results between the present method (i.e., micropump) and the previous method (i.e., syringe pump), the micropump provides consistent results when compared with the syringe pump. Thereafter, representative biophysical properties, including the RBC aggregation, are consistently obtained for specific blood samples prepared with dextran solutions ranging from 0 to 40 mg/mL. In conclusion, the present method could be considered as an effective method for quantifying the physical properties of blood samples, where the reference fluid is supplied with a cheap and portable micropump.

Computational investigation of blood cell transport in retinal microaneurysms

Microaneurysms (MAs) are one of the earliest clinically visible signs of diabetic retinopathy (DR). MA leakage or rupture may precipitate local pathology in the surrounding neural retina that impacts visual function. Thrombosis in MAs may affect their turnover time, an indicator associated with visual and anatomic outcomes in the diabetic eyes. In this work, we perform computational modeling of blood flow in microchannels containing various MAs to investigate the pathologies of MAs in DR. The particle-based model employed in this study can explicitly represent red blood cells (RBCs) and platelets as well as their interaction in the blood flow, a process that is very difficult to observe in vivo. Our simulations illustrate that while the main blood flow from the parent vessels can perfuse the entire lumen of MAs with small body-to-neck ratio (BNR), it can only perfuse part of the lumen in MAs with large BNR, particularly at a low hematocrit level, leading to possible hypoxic conditions inside MAs. We also quantify the impacts of the size of MAs, blood flow velocity, hematocrit and RBC stiffness and adhesion on the likelihood of platelets entering MAs as well as their residence time inside, two factors that are thought to be associated with thrombus formation in MAs. Our results show that enlarged MA size, increased blood velocity and hematocrit in the parent vessel of MAs as well as the RBC-RBC adhesion promote the migration of platelets into MAs and also prolong their residence time, thereby increasing the propensity of thrombosis within MAs. Overall, our work suggests that computational simulations using particle-based models can help to understand the microvascular pathology pertaining to MAs in DR and provide insights to stimulate and steer new experimental and computational studies in this area.

Dependence of cell-free-layer width on rheological parameters: Combining empirical data on flow separation at microvascular bifurcations with geometrical considerations

When blood flows through vessel networks, red blood cells (RBCs) are typically concentrated close to the vessel center line, creating a lubrication layer near the vessel wall. As RBCs bind oxygen, the width of this cell-free layer (CFL) impacts not only the blood rheology inside the vasculature, but also oxygen delivery to the tissues they perfuse and, hence, their function. Existing attempts to relate the width of the CFL to the rheological properties of the blood and the geometrical properties of the vessel are based on an analysis of the forces acting on RBCs suspended in the blood. However, the complexity of interactions in the blood makes this a challenging task. Here, we propose an alternative, two-step approach to derive such a functional relationship. First, we extend widely accepted empirical fits describing the minimum flow fraction needed for RBCs to enter a daughter vessel downstream of a microvascular bifurcation so that it depends not only on the diameter and discharge haematocrit of the parent vessel, but also on its average shear rate. Second, we propose a simple geometrical model for the minimum flow fraction based on the cross-sectional blood flow profile within the parent vessel upstream of the bifurcation-considering uniform, parabolic, and blunt velocity profiles-and derive the leading-order approximation to this model for small CFL widths. By equating the functional relationships obtained using these two approaches, we derive expressions relating the CFL width to the vessel diameter, discharge haematocrit, and mean shear rate. The resulting expressions are in good agreement with available in vivo data and represent a promising basis for future research.

Recent advances in blood rheology: a review

Due to the potential impact on the diagnosis and treatment of various cardiovascular diseases, work on the rheology of blood has significantly expanded in the last decade, both experimentally and theoretically. Experimentally, blood has been confirmed to demonstrate a variety of non-Newtonian rheological characteristics, including pseudoplasticity, viscoelasticity, and thixotropy. New rheological experiments and the development of more controlled experimental protocols on more extensive, broadly physiologically characterized, human blood samples demonstrate the sensitivity of aspects of hemorheology to several physiological factors. For example, at high shear rates the red blood cells elastically deform, imparting viscoelasticity, while at low shear rates, they form "rouleaux" structures that impart additional, thixotropic behavior. In addition to the advances in experimental methods and validated data sets, significant advances have also been made in both microscopic simulations and macroscopic, continuum, modeling, as well as novel, multiscale approaches. We outline and evaluate the most promising of these recent developments. Although we primarily focus on human blood rheology, we also discuss recent observations on variations observed across some animal species that provide some indication on evolutionary effects.

Emergent cell-free layer asymmetry and biased haematocrit partition in a biomimetic vascular network of successive bifurcations

Blood is a vital soft matter, and its normal circulation in the human body relies on the distribution of red blood cells (RBCs) at successive bifurcations. Understanding how RBCs are partitioned at bifurcations is key for the optimisation of microfluidic devices as well as for devising novel strategies for diagnosis and treatment of blood-related diseases. We report the dynamics of RBC suspensions flowing through a biomimetic vascular network incorporating three generations of microchannels and two classical types of bifurcations at the arteriole level. Our microfluidic experiments with dilute and semidilute RBC suspensions demonstrate the emergence of excessive heterogeneity of RBC concentration in downstream generations upon altering the network's outflow rates. Through parallel simulations using the immersed-boundary-lattice-Boltzmann method, we reveal that the heterogeneity is attributed to upstream perturbations in the cell-free layer (CFL) and lack of its recovery between consecutive bifurcations owing to suppressed hydrodynamic lift under reduced flow conditions. In the dilute/semidilute regime, this perturbation dominates over the effect of local fractional flow at the bifurcation and can lead to inherently unfavourable child branches that are deprived of RBCs even for equal flow split. Our work highlights the importance of CFL asymmetry cascading down a vascular network, which leads to biased phase separation that deviates from established empirical predictions.

Hemodynamic analysis for stenosis microfluidic model of thrombosis with refined computational fluid dynamics simulation

Disturbed blood flow has been increasingly recognized for its critical role in platelet aggregation and thrombosis. Microfluidics with hump shaped contractions have been developed to mimic microvascular stenosis and recapitulate the prothrombotic effect of flow disturbance. However the physical determinants of microfluidic hemodynamics are not completely defined. Here, we report a refined computational fluid dynamics (CFD) simulation approach to map the shear rate (γ) and wall shear stress (τ) distribution in the stenotic region at high accuracy. Using ultra-fine meshing with sensitivity verification, our CFD results show that the stenosis level (S) is dominant over the bulk shear rate (γ₀) and contraction angle (α) in determining γ and τ distribution at stenosis. In contrast, α plays a significant role in governing the shear rate gradient (γ^′) distribution while it exhibits subtle effects on the peak γ. To investigate the viscosity effect, we employ a Generalized Power-Law model to simulate blood flow as a non-Newtonian fluid, showing negligible difference in the γ distribution when compared with Newtonian simulation with water medium. Together, our refined CFD method represents a comprehensive approach to examine microfluidic hemodynamics in three dimensions and guide microfabrication designs. Combining this with hematological experiments promises to advance understandings of the rheological effect in thrombosis and platelet mechanobiology.

Development of margination of platelet-sized particles in red blood cell suspensions flowing through Y-shaped bifurcating microchannels

BACKGROUND

In the blood flow through microvessels, platelets exhibit enhanced concentrations in the layer free of red blood cells (cell-free layer) adjacent to the vessel wall. The motion of platelets in the cell-free layer plays an essential role in their interaction with the vessel wall, and hence it affects their functions of hemostasis and thrombosis.

OBJECTIVE

We aimed to estimate the diffusivity of platelet-sized particles in the transverse direction (the direction of vorticity) across the channel width in the cell-free layer by in vitro experiments for the microchannel flow of red blood cell (RBC) suspensions containing platelet-sized particles.

METHODS

Fluorescence microscope observations were performed to measure the transverse distribution of spherical particles immersed in RBC suspensions flowing through a Y-shaped bifurcating microchannel. We examined the development of the particle concentration profiles along the flow direction in the daughter channels, starting from asymmetric distributions with low concentrations on the inner side of the bifurcation at the inlet of the daughter channels.

RESULTS

In daughter channels of 40 μm width, reconstruction of particle margination revealed that a symmetric concentration profile was attained in ∼30 mm from the bifurcation, independent of flow rate.

CONCLUSIONS

We presented experimental evidence of particle margination developing in a bifurcating flow channel where the diffusivity of 2.9-μm diameter particles was estimated to be ∼40 μm2/s at a shear rate of 1000 s-1 and hematocrit of 0.2.

Integrating blood cell mechanics, platelet adhesive dynamics and coagulation cascade for modelling thrombus formation in normal and diabetic blood

Normal haemostasis is an important physiological mechanism that prevents excessive bleeding during trauma, whereas the pathological thrombosis especially in diabetics leads to increased incidence of heart attacks and strokes as well as peripheral vascular events. In this work, we propose a new multiscale framework that integrates seamlessly four key components of blood clotting, namely transport of coagulation factors, coagulation kinetics, blood cell mechanics and platelet adhesive dynamics, to model the development of thrombi under physiological and pathological conditions. We implement this framework to simulate platelet adhesion due to the exposure of tissue factor in a three-dimensional microchannel. Our results show that our model can simulate thrombin-mediated platelet activation in the flowing blood, resulting in platelet adhesion to the injury site of the channel wall. Furthermore, we simulate platelet adhesion in diabetic blood, and our results show that both the pathological alterations in the biomechanics of blood cells and changes in the amount of coagulation factors contribute to the excessive platelet adhesion and aggregation in diabetic blood. Taken together, this new framework can be used to probe synergistic mechanisms of thrombus formation under physiological and pathological conditions, and open new directions in modelling complex biological problems that involve several multiscale processes.

Red blood cell hitchhiking enhances the accumulation of nano- and micro-particles in the constriction of a stenosed microvessel

We investigate the circulation of nano- and micro-particles, including spherical particles and filamentous nanoworms, with red blood cells (RBCs) suspension in a constricted channel that mimics a stenosed microvessel. Through three-dimensional simulations using the immersed boundary-based Lattice Boltzmann method, the influence of channel geometries, such as the length and ratio of the constriction, on the accumulation of particles is systematically studied. Firstly, we find that the accumulation of spherical particles with 1 μm diameter in the constriction increases with the increases of both the length and ratio of the constriction. This is attributed to the interaction between spheres and RBCs. The RBCs "carry" the spheres and they accumulate inside the constriction together, due to the altered local hydrodynamics induced by the existence of the constriction. Secondly, nanoworms demonstrate higher accumulation than that of spheres inside the constriction, which is associated with the escape of nanoworms from RBC clusters and their accumulation near the wall of main channel. The accumulated near-wall nanoworms will eventually enter the constriction, thus enhancing their concentration inside the constriction. However, an exceptional case occurs in the case of constrictions with large ratio and long length. In such circumstances, the RBCs aggregate together tightly and concentrate at the center of the channel, which makes the nanoworms hardly able to escape from RBC clusters, leading to a similar accumulation of nanoworms and spheres inside the constriction. This study may provide theoretical guidance for the design of nano- and micro-particles for biomedical engineering applications, such as drug delivery systems for patients with stenosed microvessels.

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.

Effects of ankle pump exercise frequency on venous hemodynamics of the lower limb

PURPOSE

The purpose of this study was (1) to compare the effects of ankle pump exercise(APE) frequency on hemodynamics of the common femoral vein(CFV) and (2) to analyse the relationship between APE duration and lower limb fatigue.

METHODS

Twenty-seven males and thirty-three females performed APE. Among them, there were thirty participants with non-lower limb fracture (N-LLF) and thirty participants with lower limb fracture (LLF).The colour doppler ultrasound was used to record the time-averaged mean velocity (TAMV) of common femoral venous flow when the participants at rest and move at different frequencies of 6 times/min, 10 times/min, 30 times/min, 60 times/min. The ratings of perceived exertion (RPE) was used to assess the lower limb fatigue of the participants when performing APE at the frequency of 60 times/min after 1 min, 2 mins, 3 mins, 4 mins and 5 mins.

RESULTS

With the increase of frequency, TAMV increased significantly in both the participants with N-LLF and LLF (p < 0.01). The TAMV were 19.82±3.86, 33.78±8.76, 37.06±8.67, 43.82±10.40, 52.18±10.53, respectively in the participants with N-LLF and 16.98±3.01, 22.20±4.96, 24.01±5.78, 29.20±7.05, 35.75±9.28, respectively in the injured limb of patients with LLF when at rest and moving at the frequency of 6 times/min, 10 times/min, 30 times/min, 60 times/min. There was a positive correlation (p < 0.01) between lower limb fatigue and exercise duration.When the RPE was 16 points (the corresponding fatigue degree is "tired"), the exercise duration of the participants with N-LLF was 3 mins and that of the injured limb of patients with LLF was 2 mins.

CONCLUSIONS

Both fast and slow-frequency APE can promote venous blood return in the lower limb. Despite of the equivalent APE duration, fast-frequency APE can promote venous blood return more effectively. When the frequency of APE was 60 times/min, participants with N-LLF can exercise for 3 mins, and the injured limb of patients with LLF can exercise for 2 mins.

Spatiotemporal Dynamics of Dilute Red Blood Cell Suspensions in Low-Inertia Microchannel Flow

Microfluidic technologies are commonly used for the manipulation of red blood cell (RBC) suspensions and analyses of flow-mediated biomechanics. To enhance the performance of microfluidic devices, understanding the dynamics of the suspensions processed within is crucial. We report novel, to our knowledge, aspects of the spatiotemporal dynamics of RBC suspensions flowing through a typical microchannel at low Reynolds number. Through experiments with dilute RBC suspensions, we find an off-center two-peak (OCTP) profile of cells contrary to the centralized distribution commonly reported for low-inertia flows. This is reminiscent of the well-known "tubular pinch effect," which arises from inertial effects. However, given the conditions of negligible inertia in our experiments, an alternative explanation is needed for this OCTP profile. Our massively parallel simulations of RBC flow in real-size microfluidic dimensions using the immersed-boundary-lattice-Boltzmann method confirm the experimental findings and elucidate the underlying mechanism for the counterintuitive RBC pattern. By analyzing the RBC migration and cell-free layer development within a high-aspect-ratio channel, we show that such a distribution is co-determined by the spatial decay of hydrodynamic lift and the global deficiency of cell dispersion in dilute suspensions. We find a cell-free layer development length greater than 46 and 28 hydraulic diameters in the experiment and simulation, respectively, exceeding typical lengths of microfluidic designs. Our work highlights the key role of transient cell distribution in dilute suspensions, which may negatively affect the reliability of experimental results if not taken into account.

Flow-induced segregation and dynamics of red blood cells in sickle cell disease

Blood flow in sickle cell disease (SCD) can substantially differ from normal blood flow due to significant alterations in the physical properties of the red blood cells (RBCs). Chronic complications, such as inflammation of the endothelial cells lining blood vessel walls, are associated with SCD, for reasons that are unclear. Here, detailed boundary integral simulations are performed to investigate an idealized model flow flow in SCD, a binary suspension of flexible biconcave discoidal fluid-filled capsules and stiff curved prolate capsules that represent healthy and sickle RBCs, respectively, subjected to pressure-driven flow in a planar slit. The stiff component is dilute. The key observation is that, unlike healthy RBCs that concentrate around the center of the channel and form an RBC-depleted layer (i.e. cell-free layer) next to the walls, sickle cells are largely drained from the bulk of the suspension and aggregate inside the cell-free layer, displaying strong margination. These cells are found to undergo a rigid-body-like rolling orbit near the walls. A binary suspension of flexible biconcave discoidal capsules and stiff straight (non-curved) prolate capsules is also considered for comparison, and the curvature of the stiff component is found to play a minor role in the behavior. Additionally, by considering a mixture of flexible and stiff biconcave discoids, we reveal that rigidity difference by itself is sufficient to induce the segregation behavior in a binary suspension. Furthermore, the additional shear stress on the walls induced by the presence of cells is computed for the various cases. Compared to the small fluctuations in wall shear stress for a suspension of healthy RBCs, large local peaks in wall shear stress are observed for the binary suspensions, due to the proximity of the marginated stiff cells to the walls. This effect is most marked for the straight prolate capsules. As endothelial cells are known to mechanotransduce physical forces such as aberrations in shear stress and convert them to physiological processes such as activation of inflammatory signals, these results may aid in understanding mechanisms for endothelial dysfunction associated with SCD.

A Molecular Communications System for Live Detection of Hyperviscosity Syndrome

The hyperviscosity syndrome consists of an excessive viscosity of the blood. It has to be timely detected, since it significantly increases the risk of many cardiovascular diseases. It can manifest essentially in two different ways, depending on the capability of the cardiocirculatory system to compensate the increased resistance experienced by the blood in vessels. In case of compensation, the blood flow velocity remains nearly unaltered; otherwise a significant decrease can occur. In this paper, we propose a blood viscosity monitoring system based on molecular communications. It could be used in risky patients for continuous monitoring and can detect blood hyperviscosity both in case of compensation, even partial, and in case of uncompensated cardiocirculatory systems. This proposal results to be both simple and effective for real time estimating the blood viscosity level, without any need of resorting to invasive blood analyses. The main idea is to leverage the statistics of particle dispersion and absorption in a blood vessel as a function of blood viscosity. Viscosity is estimated by processing the statistics of the received signal molecules, which are easy to collect. These features open the way to the implementation of the proposed solution in a simple and cheap device.

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.

Reduced-order modeling of hemodynamics across macroscopic through mesoscopic circulation scales

We propose a hemodynamic reduced-order model bridging macroscopic and mesoscopic blood flow circulation scales from arteries to capillaries. In silico tree-like vascular geometries, mathematically described by graphs, are synthetically generated by means of stochastic growth algorithms constrained by statistical morphological and topological principles. Scale-specific pruning gradation of the tree is then proposed in order to fit computational budget requirement. Different compliant structural models with respect to pressure loads are used depending on vessel walls thicknesses and structures, which vary considerably from macroscopic to mesoscopic circulation scales. Nonlinear rheological properties of blood are also included, and microcirculation network responses are computed for different rheologies. Numerical results are in very good agreement with available experimental measurements. The computational model captures the dynamic transition between large- to small-scale flow pulsatility speeds and magnitudes and wall shear stresses, which have wide-ranging physiological influences.

Shear-induced non-monotonic viscosity dependence for model red blood cell suspensions in microvessels

The cell-free layer thickness of an aggregating red blood cell (RBC) suspension in a rectangular microchannel is investigated by hybrid fluid-particle numerical modeling. Several factors affect the suspension viscosity, cell-free layer thickness, and the cell aggregate distribution. These include the hematocrit, vessel size, red cell stiffness, aggregation interaction, and shear rate. In particular, the effect of the shear rate on the cell-free layer thickness is controversial. We found that the suspension viscosity increases along with a decrease in the cell-free layer thickness as the shear rate increases for aggregating model RBCs at low shear rates. At moderate to high shear rates, the cell-free layer thickness increases with the increasing shear rate from medium to strong shear flow for both 10% and 20% red blood cell suspensions.

Identifying the start of a platelet aggregate by the shear rate and the cell-depleted layer

Computer simulations were performed to study the transport of red blood cells and platelets in high shear flows, mimicking earlier published in vitro experiments in microfluidic devices with high affinity for platelet aggregate formation. The goal is to understand and predict where thrombus formation starts. Additionally, the need of cell-based modelling in these microfluidic devices is demonstrated by comparing our results with macroscopic models, wherein blood is modelled as a continuous fluid. Hemocell, a cell-based blood flow simulation framework is used to investigate the transport physics in the microfluidic devices. The simulations show an enlarged cell-depleted layer at the site where a platelet aggregate forms in the experiments. In this enlarged cell-depleted layer, the probability to find a platelet is higher than in the rest of the microfluidic device. In addition, the shear rates are sufficiently high to allow for the von Willebrand factor to elongate in this region. We hypothesize that the enlarged cell-depleted layer combined with a sufficiently large platelet flux and sufficiently high shear rates result in an haemodynamic environment that is a preferred location for initial platelet aggregation.

Accounting for Tube Hematocrit in Modeling of Blood Flow in Cerebral Capillary Networks

The aim of this paper consists in the derivation of an analytic formula for the hydraulic resistance of capillaries, taking into account the tube hematocrit level. The consistency of the derived formula is verified using Finite Element simulations. Such an effective formula allows for assigning resistances, depending on the hematocrit level, to the edges of networks modeling biological capillary systems, which extends our earlier models of blood flow through large capillary networks. Numerical simulations conducted for large capillary networks with random topologies demonstrate the importance of accounting for the hematocrit level for obtaining consistent results.

Antimargination of Microparticles and Platelets in the Vicinity of Branching Vessels

We investigate the margination of microparticles/platelets in blood flow through complex geometries typical for in vivo vessel networks: a vessel confluence and a bifurcation. Using three-dimensional lattice Boltzmann simulations, we confirm that behind the confluence of two vessels, a cell-free layer devoid of red blood cells develops in the channel center. Despite its small size of roughly 1 μm, this central cell-free layer persists for up to 100 μm after the confluence. Most importantly, we show from simulations that this layer also contains a significant amount of microparticles/platelets and validate this result by in vivo microscopy in mouse venules. At bifurcations, however, a similar effect does not appear, and margination is largely unaffected by the geometry. This antimargination toward the vessel center after a confluence may explain earlier in vivo observations, which found that platelet concentrations near the vessel wall are seen to be much higher on the arteriolar side (containing bifurcations) than on the venular side (containing confluences) of the vascular system.

Experimental and Numerical Study of Blood Flow in μ-vessels: Influence of the Fahraeus–Lindqvist Effect

The study of hemodynamics is particularly important in medicine and biomedical engineering as it is crucial for the design of new implantable devices and for understanding the mechanism of various diseases related to blood flow. In this study, we experimentally identify the cell free layer (CFL) width, which is the result of the Fahraeus–Lindqvist effect, as well as the axial velocity distribution of blood flow in microvessels. The CFL extent was determined using microscopic photography, while the blood velocity was measured by micro-particle image velocimetry (μ-PIV). Based on the experimental results, we formulated a correlation for the prediction of the CFL width in small caliber (D < 300 μm) vessels as a function of a modified Reynolds number (Re∞) and the hematocrit (Hct). This correlation along with the lateral distribution of blood viscosity were used as input to a “two-regions” computational model. The reliability of the code was checked by comparing the experimentally obtained axial velocity profiles with those calculated by the computational fluid dynamics (CFD) simulations. We propose a methodology for calculating the friction loses during blood flow in μ-vessels, where the Fahraeus–Lindqvist effect plays a prominent role, and show that the pressure drop may be overestimated by 80% to 150% if the CFL is neglected.

Effects of red blood cell aggregation on the blood flow in a symmetrical stenosed microvessel

In order to figure out whether red blood cell (RBC) aggregation is beneficial or deleterious for the blood flow through a stenosis, fluid mechanics of a microvascular stenosis was examined through simulating the dynamics of deformable red blood cells suspended in plasma using dissipative particle dynamics. The spatial variation in time-averaged cell-free layer (CFL) thickness and velocity profiles indicated that the blood flow exhibits asymmetry along the flow direction. The RBC accumulation occurs upstream the stenosis, leading to a thinner CFL and reduced flow velocity. Therefore, the emergence of stenosis produces an increased blood flow resistance. In addition, an enhanced Fahraeus-Lindqvist effect was observed in the presence of the stenosis. Finally, the effect of RBC aggregation combined with decreased stenosis on the blood flow was investigated. The findings showed that when the RBC clusters pass through the stenosis with a throat comparable to the RBC core in diameter, the blood flow resistance decreases with increasing intercellular interaction strength. But if the RBC core is larger and even several times than the throat, the blood flow resistance increases largely under strong RBC aggregation, which may contribute to the mechanism of the microthrombus formation.

Microfluidic blood vasculature replicas using backside lithography

Blood vessels in living tissues are an organized and hierarchical network of arteries, arterioles, capillaries, veinules and veins. Their sizes, lengths, shapes and connectivity are set up for an optimum perfusion of the tissues in which they deploy. In order to study the hemodynamics and hemophysics of blood flows and also to investigate artificial vasculature for organs on a chip, it is essential to reproduce most of these geometric features. Common microfluidic techniques produce channels with a uniform height and a rectangular cross section that do not capture the size hierarchy observed in vivo. This paper presents a new single-mask photolithography process using an optical diffuser to produce a backside exposure leading to microchannels with both a rounded cross section and a direct proportionality between local height and local width, allowing a one-step design of intrinsically hierarchical networks.

The Margination of Particles in Areas of Constricted Blood Flow

Stroke is a leading cause of death globally and is caused by stenoses, abnormal narrowings of blood vessels. Recently, there has been an increased interest in shear-activated particle clusters for the treatment of stenosis, but there is a lack of literature investigating the impact of different stenosis geometries on particle margination. Margination refers to the movement of particles toward the blood vessel wall and is desirable for drug delivery. The current study investigated ten different geometries and their effects on margination. Microfluidic devices with a constricted area were fabricated to mimic a stenosed blood vessel with different extent of occlusion, constricted length, and eccentricity (gradualness of the constriction and expansion). Spherical fluorescent particles with a diameter of 2.11 μm were suspended in blood and tracked as they moved into, through, and out of the constricted area. A margination parameter, M, was used to quantify margination based on the particle distribution after velocity normalization. Experimental results suggested that a constriction leads to an enhanced margination, whereas an expansion is responsible for a decrease in margination. Further, margination was found to increase with increasing percent occlusion and constriction length, likely a result of higher shear rate and longer residence time, respectively. Margination decreases as the stenosis geometry becomes more gradual (eccentricity increases) with the exception of a sudden constriction/expansion geometry. The findings demonstrate the importance of geometric effects on margination and call for detailed numerical modeling and geometric characterization of the stenosed areas to fully understand the underlying physics.