Identification of geometric parameters influencing the flow-induced vibration of a two-layer self-oscillating computational vocal fold model

Simplified models have been used to simulate and study the flow-induced vibrations of the human vocal folds. While it is clear that the models' responses are sensitive to geometry, it is not clear how and to what extent specific geometric features influence model motion. In this study geometric features that played significant roles in governing the motion of a two-layer (body-cover), two-dimensional, finite element vocal fold model were identified. The model was defined using a flow solver based on the viscous, unsteady, Navier-Stokes equations and a solid solver that allowed for large strain and deformation. A screening-type design-of-experiments approach was used to identify the relative importance of 13 geometric parameters. Five output measures were analyzed to assess the magnitude of each geometric parameter's effect on the model's motion. The measures related to frequency, glottal width, flow rate, intraglottal angle, and intraglottal phase delay. The most significant geometric parameters were those associated with the cover--primarily the pre-phonatory intraglottal angle--as well as the body inferior angle. Some models exhibited evidence of improved model motion, including mucosal wave-like motion and alternating convergent-divergent glottal profiles, although further improvements are still needed to more closely mimic human vocal fold motion.