High-volume fraction simulations of two-dimensional vesicle suspensions

Bryan Quaife,George Biros

DOI: https://doi.org/10.1016/j.jcp.2014.06.013

2014-04-15

Abstract:We consider numerical algorithms for the simulation of the rheology of two-dimensional vesicles suspended in a viscous Stokesian fluid. The vesicle evolution dynamics is governed by hydrodynamic and elastic forces. The elastic forces are due to local inextensibility of the vesicle membrane and resistance to bending. Numerically resolving vesicle flows poses several challenges. For example, we need to resolve moving interfaces, address stiffness due to bending, enforce the inextensibility constraint, and efficiently compute the (non-negligible) long-range hydrodynamic interactions. Our method is based on the work of {\em Rahimian, Veerapaneni, and Biros, "Dynamic simulation of locally inextensible vesicles suspended in an arbitrary two-dimensional domain, a boundary integral method", Journal of Computational Physics, 229 (18), 2010}. It is a boundary integral formulation of the Stokes equations coupled to the interface mass continuity and force balance. We extend the algorithms presented in that paper to increase the robustness of the method and enable simulations with concentrated suspensions. In particular, we propose a scheme in which both intra-vesicle and inter-vesicle interactions are treated semi-implicitly. In addition we use special integration for near-singular integrals and we introduce a spectrally accurate collision detection scheme. We test the proposed methodologies on both unconfined and confined flows for vesicles whose internal fluid may have a viscosity contrast with the bulk medium. Our experiments demonstrate the importance of treating both intra-vesicle and inter-vesicle interactions accurately.

Numerical Analysis

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the numerical algorithm problem in the rheological simulation of two - dimensional vesicle suspensions at high volume fractions. Specifically, the evolution dynamics of vesicles in a viscous Stokes fluid is affected by hydrodynamic and elastic forces. The elastic force mainly originates from the local inextensibility and bending resistance of the vesicle membrane. Numerically resolving vesicle flow faces multiple challenges, such as the need to resolve moving interfaces, handle the stiffness caused by bending, enforce inextensibility constraints, and efficiently calculate long - range hydrodynamic interactions.

The main contributions of the paper are:
1. **Implicit Inter - Vesicle Interactions**: A new semi - implicit time - stepping scheme is proposed, which can handle both inter - vesicle interactions and self - interactions simultaneously, thereby improving the stability of the method.
2. **Near - Singular Integrals**: When vesicles are close to each other, the hydrodynamic interaction terms need to calculate near - singular integrals. The paper introduces an efficient fifth - order - accurate near - singular integral method.
3. **Spectral Collision Detection**: A spectral collision detection scheme based on potential theory is proposed, which is suitable for smooth geometries and can achieve linear - time complexity in the number of computational degrees of freedom.
4. **Pressure and Stress Field Calculation**: A method for calculating the average pressure and stress in the region of interest is provided.

Through these improvements, the paper aims to improve the efficiency and accuracy of the simulation of high - concentration vesicle suspensions.