Tension of red blood cell membrane in simple shear flow

Effects of membrane viscoelasticity on the red blood cell dynamics in a microcapillary

The mechanical properties of red blood cells (RBCs) play key roles in their biological functions in microcirculation. In particular, RBCs must deform significantly to travel through microcapillaries with sizes comparable with or even smaller than their own. Although the dynamics of RBCs in microcapillaries have received considerable attention, the effect of membrane viscoelasticity has been largely overlooked. In this work, we present a computational study based on the boundary integral method and thin-shell mechanics to examine how membrane viscoelasticity influences the dynamics of RBCs flowing through straight and constricted microcapillaries. Our results reveal that the cell with a viscoelastic membrane undergoes substantially different motion and deformation compared with results based on a purely elastic membrane model. Comparisons with experimental data also suggest the importance of accounting for membrane viscoelasticity to properly capture the transient dynamics of an RBC flowing through a microcapillary. Taken together, these findings demonstrate the significant effects of membrane viscoelasticity on RBC dynamics in different microcapillary environments. The computational framework also lays the groundwork for more accurate quantitative modeling of the mechanical response of RBCs in their mechanotransduction process in subsequent investigations.

Structural Configuration of Blood Cell Membranes Determines Their Nonlinear Deformation Properties

The ability of neutrophils and red blood cells (RBCs) to undergo significant deformations is a key to their normal functioning. Disruptions of these processes can lead to pathologies. This work studied the influence of structural configuration rearrangements of membranes after exposure to external factors on the ability of native membranes of neutrophils and RBCs to undergo deep deformation. The rearrangement of the structural configuration of neutrophil and RBC membranes under the influence of cytological fixatives caused nonlinear deformation phenomena. There were an increase in Young's modulus, a decrease in the depth of homogeneous bending, and a change in the distance between cytoskeletal junctions. Based on the results of the analysis of experimental data, a mathematical model was proposed that describes the process of deep bending of RBСs and neutrophil membranes.

Nonlinear large deformation of a spherical red blood cell induced by ultrasonic standing wave

A computational model is developed to investigate the nonlinear static deformation of a spherical (osmotically swollen) red blood cell (RBC) induced by ultrasonic standing wave. The ultrasonic standing wave can generate steady acoustic radiation stress to deform the cell, and in turn, the deformed cell reshapes the acoustic field. This is a real-time coupling problem between the acoustic field and the mechanical field. In the computational model, the acoustic radiation stress acting on the RBC membrane is modeled by adopting the nonviscous momentum flux theory. The RBC membrane is modeled as a hyperelastic shell considering the in-plane elasticity, bending elasticity, and surface tension of the membrane. The volume conservation constraint of the membrane sealing fluid is applied to ensure the osmotic balance of the membrane. To address this real-time coupling problem, the computational model is implemented by a finite element method algorithm. The numerical results are compared with the existing theoretical model and experimental data, and the strain hardening trend of the experimental data is successfully predicted, which verifies the accuracy and effectiveness of the computational model. The computational model can accurately extract the mechanical properties of cells from acoustic deformation experiments, which is helpful for the diagnosis of some human diseases.

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.

Numerical Simulations of Red-Blood Cells in Fluid Flow: A Discrete Multiphysics Study

In this paper, we present a methodological study of modelling red blood cells (RBCs) in shear-induced flows based on the discrete multiphysics (DMP) approach. The DMP is an alternative approach from traditional multiphysics based on meshless particle-based methods. The proposed technique has been successful in modelling multiphysics and multi-phase problems with large interfacial deformations such as those in biological systems. In this study, we present the proposed method and introduce an accurate geometrical representation of the RBC. The results were validated against available data in the literature. We further illustrate that the proposed method is capable of modelling the rupture of the RBC membrane with minimum computational difficulty.

Numerical simulations of capsule deformation using a dual time-stepping lattice Boltzmann method

