An immersed boundary method for rigid bodies

Bakytzhan Kallemov, Amneet Pal Singh Bhalla, Boyce E. Griffith, Aleksandar Donev

Research output: Contribution to journalArticle

Abstract

We develop an immersed boundary (IB) method for modeling flows around fixed or moving rigid bodies that is suitable for a broad range of Reynolds numbers, including steady Stokes flow. The spatio-temporal discretization of the fluid equations is based on a standard staggered-grid approach. Fluid-body interaction is handled using Peskin's IB method; however, unlike existing IB approaches to such problems, we do not rely on penalty or fractional-step formulations. Instead, we use an unsplit scheme that ensures the no-slip constraint is enforced exactly in terms of the Lagrangian velocity field evaluated at the IB markers. Fractionalstep approaches, by contrast, can impose such constraints only approximately, which can lead to penetration of the flow into the body, and are inconsistent for steady Stokes flow. Imposing no-slip constraints exactly requires the solution of a large linear system that includes the fluid velocity and pressure as well as Lagrange multiplier forces that impose the motion of the body. The principal contribution of this paper is that it develops an efficient preconditioner for this exactly constrained IB formulation which is based on an analytical approximation to the Schur complement. This approach is enabled by the near translational and rotational invariance of Peskin's IB method. We demonstrate that only a few cycles of a geometric multigrid method for the fluid equations are required in each application of the preconditioner, and we demonstrate robust convergence of the overall Krylov solver despite the approximations made in the preconditioner. We empirically observe that to control the condition number of the coupled linear system while also keeping the rigid structure impermeable to fluid, we need to place the immersed boundary markers at a distance of about two grid spacings, which is significantly larger from what has been recommended in the literature for elastic bodies. We demonstrate the advantage of our monolithic solver over split solvers by computing the steady state flow through a two-dimensional nozzle at several Reynolds numbers. We apply the method to a number of benchmark problems at zero and finite Reynolds numbers, and we demonstrate first-order convergence of the method to several analytical solutions and benchmark computations.

Original languageEnglish (US)
Pages (from-to)79-141
Number of pages63
JournalCommunications in Applied Mathematics and Computational Science
Volume11
Issue number1
DOIs
StatePublished - 2016

Fingerprint

Immersed Boundary Method
Immersed Boundary
Rigid Body
Fluid
Reynolds number
Fluids
Preconditioner
Linear systems
Stokes Flow
Steady Flow
Slip
Demonstrate
Rigid structures
Linear Systems
Lagrange multipliers
Body fluids
Benchmark
Fractional Step
Invariance
Staggered Grid

Keywords

  • Fluid-structure interaction
  • Immersed boundary method
  • Rigid body

ASJC Scopus subject areas

  • Applied Mathematics
  • Computational Theory and Mathematics
  • Computer Science Applications

Cite this

An immersed boundary method for rigid bodies. / Kallemov, Bakytzhan; Bhalla, Amneet Pal Singh; Griffith, Boyce E.; Donev, Aleksandar.

In: Communications in Applied Mathematics and Computational Science, Vol. 11, No. 1, 2016, p. 79-141.

Research output: Contribution to journalArticle

Kallemov, Bakytzhan ; Bhalla, Amneet Pal Singh ; Griffith, Boyce E. ; Donev, Aleksandar. / An immersed boundary method for rigid bodies. In: Communications in Applied Mathematics and Computational Science. 2016 ; Vol. 11, No. 1. pp. 79-141.
@article{5bd014bee8d74220b0d8dcf83c148cef,
title = "An immersed boundary method for rigid bodies",
abstract = "We develop an immersed boundary (IB) method for modeling flows around fixed or moving rigid bodies that is suitable for a broad range of Reynolds numbers, including steady Stokes flow. The spatio-temporal discretization of the fluid equations is based on a standard staggered-grid approach. Fluid-body interaction is handled using Peskin's IB method; however, unlike existing IB approaches to such problems, we do not rely on penalty or fractional-step formulations. Instead, we use an unsplit scheme that ensures the no-slip constraint is enforced exactly in terms of the Lagrangian velocity field evaluated at the IB markers. Fractionalstep approaches, by contrast, can impose such constraints only approximately, which can lead to penetration of the flow into the body, and are inconsistent for steady Stokes flow. Imposing no-slip constraints exactly requires the solution of a large linear system that includes the fluid velocity and pressure as well as Lagrange multiplier forces that impose the motion of the body. The principal contribution of this paper is that it develops an efficient preconditioner for this exactly constrained IB formulation which is based on an analytical approximation to the Schur complement. This approach is enabled by the near translational and rotational invariance of Peskin's IB method. We demonstrate that only a few cycles of a geometric multigrid method for the fluid equations are required in each application of the preconditioner, and we demonstrate robust convergence of the overall Krylov solver despite the approximations made in the preconditioner. We empirically observe that to control the condition number of the coupled linear system while also keeping the rigid structure impermeable to fluid, we need to place the immersed boundary markers at a distance of about two grid spacings, which is significantly larger from what has been recommended in the literature for elastic bodies. We demonstrate the advantage of our monolithic solver over split solvers by computing the steady state flow through a two-dimensional nozzle at several Reynolds numbers. We apply the method to a number of benchmark problems at zero and finite Reynolds numbers, and we demonstrate first-order convergence of the method to several analytical solutions and benchmark computations.",
keywords = "Fluid-structure interaction, Immersed boundary method, Rigid body",
author = "Bakytzhan Kallemov and Bhalla, {Amneet Pal Singh} and Griffith, {Boyce E.} and Aleksandar Donev",
year = "2016",
doi = "10.2140/camcos.2016.11.79",
language = "English (US)",
volume = "11",
pages = "79--141",
journal = "Communications in Applied Mathematics and Computational Science",
issn = "1559-3940",
publisher = "Mathematical Sciences Publishers",
number = "1",

}

