An investigation of jet trajectory in flow through scaled vocal fold models with asymmetric glottal passages

Pulsatile two-dimensional flow through asymmetric static divergent models of the human vocal folds is investigated. Included glottal divergence angles are varied between 10° and 30°, with asymmetry angles between the vocal fold pairs ranging from 5° to 15°. The model glottal configurations represent asymmetries that arise during a phonatory cycle due to voice disorders. The flow is scaled to physiological values of Reynolds, Strouhal, and Euler numbers. Data are acquired in the anterior–posterior mid-plane of the vocal fold models using phase-averaged Particle Image Velocimetry (PIV) acquired at ten discrete locations in a phonatory cycle. Glottal jet stability arising from the vocal fold asymmetries is investigated and compared to previously reported work for symmetric vocal fold passages. Jet stability is enhanced with an increase in the included divergence angle, and the glottal asymmetry. Concurrently, the bi-modal jet trajectory and flow unsteadiness diminishes. Consistent with previous findings, the flow attachment due to the Coanda effect occurs when the acceleration of the forcing function is zero.     

Three-dimensional laryngeal flow fields induced by a model vocal fold polyp

Pathological laryngeal flow fields are investigated in a dynamically-driven, scaled-up model of the vocal folds. Disruption of the flow field due to the presence of a geometric protuberance, representative of a sessile unilateral polyp, is investigated in both the streamwise and transverse flow directions using phase-averaged particle image velocimetry. It is shown that the protuberance disrupts the normal flow behavior of the glottal jet throughout the phonatory cycle. During the divergent portions of the glottal cycle, the flow is characterized by the formation of hairpin vortices downstream of the protuberance. The protuberance also introduces significant velocity gradients in the anterior-posterior direction, which cover ∼30 − 40% of the vocal fold length. It is proposed that the disruption of the normal velocity behavior owing to the presence of a polyp will adversely impact the aerodynamic loadings that drive vocal fold motion, contributing to the temporal and spatial vocal fold asymmetries that are clinically-observed in patients with unilateral polyps.     

An investigation of asymmetric flow features in a scaled-up driven model of the human vocal folds

Flow through a driven, 7.5 times life-size vocal fold model was investigated at corresponding life-size flow rates of Q _mean  = 89.1 ml/s, 159.4 ml/s, and 253.0 ml/s. The flow was scaled to match physiological values for Reynolds, Strouhal, and Euler numbers. The models were driven at a life-size frequency of 94 Hz. Particle image velocimetry (PIV) data were acquired in the anterior–posterior midplane of the glottis, and the unsteady transglottal pressure drop across the vocal folds was simultaneously measured. Flow and pressure data were obtained at four discrete instances during the closing phases of the phonatory cycle for which t/T _open  = 0.60, 0.70, 0.80, and 0.90. The glottal jet trajectory exhibited a bimodal distribution of flow attachment between the two medial surfaces of the glottis. Vortex shedding at the trailing edge separation point generated instabilities in the shear layer, which caused large oscillations in the glottal jet orientation downstream of the glottal exit. The development of the Coanda effect during the glottal cycle was found to have minimal impact on the transglottal pressure drop, suggesting that flow orientation does not directly influence the dipole sound source. The change in transglottal pressure drop as a result of jet trajectory was less than 2% for all three investigated flow rates.     

A reduced-order flow model for vocal fold vibration: From idealized to subject-specific models

We present a reduced-order model for fluid–structure interaction (FSI) simulation of vocal fold vibration during phonation. This model couples the three-dimensional (3D) tissue mechanics and a one-dimensional (1D) flow model that is derived from the momentum and mass conservation equations for the glottal airflow. The effects of glottal entrance and pressure loss in the glottis are incorporated in the flow model. We consider both idealized vocal fold geometries and subject-specific anatomical geometries segmented from the MRI images of rabbits. For the idealized vocal fold geometries, we compare the simulation results from the 1D/3D hybrid FSI model with those from the full 3D FSI simulation based on an immersed-boundary method. For the subject-specific geometries, we incorporate previously estimated tissue properties for individual samples and compare the results with those from the high-speed imaging experiment of in vivo phonation. In both setups, the comparison shows good agreement in the vibration frequency, amplitude, phase delay, and deformation pattern of the vocal fold, which suggests potential application of the present approach for future patient-specific modeling.     