In this work a quasisteady, dual time-stepping lattice Boltzmann method is proposed for simulation of capsule deformation. At each time step the steady-state lattice Boltzmann equation is solved using the full approximation storage multigrid scheme for nonlinear equations. The capsule membrane is modeled as an infinitely thin shell suspended in an ambient fluid domain with the fluid structure interaction computed using the immersed boundary method. A finite element method is used to compute the elastic forces exerted by the capsule membrane. Results for a wide range of parameters and initial configurations are presented. The proposed method is found to reduce the computational time by a factor of ten.

Numerical Model for the Determination of Erythrocyte Mechanical Properties and Wall Shear Stress in vivo From Intravital Microscopy

The mechanical properties and deformability of Red Blood Cells (RBCs) are important determinants of blood rheology and microvascular hemodynamics. The objective of this study is to quantify the mechanical properties and wall shear stress experienced by the RBC membrane during capillary plug flow in vivo utilizing high speed video recording from intravital microscopy, biomechanical modeling, and computational methods. Capillaries were imaged in the rat cremaster muscle pre- and post-RBC transfusion of stored RBCs for 2-weeks. RBC membrane contours were extracted utilizing image processing and parametrized. RBC parameterizations were used to determine updated deformation gradient and Lagrangian Green strain tensors for each point along the parametrization and for each frame during plug flow. The updated Lagrangian Green strain and Displacement Gradient tensors were numerically fit to the Navier-Lame equations along the parameterized boundary to determined Lame's constants. Mechanical properties and wall shear stress were determined before and transfusion, were grouped in three populations of erythrocytes: native cells (NC) or circulating cells before transfusion, and two distinct population of cells after transfusion with stored cells (SC1 and SC2). The distinction, between the heterogeneous populations of cells present after the transfusion, SC1 and SC2, was obtained through principle component analysis (PCA) of the mechanical properties along the membrane. Cells with the first two principle components within 3 standard deviations of the mean, were labeled as SC1, and those with the first two principle components greater than 3 standard deviations from the mean were labeled as SC2. The calculated shear modulus average was 1.1±0.2, 0.90±0.15, and 12 ± 8 MPa for NC, SC1, and SC2, respectively. The calculated young's modulus average was 3.3±0.6, 2.6±0.4, and 32±20 MPa for NC, SC1, and SC2, respectively. o our knowledge, the methods presented here are the first estimation of the erythrocyte mechanical properties and shear stress in vivo during capillary plug flow. In summary, the methods introduced in this study may provide a new avenue of investigation of erythrocyte mechanics in the context of hematologic conditions that adversely affect erythrocyte mechanical properties.

Estimation of membrane bending modulus of stiffness tuned human red blood cells from micropore filtration studies

Human red blood cells (RBCs) need to deform in order to pass through capillaries in human vasculature with diameter smaller than that of the RBC. An altered RBC cell membrane stiffness (CMS), thereby, is likely to have consequences on their flow rate. RBC CMS is known to be affected by several commonly encountered disease conditions. This study was carried out to investigate whether an increase in RBC CMS, to the extent seen in such commonly encountered medical conditions, affects the RBC flow rate through channels with diameters comparable to that of the RBC. To do this, we use RBCs extracted from a healthy individual with no known medical conditions and treated with various concentrations of Bovine Serum Albumin (BSA). We study their flow through polycarbonate membranes with pores of diameter 5μm and 8μm which are smaller than and comparable to the RBC diameter respectively. The studies are carried out at constant hematocrit and volumetric flow rate. We find that when the diameter of the capillary is smaller than that of the RBC, the flow rate of the RBCs is lowered as the concentration of BSA is increased while the reverse is true when the diameter is comparable to that of the RBC. We confirm that this is a consequence of altered CMS of the RBCs from their reorientation dynamics in an Optical Tweezer. We find that a treatment with 0.50mg/ml BSA mimics the situation for RBCs extracted from a healthy individual while concentrations higher than 0.50mg/ml elevate the RBC CMS across a range expected for individuals with a condition of hyperglycemia. Using a simple theoretical model of the RBC deformation process at the entry of a narrow channel, we extract the RBC membrane bending modulus from their flow rate.

Deformation of a Red Blood Cell in a Narrow Rectangular Microchannel

