A low-dimensional model for the red blood cell

Unresolved RBCs: An upscaling strategy for the CFD-DEM simulation of blood flow with deformable cells

Numerical simulation of blood flow is a challenging topic due to the multiphase nature of this biological fluid. The choice of a specific method among the ones available in literature is often motivated by the physical scale of interest. Single-phase approximation allows for lower computational time, but does not consider this multiphase nature. Cell-level simulation, on the other hand, requires high computational resources and is limited to small scales. This work proposes a scale-up approach for cell-level simulation of blood flow, in the framework of unresolved CFD-DEM technique. This method offers the possibility to simulate hundreds of thousands of particles with limited computational effort, but requires specific models for fluid-particle interactions. Regarding blood flow, drag and lift force acting on the red blood cells (RBCs) are responsible for several macroscopic blood characteristics. Despite several correlations available for drag and lift force acting on rigid particles, specific force models for the simulation of deformable particles compatible with RBCs physics are missing. This study employs data obtained from cell-level simulations to derive equations then used in unresolved simulation of RBCs. The strategy followed during the modeling phase is presented, together with the model verification and validation. This approach returns satisfying results when used to simulate blood flow in large-scale channels. Up to half a million RBCs are considered, and computational effort is reported to allow a comparison with other existing methods. Future perspectives include further improvement of the model, such as a deeper understanding of particle-particle interactions.

Dynamic cellular responses to gravitational forces: Exploring the impact on white blood cell(s)

In recent years, the allure of space exploration and human spaceflight has surged, yet the effects of microgravity on the human body remain a significant concern. Immune and red blood cells rely on hematic or lymphatic streams as their primary means of transportation, posing notable challenges under microgravity conditions. This study sheds light on the intricate dynamics of cell behavior when suspended in bio-fluid under varying gravitational forces. Utilizing the dissipative particle dynamics approach, blood and white blood cells were modeled, with gravity applied as an external force along the vertical axis, ranging from 0 to 2 g in parameter sweeps. The results revealed discernible alterations in the cell shape and spatial alignment in response to gravity, quantified through metrics such as elongation and deformation indices, pitch angle, and normalized center of mass. Statistical analysis using the Mann-Whitney U test underscored clear distinctions between microgravity (<1 g) and hypergravity (>1 g) samples compared to normal gravity (1 g). Furthermore, the examination of forces exerted on the solid, including drag, shear stress, and solid forces, unveiled a reduction in the magnitude as the gravitational force increased. Additional analysis through dimensionless numbers unveiled the dominance of capillary and gravitational forces, which impacted cell velocity, leading to closer proximity to the wall and heightened viscous interaction with surrounding fluid particles. These interactions prompted shape alterations and reduced white blood cell area while increasing red blood cells. This study represents an effort in comprehending the effects of gravity on blood cells, offering insights into the intricate interplay between cellular dynamics and gravitational forces.

Numerical investigation on red blood cell flow based on unstructured grid

Prediction of blood cell flow is known as the difficult research by reason of the complexity of blood vessel. In this study, considering the complex structure of blood vessels, a mechanical model for red blood cell (RBC) based on unstructured grid has been established to study the flow characteristics of RBCs in complex blood vessels. In the model, the strain-energy function by Skalak is employed to model the shear elasticity and surface-area conservation of the membrane, and the hinge spring is used to describe the forces originating from local bending of the membrane. The immersed boundary method is utilized to couple the interphase force. Using the model, the stretching test of RBC is compared with the experiment data, and the good agreement verified the validation of the present model. The morphology of red blood cell and the blood viscosity in micro-vessel are studied. RBCs move with a symmetric shape (parachute shape) in small blood vessels, and the buckling instability is observed when the RBC flow slowly through a micro-vessel or a converging-diverging capillary. When the vessel diameter is around 10 μm, the reverse Fahraeus-Lindqvist effect is presented. The blood apparent viscosity shows linear increase with the blood hematocrit. In addition, Malaria infection can make the RBC deformability decreased and the blood apparent viscosity increased.

