More Web Proxy on the site http://driver.im/

article

A spectral boundary integral method for flowing blood cells

Authors:

Amir H. G. Isfahani,

Jonathan B. FreundAuthors Info & Claims

Journal of Computational Physics, Volume 229, Issue 10

Pages 3726 - 3744

https://doi.org/10.1016/j.jcp.2010.01.024

Published: 01 May 2010 Publication History

Abstract

A spectral boundary integral method for simulating large numbers of blood cells flowing in complex geometries is developed and demonstrated. The blood cells are modeled as finite-deformation elastic membranes containing a higher viscosity fluid than the surrounding plasma, but the solver itself is independent of the particular constitutive model employed for the cell membranes. The surface integrals developed for solving the viscous flow, and thereby the motion of the massless membrane, are evaluated using an O(NlogN) particle-mesh Ewald (PME) approach. The cell shapes, which can become highly distorted under physiologic conditions, are discretized with spherical harmonics. The resolution of these global basis functions is, of course, excellent, but more importantly they facilitate an approximate de-aliasing procedure that stabilizes the simulations without adding any numerical dissipation or further restricting the permissible numerical time step. Complex geometry no-slip boundaries are included using a constraint method that is coupled into an implicit system that is solved as part of the time advancement routine. The implementation is verified against solutions for axisymmetric flows reported in the literature, and its accuracy is demonstrated by comparison against exact solutions for relaxing surface deformations. It is also used to simulate flow of blood cells at 30% volume fraction in tubes between 4.9 and 16.9@mm in diameter. For these, it is shown to reproduce the well-known non-monotonic dependence of the effective viscosity on the tube diameter.

References

[1]

Abbitt, K.B. and Nash, G.B., Rheological properties of the blood influencing selectin-mediated adhesion of flowing leukocytes. Am. J. Physiol. Heart Circ. Physiol. v285. H229

[2]

J.C. Admas, P.N. Swarztrauber, SPHEREPAKC 2.0: a model development facility, Tech. rep., NCAR, NCAR/TN-436-STR, 1997.

[3]

Bagchi, P., Mesoscale simulation of blood flow in small vessels. Biophs. J. v92 i6. 1858-1877.

[4]

Bagchi, P., Johnson, P.C. and Popel, A.S., Computational fluid dynamic simulation of aggregation of deformable cells in a shear flow. J. Bio. Eng. v127. 1070-1080.

[5]

Bruno, O.P. and Kunyansky, L.A., A fast, high-order algorithm for the solution of surface scattering problems: Basic implementation, tests, and applications. J. Comput. Phys. v169. 80-110.

Digital Library

[6]

Canuto, C., Hussaini, M.Y., Quarteroni, A. and Zang, T.A., Spectral Methods in Fluid Dynamics. 1987. Springer-Verlag, Berlin.

[7]

Dao, M., Li, J. and Suresh, S., Molecularly based analysis of deformation of spectrin network and human erythrocyte. Mater. Sci. Eng. C. v26. 1232-1244.

[8]

Deserno, M. and Holm, C., How to mesh up Ewald sums. I. A theoretical and numerical comparison of various particle mesh routines. J. Chem. Phys. v109 i18. 7678-7693.

[9]

Deserno, M. and Holm, C., How to mesh up Ewald sums. II. An accurate error estimate for the particle-particle-particle-mesh algorithm. J. Chem. Phys. v109 i18. 7694-7701.

[10]

Dintenfass, L., Internal viscosity of the red cell and a blood viscosity equation. Nature. v219. 956-958.

[11]

Duffy, M.G., Quadrature over a pyramid or cube of integrands with a singularity at a vertex. SIAM J. Numer. Anal. v19 i6. 1260-1262.

[12]

Fåhræus, R. and Lindqvist, T., The viscosity of blood in narrow capillary tubes. Am. J. Physiol. v96. 562-568.

[13]

Firrel, J.C. and Lipowsky, H.H., Leukocyte margination and deformation in mesnenteric venules of rate. Am. J. Physiol. Heart Circ. Physiol. v256. H1667

[14]

Freund, J.B., Leukocyte margination in a model microvessel. Phys. Fluids. v19. 023301

