Glottal flow through a two-mass model: comparison of Navier-Stokes solutions with simplified models

A new numerical model of the vocal folds is presented based on the well-known two-mass models of the vocal folds. The two-mass model is coupled to a model of glottal airflow based on the incompressible Navier-Stokes equations. Glottal waves are produced using different initial glottal gaps and different subglottal pressures. Fundamental frequency, glottal peak flow, and closed phase of the glottal waves have been compared with values known from the literature. The phonation threshold pressure was determined for different initial glottal gaps. The phonation threshold pressure obtained using the flow model with Navier-Stokes equations corresponds better to values determined in normal phonation than the phonation threshold pressure obtained using the flow model based on the Bernoulli equation. Using the Navier-Stokes equations, an increase of the subglottal pressure causes the fundamental frequency and the glottal peak flow to increase, whereas the fundamental frequency in the Bernoulli-based model does not change with increasing pressure.     

The effects of rehydration on phonation in excised canine larynges

Experiments using excised canine larynges were conducted to study the restoration of vocal efficiency in dehydrated larynges. Excised larynges were dehydrated with warm, dry air to the point that airflow through the approximated vocal folds would not entrain the folds to produce phonation. The dehydrated vocal folds were then bathed in a saline solution. The rehydrated larynges were then remounted on the bench apparatus that enabled phonation with a constant humidified airflow, and measurements were made of phonation threshold pressure, glottal airflow, and amplitude. Hydration resulted in significantly increased efficiency and decrease in phonation threshold pressure. The findings confirm clinical impressions that hydration is critical in the physiology of normal phonation.     

Optimal glottal configuration for ease of phonation

Recent experimental studies have shown the existence of optimalvalues of the glottal width and convergence angle, at which the phonation threshold pressure is minimum. These results indicate the existence of an optimal glottal configuration for ease of phonation, not predicted by the previous theory. In this paper, the origin of the optimal configuration is investigated using a low dimensional mathematical model of the vocal fold. Two phenomena of glottal aerodynamics are examined: pressure losses due to air viscosity, and air flow separation from a divergent glottis. The optimal glottal configuration seems to be a consequence of the combined effect of both factors. The results agree with the experimental data, showing that the phonation threshold pressure is minimum when the vocal folds are slightly separated in a near rectangular glottis.     

Phonation Threshold Power in Ex Vivo Laryngeal Models

This study hypothesized that phonation threshold power is measureable and sensitive to changes in the biomechanical properties of the vocal folds. Phonation threshold power was measured in three sample populations of 10 excised canine larynges treated with variable posterior glottal gap, variable bilateral vocal fold elongation, and variable vocal fold lesioning. Posterior glottal gap varied from 0 to 4 mm in 0.5 mm intervals. Bilateral vocal fold elongation varied from 0% to 20% in 5% intervals. Vocal fold lesion treatments included unilateral and bilateral vocal fold lesion groups. Each treatment was investigated independently in a sample population of 10 excised canine larynges. Linear regression analysis indicated that phonation threshold power was sensitive to posterior glottal gap (R² = 0.298, P < 0.001) and weakly to vocal fold elongation (R² = 0.052, P = 0.003). A one-way repeated measures analysis of variance indicated that phonation threshold power was sensitive to the presence of lesions (P < 0.001). Theoretical and experimental evidence presented here suggests that phonation threshold power could be used as a broad screening parameter sensitive to certain changes in the biomechanical properties of the larynx. It has not yet been measured in humans, but because it has the potential to represent the airflow-tissue energy transfer more completely than the phonation threshold pressure or flow alone, it may be a more useful parameter than these and could be used to indicate that laryngeal health is likely abnormal.     

Influence of supraglottal structures on the glottal jet exiting a two-layer synthetic, self-oscillating vocal fold model

A synthetic two-layer, self-oscillating, life-size vocal fold model was used to study the influence of the vocal tract and false folds on the glottal jet. The model vibrated at frequencies, pressures, flow rates, and amplitudes consistent with human phonation, although some differences in behavior between the model and the human vocal folds are noted. High-speed images of model motion and flow visualization were acquired. Phase-locked ensemble-averaged glottal jet velocity measurements using particle image velocimetry (PIV) were acquired with and without an idealized vocal tract, with and without false folds. PIV data were obtained with varying degrees of lateral asymmetric model positioning. Glottal jet velocity magnitudes were consistent with those measured using excised larynges. A starting vortex was observed in all test cases. The false folds interfered with the starting vortex, and in some cases vortex shedding from the false folds was observed. In asymmetric cases without false folds, the glottal jet tended to skew toward the nearest wall; with the false folds, the opposite trend was observed. rms velocity calculations showed the jet shear layer and laminar core. The rms velocities were higher in the vocal tract cases compared to the open jet and false fold cases.     

