A fast algorithm for simulating vesicle flows in three dimensions

Shravan K. Veerapaneni, Abtin Rahimian, George Biros, Denis Zorin

Research output: Contribution to journalArticle

Abstract

Vesicles are locally-inextensible fluid membranes that can sustain bending. In this paper, we extend the study of Veerapaneni et al. [S.K. Veerapaneni, D. Gueyffier, G. Biros, D. Zorin, A numerical method for simulating the dynamics of 3D axisymmetric vesicles suspended in viscous flows, Journal of Computational Physics 228 (19) (2009) 7233-7249] to general non-axisymmetric vesicle flows in three dimensions.Although the main components of the algorithm are similar in spirit to the axisymmetric case (spectral approximation in space, semi-implicit time-stepping scheme), important new elements need to be introduced for a full 3D method. In particular, spatial quantities are discretized using spherical harmonics, and quadrature rules for singular surface integrals need to be adapted to this case; an algorithm for surface reparameterization is needed to ensure stability of the time-stepping scheme, and spectral filtering is introduced to maintain reasonable accuracy while minimizing computational costs. To characterize the stability of the scheme and to construct preconditioners for the iterative linear system solvers used in the semi-implicit time-stepping scheme, we perform a spectral analysis of the evolution operator on the unit sphere.By introducing these algorithmic components, we obtain a time-stepping scheme that circumvents the stability constraint on the time-step and achieves spectral accuracy in space. We present results to analyze the cost and convergence rates of the overall scheme. To illustrate the applicability of the new method, we consider a few vesicle-flow interaction problems: a single vesicle in relaxation, sedimentation, shear flows, and many-vesicle flows.

Original languageEnglish (US)
Pages (from-to)5610-5634
Number of pages25
JournalJournal of Computational Physics
Volume230
Issue number14
DOIs
StatePublished - Jun 20 2011

Fingerprint

Flow interactions
Viscous flow
Shear flow
Sedimentation
Spectrum analysis
Linear systems
Costs
Numerical methods
costs
Physics
Membranes
viscous flow
spherical harmonics
linear systems
Fluids
shear flow
quadratures
spectrum analysis
membranes
operators

Keywords

  • Boundary integral methods
  • Fast algorithms
  • High-order methods
  • Vesicle simulations

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy (miscellaneous)

Cite this

A fast algorithm for simulating vesicle flows in three dimensions. / Veerapaneni, Shravan K.; Rahimian, Abtin; Biros, George; Zorin, Denis.

In: Journal of Computational Physics, Vol. 230, No. 14, 20.06.2011, p. 5610-5634.

Research output: Contribution to journalArticle

Veerapaneni, Shravan K. ; Rahimian, Abtin ; Biros, George ; Zorin, Denis. / A fast algorithm for simulating vesicle flows in three dimensions. In: Journal of Computational Physics. 2011 ; Vol. 230, No. 14. pp. 5610-5634.
@article{baa39ddab31b46eda2de44ea23c90ffb,
title = "A fast algorithm for simulating vesicle flows in three dimensions",
abstract = "Vesicles are locally-inextensible fluid membranes that can sustain bending. In this paper, we extend the study of Veerapaneni et al. [S.K. Veerapaneni, D. Gueyffier, G. Biros, D. Zorin, A numerical method for simulating the dynamics of 3D axisymmetric vesicles suspended in viscous flows, Journal of Computational Physics 228 (19) (2009) 7233-7249] to general non-axisymmetric vesicle flows in three dimensions.Although the main components of the algorithm are similar in spirit to the axisymmetric case (spectral approximation in space, semi-implicit time-stepping scheme), important new elements need to be introduced for a full 3D method. In particular, spatial quantities are discretized using spherical harmonics, and quadrature rules for singular surface integrals need to be adapted to this case; an algorithm for surface reparameterization is needed to ensure stability of the time-stepping scheme, and spectral filtering is introduced to maintain reasonable accuracy while minimizing computational costs. To characterize the stability of the scheme and to construct preconditioners for the iterative linear system solvers used in the semi-implicit time-stepping scheme, we perform a spectral analysis of the evolution operator on the unit sphere.By introducing these algorithmic components, we obtain a time-stepping scheme that circumvents the stability constraint on the time-step and achieves spectral accuracy in space. We present results to analyze the cost and convergence rates of the overall scheme. To illustrate the applicability of the new method, we consider a few vesicle-flow interaction problems: a single vesicle in relaxation, sedimentation, shear flows, and many-vesicle flows.",
keywords = "Boundary integral methods, Fast algorithms, High-order methods, Vesicle simulations",
author = "Veerapaneni, {Shravan K.} and Abtin Rahimian and George Biros and Denis Zorin",
year = "2011",
month = "6",
day = "20",
doi = "10.1016/j.jcp.2011.03.045",
language = "English (US)",
volume = "230",
pages = "5610--5634",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "14",

}

