Blood flow and cell-free layer in microvessels

Non-Newtonian flow of blood in arterioles: consequences for wall shear stress measurements

OBJECTIVE

Our primary goal is to investigate the effects of non-Newtonian blood properties on wall shear stress in microvessels. The secondary goal is to derive a correction factor for the Poiseuille-law-based indirect measurements of wall shear stress.

METHODS

The flow is assumed to exhibit two distinct, immiscible and homogeneous fluid layers: an inner region densely packed with RBCs, and an outer cell-free layer whose thickness depends on discharge hematocrit. The cell-free layer is assumed to be Newtonian, while rheology of the RBC-rich core is modeled using the Quemada constitutive law.

RESULTS

Our model provides a realistic description of experimentally observed blood velocity profiles, tube hematocrit, core hematocrit, and apparent viscosity over a wide range of vessel radii and discharge hematocrits.

CONCLUSIONS

Our analysis reveals the importance of incorporating this complex blood rheology into estimates of WSS in microvessels. The latter is accomplished by specifying a correction factor, which accounts for the deviation of blood flow from the Poiseuille law.

Cystic fibrosis transmembrane conductance regulator is involved in polyphenol-induced swelling of the endothelial glycocalyx

Previous studies show that polyphenol-rich compounds can induce a swelling of the endothelial glycocalyx (eGC). Our goal was to reveal the mechanism behind the eGC-swelling. As polyphenols are potent modulators of fibrosis transmembrane conductance regulator (CFTR) Cl(-) channel, the hypothesis was tested whether polyphenol-induced increase in CFTR activity is responsible for the eGC-swelling. The impact of the polyphenols resveratrol, (-)-epicatechin, and quercetin on nanomechanics of living endothelial GM7373 cells was monitored by AFM-nanoindentation. The tested polyphenols lead to eGC-swelling with a simultaneous decrease in cortical stiffness. EGC-swelling, but not the change in cortical stiffness, was prevented by the inhibition of CFTR. Polyphenol-induced eGC-swelling could be mimicked by cytochalasin D, an actin-depolymerizing agent. Thus, in the vascular endothelium, polyphenols induce eGC-swelling by softening cortical actin and activating CFTR. Our findings imply that CFTR plays an important role in the maintenance of vascular homeostasis and may explain the vasoprotective properties of polyphenols.

FROM THE CLINICAL EDITOR

Many vascular problems clinically can be attributed to a dysregulation of endothelial glycocalyx (eGC). The underlying mechanism however remains unclear. In this article, the authors used nanoindentation and showed that polyphenols could swell the endothelial glycocalyx and alter its function. This investigative method can lead to further mechanistic studies of other molecular pathways.

Microvascular blood flow resistance: Role of red blood cell migration and dispersion

Microvascular blood flow resistance has a strong impact on cardiovascular function and tissue perfusion. The flow resistance in microcirculation is governed by flow behavior of blood through a complex network of vessels, where the distribution of red blood cells across vessel cross-sections may be significantly distorted at vessel bifurcations and junctions. In this paper, the development of blood flow and its resistance starting from a dispersed configuration of red blood cells is investigated in simulations for different hematocrit levels, flow rates, vessel diameters, and aggregation interactions between red blood cells. Initially dispersed red blood cells migrate toward the vessel center leading to the formation of a cell-free layer near the wall and to a decrease of the flow resistance. The development of cell-free layer appears to be nearly universal when scaled with a characteristic shear rate of the flow. The universality allows an estimation of the length of a vessel required for full flow development, lc ≲ 25D, for vessel diameters in the range 10 μm < D < 100 μm. Thus, the potential effect of red blood cell dispersion at vessel bifurcations and junctions on the flow resistance may be significant in vessels which are shorter or comparable to the length lc. Aggregation interactions between red blood cells generally lead to a reduction of blood flow resistance. The simulations are performed using the same viscosity for both external and internal fluids and the RBC membrane viscosity is not considered; however, we discuss how the viscosity contrast may affect the results. Finally, we develop a simple theoretical model which is able to describe the converged cell-free-layer thickness at steady-state flow with respect to flow rate. The model is based on the balance between a lift force on red blood cells due to cell-wall hydrodynamic interactions and shear-induced effective pressure due to cell-cell interactions in flow. We expect that these results can also be used to better understand the flow behavior of other suspensions of deformable particles such as vesicles, capsules, and cells.

