Flow visualization and pressure distributions in a model of the glottis with a symmetric and oblique divergent angle of 10 degrees

Modeling the human larynx can provide insights into the nature of the flow and pressures within the glottis. In this study, the intraglottal pressures and glottal jet flow were studied for a divergent glottis that was symmetric for one case and oblique for another. A Plexiglas model of the larynx (7.5 times life size) with interchangeable vocal folds was used. Each vocal fold had at least 11 pressure taps. The minimal glottal diameter was held constant at 0.04 cm. The glottis had an included divergent angle of 10 degrees. In one case the glottis was symmetric. In the other case, the glottis had an obliquity of 15 degrees. For each geometry, transglottal pressure drops of 3, 5, 10, and 15 cm H2O were used. Pressure distribution results, suggesting significantly different cross-channel pressures at glottal entry for the oblique case, replicate the data in another study by Scherer et al. [J. Acoust. Soc. Am. 109, 1616-1630 (2001b)]. Flow visualization using a LASER sheet and seeded airflow indicated separated flow inside the glottis. Separation points did not appear to change with flow for the symmetric glottis, but for the oblique glottis moved upstream on the divergent glottal wall as flow rate increased. The outgoing glottal jet was skewed off-axis for both the symmetric and oblique cases. The laser sheet showed asymmetric circulating regions in the downstream region. The length of the laminar core of the glottal jet was less than approximately 0.6 cm, and decreased in length as flow increased. The results suggest that the glottal obliquity studied here creates significantly different driving forces on the two sides of the glottis (especially at the entrance to the glottis), and that the skewed glottal jet characteristics need to be taken into consideration for modeling and aeroacoustic purposes.     

The effect of glottal angle on intraglottal pressure

Intraglottal pressure distributions depend upon glottal shape, size, and diameter. This study reports the effects of varying glottal angle on intraglottal and transglottal pressures using a three-dimensional Plexiglas model with a glottis having nine symmetric glottal angles and a constant minimal glottal diameter of 0.06 cm. The empirical data were supported by computational results using FLUENT. The results suggested that (1) the greater the convergent glottal angle, the greater outward driving forces (higher intraglottal pressures) on the vocal folds; (2) flow resistance was greatest for the uniform glottis, and least for the 10 degrees divergent glottis; (3) the greatest negative pressure in the glottis and therefore the greatest pressure recovery for diverging glottal shapes occurred for an angle of 10 degrees; (4) the smaller the convergent angle, the greater the flow resistance; (5) FLUENT was highly accurate in predicting the empirical pressures of this model; (6) flow separation locations (given by FLUENT) for the divergent glottis moved upstream for larger flows and larger glottal angles. The results suggest that phonatory efficiency related to aerodynamics may be enhanced with vocal fold oscillations that include large convergent angles during glottal opening and small (5 degrees - 10 degrees) divergent angles during glottal closing.     

The effect of entrance radii on intraglottal pressure distributions in the divergent glottis

Modeling laryngeal aerodynamics requires specification of the glottal geometry. Changing the glottal exit radius alters the intraglottal pressure distributions in the converging glottis [Scherer et al., J. Acoust. Soc. Am. 110, 2267-2269 (2001)]. This study examined the effects of the glottal entrance radius on the intraglottal pressure distributions for divergent angles of 5°, 10°, 20°, 30°, and 40°. Glottal airflow and minimal glottal diameter were held constant at 73.2 cm(3)/s and 0.02 cm, respectively. The computational code FLUENT was used to obtain the pressure distributions. Results suggest that a smaller glottal entrance radius tends to (1) lower the transglottal pressure (reduce glottal flow resistance), although this is angle dependent, (2) make the pressure dip near the glottal entrance more negative in value, (3) increase the slope of the pressure distribution just upstream of the glottal entrance, and (4) make the initial pressure recovery (rise) in the glottis steeper. A general empirical equation for transglottal pressure as a function of radius, angle, and separation point location is offered. These results suggest that glottal entrance curvature for the divergent glottis significantly affects the driving pressures on the vocal folds, and needs to be well specified when building computational and physical models.     

Unsteady behavior of flow in a scaled-up vocal folds model