On the jet formation through a leaky glottis

The present experimental study aims at analyzing the jet formation of the glottal jet flow using a model of a leaky glottis. Experiments were performed in a flow channel with dynamic models of the vocal folds in order to measure the glottal waveform and the velocity distribution in the supraglottal region using High Speed Particle Image Velocimetry (HSPIV). Proper Orthogonal Decomposition (POD) of the vortex Q-criterion was performed in order to detect the energetically most significant large-scale vortex structures and their appearance in the jet flow. The spectral analysis of the glottal waveforms results in an increased spectral decay with more prominent peaks at higher frequencies in the case of a leaky glottis compared to a completely closing reference case. Vortex induced fluctuation frequencies have similar spectral content in both cases as they appear as trains of vortex packets in a regular manner over the glottal cycle. However, when removing the false vocal folds in the leaky glottis model, coherence of vortex generation is lost over the motion cycle. Thus, the presence of the false vocal folds retains most of the vortex induced characteristics in the source spectrum even when the vocal folds do not close fully.     

合成射流对失速状态下翼型大分离流动控制的试验研究

为研究低速状态合成射流在抑制翼型气流分离和推迟失速方面的控制机理, 开展了NACA0021 翼型失速特性射流控制的风洞试验研究. 通过系统性的模型测力、翼型瞬态流场粒子图像测速和边界层速度测定的对比试验, 深入探索了合成射流各参数对翼型失速特性控制效果的影响规律. 试验结果表明射流偏角在翼型升力和失速迎角控制方面的效果对射流动量系数较为敏感: 当动量系数较大时, 近切向射流的控制效果更好. 射流动量系数为0.033 时, 偏角30°的射流使得翼型升力系数峰值提高23.56%, 失速迎角增大5°; 而动量系数较小时, 偏角较大的射流能够获得最佳控制效果. 射流动量系数为0.0026 时, 法向射流对翼型最大升力系数控制效果最好(提高9.2%).      

High-speed PIV measurements of the flow downstream of a dynamic mechanical model of the human vocal folds

The objective of the present study is the detailed analysis of the unsteady vortex dynamics downstream of the human glottis during phonation at typical fundamental frequencies of the male voice of about 120 Hz. A hydraulic respiratory mock circuit has been built, including a factor of three up-scaled realistic dynamic model of the vocal folds. Time-resolving flow measurements were carried out downstream of the glottis by means of high-speed particle image velocimetry (PIV). The function of the human glottis is reproduced by two counter-rotating cams, each of which is covered with a stretched silicone membrane. The three-dimensional (3-D) geometry of the cams is designed such that the rotation leads to a realistic time-varying motion and profile of the glottis and waveform of the glottal cycle. Using high-speed PIV, the velocity field is captured with high spatial and temporal resolution to investigate the unsteady vortex dynamics of the cyclic jet-like flow in the vocal tract. The results help us to understand the vorticity interaction within the pulsating jet and, consequently, the generated sound in a human voice. In addition, changing the 3-D contours of the cams enables us to investigate basic pathological differences of the glottis function and the resulting alterations of the velocity and vorticity field in the vocal tract. The results are presented for typical physiological flow conditions in the human glottis. The frequencies of periodic vortex structures generated downstream of the glottis are fivefold higher than the fundamental frequency of the vocal folds oscillation. The highest vorticity fluctuations have a phase shift of 35% relative to the opening of the glottis. Finally, the flow field in the vocal tract is identified to be highly three-dimensional.     

A three-dimensional study of the glottal jet

This work builds upon the efforts to characterize the three-dimensional features of the glottal jet during vocal fold vibration. The study uses a Stereoscopic Particle Image Velocimetry setup on a self-oscillating physical model of the vocal folds with a uniform vocal tract. Time averages are documented and analyzed within the framework given by observations reported for jets exiting elongated nozzles. Phase averages are locked to the audio signal and used to obtain a volumetric reconstruction of the jet. From this reconstruction, the intra-cycle dynamics of the jet axis switching is disclosed.     

Measurement of flow separation in a human vocal folds model

