We describe here a rigorous and accurate model for the simulation of three-dimensional deformable particles (DPs). The method is very versatile, easily simulating various types of deformable particles such as vesicles, capsules, and biological cells. Each DP is resolved explicitly and advects within the surrounding Newtonian fluid. The DPs have a preferred rest shape (e.g., spherical for vesicles, or biconcave for red blood cells). The model uses a classic hybrid system: an Eulerian approach is used for the Navier-Stokes solver (the lattice Boltzmann method) and a Lagrangian approach for the evolution of the DP mesh. Coupling is accomplished through the lattice Boltzmann velocity field, which transmits force to the membranes of the DPs. The novelty of this method resides in its ability (by design) to simulate a large number of DPs within the bounds of current computational limitations: our simple and efficient approach is to (i) use the lattice Boltzmann method because of its acknowledged efficiency at low Reynolds number and its ease of parallelization, and (ii) model the DP dynamics using a coarse mesh (approximately 500 nodes) and a spring model constraining (if necessary) local area, total area, cell volume, local curvature, and local primary stresses. We show that this approach is comparable to the more common-yet numerically expensive-approach of membrane potential function, through a series of quantitative comparisons. To demonstrate the capabilities of the model, we simulate the flow of 200 densely packed red blood cells-a computationally challenging task. The model is very efficient, requiring of the order of minutes for a single DP in a 50 μm×40 μm×40 μm simulation domain and only hours for 200 DPs in 80 μm×30 μm×30 μm. Moreover, the model is highly scalable and efficient compared to other models of blood cells in flow, making it an ideal and unique tool for studying blood flow in microvessels or vesicle or capsule flow (or a mixture of different particles). In addition to directly predicting fluid dynamics in complex suspension in any geometry, the model allows determination of accurate, empirical rules which may improve existing macroscopic, continuum models.

For more than a century, pioneering researchers have been using novel experimental and computational approaches to probe the mysteries of blood flow. Thanks to their efforts, we know that blood cells generally prefer to migrate to the axis of flow, that red and white cells segregate in flow, and that cell deformability and their tendency to reversibly aggregate contribute to the non-Newtonian nature of this unique fluid. All of these properties have beneficial physiological consequences, allowing blood to perform a variety of critical functions. Our current understanding of these unusual flow properties of blood have been made possible by the ingenuity and diligence of a number of researchers, including Harry Goldsmith, who developed novel technologies to visualize and quantify the flow of blood at the level of individual cells. Here we summarize efforts in our lab to continue this tradition and to further our understanding of how blood cells interact with each other and with the blood vessel wall.

No abstract