Numerical simulation of the motion of red blood cells and vesicles in microfluidic flows

We study the mathematical modeling and numerical simulation of the motion of red blood cells (RBC) and vesicles subject to an external incompressible flow in a microchannel. RBC and vesicles are viscoelastic bodies consisting of a deformable elastic membrane enclosing an incompressible fluid. We provide an extension of the finite element immersed boundary method by Boffi and Gastaldi (Comput Struct 81:491–501, 2003), Boffi et al. (Math Mod Meth Appl Sci 17:1479–1505, 2007), Boffi et al. (Comput Struct 85:775–783, 2007) based on a model for the membrane that additionally accounts for bending energy and also consider inflow/outflow conditions for the external fluid flow. The stability analysis requires both the approximation of the membrane by cubic splines (instead of linear splines without bending energy) and an upper bound on the inflow velocity. In the fully discrete case, the resulting CFL-type condition on the time step size is also more restrictive. We perform numerical simulations for various scenarios including the tank treading motion of vesicles in microchannels, the behavior of ‘healthy’ and ‘sick’ RBC which differ by their stiffness, and the motion of RBC through thin capillaries. The simulation results are in very good agreement with experimentally available data.

[1]  G. Cokelet,et al.  Rheology and hemodynamics. , 1980, Annual review of physiology.

[2]  Charles S. Peskin,et al.  Improved Volume Conservation in the Computation of Flows with Immersed Elastic Boundaries , 1993 .

[3]  L. Heltai,et al.  On the hyper-elastic formulation of the immersed boundary method , 2008 .

[4]  Charles S. Peskin,et al.  Fluid Flow in Collapsible Elastic Tubes: A Three-Dimensional Numerical Model , 2001 .

[5]  Hiroshi Noguchi,et al.  Swinging and tumbling of fluid vesicles in shear flow. , 2007, Physical review letters.

[6]  Michael H.G. Duits,et al.  Deformation of giant lipid bilayer vesicles in a shear flow , 1996 .

[7]  Victor Steinberg,et al.  Orientation and dynamics of a vesicle in tank-treading motion in shear flow. , 2005, Physical review letters.

[8]  T. Biben,et al.  Tumbling of vesicles under shear flow within an advected-field approach. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[9]  Gunther Guetens,et al.  Human Red Blood Cells: Rheological Aspects, Uptake, and Release of Cytotoxic Drugs , 2004, Critical reviews in clinical laboratory sciences.

[10]  Tsorng-Whay Pan,et al.  International Journal of C 2009 Institute for Scientific Numerical Analysis and Modeling Computing and Information Dynamical Simulation of Red Blood Cell Rheology in Microvessels Tsorng-whay Pan and Tong Wang , 2022 .

[11]  David Farrell,et al.  Immersed finite element method and its applications to biological systems. , 2006, Computer methods in applied mechanics and engineering.

[12]  Thomas Podgorski,et al.  Rheology of dilute suspensions of vesicles and red blood cells , 2007 .

[13]  C. Peskin,et al.  A three-dimensional computational method for blood flow in the heart. 1. Immersed elastic fibers in a viscous incompressible fluid , 1989 .

[14]  A Chabanel,et al.  Influence of cholesterol content on red cell membrane viscoelasticity and fluidity. , 1983, Biophysical journal.

[15]  Cherl‐Ho Lee,et al.  Rheological properties of erythrocytes from male hypercholesterolemia. , 2004, Microvascular research.

[16]  Hiroshi Noguchi,et al.  Dynamics of fluid vesicles in flow through structured microchannels , 2008, 0811.0862.

[17]  S. Chien,et al.  Red cell rheology in stomatocyte-echinocyte transformation: roles of cell geometry and cell shape. , 1986, Blood.

[18]  Y. C. Fung,et al.  Improved measurements of the erythrocyte geometry. , 1972, Microvascular research.

[19]  H. Noguchi,et al.  Shape transitions of fluid vesicles and red blood cells in capillary flows. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Aleksander S Popel,et al.  Red blood cell aggregation and dissociation in shear flows simulated by lattice Boltzmann method. , 2008, Journal of biomechanics.

