Boundary integral method for the flow of vesicles with viscosity contrast in three dimensions

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

Research output: Contribution to journalArticle

Abstract

We propose numerical algorithms for the simulation of the dynamics of three-dimensional vesicles suspended in viscous Stokesian fluid. Our method is an extension of our previous work (S.K. Veerapaneni et al., 2011) [37] to flows with viscosity contrast. This generalization requires a change in the boundary integral formulation of the solution, in which a double-layer Stokes integral is introduced, and leads to changes in the fluid dynamics due to the viscosity contrast of the vesicles, which can no longer be efficiently resolved with existing algorithms. In this paper we describe the algorithms needed to handle flows with viscosity contrast accurately and efficiently. We show that a globally semi-implicit method does not have any time-step stability constraint for flows with single and multiple vesicles with moderate viscosity contrast and the computational cost per simulation unit time is comparable to or less than that of an explicit scheme. Automatic oversampling adaptation enables us to achieve high accuracy with very low spectral resolution. We conduct numerical experiments to investigate the stability, accuracy, and the computational cost of the algorithms. Overall, our method achieves several orders of magnitude speed-up compared to the standard explicit schemes.

Original languageEnglish (US)
Pages (from-to)766-786
Number of pages21
JournalJournal of Computational Physics
Volume298
DOIs
StatePublished - 2015

Fingerprint

boundary integral method
Viscosity
viscosity
costs
Spectral resolution
viscous fluids
fluid dynamics
Fluid dynamics
spectral resolution
Costs
simulation
formulations
Fluids
Experiments

Keywords

  • Boundary integral
  • Galerkin method
  • Implicit time stepping
  • Stokes flow
  • Viscosity contrast

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy (miscellaneous)

Cite this

Boundary integral method for the flow of vesicles with viscosity contrast in three dimensions. / Rahimian, Abtin; Veerapaneni, Shravan K.; Zorin, Denis; Biros, George.

In: Journal of Computational Physics, Vol. 298, 2015, p. 766-786.

Research output: Contribution to journalArticle

Rahimian, Abtin ; Veerapaneni, Shravan K. ; Zorin, Denis ; Biros, George. / Boundary integral method for the flow of vesicles with viscosity contrast in three dimensions. In: Journal of Computational Physics. 2015 ; Vol. 298. pp. 766-786.
@article{223260c3a154467c86cecdf93e4ddc08,
title = "Boundary integral method for the flow of vesicles with viscosity contrast in three dimensions",
abstract = "We propose numerical algorithms for the simulation of the dynamics of three-dimensional vesicles suspended in viscous Stokesian fluid. Our method is an extension of our previous work (S.K. Veerapaneni et al., 2011) [37] to flows with viscosity contrast. This generalization requires a change in the boundary integral formulation of the solution, in which a double-layer Stokes integral is introduced, and leads to changes in the fluid dynamics due to the viscosity contrast of the vesicles, which can no longer be efficiently resolved with existing algorithms. In this paper we describe the algorithms needed to handle flows with viscosity contrast accurately and efficiently. We show that a globally semi-implicit method does not have any time-step stability constraint for flows with single and multiple vesicles with moderate viscosity contrast and the computational cost per simulation unit time is comparable to or less than that of an explicit scheme. Automatic oversampling adaptation enables us to achieve high accuracy with very low spectral resolution. We conduct numerical experiments to investigate the stability, accuracy, and the computational cost of the algorithms. Overall, our method achieves several orders of magnitude speed-up compared to the standard explicit schemes.",
keywords = "Boundary integral, Galerkin method, Implicit time stepping, Stokes flow, Viscosity contrast",
author = "Abtin Rahimian and Veerapaneni, {Shravan K.} and Denis Zorin and George Biros",
year = "2015",
doi = "10.1016/j.jcp.2015.06.017",
language = "English (US)",
volume = "298",
pages = "766--786",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - Boundary integral method for the flow of vesicles with viscosity contrast in three dimensions

AU - Rahimian, Abtin

AU - Veerapaneni, Shravan K.

AU - Zorin, Denis

AU - Biros, George

PY - 2015

Y1 - 2015

N2 - We propose numerical algorithms for the simulation of the dynamics of three-dimensional vesicles suspended in viscous Stokesian fluid. Our method is an extension of our previous work (S.K. Veerapaneni et al., 2011) [37] to flows with viscosity contrast. This generalization requires a change in the boundary integral formulation of the solution, in which a double-layer Stokes integral is introduced, and leads to changes in the fluid dynamics due to the viscosity contrast of the vesicles, which can no longer be efficiently resolved with existing algorithms. In this paper we describe the algorithms needed to handle flows with viscosity contrast accurately and efficiently. We show that a globally semi-implicit method does not have any time-step stability constraint for flows with single and multiple vesicles with moderate viscosity contrast and the computational cost per simulation unit time is comparable to or less than that of an explicit scheme. Automatic oversampling adaptation enables us to achieve high accuracy with very low spectral resolution. We conduct numerical experiments to investigate the stability, accuracy, and the computational cost of the algorithms. Overall, our method achieves several orders of magnitude speed-up compared to the standard explicit schemes.

AB - We propose numerical algorithms for the simulation of the dynamics of three-dimensional vesicles suspended in viscous Stokesian fluid. Our method is an extension of our previous work (S.K. Veerapaneni et al., 2011) [37] to flows with viscosity contrast. This generalization requires a change in the boundary integral formulation of the solution, in which a double-layer Stokes integral is introduced, and leads to changes in the fluid dynamics due to the viscosity contrast of the vesicles, which can no longer be efficiently resolved with existing algorithms. In this paper we describe the algorithms needed to handle flows with viscosity contrast accurately and efficiently. We show that a globally semi-implicit method does not have any time-step stability constraint for flows with single and multiple vesicles with moderate viscosity contrast and the computational cost per simulation unit time is comparable to or less than that of an explicit scheme. Automatic oversampling adaptation enables us to achieve high accuracy with very low spectral resolution. We conduct numerical experiments to investigate the stability, accuracy, and the computational cost of the algorithms. Overall, our method achieves several orders of magnitude speed-up compared to the standard explicit schemes.

KW - Boundary integral

KW - Galerkin method

KW - Implicit time stepping

KW - Stokes flow

KW - Viscosity contrast

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

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

U2 - 10.1016/j.jcp.2015.06.017

DO - 10.1016/j.jcp.2015.06.017

M3 - Article

AN - SCOPUS:84955275278

VL - 298

SP - 766

EP - 786

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -