Inversion of hematocrit partition at microfluidic bifurcations

Partitioning of red blood cells (RBCs) at the level of bifurcations in the microcirculatory system affects many physiological functions yet it remains poorly understood. We address this problem by using T-shaped microfluidic bifurcations as a model. Our computer simulations and in vitro experiments reveal that the hematocrit ($\phi_0$) partition depends strongly on RBC deformability, as long as $\phi_0 <20$% (within the normal range in microcirculation), and can even lead to complete deprivation of RBCs in a child branch. Furthermore, we discover a deviation from the Zweifach-Fung effect which states that the child branch with lower flow rate recruits less RBCs than the higher flow rate child branch. At small enough $\phi_0$, we get the inverse scenario, and the hematocrit in the lower flow rate child branch is even higher than in the parent vessel. We explain this result by an intricate up-stream RBC organization and we highlight the extreme dependence of RBC transport on geometrical and cell mechanical properties. These parameters can lead to unexpected behaviors with consequences on the microcirculatory function and oxygen delivery in healthy and pathological conditions.Comment: 16 page

A simplified particulate model for coarse-grained hemodynamics simulations

Human blood flow is a multi-scale problem: in first approximation, blood is a dense suspension of plasma and deformable red cells. Physiological vessel diameters range from about one to thousands of cell radii. Current computational models either involve a homogeneous fluid and cannot track particulate effects or describe a relatively small number of cells with high resolution, but are incapable to reach relevant time and length scales. Our approach is to simplify much further than existing particulate models. We combine well established methods from other areas of physics in order to find the essential ingredients for a minimalist description that still recovers hemorheology. These ingredients are a lattice Boltzmann method describing rigid particle suspensions to account for hydrodynamic long range interactions and---in order to describe the more complex short-range behavior of cells---anisotropic model potentials known from molecular dynamics simulations. Paying detailedness, we achieve an efficient and scalable implementation which is crucial for our ultimate goal: establishing a link between the collective behavior of millions of cells and the macroscopic properties of blood in realistic flow situations. In this paper we present our model and demonstrate its applicability to conditions typical for the microvasculature.Comment: 12 pages, 11 figure

Blood flow in microvascular networks

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.Simulation of blood presents a very complex haemodynamics problem especially in relation to the understanding of atherogenesis. In many simulations, blood has been treated as a single-phase homogeneous fluid, a classical approach that does not account for the presence of red blood cells (RBCs). Although this approach provides satisfactory tools to describe certain aspects of blood flow in large arteries, it fails to give an adequate representation of the flow field in the vessels of smaller diameter where the size of the RBC becomes significant relative to vessel diameter. So, this article is concerned with the study of non-Newtonian blood flow in microvascular networks with the intention of developing a new cell depletion layer model to represent the behaviour of RBCs through bifurcating networks. The model is tested in a microvascular network constructed possessing realistic bifurcation features, with controlled dimensions and angles. The RBC depletion model treats blood as two continuum layers, with a central, non-Newtonian core region of concentrated red cell suspension that is surrounded by a layer of plasma (Newtonian fluid) adjacent to the vessel wall. In the central core region, blood is described by Quemada's non-Newtonian rheological model. Geometry differences are shown to significantly affect flow rates, haematocrit distributions and the corresponding cell depletion layers

Red blood cell lingering modulates hematocrit distribution in the microcirculation

The distribution of red blood cells (RBCs) in the microcirculation determines the oxygen delivery and solute transport to tissues. This process relies on the partitioning of RBCs at successive bifurcations throughout the microvascular network, and it has been known since the last century that RBCs partition disproportionately to the fractional blood flow rate, therefore leading to heterogeneity of the hematocrit (i.e., volume fraction of RBCs in blood) in microvessels. Usually, downstream of a microvascular bifurcation, the vessel branch with a higher fraction of blood flow receives an even higher fraction of RBC flux. However, both temporal and time-average deviations from this phase-separation law have been observed in recent studies. Here, we quantify how the microscopic behavior of RBC lingering (i.e., RBCs temporarily residing near the bifurcation apex with diminished velocity) influences their partitioning, through combined in vivo experiments and in silico simulations. We developed an approach to quantify the cell lingering at highly confined capillary-level bifurcations and demonstrate that it correlates with deviations of the phase-separation process from established empirical predictions by Pries et al. Furthermore, we shed light on how the bifurcation geometry and cell membrane rigidity can affect the lingering behavior of RBCs; e.g., rigid cells tend to linger less than softer ones. Taken together, RBC lingering is an important mechanism that should be considered when studying how abnormal RBC rigidity in diseases such as malaria and sickle-cell disease could hinder the microcirculatory blood flow or how the vascular networks are altered under pathological conditions (e.g., thrombosis, tumors, aneurysm)

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

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.National Institutes of Health (U.S.)Singapore-MIT Alliance for Research and Technology (SMART)United States. Dept. of Energy. Collaboratory on Mathematics for Mesoscopic Modeling of MaterialsUnited States. Dept. of Energy (INCITE Award

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 3D 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 one micrometer, this central cell-free layer persists for up to 100 $\mu$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 mouce venules. At bifurcations, however, a similar effect does not appear and margination is largely unaffected by the geometry. This anti-margination towards the vessel center after a confluence may explain in vivo observations by Woldhuis et al. [Am. J. Physiol. 262, H1217 (1992)] where 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