Immersed boundary method for variable viscosity and variable density problems using fast constant-coefficient linear solvers i: Numerical method and results

Thomas G. Fai, Boyce E. Griffith, Yoichiro Mori, Charles Peskin

Research output: Contribution to journalArticle

Abstract

We present a general variable viscosity and variable density immersed boundary method that is first-order accurate in the variable density case and, for problems possessing sufficient regularity, second-order accurate in the constant density case. The viscosity and density are considered material properties and are defined by a dynamically updated tesselation. Empirical convergence rates are reported for a test problem of a two-dimensional viscoelastic shell with spatially varying material properties. The reduction to first-order accuracy in the variable density case can be avoided by using an iterative scheme, although this approach may not be efficient enough for practical use. In our time-stepping scheme, both the inertial and viscous terms are split into two parts: a constant-coefficient part that is treated implicitly, and a variable-coefficient part that is treated explicitly. This splitting allows the resulting equations to be solved efficiently using fast constantcoefficient linear solvers, and in this work, we use solvers based on the fast Fourier transform. As an application of this method, we perform fully three-dimensional, two-phase simulations of red blood cells accounting for variable viscosity and variable density. We study the behavior of red cells during shear flow and during capillary flow.

Original languageEnglish (US)
JournalSIAM Journal on Scientific Computing
Volume35
Issue number5
DOIs
StatePublished - 2013

Fingerprint

Variable Viscosity
Immersed Boundary Method
Numerical methods
Numerical Methods
Viscosity
Numerical Results
Materials properties
Coefficient
Cells
Capillary flow
Shear flow
Fast Fourier transforms
Material Properties
Blood
First-order
Red Blood Cells
Time Stepping
Fast Fourier transform
Shear Flow
Iterative Scheme

Keywords

  • Immersed boundary method
  • Red blood cells
  • Variable viscosity

ASJC Scopus subject areas

  • Applied Mathematics
  • Computational Mathematics

Cite this

Immersed boundary method for variable viscosity and variable density problems using fast constant-coefficient linear solvers i : Numerical method and results. / Fai, Thomas G.; Griffith, Boyce E.; Mori, Yoichiro; Peskin, Charles.

In: SIAM Journal on Scientific Computing, Vol. 35, No. 5, 2013.

Research output: Contribution to journalArticle

@article{e40799897e724154b3f873b51f7e1c2e,
title = "Immersed boundary method for variable viscosity and variable density problems using fast constant-coefficient linear solvers i: Numerical method and results",
abstract = "We present a general variable viscosity and variable density immersed boundary method that is first-order accurate in the variable density case and, for problems possessing sufficient regularity, second-order accurate in the constant density case. The viscosity and density are considered material properties and are defined by a dynamically updated tesselation. Empirical convergence rates are reported for a test problem of a two-dimensional viscoelastic shell with spatially varying material properties. The reduction to first-order accuracy in the variable density case can be avoided by using an iterative scheme, although this approach may not be efficient enough for practical use. In our time-stepping scheme, both the inertial and viscous terms are split into two parts: a constant-coefficient part that is treated implicitly, and a variable-coefficient part that is treated explicitly. This splitting allows the resulting equations to be solved efficiently using fast constantcoefficient linear solvers, and in this work, we use solvers based on the fast Fourier transform. As an application of this method, we perform fully three-dimensional, two-phase simulations of red blood cells accounting for variable viscosity and variable density. We study the behavior of red cells during shear flow and during capillary flow.",
keywords = "Immersed boundary method, Red blood cells, Variable viscosity",
author = "Fai, {Thomas G.} and Griffith, {Boyce E.} and Yoichiro Mori and Charles Peskin",
year = "2013",
doi = "10.1137/120903038",
language = "English (US)",
volume = "35",
journal = "SIAM Journal of Scientific Computing",
issn = "1064-8275",
publisher = "Society for Industrial and Applied Mathematics Publications",
number = "5",

}

TY - JOUR

T1 - Immersed boundary method for variable viscosity and variable density problems using fast constant-coefficient linear solvers i

T2 - Numerical method and results

AU - Fai, Thomas G.

AU - Griffith, Boyce E.

AU - Mori, Yoichiro

AU - Peskin, Charles

PY - 2013

Y1 - 2013

N2 - We present a general variable viscosity and variable density immersed boundary method that is first-order accurate in the variable density case and, for problems possessing sufficient regularity, second-order accurate in the constant density case. The viscosity and density are considered material properties and are defined by a dynamically updated tesselation. Empirical convergence rates are reported for a test problem of a two-dimensional viscoelastic shell with spatially varying material properties. The reduction to first-order accuracy in the variable density case can be avoided by using an iterative scheme, although this approach may not be efficient enough for practical use. In our time-stepping scheme, both the inertial and viscous terms are split into two parts: a constant-coefficient part that is treated implicitly, and a variable-coefficient part that is treated explicitly. This splitting allows the resulting equations to be solved efficiently using fast constantcoefficient linear solvers, and in this work, we use solvers based on the fast Fourier transform. As an application of this method, we perform fully three-dimensional, two-phase simulations of red blood cells accounting for variable viscosity and variable density. We study the behavior of red cells during shear flow and during capillary flow.

AB - We present a general variable viscosity and variable density immersed boundary method that is first-order accurate in the variable density case and, for problems possessing sufficient regularity, second-order accurate in the constant density case. The viscosity and density are considered material properties and are defined by a dynamically updated tesselation. Empirical convergence rates are reported for a test problem of a two-dimensional viscoelastic shell with spatially varying material properties. The reduction to first-order accuracy in the variable density case can be avoided by using an iterative scheme, although this approach may not be efficient enough for practical use. In our time-stepping scheme, both the inertial and viscous terms are split into two parts: a constant-coefficient part that is treated implicitly, and a variable-coefficient part that is treated explicitly. This splitting allows the resulting equations to be solved efficiently using fast constantcoefficient linear solvers, and in this work, we use solvers based on the fast Fourier transform. As an application of this method, we perform fully three-dimensional, two-phase simulations of red blood cells accounting for variable viscosity and variable density. We study the behavior of red cells during shear flow and during capillary flow.

KW - Immersed boundary method

KW - Red blood cells

KW - Variable viscosity

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

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

U2 - 10.1137/120903038

DO - 10.1137/120903038

M3 - Article

VL - 35

JO - SIAM Journal of Scientific Computing

JF - SIAM Journal of Scientific Computing

SN - 1064-8275

IS - 5

ER -