Exact Singularity Subtraction from Boundary Integral Equations in Modeling Vesicles and Red Blood Cells

The study of vesicles, capsules and red blood cells (RBCs) under flow is a field of active research, belonging to the general problematic of fluid/structure interactions. Here, we are interested in modeling vesicles, capsules and RBCs using a boundary integral formulation, and focus on exact singularity subtractions of the kernel of the integral equations in 3D. In order to increase the precision of singular and near-singular integration, we propose here a refinement procedure in the vicinity of the pole of the Green-Oseen kernel. The refinement is performed homogeneously everywhere on the source surface in order to reuse the additional quadrature nodes when calculating boundary integrals in multiple target points. We also introduce a multi-level look-up algorithm in order to select the additional quadrature nodes in vicinity of the pole of the Green-Oseen kernel. The expected convergence rate of the proposed algorithm is of order$\mathcal{O}(1/N^2)$ while the computational complexity is of order$\mathcal{O}$($N^2$ln$N$), where $N$ is the number of degrees of freedom used for surface discretization. Several numerical tests are presented to demonstrate the convergence and the efficiency of the method.     

Lattice Boltzmann Study on the Motion of Dual Droplets in Microchannels With Contact Angle Hysteresis北大核心CSCD

接触角滞后表现为流体在非理想固体表面上运动时前进接触角和后退接触角不同,是两相流体在润湿表面上流动的重要现象.该文采用改进的伪势格子Boltzmann(LB)多组分模型,并与几何润湿边界条件相结合,研究了两个液滴在具有接触角滞后性微通道表面上的运动行为,主要研究了通道内特征数、通道表面性质以及液滴初始参数的影响.研究结果表明：毛细数的增大有助于液滴的移动,然而并不利于液滴的排出,且毛细数的增加对上游液滴的影响大于其对下游液滴的影响;另一方面,接触角滞后性窗口越大,液滴运动和形变更迟缓,但形变程度更明显,两液滴更早地发生合并,但更晚地排出管道;液滴间距的增加使液滴的运动行为在不同阶段表现为不同的模式,但都导致通道中残留小液滴,使得液滴排出通道的时间增加.研究结果还表明：上游液滴和下游液滴的相对尺寸差距越大,越不利于液滴排出管道.     

Numerical Simulation of Red Blood Cell Suspensions Behind a Moving Interface in a Capillary

Computational modeling and simulation are presented on the motion of red blood cells behind a moving interface in a capillary. The methodology is based on an immersed boundary method and the skeleton structure of the red blood cell (RBC) membrane is modeled as a spring network. As by the nature of the problem, the computational domain is moving with either a designated RBC or an interface in an infinitely long two-dimensional channel with an undisturbed flow field in front of the computational domain. The tanking-treading and the inclination angle of a cell in a simple shear flow are briefly discussed for the validation purpose. We then present and discuss the results of the motion of red blood cells behind a moving interface in a capillary, which show that the RBCs with higher velocity than the interface speed form a concentrated slug behind the moving interface.     

A DLM/FD/IB Method for Simulating Compound Cell Interacting with Red Blood Cells in a Microchannel

In this article, a computational model and related methodologies have been tested for simulating the motion of a malaria infected red blood cell (iRBC for short) in Poiseuille flow at low Reynolds numbers. Besides the deformability of the red blood cell membrane, the migration of a neutrally buoyant particle (used to model the malaria parasite inside the membrane) is another factor to determine the iRBC motion. Typically an iRBC oscillates in a Poiseuille flow due to the competition between these two factors. The interaction of an iRBC and several RBCs in a narrow channel shows that, at lower flow speed, the iRBC can be easily pushed toward the wall and stay there to block the channel. But, at higher flow speed, RBCs and iRBC stay in the central region of the channel since their migrations are dominated by the motion of the RBC membrane.