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.     

Direct-numerical simulation of the glottal jet and vocal-fold dynamics in a three-dimensional laryngeal model

An immersed-boundary method based flow solver coupled with a finite-element solid dynamics solver is employed in order to conduct direct-numerical simulations of phonatory dynamics in a three-dimensional model of the human larynx. The computed features of the glottal flow including mean and peak flow rates, and the open and skewness quotients are found to be within the normal physiological range. The flow-induced vibration pattern shows the classical "convergent-divergent" glottal shape, and the vibration amplitude is also found to be typical for human phonation. The vocal fold motion is analyzed through the method of empirical eigenfunctions and this analysis indicates a 1:1 modal entrainment between the "adduction-abduction" mode and the "mucosal wave" mode. The glottal jet is found to exhibit noticeable cycle-to-cycle asymmetric deflections and the mechanism underlying this phenomenon is examined.     

Flow separation in a computational oscillating vocal fold model

A finite-volume computational model that solves the time-dependent glottal airflow within a forced-oscillation model of the glottis was employed to study glottal flow separation. Tracheal input velocity was independently controlled with a sinusoidally varying parabolic velocity profile. Control parameters included flow rate (Reynolds number), oscillation frequency and amplitude of the vocal folds, and the phase difference between the superior and inferior glottal margins. Results for static divergent glottal shapes suggest that velocity increase caused glottal separation to move downstream, but reduction in velocity increase and velocity decrease moved the separation upstream. At the fixed frequency, an increase of amplitude of the glottal walls moved the separation further downstream during glottal closing. Increase of Reynolds number caused the flow separation to move upstream in the glottis. The flow separation cross-sectional ratio ranged from approximately 1.1 to 1.9 (average of 1.47) for the divergent shapes. Results suggest that there may be a strong interaction of rate of change of airflow, inertia, and wall movement. Flow separation appeared to be "delayed" during the vibratory cycle, leading to movement of the separation point upstream of the glottal end only after a significant divergent angle was reached, and to persist upstream into the convergent phase of the cycle.     

A theoretical model of the pressure field arising from asymmetric intraglottal flows applied to a two-mass model of the vocal folds

A theoretical flow solution is presented for predicting the pressure distribution along the vocal fold walls arising from asymmetric flow that forms during the closing phases of speech. The resultant wall jet was analyzed using boundary layer methods in a non-inertial reference frame attached to the moving wall. A solution for the near-wall velocity profiles on the flow wall was developed based on a Falkner-Skan similarity solution and it was demonstrated that the pressure distribution along the flow wall is imposed by the velocity in the inviscid core of the wall jet. The method was validated with experimental velocity data from 7.5 times life-size vocal fold models, acquired for varying flow rates and glottal divergence angles. The solution for the asymmetric pressures was incorporated into a widely used two-mass model of vocal fold oscillation with a coupled acoustical model of sound propagation. Asymmetric pressure loading was found to facilitate glottal closure, which yielded only slightly higher values of maximum flow declination rate and radiated sound, and a small decrease in the slope of the spectral tilt. While the impact on symmetrically tensioned vocal folds was small, results indicate the effect becomes more significant for asymmetrically tensioned vocal folds.     

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.     

Impact of wall rotation on supraglottal jet stability in voiced speech

Supraglottal jet variability was investigated in a scaled-up flow facility incorporating driven vocal fold models with and without wall rotation. Principle component analysis was performed on the experimental supraglottal flow fields to ascertain the roll of glottal wall motion on the development of the supraglottal jet. It is shown that intraglottal flow asymmetries that develop due to wall rotation are not the primary mechanism for generating large-scale cycle-to-cycle deflection of the supraglottal jet. However, wall rotation does decrease the energy content of the first mode, redistributing it to the higher modes through an increase in unstructured flow variability.     

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.     

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.     

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.     

Asymmetric glottal jet deflection: differences of two- and three-dimensional models

Flow is studied through a channel with an oscillating orifice mimicking the motion of the glottal-gap during phonation. Simulations with prescribed flow and wall-motion are carried out for different orifice geometries, a 2D slit-like and a 3D lens-like one. Although the jet emerges from a symmetric orifice a significant deflection occurs in case of the slit-like geometry, contrary to the 3D lens-like one. The results demonstrate the dependency of jet entrainment and vortex dynamics on the orifice geometry and the interpretation of asymmetric jet deflection with regard to the relevance of the Coanda effect in the process of human phonation.     

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.     

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.     

Velocity distributions in glottal models

