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.     

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.     

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 closure, transglottal airflow, and voice quality in healthy middle-aged women

Seventeen healthy women, 45 to 61 years old, were examined using videofiberstroboscopy during phonation at three loudness levels. Two phoniatricians evaluated glottal closure using category and ratio scales. Transglottal airflow was studied by inverse filtering of the oral airflow signal recorded in a flow mask (Glottal Enterprises System) during the spoken phrase /ba:pa:pa:pa:p/ at three loudness levels. Subglottal pressure was estimated from the intraoral pressure during p occlusion. Running speech and the repeated /pa:/ syllables were perceptually evaluated by three speech pathologists regarding breathiness, hypo-, and hyperfunction, using continuous scales. Incomplete glottal closure was found in 35 of 46 phonations (76%). The degree of glottal closure increased significantly with raised loudness. Half of the women closed the glottis completely during loud phonation. Posterior glottal chink (PGC) was the most common gap configuration and was found in 28 of 46 phonations (61%). One third of the PGCs were in the cartilaginous glottis (PGCc) only. Two thirds extended into the membranous portion (PGCm); most of these occurred during soft phonation. Peak flow, peak-to-peak (AC) flow, and the maximum rate of change for the flow in the closing phase increased significantly with raised loudness. Minimum flow decreased significantly from normal to loud voice. Breathiness decreased with increased loudness. The results suggest that the incomplete closure patterns PGCc and PGCm during soft phonation ought primarily to be regarded as normal for Swedish women in this age group.     

Comparing turbulence models for flow through a rigid glottal model

Flow through a rigid model of the human vocal tract featuring a divergent glottis was numerically modeled using the Reynolds-averaged Navier-Stokes approach. A number of different turbulence models, available in a widely used commercial computational fluid dynamics code, were tested to determine their ability to capture various flow features recently observed in laboratory experiments and large eddy simulation studies. The study reveals that results from unsteady simulations employing the k-omega shear stress transport model were in much better agreement with previous measurements and predictions with regard to the ability to predict glottal jet skewing due to the Coanda effect and the intraglottal pressure distribution or related skin friction coefficient, than either steady or unsteady simulations using the Spalart-Allmaras model or any other two-equation turbulence model investigated in this study.     

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.     

Numerical simulation of turbulence transition and sound radiation for flow through a rigid glottal model

Large eddy simulation (LES)-based computational aeroacoustics techniques were applied to a static model of the human glottis, idealized here as a planar channel with an orifice, to study flow-acoustic interactions related to speech. Rigid models of both converging and diverging glottal passages, each featuring a 20 deg included angle and a minimal glottal diameter of 0.04 cm, with an imposed transglottal pressure of 15 cm H2O, were studied. The Favre-filtered compressible Navier-Stokes equations were integrated for this low-Mach-number flow using an additive semi-implicit Runge-Kutta method and a high-order compact finite-difference scheme with characteristic-based nonreflecting boundary conditions and a multiblock approach. Flow asymmetries related to the Coanda effect and transition to turbulence, as well as the far-field sound, were captured. Acoustic-analogy-based far-field sound predictions were compared with direct simulations and showed that dipole sources, arising from unsteady flow forces exerted on the glottal walls, are primarily responsible for the tonal sound observed in the divergent glottis case.     

Investigation of a glottal related harmonics-to-noise ratio and spectral tilt as indicators of glottal noise in synthesized and human voice signals

The harmonics-to-noise ratio (HNR) of the voiced speech signal has implicitly been used to infer information regarding the turbulent noise level at the glottis. However, two problems exist for inferring glottal noise attributes from the HNR of the speech wave form: (i) the measure is fundamental frequency (f0) dependent for equal levels of glottal noise, and (ii) any deviation from signal periodicity affects the ratio, not just turbulent noise. An alternative harmonics-to-noise ratio formulation [glottal related HNR (GHNR')] is proposed to overcome the former problem. In GHNR' a mean over the spectral range of interest of the HNRs at specific harmonic/between-harmonic frequencies (expressed in linear scale) is calculated. For the latter issue [(ii)] two spectral tilt measures are shown, using synthesis data, to be sensitive to glottal noise while at the same time being comparatively insensitive to other glottal aperiodicities. The theoretical development predicts that the spectral tilt measures reduce as noise levels increase. A conventional HNR estimator, GHNR' and two spectral tilt measures are applied to a data set of 13 pathological and 12 normal voice samples. One of the tilt measures and GHNR' are shown to provide statistically significant differentiating power over a conventional HNR estimator.     

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.     

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.     

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.     

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.     

A preliminary study of particle velocity during phonation in an in vivo canine model

The particle velocity across the glottis was measured with simultaneous electroglottography, photoglottography, and subglottic pressure in an in vivo canine model of phonation. A constant temperature anemometer measured flow velocity at five midline anterior to posterior glottal positions. Tracheal input air flow was varied in five steps from 175 to 500 cc/s, while vocal fold approximation was achieved by constant electrical stimulation of the laryngeal nerves. For all levels of air flow, a decreasing peak velocity gradient was observed from the anterior commissure to the arytenoids. Time-varying features of the flow velocity are discussed in relation to glottal vibratory events and aerodynamics.     

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.     

腭裂语音中齿龈塞音的声门代偿现象声学分析与判定

对齿龈塞音在腭裂语音中的声门塞音代偿现象进行了声学分析,计算频谱分布的多阶统计量—谱矩,并将代偿塞音和正常塞音进行对比。结果显示声门塞音爆破段的第一阶谱矩即频谱质心的频率位置比正常塞音低,因为声门塞音的阻塞部位在声门,导致声道腔体偏长从而共振频率偏低。还观察到声门塞音的第二阶谱矩即标准偏差偏高,说明其谱能量分布比正常塞音更加分散。声门塞音的第三阶谱矩即偏度大多为正值,反映了声门塞音功率谱的非对称性且大头朝向低频区而长尾朝向高频区。采用逻辑回归模型进行样本分类,通过交叉验证选出最优的四阶谱矩作为模型自变量,分类正确率为89.7%。结合塞音爆破时刻自动检测,实现了音节/di/的声门塞音客观判定。      

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.     

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.     

Laryngeal vibrations: a comparison between high-speed filming and glottographic techniques

This study was designed to compare information on laryngeal vibrations obtained by high-speed filming, photoglottography (PGG), and electroglottography (ECG). Simultaneous glottographic signals and high-speed films were obtained from two subjects producing steady phonation. Measurements of glottal width were made at three points along the glottis in the anterior--posterior dimension and aligned with the other records. Results indicate that PGG and film measurements give essentially the same information for peak glottal opening and glottal closure. The EGG signal appears to reliably indicate vocal-fold contact. Together, PGG and EGG may provide much of the information obtained from high-speed filming as well as potentially detect horizontal phase differences during opening and closing.     

Broadband sound generation by confined pulsating jets in a mechanical model of the human larynx

Experiments were performed to study the production of broadband sound in confined pulsating jets through orifices with a time-varying area. The goal was to better understand broadband sound generation at the human glottis during voicing. The broadband component was extracted from measured sound signals by the elimination of the periodic component through ensemble averaging. Comparisons were made between the probability density functions of the broadband sound in pulsating jets and of comparable stationary jets. The results indicate that the quasi-steady approximation may be valid for the broadband component when the turbulence is well established and the turbulence kinetic energy is comparatively large. A wavelet analysis of the broadband sound showed that random sound production was modulated at the driving frequency. Two distinct sound production peaks were observed during one cycle, presumably associated firstly with jet formation and secondly with flow deceleration during orifice closing. Most high-frequency sound was produced during the closing phase. Deviations from quasi-steady behavior were observed. As the driving frequency increased, sound production during the opening phase was reduced, possibly due to the shorter time available for turbulence to develop. These results may be useful for better quality voice synthesis.     

Laryngeal Mechanisms During Human 4-kHz Vocalization Studied With CT, Videostroboscopy, and Color Doppler Imaging

The objective of this study was to investigate the underlying laryngeal mechanisms during the specific human 4-kHz vocalization. The laryngeal configuration during this vocalization was measured using high-resolution computerized tomographic scan and videostrobolaryngoscopy. The color Doppler imaging (CDI) of medical ultrasound was used to detect the vibrations of glottal and supraglottal mucosa. During the 4-kHz vocalization, the ventricular folds were adducted in the shape of a bimodal chink and the vocal folds were shaped as a "V" with an opening at the posterior glottis. In the coronal view, the laryngeal ventricles had collapsed and a divergent shaped conduit was observed at the posterior portion of the larynx. The surface mucosa vibration detected by CDI was noted over the bilateral ventricular folds and aryepiglottic folds. The vibration displacement was estimated to be on the order of 0.1mm. This vibration amplitude was too small to be detected in videostrobolaryngoscopy. The laryngeal configuration and CDI data suggested a diffuser jet with periodic vorticity bursts in the larynx producing 4 kHz voice.