Cell-scale hemolysis evaluation of intervenient ventricular assist device based on dissipative particle dynamics

Most of the existing hemolysis mechanism studies are carried out on the macro flow scale. They assume that the erythrocyte membranes with different loads will suffer the same damage, which obviously has limitations. Thus, exploring the hemolysis mechanism through the macroscopic flow field information is a tough challenge. In order to further understand the non-physiological shear hemolysis phenomenon at the cell scale, this study used the coarse-grained erythrocytes damage model at the mesoscopic scale based on the transport dissipative particle dynamics (tDPD) method. Combined with computational fluid dynamics the hemolysis of scalarized shear stress (τ) in the clearance of "Impella 5.0" was evaluated under the Lagrange perspective and Euler perspective. The results from the Lagrange perspective showed that the change rate of scaled shear stress (τ˙) was the most critical factor in damaging RBCs in the rotor region of "Impella 5.0"and other transvalvular micro-axial blood pumps. Then, we propose a dimensionless number Dk with time integration based on τ˙ to evaluate hemolysis. The Dissipative particle dynamics simulation results are consistent with the Dk evaluation results, so τ˙ may be an important factor in the hemolysis of VADs. Finally, we tested the hemolysis of 30% hematocrit whole blood in the "Impella 5.0" shroud clearance from the Euler perspective. Relevant results indicate that because of the wall effect, the RBCs near the impeller side are more prone to damage, and most of the cytoplasm is also gathered at the rotor side.

Hemolysis prediction in bio-microfluidic applications using resolved CFD-DEM simulations

BACKGROUND AND OBJECTIVE

Hemolysis, namely hemoglobin leakage from red blood cells (RBCs), is one of the major sources of incorrect results in clinical tests, especially when passive microfluidics is involved. This is due to small characteristic dimensions which could cause strong RBCs deformation. Prediction of hemolysis is essential in the design and optimization of lab-on-a-chip devices for cell sorting and plasma separation. The aim of this work is to provide a numerical simulation tool this purpose applicable to real-scale bio-microfluidic devices with affordable computational cost.

METHODS

Blood is modelled as a suspension of biological cells, mainly RBCs, in liquid plasma assumed as a Newtonian, incompressible carrier fluid. Therefore, the physics of cells and carrier fluid is coupled by means of an immersed boundary concept known as resolved CFD-DEM. In this approach, the Navier-Stokes equations are numerically solved through a finite volume method with an additional penalty term to account for the presence of RBCs. RBCs' positions and velocities are updated by solving Newton and Euler equations for conservation of linear and angular momentum. To model the RBCs deformation, a reduced-order model is employed, where each RBC is represented by a clump of overlapping rigid spheres connected by fictional numerical bonds, whose properties are tuned to reproduce the ones of RBCs viscoelastic membrane. This coupled approach allows access to cell-level information and facilitates the usage of strain-based hemolysis models.

RESULTS

Different micro-channel geometries and blood hematocrits are simulated, to explore the influence of these factors on RBCs damage. Statistical analysis is performed to extract relevant biophysical quantities from numerical simulations such as hemolysis index distribution at the channel exit. Finally, the effect of carrier fluid viscosity is studied in relation to cell-cell interactions.

CONCLUSIONS

Simulation results show that hemolysis occurrence is almost independent of the hematocrit values in the microchannel, implying the possibility to speed up calculation using low hematocrit values. Nevertheless, using whole blood viscosity for the carrier fluid overestimates the value of the hemolysis index by almost one order of magnitude.

Morphology of depletant-induced erythrocyte aggregates