Synergy between shear-induced migration and secondary flows on red blood cells transport in arteries: considerations on oxygen transport

Shear-induced migration of red blood cells (RBCs) is a well-known phenomenon characterizing blood flow in the small vessels (micrometre to millimetre size) of the cardiovascular system. In large vessels, like the abdominal aorta and the carotid artery (millimetre to centimetre size), the extent of this migration and its interaction with secondary flows has not been fully elucidated. RBC migration exerts its influence primarily on platelet concentration, oxygen transport and oxygen availability at the luminal surface, which could influence vessel wall disease processes in and adjacent to the intima. Phillips' shear-induced particle migration model, coupled to the Quemada viscosity model, was employed to simulate the macroscopic behaviour of RBCs in four patient-specific geometries: a normal abdominal aorta, an abdominal aortic aneurysm (AAA), a normal carotid bifurcation and a stenotic carotid bifurcation. Simulations show a migration of RBCs from the near-wall region with a lowering of wall haematocrit (volume fraction of RBCs) on the posterior side of the normal aorta and on the lateral-external side of the iliac arteries. A marked migration is observed on the outer wall of the carotid sinus, along the common carotid artery and in the carotid stenosis. No significant migration is observed in the AAA. The spatial and temporal patterns of wall haematocrit are correlated with the near-wall shear layer and with the secondary flows induced by the vessel curvature. In particular, secondary flows accentuate the initial lowering in RBC near-wall concentration by convecting RBCs from the inner curvature side to the outer curvature side. The results reinforce data in literature showing a decrease in oxygen partial pressure on the inner curvature wall of the carotid sinus induced by the presence of secondary flows. The lowering of wall haematocrit is postulated to induce a decrease in oxygen availability at the luminal surface through a diminished concentration of oxyhaemoglobin, hence contributing, with the reported lowered oxygen partial pressure, to local hypoxia.

Cell-free layer development process in the entrance region of microvessels

We simulated red blood cell flows through a finite length channel with a two-dimensional immersed boundary lattice Boltzmann model. The local instantaneous variation in wall-cell distance has been examined in details, and a nominal cell-free layer (CFL) thickness has been proposed. The CFL development process along the channel has been then analyzed, showing that the CFL thickness profile can be basically split into two regimes: the initial rapid increase due to cell migration and the later gradual growth due to cell reorganization. Effects of various hemorheological factors, such as rigidity, aggregation, hematocrit, and channel width, have also been investigated. The development length of the CFL to 90% of its final width ranges from 150 to 300 μm, and the development length is sensitive to changes in hemorheological conditions. The correlation between the CFL features and hemorheological parameters has also been explored. The simulation results have been compared to available experimental studies, and qualitative agreement has been noticed. In spite of the model limitations, this study reveals the complexity of CFL development process, and it could be useful for better understanding relevant processes and phenomena in the microcirculation.

Hydrodynamic characteristics of a membrane oxygenator: modeling of pressure-flow characteristics and their influence on apparent viscosity

The viscosity obtained from pressure-flow characteristics of an oxygenator may help to detect factors that change oxygenator resistance. The objective of this study was to model pressure-flow characteristics of a membrane oxygenator with an integrated arterial filter and to quantify their influence on apparent viscosity of non-Newtonian fluids. One Newtonian fluid (glycerin solution) and two non-Newtonian fluids (whole bovine blood and a human red blood cell suspension) were perfused through an oxygenator and their pressure-flow characteristics examined systematically. Four resistance parameters for the pressure gradient characteristics approximation equation were obtained by the least squares method from the relational expression of pressure-flow characteristics and viscosity. For all three fluids, a non-linear flow to pressure change was observed with a coefficient of determination of almost 1 by exponential approximation. The glycerin solution had a higher pressure gradient (10-70%) than the other fluids; the apparent viscosity of the non-Newtonian fluids was around 35% lower than the static one measured by a torsional oscillation viscometer. Overall, our study demonstrated that the influence on the apparent viscosity of non-Newtonian fluids can be quantified by pressure gradient differences in a membrane oxygenator with an integrated arterial filter.

Two-dimensional simulation of red blood cell motion near a wall under a lateral force

