Premium
Three‐dimensional numerical simulation of red blood cell motion in Poiseuille flows
Author(s) -
Shi Lingling,
Pan TsorngWhay,
Glowinski Roland
Publication year - 2014
Publication title -
international journal for numerical methods in fluids
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.938
H-Index - 112
eISSN - 1097-0363
pISSN - 0271-2091
DOI - 10.1002/fld.3939
Subject(s) - hagen–poiseuille equation , reynolds number , mechanics , flow (mathematics) , immersed boundary method , physics , tube (container) , motion (physics) , position (finance) , capillary action , geometry , materials science , classical mechanics , boundary (topology) , mathematics , mathematical analysis , thermodynamics , turbulence , composite material , finance , economics
SUMMARY An immersed boundary method based on an FEM has been successfully combined with an elastic spring network model for simulating the dynamical behavior of a red blood cell (RBC) in Poiseuille flows. This elastic spring network preserves the biconcave shape of the RBC in the sense that after the removal of the body force for driving the Poiseuille flow, the RBC with its typical parachute shape in a tube does restore its biconcave resting shape. As a benchmark test, the relationship between the deformation index and the capillary number of the RBCs flowing through a narrow cylindrical tube has been validated. For the migration properties of a single cell in a slit Poiseuille flow, a slipper shape accompanied by a cell membrane tank‐treading motion is obtained for Re ⩾ 0 . 03 , and the cell mass center is away from the center line of the channel due to its asymmetric slipper shape. For the lower Re ⩽0.0137, an RBC with almost undeformed biconcave shape has a tumbling motion. A transition from tumbling to tank‐treading happens at the Reynolds number between 0.0137 and 0.03. In slit Poiseuille flow, the RBC can also exhibit a rolling motion like a wheel during the migration when the cell is released in the fluid flow with φ = π /2 and θ = π /2 (see Figure 12 for the definition of φ and θ ). The lower the Reynolds number, the longer the rolling motion lasts; but the equilibrium shape and position are independent from the cell initial position in the channel. Copyright © 2014 John Wiley & Sons, Ltd.