[15]

Hashimoto, H., On the periodic fundamental solutions of the stokes equations and their application to viscous flow past a cubic array of spheres. J. Fluid Mech. v5. 317-328.

[16]

Hebecker, F.K., Efficient boundary element methods for three-dimensional exterior viscous flows. Numer. Methods Part. Diff. Equations. v2. 273-297.

[17]

Johnston, P.R. and Elliott, D., A sinh transformation for evaluating nearly singular boundary element integrals. Int. J. Numer. Methods Eng. v65. 564-678.

[18]

Kim, S. and Karrila, S.J., Microhydrodynamics: Principles and Selected Applications. 1991. Butterworth Heinemann, Stoneham, MA.

[19]

Kravchenko, A.G. and Moin, P., On the effect of numerical errors in large eddy simulations of turbulent flows. J. Comput. Phys. v131. 310-322.

Digital Library

[20]

Kress, R., Linear Integral Equations. 1989. Springer-Verlag.

[21]

Lee, L. and LeVeque, R.J., An immersed interface method for incompressible Navier-Stokes equations. SIAM J. Sci. Comput. v25 i3. 832-856.

[22]

Li, J., Dao, M., Lim, C.T. and Suresh, S., Spectrin-level modeling of the cytoskeleton and optical tweezers stretching of the erythrocyte. Biophys. J. v88. 3707-3719.

[23]

Liu, W.K., Liu, Y., Farrell, D., Zhang, L., Wang, X.S., Fukui, Y., Patankar, N., Zhang, Y. and Bajaj, C., Immersed finite element method and its applications to biological systems. Comput. Methods Appl. Mech. Eng. v195. 1722-1749.

[24]

MacMeccan, R.M., Clausen, J.R., Neitzel, G.P. and Aidun, C.K., Simulating deformable particle suspensions using a coupled lattice-Boltzmann and finite-element method. J. Fluid Mech. v618. 13-39.

[25]

Mityushev, V. and Adler, P.M., Longitudinal permeability of spatially periodic rectangular arrays of circular cylinders. I. A single cylinder in the unit cell. Z. Angew. Math. Mech. v82 i5. 335-345.

[26]

Muldowney, G.P. and Higdon, J.J.L., A spectral boundary element approach to three-dimensional Stokes flow. J. Fluid Mech. v298. 167-192.

[27]

Orszag, S.A., On the elimination of aliasing in finite difference schemes by filtering high-wavenumber components. J. Atoms. Sci. v28. 1074

[28]

Peskin, C.S., Flow patterns around heart valves: a numerical method. J. Comput. Phys. v10 i2. 252-271.

[29]

Pivkin, I.V. and Karniadakis, P.E., Accurate coarse-grained modeling of red blood cells. Phys. Rev. Lett. v101. 118105

[30]

Pogorzelski, W., Integral Equations and Their Applications. 1966. Pergamon Press.

[31]

Power, H. and Wrobel, L.C., Boundary Integral Methods in Fluid Mechanics. 1994. Computational Mechanics Publications.

[32]

Pozrikidis, C., Boundary Integral and Singularity Methods for Linearized Viscous Flow. 1992. Cambridge University Press, Cambridge.

[33]

Pozrikidis, C., Numerical simulation of the flow-induced deformation of red blood cells. Ann. Biomed. Eng. v31. 1194-1205.

[34]

Pozrikidis, C., Axisymmetric motion of a file of red blood cells through capillaries. Phys. Fluids. v17. 031503

[35]

Pries, A.R., Ley, K., Clasassen, M. and Gaehtgens, P., Red cell distribution at microvascular bifurcations. Microvasc. Res. v38. 81-101.

[36]

Pries, A.R., Neuhaus, D. and Gaehtgens, P., Blood viscosity in tube flow: dependence on diameter and hematocrit. Am. J. Physiol. Heart Circ. Physiol. v263. H1770-H1778.

[37]

Rallison, J.M. and Acrivos, A., A numerical study of the deformation and burst of a viscous drop in an extensional flow. J. Fluid Mech. v89. 191-200.

[38]

Rochal, S.B., Lorman, V.L. and Mennessier, G., Viscoelastic dynamics of spherical composite vesicles. Phys. Rev. E. v71. 021905

