Simulation of a flapping flexible filament in a flowing soap film by the immersed boundary method

Luoding Zhu, Charles Peskin

Research output: Contribution to journalArticle

Abstract

This paper reports the computer simulation of a flapping flexible filament in a flowing soap film using the immersed boundary method. Our mathematical formulation includes filament mass and elasticity, gravity, air resistance, and the two wires that bound the flowing soap film. The incompressible viscous Navier-Stokes equations, which are used to describe the motion of the soap film and filament in our formulation, are discretized on a fixed uniform Eulerian lattice while the filament equations are discretized on a moving Lagrangian array of points which do not necessarily coincide with the fixed Eulerian mesh points of the fluid computation. The interaction between the filament and the soap film is handled by a smoothed approximation to the Dirac delta function. This delta function approximation is used not only to interpolate the fluid velocity and to apply force to the fluid (as is commonly done in immersed boundary computations), but also to handle the mass of the filament, which is represented in our calculation as delta function layer of fluid mass density supported along the immersed filament. Because of this nonuniform density, we need to use a multigrid method for solving the discretized fluid equations. This replaces the FFT-based method that is commonly used in the uniform-density case. Our main results are as follows. (i) The sustained flapping of the filament only occurs when filament mass is included in the formulation of the model; within a certain range of mass, the more the mass of the filament, the bigger the amplitude of the flapping. (ii) When the length of filament is short enough (below some critical length), the filament always approaches its straight (rest) state, in which the filament points downstream; but when the length is larger, the system is bistable, which means that it can settle into either state (rest state or sustained flapping) depending on the initial conditions. This numerical result we observed in computer simulation is the same as that of the laboratory experiment even though the Reynolds number of the computations is lower than that of the laboratory experiment by two orders of magnitude.

Original languageEnglish (US)
Pages (from-to)452-468
Number of pages17
JournalJournal of Computational Physics
Volume179
Issue number2
DOIs
StatePublished - Jul 1 2002

Fingerprint

flapping
soaps
Soaps (detergents)
filaments
Delta functions
Fluids
simulation
delta function
Computer simulation
fluids
Fast Fourier transforms
Navier Stokes equations
Elasticity
Gravitation
Reynolds number
Experiments
formulations
Wire
computerized simulation
Air

Keywords

  • Bistability
  • Computational fluid dynamics
  • Flag-in-wind
  • Flapping filament
  • Immersed boundary method
  • Multigrid method

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy(all)

Cite this

Simulation of a flapping flexible filament in a flowing soap film by the immersed boundary method. / Zhu, Luoding; Peskin, Charles.

In: Journal of Computational Physics, Vol. 179, No. 2, 01.07.2002, p. 452-468.

Research output: Contribution to journalArticle

@article{d130db4cea2d4a63b5700ba1b7d54153,
title = "Simulation of a flapping flexible filament in a flowing soap film by the immersed boundary method",
abstract = "This paper reports the computer simulation of a flapping flexible filament in a flowing soap film using the immersed boundary method. Our mathematical formulation includes filament mass and elasticity, gravity, air resistance, and the two wires that bound the flowing soap film. The incompressible viscous Navier-Stokes equations, which are used to describe the motion of the soap film and filament in our formulation, are discretized on a fixed uniform Eulerian lattice while the filament equations are discretized on a moving Lagrangian array of points which do not necessarily coincide with the fixed Eulerian mesh points of the fluid computation. The interaction between the filament and the soap film is handled by a smoothed approximation to the Dirac delta function. This delta function approximation is used not only to interpolate the fluid velocity and to apply force to the fluid (as is commonly done in immersed boundary computations), but also to handle the mass of the filament, which is represented in our calculation as delta function layer of fluid mass density supported along the immersed filament. Because of this nonuniform density, we need to use a multigrid method for solving the discretized fluid equations. This replaces the FFT-based method that is commonly used in the uniform-density case. Our main results are as follows. (i) The sustained flapping of the filament only occurs when filament mass is included in the formulation of the model; within a certain range of mass, the more the mass of the filament, the bigger the amplitude of the flapping. (ii) When the length of filament is short enough (below some critical length), the filament always approaches its straight (rest) state, in which the filament points downstream; but when the length is larger, the system is bistable, which means that it can settle into either state (rest state or sustained flapping) depending on the initial conditions. This numerical result we observed in computer simulation is the same as that of the laboratory experiment even though the Reynolds number of the computations is lower than that of the laboratory experiment by two orders of magnitude.",
keywords = "Bistability, Computational fluid dynamics, Flag-in-wind, Flapping filament, Immersed boundary method, Multigrid method",
author = "Luoding Zhu and Charles Peskin",
year = "2002",
month = "7",
day = "1",
doi = "10.1006/jcph.2002.7066",
language = "English (US)",
volume = "179",
pages = "452--468",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "2",

}