Red blood cells suspended in quiescent plasma tend to aggregate into multicellular assemblages, including linearly stacked columnar rouleaux, which can reversibly form more complex clusters or branching networks. While these aggregates play an essential role in establishing hemorheological and pathological properties, the biophysics behind their self-assembly into dynamic mesoscopic structures remains under-explored. We employ coarse-grained molecular simulations to model low-hematocrit erythrocytes subject to short-range implicit depletion forces, and demonstrate not only that depletion interactions are sufficient to account for a sudden dispersion-aggregate transition, but also that the volume fraction of depletant macromolecules controls small aggregate morphology. We observe a sudden transition from a dispersion to a linear column rouleau, followed by a slow emergence of disorderly amorphous clusters of many short rouleaux at larger volume fractions. This work demonstrates how discocyte topology and short-range, non-specific, physical interactions are sufficient to self-assemble erythrocytes into various aggregate structures, with markedly different morphologies and biomedical consequences.

The dynamics of a healthy and infected red blood cell in flow through constricted channels: A DPD simulation

Understanding the dynamics of red blood cell (RBC) motion under in silico conditions is central to the development of cost-effective diagnostic tools. Specifically, unraveling the relationship between the rheological properties and the nature of shape change in the RBC (healthy or infected) can be extremely useful. In case of malarial infection, RBC progressively loses its deformability and tends to occlude the microvessel. In the present study, detailed mesoscopic simulations are performed to investigate the deformation dynamics of an RBC in flow through a constricted channel. Specifically, the manifestation of viscous forces (through flow rates) on the passage and blockage characteristics of a healthy red blood cell (hRBC) vis-á-vis an infected red blood cell (iRBC) are investigated. A finite-sized dissipative particle dynamics framework is used to model plasma in conjunction with a discrete model for the RBC. Instantaneous wall boundary method was used to model no-slip wall boundary conditions with a good control on the near-wall density fluctuations and compressibility effects. To investigate the microvascular occlusion, the RBC motion through 2 types of constricted channels, viz, (1) a tapered microchannel and (2) a stenosed-type microchannel, were simulated. It was observed that the deformation of an infected cell was much less compared with a healthy cell, with an attendant increase in the passage time. Apart from the qualitative features, deformation indices were obtained. The deformation of hRBC was sudden, while the iRBC deformed slowly as it traversed through the constriction. For higher flow rates, both hRBC and iRBC were found to undergo severe deformation. Even under low flow rates, hRBC could easily traverse past the constricted channel. However, for sufficiently slow flow rates (eg, capillary flows), the microchannel was found to be completely blocked by the iRBC.

Synergistic Integration of Laboratory and Numerical Approaches in Studies of the Biomechanics of Diseased Red Blood Cells

In red blood cell (RBC) disorders, such as sickle cell disease, hereditary spherocytosis, and diabetes, alterations to the size and shape of RBCs due to either mutations of RBC proteins or changes to the extracellular environment, lead to compromised cell deformability, impaired cell stability, and increased propensity to aggregate. Numerous laboratory approaches have been implemented to elucidate the pathogenesis of RBC disorders. Concurrently, computational RBC models have been developed to simulate the dynamics of RBCs under physiological and pathological conditions. In this work, we review recent laboratory and computational studies of disordered RBCs. Distinguished from previous reviews, we emphasize how experimental techniques and computational modeling can be synergically integrated to improve the understanding of the pathophysiology of hematological disorders.

Computational Biomechanics of Human Red Blood Cells in Hematological Disorders

We review recent advances in multiscale modeling of the biomechanical characteristics of red blood cells (RBCs) in hematological diseases, and their relevance to the structure and dynamics of defective RBCs. We highlight examples of successful simulations of blood disorders including malaria and other hereditary disorders, such as sickle-cell anemia, spherocytosis, and elliptocytosis.

Cellular Level In-silico Modeling of Blood Rheology with An Improved Material Model for Red Blood Cells

