TY - JOUR
T1 - An immersed boundary method for simulating vesicle dynamics in three dimensions
AU - Seol, Yunchang
AU - Hu, Wei Fan
AU - Kim, Yongsam
AU - Lai, Ming Chih
N1 - Publisher Copyright:
© 2016 Elsevier Inc.
PY - 2016/10/1
Y1 - 2016/10/1
N2 - We extend our previous immersed boundary (IB) method for 3D axisymmetric inextensible vesicle in Navier-Stokes flows (Hu et al., 2014 [17]) to general three dimensions. Despite a similar spirit in numerical algorithms to the axisymmetric case, the fully 3D numerical implementation is much more complicated and is far from straightforward. A vesicle membrane surface is known to be incompressible and exhibits bending resistance. As in 3D axisymmetric case, instead of keeping the vesicle locally incompressible, we adopt a modified elastic tension energy to make the vesicle surface patch nearly incompressible so that solving the unknown tension (Lagrange multiplier for the incompressible constraint) can be avoided. Nevertheless, the new elastic force derived from the modified tension energy has exactly the same mathematical form as the original one except the different definitions of tension. The vesicle surface is discretized on a triangular mesh where the elastic tension and bending force are calculated on each vertex (Lagrangian marker in the IB method) of the triangulation. A series of numerical tests on the present scheme are conducted to illustrate the robustness and applicability of the method. We perform the convergence study for the immersed boundary forces and the fluid velocity field. We then study the vesicle dynamics in various flows such as quiescent, simple shear, and gravitational flows. Our numerical results show good agreements with those obtained in previous theoretical, experimental and numerical studies.
AB - We extend our previous immersed boundary (IB) method for 3D axisymmetric inextensible vesicle in Navier-Stokes flows (Hu et al., 2014 [17]) to general three dimensions. Despite a similar spirit in numerical algorithms to the axisymmetric case, the fully 3D numerical implementation is much more complicated and is far from straightforward. A vesicle membrane surface is known to be incompressible and exhibits bending resistance. As in 3D axisymmetric case, instead of keeping the vesicle locally incompressible, we adopt a modified elastic tension energy to make the vesicle surface patch nearly incompressible so that solving the unknown tension (Lagrange multiplier for the incompressible constraint) can be avoided. Nevertheless, the new elastic force derived from the modified tension energy has exactly the same mathematical form as the original one except the different definitions of tension. The vesicle surface is discretized on a triangular mesh where the elastic tension and bending force are calculated on each vertex (Lagrangian marker in the IB method) of the triangulation. A series of numerical tests on the present scheme are conducted to illustrate the robustness and applicability of the method. We perform the convergence study for the immersed boundary forces and the fluid velocity field. We then study the vesicle dynamics in various flows such as quiescent, simple shear, and gravitational flows. Our numerical results show good agreements with those obtained in previous theoretical, experimental and numerical studies.
KW - Immersed boundary method
KW - Incompressible membrane
KW - Navier-Stokes equations
KW - Three-dimensional vesicle
UR - http://www.scopus.com/inward/record.url?scp=84977505879&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2016.06.035
DO - 10.1016/j.jcp.2016.06.035
M3 - 期刊論文
AN - SCOPUS:84977505879
SN - 0021-9991
VL - 322
SP - 125
EP - 141
JO - Journal of Computational Physics
JF - Journal of Computational Physics
ER -