The motion of a red blood cell suspended in a linear shear flow adjacent to a fixed boundary subject to an applied lateral force directed toward the boundary is simulated. A two-dimensional model is used that represents the viscous and elastic properties of normal red blood cells. Shear rates in the range of 100 to 600 s^{-1} are considered, and the suspending medium viscosity is 1 cP. In the absence of a lateral force, the cell executes a tumbling motion. With increasing lateral force, a transition from tumbling to tank-treading is predicted. The minimum force required to ensure tank-treading increases nonlinearly with the shear rate. Transient swinging motions occur when the force is slightly larger than the transition value. The applied lateral force is balanced by a hydrodynamic lift force resulting from the positive orientation of the long axis of the cell with respect to the wall. In the case of cyclic tumbling motions, the orientation angle takes positive values through most of the cycle, resulting in lift generation. These results are used to predict the motion of a cell close to the outer edge of the cell-rich core region that is generated when blood flows in a narrow tube. In this case, the lateral force is generated by shear-induced dispersion, resulting from cell-cell interactions in a region with a concentration gradient. This force is estimated using previous data on shear-induced dispersion. The cell is predicted to execute tank-treading motions at normal physiological hematocrit levels, with the possibility of tumbling at lower hematocrit levels.

Multiscale Particle-Based Modeling of Flowing Platelets in Blood Plasma Using Dissipative Particle Dynamics and Coarse Grained Molecular Dynamics

We developed a multiscale particle-based model of platelets, to study the transport dynamics of shear stresses between the surrounding fluid and the platelet membrane. This model facilitates a more accurate prediction of the activation potential of platelets by viscous shear stresses - one of the major mechanisms leading to thrombus formation in cardiovascular diseases and in prosthetic cardiovascular devices. The interface of the model couples coarse-grained molecular dynamics (CGMD) with dissipative particle dynamics (DPD). The CGMD handles individual platelets while the DPD models the macroscopic transport of blood plasma in vessels. A hybrid force field is formulated for establishing a functional interface between the platelet membrane and the surrounding fluid, in which the microstructural changes of platelets may respond to the extracellular viscous shear stresses transferred to them. The interaction between the two systems preserves dynamic properties of the flowing platelets, such as the flipping motion. Using this multiscale particle-based approach, we have further studied the effects of the platelet elastic modulus by comparing the action of the flow-induced shear stresses on rigid and deformable platelet models. The results indicate that neglecting the platelet deformability may overestimate the stress on the platelet membrane, which in turn may lead to erroneous predictions of the platelet activation under viscous shear flow conditions. This particle-based fluid-structure interaction multiscale model offers for the first time a computationally feasible approach for simulating deformable platelets interacting with viscous blood flow, aimed at predicting flow induced platelet activation by using a highly resolved mapping of the stress distribution on the platelet membrane under dynamic flow conditions.

Particle margination and its implications on intravenous anticancer drug delivery

"Margination" refers to the movement of particles in flow toward the walls of a channel. The term was first coined in physiology for describing the behavior of white blood cells (WBCs) and platelets in blood flow. The margination of particles is desirable for anticancer drug delivery because it results in the close proximity of drug-carrying particles to the endothelium, where they can easily diffuse into cancerous tumors through the leaky vasculature. Understanding the fundamentals of margination may further lead to the rational design of particles and allow for more specific delivery of anticancer drugs into tumors, thereby increasing patient comfort during cancer treatment. This paper reviews existing theoretical and experimental studies that focus on understanding margination. Margination is a complex phenomenon that depends on the interplay between inertial, hydrodynamic, electrostatic, lift, van der Waals, and Brownian forces. Parameters that have been explored thus far include the particle size, shape, density, stiffness, shear rate, and the concentration and aggregation state of red blood cells (RBCs). Many studies suggested that there exists an optimal particle size for margination to occur, and that nonspherical particles tend to marginate better than spherical particles. There are, however, conflicting views on the effects of particle density, stiffness, shear rate, and RBCs. The limitations of using the adhesion of particles to the channel walls in order to quantify margination propensity are explained, and some outstanding questions for future research are highlighted.

Hydrodynamic lift of vesicles and red blood cells in flow--from Fåhræus & Lindqvist to microfluidic cell sorting