Many of the intriguing properties of blood originate from its cellular nature. Therefore, accurate modeling of blood flow related phenomena requires a description of the dynamics at the level of individual cells. This, however, presents several computational challenges that can only be addressed by high performance computing. We present Hemocell, a parallel computing framework which implements validated mechanical models for red blood cells and is capable of reproducing the emergent transport characteristics of such a complex cellular system. It is computationally capable of handling large domain sizes, thus it is able to bridge the cell-based micro-scale and macroscopic domains. We introduce a new material model for resolving the mechanical responses of red blood cell membranes under various flow conditions and compare it with a well established model. Our new constitutive model has similar accuracy under relaxed flow conditions, however, it performs better for shear rates over 1,500 s^-1. We also introduce a new method to generate randomized initial conditions for dense mixtures of different cell types free of initial positioning artifacts.

Hybrid smoothed dissipative particle dynamics and immersed boundary method for simulation of red blood cells in flows

In biofluid flow systems, often the flow problems of fluids of complex structures, such as the flow of red blood cells (RBCs) through complex capillary vessels, need to be considered. The smoothed dissipative particle dynamics (SDPD), a particle-based method, is one of the easy and flexible methods to model such complex structure fluids. It couples the best features of the smoothed particle hydrodynamics (SPH) and dissipative particle dynamics (DPD), with parameters having specific physical meaning (coming from SPH discretization of the Navier-Stokes equations), combined with thermal fluctuations in a mesoscale simulation, in a similar manner to the DPD. On the other hand, the immersed boundary method (IBM), a preferred method for handling fluid-structure interaction problems, has also been widely used to handle the fluid-RBC interaction in RBC simulations. In this paper, we aim to couple SDPD and IBM together to carry out the simulations of RBCs in complex flow problems. First, we develop the SDPD-IBM model in details, including the SDPD model for the evolving fluid flow, the RBC model for calculating RBC deformation force, the IBM for treating fluid-RBC interaction, and the solid boundary treatment model as well. We then conduct the verification and validation of the combined SDPD-IBM method. Finally, we demonstrate the capability of the SDPD-IBM method by simulating the flows of RBCs in rectangular, cylinder, curved, bifurcated, and constricted tubes, respectively.

Theory and algorithms to compute Helfrich bending forces: a review

Cell membranes are vital to shield a cell's interior from the environment. At the same time they determine to a large extent the cell's mechanical resistance to external forces. In recent years there has been considerable interest in the accurate computational modeling of such membranes, driven mainly by the amazing variety of shapes that red blood cells and model systems such as vesicles can assume in external flows. Given that the typical height of a membrane is only a few nanometers while the surface of the cell extends over many micrometers, physical modeling approaches mostly consider the interface as a two-dimensional elastic continuum. Here we review recent modeling efforts focusing on one of the computationally most intricate components, namely the membrane's bending resistance. We start with a short background on the most widely used bending model due to Helfrich. While the Helfrich bending energy by itself is an extremely simple model equation, the computation of the resulting forces is far from trivial. At the heart of these difficulties lies the fact that the forces involve second order derivatives of the local surface curvature which by itself is the second derivative of the membrane geometry. We systematically derive and compare the different routes to obtain bending forces from the Helfrich energy, namely the variational approach and the thin-shell theory. While both routes lead to mathematically identical expressions, so-called linear bending models are shown to reproduce only the leading order term while higher orders differ. The main part of the review contains a description of various computational strategies which we classify into three categories: the force, the strong and the weak formulation. We finally give some examples for the application of these strategies in actual simulations.

MD/DPD Multiscale Framework for Predicting Morphology and Stresses of Red Blood Cells in Health and Disease

