Shape changes and motion of a vesicle in a fluid using a lattice Boltzmann model

Huabing Li,Houhui Yi,Xiaowen Shan,Haiping Fang

DOI: https://doi.org/10.48550/arXiv.physics/0607074

2007-01-12

Abstract:We study the deformation and motion of an erythrocyte in fluid flows via a lattice Boltzmann method. To this purpose, the bending rigidity and the elastic modulus of isotropic dilation are introduced and incorporated with the lattice Boltzmann simulation, and the membrane-flow interactions on both sides of the membrane are carefully examined. We find that the static biconcave shape of an erythrocyte is quite stable and can effectively resist the pathological changes on their membrane. Further, our simulation results show that in shear flow, the erythrocyte will be highly flattened and undergo tank tread-like motion. This phenomenon has been observed by experiment very long time ago, but it has feazed the boundary integral and singularity methods up to the present. Because of its intrinsically parallel dynamics, this lattice Boltzmann method is expected to find wide applications for both single and multi-vesicles suspension as well as complex open membranes in various fluid flows for a wide range of Reynolds numbers.

Fluid Dynamics

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to study the deformation and movement of erythrocytes in fluid flow, especially the stability of their biconcave shape and the tank - tread - like motion that occurs in shear flow. Specifically:

1. **Stability of erythrocyte shape**: The paper explores whether the static biconcave shape of erythrocytes can effectively resist pathological changes on the membrane. By introducing the bending rigidity of the membrane and the elastic modulus of isotropic dilation, and combining with the lattice Boltzmann method (LBM) for simulation, the research finds that this biconcave shape is very stable.

2. **Deformation and movement in shear flow**: In shear flow, erythrocytes will become highly flattened and exhibit tank - tread - like motion. Although this phenomenon has been experimentally observed for a long time, traditional boundary integral methods and singularity methods are difficult to accurately simulate. The paper points out that due to the inherent parallel characteristics of LBM, this method can more effectively simulate this process.

3. **Broad application prospects**: The paper also emphasizes the application potential of LBM in simulating single or multiple vesicle suspensions and complex open membranes in various fluid flows, especially in different Reynolds number ranges.

### Key formulas

- **Bending energy**:
\[
\varphi_B=\frac{1}{2}k_B\int k^2 \, dl
\]
where \(k_B\) is the bending modulus, \(k\) is the curvature of the membrane, and \(l\) is the arc length.

- **Viscoelastic stress**:
\[
T_{mn} = 2\eta\dot{\varepsilon}_{mn}
\]
where \(T_{mn}\) is the viscous stress, \(\dot{\varepsilon}_{mn}\) is the strain rate, and \(\eta\) is the viscosity coefficient.

- **Isoperimetric elastic potential energy**:
\[
\varphi_k=\frac{1}{2}k_k\sum_{i = 1}^{N}(l_i - l_0)^2
\]
where \(k_k\) is the elastic coefficient of the membrane, \(l_0\) and \(l_i\) are the original length and the simulated length respectively.

- **Trans - membrane pressure difference**:
\[
\Delta\mu=RT\ln\left(\frac{p_{\text{out}}}{p_{\text{in}}}\right)
\]
where \(p_{\text{out}}\) and \(p_{\text{in}}\) are the pressures outside and inside the erythrocyte respectively, \(R\) is the gas constant, and \(T\) is the temperature.

### Summary

By introducing LBM, this paper successfully solves the problems of deformation and movement of erythrocytes in fluid flow, especially the stability of their biconcave shape and their behavior in shear flow. This method not only verifies the experimental results but also shows its broad application prospects in simulating complex fluid systems.