Hydrodynamic lift forces acting on cells and particles in fluid flow receive ongoing attention from medicine, mathematics, physics and engineering. The early findings of Fåhræus & Lindqvist on the viscosity change of blood with the diameter of capillaries motivated extensive studies both experimentally and theoretically to illuminate the underlying physics. We review this historical development that led to the discovery of the inertial and non-inertial lift forces and elucidate the origins of these forces that are still not entirely clear. Exploiting microfluidic techniques induced a tremendous amount of new insights especially into the more complex interactions between the flow field and deformable objects like vesicles or red blood cells. We trace the way from the investigation of single cell dynamics to the recent developments of microfluidic techniques for particle and cell sorting using hydrodynamic forces. Such continuous and label-free on-chip cell sorting devices promise to revolutionize medical analyses for personalized point-of-care diagnosis. We present the state-of-the-art of different hydrodynamic lift-based techniques and discuss their advantages and limitations.

White blood cell margination in microcirculation

Proper functioning of white blood cells is not possible without their ability to adhere to vascular endothelium, which may occur only if they are close enough to vessel walls. To facilitate the adhesion, white blood cells migrate toward the vessel walls in blood flow through a process called margination. The margination of white cells depends on a number of conditions including local hematocrit, flow rate, red blood cell aggregation, and the deformability of both red and white cells. To better understand the margination process of white blood cells, we employ mesoscopic hydrodynamic simulations of a three-dimensional model of blood flow, which has been previously shown to capture quantitatively realistic blood flow properties and rheology. The margination properties of white blood cells are studied for a wide range of hematocrit values and flow conditions. Efficient white blood cell margination is found in an intermediate range of hematocrit values of Ht ≈ 0.2-0.4 and at relatively low flow rates, characteristic of the venular part of microcirculation. In addition, aggregation interactions between red blood cells lead to enhanced white-blood-cell margination. This simulation study provides a quantitative description of the margination of white blood cells, and is also highly relevant for the margination of particles or cells of similar size such as circulating tumor cells.

Margination of micro- and nano-particles in blood flow and its effect on drug delivery

Drug delivery by micro- and nano-carriers enables controlled transport of pharmaceuticals to targeted sites. Even though carrier fabrication has made much progress recently, the delivery including controlled particle distribution and adhesion within the body remains a great challenge. The adhesion of carriers is strongly affected by their margination properties (migration toward walls) in the microvasculature. To investigate margination characteristics of carriers of different shapes and sizes and to elucidate the relevant physical mechanisms, we employ mesoscopic hydrodynamic simulations of blood flow. Particle margination is studied for a wide range of hematocrit values, vessel sizes, and flow rates, using two- and three-dimensional models. The simulations show that the margination properties of particles improve with increasing carrier size. Spherical particles yield slightly better margination than ellipsoidal carriers; however, ellipsoidal particles exhibit a slower rotational dynamics near a wall favoring their adhesion. In conclusion, micron-sized ellipsoidal particles are favorable for drug delivery in comparison with sub-micron spherical particles.

Size-Dependent Nanoparticle Margination and Adhesion Propensity in a Microchannel

Intravenous injection of nanoparticles as drug delivery vehicles is a common practice in clinical trials of therapeutic agents to target specific cancerous or pathogenic sites. The vascular flow dynamics of nanocarriers (NCs) in human microcapillaries play an important role in the ultimate efficacy of this drug delivery method. This article reports an experimental study of the effect of nanoparticle size on their margination and adhesion propensity in microfluidic channels of a half-elliptical cross section. Spherical polystyrene particles ranging in diameter from 60 to 970 nm were flown in the microchannels and individual particles adhered to either the top or bottom wall of the channel were imaged using fluorescence microscopy. When the number concentration of particles in the flow was kept constant, the percentage of nanoparticles adhered to the top wall increased with decreasing diameter (d), with the number of particles adhered to the top wall following a d−3 trend. When the volume concentration of particles in solution was kept constant, no discernible trend was found. This experimental finding is explained by the competition between the Brownian force promoting margination and repulsive particle–particle electrostatic forces retarding adhesion to the wall. The 970 nm particles were found to adhere to the bottom wall much more than to the top wall for each of the three physiologically relevant shear rates tested, revealing the effect of gravitational force on the large particles. These findings on the flow behavior of spherical nanoparticles in artificial microcapillaries provide further insight for the rational design of NCs for targeted cancer therapeutics.

Computational biorheology of human blood flow in health and disease

