Stefan Kniesburges

Flow-structure-acoustic interaction in a human voice model

For the investigation of the physical processes of human phonation, inhomogeneous synthetic vocal folds were developed to represent the full fluid-structure-acoustic coupling. They consisted of polyurethane rubber with a stiffness in the range of human vocal folds and were mounted in a channel, shaped like the vocal tract in the supraglottal region. This test facility permitted extensive observations of flow-induced vocal fold vibrations, the periodic flow field, and the acoustic signals in the far field of the channel. Detailed measurements were performed applying particle-image velocimetry, a laser-scanning vibrometer, a microphone, unsteady pressure sensors, and a hot-wire probe, with the aim of identifying the physical mechanisms in human phonation. The results support the existence of the Coanda effect during phonation, with the flow attaching to one vocal fold and separating from the other. This behavior is not linked to one vocal fold and changes stochastically from cycle to cycle. The oscillating flow field generates a tonal sound. The broadband noise is presumed to be caused by the interaction of the asymmetric flow with the downstream-facing surfaces of the vocal folds, analogous to trailing-edge noise.

In vitro experimental investigation of voice production.

The process of human phonation involves a complex interaction between the physical domains of structural dynamics, fluid flow, and acoustic sound production and radiation. Given the high degree of nonlinearity of these processes, even small anatomical or physiological disturbances can significantly affect the voice signal. In the worst cases, patients can lose their voice and hence the normal mode of speech communication. To improve medical therapies and surgical techniques it is very important to understand better the physics of the human phonation process. Due to the limited experimental access to the human larynx, alternative strategies, including artificial vocal folds, have been developed. The following review gives an overview of experimental investigations of artificial vocal folds within the last 30 years. The models are sorted into three groups: static models, externally driven models, and self-oscillating models. The focus is on the different models of the human vocal folds and on the ways in which they have been applied.

3D Reconstruction of Human Laryngeal Dynamics Based on Endoscopic High-Speed Recordings

Standard laryngoscopic imaging techniques provide only limited two-dimensional insights into the vocal fold vibrations not taking the vertical component into account. However, previous experiments have shown a significant vertical component in the vibration of the vocal folds. We present a 3D reconstruction of the entire superior vocal fold surface from 2D high-speed videoendoscopy via stereo triangulation. In a typical camera-laser set-up the structured laser light pattern is projected on the vocal folds and captured at 4000 fps. The measuring device is suitable for in vivo application since the external dimensions of the miniaturized set-up barely exceed the size of a standard rigid laryngoscope. We provide a conservative estimate on the resulting resolution based on the hardware components and point out the possibilities and limitations of the miniaturized camera-laser set-up. In addition to the 3D vocal fold surface, we extended previous approaches with a G2-continuous model of the vocal fold edge. The clinical applicability was successfully established by the reconstruction of visual data acquired from 2D in vivo high-speed recordings of a female and a male subject. We present extracted dynamic parameters like maximum amplitude and velocity in the vertical direction. The additional vertical component reveals deeper insights into the vibratory dynamics of the vocal folds by means of a non-invasive method. The successful miniaturization allows for in vivo application giving access to the most realistic model available and hence enables a comprehensive understanding of the human phonation process.

Influence of Vortical Flow Structures on the Glottal Jet Location in the Supraglottal Region

Within the fully coupled multiphysics phonation process, the fluid flow plays an important role for sound production. This study addresses phenomena in the flow downstream of synthetic self-oscillating vocal folds. An experimental setup consisting of devices for producing and conditioning the flow including the main test channel was applied. The supraglottal channel was designed to prevent an acoustic coupling to the vocal folds. Hence, the oscillations were aerodynamically driven. The cross-section of the supraglottal channel was systematically varied by increasing the distance between the lateral channel walls. The vocal folds consisted of silicone rubber of homogenous material distribution generating self-sustained oscillations. The airflow was visualized in the immediate supraglottal region using a laser-sheet technique and a digital high-speed camera. Furthermore, the flow was studied by measuring the static pressure distributions on both lateral supraglottal channel walls. The results clearly showed different flow characteristics depending on the supraglottal configuration. In all cases with supraglottal channel, the jet was located asymmetrical and bent in medial-lateral direction. Furthermore, the side to which the jet was deflected changed in between the consecutive cycles showing a bifurcational behavior. Previously, this phenomenon was explained by the Coanda effect. However, the present data suggest that the deflection of the jet was mainly caused by large air vortices in the supraglottal channel produced by the flow field of previous oscillations. In contrast, for the case without supraglottal channel, the air jet was found totally symmetrical stabilized by the constant pressure in the ambient region. The emitted sound signal showed additional subharmonic tonal peaks for the asymmetric flow cases, which are characteristics for diplophonia.

Dynamic vocal fold parameters with changing adduction in ex-vivo hemilarynx experiments

Ex-vivo hemilarynx experiments allow the visualization and quantification of three-dimensional dynamics of the medial vocal fold surface. For three excised human male larynges, the vibrational output, the glottal flow resistance, and the sound pressure during sustained phonation were analyzed as a function of vocal fold adduction for varying subglottal pressure. Empirical eigenfunctions, displacements, and velocities were investigated along the vocal fold surface. For two larynges, an increase of adduction level resulted in an increase of the glottal flow resistance at equal subglottal pressures. This caused an increase of lateral and vertical oscillation amplitudes and velocity indicating an improved energy transfer from the airflow to the vocal folds. In contrast, the third larynx exhibited an amplitude decrease for rising adduction accompanying reduction of the flow resistance. By evaluating the empirical eigenfunctions, this reduced flow resistance was assigned to an unbalanced oscillation pattern with predominantly lateral amplitudes. The results suggest that adduction facilitates the phonatory process by increasing the glottal flow resistance and enhancing the vibrational amplitudes. However, this interrelation only holds for a maintained balanced ratio between vertical and lateral displacements. Indeed, a balanced vertical-lateral oscillation pattern may be more beneficial to phonation than strong periodicity with predominantly lateral vibrations.