Measurements of the fluid flow through a scaled-up model of the human glottis are presented to determine whether glottal flow may be approximated as unsteady. Time- and space-resolved velocity vector fields from digital particle image velocimetry (DPIV) measurements of the flow through the gap between two moving, rigid walls are presented in four cases, over a range of Strouhal numbers: 0.010, 0.018, 0.035, 0.040, corresponding to life-scale f(0) of 30, 58, 109, and 126 Hz, respectively, at a Reynolds number of 8000. It is observed that (1) glottal flow onset is delayed after glottal opening and (2) glottal flow shutoff occurs prior to closure. A comparison between flow through a fully open, nonmoving glottis and that through the moving vocal folds shows a marked difference in spatial structure of the glottal jet. The following features of the flow are seen to exhibit strong dependence on cycle frequency: (a) glottal exit plane velocity, (b) volume flow, (c) vortex shedding rates, and (d) vortex amplitude. Vortex shedding appears to be a factor both in controlling flow resistance and in cycle-to-cycle volume flow variations. All these observations strongly suggest that glottal flow is inherently unsteady.     

Theoretical assessment of unsteady aerodynamic effects in phonation

This paper ranks the importance of unsteady aerodynamic mechanisms in glottal flow. Particular emphasis is given to separation point motion, acceleration of glottal airflow by vocal fold motion, and viscous blockage. How nondimensional parameters such as the Reynolds, Strouhal, and Womersley numbers help in this ranking is also addressed. An equation of motion is derived which includes terms explicitly describing the effects of interest, assuming (1) a symmetrical glottis, (2) zero pressure recovery downstream of the vocal folds, and (3) a quasisteady glottal jet. Estimating the order of magnitude of the terms in this equation, it is shown that the flow is characterized by two temporal regimes: (1) a flow initiation/shutoff regime where local unsteady acceleration and wall motion dominate, and (2) a "quasisteady" regime where the flow is dominated by convective acceleration. In the latter case, separation point motion and viscous blockage are shown to be out of phase with motion of the vocal folds, thereby impacting the shape of the glottal volume flow waveform. The analysis suggests that glottal flow may be considered quasisteady only insofar as traditional assumptions concerning glottal jet behavior can be confirmed.     

Pressure and velocity profiles in a static mechanical hemilarynx model

This study examined pressure and velocity profiles in a hemilarynx mechanical model of phonation. The glottal section had parallel walls and was fabricated from hard plastic. Twelve pressure taps were created in the vocal fold surface and connected to a differential pressure transducer through a pressure switch. The glottal gap was measured with feeler gauges and the uniform glottal duct was verified by use of a laser system. Eight pressure transducers were placed in the flat wall opposite the vocal fold. Hot-wire anemometry was used to obtain velocity profiles upstream and downstream of the glottis. The results indicate that the pressure distribution on the vocal fold surface was consistent with pressure change along a parallel duct, whereas the pressures on the opposite flat wall typically were lower (by 8%-40% of the transglottal pressure just past mid-glottis). The upstream velocity profiles were symmetric regardless of the constriction shape and size. The jet flow downstream of the glottis was turbulent even for laminar upstream conditions. The front of the jet was consistently approximately 1.5 mm from the flat wall for glottal gaps of 0.4, 0.8 and 1.2 mm. The turbulence intensity also remained approximately at the same location of about 4 mm from the flat wall for the two larger gaps.     

Gender recognition from vocal source

Efficiency of automatic recognition of male and female voices based on solving the inverse problem for glottis area dynamics and for waveform of the glottal airflow volume velocity pulse is studied. The inverse problem is regularized through the use of analytical models of the voice excitation pulse and of the dynamics of the glottis area, as well as the model of one-dimensional glottal airflow. Parameters of these models and spectral parameters of the volume velocity pulse are considered. The following parameters are found to be most promising: the instant of maximum glottis area, the maximum derivative of the area, the slope of the spectrum of the glottal airflow volume velocity pulse, the amplitude ratios of harmonics of this spectrum, and the pitch. On the plane of the first two main components in the space of these parameters, an almost twofold decrease in the classification error relative to that for the pitch alone is attained. The male voice recognition probability is found to be 94.7%, and the female voice recognition probability is 95.9%.     

