The mechanisms of subharmonic tone generation in a synthetic larynx model

Synthetic, self-oscillating vocal fold models for voice production researcha)

Sound for the human voice is produced by vocal fold flow-induced vibration and involves a complex coupling between flow dynamics, tissue motion, and acoustics. Over the past three decades, synthetic, self-oscillating vocal fold models have played an increasingly important role in the study of these complex physical interactions. In particular, two types of models have been established: "membranous" vocal fold models, such as a water-filled latex tube, and "elastic solid" models, such as ultrasoft silicone formed into a vocal fold-like shape and in some cases with multiple layers of differing stiffness to mimic the human vocal fold tissue structure. In this review, the designs, capabilities, and limitations of these two types of models are presented. Considerations unique to the implementation of elastic solid models, including fabrication processes and materials, are discussed. Applications in which these models have been used to study the underlying mechanical principles that govern phonation are surveyed, and experimental techniques and configurations are reviewed. Finally, recommendations for continued development of these models for even more lifelike response and clinical relevance are summarized.

Restraining vocal fold vertical motion reduces source-filter interaction in a two-mass model

Previous experimental studies suggested that restraining the vocal fold vertical motion may reduce the coupling strength between the voice source and vocal tract. In this study, the effects of vocal fold vertical motion on source-filter interaction were systematically examined in a two-dimensional two-mass model coupled to a compressible flow simulation. The results showed that when allowed to move vertically, the vocal folds exhibited subharmonic vibration due to entrainment to the first vocal tract acoustic resonance. Restraining the vertical motion suppressed this entrainment. This indicates that the vertical mobility of the vocal folds may play a role in regulating source-filter interaction.

An Investigation of Acoustic Back-Coupling in Human Phonation on a Synthetic Larynx Model