TY - JOUR

T1 - A fast algorithm for simulating vesicle flows in three dimensions

AU - Veerapaneni, Shravan K.

AU - Rahimian, Abtin

AU - Biros, George

AU - Zorin, Denis

PY - 2011/6/20

Y1 - 2011/6/20

N2 - Vesicles are locally-inextensible fluid membranes that can sustain bending. In this paper, we extend the study of Veerapaneni et al. [S.K. Veerapaneni, D. Gueyffier, G. Biros, D. Zorin, A numerical method for simulating the dynamics of 3D axisymmetric vesicles suspended in viscous flows, Journal of Computational Physics 228 (19) (2009) 7233-7249] to general non-axisymmetric vesicle flows in three dimensions.Although the main components of the algorithm are similar in spirit to the axisymmetric case (spectral approximation in space, semi-implicit time-stepping scheme), important new elements need to be introduced for a full 3D method. In particular, spatial quantities are discretized using spherical harmonics, and quadrature rules for singular surface integrals need to be adapted to this case; an algorithm for surface reparameterization is needed to ensure stability of the time-stepping scheme, and spectral filtering is introduced to maintain reasonable accuracy while minimizing computational costs. To characterize the stability of the scheme and to construct preconditioners for the iterative linear system solvers used in the semi-implicit time-stepping scheme, we perform a spectral analysis of the evolution operator on the unit sphere.By introducing these algorithmic components, we obtain a time-stepping scheme that circumvents the stability constraint on the time-step and achieves spectral accuracy in space. We present results to analyze the cost and convergence rates of the overall scheme. To illustrate the applicability of the new method, we consider a few vesicle-flow interaction problems: a single vesicle in relaxation, sedimentation, shear flows, and many-vesicle flows.

AB - Vesicles are locally-inextensible fluid membranes that can sustain bending. In this paper, we extend the study of Veerapaneni et al. [S.K. Veerapaneni, D. Gueyffier, G. Biros, D. Zorin, A numerical method for simulating the dynamics of 3D axisymmetric vesicles suspended in viscous flows, Journal of Computational Physics 228 (19) (2009) 7233-7249] to general non-axisymmetric vesicle flows in three dimensions.Although the main components of the algorithm are similar in spirit to the axisymmetric case (spectral approximation in space, semi-implicit time-stepping scheme), important new elements need to be introduced for a full 3D method. In particular, spatial quantities are discretized using spherical harmonics, and quadrature rules for singular surface integrals need to be adapted to this case; an algorithm for surface reparameterization is needed to ensure stability of the time-stepping scheme, and spectral filtering is introduced to maintain reasonable accuracy while minimizing computational costs. To characterize the stability of the scheme and to construct preconditioners for the iterative linear system solvers used in the semi-implicit time-stepping scheme, we perform a spectral analysis of the evolution operator on the unit sphere.By introducing these algorithmic components, we obtain a time-stepping scheme that circumvents the stability constraint on the time-step and achieves spectral accuracy in space. We present results to analyze the cost and convergence rates of the overall scheme. To illustrate the applicability of the new method, we consider a few vesicle-flow interaction problems: a single vesicle in relaxation, sedimentation, shear flows, and many-vesicle flows.

KW - Boundary integral methods

KW - Fast algorithms

KW - High-order methods

KW - Vesicle simulations

UR - http://www.scopus.com/inward/record.url?scp=79955741985&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=79955741985&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2011.03.045

DO - 10.1016/j.jcp.2011.03.045

M3 - Article

AN - SCOPUS:79955741985

VL - 230

SP - 5610

EP - 5634

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 14

ER -