Two-dimensional model of vocal fold vibration for sound synthesis of voice and soprano singing

Voiced sounds were simulated with a computer model of the vocal fold composed of a single mass vibrating both parallel and perpendicular to the airflow. Similarities with the two-mass model are found in the amplitudes of the glottal area and the glottal volume flow velocity, the variation in the volume flow waveform with the vocal tract shape, and the dependence of the oscillation amplitude upon the average opening area of the glottis, among other similar features. A few dissimilarities are also found in the more symmetric glottal and volume flow waveforms in the rising and falling phases. The major improvement of the present model over the two-mass model is that it yields a smooth transition between oscillations with an inductive load and a capacitive load of the vocal tract with no sudden jumps in the vibration frequency. Self-excitation is possible both below and above the first formant frequency of the vocal tract. By taking advantage of the wider continuous frequency range, the two-dimensional model can successfully be applied to the sound synthesis of a high-pitched soprano singing, where the fundamental frequency sometimes exceeds the first formant frequency.     

Experimental investigation of the influence of a posterior gap on glottal flow and sound

The influence of a posterior gap on the airflow through the human glottis was investigated using a driven synthetic model. Instantaneous orifice discharge coefficient of a glottal shaped orifice was obtained from the time-varying orifice area and the velocity distribution of the pulsated jet measured on the axial plane using a single hot-wire probe. Instantaneous orifice discharge coefficient values were found to undergo a cyclic hysteresis loop when plotted versus Reynolds number and time, indicating a pressure head increase and a net energy transfer from the air flow to the orifice wall. The net energy transferred was estimated to be around 10% of the value presumably required to achieve self-sustained oscillation. The radiated sound pressure was measured to characterize the influence of the minimum flow through the posterior gap on the broadband component of the radiated sound. The presence of a posterior gap was found to significantly increase the broadband sound level produced over the frequency range in which human hearing is most sensitive.     

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.     

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.     

Intraglottal pressure profiles for a symmetric and oblique glottis with a divergence angle of 10 degrees

Human phonation does not always involve symmetric motions of the two vocal folds. Asymmetric motions can create slanted or oblique glottal angles. This study reports intraglottal pressure profiles for a Plexiglas model of the larynx with a glottis having a 10-degree divergence angle and either a symmetric orientation or an oblique angle of 15 degrees. For the oblique glottis, one side was divergent and the other convergent. The vocal fold surfaces had 14 pressure taps. The minimal glottal diameter was held constant at 0.04 cm. Results indicated that for either the symmetric or oblique case, the pressure profiles were different on the two sides of the glottis except for the symmetric geometry for a transglottal pressure of 3 cm H2O. For the symmetric case, flow separation created lower pressures on the side where the flow stayed attached to the wall, and the largest pressure differences between the two sides of the channel were 5%-6% of the transglottal pressure. For the oblique case, pressures were lower on the divergent glottal side near the glottal entry and exit, and the cross-channel pressures at the glottis entrance differed by 27% of the transglottal pressure. The empirical pressure distributions were supported by computational results. The observed aerodynamic asymmetries could be a factor contributing to normal jitter values and differences in vocal fold phasing.     

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.     

Characterizing glottal jet turbulence

Air pressure associated with airflow from the lungs drives the vocal folds into oscillation and allows the air to exit the glottis as a turbulent jet, even though laminar flow may enter the glottis from the trachea. The separation of the turbulence from the deterministic portion of the glottal jet was investigated in the excised canine larynx model. The present study is methodological in that the main goal was to examine three methods of obtaining reasonable representations of both the deterministic signal and the residual turbulence portion: (a) smoothing, (b) wavelet denoising, and (c) ensemble averaging. Ensemble averaging resulted in a deterministic signal that disregarded gross cyclic alterations while exaggerating the turbulence intensity. Wavelet denoising can perform an excellent analysis and synthesis of the glottal velocity, but was problematic in determining which levels of analysis to choose to represent both the deterministic and turbulence appropriately. Smoothing appeared to be the most appropriate for phonation velocities because it preserved gross cyclic variations important to perturbations and modulations, while extracting turbulence at what appears to be reasonable levels.     

The occurrence of the Coanda effect in pulsatile flow through static models of the human vocal folds