[21]  Gerhard Gompper,et al.  Low-Reynolds-number hydrodynamics of complex fluids by multi-particle-collision dynamics , 2004 .

[22]  J S Owen,et al.  Erythrocyte echinocytosis in liver disease. Role of abnormal plasma high density lipoproteins. , 1985, The Journal of clinical investigation.

[23]  Marina V Kameneva,et al.  Development of Standard Tests to Examine Viscoelastic Properties of Blood of Experimental Animals for Pediatric Mechanical Support Device Evaluation , 2006, ASAIO journal.

[24]  C. Pozrikidis.  Axisymmetric motion of a file of red blood cells through capillaries , 2005 .

[25]  T Biben,et al.  Analytical analysis of a vesicle tumbling under a shear flow. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[26]  Gerhard Gompper,et al.  Multiparticle collision dynamics modeling of viscoelastic fluids. , 2008, The Journal of chemical physics.

[27]  G. Thurston,et al.  Clinical blood rheology. , 1979, Biorheology.

[28]  Chaouqi Misbah,et al.  Rheology and dynamics of vesicle suspension in comparison with droplet emulsion , 2008 .

[29]  N. Mohandas,et al.  Shear-Response of the Spectrin Dimer-Tetramer Equilibrium in the Red Blood Cell Membrane* , 2002, The Journal of Biological Chemistry.

[30]  H Schmid-Schönbein,et al.  The red cell as a fluid droplet: tank tread-like motion of the human erythrocyte membrane in shear flow. , 1978, Science.

[31]  Aleksander S Popel,et al.  Computational fluid dynamic simulation of aggregation of deformable cells in a shear flow. , 2005, Journal of biomechanical engineering.

[32]  D. Braasch,et al.  Red Cell Deformability and Capillary Blood Flow , 1971 .

[33]  Colette Lartigue,et al.  Tank treading and unbinding of deformable vesicles in shear flow: determination of the lift force. , 2002, Physical review letters.

[34]  N. Mohandas,et al.  Erythrocyte membrane deformability and stability: two distinct membrane properties that are independently regulated by skeletal protein associations , 1986, The Journal of cell biology.

[35]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[36]  C. Pozrikidis,et al.  Modeling and Simulation of Capsules and Biological Cells , 2003 .

[37]  Udo Seifert.  Hydrodynamic lift on bound vesicles , 1999 .

[38]  M. Abkarian,et al.  Dynamics of vesicles in a wall-bounded shear flow. , 2005, Biophysical journal.

[39]  Udo Seifert.  Fluid membranes in hydrodynamic flow fields: Formalism and an application to fluctuating quasispherical vesicles in shear flow , 1999 .

[40]  T. Biben,et al.  Steady to unsteady dynamics of a vesicle in a flow. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[41]  G B Thurston,et al.  Intrinsic viscoelasticity of blood cell suspensions: effects of erythrocyte deformability. , 1987, Biorheology.

[42]  Philippe G. Ciarlet,et al.  The finite element method for elliptic problems , 2002, Classics in applied mathematics.

[43]  Thomas Podgorski,et al.  Tumbling of lipid vesicles, enclosing a viscous fluid, under a shear flow , 2005 .

[44]  A Chabanel,et al.  Increased resistance to membrane deformation of shape-transformed human red blood cells. , 1987, Blood.

[45]  C. Peskin.  Numerical analysis of blood flow in the heart , 1977 .

[46]  L. Heltai,et al.  A finite element approach to the immersed boundary method , 2003 .

[47]  Prosenjit Bagchi,et al.  Mesoscale simulation of blood flow in small vessels. , 2007, Biophysical journal.

[48]  G. Thurston,et al.  Effects of viscoelasticity of blood on wave propagation in the circulation. , 1976, Journal of biomechanics.

[49]  Hiroshi Noguchi,et al.  Fluid vesicles with viscous membranes in shear flow. , 2004, Physical review letters.

[50]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

[51]  P. Grisvard.  Elliptic Problems in Nonsmooth Domains , 1985 .

[52]  G. Thurston,et al.  Viscoelastic properties of blood and blood analogs , 1996 .

[53]  L. Heltai,et al.  Mathematical Models and Methods in Applied Sciences Vol. 17, No. 10 (17) 1479–1505 c ○ World Scientific Publishing Company NUMERICAL STABILITY OF THE FINITE ELEMENT IMMERSED BOUNDARY METHOD , 2005 .

[54]  Seifert,et al.  Fluid Vesicles in Shear Flow. , 1996, Physical review letters.

[55]  M. Dupin,et al.  Modeling the flow of dense suspensions of deformable particles in three dimensions. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[56]  Chaouqi Misbah,et al.  Rheology of a dilute suspension of vesicles. , 2007, Physical review letters.

[57]  R. A. Cooper,et al.  Anemia with spur cells: a red cell defect acquired in serum and modified in the circulation. , 1969, The Journal of clinical investigation.

[58]  C. Pozrikidis,et al.  Finite deformation of liquid capsules enclosed by elastic membranes in simple shear flow , 1995, Journal of Fluid Mechanics.

[59]  Daniel T Chiu,et al.  A microfluidic model for single-cell capillary obstruction by Plasmodium falciparum-infected erythrocytes , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[60]  R M Hochmuth,et al.  Erythrocyte membrane elasticity and viscosity. , 1987, Annual review of physiology.

[61]  H. Chmiel,et al.  The determination of blood viscoelasticity in clinical hemorheology. , 1990 .

[62]  G. Thurston,et al.  Viscoelasticity of human blood. , 1972, Biophysical journal.

[63]  Thomas Podgorski,et al.  Micro-macro link in rheology of erythrocyte and vesicle suspensions. , 2008, Biophysical journal.

[64]  C. Pozrikidis,et al.  Effect of membrane bending stiffness on the deformation of capsules in simple shear flow , 2001, Journal of Fluid Mechanics.

[65]  Thomas Podgorski,et al.  Dynamics and rheology of a dilute suspension of vesicles: higher-order theory. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[66]  H. Schmid-schönbein,et al.  Tank Tread Motion of Red Cell Membranes in Viscometric Flow: Behavior of Intracellular and Extracellular Markers (with Film) , 1978 .

[67]  H Kiesewetter,et al.  Basic phenomena of red blood cell rouleaux formation. , 1999, Biorheology.

[68]  Shu Chien,et al.  Shear Dependence of Effective Cell Volume as a Determinant of Blood Viscosity , 1970, Science.

[69]  H Minamitani,et al.  Direct measurement of erythrocyte deformability in diabetes mellitus with a transparent microchannel capillary model and high-speed video camera system. , 2001, Microvascular research.

[70]  Shu Chien,et al.  Handbook of Bioengineering , 1986 .

[71]  Shigeo Wada,et al.  Simulation Study on Effects of Hematocrit on Blood Flow Properties Using Particle Method , 2006 .

[72]  A. Popel,et al.  Large deformation of red blood cell ghosts in a simple shear flow. , 1998, Physics of fluids.

[73]  Yaling Liu,et al.  Rheology of red blood cell aggregation by computer simulation , 2006, J. Comput. Phys..

[74]  R. Skalak,et al.  Motion of a tank-treading ellipsoidal particle in a shear flow , 1982, Journal of Fluid Mechanics.

[75]  Marina V. Kameneva,et al.  RHEOLOGICAL PROPERTIES OF BOVINE BLOOD ASSESSED BY A VISCOELASTICITY ANALYZER. POTENTIAL APPLICATIONS IN TESTING OF HEART-ASSIST DEVICE BIOCOMPATIBILITY , 2006 .

[76]  Luca Heltai,et al.  On the CFL condition for the finite element immersed boundary method , 2007 .

[77]  R Josephs,et al.  Remodeling the shape of the skeleton in the intact red cell. , 1996, Biophysical journal.

[78]  S. Chien.  Red cell deformability and its relevance to blood flow. , 1987, Annual review of physiology.

[79]  T. Secomb,et al.  Analysis of red blood cell motion through cylindrical micropores: effects of cell properties. , 1996, Biophysical journal.