An Immersed Boundary method with divergence-free velocity interpolation and force spreading

Yuanxun Bao, Aleksandar Donev, Boyce E. Griffith, David M. McQueen, Charles Peskin

Research output: Contribution to journalArticle

Abstract

The Immersed Boundary (IB) method is a mathematical framework for constructing robust numerical methods to study fluid–structure interaction in problems involving an elastic structure immersed in a viscous fluid. The IB formulation uses an Eulerian representation of the fluid and a Lagrangian representation of the structure. The Lagrangian and Eulerian frames are coupled by integral transforms with delta function kernels. The discretized IB equations use approximations to these transforms with regularized delta function kernels to interpolate the fluid velocity to the structure, and to spread structural forces to the fluid. It is well-known that the conventional IB method can suffer from poor volume conservation since the interpolated Lagrangian velocity field is not generally divergence-free, and so this can cause spurious volume changes. In practice, the lack of volume conservation is especially pronounced for cases where there are large pressure differences across thin structural boundaries. The aim of this paper is to greatly reduce the volume error of the IB method by introducing velocity-interpolation and force-spreading schemes with the properties that the interpolated velocity field in which the structure moves is at least C1 and satisfies a continuous divergence-free condition, and that the force-spreading operator is the adjoint of the velocity-interpolation operator. We confirm through numerical experiments in two and three spatial dimensions that this new IB method is able to achieve substantial improvement in volume conservation compared to other existing IB methods, at the expense of a modest increase in the computational cost. Further, the new method provides smoother Lagrangian forces (tractions) than traditional IB methods. The method presented here is restricted to periodic computational domains. Its generalization to non-periodic domains is important future work.

Original languageEnglish (US)
Pages (from-to)183-206
Number of pages24
JournalJournal of Computational Physics
Volume347
DOIs
StatePublished - Oct 15 2017

Fingerprint

interpolation
Interpolation
divergence
Delta functions
Conservation
Fluids
Mathematical operators
conservation
delta function
fluids
velocity distribution
Numerical methods
operators
integral transformations
traction
viscous fluids
Costs
Experiments
costs
formulations

Keywords

  • Fluid–structure interaction
  • Force spreading
  • Immersed Boundary method
  • Incompressible flow
  • Velocity interpolation
  • Volume conservation

ASJC Scopus subject areas

  • Physics and Astronomy (miscellaneous)
  • Computer Science Applications

Cite this

An Immersed Boundary method with divergence-free velocity interpolation and force spreading. / Bao, Yuanxun; Donev, Aleksandar; Griffith, Boyce E.; McQueen, David M.; Peskin, Charles.

In: Journal of Computational Physics, Vol. 347, 15.10.2017, p. 183-206.

Research output: Contribution to journalArticle