In the human phonation process, acoustic standing waves in the vocal tract can influence the fluid flow through the glottis as well as vocal fold oscillation. To investigate the amount of acoustic back-coupling, the supraglottal flow field has been recorded via high-speed particle image velocimetry (PIV) in a synthetic larynx model for several configurations with different vocal tract lengths. Based on the obtained velocity fields, acoustic source terms were computed. Additionally, the sound radiation into the far field was recorded via microphone measurements and the vocal fold oscillation via high-speed camera recordings. The PIV measurements revealed that near a vocal tract resonance frequency fR, the vocal fold oscillation frequency fo (and therefore also the flow field's fundamental frequency) jumps onto fR. This is accompanied by a substantial relative increase in aeroacoustic sound generation efficiency. Furthermore, the measurements show that fo-fR-coupling increases vocal efficiency, signal-to-noise ratio, harmonics-to-noise ratio and cepstral peak prominence. At the same time, the glottal volume flow needed for stable vocal fold oscillation decreases strongly. All of this results in an improved voice quality and phonation efficiency so that a person phonating with fo-fR-coupling can phonate longer and with better voice quality.

Influence of Reduced Saliva Production on Phonation in Patients With Ectodermal Dysplasia

OBJECTIVE

Patients with ectodermal dysplasia (ED) suffer from an inherited disorder in the development of the ectodermal structures. Besides the main symptoms, i.e. significantly reduced formation/expression of teeth, hair and sweat glands, a decreased saliva production is objectively accounted. In addition to difficulties with chewing/swallowing, ED patients frequently report on the subjective impression of rough and hoarse voices. A correlation between the reduced production of saliva and an affliction of the voice has not yet been investigated objectively for this rare disease.

METHODS

Following an established measurement protocol, a study has been conducted on 31 patients with ED and 47 controls (no ED, healthy voice). Additionally, the vocal fold oscillations were recorded by high-speed videoendoscopy (HSV@4 kHz). The glottal area waveform was determined by segmentation and objective glottal dynamic parameters were calculated. The generated acoustic signal was evaluated by objective and subjective measures. The individual impairment was documented by a standardized questionnaire (VHI). Additionally, the amount of generated saliva was measured for a defined period of time.

RESULTS

ED patients displayed a significantly reduced saliva production compared to the control group. Furthermore, the auditory-perceptual evaluation yielded significantly higher ratings for breathiness and hoarseness in the voices of male ED patients compared to male controls. The majority of male ED patients (67%) indicated at least minor impairment in the self-evaluation. Objective acoustic measures like Jitter and Shimmer confirmed the decreased acoustic quality in male ED patients, whereas none of the investigated HSV parameters showed significant differences between the test groups. Statistical analysis did not confirm a statistically significant correlation between reduced voice quality and amount of saliva.

CONCLUSIONS

An objective impairment of the acoustic outcome was demonstrated for male ED patients. However, the vocal folds dynamics in HSV recordings seem unaffected.

Effect of Ligament Fibers on Dynamics of Synthetic, Self-Oscillating Vocal Folds in a Biomimetic Larynx Model

Synthetic silicone larynx models are essential for understanding the biomechanics of physiological and pathological vocal fold vibrations. The aim of this study is to investigate the effects of artificial ligament fibers on vocal fold vibrations in a synthetic larynx model, which is capable of replicating physiological laryngeal functions such as elongation, abduction, and adduction. A multi-layer silicone model with different mechanical properties for the musculus vocalis and the lamina propria consisting of ligament and mucosa was used. Ligament fibers of various diameters and break resistances were cast into the vocal folds and tested at different tension levels. An electromechanical setup was developed to mimic laryngeal physiology. The measurements included high-speed video recordings of vocal fold vibrations, subglottal pressure and acoustic. For the evaluation of the vibration characteristics, all measured values were evaluated and compared with parameters from ex and in vivo studies. The fundamental frequency of the synthetic larynx model was found to be approximately 200-520 Hz depending on integrated fiber types and tension levels. This range of the fundamental frequency corresponds to the reproduction of a female normal and singing voice range. The investigated voice parameters from vocal fold vibration, acoustics, and subglottal pressure were within normal value ranges from ex and in vivo studies. The integration of ligament fibers leads to an increase in the fundamental frequency with increasing airflow, while the tensioning of the ligament fibers remains constant. In addition, a tension increase in the fibers also generates a rise in the fundamental frequency delivering the physiological expectation of the dynamic behavior of vocal folds.

Influence of Perspective Distortion in Laryngoscopy

OBJECTIVE

An experiment with controllable boundaries was designed to assess the influence of the recording angle and distance on two-dimensional (2D) imaging in laryngoscopy and resulting 2D parameter calculation derived from the glottal area waveform (GAW).

METHOD

Two high-speed camera setups were used to synchronously record an oscillating synthetic vocal fold (VF) model, simulating a high-speed videoendoscopy. One camera recorded at variable lateral recording angles and a reference camera in superior perspective. This was performed at different physiological recording distances and for two oscillation modes (with/without contacting VFs). The GAW was derived from the segmented glottis, and two parameters each for the categories of symmetry, periodicity, and closure were calculated, as well as two derivative measures. The percentage difference between the variable and reference camera value pairs was calculated, and the angle and height dependencies were quantified using linear regression.

RESULTS

The visual perception of a laryngoscopy was found to be influenced by the lateral recording angle, which may lead to misinterpretation of VF symmetry among inexperienced observers. The strongest influence of recording angle was observed for symmetry parameters, the strongest being the Amplitude Symmetry Index with up to 2.6%/° (p < .05). A dependence on the recording distance was only found for the Maximum Area Declination Rate.

CONCLUSIONS

The recording angle in 2D laryngoscopy should be carefully considered during visual inspection of the VF dynamics. Most of the investigated objective parameters were unaffected by the examined perspective distortion. However, especially left-right symmetry measures should only be used under controlled boundary conditions to avoid misdiagnosis and misinterpretation.

SUPPLEMENTAL MATERIAL

https://doi.org/10.23641/asha.23961183.

The mechanisms of harmonic sound generation during phonation: A multi-modal measurement-based approach

Sound generation during voiced speech remains an open research topic because the underlying process within the human larynx is hardly accessible for direct measurements. In the present study, harmonic sound generation during phonation was investigated with a model that replicates the fully coupled fluid-structure-acoustic interaction (FSAI). The FSAI was captured using a multi-modal approach by measuring the flow and acoustic source fields based on particle image velocimetry, as well as the surface velocity of the vocal folds based on laser vibrometry and high-speed imaging. Strong harmonic sources were localized near the glottis, as well as further downstream, during the presence of the supraglottal jet. The strongest harmonic content of the vocal fold surface motion was verified for the area near the glottis, which directly interacts with the glottal jet flow. Also, the acoustic back-coupling of the formant frequencies onto the harmonic oscillation of the vocal folds was verified. These findings verify that harmonic sound generation is the result of a strong interrelation between the vocal fold motion, modulated flow field, and vocal tract geometry.

3D-FV-FE Aeroacoustic Larynx Model for Investigation of Functional Based Voice Disorders

For the clinical analysis of underlying mechanisms of voice disorders, we developed a numerical aeroacoustic larynx model, called simVoice, that mimics commonly observed functional laryngeal disorders as glottal insufficiency and vibrational left-right asymmetries. The model is a combination of the Finite Volume (FV) CFD solver Star-CCM+ and the Finite Element (FE) aeroacoustic solver CFS++. simVoice models turbulence using Large Eddy Simulations (LES) and the acoustic wave propagation with the perturbed convective wave equation (PCWE). Its geometry corresponds to a simplified larynx and a vocal tract model representing the vowel /a/. The oscillations of the vocal folds are externally driven. In total, 10 configurations with different degrees of functional-based disorders were simulated and analyzed. The energy transfer between the glottal airflow and the vocal folds decreases with an increasing glottal insufficiency and potentially reflects the higher effort during speech for patients being concerned. This loss of energy transfer may also have an essential influence on the quality of the sound signal as expressed by decreasing sound pressure level (SPL), Cepstral Peak Prominence (CPP), and Vocal Efficiency (VE). Asymmetry in the vocal fold oscillations also reduces the quality of the sound signal. However, simVoice confirmed previous clinical and experimental observations that a high level of glottal insufficiency worsens the acoustic signal quality more than oscillatory left-right asymmetry. Both symptoms in combination will further reduce the quality of the sound signal. In summary, simVoice allows for detailed analysis of the origins of disordered voice production and hence fosters the further understanding of laryngeal physiology, including occurring dependencies. A current walltime of 10 h/cycle is, with a prospective increase in computing power, auspicious for a future clinical use of simVoice.

Hybrid aeroacoustic approach for the efficient numerical simulation of human phonation

A hybrid aeroacoustic approach was developed for the efficient numerical computation of human phonation. In the first step, an incompressible flow simulation on a three-dimensional (3 D) computational grid, which is capable of resolving all relevant turbulent scales, is performed using STARCCM+ and finite volume method. In the second step, the acoustic source terms on the flow grid are computed and a conservative interpolation to the acoustic grid is performed. Finally, the perturbed convective wave equation is solved to obtain the acoustic field in 3 D with the finite element solver CFS++. Thereby, the conservative transformation of the acoustic sources from the flow grid to the acoustic grid is a key step to allow coarse acoustic grids without reducing accuracy. For this transformation, two different interpolation strategies are compared and grid convergence is assessed. Overall, 16 simulation setups are compared. The initial (267 000 degrees of freedom) and the optimized (21 265 degrees of freedom) simulation setup were validated by measurements of a synthetic larynx model. To conclude, the total computational time of the acoustic simulation is reduced by 95% compared to the initial simulation setup without a significant reduction of accuracy, being 7%, in the frequency range of interest.

Aerodynamic impact of the ventricular folds in computational larynx models

Ventricular folds (VeFs) act as passive, non-moving structures during normal phonation. According to the literature, VeFs potentially aid the flow-driven oscillations of the vocal folds (VFs) that produce the primary sound of human phonation. In this study, large eddy simulations were performed to analyze this influence in a numerical model with imposed VF motion as measured experimentally from a synthetic silicone vocal fold model. Model configurations with and without VeFs were considered. Furthermore, configurations with rectangular and elliptical glottis shapes were simulated to investigate the effects of three-dimensional glottal jet evolutions. Results showed that VeFs increased flow rate and transglottal pressure difference by a decrease in the pressure level in the ventricles immediately downstream of the VFs. This led to an increase in the glottal flow resistance, increased energy transfer rate between the flow and VFs, and a simultaneous decrease in the laryngeal flow resistance, which shows a higher amount of kinetic energy in the glottal flow. This enhancement was more pronounced in the rectangular glottis and varied with the subglottal pressure and VeF gap size.

Vibratory Characteristics of Diplophonia Studied by High Speed Video and Vibrogram Analysis

Diplophonia can occur in patients with polyps, atrophy, paralysis, or scars. Its vibratory patterns have not been well characterized. High-speed video (HSV) analysis can contribute to their understanding. Twenty subjects with a diplophonic voice quality were studied by HSV. Diplophonia was due to medical causes including vocal fold paresis (n = 7), vocal atrophy (n = 5), polyps (n = 5), and scars/sulci (n = 3). The HSV was analyzed using a multislice digital videokymography (DKG). The DKG tracing was analyzed qualitatively and then transformed into a vibrogram waveform signal for frequency analysis. RESULTS: Vibratory abnormalities seen on HSVs explained the diplophonia. Subharmonics to the fundamental frequency can be visualized by DKG. None could be resolved by stroboscopy. One can stratify diplophonia as symmetric or asymmetric based on the involvement of one or both vocal folds. Scars and atrophy showed symmetric subharmonic production with ectopic beats every 4-10 beats. Some subjects showed anterior and posterior independent vocal fold oscillators. Asymmetric causes of diplophonia are common in patients with paralysis. Two different oscillation frequencies of each vocal fold generate in and then out of phase interaction between the two sides. Vibrogram analysis documents the frequent presence of interharmonic energy peaks above the dominant fundamental frequency. Eighteen of the 20 subjects have obvious subharmonic peaks. CONCLUSION: Patients with diplophonia have vibratory abnormalities arising from the vocal folds. HSV and vibrogram analysis followed by frequency analysis of the vibrogram can resolve vibratory abnormality into symmetric versus asymmetric causes and can document the type of vibratory abnormality.

Computational Models of Laryngeal Aerodynamics: Potentials and Numerical Costs

Human phonation is based on the interaction between tracheal airflow and laryngeal dynamics. This fluid-structure interaction is based on the energy exchange between airflow and vocal folds. Major challenges in analyzing the phonatory process in-vivo are the small dimensions and the poor accessibility of the region of interest. For improved analysis of the phonatory process, numerical simulations of the airflow and the vocal fold dynamics have been suggested. Even though most of the models reproduced the phonatory process fairly well, development of comprehensive larynx models is still a subject of research. In the context of clinical application, physiological accuracy and computational model efficiency are of great interest. In this study, a simple numerical larynx model is introduced that incorporates the laryngeal fluid flow. It is based on a synthetic experimental model with silicone vocal folds. The degree of realism was successively increased in separate computational models and each model was simulated for 10 oscillation cycles. Results show that relevant features of the laryngeal flow field, such as glottal jet deflection, develop even when applying rather simple static models with oscillating flow rates. Including further phonatory components such as vocal fold motion, mucosal wave propagation, and ventricular folds, the simulations show phonatory key features like intraglottal flow separation and increased flow rate in presence of ventricular folds. The simulation time on 100 CPU cores ranged between 25 and 290 hours, currently restricting clinical application of these models. Nevertheless, results show high potential of numerical simulations for better understanding of phonatory process.