[39]

Roma, A.M., Peskin, C.S. and Berger, M.J., An adaptive version of the immersed boundary method. J. Comput. Phys. v153 i2. 509-534.

[40]

Saad, Y. and Schultz, M., A generalized minimal residual algorithm for solving nonsymmetric linear systems. SIAM J. Sci. Stat. Comput. v7. 856-869.

Digital Library

[41]

Saintillan, D., Darve, E. and Shaqfeh, E.S.G., A smooth particle-mesh Ewald algorithm for Stokes suspension simulations: the sedimentation of fibers. Phys. Fluids. v17. 033301

[42]

Silva, J.J.R., Power, H. and Wrobel, L.C., A hypersingular integral equation formulation for Stokes' flow in ducts. Eng. Anal. Bound. Element. v12. 185-193.

[43]

Skalak, R., Tözeren, A., Zarda, P.R. and Chien, S., Strain energy function of red blood cell membranes. Biophys. J. v13. 245

[44]

Swarztrauber, P.N. and Spotz, W.F., Generalized discrete spherical harmonic transforms. J. Comput. Phys. v159. 213-230.

Digital Library

[45]

Wang, H., Lei, T., Li, J., Huang, J. and Yao, Z., A parallel fast multipole accelerated integral equation scheme for 3D Stokes equations. Int. J. Numer. Methods Eng. v70. 812-839.

[46]

Ying, L., Biros, G. and Zorin, D., A high-order 3D boundary integral equation solver for elliptic PDEs in smooth domains. J. Comput. Phys. v219. 247-275.

Digital Library

[47]

Yoshida, K., Nishimura, N. and Kobayashi, S., Application of new fast multipole boundary integral equation method to crack problems in 3D. Eng. Anal. Bound. Elements. v25. 239-247.

[48]

Zinchenko, A.Z. and Davis, R.H., Algorithm for direct numerical simulation of emulsion flow through a granular material. J. Comput. Phys. v227. 7841-7888.

Digital Library

Cited By

Huet DWachs A(2023)A Cartesian-octree adaptive front-tracking solver for immersed biological capsules in large complex domainsJournal of Computational Physics10.1016/j.jcp.2023.112424492:COnline publication date: 18-Oct-2023
https://dl.acm.org/doi/10.1016/j.jcp.2023.112424
Sorgentone CTornberg A(2023)Estimation of quadrature errors for layer potentials evaluated near surfaces with spherical topologyAdvances in Computational Mathematics10.1007/s10444-023-10083-749:6Online publication date: 23-Nov-2023
https://dl.acm.org/doi/10.1007/s10444-023-10083-7
Rydquist GEsmaily M(2022)A cell-resolved, Lagrangian solver for modeling red blood cell dynamics in macroscale flowsJournal of Computational Physics10.1016/j.jcp.2022.111204461:COnline publication date: 15-Jul-2022
https://dl.acm.org/doi/10.1016/j.jcp.2022.111204
Show More Cited By

A spectral boundary integral method for flowing blood cells

Recommendations

Studying the blood plasma flow past a red blood cell, with the mathematical method of Kelvin's transformation
BIBE '12: Proceedings of the 2012 IEEE 12th International Conference on Bioinformatics & Bioengineering (BIBE)

A mathematical tool, namely the Kelvin transformation, has been employed in order to derive analytical expressions for important hydrodynamic quantities, aiming to the understanding and the study of the blood plasma flow past a Red Blood Cell (RBC). ...
Accelerated boundary integral method for multiphase flow in non-periodic geometries

An accelerated boundary integral method for Stokes flow of a suspension of deformable particles is presented for an arbitrary domain and implemented for the important case of a planar slit geometry. The computational complexity of the algorithm scales ...
Dynamic simulation of locally inextensible vesicles suspended in an arbitrary two-dimensional domain, a boundary integral method

We consider numerical algorithms for the simulation of hydrodynamics of two-dimensional vesicles suspended in a viscous Stokesian fluid. The motion of vesicles is governed by the interplay between hydrodynamic and elastic forces. Continuum models of ...

Comments

