Application of Chimera grid to modelling cell motion and aggregation in a narrow tube

Biomechanics in thrombus formation from direct cellular simulations

Numerically reproducing the process of thrombus formation is highly desired for understanding its mechanism but still remains challenging due to the polydisperse feature of blood components and their multiple biochemical or biomechanical behaviors involved. We numerically implemented a simplified version of the process from the perspective of biomechanics, using a mesoscale particle-based method, smoothed dissipative particle dynamics-immersed boundary method. This version covers the adhesion and aggregation of platelets (PLTs), the deformation and aggregation of red blood cells (RBCs), and the interaction between PLTs and RBCs, as well as the blockage of microvessels. Four critical factors that can affect thrombus formation were investigated: the velocity of blood flow, the adhesive ability of PLTs, the interaction strength between PLTs and RBCs, and the deformability of RBCs. Increasing the velocity of blood flow was found to be the most effective way to reduce the microvessel blockage, and reducing the adhesive ability of PLTs is also a direct and efficient way. However, decreasing the interaction strength between PLTs and RBCs sometimes does not alleviate thrombus formation, and similarly, increasing the deformability of RBCs does not have a significant improvement for the severely blocked microvessel. These results imply that maintaining high-rate blood flow plays a crucial role in the prevention and treatment of thrombosis, which is even more effective than antiplatelet or anticoagulant drugs. The drugs or treatments concentrating on reducing the PLT-RBC interaction or softening the RBCs may not have a significant effect on the thrombosis.

Margination and adhesion dynamics of tumor cells in a real microvascular network

In tumor metastasis, the margination and adhesion of tumor cells are two critical and closely related steps, which may determine the destination where the tumor cells extravasate to. We performed a direct three-dimensional simulation on the behaviors of the tumor cells in a real microvascular network, by a hybrid method of the smoothed dissipative particle dynamics and immersed boundary method (SDPD-IBM). The tumor cells are found to adhere at the microvascular bifurcations more frequently, and there is a positive correlation between the adhesion of the tumor cells and the wall-directed force from the surrounding red blood cells (RBCs). The larger the wall-directed force is, the closer the tumor cells are marginated towards the wall, and the higher the probability of adhesion behavior happen is. A relatively low or high hematocrit can help to prevent the adhesion of tumor cells, and similarly, increasing the shear rate of blood flow can serve the same purpose. These results suggest that the tumor cells may be more likely to extravasate at the microvascular bifurcations if the blood flow is slow and the hematocrit is moderate.

The key events of thrombus formation: platelet adhesion and aggregation

Thrombus formation is a complex, dynamic and multistep process, involving biochemical reactions, mechanical stimulation, hemodynamics, and so on. In this study, we concentrate on its two crucial steps: (i) platelets adhered to a vessel wall, or simply platelet adhesion, and (ii) platelets clumping and arrested to the adherent platelets, named platelet aggregation. We report the first direct simulation of three modes of platelet adhesion, detachment, rolling adhesion and firm adhesion, as well as the formation, disintegration, arrestment and consolidation of platelet plugs. The results show that the bond dissociation in the detachment mode is mainly attributed to a high probability of rupturing bonds, such that any existing bond can be quickly ruptured and all bonds would be completely broken. In the rolling adhesion, however, it is mainly attributed to the strong traction from the shear flow or erythrocytes, causing that the bonds are ruptured at the trailing edge of the platelet. The erythrocytes play an important role in platelet activities, such as the formation, disintegration, arrestment and consolidation of platelet plugs. They exert an aggregate force on platelets, a repulsion at a near distance but an attraction at a far distance to the platelets. This aggregate force can promote platelets to form a plug and/or bring along a part of a platelet plug causing its disintegration. It also greatly influences the arrestment and consolidation of platelet plugs, together with the adhesive force from the thrombus.

Red blood cell distribution in a microvascular network with successive bifurcations