Healthy red blood cells (RBCs) have remarkable deformability, squeezing through narrow capillaries as small as 3 microns in diameter without any damage. However, in many hematological disorders the spectrin network and lipid bilayer of diseased RBCs may be significantly altered, leading to impaired functionality including loss of deformability. We employ a two-component whole-cell multiscale model to quantify the biomechanical characteristics of the healthy and diseased RBCs, including Plasmodium falciparum-infected RBCs (Pf-RBCs) and defective RBCs in hereditary disorders, such as spherocytosis and elliptocytosis. In particular, we develop a two-step multiscale framework based on coarse-grained molecular dynamics (CGMD) and dissipative particle dynamics (DPD) to predict the static and dynamic responses of RBCs subject to tensile forcing, using experimental information only on the structural defects in the lipid bilayer, cytoskeleton, and their interaction. We first employ CGMD on a small RBC patch to compute the shear modulus, bending stiffness, and network parameters, which are subsequently used as input to a whole-cell DPD model to predict the RBC shape and corresponding stress field. For Pf-RBCs at trophozoite and schizont stages, the presence of cytoadherent knobs elevates the shear response in the lipid bilayer and stiffens the RBC membrane. For RBCs in spherocytosis and elliptocytosis, the bilayer-cytoskeleton interaction is weakened, resulting in substantial increase of the tensile stress in the lipid bilayer. Furthermore, we investigate the transient behavior of stretching deformation and shape relaxation of the normal and defective RBCs. Different from the normal RBCs possessing high elasticity, our simulations reveal that the defective RBCs respond irreversibly, i.e., they lose their ability to recover the normal biconcave shape in successive loading cycles of stretching and relaxation. Our findings provide fundamental insights into the microstructure and biomechanics of RBCs, and demonstrate that the two-step multiscale framework presented here can be used effectively for in silico studies of hematological disorders based on first principles and patient-specific experimental input at the protein level.

Dynamic and rheological properties of soft biological cell suspensions

Quantifying dynamic and rheological properties of suspensions of soft biological particles such as vesicles, capsules, and red blood cells (RBCs) is fundamentally important in computational biology and biomedical engineering. In this review, recent studies on dynamic and rheological behavior of soft biological cell suspensions by computer simulations are presented, considering both unbounded and confined shear flow. Furthermore, the hemodynamic and hemorheological characteristics of RBCs in diseases such as malaria and sickle cell anemia are highlighted.

Flow structures and red blood cell dynamics in arteriole of dilated or constricted cross section

Vessel with 'circular' or 'star-shaped' cross sections are studied, representing respectively dilated or constricted cases where endothelial cells smoothly line or bulge into the lumen. Computational haemodynamics simulations are carried out on idealised periodic arteriole-sized vessels, with red blood cell 'tube' hematocrit value=24%. A further simulation of a single red blood cell serves for comparison purposes. The bulk motion of the red blood cells reproduces well-known effects, including the presence of a cell-free layer and the apparent shear-thinning non-Newtonian rheology. The velocity flow field is analysed in a Lagrangian reference frame, relative to any given red blood cell, hence removing the bulk coaxial motion and highlighting instead the complex secondary flow patterns. An aggregate formation becomes apparent, continuously rearranging and dynamic, brought about by the inter-cellular fluid mechanics interactions and the deformability properties of the cells. The secondary flow field induces a vacillating radial migration of the red blood cells. At different radial locations, the red blood cells express different residence times, orientation and shape. The shear stresses exerted by the flow on the vessel wall are influenced by the motion of red blood cells, despite the presence of the cell-free layer. Spatial (and temporal) variations of wall shear stress patters are observed, especially for the 'circular' vessel. The 'star-shaped' vessel bears considerable stress at the protruding endothelial cell crests, where the stress vectors are coaxially aligned. The bulging endothelial cells hence regularise the transmission of stresses on the vessel wall.

Hydrodynamic interactions of deformable polymeric nanocarriers and the effect of crosslinking

We report theoretical as well as numerical investigations of deformable nanocarriers (NCs) under physiologically relevant flow conditions. Specifically, to model the deformable lysozyme-core/dextran-shell crosslinked polymer based NC with internal nanostructure and subject it to external hydrodynamic shear, we have introduced a coarse-grained model for the NC and have adopted a Brownian dynamics framework, which incorporates hydrodynamic interactions, in order to describe the static and dynamic properties of the NC. In order to represent the fluidity of the polymer network in the dextran brush-like corona, we coarse-grain the structure of the NC based on the hypothesis that Brownian motion, polymer melt reptations, and crosslinking density dominate their structure and dynamics. In our model, we specify a crosslinking density and employ the simulated annealing protocol to mimic the experimental synthesis steps in order to obtain the appropriate internal structure of the core-shell polymer. We then compute the equilibrium as well as steady shear rheological properties as functions of the Péclet number and the crosslinking density, in the presence of hydrodynamic interactions. We find that with increasing crosslinking, the stiffness of the nanocarrier increases, the radius of gyration decreases, and as a consequence the self-diffusivity increases. The nanocarrier under shear deforms and orients along the direction of the applied shear and we find that the orientation and deformation under shear are dependent on the shear rate and the crosslinking density. We compare various dynamic properties of the NC as a function of the shear force, such as orientation, deformation, intrinsic stresses etc., with previously reported computational and experimental results of other model systems. The computational approach described here serves as a powerful tool for the rational design of NCs by taking both the physiological as well as the hydrodynamic environments into consideration. Development of such models is essential in order to gain useful insights that may be translated into the optimal design of NCs for diagnostic as well as targeted drug delivery applications.

Probing red blood cell mechanics, rheology and dynamics with a two-component multi-scale model

This study is partially motivated by the validation of a new two-component multi-scale cell model we developed recently that treats the lipid bilayer and the cytoskeleton as two distinct components. Here, the whole cell model is validated and compared against several available experiments that examine red blood cell (RBC) mechanics, rheology and dynamics. First, we investigated RBC deformability in a microfluidic channel with a very small cross-sectional area and quantified the mechanical properties of the RBC membrane. Second, we simulated twisting torque cytometry and compared predicted rheological properties of the RBC membrane with experimental measurements. Finally, we modelled the tank-treading (TT) motion of a RBC in a shear flow and explored the effect of channel width variation on the TT frequency. We also investigated the effects of bilayer-cytoskeletal interactions on these experiments and our simulations clearly indicated that they play key roles in the determination of cell membrane mechanical, rheological and dynamical properties. These simulations serve as validation tests and moreover reveal the capabilities and limitations of the new whole cell model.

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.

Temporal Multiscale Approach for Nanocarrier Motion with Simultaneous Adhesion and Hydrodynamic Interactions in Targeted Drug Delivery

We present a fluctuating hydrodynamics approach and a hybrid approach combining fluctuating hydrodynamics with generalized Langevin dynamics to resolve the motion of a nanocarrier when subject to both hydrodynamic interactions and adhesive interactions. Specifically, using these approaches, we compute equilibrium probability distributions at constant temperature as well as velocity autocorrelation functions of the nanocarrier subject to thermal motion in a quiescent Newtonian fluid medium, when tethered by a harmonic spring force mimicking a tether due to a single receptor-ligand bond. We demonstrate that the thermal equipartition of translation, rotation, and spring degrees of freedom are preserved by our formalism while simultaneously resolving the nature of the hydrodynamic correlations. Additionally, we evaluate the potential of mean force (or free energy density) along a specified reaction coordinate to faciltate extensive conformational sampling of the nanocarrier motion. We show that our results are in excellent agreement with analytical results and Monte Carlo simulations, thereby validating our methodologies. The frameworks we have presented provide a comprehensive platform for temporal multiscale modeling of hydrodynamic and microscopic interactions mediating nanocarrier motion and adhesion in vascular targeted drug delivery.

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.

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.

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.

Predicting the morphology of sickle red blood cells using coarse-grained models of intracellular aligned hemoglobin polymers

Sickle red blood cells (SS-RBCs) exhibit heterogeneous cell morphologies (sickle, holly leaf, granular, etc.) in the deoxygenated state due to the polymerization of the sickle hemoglobin. Experimental evidence points to a close relationship between SS-RBC morphology and intracellular aligned hemoglobin polymers. Here, we develop a coarse-grained (CG) stochastic model to represent the growth of the intracellular aligned hemoglobin polymer domain. The CG model is calibrated based on the mechanical properties (Young's modulus, bending rigidity) of the sickle hemoglobin fibers reported in experiments. The process of the cell membrane transition is simulated for physiologic aligned hemoglobin polymer configurations and mean corpuscular hemoglobin concentration. Typical SS-RBC morphologies observed in experiments can be obtained from the current model as a result of the intracellular aligned hemoglobin polymer development without introducing any further ad hoc assumptions. It is found that the final shape of SS-RBCs is primarily determined by the angular width of the aligned hemoglobin polymer domain, but it also depends, to a lesser degree, on the polymer growth rate and the cell membrane rigidity. Cell morphologies are quantified by structural shape factors, which agree well with experimental results from medical images.

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.

Predicting human blood viscosity in silico

The viscosity of blood has long been used as an indicator in the understanding and treatment of disease, and the advent of modern viscometers allows its measurement with ever-improving clinical convenience. However, these advances have not been matched by theoretical developments that can yield a quantitative understanding of blood's microrheology and its possible connection to relevant biomolecules (e.g., fibrinogen). Using coarse-grained molecular dynamics and two different red blood cell models, we accurately predict the dependence of blood viscosity on shear rate and hematocrit. We explicitly represent cell-cell interactions and identify the types and sizes of reversible rouleaux structures that yield a tremendous increase of blood viscosity at low shear rates. We also present the first quantitative estimates of the magnitude of adhesive forces between red cells. In addition, our simulations support the hypothesis, previously deduced from experiments, of yield stress as an indicator of cell aggregation. This non-Newtonian behavior is analyzed and related to the suspension's microstructure, deformation, and dynamics of single red blood cells. The most complex cell dynamics occurs in the intermediate shear rate regime, where individual cells experience severe deformation and transient folded conformations. The generality of these cell models together with single-cell measurements points to the future prediction of blood-viscosity anomalies and the corresponding microstructures associated with various diseases (e.g., malaria, AIDS, and diabetes mellitus). The models can easily be adapted to tune the properties of a much wider class of complex fluids including capsule and vesicle suspensions.

Predicting dynamics and rheology of blood flow: A comparative study of multiscale and low-dimensional models of red blood cells

We compare the predictive capability of two mathematical models for red blood cells (RBCs) focusing on blood flow in capillaries and arterioles. Both RBC models as well as their corresponding blood flows are based on the dissipative particle dynamics (DPD) method, a coarse-grained molecular dynamics approach. The first model employs a multiscale description of the RBC (MS-RBC), with its membrane represented by hundreds or even thousands of DPD-particles connected by springs into a triangular network in combination with out-of-plane elastic bending resistance. Extra dissipation within the network accounts for membrane viscosity, while the characteristic biconcave RBC shape is achieved by imposition of constraints for constant membrane area and constant cell volume. The second model is based on a low-dimensional description (LD-RBC) constructed as a closed torus-like ring of only 10 large DPD colloidal particles. They are connected into a ring by worm-like chain (WLC) springs combined with bending resistance. The LD-RBC model can be fitted to represent the entire range of nonlinear elastic deformations as measured by optical-tweezers for healthy and for infected RBCs in malaria. MS-RBCs suspensions model the dynamics and rheology of blood flow accurately for any vessel size but this approach is computationally expensive for vessel diameters above 100μm. Surprisingly, the much more economical suspensions of LD-RBCs also capture the blood flow dynamics and rheology accurately except for small-size vessels comparable to RBC diameter. In particular, the LD-RBC suspensions are shown to properly capture the experimental data for the apparent viscosity of blood and its cell-free layer (CFL) in tube flow. Taken together, these findings suggest a hierarchical approach in modeling blood flow in the arterial tree, whereby the MS-RBC model should be employed for capillaries and arterioles below 100μm, the LD-RBC model for arterioles, and the continuum description for arteries.