Validation of a simple method for representing spheres and slender bodies in an immersed boundary method for Stokes flow on an unbounded domain

Thomas T. Bringley, Charles Peskin

Research output: Contribution to journalArticle

Abstract

We test the efficacy of using a single Lagrangian point to represent a sphere, and a one-dimensional array of such points to represent a slender body, in a new immersed boundary method for Stokes flow. A numerical parameter, the spacing of the Eulerian grid, is used to determine the effective radius of the immersed sphere or slender body. Such representations are much less expensive computationally than those with two or three-dimensional meshes of Lagrangian points. To perform this test, we develop a numerical method to solve the discretized Stokes equations on an unbounded Eulerian grid which contains an arbitrary configuration of Lagrangian points that apply force to the fluid and that move with the fluid. We compare results computed with this new immersed boundary method to known results for spheres and rigid cylinders in Stokes flow in R3. We find that, for certain choices of parameters, the interactions with the fluid of a single Lagrangian point accurately replicate those of a sphere of some particular radius, independent of the location of the point with respect to the Eulerian grid. The interactions of a linear array of Lagrangian points, for certain choices of parameters, accurately replicate those of a cylinder of some particular radius, independent of the position and orientation of the array with respect to the Eulerian grid. The effective radius of the sphere and the effective radius of the cylinder turn out to be related in a simple and natural way. Our results suggest recipes for choosing parameters that should be useful to practitioners. One surprising result is that one must not use too many Lagrangian points in an array. Another is that the approximate delta functions traditionally used in the immersed boundary method perform much better than higher order delta functions with the same support.

Original languageEnglish (US)
Pages (from-to)5397-5425
Number of pages29
JournalJournal of Computational Physics
Volume227
Issue number11
DOIs
StatePublished - May 10 2008

Fingerprint

slender bodies
Stokes flow
Delta functions
radii
grids
delta function
Fluids
fluids
linear arrays
Numerical methods
mesh
spacing
interactions
configurations

Keywords

  • Immersed boundary method
  • Numerical methods
  • Slender-body theory
  • Stokes equations

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy(all)

Cite this

Validation of a simple method for representing spheres and slender bodies in an immersed boundary method for Stokes flow on an unbounded domain. / Bringley, Thomas T.; Peskin, Charles.

In: Journal of Computational Physics, Vol. 227, No. 11, 10.05.2008, p. 5397-5425.

Research output: Contribution to journalArticle

@article{dda783a958564284a154b168146bed39,
title = "Validation of a simple method for representing spheres and slender bodies in an immersed boundary method for Stokes flow on an unbounded domain",
abstract = "We test the efficacy of using a single Lagrangian point to represent a sphere, and a one-dimensional array of such points to represent a slender body, in a new immersed boundary method for Stokes flow. A numerical parameter, the spacing of the Eulerian grid, is used to determine the effective radius of the immersed sphere or slender body. Such representations are much less expensive computationally than those with two or three-dimensional meshes of Lagrangian points. To perform this test, we develop a numerical method to solve the discretized Stokes equations on an unbounded Eulerian grid which contains an arbitrary configuration of Lagrangian points that apply force to the fluid and that move with the fluid. We compare results computed with this new immersed boundary method to known results for spheres and rigid cylinders in Stokes flow in R3. We find that, for certain choices of parameters, the interactions with the fluid of a single Lagrangian point accurately replicate those of a sphere of some particular radius, independent of the location of the point with respect to the Eulerian grid. The interactions of a linear array of Lagrangian points, for certain choices of parameters, accurately replicate those of a cylinder of some particular radius, independent of the position and orientation of the array with respect to the Eulerian grid. The effective radius of the sphere and the effective radius of the cylinder turn out to be related in a simple and natural way. Our results suggest recipes for choosing parameters that should be useful to practitioners. One surprising result is that one must not use too many Lagrangian points in an array. Another is that the approximate delta functions traditionally used in the immersed boundary method perform much better than higher order delta functions with the same support.",
keywords = "Immersed boundary method, Numerical methods, Slender-body theory, Stokes equations",
author = "Bringley, {Thomas T.} and Charles Peskin",
year = "2008",
month = "5",
day = "10",
doi = "10.1016/j.jcp.2008.01.048",
language = "English (US)",
volume = "227",
pages = "5397--5425",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "11",

}

TY - JOUR

T1 - Validation of a simple method for representing spheres and slender bodies in an immersed boundary method for Stokes flow on an unbounded domain

AU - Bringley, Thomas T.

AU - Peskin, Charles

PY - 2008/5/10

Y1 - 2008/5/10

N2 - We test the efficacy of using a single Lagrangian point to represent a sphere, and a one-dimensional array of such points to represent a slender body, in a new immersed boundary method for Stokes flow. A numerical parameter, the spacing of the Eulerian grid, is used to determine the effective radius of the immersed sphere or slender body. Such representations are much less expensive computationally than those with two or three-dimensional meshes of Lagrangian points. To perform this test, we develop a numerical method to solve the discretized Stokes equations on an unbounded Eulerian grid which contains an arbitrary configuration of Lagrangian points that apply force to the fluid and that move with the fluid. We compare results computed with this new immersed boundary method to known results for spheres and rigid cylinders in Stokes flow in R3. We find that, for certain choices of parameters, the interactions with the fluid of a single Lagrangian point accurately replicate those of a sphere of some particular radius, independent of the location of the point with respect to the Eulerian grid. The interactions of a linear array of Lagrangian points, for certain choices of parameters, accurately replicate those of a cylinder of some particular radius, independent of the position and orientation of the array with respect to the Eulerian grid. The effective radius of the sphere and the effective radius of the cylinder turn out to be related in a simple and natural way. Our results suggest recipes for choosing parameters that should be useful to practitioners. One surprising result is that one must not use too many Lagrangian points in an array. Another is that the approximate delta functions traditionally used in the immersed boundary method perform much better than higher order delta functions with the same support.

AB - We test the efficacy of using a single Lagrangian point to represent a sphere, and a one-dimensional array of such points to represent a slender body, in a new immersed boundary method for Stokes flow. A numerical parameter, the spacing of the Eulerian grid, is used to determine the effective radius of the immersed sphere or slender body. Such representations are much less expensive computationally than those with two or three-dimensional meshes of Lagrangian points. To perform this test, we develop a numerical method to solve the discretized Stokes equations on an unbounded Eulerian grid which contains an arbitrary configuration of Lagrangian points that apply force to the fluid and that move with the fluid. We compare results computed with this new immersed boundary method to known results for spheres and rigid cylinders in Stokes flow in R3. We find that, for certain choices of parameters, the interactions with the fluid of a single Lagrangian point accurately replicate those of a sphere of some particular radius, independent of the location of the point with respect to the Eulerian grid. The interactions of a linear array of Lagrangian points, for certain choices of parameters, accurately replicate those of a cylinder of some particular radius, independent of the position and orientation of the array with respect to the Eulerian grid. The effective radius of the sphere and the effective radius of the cylinder turn out to be related in a simple and natural way. Our results suggest recipes for choosing parameters that should be useful to practitioners. One surprising result is that one must not use too many Lagrangian points in an array. Another is that the approximate delta functions traditionally used in the immersed boundary method perform much better than higher order delta functions with the same support.

KW - Immersed boundary method

KW - Numerical methods

KW - Slender-body theory

KW - Stokes equations

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

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

U2 - 10.1016/j.jcp.2008.01.048

DO - 10.1016/j.jcp.2008.01.048

M3 - Article

AN - SCOPUS:42649127705

VL - 227

SP - 5397

EP - 5425

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 11

ER -