The deformability of a red blood cell (RBC) is one of the most important biological parameters affecting blood flow, both in large arteries and in the microcirculation, and hence it can be used to quantify the cell state. Despite numerous studies on the mechanical properties of RBCs, including cell rigidity, much is still unknown about the relationship between deformability and the configuration of flowing cells, especially in a confined rectangular channel. Recent computer simulation techniques have successfully been used to investigate the detailed behavior of RBCs in a channel, but the dynamics of a translating RBC in a narrow rectangular microchannel have not yet been fully understood. In this study, we numerically investigated the behavior of RBCs flowing at different velocities in a narrow rectangular microchannel that mimicked a microfluidic device. The problem is characterized by the capillary number C a , which is the ratio between the fluid viscous force and the membrane elastic force. We found that confined RBCs in a narrow rectangular microchannel maintained a nearly unchanged biconcave shape at low C a , then assumed an asymmetrical slipper shape at moderate C a , and finally attained a symmetrical parachute shape at high C a . Once a RBC deformed into one of these shapes, it was maintained as the final stable configurations. Since the slipper shape was only found at moderate C a , measuring configurations of flowing cells will be helpful to quantify the cell state.

Passive swimming of a microcapsule in vertical fluid oscillation

The artificial microswimmer is a cutting-edge technology with applications in drug delivery and micro-total-analysis systems. The flow field around a microswimmer can be regarded as Stokes flow, in which reciprocal body deformation cannot induce migration. In this study, we propose a microcapsule swimmer that undergoes amoeboidlike shape deformations under fluid oscillation conditions. This is a study on the propulsion principle using a capsule with a solid membrane, and one of only a few studies using fluid oscillation. The microswimmer consists of an elastic capsule containing fluid and a rigid sphere. Opposing forces are generated when fluid oscillations are applied, because the densities of the internal fluid and sphere are different. The opposing forces induce nonreciprocal body deformation, which leads to migration of the microswimmer under Stokes flow conditions. Using numerical simulations, we found that the microswimmer propels itself in one of two modes, i.e., stroke swimming or drag swimming. We discuss the feasibility of the proposed microswimmer and show that the most efficient swimmer can migrate tens of micrometers per second. These findings pave the way for future artificial microswimmer designs.

Simulation of the nodal flow of mutant embryos with a small number of cilia: comparison of mechanosensing and vesicle transport hypotheses

Left-right (L-R) asymmetry in the body plan is determined by nodal flow in vertebrate embryos. Shinohara et al. (Shinohara K et al. 2012 Nat. Commun.3, 622 (doi:10.1038/ncomms1624)) used Dpcd and Rfx3 mutant mouse embryos and showed that only a few cilia were sufficient to achieve L-R asymmetry. However, the mechanism underlying the breaking of symmetry by such weak ciliary flow is unclear. Flow-mediated signals associated with the L-R asymmetric organogenesis have not been clarified, and two different hypotheses-vesicle transport and mechanosensing-are now debated in the research field of developmental biology. In this study, we developed a computational model of the node system reported by Shinohara et al. and examined the feasibilities of the two hypotheses with a small number of cilia. With the small number of rotating cilia, flow was induced locally and global strong flow was not observed in the node. Particles were then effectively transported only when they were close to the cilia, and particle transport was strongly dependent on the ciliary positions. Although the maximum wall shear rate was also influenced by ciliary position, the mean wall shear rate at the perinodal wall increased monotonically with the number of cilia. We also investigated the membrane tension of immotile cilia, which is relevant to the regulation of mechanotransduction. The results indicated that tension of about 0.1 μN m-1 was exerted at the base even when the fluid shear rate was applied at about 0.1 s-1. The area of high tension was also localized at the upstream side, and negative tension appeared at the downstream side. Such localization may be useful to sense the flow direction at the periphery, as time-averaged anticlockwise circulation was induced in the node by rotation of a few cilia. Our numerical results support the mechanosensing hypothesis, and we expect that our study will stimulate further experimental investigations of mechanotransduction in the near future.

Relationship between transit time and mechanical properties of a cell through a stenosed microchannel

The changes in the mechanical properties of a cell are not only the cause of some diseases, but can also be a biomarker for some disease states. In recent times, microfluidic devices with built-in constrictions have been widely used to measure these changes. The transit time in such devices, defined as the time that a cell takes to pass through a constriction, has been found to be a crucial factor associated with the cell mechanical properties. Here, we use smoothed dissipative particle dynamics (SDPD), a particle-based numerical method, to explore the relationship between the transit time and mechanical properties of a cell. Three expressions of the transit time are developed from our simulation data, with respect to the stenosed size of constrictions, the shear modulus and bending modulus of cells, respectively. We show that a convergent constriction (the inlet is wider than the outlet), and a sharp-corner constriction (the constriction outlet is narrow) are better in identifying the differences in the transit time of cells. Moreover, the transit time increases and gradually approaches a constant as the shear modulus of cells increases, but increases first and then decreases as the bending modulus increases. These results suggest that the mechanical properties of cells can indeed be measured by analyzing their transit time, based on the recommended microfluidic device.

Off-center motion of a trapped elastic capsule in a microfluidic channel with a narrow constriction

Owing to their significance in capsule-related engineering and biomedical applications, a number of studies have considered the dynamics of elastic capsules flowing in constricted microchannels. However, these studies have focused on capsules moving along the channel centerline. In the present study, we numerically investigate the transient motion of an elastic capsule in a microfluidic channel with a rectangular constriction, which is initially trapped at the constriction inlet while off the channel centerline (i.e., on the channel bottom-wall). Under the push of the surrounding flow, the capsule can squeeze into the constriction, but only if the capsule deformability or the constriction size is sufficiently large. We find that the critical capillary number leading to the penetration of the capsule into the constriction is larger for off-centerline capsules compared to centered capsules. The centered capsule is stationary at the steady state when it remains stuck at the constriction; in contrast, the off-centerline capsule is not stationary but exhibits a tank-treading motion, i.e., its overall shape maintains a nonspherical shape with a protrusion into the constriction while its membrane exhibits a continuous rotation. Further, we examine the dependence of the capsule motion type, capsule deformation degree and membrane tension distribution on the capillary number (measuring the effects of flow strength and membrane mechanics) and constriction geometries (including the constriction height and width). Finally, we discuss the mechanism governing the capsule motion by analyzing the hydrodynamic forces acting on the capsule. The shear force acting on the capsule top owing to the fluid flow in the gap between the capsule top and the channel top-wall is the main source inducing the membrane tank-treading rotation.

Trapping and patterning of large particles and cells in a 1D ultrasonic standing wave

The use of ultrasound for trapping and patterning particles or cells in microfluidic systems is usually confined to particles which are considerably smaller than the acoustic wavelength. In this regime, the primary forces result in particle clustering at certain locations in the sound field, whilst secondary forces, those arising due to particle-particle interaction forces, assist this clustering process. Using a wavelength closer to the size of the particles allows one particle to be held at each primary force minimum. However, to achieve this, the influence of secondary forces needs to be carefully studied, as inter-particle attraction is highly undesirable. Here, we study the effect of particle size and material properties on both the primary and secondary acoustic forces as the particle diameter is increased towards the wavelength of the 1-dimensional axisymmetric ultrasonic field. We show that the resonance frequencies of the solid sphere have an important role in the resulting secondary forces which leads to a narrow band of frequencies that allow the patterning of large particles in a 1-D array. Knowledge regarding the naturally existent secondary forces would allow for system designs enabling single cell studies to be conducted in a biologically safe manner.

Optimal Branching Structure of Fluidic Networks with Permeable Walls