Bao, Yuanxun ; Donev, Aleksandar ; Griffith, Boyce E. ; McQueen, David M. ; Peskin, Charles. / An Immersed Boundary method with divergence-free velocity interpolation and force spreading. In: Journal of Computational Physics. 2017 ; Vol. 347. pp. 183-206.
@article{f3c282bfb7eb42ad94e4ee0733fab401,
title = "An Immersed Boundary method with divergence-free velocity interpolation and force spreading",
abstract = "The Immersed Boundary (IB) method is a mathematical framework for constructing robust numerical methods to study fluid–structure interaction in problems involving an elastic structure immersed in a viscous fluid. The IB formulation uses an Eulerian representation of the fluid and a Lagrangian representation of the structure. The Lagrangian and Eulerian frames are coupled by integral transforms with delta function kernels. The discretized IB equations use approximations to these transforms with regularized delta function kernels to interpolate the fluid velocity to the structure, and to spread structural forces to the fluid. It is well-known that the conventional IB method can suffer from poor volume conservation since the interpolated Lagrangian velocity field is not generally divergence-free, and so this can cause spurious volume changes. In practice, the lack of volume conservation is especially pronounced for cases where there are large pressure differences across thin structural boundaries. The aim of this paper is to greatly reduce the volume error of the IB method by introducing velocity-interpolation and force-spreading schemes with the properties that the interpolated velocity field in which the structure moves is at least C1 and satisfies a continuous divergence-free condition, and that the force-spreading operator is the adjoint of the velocity-interpolation operator. We confirm through numerical experiments in two and three spatial dimensions that this new IB method is able to achieve substantial improvement in volume conservation compared to other existing IB methods, at the expense of a modest increase in the computational cost. Further, the new method provides smoother Lagrangian forces (tractions) than traditional IB methods. The method presented here is restricted to periodic computational domains. Its generalization to non-periodic domains is important future work.",
keywords = "Fluid–structure interaction, Force spreading, Immersed Boundary method, Incompressible flow, Velocity interpolation, Volume conservation",
author = "Yuanxun Bao and Aleksandar Donev and Griffith, {Boyce E.} and McQueen, {David M.} and Charles Peskin",
year = "2017",
month = "10",
day = "15",
doi = "10.1016/j.jcp.2017.06.041",
language = "English (US)",
volume = "347",
pages = "183--206",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - An Immersed Boundary method with divergence-free velocity interpolation and force spreading

AU - Bao, Yuanxun

AU - Donev, Aleksandar

AU - Griffith, Boyce E.

AU - McQueen, David M.

AU - Peskin, Charles

PY - 2017/10/15

Y1 - 2017/10/15

N2 - The Immersed Boundary (IB) method is a mathematical framework for constructing robust numerical methods to study fluid–structure interaction in problems involving an elastic structure immersed in a viscous fluid. The IB formulation uses an Eulerian representation of the fluid and a Lagrangian representation of the structure. The Lagrangian and Eulerian frames are coupled by integral transforms with delta function kernels. The discretized IB equations use approximations to these transforms with regularized delta function kernels to interpolate the fluid velocity to the structure, and to spread structural forces to the fluid. It is well-known that the conventional IB method can suffer from poor volume conservation since the interpolated Lagrangian velocity field is not generally divergence-free, and so this can cause spurious volume changes. In practice, the lack of volume conservation is especially pronounced for cases where there are large pressure differences across thin structural boundaries. The aim of this paper is to greatly reduce the volume error of the IB method by introducing velocity-interpolation and force-spreading schemes with the properties that the interpolated velocity field in which the structure moves is at least C1 and satisfies a continuous divergence-free condition, and that the force-spreading operator is the adjoint of the velocity-interpolation operator. We confirm through numerical experiments in two and three spatial dimensions that this new IB method is able to achieve substantial improvement in volume conservation compared to other existing IB methods, at the expense of a modest increase in the computational cost. Further, the new method provides smoother Lagrangian forces (tractions) than traditional IB methods. The method presented here is restricted to periodic computational domains. Its generalization to non-periodic domains is important future work.

AB - The Immersed Boundary (IB) method is a mathematical framework for constructing robust numerical methods to study fluid–structure interaction in problems involving an elastic structure immersed in a viscous fluid. The IB formulation uses an Eulerian representation of the fluid and a Lagrangian representation of the structure. The Lagrangian and Eulerian frames are coupled by integral transforms with delta function kernels. The discretized IB equations use approximations to these transforms with regularized delta function kernels to interpolate the fluid velocity to the structure, and to spread structural forces to the fluid. It is well-known that the conventional IB method can suffer from poor volume conservation since the interpolated Lagrangian velocity field is not generally divergence-free, and so this can cause spurious volume changes. In practice, the lack of volume conservation is especially pronounced for cases where there are large pressure differences across thin structural boundaries. The aim of this paper is to greatly reduce the volume error of the IB method by introducing velocity-interpolation and force-spreading schemes with the properties that the interpolated velocity field in which the structure moves is at least C1 and satisfies a continuous divergence-free condition, and that the force-spreading operator is the adjoint of the velocity-interpolation operator. We confirm through numerical experiments in two and three spatial dimensions that this new IB method is able to achieve substantial improvement in volume conservation compared to other existing IB methods, at the expense of a modest increase in the computational cost. Further, the new method provides smoother Lagrangian forces (tractions) than traditional IB methods. The method presented here is restricted to periodic computational domains. Its generalization to non-periodic domains is important future work.

KW - Fluid–structure interaction

KW - Force spreading

KW - Immersed Boundary method

KW - Incompressible flow

KW - Velocity interpolation

KW - Volume conservation

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

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

U2 - 10.1016/j.jcp.2017.06.041

DO - 10.1016/j.jcp.2017.06.041

M3 - Article

VL - 347

SP - 183

EP - 206

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -