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

2017 ◽  
Vol 95 (6) ◽  
Author(s):  
Ting Ye ◽  
Nhan Phan-Thien ◽  
Chwee Teck Lim ◽  
Lina Peng ◽  
Huixin Shi
Author(s):  
Orest Shardt ◽  
Jos Derksen

We describe the direct simulation of high-solids-fraction suspensions of non-spherical rigid particles that are slightly denser than the fluid. The lattice-Boltzmann method is used to solve the flow of the interstitial Newtonian fluid, and the immersed boundary method is used to enforce a no-slip boundary condition at the surface of each particle. The surface points for the immersed boundary method are also employed for collision handling by applying repulsive forces between nearby surface points. Due to the finite number of these points, the method simulates rough surface collisions. We also discuss methods for integrating the equations of particle motion at low density ratios and propose a method with improved accuracy. Rigid particles shaped like red blood cells were simulated. Simulations of a single particle showed that the particle settles in its original orientation when the Reynolds number is low (1.2) but flips to a higher drag orientation when the Reynolds number is higher (7.3). A simulation with a 45% solids volume fraction and a low solid over fluid density ratio showed the possibility of simulating blood as it is found in the body. A simulation at a lower solids volume fraction (35%) was used to compare the results with the erythrocyte sedimentation rate (ESR), a common blood test. The sedimentation rate was estimated as 0.2 mm/hr, which is an order of magnitude lower than a typical ESR of about 6 mm/hr for a healthy adult. The most likely reasons for the discrepancy are the omission of agglomeration-inducing inter-cellular forces from the simulations and the treatment of the red blood cells as rigid particles.


Author(s):  
Danny Bluestein ◽  
João S. Soares ◽  
Peng Zhang ◽  
Chao Gao ◽  
Seetha Pothapragada ◽  
...  

The coagulation cascade of blood may be initiated by flow induced platelet activation, which prompts clot formation in prosthetic cardiovascular devices and arterial disease processes. While platelet activation may be induced by biochemical agonists, shear stresses arising from pathological flow patterns enhance the propensity of platelets to activate and initiate the intrinsic pathway of coagulation, leading to thrombosis. Upon activation platelets undergo complex biochemical and morphological changes: organelles are centralized, membrane glycoproteins undergo conformational changes, and adhesive pseudopods are extended. Activated platelets polymerize fibrinogen into a fibrin network that enmeshes red blood cells. Activated platelets also cross-talk and aggregate to form thrombi. Current numerical simulations to model this complex process mostly treat blood as a continuum and solve the Navier-Stokes equations governing blood flow, coupled with diffusion-convection-reaction equations. It requires various complex constitutive relations or simplifying assumptions, and is limited to μm level scales. However, molecular mechanisms governing platelet shape change upon activation and their effect on rheological properties can be in the nm level scales. To address this challenge, a multiscale approach which departs from continuum approaches, may offer an effective means to bridge the gap between macroscopic flow and cellular scales. Molecular dynamics (MD) and dissipative particle dynamics (DPD) methods have been employed in recent years to simulate complex processes at the molecular scales, and various viscous fluids at low-to-high Reynolds numbers at mesoscopic scales. Such particle methods possess important properties at the mesoscopic scale: complex fluids with heterogeneous particles can be modeled, allowing the simulation of processes which are otherwise very difficult to solve by continuum approaches. It is becoming a powerful tool for simulating complex blood flow, red blood cells interactions, and platelet-mediated thrombosis involving platelet activation, aggregation, and adhesion.


2016 ◽  
Vol 13 (04) ◽  
pp. 1641003
Author(s):  
L. W. Zhou ◽  
Yu-Qian Zhang ◽  
Xiao-Long Deng ◽  
M. B. Liu

This paper presents a dissipative particle dynamics (DPDs) method for investigating the movement and deformation of biconcave shape red blood cells (RBCs) with the worm-like chain (WLC) bead spring. First, the stretching of a RBC is modeled and the obtained shape evolution of the cell agrees well with experimental results. Second, the movement and deformation of a RBC in shear flows are investigated and three typical modes (tumbling, intermittent and tank-treading) are observed. Lastly, an illustrating example of multi-RBCs in Poiseuille flow in a tube is simulated. We conclude that the presented DPD method with WLC spring can effectively model the movement and deformation of bioconcave cells.


2017 ◽  
Vol 217 ◽  
pp. 171-179 ◽  
Author(s):  
Ansel L. Blumers ◽  
Yu-Hang Tang ◽  
Zhen Li ◽  
Xuejin Li ◽  
George E. Karniadakis

2009 ◽  
Vol 131 (7) ◽  
Author(s):  
Xiaobo Gong ◽  
Kazuyasu Sugiyama ◽  
Shu Takagi ◽  
Yoichiro Matsumoto

The deformation of multiple red blood cells in a capillary flow was studied numerically. The immersed boundary method was used for the fluid red blood cells interaction. The membrane of the red blood cell was modeled as a hyperelastic thin shell. The numerical results show that the apparent viscosity in the capillary flow is more sensitive to the change of shear coefficient of the membrane than the bending coefficient and surface dilation coefficient, and the increase in the shear coefficient results in an increase in the pressure drop in the blood flow in capillary vessels in order to sustain the same flux rate of red blood cells.


Sign in / Sign up

Export Citation Format

Share Document