Research Papers

A Coupled Sharp-Interface Immersed Boundary-Finite-Element Method for Flow-Structure Interaction With Application to Human Phonation

[+] Author and Article Information
X. Zheng, Q. Xue

Department of Mechanical Engineering, Johns Hopkins University, Baltimore, MD 21218

R. Mittal1

Department of Mechanical Engineering, Johns Hopkins University, Baltimore, MD 21218mittal@jhu.edu

S. Beilamowicz

Division of Otolaryngology, George Washington University, Washington, DC 20052


Corresponding author.

J Biomech Eng 132(11), 111003 (Oct 15, 2010) (12 pages) doi:10.1115/1.4002587 History: Received September 24, 2009; Revised September 02, 2010; Posted September 21, 2010; Published October 15, 2010; Online October 15, 2010

A new flow-structure interaction method is presented, which couples a sharp-interface immersed boundary method flow solver with a finite-element method based solid dynamics solver. The coupled method provides robust and high-fidelity solution for complex flow-structure interaction (FSI) problems such as those involving three-dimensional flow and viscoelastic solids. The FSI solver is used to simulate flow-induced vibrations of the vocal folds during phonation. Both two- and three-dimensional models have been examined and qualitative, as well as quantitative comparisons, have been made with established results in order to validate the solver. The solver is used to study the onset of phonation in a two-dimensional laryngeal model and the dynamics of the glottal jet in a three-dimensional model and results from these studies are also presented.

Copyright © 2010 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Figure 1

2D schematic describing ghost-cell methodology

Grahic Jump Location
Figure 2

2D illustration of contact force. The upper vocal fold intrudes into the center contact plane.

Grahic Jump Location
Figure 3

Schematic flow past an elastically mounted sphere

Grahic Jump Location
Figure 4

(a) A coronal view of CT scan of human larynx and the current flow domain that attempts to match the key geometrical features in the CT scan. (b) Three-layer vocal fold inner structure inside an idealized geometric VF model based on CT scan and triangular elements used in the current solver.

Grahic Jump Location
Figure 5

Grid used in the current simulation. The figure shows every three grid points in each direction.

Grahic Jump Location
Figure 6

Variation of fundamental frequency with subglottal pressure

Grahic Jump Location
Figure 7

Time variation of two-dimensional glottal volume fluxes in the stationary vibration stage

Grahic Jump Location
Figure 8

Contours of spanwise vorticity during vocal fold vibration: (a) 0.3384 s, (b) 0.3398 s, (c) 0.3409 s, and (d) 0.3454 s

Grahic Jump Location
Figure 9

Contours of stresses (kPa) in the vocal folds during the open and closed phases of the vibration cycle

Grahic Jump Location
Figure 10

Flow domain in the 3D flow-structure simulation and finite-element mesh for the true vocal folds

Grahic Jump Location
Figure 11

Time variation of three-dimensional glottal volume flux in the stationary vibration stage

Grahic Jump Location
Figure 12

Isosurface of swirl strength at six different time instants over one vocal fold vibration cycle: (a) 0.0245 s, (b) 0.02625 s, (c) 0.02695 s, (d) 0.02765 s, (e) 0.02800 s, and (f) 0.02870 s

Grahic Jump Location
Figure 13

Isosurface of the turbulent kinetic energy corresponding to a value of KT=0.15 and a contour of turbulent kinetic energy at the center plane (z=0.75 cm)

Grahic Jump Location
Figure 14

(a) Five locations along the center line in X-Y plane chosen to perform the flow spectrum analysis. ((b)–(f)) Span-averaged streamwise velocity spectra: (b) point 1 (x=3.025 cm), (c) point 2 (x=3.175 cm), (d) point 3 (x=3.325 cm), (e) point 4 (x=3.575 cm), and (f) point 5 (x=3.975 cm). Dash-dot line -⋅-⋅-⋅-⋅-⋅- corresponds to k−5/3.



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In