Hematologic disorders arising from infectious diseases, hereditary factors and environmental influences can lead to, and can be influenced by, significant changes in the shape, mechanical and physical properties of red blood cells (RBCs), and the biorheology of blood flow. Hence, modeling of hematologic disorders should take into account the multiphase nature of blood flow, especially in arterioles and capillaries. We present here an overview of a general computational framework based on dissipative particle dynamics (DPD) which has broad applicability in cell biophysics with implications for diagnostics, therapeutics and drug efficacy assessments for a wide variety of human diseases. This computational approach, validated by independent experimental results, is capable of modeling the biorheology of whole blood and its individual components during blood flow so as to investigate cell mechanistic processes in health and disease. DPD is a Lagrangian method that can be derived from systematic coarse-graining of molecular dynamics but can scale efficiently up to arterioles and can also be used to model RBCs down to the spectrin level. We start from experimental measurements of a single RBC to extract the relevant biophysical parameters, using single-cell measurements involving such methods as optical tweezers, atomic force microscopy and micropipette aspiration, and cell-population experiments involving microfluidic devices. We then use these validated RBC models to predict the biorheological behavior of whole blood in healthy or pathological states, and compare the simulations with experimental results involving apparent viscosity and other relevant parameters. While the approach discussed here is sufficiently general to address a broad spectrum of hematologic disorders including certain types of cancer, this paper specifically deals with results obtained using this computational framework for blood flow in malaria and sickle cell anemia.

Blood viscosity in microvessels: experiment and theory

The apparent viscosity of blood flowing through narrow glass tubes decreases strongly with decreasing tube diameter over the range from about 300 μm to about 10 μm. This phenomenon, known as the Fåhraeus-Lindqvist effect, occurs because blood is a concentrated suspension of deformable red blood cells with a typical dimension of about 8 μm. Most of the resistance to blood flow through the circulatory system resides in microvessels with diameters in this range. Apparent viscosity of blood in microvessels in vivo has been found to be significantly higher than in glass tubes with corresponding diameters. Here we review experimental observations of blood's apparent viscosity in vitro and in vivo, and progress towards a quantitative theoretical understanding of the mechanisms involved.

Multiscale modeling of blood flow: from single cells to blood rheology

Mesoscale simulations of blood flow, where the red blood cells are described as deformable closed shells with a membrane characterized by bending rigidity and stretching elasticity, have made much progress in recent years to predict the flow behavior of blood cells and other components in various flows. To numerically investigate blood flow and blood-related processes in complex geometries, a highly efficient simulation technique for the plasma and solutes is essential. In this review, we focus on the behavior of single and several cells in shear and microcapillary flows, the shear-thinning behavior of blood and its relation to the blood cell structure and interactions, margination of white blood cells and platelets, and modeling hematologic diseases and disorders. Comparisons of the simulation predictions with existing experimental results are made whenever possible, and generally very satisfactory agreement is obtained.

Blood flow in small tubes: quantifying the transition to the non-continuum regime

In small vessels blood is usually treated as a Newtonian fluid down to diameters of ~200 μm. We investigate the flow of red blood cell (RBC) suspensions driven through small tubes (diameters 10-150 μm) in the range marking the transition from arterioles and venules to the largest capillary vessels. The results of the simulations combined with previous simulations of uniform shear flow and experimental data show that for diameters less than ~100 μm the suspension's stress cannot be described as a continuum, even a heterogeneous one. We employ the dissipative particle dynamics (DPD) model, which has been successfully used to predict human blood bulk viscosity in homogeneous shear flow. In tube flow the cross-stream stress gradient induces an inhomogeneous distribution of RBCs featuring a centreline cell density peak, and a cell-free layer (CFL) next to the wall. For a neutrally buoyant suspension the imposed linear shear-stress distribution together with the differentiable velocity distribution allow the calculation of the local viscosity across the tube section. The viscosity across the section as a function of the strain rate is found to be essentially independent of tube size for the larger diameters and is determined by the local haematocrit (H) and shear rate. Other RBC properties such as asphericity, deformation, and cell-flow orientation exhibit similar dependence for the larger tube diameters. As the tube size decreases below ~100 μm in diameter, the viscosity in the central region departs from the large-tube similarity function of the shear rate, since H increases significantly towards the centreline. The dependence of shear stress on tube size, in addition to the expected local shear rate and local haematocrit, implies that blood flow in small tubes cannot be described as a heterogeneous continuum. Based on the analysis of the DPD simulations and on available experimental results, we propose a simple velocity-slip model that can be used in conjunction with continuum-based simulations.

Continuum- and particle-based modeling of shapes and dynamics of red blood cells in health and disease

