“…While in the particle-based methods, such as dissipative particle dynamics (DPD), blood cell models are constructed using DPD particles and thus are naturally assimilated with the background flow, continuum-based RBC models often implemented a boundary integral algorithm or the immersed boundary method (IBM) to couple RBC models with the background flow, which are solved using different solvers, including the finite volume method, the finite element method and the lattice Boltzmann method [ 31 , 32 , 33 , 34 , 35 , 36 , 37 , 38 , 39 , 40 , 41 ]. Motivated by experimental studies, recent progress in computational modeling has enabled simulations of fluid dynamics and cell aggregation dynamics under physiological and pathological states [ 10 , 27 , 42 , 43 , 44 , 45 , 46 , 47 , 48 , 49 , 50 ], which has provided insight into the pathogenesis of the disease as well as facilitated the development of therapeutic treatments [ 51 , 52 ]. For example, RBC models developed using DPD are widely applied to simulate the deformation and aggregation of healthy RBC doublets [ 42 ] together with their effects on the blood cell dynamics in stenosed microvessels [ 43 ].…”