The paper provides experimental data on flow separation from a model of the human vocal folds. Data were measured on a four times scaled physical model, where one vocal fold was fixed and the other oscillated due to fluid–structure interaction. The vocal folds were fabricated from silicone rubber and placed on elastic support in the wall of a transparent wind tunnel. A PIV system was used to visualize the flow fields immediately downstream of the glottis and to measure the velocity fields. From the visualizations, the position of the flow separation point was evaluated using a semiautomatic procedure and plotted for different airflow velocities. The separation point position was quantified relative to the orifice width separately for the left and right vocal folds to account for flow asymmetry. The results indicate that the flow separation point remains close to the narrowest cross-section during most of the vocal fold vibration cycle, but moves significantly further downstream shortly prior to and after glottal closure.     

超声速气流中雾化燃料喷射的三维数值研究

首次用双流体模型对雾化燃料在扩张形超燃室中沿九喷嘴顺流喷射的混合问题进行了数值研究。气相用迎风 TVD格式求解三维全 Navier- Stokes方程 ,液相用预估、校正 NND格式求解三维 Euler方程。相间相互作用的常微分方程用预估、校正Runge- Kutta法求解。用三维 Poisson方程生成网格。结果表明 :气相较液相的扩散效果好 ,小直径液滴的扩散效果好。相间速度滑移、改变气相喷射压力和喷射速度对液相扩散的贡献不大 ,但调整喷射角度会明显地增强液相的扩散、混合 ,本文结果未出现阻塞。     

Modified flapping jet for increased jet spreading using synthetic jets

The present paper is an experimental investigation, using a PIV system, on modified rectangular jet flow co-flowing with a pair of synthetic jets placed symmetrically with respect to the geometric centerline of the main flow. The objective was to determine the optimal forcing conditions that would result in jet spreading beyond what would be obtained in a simple flapped jet. The main jet had an exit Re_h = 36,000, based on the slot height, h. The synthetic jets were operated in a periodic manner with a periodic momentum coefficient of about 3.3% and at a frequency of the main jet preferred mode. A short, wide angle diffuser of half angle of about 45° was attached to the main jet. Generally for the vectored jet, much of the flow features found here resembled those reported in the literature except that the deflection angle in this study increased with downstream distances inside the diffuser and then remained roughly unchanged thereafter. Larger jet spreading was achieved when the main jet was subjected to simultaneous actuation of the synthetic jets but the flow did not achieve the initial jet spreading that was observed in the vectored jet. Further jet spreading was achieved when the synthetic jets were alternately actuated in which each synthetic jet was actuated for a number of cycles before switching. This technique allowed the jet to flap across the flow between transverse positions larger than what would be obtained in a simple flip-flop jet. Under the present flow geometry and Reynolds number, it was found that when the ratio f_s/f_al, where f_s is the synthetic actuation frequency and f_al is the alternating frequency, was larger than 10, the mean streamwise velocity of the main jet had two peaks symmetrically placed with respect to the jet axis and the jet had the appearance of flowing into two streams each moving nearly parallel to the diffuser wall. For a value of f_s/f_al of about 10, the optimal value in this study, the desired flow properties were achieved in that, the mean velocity was nearly uniform with an increase in the jet width compared to the simultaneous actuations, and the jet flapping was more effective in redistributing and homogenizing the turbulent kinetic energy across the main jet.     

Interaction of a transverse gas jet with supersonic flow in a dihedral

The results of an experimental investigation of supersonic flow at the Mach number M₁=3 past a transverse gas jet flowing from an orifice in the edge of a dihedral with a linear angle of 90° are presented. The Mach number of the jet was varied from 1 to 3, and the ratio of the total pressure in the jet to the free stream pressure from 90 to 760. Visualization of the flow near the faces of the dihedral revealed the existence of internal lines of flow convergence and divergence in the region of three-dimensional separated flow, which indicates the presence of complex vortex structures. The dependence of the dimensions of the separated flow zone and the characteristic pressures in it on the jet parameters is explored.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 6, pp. 161–165, November–December, 1994.     

The effect of density ratio on the near field of a naturally occurring oscillating jet

