Simulation of unsteady compressible flow in a channel with vibrating walls - Influence of the frequency

Abstract This study deals with the numerical solution of a 2D unsteady flow of a compressible viscous fluid in a channel for low inlet airflow velocity. The unsteadiness of the flow is caused by a prescribed periodic motion of a part of the channel wall with large amplitudes, nearly closing the channel during oscillations. The channel is a simplified model of the glottal space in the human vocal tract and the flow can represent a model of airflow coming from the trachea, through the glottal region with periodically vibrating vocal folds to the human vocal tract. The flow is described by the system of Navier–Stokes equations for laminar flows. The numerical solution is implemented using the finite volume method (FVM) and the predictor–corrector MacCormack scheme with Jameson artificial viscosity using a grid of quadrilateral cells. Due to the motion of the grid, the basic system of conservation laws is considered in the Arbitrary Lagrangian–Eulerian (ALE) form. The authors present the numerical simulations of flow fields in the channel, acquired from a program developed exclusively for this purpose. The numerical results for unsteady flows in the channel are presented for inlet Mach number M ∞  = 0.012, Reynolds number Re ∞  = 4.5 × 10 3 and the wall motion frequency 20 and 100 Hz.

[1]  Jiří Fürst,et al.  Finite Volume Solution of 2D and 3D Euler and Navier-Stokes Equations , 2001 .

[2]  Jiří Fürst,et al.  An unsteady numerical solution of viscous compressible flows in a channel , 2006 .

[3]  H. K. Schutte,et al.  Glottal flow through a two-mass model: comparison of Navier-Stokes solutions with simplified models. , 2002, Journal of the Acoustical Society of America.

[4]  Manfred Kaltenbacher,et al.  Human phonation analysis by 3d aero-acoustic computation , 2009 .

[5]  van Rr René Hassel,et al.  Theoretical and experimental study of quasisteady‐flow separation within the glottis during phonation. Application to a modified two‐mass model , 1994 .

[6]  A. Alwan,et al.  A contribution to simulating a three-dimensional larynx model using the finite element method. , 2003, The Journal of the Acoustical Society of America.

[7]  Zhaoyan Zhang,et al.  Coherent structures of the near field flow in a self-oscillating physical model of the vocal folds. , 2007, The Journal of the Acoustical Society of America.

[8]  Jaromír Horácek,et al.  Numerical simulation of airflow through the oscillating glottis , 2007, MAVEBA.

[9]  Petr Šidlof,et al.  Numerical simulation of self-oscillations of human vocal folds with Hertz model of impact forces , 2005 .

[10]  P. Šidlof Fluid-structure interaction in human vocal folds , 2007 .

[11]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[12]  K. Kozel,et al.  Remark on numerical simulation of 2D unsteady transonic flows , 1993 .

[13]  Ingo R. Titze,et al.  Combined simulation of two‐dimensional airflow and vocal‐fold vibration , 1997 .