Pulsatile flow through a one-sided diffuser and static divergent vocal-fold models is investigated to ascertain the relevance of viscous-driven flow asymmetries in the larynx. The models were 7.5 times real size, and the flow was scaled to match Reynolds and Strouhal numbers, as well as the translaryngeal pressure drop. The Reynolds number varied from 0-2000, for flow oscillation frequencies corresponding to 100 and 150 Hz life-size. Of particular interest was the development of glottal flow skewing by attachment to the bounding walls, or Coanda effect, in a pulsatile flow field, and its impact on speech. The vocal folds form a divergent passage during phases of the phonation cycle when viscous effects such as flow separation are important. It was found that for divergence angles of less than 20 degrees, the attachment of the flow to the vocal-fold walls occurred when the acceleration of the forcing function was zero, and the flow had reached maximum velocity. For a divergence angle of 40 degrees, the fully separated central jet never attached to the vocal-fold walls. Inferences are made regarding the impact of the Coanda effect on the sound source contribution in speech.     

Shaping function models of the phonatory excitation signal

The phonatory excitation signal is the acoustic signal that is generated at the glottis by the vibrating vocal folds and pulsatile airflow. A shaping function model is a nonlinear memoryless input-output characteristic that transforms a simple harmonic into the desired output. The model can be fitted linearly to observed or simulated template cycles. The instantaneous values of the excitation cycle centroid, amplitude as well as length, and the cues for phonatory identity are set via distinct parameters. The synthetic phonatory excitation signal is zero on average, as well as identically zero when the glottal airflow rate is constant.     

Physiologic and acoustic differences between male and female voices

Comparison is drawn between male and female larynges on the basis of overall size, vocal fold membranous length, elastic properties of tissue, and prephonatory glottal shape. Two scale factors are proposed that are useful for explaining differences in fundamental frequency, sound power, mean airflow, and glottal efficiency. Fundamental frequency is scaled primarily according to the membranous length of the vocal folds (scale factor of 1.6), whereas mean airflow, sound power, glottal efficiency, and amplitude of vibration include another scale factor (1.2) that relates to overall larynx size. Some explanations are given for observed sex differences in glottographic waveforms. In particular, the simulated (computer-modeled) vocal fold contact area is used to infer male-female differences in the shape of the glottis. The female glottis appears to converge more linearly (from bottom to top) than the male glottis, primarily because of medial surface bulging of the male vocal folds.     

On the role of glottis-interior sources in the production of voiced sound

The voice source is dominated by aeroacoustic sources downstream of the glottis. In this paper an investigation is made of the contribution to voiced speech of secondary sources within the glottis. The acoustic waveform is ultimately determined by the volume velocity of air at the glottis, which is controlled by vocal fold vibration, pressure forcing from the lungs, and unsteady backreactions from the sound and from the supraglottal air jet. The theory of aerodynamic sound is applied to study the influence on the fine details of the acoustic waveform of "potential flow" added-mass-type glottal sources, glottis friction, and vorticity either in the glottis-wall boundary layer or in the portion of the free jet shear layer within the glottis. These sources govern predominantly the high frequency content of the sound when the glottis is near closure. A detailed analysis performed for a canonical, cylindrical glottis of rectangular cross section indicates that glottis-interior boundary/shear layer vortex sources and the surface frictional source are of comparable importance; the influence of the potential flow source is about an order of magnitude smaller.     

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.     

A computational study of asymmetric glottal jet deflection during phonation

Two-dimensional numerical simulations are used to explore the mechanism for asymmetric deflection of the glottal jet during phonation. The model employs the full Navier-Stokes equations for the flow but a simple laryngeal geometry and vocal-fold motion. The study focuses on the effect of Reynolds number and glottal opening angle with a particular emphasis on examining the importance of the so-called "Coanda effect" in jet deflection. The study indicates that the glottal opening angle has no substantial effect on glottal jet deflection. Deflection in the glottal jet is always preceded by large-scale asymmetry in the downstream portion of the glottal jet. A detailed analysis of the velocity and vorticity fields shows that these downstream asymmetric vortex structures induce a flow at the glottal exit which is the primary driver for glottal jet deflection.