The triangular oscillating jet nozzle generates a triangular jet partially confined within an axi-symmetric chamber to produce a large scale flow oscillation that has application in thermal processes. Particle image velocimetry and oscillation frequency measurements were conducted to investigate the influence of the jet fluid to ambient fluid density ratio on the resulting oscillating flow. The investigation was conducted with a jet momentum flux of 0.06 kg m s⁻² (Re = 7.3−47.2 × 10³) and density ratios ranging from 0.2 to 5.0. The initial spread and decay of the emerging jet was found to depend upon the density ratio but in a more complex way than does an unconfined jet. Both the spread and decay are strongly influenced by the instantaneous angle of jet deflection, with greater deflection leading to increased spreading and decay of the jet. Decreasing the density ratio below unity results in a rapid decrease in the deflection angle, while increasing the density ratio above unity results in an increase in the deflection angle, albeit with less sensitivity. The frequency of oscillation was also shown to depend on the density ratio with an increase in the density ratio causing a decrease in the dominant oscillation frequency.     

The formation of vortex rings in a strongly forced round jet

The periodic formation of vortex rings in the developing region of a round jet subjected to high-amplitude acoustic forcing is investigated with High-Speed Particle Image Velocimetry. Harmonic velocity oscillations ranging from 20 to 120% of the mean exit velocity of the jet was achieved at several forcing frequencies determined by the acoustic response of the system. The time-resolved history of the formation process and circulation of the vortex rings are evaluated as a function of the forcing conditions. Overall, high-amplitude forcing causes the shear layers of the jet to breakup into a train of large-scale vortex rings, which share many of the features of starting jets. Features of the jet breakup such as the roll-up location and vortex size were found to be both amplitude and frequency dependent. A limiting time-scale of t/T ≈ 0.33 based on the normalized forcing period was found to restrict the growth of a vortex ring in terms of its circulation for any given arrangement of jet forcing conditions. In sinusoidally forced jets, this time-scale corresponds to a kinematic constraint where the translational velocity of the vortex ring exceeds the shear layer velocity that imposes pinch-off. This kinematic constraint results from the change in sign in the jet acceleration between t = 0 and t = 0.33T. However, some vortex rings were observed to pinch-off before t = 0.33T suggesting that they had acquired their maximum circulation. By invoking the slug model approximations and defining the slug parameters based on the experimentally obtained time- and length-scales, an analytical model based on the slug and ring energies revealed that the formation number for a sinusoidally forced jet is L/D ≈ 4 in agreement with the results of Gharib et al. (J Fluid Mech 360:121–140, 1998).     

基于主动流动控制的射流矢量偏转技术

本文介绍了一种基于主动流动控制技术的射流矢量偏转新方法和控制思路。通过在主射流出口两侧加装斜置扩张固壁板来降低射流两侧与固壁边界间的流体压力,将射流偏转由"不敏感-难控"转变成"敏感-易控",再在固壁板布置自行研制的斜出口合成射流激励器对主射流进行比例偏转控制。实验结果表明,射流最大偏转角可达15°。此外还研究了激励位置角度、激励频率、激励电压不同工作参数对射流矢量偏转控制的影响,实现了主射流偏转角的比例控制。当合成射流与主射流动量比为1∶43时,主射流偏转角可达13°,合成射流激励器消耗的能量为1.5W,初步实现了以小的能量消耗获取高的控制效益。     

Kinematic analysis of disturbed flow past a plate in a channel with porous walls

Analytic expressions for the complex flow potential are obtained in the linear formulation in the neighborhood of a plate at a small angle of incidence and near porous channel walls. The general solution includes the limiting cases of a plate in a channel with impermeable walls and in a jet. Numerical results concerning the effect of porosity on the flow geometry in the neighborhood of the plate and the channel walls are presented. The disturbed-flow velocity distributions along the channel walls and the flow rate of the fluid sinking at infinity are obtained.Cheboksary. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 13–19, July–August, 1995.     

Influence of subglottic stenosis on the flow-induced vibration of a computational vocal fold model

