Numerical simulation of turbulence transition and sound radiation for flow through a rigid glottal model

Large eddy simulation (LES)-based computational aeroacoustics techniques were applied to a static model of the human glottis, idealized here as a planar channel with an orifice, to study flow-acoustic interactions related to speech. Rigid models of both converging and diverging glottal passages, each featuring a 20 deg included angle and a minimal glottal diameter of 0.04 cm, with an imposed transglottal pressure of 15 cm H2O, were studied. The Favre-filtered compressible Navier-Stokes equations were integrated for this low-Mach-number flow using an additive semi-implicit Runge-Kutta method and a high-order compact finite-difference scheme with characteristic-based nonreflecting boundary conditions and a multiblock approach. Flow asymmetries related to the Coanda effect and transition to turbulence, as well as the far-field sound, were captured. Acoustic-analogy-based far-field sound predictions were compared with direct simulations and showed that dipole sources, arising from unsteady flow forces exerted on the glottal walls, are primarily responsible for the tonal sound observed in the divergent glottis case.