Effect of inferior surface angle on the self-oscillation of a computational vocal fold model

Geometry of the human vocal folds strongly influences their oscillatory motion. While the effect of intraglottal geometry on phonation has been widely investigated, the study of the geometry of the inferior surface of the vocal folds has been limited. In this study the way in which the inferior vocal fold surface angle affects vocal fold vibration was explored using a two-dimensional, self-oscillating finite element vocal fold model. The geometry was parameterized to create models with five different inferior surface angles. Four of the five models exhibited self-sustained oscillations. Comparisons of model motion showed increased vertical displacement and decreased glottal width amplitude with decreasing inferior surface angle. In addition, glottal width and air flow rate waveforms changed as the inferior surface angle was varied. Structural, rather than aerodynamic, effects are shown to be the cause of the changes in model response as the inferior surface angle was varied. Supporting data including glottal pressure distribution, average intraglottal pressure, energy transfer, and flow separation point locations are discussed, and suggestions for future research are given.     

Dynamic glottal pressures in an excised hemilarynx model

During phonation, air pressures act upon the vocal folds to help maintain their oscillation. The air pressures vary dynamically along the medial surface of the vocal folds, although no live human or excised studies have shown how those pressure profiles vary in time. The purpose of this study was to examine time-dependent glottal pressure profiles using a canine hemilarynx approach. The larynx tissue was cut in the midsaggital plane from the top to about 5 mm below the vocal folds. The right half was replaced with a Plexiglas pane with imbedded pressure taps. Simultaneous recordings were made of glottal pressure signals, subglottal pressure, particle velocity, and average airflow at various levels of adduction. The data indicate that the pressures in the glottis (on the Plexiglas) vary both vertically and longitudinally throughout the phonatory cycle. Pressures vary most widely near the location of maximum vibratory amplitude, and can include negative pressures during a portion of the cycle. Pressures anterior and posterior to the maximum amplitude location may have less variation and may remain positive throughout the cycle, giving rise to a new concept called dynamic bidirectional pressure gradients in the glottis. This is an important concept that may relate strongly to tissue health as well as basic oscillatory mechanics.     

Effect of a Single Injection of Basic Fibroblast Growth Factor into the Vocal Folds: A 36-Month Clinical Study

ObjectiveA single injection of basic fibroblast growth factor (bFGF) into the vocal folds of patients with glottal insufficiency has been shown to be effective for a few years. However, the long-term therapeutic effect of a single injection of bFGF into the vocal folds has yet to be demonstrated. In this study, the therapeutic effect of a single injection of bFGF into the vocal folds was investigated over several years by monitoring patients for 36 months following this treatment.MethodsNineteen patients with glottal insufficiency received injections of bFGF diluted to 20 μg/mL in the superficial layer of the lamina propria of the bilateral vocal folds. The following parameters were evaluated at preinjection baseline and 6, 12, 18, 24, and 36 months later, and statistical comparisons were performed. The parameters evaluated were: the Grade, Rough, Breathy, Asthenic, and Strained (GRBAS) scale score; maximum phonation time; acoustic analysis; and glottal wave analysis (GWA) and kymograph edge analysis (KEA) using high-speed digital imaging (HSDI). The amplitude perturbation quotient (APQ) and period perturbation quotient (PPQ) were measured by acoustic analysis. The mean minimum glottal area during vocalization and mean minimum distance between the vocal folds were measured by GWA. The amplitudes of the bilateral vocal folds were measured by KEA.ResultsPostinjection, the GRBAS scale score decreased from 6 months after injection, and maximum phonation time was prolonged. The mean minimum glottal area during vocalization and the mean minimum distance between the vocal folds calculated by GWA of HSDI decreased significantly after 6 months. These effects persisted until 36 months postinjection. APQ and PPQ derived from acoustic analysis tended to decrease, but not significantly. There was no clear change in the amplitudes of the bilateral vocal folds calculated by KEA of HSDI before and after injection.ConclusionsThese results suggest that the effects of a single injection of bFGF into the vocal folds persist for 36 months.     

Pressure-flow relationships during phonation in the canine larynx

Measurements of air pressure and flow were made using an in vivo canine model of the larynx. Subglottic pressures at varying flow rates were taken during phonation induced by laryngeal nerve stimulation. Results showed that during constant vocal fold stiffness, subglottic pressure rose slightly with increased air flow. The larynx in the in vivo canine model exhibited a flow-dependent decrease in laryngeal airway resistance. Increasing flow rate was associated with an increase in frequency of phonation and open quotient, as measured glottographically. Results from this experiment were compared with a theoretical two-mass model of the larynx and other theoretical models of phonation. The influence of aerodynamic forces on glottal vibration is explained by increased lateral excursion of the vocal folds during the open interval and shortening of the closed interval during the glottal cycle.     

