Finite element approximation of flow induced vibrations of human vocal folds model: Effects of inflow boundary conditions and the length of subglottal and supraglottal channel on phonation onset

The paper presents the numerical analysis of interaction of the vibrating simplified human vocal folds model with the incompressible viscous airflow in a channel modeling simplified subglottal and supraglottal spaces. The flow in the considered 2D computational fluid domain is governed by the Navier–Stokes equations written in the Arbitrary Lagrangian–Eulerian form. The stabilized finite element method is applied for numerical approximation and the choice of boundary conditions and their implementation is discussed. For the considered model problem the prescribed inlet velocity and prescribed pressure difference formulations were numerically analyzed. The prescribed inlet velocity formulation was successful in predicting of the flutter velocity value, whereas the prescribed pressure difference gave nonphysical results. Finally a modified inlet boundary condition motivated by the penalization approach is suggested. It is shown that this approach gives possibilities to optimize the inlet boundary condition related to a physical reality by changing smoothly the penalty parameter in the interval between the two extremes and to treat the complete closures of the channel.

[1]  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.

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

[3]  C. Simerska,et al.  Modelling of aeroelastic behaviour of a plate in channel , 1996 .

[4]  John W. Barrett,et al.  Finite element approximation of the Dirichlet problem using the boundary penalty method , 1986 .

[5]  I. Babuska.  The Finite Element Method with Penalty , 1973 .

[6]  Jack J Jiang,et al.  Asymmetric airflow and vibration induced by the Coanda effect in a symmetric model of the vocal folds. , 2007, The Journal of the Acoustical Society of America.

[7]  Qian Xue,et al.  Subject-specific computational modeling of human phonation. , 2014, The Journal of the Acoustical Society of America.

[8]  I. Titze.  The myoelastic aerodynamic theory of phonation , 2006 .

[9]  J. Szmelter.  Incompressible flow and the finite element method , 2001 .

[10]  Tomáš Vampola,et al.  Development of Self-oscillating Human Vocal Folds Prosthesis , 2016 .

[11]  Fumio Inada,et al.  Flow Induced Vibrations : Classifications and Lessons from Practical Experiences , 2008 .

[12]  Maxim A. Olshanskii,et al.  Stabilized finite element schemes with LBB-stable elements for incompressible flows , 2005 .

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

[14]  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 .

[15]  Thomas J. R. Hughes,et al.  Weak imposition of Dirichlet boundary conditions in fluid mechanics , 2007 .

[16]  Thomas J. R. Hughes,et al.  An arbitrary Lagrangian-Eulerian finite rigid element method for interaction of fluid and a rigid body , 1992 .

[17]  Miloslav Feistauer,et al.  Application of a Stabilized FEM to Problems of Aeroelasticity , 2004 .

[18]  J. Flanagan,et al.  Synthesis of voiced sounds from a two-mass model of the vocal cords , 1972 .

[19]  Rolf Rannacher,et al.  ARTIFICIAL BOUNDARIES AND FLUX AND PRESSURE CONDITIONS FOR THE INCOMPRESSIBLE NAVIER–STOKES EQUATIONS , 1996 .

[20]  Petr Sváček,et al.  Numerical Simulation of Glottal Flow in Interaction with Self Oscillating Vocal Folds: Comparison of Finite Element Approximation with a Simplified Model , 2012 .

[21]  R. Mittal,et al.  A Computational Study of the Effect of False Vocal Folds on Glottal Flow and Vocal Fold Vibration During Phonation , 2009, Annals of Biomedical Engineering.

[22]  P. Tallec,et al.  Fluid structure interaction with large structural displacements , 2001 .

[23]  P. Fabrie,et al.  EFFECTIVE DOWNSTREAM BOUNDARY CONDITIONS FOR INCOMPRESSIBLE NAVIER-STOKES EQUATIONS , 1994 .

[24]  Vít Dolejsí.  Anisotropic mesh adaptation technique for viscous flow simulation , 2001, J. Num. Math..

[25]  Isao Tokuda,et al.  Mathematical Models and Numerical Schemes for the Simulation of Human Phonation , 2011 .

[26]  Jaromír Horáček,et al.  AEROELASTIC MODEL OF VOCAL-FOLD-SHAPED VIBRATING ELEMENT FOR STUDYING THE PHONATION THRESHOLD , 2002 .

[27]  R. Codina.  Stabilization of incompressibility and convection through orthogonal sub-scales in finite element methods , 2000 .

[28]  Miloslav Feistauer,et al.  Mathematical Methods in Fluid Dynamics , 1993 .

[29]  Rajat Mittal,et al.  An immersed-boundary method for flow-structure interaction in biological systems with application to phonation , 2008, J. Comput. Phys..

[30]  J. Horáček,et al.  Numerical simulation of flow induced airfoil vibrations with large amplitudes , 2007, Programs and Algorithms of Numerical Mathematics 19.

[31]  M. Döllinger,et al.  A 2D finite-element scheme for fluid–solid–acoustic interactions and its application to human phonation , 2009 .

[32]  Luc Mongeau,et al.  Aerodynamic transfer of energy to the vocal folds. , 2005, The Journal of the Acoustical Society of America.