Nonproportional RBC distribution is an important characteristic in microvascular networks, which can result in heterogeneity of oxygen supply that may cause ischemic death in severe cases. In this paper, we perform three-dimensional numerical simulations of a large number of RBCs in a microvascular network, by using a hybrid method of smoothed dissipative particle dynamic and immersed boundary method. The distribution of multiple RBCs in a T-bifurcation is first simulated as a validation study, and a reasonable agreement is observed both qualitatively and quantitatively on the RBC flux between our results and the previously published numerical and empirical results. Next, the distribution of a large number of RBCs in a microvascular network is investigated, including the effects of cell deformability, aggregation and tube hematocrit. The simulation results indicate that decreased deformability and increased aggregation strength have a similar effect on the RBC distribution: the large RBC flux becomes larger, but the small becomes smaller. A high hematocrit also causes a similar phenomenon that the RBCs are more apt to flow into a high RBC-flux branch, because they are arranged compactly into a rouleaux and difficultly broken up at a high hematocrit. These results imply that lower cell deformability, stronger aggregation or higher tube hematocrit would be conducive to the phase separation of hematocrit and plasma skimming processes in microcirculation.

Numerical design of a microfluidic chip for probing mechanical properties of cells

Microfluidic chips have been widely used to probe the mechanical properties of cells, which are recognized as a promising label-free biomarker for some diseases. In our previous work (Ye et al., 2018), we have studied the relationships between the transit time and the mechanical properties of a cell flowing through a microchannel with a single constriction, which potentially forms a basis for a microfluidic chip to measure cell's mechanical properties. Here, we investigate this microfluidic chip design and examine its potential in performances. We first develop the simultaneous dependence of the transit time on both the shear and bending moduli of a cell, and then examine the chip sensitivity with respect to the cell mechanical properties while serializing a single constriction along the flow direction. After that, we study the effect of the flow velocity on the transit time, and also test the chip's ability to identify heterogeneous cells with different mechanical properties. The results show that the microfluidic chip designed is capable of identifying heterogeneous cells, even when only one unhealthy cell is included. The serialization of chip can greatly increase the chip sensitivity with respect to the mechanical properties of cells. The flow with a higher velocity helps in not only promoting the chip throughput, but also in providing more accurate transit time measurements, because the cell prefers a symmetric deformation under a high velocity.

A review of numerical methods for red blood cell flow simulation

In this review, we provide an overview of the simulation techniques employed for modelling the flow of red blood cells (RBCs) in blood plasma. The scope of this review omits the fluid modelling aspect while focusing on other key components in the RBC-plasma model such as (1) describing the RBC deformation with shell-based and spring-based RBC models, (2) constitutive models for RBC aggregation based on bridging theory and depletion theory and (3) additional strategies required for completing the RBC-plasma flow model. These include topics such as modelling fluid-structure interaction with the immersed boundary method and boundary integral method, and updating the variations in multiphase fluid property through the employment of index field methods. Lastly, we summarily discuss the current state and aims of RBC modelling and suggest some research directions for the further development of this field of modelling.

AN EFFICIENT RED BLOOD CELL MODEL IN THE FRAME OF IB-LBM AND ITS APPLICATION

A two-dimensional red blood cell (RBC) membrane model based on elastic and Euler–Bernoulli beam theories is introduced in the frame of immersed boundary-lattice Boltzmann method (IB-LBM). The effect of the flexible membrane is handled by the immersed boundary method in which the stress exerted by the RBC on the ambient fluid is spread onto the collocated grid points near the boundary. The fluid dynamics is obtained by solving the discrete lattice Boltzmann equation. A "ghost shape", to which the RBC returns when restoring, is introduced by prescribing a bending force along the boundary. Numerical examples involving tumbling, tank-treading and RBC aggregation in shear flow and deformation and restoration in poiseuille flow are presented to verify the method and illustrate its efficiency. As an application of the present method, a ten-RBC colony being compressed through a stenotic microvessel is studied focusing the cell–cell interaction strength. Quantitative comparisons of the pressure and velocity on specified microvessel interfaces are made between each aggregation case. It reveals that the stronger aggregation may lead to more resistance against blood flow and result in higher pressure difference at the stenosis.