We review recent advances in multiscale modeling of the mechanics of healthy and diseased red blood cells (RBCs), and blood flow in the microcirculation. We cover the traditional continuum-based methods but also particle-based methods used to model both the RBCs and the blood plasma. We highlight examples of successful simulations of blood flow including malaria and sickle cell anemia.

Spring-network-based model of a red blood cell for simulating mesoscopic blood flow

We developed a mechanical model of a red blood cell (RBC) that is capable of expressing its characteristic behaviors in shear flows. The RBC was modeled as a closed shell membrane consisting of spring networks in the framework of the energy minimum concept. The fluid forces acting on RBCs were modeled from Newton's viscosity law and the conservation of momentum. In a steady shear flow, the RBC model exhibited various behaviors, depending on the shear rate; it tumbled, tank-treaded, or both. The transition from tumbling to tank-treading occurred at a shear rate of 20 s( - 1). The simulation of an RBC in steady and unsteady parallel shear flows (Couette flows) showed that the deformation parameters of the RBC were consistent with experimental results. The RBC in Poiseuille flow migrated radially towards the central axis of the flow channel. Axial migration became faster with an increase in the viscosity of the media, qualitatively consistent with experimental results. These results demonstrate that the proposed model satisfies the essential conditions for simulating RBC behavior in blood flow. Finally, a large-scale RBC flow simulation was implemented to show the capability of the proposed model for analyzing the mesoscopic nature of blood flow.

Two-phase model for prediction of cell-free layer width in blood flow

This study aimed to develop a numerical model capable of predicting changes in the cell-free layer (CFL) width in narrow tubes with consideration of red blood cell aggregation effects. The model development integrates to empirical relations for relative viscosity (ratio of apparent viscosity to medium viscosity) and core viscosity measured on independent blood samples to create a continuum model that includes these two regions. The constitutive relations were derived from in vitro experiments performed with three different glass-capillary tubes (inner diameter=30, 50 and 100 μm) over a wide range of pseudoshear rates (5-300 s(-1)). The aggregation tendency of the blood samples was also varied by adding Dextran 500 kDa. Our model predicted that the CFL width was strongly modulated by the relative viscosity function. Aggregation increased the width of CFL, and this effect became more pronounced at low shear rates. The CFL widths predicted in the present study at high shear conditions were in agreement with those reported in previous studies. However, unlike previous multi-particle models, our model did not require a high computing cost, and it was capable of reproducing results for a thicker CFL width at low shear conditions, depending on aggregating tendency of the blood.

Quantification of red blood cell deformation at high-hematocrit blood flow in microvessels

The deformation of red blood cells in microvessels was investigated numerically for various vessel diameters, hematocrits, and shear rates. We simulated blood flow in circular channels with diameters ranging from 9 to 50 μm, hematocrits from 20% to 45%, and shear rates from 20 to 150 s(-1) using a particle-based model with parallel computing. The apparent viscosity predicted by the simulation was in good agreement with previous experimental results. We quantified the deformation of red blood cells as a function of radial position. The numerical results demonstrated that because of the shape transition in response to local shear stress and the wall effect, the radial variation of red blood cell deformation in relatively large microvessels could be classified into three different regions: near-center, middle, and near-wall regions. Effects of the local shear stress and wall varied with vessel diameter, hematocrit, and shear rate.

Motion of red blood cells near microvessel walls: effects of a porous wall layer

A two-dimensional model is used to simulate the motion and deformation of a single mammalian red blood cell (RBC) flowing close to the wall of a microvessel, taking into account the effects of a porous endothelial surface layer (ESL) lining the vessel wall. Migration of RBCs away from the wall leads to the formation of a cell-depleted layer near the wall, which has a large effect on the resistance to blood flow in microvessels. The objective is to examine the mechanical factors causing this migration, including the effects of the ESL. The vessel is represented as a straight parallel-sided channel. The RBC is represented as a set of interconnected viscoelastic elements, suspended in plasma, a Newtonian fluid. The ESL is represented as a porous medium, and plasma flow in the layer is computed using the Brinkman approximation. It is shown that an initially circular cell positioned close to the ESL in a shear flow is deformed into an asymmetric shape. This breaking of symmetry leads to migration away from the wall. With increasing hydraulic resistivity of the layer, the rate of lateral migration increases. It is concluded that mechanical interactions of RBCs flowing in microvessels with a porous wall layer may reduce the rate of lateral migration and hence reduce the width of the cell-depleted zone external to the ESL, relative to the cell-depleted zone that would be formed if the interface between the ESL and free-flowing plasma were replaced by an impermeable boundary.