Please enable JavaScript to view thecomments powered by Disqus.

Information & Contributors

Information

Published In

cover image Journal of Computational Physics

Journal of Computational Physics Volume 229, Issue 10

May, 2010

448 pages

ISSN:0021-9991

Issue’s Table of Contents

Copyright © Elsevier Inc. © 2010.

Publisher

Academic Press Professional, Inc.

United States

Publication History

Published: 01 May 2010

Author Tags

Qualifiers

Article

Contributors

Other Metrics

View Article Metrics

Bibliometrics & Citations

Bibliometrics

Article Metrics

23
Total Citations
View Citations
0
Total Downloads

Downloads (Last 12 months)0
Downloads (Last 6 weeks)0

Reflects downloads up to 03 Jan 2025

Other Metrics

View Author Metrics

Citations

Cited By

Huet DWachs A(2023)A Cartesian-octree adaptive front-tracking solver for immersed biological capsules in large complex domainsJournal of Computational Physics10.1016/j.jcp.2023.112424492:COnline publication date: 18-Oct-2023
https://dl.acm.org/doi/10.1016/j.jcp.2023.112424
Sorgentone CTornberg A(2023)Estimation of quadrature errors for layer potentials evaluated near surfaces with spherical topologyAdvances in Computational Mathematics10.1007/s10444-023-10083-749:6Online publication date: 23-Nov-2023
https://dl.acm.org/doi/10.1007/s10444-023-10083-7
Rydquist GEsmaily M(2022)A cell-resolved, Lagrangian solver for modeling red blood cell dynamics in macroscale flowsJournal of Computational Physics10.1016/j.jcp.2022.111204461:COnline publication date: 15-Jul-2022
https://dl.acm.org/doi/10.1016/j.jcp.2022.111204
Lu LMorse MRahimian AStadler GZorin DTaufer MBalaji PPeña A(2019)Scalable simulation of realistic volume fraction red blood cell flows through vascular networksProceedings of the International Conference for High Performance Computing, Networking, Storage and Analysis10.1145/3295500.3356203(1-30)Online publication date: 17-Nov-2019
https://dl.acm.org/doi/10.1145/3295500.3356203
Wala MKlöckner A(2019)A fast algorithm for Quadrature by Expansion in three dimensionsJournal of Computational Physics10.1016/j.jcp.2019.03.024388:C(655-689)Online publication date: 1-Jul-2019
https://dl.acm.org/doi/10.1016/j.jcp.2019.03.024
Maxian OKassen AStrychalski W(2018)A continuous energy-based immersed boundary method for elastic shellsJournal of Computational Physics10.1016/j.jcp.2018.05.045371:C(333-362)Online publication date: 15-Oct-2018
https://dl.acm.org/doi/10.1016/j.jcp.2018.05.045
Siegel MTornberg A(2018)A local target specific quadrature by expansion method for evaluation of layer potentials in 3DJournal of Computational Physics10.1016/j.jcp.2018.03.006364:C(365-392)Online publication date: 1-Jul-2018
https://dl.acm.org/doi/10.1016/j.jcp.2018.03.006
Deng YRong JYe WGray L(2018)An efficient grid-based direct-volume integration BEM for 3D geometrically nonlinear elasticityComputational Mechanics10.1007/s00466-017-1515-z62:4(603-616)Online publication date: 1-Oct-2018
https://dl.acm.org/doi/10.1007/s00466-017-1515-z
Lu LRahimian AZorin D(2017)Contact-aware simulations of particulate Stokesian suspensionsJournal of Computational Physics10.1016/j.jcp.2017.06.039347:C(160-182)Online publication date: 15-Oct-2017
https://dl.acm.org/doi/10.1016/j.jcp.2017.06.039
Boedec GLeonetti MJaeger M(2017)Isogeometric FEM-BEM simulations of drop, capsule and vesicle dynamics in Stokes flowJournal of Computational Physics10.1016/j.jcp.2017.04.024342:C(117-138)Online publication date: 1-Aug-2017
https://dl.acm.org/doi/10.1016/j.jcp.2017.04.024
Show More Cited By

View Options

View options

Media

Figures

Other

Tables

View Issue’s Table of Contents