Vocal fold and ventricular fold vibration in period-doubling phonation: physiological description and aerodynamic modeling

Occurrences of period-doubling are found in human phonation, in particular for pathological and some singing phonations such as Sardinian A Tenore Bassu vocal performance. The combined vibration of the vocal folds and the ventricular folds has been observed during the production of such low pitch bass-type sound. The present study aims to characterize the physiological correlates of this acoustical production and to provide a better understanding of the physical interaction between ventricular fold vibration and vocal fold self-sustained oscillation. The vibratory properties of the vocal folds and the ventricular folds during phonation produced by a professional singer are analyzed by means of acoustical and electroglottographic signals and by synchronized glottal images obtained by high-speed cinematography. The periodic variation in glottal cycle duration and the effect of ventricular fold closing on glottal closing time are demonstrated. Using the detected glottal and ventricular areas, the aerodynamic behavior of the laryngeal system is simulated using a simplified physical modeling previously validated in vitro using a larynx replica. An estimate of the ventricular aperture extracted from the in vivo data allows a theoretical prediction of the glottal aperture. The in vivo measurements of the glottal aperture are then compared to the simulated estimations.     

The membranous contact quotient: A new phonatory measure of glottal competence

The membranous contact quotient (MCQ) is introduced as a measure of dynamic glottal competence. It is defined as the ratio of the membranous contact glottis (the anterior-posterior length of contact between the two membranous vocal folds) and the membranous vocal fold length. An elliptical approximation to the vocal fold contour during phonation was used to predict MCQ values as a function of vocal process gap (adduction), maximum glottal width, and membranous glottal length. MCQ is highly dependent on the vocal process gap and the maximum glottal width, but not on vocal fold length. Five excised larynges were used to obtain MCQ data for a wide range of vocal process gaps and maximum glottal widths. Predicted and measured MCQ values had a correlation of 0.93, with an average absolute difference of 9.6% (SD = 10.5%). The model is better at higher values of MCQ. The theory for MCQ is also expressed as a function of vocal process gap and subglottal pressure to suggest production control potential. The MCQ measure is obtainable with the use of stroboscopy and appears to be a potentially useful clinical measure.     

Phonation threshold pressure and onset frequency in a two-layer physical model of the vocal folds

The influence of vocal fold geometry and stiffness on phonation onset was experimentally investigated using a body-cover physical model of the vocal folds. Results showed that a lower phonation threshold pressure and phonation onset frequency can be achieved by reducing body-layer or cover-layer stiffness, reducing medial surface thickness, or increasing cover-layer depth. Increasing body-layer stiffness also restricted vocal fold motion to the cover layer and reduced prephonatory glottal opening. Excitation of anterior-posterior modes was also observed, particularly for large values of the body-cover stiffness ratio. The results of this study were also discussed in relation to previous theoretical and experimental studies.     

The effect of viscosity changes in the vocal folds on the range of oscillation

Changes in vocal fold oscillation threshold pressure were induced in excised canine larynges by experimentally causing fluid movement into and out of the vocal folds. The transport was facilitated by exposing the vocal folds to various osmotic solutions, and it was assumed that changes in hydration caused changes in the internal tissue viscosity. A range of oscillation threshold pressures was measured for each condition of hydration by varying length and glottal width. The oscillation threshold pressure shifted as predicted. Decreased hydration (increased viscosity) raised the threshold of oscillation, and increased hydration (decreased viscosity) lowered the threshold of oscillation. This apparently represents the first in vitro model for the study of the effect of viscosity changes of the internal environment of the vocal folds on phonation.     

Simulation of vocal fold impact pressures with a self-oscillating finite-element model

Vocal fold impact pressures were studied using a self-oscillating finite-element model capable of simulating vocal fold vibration and airflow. The calculated airflow pressure is applied on the vocal fold as the driving force. The airflow region is then adjusted according to the calculated vocal fold displacement. The interaction between airflow and the vocal folds produces a self-oscillating solution. Lung pressures between 0.2 and 2.5 kPa were used to drive this self-oscillating model. The spatial distribution of the impact pressure was studied. Studies revealed that the tissue collision during phonation produces a very large impact pressure which correlates with the lung pressure and glottal width. Larger lung pressure and a narrower glottal width increase the impact pressure. The impact pressure was found to be roughly the square root of lung pressure. In the inferior-superior direction, the maximum impact pressure is related to the narrowest glottis. In the anterior-posteriorfirection, the greatest impact pressure appears at the midpoint of the vocal fold. The match between our numerical simulations and clinical observations suggests that this self-oscillating finite-element model might be valuable for predicting mechanical trauma of the vocal folds.     

On the difference between negative damping and eigenmode synchronization as two phonation onset mechanisms

