Direct-numerical simulation of the glottal jet and vocal-fold dynamics in a three-dimensional laryngeal model

An immersed-boundary method based flow solver coupled with a finite-element solid dynamics solver is employed in order to conduct direct-numerical simulations of phonatory dynamics in a three-dimensional model of the human larynx. The computed features of the glottal flow including mean and peak flow rates, and the open and skewness quotients are found to be within the normal physiological range. The flow-induced vibration pattern shows the classical "convergent-divergent" glottal shape, and the vibration amplitude is also found to be typical for human phonation. The vocal fold motion is analyzed through the method of empirical eigenfunctions and this analysis indicates a 1:1 modal entrainment between the "adduction-abduction" mode and the "mucosal wave" mode. The glottal jet is found to exhibit noticeable cycle-to-cycle asymmetric deflections and the mechanism underlying this phenomenon is examined.