The effect of red blood cell aggregation on velocity and cell-depleted layer characteristics of blood in a bifurcating microchannel

Red blood cell (RBC) aggregation is a multifaceted phenomenon, and whether it is generally beneficial or deleterious remains unclear. In order to better understand its effect on microvascular blood flow, the phenomenon must be studied in complex geometries, as it is strongly dependent on time, flow, and geometry. The cell-depleted layer (CDL) which forms at the walls of microvessels has been observed to be enhanced by aggregation; however, details of the characteristics of the CDL in complex regions, such as bifurcations, require further investigation. In the present study, a microchannel with a T-junction was used to analyze the influence of aggregation on the flow field and the CDL. Micro-PIV using RBCs as tracers provided high resolution cell velocity data. CDL characteristics were measured from the same data using a newly developed technique based on motion detection. Skewed and sharpened velocity profiles in the daughter branches were observed, contrary to the behavior of a continuous Newtonian fluid. RBC aggregation was observed to increase the skewness, but decrease the sharpening, of the velocity profiles in the daughter branches. The CDL width was found to be significantly greater, with a wider distribution, in the presence of aggregation and the mean width increased proportionally with the reciprocal of the fraction of flow entering the daughter branch. Aggregation also significantly increased the roughness of the interface between the CDL and the RBC core. The present results provide further insight into how RBC aggregation may affect the flow in complex geometries, which is of importance in both understanding its functions invivo, and utilizing it as a tool in microfluidic devices.

A microfluidics approach towards high-throughput pathogen removal from blood using margination

Sepsis is an adverse systemic inflammatory response caused by microbial infection in blood. This paper reports a simple microfluidic approach for intrinsic, non-specific removal of both microbes and inflammatory cellular components (platelets and leukocytes) from whole blood, inspired by the invivo phenomenon of leukocyte margination. As blood flows through a narrow microchannel (20 × 20 µm), deformable red blood cells (RBCs) migrate axially to the channel centre, resulting in margination of other cell types (bacteria, platelets, and leukocytes) towards the channel sides. By using a simple cascaded channel design, the blood samples undergo a 2-stage bacteria removal in a single pass through the device, thereby allowing higher bacterial removal efficiency. As an application for sepsis treatment, we demonstrated separation of Escherichia coli and Saccharomyces cerevisiae spiked into whole blood, achieving high removal efficiencies of ∼80% and ∼90%, respectively. Inflammatory cellular components were also depleted by >80% in the filtered blood samples which could help to modulate the host inflammatory response and potentially serve as a blood cleansing method for sepsis treatment. The developed technique offers significant advantages including high throughput (∼1 ml/h per channel) and label-free separation which allows non-specific removal of any blood-borne pathogens (bacteria and fungi). The continuous processing and collection mode could potentially enable the return of filtered blood back to the patient directly, similar to a simple and complete dialysis circuit setup. Lastly, we designed and tested a larger filtration device consisting of 6 channels in parallel (∼6 ml/h) and obtained similar filtration performances. Further multiplexing is possible by increasing channel parallelization or device stacking to achieve higher throughput comparable to convectional blood dialysis systems used in clinical settings.

Blood-plasma separation in Y-shaped bifurcating microfluidic channels: a dissipative particle dynamics simulation study

The motion of a suspension of red blood cells (RBCs) flowing in a Y-shaped bifurcating microfluidic channel is investigated using a validated low-dimensional RBC model based on dissipative particle dynamics. Specifically, the RBC is represented as a closed torus-like ring of ten colloidal particles, which leads to efficient simulations of blood flow in microcirculation over a wide range of hematocrits. Adaptive no-slip wall boundary conditions were implemented to model hydrodynamic flow within a specific wall structure of diverging three-dimensional microfluidic channels, paying attention to controlling density fluctuations. Plasma skimming and the all-or-nothing phenomenon of RBCs in a bifurcating microfluidic channel have been investigated in our simulations for healthy and diseased blood, including the size of a cell-free layer on the daughter branches. The feed hematocrit level in the parent channel has considerable influence on blood-plasma separation. Compared to the blood-plasma separation efficiencies of healthy RBCs, malaria-infected stiff RBCs (iRBCs) have a tendency to travel into the low flow-rate daughter branch because of their different initial distribution in the parent channel. Our simulation results are consistent with previously published experimental results and theoretical predictions.