Computational modeling of vibration-induced systemic hydration of vocal folds over a range of phonation conditions

Predicting phonation conditions that are benign to voice health remains a biomechanically relevant problem. Our objective is to provide insight into vocal fold (VF) hydration based on continuum-based VF models that are able to compute VF stresses during phonation and a scheme for the extraction and generalization of such computational data based on the principle of linear superposition. Because VF tissue is poroelastic, spatial gradients of VF hydrostatic stresses computed for a given phonation condition determine VF interstitial fluid flow. The present approach transforms, based on linear superposition principles, the computed interstitial fluid velocities at the particular phonation to those at an arbitrary phonation condition. Intersititial fluid flow characteristics for a range of phonation conditions are compared. For phonation conditions with no or moderate collision, no dehydration per vibration cycle is predicted throughout the VF. For more severe collision conditions, tissue dehydration is restricted to a region close to the glottal surface. Interstitial fluid displacement in the VF is found to be heterogeneous and strongly dependent on the phonation condition. A phonation condition is found to exist for which dehydration peaks. The proposed method significantly expands the scope and relevance of conducting isolated numerical simulations of VF vibration.