Velocity distributions within three models of the human larynx, namely, a rigid plexiglas model, an excised canine larynx, and a computational model are investigated with experimental and theoretical analyses. A plexiglas wind tunnel with interchangeable glottal constrictions was used as a two-dimensional steady-flow model to measure velocity and pressure for various glottal shapes. A canine excised larynx was used as a prototype pulsatile flow model to study pressure and velocity variations during phonation. Results of the plexiglas modelling indicated a parabolic laminar velocity profile upstream of the glottal constriction and turbulent and asymmetric velocity profiles downstream of the glottal constriction. The time-averaged velocities of the excised larynx had similarities with the plexiglas model results, and instabilities and asymmetries were also demonstrated by the computational method.     

Snapshot POD analysis of transient flow in the pilot stage of a jet pipe servo valve

The spatial evolution of a turbulent flow in the pilot stage of a jet pipe servo valve at the inlet pressure and deflection angle of the jet pipe is investigated using a large eddy simulation (LES). The pressure of the same flow field is measured by a high frequency dynamic pressure sensor in the experiments and is compared with the results of the LES, as well as their root-mean-square (RMS) and fast Fourier transform (FFT) results. The results of experiments and LES are in good agreement, indicating that LES is able to predict the flow dynamics. Velocity datasets based on LES are utilised to conduct the snapshot proper orthogonal decomposition (snapshot POD) technique. The snapshot POD analysis results of the first 4 modes show a full ability to directly visualise details of the coherent structures. The influences of the inlet pressure and deflection angle of the jet pipe are also discussed. Under different inlet pressures, the velocity eigenfunctions of the first mode are similar, while the locations and strengths of the vortices in high modes are different. The Lamb-Oseen vortices that affect the trajectory of jet streams are observed in the vicinity of the entrances of receiver channels only in the first mode, and several spindly vortices appear in the region of ?5?y/n?相似文献   

A numerical study of rotation effects on jets effusing from inclined holes into a cross-flow

In this work, the complexity of the flow field arising from the impact of the interaction of coolant jets with a hot cross-flow under rotation conditions was numerically simulated using large eddy simulation with artificial inflow boundary condition. The finite-volume method and the unsteady PISO (Pressure Implicit with Splitting of Operators) algorithm were applied on a non-uniform staggered grid. The simulations were performed for four different values of rotation number (Ro) of 0.0, 0.03021, 0.06042, and 0.12084, a jet Reynolds number of 4700, based on the hole width and the jet exit velocity. The air jet was injected at 30° and 90° in the streamwise direction with a density ratio of 1.04 and a velocity ratio of 0.5. The flow fields of the present study were compared with experimental data in order to validate the reliability of the LES technique. It was shown that the rotation has a strong impact on the jet trajectory behaviour and the film cooling effectiveness. The film trajectory always inclines centrifugally. Under rotating conditions, the film trajectory departs from the centreline to the left boundary. The deflection becomes greater as Ro increases. Furthermore, it was also found that the injection angle has a strong impact on separation and reattachment behaviour as well as the strength of the penetration into the cross-flow. As it increases, the distribution of the film cooling downstream the jet exit is more non-uniform and the film cooling effectiveness level slightly decreases.     

介质阻挡体放电对甲烷扩散火焰特性的影响

本文在高频交流激励模式下,采用同轴圆柱构型激励器,开展了介质阻挡体放电对空气/甲烷同轴剪切扩散火焰燃烧特性影响实验研究。激励器敷设在外喷嘴环缝以电离空气,采用纹影系统和B型热电偶分别获取流场形态和火焰温度,激励频率为8 kHz,通过改变气体流量和放电电压,分析了不同工况下射流流场、火焰结构和火焰温度在等离子体作用下的变化规律。结果表明:等离子体气动效应能有效增强射流湍流强度,强化空气/甲烷掺混,增大射流角,并随激励电压提高作用效果逐渐增强,实验中未形成明显扩张流动的初始射流在放电电压30 kV时其射流角最大为23.5°。贫燃条件下等离子体激励会改善火焰形态,增强燃烧稳定性,并在流量较低时缩短火焰长度。此外,富燃火焰下游温度会随着激励强度增大不断升高,而贫燃火焰下游温度变化受上游燃烧强度影响存在升高和降低两种情况。     

Caractérisation expérimentale et numérique des transferts de chaleur à travers un jet d'air plan horizontal

This study reports on heat transfer through a plan horizontal air jet. The jet maintained at a higher temperature than the ambience allows a thermal sealing at the upper side of the opening case.The thermal and hydrodynamic measurements taken on the experimental device and the numerical simulations carried out by using the Phoenics package distinguish two types of flow. For the small velocities of the jet or the high differences of temperature, the natural convection predominates: the hot jet diffuses upwards and can not seal the case. Inversely, for the high values of the jet velocities, the forced convection predominates: the case is correctly sealed but the power consumption increases like the flowrate. A dimensionless number called deflection modulus allows to characterize the transition from one flow pattern to another and it also allows to minimize the energy consumption.     

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.     

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.