TY - JOUR

T1 - Simulation of a flapping flexible filament in a flowing soap film by the immersed boundary method

AU - Zhu, Luoding

AU - Peskin, Charles

PY - 2002/7/1

Y1 - 2002/7/1

N2 - This paper reports the computer simulation of a flapping flexible filament in a flowing soap film using the immersed boundary method. Our mathematical formulation includes filament mass and elasticity, gravity, air resistance, and the two wires that bound the flowing soap film. The incompressible viscous Navier-Stokes equations, which are used to describe the motion of the soap film and filament in our formulation, are discretized on a fixed uniform Eulerian lattice while the filament equations are discretized on a moving Lagrangian array of points which do not necessarily coincide with the fixed Eulerian mesh points of the fluid computation. The interaction between the filament and the soap film is handled by a smoothed approximation to the Dirac delta function. This delta function approximation is used not only to interpolate the fluid velocity and to apply force to the fluid (as is commonly done in immersed boundary computations), but also to handle the mass of the filament, which is represented in our calculation as delta function layer of fluid mass density supported along the immersed filament. Because of this nonuniform density, we need to use a multigrid method for solving the discretized fluid equations. This replaces the FFT-based method that is commonly used in the uniform-density case. Our main results are as follows. (i) The sustained flapping of the filament only occurs when filament mass is included in the formulation of the model; within a certain range of mass, the more the mass of the filament, the bigger the amplitude of the flapping. (ii) When the length of filament is short enough (below some critical length), the filament always approaches its straight (rest) state, in which the filament points downstream; but when the length is larger, the system is bistable, which means that it can settle into either state (rest state or sustained flapping) depending on the initial conditions. This numerical result we observed in computer simulation is the same as that of the laboratory experiment even though the Reynolds number of the computations is lower than that of the laboratory experiment by two orders of magnitude.

AB - This paper reports the computer simulation of a flapping flexible filament in a flowing soap film using the immersed boundary method. Our mathematical formulation includes filament mass and elasticity, gravity, air resistance, and the two wires that bound the flowing soap film. The incompressible viscous Navier-Stokes equations, which are used to describe the motion of the soap film and filament in our formulation, are discretized on a fixed uniform Eulerian lattice while the filament equations are discretized on a moving Lagrangian array of points which do not necessarily coincide with the fixed Eulerian mesh points of the fluid computation. The interaction between the filament and the soap film is handled by a smoothed approximation to the Dirac delta function. This delta function approximation is used not only to interpolate the fluid velocity and to apply force to the fluid (as is commonly done in immersed boundary computations), but also to handle the mass of the filament, which is represented in our calculation as delta function layer of fluid mass density supported along the immersed filament. Because of this nonuniform density, we need to use a multigrid method for solving the discretized fluid equations. This replaces the FFT-based method that is commonly used in the uniform-density case. Our main results are as follows. (i) The sustained flapping of the filament only occurs when filament mass is included in the formulation of the model; within a certain range of mass, the more the mass of the filament, the bigger the amplitude of the flapping. (ii) When the length of filament is short enough (below some critical length), the filament always approaches its straight (rest) state, in which the filament points downstream; but when the length is larger, the system is bistable, which means that it can settle into either state (rest state or sustained flapping) depending on the initial conditions. This numerical result we observed in computer simulation is the same as that of the laboratory experiment even though the Reynolds number of the computations is lower than that of the laboratory experiment by two orders of magnitude.

KW - Bistability

KW - Computational fluid dynamics

KW - Flag-in-wind

KW - Flapping filament

KW - Immersed boundary method

KW - Multigrid method

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

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

U2 - 10.1006/jcph.2002.7066

DO - 10.1006/jcph.2002.7066

M3 - Article

VL - 179

SP - 452

EP - 468

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 2

ER -