Biological and engineering studies of Hess-Murray's law are focused on assemblies of tubes with impermeable walls. Blood vessels and airways have permeable walls to allow the exchange of fluid and other dissolved substances with tissues. Should Hess-Murray's law hold for bifurcating systems in which the walls of the vessels are permeable to fluid? This paper investigates the fluid flow in a porous-walled T-shaped assembly of vessels. Fluid flow in this branching flow structure is studied numerically to predict the configuration that provides greater access to the flow. Our findings indicate, among other results, that an asymmetric flow (i.e., breaking the symmetry of the flow distribution) may occur in this symmetrical dichotomous system. To derive expressions for the optimum branching sizes, the hydraulic resistance of the branched system is computed. Here we show the T-shaped assembly of vessels is only conforming to Hess-Murray's law optimum as long as they have impervious walls. Findings also indicate that the optimum relationship between the sizes of parent and daughter tubes depends on the wall permeability of the assembled tubes. Our results agree with analytical results obtained from a variety of sources and provide new insights into the dynamics within the assembly of vessels.

SPH-DEM approach to numerically simulate the deformation of three-dimensional RBCs in non-uniform capillaries

Background

Blood continuously flows through the blood vessels in the human body. When blood flows through the smallest blood vessels, red blood cells (RBCs) in the blood exhibit various types of motion and deformed shapes. Computational modelling techniques can be used to successfully predict the behaviour of the RBCs in capillaries. In this study, we report the application of a meshfree particle approach to model and predict the motion and deformation of three-dimensional RBCs in capillaries.

Methods

An elastic spring network based on the discrete element method (DEM) is employed to model the three-dimensional RBC membrane. The haemoglobin in the RBC and the plasma in the blood are modelled as smoothed particle hydrodynamics (SPH) particles. For validation purposes, the behaviour of a single RBC in a simple shear flow is examined and compared against experimental results. Then simulations are carried out to predict the behaviour of RBCs in a capillary; (i) the motion of five identical RBCs in a uniform capillary, (ii) the motion of five identical RBCs with different bending stiffness (K_b) values in a stenosed capillary, (iii) the motion of three RBCs in a narrow capillary. Finally five identical RBCs are employed to determine the critical diameter of a stenosed capillary.

Results

Validation results showed a good agreement with less than 10% difference. From the above simulations, the following results are obtained; (i) RBCs exhibit different deformation behaviours due to the hydrodynamic interaction between them. (ii) Asymmetrical deformation behaviours of the RBCs are clearly observed when the bending stiffness (K_b) of the RBCs is changed. (iii) The model predicts the ability of the RBCs to squeeze through smaller blood vessels. Finally, from the simulations, the critical diameter of the stenosed section to stop the motion of blood flow is predicted.

Conclusions

A three-dimensional spring network model based on DEM in combination with the SPH method is successfully used to model the motion and deformation of RBCs in capillaries. Simulation results reveal that the condition of blood flow stopping depends on the pressure gradient of the capillary and the severity of stenosis of the capillary. In addition, this model is capable of predicting the critical diameter which prevents motion of RBCs for different blood pressures.

Dynamics of flexible fibers and vesicles in Poiseuille flow at low Reynolds number

The dynamics of flexible fibers and vesicles in unbounded planar Poiseuille flow at low Reynolds number is shown to exhibit similar basic features, when their equilibrium (moderate) aspect ratio is the same and vesicle viscosity contrast is relatively high. Tumbling, lateral migration, accumulation and shape evolution of these two types of flexible objects are analyzed numerically. The linear dependence of the accumulation position on relative bending rigidity, and other universal scalings are derived from the local shear flow approximation.

Angle of inclination of tank-treading red cells: dependence on shear rate and suspending medium

Red cells suspended in solutions much more viscous than blood plasma assume an almost steady-state orientation when sheared above a threshold value of shear rate. This orientation is a consequence of the motion of the membrane around the red cell called tank-treading. Observed along the undisturbed vorticity of the shear flow, tank-treading red cells appear as slender bodies. Their orientation can be quantified as an angle of inclination (θ) of the major axis with respect to the undisturbed flow direction. We measured θ using solution viscosities (η0) and shear rates (γ˙) covering one and three orders of magnitude, respectively. At the lower values of η0, θ was almost independent of γ˙. At the higher values of η0, θ displayed a maximum at intermediate shear rates. The respective maximal values of θ increased by ∼10° from 10.7 to 104 mPas. After accounting for the absent membrane viscosity in models by using an increased cytoplasmic viscosity, their predictions of θ agree qualitatively with our data. Comparison of the observed variation of θ at constant γ˙ with model results suggests a change in the reference configuration of the shear stiffness of the membrane.