Evaluation of Analytical Modeling Functions for the Phonation Onset Process

The human voice originates from oscillations of the vocal folds in the larynx. The duration of the voice onset (VO), called the voice onset time (VOT), is currently under investigation as a clinical indicator for correct laryngeal functionality. Different analytical approaches for computing the VOT based on endoscopic imaging were compared to determine the most reliable method to quantify automatically the transient vocal fold oscillations during VO. Transnasal endoscopic imaging in combination with a high-speed camera (8000 fps) was applied to visualize the phonation onset process. Two different definitions of VO interval were investigated. Six analytical functions were tested that approximate the envelope of the filtered or unfiltered glottal area waveform (GAW) during phonation onset. A total of 126 recordings from nine healthy males and 210 recordings from 15 healthy females were evaluated. Three criteria were analyzed to determine the most appropriate computation approach: (1) reliability of the fit function for a correct approximation of VO; (2) consistency represented by the standard deviation of VOT; and (3) accuracy of the approximation of VO. The results suggest the computation of VOT by a fourth-order polynomial approximation in the interval between 32.2 and 67.8% of the saturation amplitude of the filtered GAW.

Numerical Simulation of Fluid-Structure- and Fluid-Structure-Acoustic Interaction Based on a Partitioned Coupling Scheme

A coupling scheme for the simulation of problems involving fluid-structure and fluid-structure-acoustic interactions is described. For the simulation of fluid-structure interactions an implicit coupling between the flow field and the mechanical structure is applied. In order to determine the acoustic radiation due to mechanical vibrations, a direct coupling between the mechanical structure and the acoustic field is employed. At the structure-acoustic interface non-matching grids are realized with mortar elements. The exchange of physical quantities between the structured mesh of the flow computation and the unstructured mesh of the structural computation is performed by MpCCI. Applications to the flow around an oscillating rigid body and the flow over a flexible plate are presented.

Automated setup for ex vivo larynx experiments

Ex vivo larynx experiments are limited in time due to degeneration of the laryngeal tissues. In order to acquire a significant and comparable amount of data, automatization of current manual experimental procedures is desirable. A computer controlled, electro-mechanical setup was developed for time-dependent variation of specific physiological parameters, including adduction and elongation level of the vocal folds and glottal flow. The setup offers a standardized method to induce defined forces on the laryngeal cartilages. Furthermore, phonation onset is detected automatically and the subsequent measurement procedure is automated and standardized to improve the efficiency of the experimental process. The setup was validated using four ex vivo porcine larynges, whereas each validation measurement series was executed with one separate larynx. Altogether 31 single measurements were undertaken, which can be summed up to a total experimental time of about 4 min. Vocal fold elongation and adduction lead both to an increase in fundamental frequency and subglottal pressure. Measurement procedures like applying defined subglottal pressure steps and onset-offset detection were reliably executed. The setup allows for a computer-based parameter control, which enables fast experimental execution over a wide range of laryngeal configurations. This maximizes the number of measurements and reduces personal effort compared with manual procedures.

Effect of the ventricular folds in a synthetic larynx model

Within the human larynx, the ventricular folds serve primarily as a protecting valve during swallowing. They are located directly above the sound-generating vocal folds. During normal phonation, the ventricular folds are passive structures that are not excited to periodical oscillations. However, the impact of the ventricular folds on the phonation process has not yet been finally clarified. An experimental synthetic human larynx model was used to investigate the effect of the ventricular folds on the phonation process. The model includes self-oscillating vocal fold models and allows the comparison of the pressure distribution at multiple locations in the larynx for configurations with and without ventricular folds. The results indicate that the ventricular folds increase the efficiency of the phonation process by reducing the phonation threshold level of the pressure below the vocal folds. Two effects caused by the ventricular folds could be identified as reasons: (1) a decrease in the mean pressure level in the region between vocal and ventricular folds (ventricles) and (2) an increase in the glottal flow resistance. The reason for the first effect is a reduction of the pressure level in the ventricles due to the jet entrainment and the low static pressure in the glottal jet. The second effect results from an increase in the glottal flow resistance that enhances the aerodynamic energy transfer into the vocal folds. This effect reduces the onset threshold of the pressure difference across the glottis.

The mechanisms of subharmonic tone generation in a synthetic larynx model

The sound spectra obtained in a synthetic larynx exhibited subharmonic tones that are characteristic for diplophonia. Although the generation of subharmonics is commonly associated with asymmetrically oscillating vocal folds, the synthetic elastic vocal folds showed symmetrical oscillations. The amplitudes of the subharmonics decreased with an increasing lateral diameter of the supraglottal channel, which indicates a strong dependence of the supraglottal boundary conditions. Investigations of the supraglottal flow field revealed small cycle-to-cycle variations of the static pressure in the region of the pulsatile glottal jet as the origin of the first subharmonic tone. It is located at half the fundamental frequency of the vocal fold oscillation. A principle component analysis of the supraglottal flow field with the fully developed glottal jet revealed a large recirculation area in the second spatial eigenvector which deflected the glottal jet slightly in a perpendicular direction of the jet axis. The rotation direction of the recirculation area changed with different oscillation cycles between clockwise and counterclockwise. As both directions were uniformly distributed across all acquired oscillation cycles, a cycle-wise change can be assumed. It is concluded that acoustic subharmonics are generated by small fluctuations of the glottal jet location favored by small lateral diameters of the supraglottal channel.