Computational fluid dynamics of aggregating red blood cells in postcapillary venules

Aggregate formation of red blood cells (RBCs) in a postcapillary venular bifurcation is investigated with three-dimensional computer simulations using the Chimera grid method. Interaction energy between the RBCs is modelled by a depletion interaction theory; RBCs are modelled as rigid oblate ellipsoids. The cell-cell interactions of RBCs are strongly dependent on vessel geometry and shear rates. The experimental data on vessel geometry, pseudoshear rates, and Dextran concentration obtained in our previous in vivo RBC aggregation study in postcapillary venules of the rat spinotrapezius muscle were used to simulate RBC aggregation. The computational results were compared to the experimental results from the in vivo study. The results show that cells have a larger tendency to form an aggregate under reduced flows. Aggregate formation also depends on the angle and location of the cells before they enter the bifurcation region. Comparisons with experimental data are discussed.

Numerical simulation of rheology of red blood cell rouleaux in microchannels

An elastic spring model is applied to simulate the skeletal structure of the red blood cell (RBC) membrane and to study the dynamical behaviors of the red blood cell rouleaux (aggregates) in microchannels. The biconcave shape of RBCs in static plasma and the tank-treading phenomenon of single RBCs in simple shear flows have been successfully captured using this model. The aggregation and dissociation of RBCs with different deformability have been investigated in both shear and Poiseuille flows by taking into consideration the rheology of the cells and the intercellular interaction kinetics. It is found that the equilibrium configuration of the rouleaux formed under no-flow condition, the motion of the rouleaux in the flows, and the rheological behavior of individual cells in the rouleaux is closely related to the intercellular interaction strength, hydrodynamic viscous forces, and the deformability of the cell membrane.

Effects of erythrocyte deformability and aggregation on the cell free layer and apparent viscosity of microscopic blood flows

Concentrated erythrocyte (i.e., red blood cell) suspensions flowing in microchannels have been simulated with an immersed-boundary lattice Boltzmann algorithm, to examine the cell layer development process and the effects of cell deformability and aggregation on hemodynamic and hemorheological behaviors. The cells are modeled as two-dimensional deformable biconcave capsules and experimentally measured cell properties have been utilized. The aggregation among cells is modeled by a Morse potential. The flow development process demonstrates how red blood cells migrate away from the boundary toward the channel center, while the suspending plasma fluid is displaced to the cell free layer regions left by the migrating cells. Several important characteristics of microscopic blood flows observed experimentally have been well reproduced in our model, including the cell free layer, blunt velocity profile, changes in apparent viscosity, and the Fahraeus effect. We found that the cell free layer thickness increases with both cell deformability and aggregation strength. Due to the opposing effects of the cell free layer lubrication and the high viscosity of cell-concentrated core, the influence of aggregation is complex but the lubrication effect appears to dominate, causing the relative apparent viscosity to decrease with aggregation. It appears therefore that the immersed-boundary lattice Boltzmann numerical model may be useful in providing valuable information on microscopic blood flows in various microcirculation situations.

Red blood cell aggregation and dissociation in shear flows simulated by lattice Boltzmann method

In this paper we develop a lattice Boltzmann algorithm to simulate red blood cell (RBC) behavior in shear flows. The immersed boundary method is employed to incorporate the fluid-membrane interaction between the flow field and deformable cells. The cell membrane is treated as a neo-Hookean viscoelastic material and a Morse potential is adopted to model the intercellular interaction. Utilizing the available mechanical properties of RBCs, multiple cells have been studied in shear flows using a two-dimensional approximation. These cells aggregate and form a rouleau under the action of intercellular interaction. The equilibrium configuration is related to the interaction strength. The end cells exhibit concave shapes under weak interaction and convex shapes under strong interaction. In shear flows, such a rouleau-like aggregate will rotate or be separated, depending on the relative strengths of the intercellular interaction and hydrodynamic viscous forces. These behaviors are qualitatively similar to experimental observations and show the potential of this numerical scheme for future studies of blood flow in microvessels.