Deformation of a micro-torque swimmer

The membrane tension of some kinds of ciliates has been suggested to regulate upward and downward swimming velocities under gravity. Despite its biological importance, deformation and membrane tension of a ciliate have not been clarified fully. In this study, we numerically investigated the deformation of a ciliate swimming freely in a fluid otherwise at rest. The cell body was modelled as a capsule with a hyperelastic membrane enclosing a Newtonian fluid. Thrust forces due to the ciliary beat were modelled as torques distributed above the cell body. The effects of membrane elasticity, the aspect ratio of the cell's reference shape, and the density difference between the cell and the surrounding fluid were investigated. The results showed that the cell deformed like a heart shape, when the capillary number was sufficiently large. Under the influence of gravity, the membrane tension at the anterior end decreased in the upward swimming while it increased in the downward swimming. Moreover, gravity-induced deformation caused the cells to move gravitationally downwards or upwards, which resulted in a positive or negative geotaxis-like behaviour with a physical origin. These results are important in understanding the physiology of a ciliate's biological responses to mechanical stimuli.

Hemodynamics in the microcirculation and in microfluidics

Hemodynamics in microcirculation is important for hemorheology and several types of circulatory disease. Although hemodynamics research has a long history, the field continues to expand due to recent advancements in numerical and experimental techniques at the micro-and nano-scales. In this paper, we review recent computational and experimental studies of blood flow in microcirculation and microfluidics. We first focus on the computational studies of red blood cell (RBC) dynamics, from the single cellular level to mesoscopic multiple cellular flows, followed by a review of recent computational adhesion models for white blood cells, platelets, and malaria-infected RBCs, in which the cell adhesion to the vascular wall is essential for cellular function. Recent developments in optical microscopy have enabled the observation of flowing blood cells in microfluidics. Experimental particle image velocimetry and particle tracking velocimetry techniques are described in this article. Advancements in micro total analysis system technologies have facilitated flowing cell separation with microfluidic devices, which can be used for biomedical applications, such as a diagnostic tool for breast cancer or large intestinal tumors. In this paper, cell-separation techniques are reviewed for microfluidic devices, emphasizing recent advances and the potential of this fast-evolving research field in the near future.

Numerical analysis of a red blood cell flowing through a thin micropore

Red blood cell (RBC) deformability plays a key role in microcirculation, especially in vessels that have diameters even smaller than the nominal cell size. In this study, we numerically investigate the dynamics of an RBC in a thin micropore. The RBC is modeled as a capsule with a thin hyperelastic membrane. In a numerical simulation, we employ a boundary element method for fluid mechanics and a finite element method for membrane mechanics. The resulting RBC deformation towards the flow direction is suppressed considerably by increased cytoplasm viscosity, whereas the gap between the cell membrane and solid wall becomes smaller with higher cytoplasm viscosity. We also measure the transit time of the RBC and find that nondimensional transit time increases nonlinearly with respect to the viscosity ratio, whereas it is invariant to the capillary number. In conclusion, cytoplasmic viscosity plays a key role in the dynamics of an RBC in a thin pore. The results of this study will be useful for designing a microfluidic device to measure cytoplasmic viscosity.

Membrane compression in tumbling and vacillating-breathing regimes for quasispherical vesicles

We derive some analytical results of a well-known model for quasispherical vesicles in a linear shear flow at low deformability. Attention is focussed on the oscillatory regimes: the tumbling (TB) mode, vacillating-breathing (VB) mode, and the transition from vacillating-breathing to tumbling, depending on a control parameter Γ. It is shown that, during the VB-to-TB transition (Γ=1), the vesicle momentarily attains its maximal extension in the vorticity direction and transits through a circular profile in the shear plane for which the radius is exactly determined. In addition, we provide an explicit analytical expression for the effective membrane tension for different types of motions. We find a critical bending number below which the membrane undergoes compression at each instant and show that, during the VB-to-TB transition, a fourth-order membrane deformation is possible.