Negative damping and eigenmode synchronization as two different mechanisms of phonation onset are distinguished. Although both mechanisms lead to a favorable phase relationship between the flow pressure and the vocal fold motion as required for a net energy transfer into the vocal folds, the underlying mechanisms for this favorable phase relationship are different. The negative damping mechanism relies on glottal aerodynamics or acoustics to establish before onset and maintain after onset the favorable phase relationship, and therefore has minimum requirements on vocal fold geometry and biomechanics. A single degree-of-freedom vocal fold model is all that is needed for self-oscillation in the presence of a negative damping mechanism. In contrast, the mechanism of eigenmode synchronization critically depends on the geometrical and biomechanical properties of the vocal folds (at least 2-degrees-of-freedom are required), and has little requirement on the glottal aerodynamics other than flow separation. The favorable phase relation is established once synchronization occurs, regardless of the phase relationship imposed by glottal aerodynamics before onset. Unlike that of the negative damping mechanism, initiation of eigenmode synchronization requires neither a velocity-dependent flow pressure nor an alternating convergent-divergent glottis. The clinical implications of the distinctions between these two mechanisms are discussed.     

In Vivo Quantification of the Intraglottal Pressure: Modal Phonation and Voice Onset

Intraglottal pressure is the driving force of vocal fold vibration. Its time course during the open phase of the vibratory cycle is essential in the mechanics of phonation, but measuring it directly is difficult and may hinder spontaneous voicing. However, it can be computed from the in vivo measured transglottal flow and glottal area (hence the air particle velocity) on the basis of the Bernoulli energy law and the interaction with the inertance of the vocal tract. As to sustained modal phonation, calculations are presented for the two possible shapes of glottal duct: convergent and divergent, including absolute calibration in order to obtain quantitative physical values. Whatever the glottal duct configuration, the calculations based on measured values of glottal area and air flow show that the integrated intraglottal pressure during the opening phase systematically exceeds that during the closing phase, which is the basic condition for sustaining vocal fold oscillation. The key point is that the airflow curve is skewed to the right relative to the glottal area curve. The skewing results from air compressibility and vocal tract inertance. The intraglottal pressure becomes negative during the closing phase. As to the soft (or physiological) voice onset, a similar approach shows that the integrated pressure differences (opening phase − closing phase) actually increase as the onset progresses, and this applies to the results based on Bernoulli's energy law as well as to those based on the interaction with the inertance of the vocal tract. Furthermore and similarly, the phase lead of the pressure wave with respect to the glottal opening progressively increases. The underlying explanation lies in the progressively increasing skewing of the airflow curve to the right with respect to the glottal area curve.     

Theoretical simulation and experimental validation of inverse quasi-one-dimensional steady and unsteady glottal flow models

In physical modeling of phonation, the pressure drop along the glottal constriction is classically assessed with the glottal geometry and the subglottal pressure as known input parameters. Application of physical modeling to study phonation abnormalities and pathologies requires input parameters related to in vivo measurable quantities commonly corresponding to the physical model output parameters. Therefore, the current research presents the inversion of some popular simplified flow models in order to estimate the subglottal pressure, the glottal constriction area, or the separation coefficient inherent to the simplified flow modeling for steady and unsteady flow conditions. The inverse models are firstly validated against direct simulations and secondly against in vitro measurements performed for different configurations of rigid vocal fold replicas mounted in a suitable experimental setup. The influence of the pressure corrections related to viscosity and flow unsteadiness on the flow modeling is quantified. The inversion of one-dimensional glottal flow models including the major viscous effects can predict the main flow quantities with respect to the in vitro measurements. However, the inverse model accuracy is strongly dependent on the pertinence of the direct flow modeling. The choice of the separation coefficient is preponderant to obtain pressure predictions relevant to the experimental data.     

Phonation threshold pressure: comparison of calculations and measurements taken with physical models of the vocal fold mucosa

In an important paper on the physics of small amplitude oscillations, Titze showed that the essence of the vertical phase difference, which allows energy to be transferred from the flowing air to the motion of the vocal folds, could be captured in a surface wave model, and he derived a formula for the phonation threshold pressure with an explicit dependence on the geometrical and biomechanical properties of the vocal folds. The formula inspired a series of experiments [e.g., R. Chan and I. Titze, J. Acoust. Soc. Am 119, 2351-2362 (2006)]. Although the experiments support many aspects of Titze's formula, including a linear dependence on the glottal half-width, the behavior of the experiments at the smallest values of this parameter is not consistent with the formula. It is shown that a key element for removing this discrepancy lies in a careful examination of the properties of the entrance loss coefficient. In particular, measurements of the entrance loss coefficient at small widths done with a physical model of the glottis (M5) show that this coefficient varies inversely with the glottal width. A numerical solution of the time-dependent equations of the surface wave model shows that adding a supraglottal vocal tract lowers the phonation threshold pressure by an amount approximately consistent with Chan and Titze's experiments.