TY - JOUR

T1 - An immersed boundary method for rigid bodies

AU - Kallemov, Bakytzhan

AU - Bhalla, Amneet Pal Singh

AU - Griffith, Boyce E.

AU - Donev, Aleksandar

PY - 2016

Y1 - 2016

N2 - We develop an immersed boundary (IB) method for modeling flows around fixed or moving rigid bodies that is suitable for a broad range of Reynolds numbers, including steady Stokes flow. The spatio-temporal discretization of the fluid equations is based on a standard staggered-grid approach. Fluid-body interaction is handled using Peskin's IB method; however, unlike existing IB approaches to such problems, we do not rely on penalty or fractional-step formulations. Instead, we use an unsplit scheme that ensures the no-slip constraint is enforced exactly in terms of the Lagrangian velocity field evaluated at the IB markers. Fractionalstep approaches, by contrast, can impose such constraints only approximately, which can lead to penetration of the flow into the body, and are inconsistent for steady Stokes flow. Imposing no-slip constraints exactly requires the solution of a large linear system that includes the fluid velocity and pressure as well as Lagrange multiplier forces that impose the motion of the body. The principal contribution of this paper is that it develops an efficient preconditioner for this exactly constrained IB formulation which is based on an analytical approximation to the Schur complement. This approach is enabled by the near translational and rotational invariance of Peskin's IB method. We demonstrate that only a few cycles of a geometric multigrid method for the fluid equations are required in each application of the preconditioner, and we demonstrate robust convergence of the overall Krylov solver despite the approximations made in the preconditioner. We empirically observe that to control the condition number of the coupled linear system while also keeping the rigid structure impermeable to fluid, we need to place the immersed boundary markers at a distance of about two grid spacings, which is significantly larger from what has been recommended in the literature for elastic bodies. We demonstrate the advantage of our monolithic solver over split solvers by computing the steady state flow through a two-dimensional nozzle at several Reynolds numbers. We apply the method to a number of benchmark problems at zero and finite Reynolds numbers, and we demonstrate first-order convergence of the method to several analytical solutions and benchmark computations.

AB - We develop an immersed boundary (IB) method for modeling flows around fixed or moving rigid bodies that is suitable for a broad range of Reynolds numbers, including steady Stokes flow. The spatio-temporal discretization of the fluid equations is based on a standard staggered-grid approach. Fluid-body interaction is handled using Peskin's IB method; however, unlike existing IB approaches to such problems, we do not rely on penalty or fractional-step formulations. Instead, we use an unsplit scheme that ensures the no-slip constraint is enforced exactly in terms of the Lagrangian velocity field evaluated at the IB markers. Fractionalstep approaches, by contrast, can impose such constraints only approximately, which can lead to penetration of the flow into the body, and are inconsistent for steady Stokes flow. Imposing no-slip constraints exactly requires the solution of a large linear system that includes the fluid velocity and pressure as well as Lagrange multiplier forces that impose the motion of the body. The principal contribution of this paper is that it develops an efficient preconditioner for this exactly constrained IB formulation which is based on an analytical approximation to the Schur complement. This approach is enabled by the near translational and rotational invariance of Peskin's IB method. We demonstrate that only a few cycles of a geometric multigrid method for the fluid equations are required in each application of the preconditioner, and we demonstrate robust convergence of the overall Krylov solver despite the approximations made in the preconditioner. We empirically observe that to control the condition number of the coupled linear system while also keeping the rigid structure impermeable to fluid, we need to place the immersed boundary markers at a distance of about two grid spacings, which is significantly larger from what has been recommended in the literature for elastic bodies. We demonstrate the advantage of our monolithic solver over split solvers by computing the steady state flow through a two-dimensional nozzle at several Reynolds numbers. We apply the method to a number of benchmark problems at zero and finite Reynolds numbers, and we demonstrate first-order convergence of the method to several analytical solutions and benchmark computations.

KW - Fluid-structure interaction

KW - Immersed boundary method

KW - Rigid body

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

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

U2 - 10.2140/camcos.2016.11.79

DO - 10.2140/camcos.2016.11.79

M3 - Article

VL - 11

SP - 79

EP - 141

JO - Communications in Applied Mathematics and Computational Science

JF - Communications in Applied Mathematics and Computational Science

SN - 1559-3940

IS - 1

ER -