The effect of subglottic stenosis on vocal fold vibration is investigated. An idealized stenosis is defined, parameterized, and incorporated into a two-dimensional, fully coupled finite element model of the vocal folds and laryngeal airway. Flow-induced responses of the vocal fold model to varying severities of stenosis are compared. The model vibration was not appreciably affected by stenosis severities of up to 60% occlusion. Model vibration was altered by stenosis severities of 90% or greater, evidenced by decreased superior model displacement, glottal width amplitude, and flow rate amplitude. Predictions of vibration frequency and maximum flow declination rate were also altered by high stenosis severities. The observed changes became more pronounced with increasing stenosis severity and inlet pressure, and the trends correlated well with flow resistance calculations. Flow visualization was used to characterize subglottal flow patterns in the space between the stenosis and the vocal folds. Underlying mechanisms for the observed changes, possible implications for human voice production, and suggestions for future work are discussed.     

Characteristics of strongly-forced turbulent jets and non-premixed jet flames

Previous researchers have demonstrated that strong pulsations of the fuel flow rate can significantly reduce the flame length and luminosity of laminar/transitional non-premixed jet flames. The physical mechanisms responsible for these changes are investigated experimentally in acoustically-forced jet flows where the peak velocity fluctuations are up to eight times the mean flow velocity. Both reacting and non-reacting flows were studied and Reynolds numbers, based on the mean flow properties, ranged from 800 to 10,000 (corresponding to peak Reynolds numbers of 1,450–23,000), and forcing frequencies ranged from 290 to 1,140 Hz. Both the first and second organ-pipe resonance modes of the fuel delivery tube were excited to obtain these frequencies. An analysis of the acoustic forcing characteristics within the resonance tube is provided in order to understand the source of the high amplitude forcing. Flow visualization of jets with first resonant forcing confirms the presence of large-scale coherent vortices and strong reverse flow near the exit of the fuel tube. With second-resonant forcing, however, vortices are not emitted from the tube as they are drawn back into the fuel tube before they can fully form. Increased fine-scale turbulence is associated with both resonant cases, but particularly at second resonance. The power spectra of the velocity fluctuations for a resonantly pulsed jet show the presence of an inertial subrange indicating that the flow becomes fully turbulent even for mean-Reynolds-number jets that are nominally laminar. It is shown that these pulsed jet flows exhibit strong similarities to synthetic jets and that the Strouhal number, based on the maximum velocity at the fuel tube exit, is the dominant parameter for scaling these flows. The Strouhal number determines the downstream location where the coherent vortices breakdown, and is found to provide better collapse of flame length data (both current and previous) than other parameters that have been used in the literature.     

Mixed convection flow investigation in a rectangular horizontal tube stenosis via liquid crystal thermography and planar PIV

The temperature and velocity field in a horizontal convergent-divergent rectangular channel heated from below is studied experimentally for a Reynolds range 8-120, Grashof numbers from 0.44 × 10⁵ to 2.56 × 10⁵ and Richardson numbers from 3 to 4000, using water as working fluid. The duct aspect ratio (width/height) varies from 1 at its inlet to 2.28 at the stenosis neck, and both the upper and bottom walls are tilted with an angle of ±15.7° with respect to the horizontal. The temperature of the bottom wall is kept constant above that of the issuing fluid. The temperature field is recorded by liquid crystals in the vertical mid plane whereas the velocity field is measured in this plane as well as in four cross sections of the divergent passage by planar PIV, revealing the characteristics of the secondary velocity field. For all the examined cases the flow in the convergent passage is free of thermal plumes, and the thermal boundary layer is thin. In contrast, the divergent passage is characterized by a thermal plume which is shifted upstream with increasing Gr or reducing Re. Both transversal and longitudinal rolls emerge in this diffuser the strength of which depend on Re and Gr influencing the streamwise distribution of Nusselt which for low Re presents a minimum.     

Rotation of the thrust vector of a two-dimensional nozzle by means of displacing the critical surface

The method of thrust vector deflection by means of displacing the critical surface of a nozzle is numerically and experimentally investigated. The displacement is realized at the expense of extending rotatable walls, or deflectors, into the flow; in this case, one of the edges of the critical surface is displaced from the throat onto the deflector surface. Two nozzle configurations, with short and long deflectors, are studied. The thrust vector deflection angles and the nozzle thrust coefficients are determined in the thrust vector deflection regime. For the configuration with long deflectors the angle of rotation of the thrust vector amounts to 30° and is determined by the effect of jet ejection toward the opposite wall, similar with the Coanda effect.