A numerical method for simulating the dynamics of 3D axisymmetric vesicles suspended in viscous flows

We extend [Shravan K. Veerapaneni, Denis Gueyffier, Denis Zorin, George Biros, A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D, Journal of Computational Physics 228(7) (2009) 2334-2353] to the case of three-dimensional axisymmetric vesicles of spherical or toroidal topology immersed in viscous flows. Although the main components of the algorithm are similar in spirit to the 2D case-spectral approximation in space, semi-implicit time-stepping scheme-the main differences are that the bending and viscous force require new analysis, the linearization for the semi-implicit schemes must be rederived, a fully implicit scheme must be used for the toroidal topology to eliminate a CFL-type restriction and a novel numerical scheme for the evaluation of the 3D Stokes single layer potential on an axisymmetric surface is necessary to speed up the calculations. By introducing these novel components, we obtain a time-scheme that experimentally is unconditionally stable, has low cost per time step, and is third-order accurate in time. We present numerical results to analyze the cost and convergence rates of the scheme. To verify the solver, we compare it to a constrained variational approach to compute equilibrium shapes that does not involve interactions with a viscous fluid. To illustrate the applicability of method, we consider a few vesicle-flow interaction problems: the sedimentation of a vesicle, interactions of one and three vesicles with a background Poiseuille flow.

[1]  C. Pozrikidis,et al.  Interfacial dynamics for Stokes flow , 2001 .

[2]  U. Seifert,et al.  Influence of shear flow on vesicles near a wall: A numerical study. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[3]  Thomas Podgorski,et al.  Deformation of vesicles flowing through capillaries , 2004 .

[4]  Seifert,et al.  Fluid Vesicles in Shear Flow. , 1996, Physical review letters.

[5]  Lin Ma,et al.  Viscous regularization and r-adaptive remeshing for finite element analysis of lipid membrane mechanics , 2007, J. Comput. Phys..

[6]  C. Pozrikidis.  Axisymmetric motion of a file of red blood cells through capillaries , 2005 .

[7]  M I Bloor,et al.  Method for efficient shape parametrization of fluid membranes and vesicles. , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[8]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[9]  C. Pozrikidis,et al.  Effect of membrane bending stiffness on the deformation of capsules in simple shear flow , 2001, Journal of Fluid Mechanics.

[10]  M. Kropinski.  An efficient numerical method for studying interfacial motion in two-dimensional creeping flows , 2001 .

[11]  Leslie Greengard,et al.  A fast multipole method for the three-dimensional Stokes equations , 2008, J. Comput. Phys..

[12]  F. Campelo,et al.  Dynamic model and stationary shapes of fluid vesicles , 2006, The European physical journal. E, Soft matter.

[13]  Luiz C. Wrobel,et al.  Boundary Integral Methods in Fluid Mechanics , 1995 .

[14]  Michael Shelley,et al.  Simulating the dynamics and interactions of flexible fibers in Stokes flows , 2004 .

[15]  L. Trefethen.  Spectral Methods in MATLAB , 2000 .

[16]  Andrei Ludu.  Nonlinear Waves and Solitons on Contours and Closed Surfaces , 2007 .

[17]  Seifert,et al.  Shape transformations of vesicles: Phase diagram for spontaneous- curvature and bilayer-coupling models. , 1991, Physical review. A, Atomic, molecular, and optical physics.

[18]  V. Rovenski,et al.  Differential Geometry of Curves and Surfaces: A Concise Guide , 2005 .

[19]  Q. Du,et al.  A phase field approach in the numerical study of the elastic bending energy for vesicle membranes , 2004 .

[20]  Jian Zhang,et al.  Adaptive Finite Element Method for a Phase Field Bending Elasticity Model of Vesicle Membrane Deformations , 2008, SIAM J. Sci. Comput..

[21]  George Biros,et al.  A boundary integral method for simulating the dynamics of inextensible vesicles suspended in a viscous fluid in 2D , 2009, J. Comput. Phys..

[22]  Alex Solomonoff,et al.  Accuracy Enhancement for Higher Derivatives using Chebyshev Collocation and a Mapping Technique , 1997, SIAM J. Sci. Comput..

[23]  Udo Seifert,et al.  Configurations of fluid membranes and vesicles , 1997 .

[24]  Steven J. Ruuth,et al.  Implicit-explicit methods for time-dependent partial differential equations , 1995 .

[25]  Uri M. Ascher,et al.  Computer methods for ordinary differential equations and differential-algebraic equations , 1998 .

[26]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[27]  Feng Feng,et al.  Finite element modeling of lipid bilayer membranes , 2006, J. Comput. Phys..

[28]  Bradley K. Alpert,et al.  Hybrid Gauss-Trapezoidal Quadrature Rules , 1999, SIAM J. Sci. Comput..