A high-resolution fast boundary-integral method is developed for simulating red blood cell motion in complex geometries. The algorithm employs a particle-mesh-Ewald method (PME) to achieve computational efficiency superior to that of standard boundary-element implementations. The computational expense scales with O(N logN), where N is the number of collocation points distributed over the cell membranes. The no-slip boundary condition on vessel walls with complex geometry is enforced implicitly in terms of a linear system for unknown forces required to stop the flow at the walls. In the numerical implementation, cell shapes are represented by spherical harmonic functions interpolated through collocation points. Because the resolution of the global spectral basis functions is perfect in some sense, accurate solutions can be obtained efficiently and error-free interpolation can be achieved for switching resolutions. The procedure facilitates the control of aliasing error and circumvents explicit filtering and implicit numerical dissipation; both would degrade the numerical accuracy. The resolution of the scheme can be set based on desired accuracy, unconstrained by stability considerations. The overall approach, motivation, advantages, and drawbacks of the numerical scheme are discussed in detail. The solver is demonstrated for case studies involving the motion of a red blood cell through a constriction, the computation of the effective blood viscosity in tube flow, the transport of a leukocyte in a microvessel, and flow in a model network. Considering the simplicity of the finite-deformation elastic constitutive model used to describe the cell membrane, the calculated effective viscosity reproduces remarkably well the well-known non-monotonic dependence on vessel diameter.
ASJC Scopus subject areas
- Biochemistry, Genetics and Molecular Biology(all)
- Physics and Astronomy(all)