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

Asymmetric triangular body-cover model of the vocal folds with bilateral intrinsic muscle activation

Many voice disorders are linked to imbalanced muscle activity and known to exhibit asymmetric vocal fold vibration. However, the relation between imbalanced muscle activation and asymmetric vocal fold vibration is not well understood. This study introduces an asymmetric triangular body-cover model of the vocal folds, controlled by the activation of bilateral intrinsic laryngeal muscles, to investigate the effects of muscle imbalance on vocal fold oscillation. Various scenarios were considered, encompassing imbalance in individual muscles and muscle pairs, as well as accounting for asymmetry in lumped element parameters. Measurements of amplitude and phase asymmetries were employed to match the oscillatory behavior of two pathological cases: unilateral paralysis and muscle tension dysphonia. The resulting simulations exhibit muscle imbalance consistent with expectations in the composition of these voice disorders, yielding asymmetries exceeding 30% for paralysis and below 5% for dysphonia. This underscores the relevance of muscle imbalance in representing phonatory scenarios and its potential for characterizing asymmetry in vocal fold vibration.

Effects of implant and vocal fold stiffness on voice production after medialization laryngoplasty in an MRI-based vocal fold model

Medialization laryngoplasty is one of the primary surgical interventions in the treatment of glottal insufficiency due to vocal fold paralysis, paresis, or atrophy. During the surgery, an implant is laterally inserted into the larynx to medialize the affected vocal fold toward glottal midline, with the goal of improving glottal closure during phonation and voice production efficiency. While implants of different materials and geometry designs have been used, the effect of implant design on the voice outcome remains unclear. In this simulation study, the effect of implant stiffness was investigated in an MRI-based model of the vocal folds after medialization laryngoplasty. The results showed that implant stiffness had a significant impact on the phonation threshold pressure, glottal area waveform, and fundamental frequency, but only small effect on the closed quotient and other acoustic measures of the produced voice. The effect of implant stiffness also exhibited variability, depending on the stiffness conditions of the vocal fold and paraglottic tissues, indicating that individual differences need to be considered during the planning of medialization laryngoplasty.

Fluid-Structure Interaction Analysis of Aerodynamic and Elasticity Forces During Vocal Fold Vibration

The effect of the intraglottal vortices on the glottal flow waveform was explored using flow-structure-interaction (FSI) modeling. These vortices form near the superior aspect of the vocal folds during the closing phase of the folds' vibration. The geometry of the vocal fold was based on the well-known M5 model. The model did not include a vocal tract to remove its inertance effect on the glottal flow. Material properties for the cover and body layers of the folds were set using curve fit to experimental data of tissue elasticity. A commercially available FSI solver was used to perform simulations at low and high values of subglottal input pressure. Validation of the FSI results showed a good agreement for the glottal flow and the vocal fold displacement data with measurements taken in the excised canine larynx model. The simulations result further support the hypothesis that intraglottal vortices can affect the glottal flow waveform, specifically its maximum flow declination rate (MFDR). It showed that MFDR occurs at the same phase when the highest intraglottal vortical strength and the negative pressure occur. It also showed that when MFDR occurs, the magnitude of the aerodynamic force acting on the glottal wall is greater than the elastic recoil force predicted in the tissue. These findings are significant because nearly all theoretical and computational models that study the vocal fold vibrations mechanism do not consider the intraglottal negative pressure caused by the vortices as an additional closing force acting on the folds.

Flow-induced self-sustained oscillations in a straight channel with rigid walls and elastic supports

This work considers the two-dimensional flow field of an incompressible viscous fluid in a parallel-sided channel. In our study, one of the walls is fixed whereas the other one is elastically mounted, and sustained oscillations are induced by the fluid motion. The flow that forces the wall movement is produced as a consequence that one of the ends of the channel is pressurized, whereas the opposite end is at atmospheric pressure. The study aims at reducing the complexity of models for several physiological systems in which fluid-structure interaction produces large deformation of the wall. We report the experimental results of the observed self-sustained oscillations. These oscillations occur at frequencies close to the natural frequency of the system. The vertical motion is accompanied by a slight trend to rotate the moving mass at intervals when the gap height is quite narrow. We propose a simplified analytical model to explore the conditions under which this motion is possible. The analytical approach considers asymptotic solutions of the Navier-Stokes equation with a perturbation technique. The comparison between the experimental pressure measured at the midlength of the channel and the analytical result issued with a model neglecting viscous effects shows a very good agreement. Also, the rotating trend of the moving wall can be explained in terms of the quadratic dependence of the pressure with the streamwise coordinate that is predicted by this simplified model.

Subject-Specific Computational Fluid-Structure Interaction Modeling of Rabbit Vocal Fold Vibration

A full three-dimensional (3D) fluid-structure interaction (FSI) study of subject-specific vocal fold vibration is carried out based on the previously reconstructed vocal fold models of rabbit larynges. Our primary focuses are the vibration characteristics of the vocal fold, the unsteady 3D flow field, and comparison with a recently developed 1D glottal flow model that incorporates machine learning. The 3D FSI model applies strong coupling between the finite-element model for the vocal fold tissue and the incompressible Navier-Stokes equation for the flow. Five different samples of the rabbit larynx, reconstructed from the magnetic resonance imaging (MRI) scans after the in vivo phonation experiments, are used in the FSI simulation. These samples have distinct geometries and a different inlet pressure measured in the experiment. Furthermore, the material properties of the vocal fold tissue were determined previously for each individual sample. The results demonstrate that the vibration and the intraglottal pressure from the 3D flow simulation agree well with those from the 1D flow model based simulation. Further 3D analyses show that the inferior and supraglottal geometries play significant roles in the FSI process. Similarity of the flow pattern with the human vocal fold is discussed. This study supports the effective usage of rabbit larynges to understand human phonation and will help guide our future computational studies that address vocal fold disorders.

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.

The influence of flow model selection on finite element model parameter estimation using Bayesian inference

Recently, Bayesian estimation coupled with finite element modeling has been demonstrated as a viable tool for estimating vocal fold material properties from kinematic information obtained via high-speed video recordings. In this article, the sensitivity of the parameter estimations to the employed fluid model is explored by considering Bernoulli and one-dimensional viscous fluid flow models. Simulation results indicate that prescribing an ad hoc separation location for the Bernoulli flow model can lead to large estimate biases, whereas including the separation location as an estimated parameter leads to results comparable to that of the viscous fluid flow model.

A one-dimensional flow model enhanced by machine learning for simulation of vocal fold vibration

A one-dimensional (1D) unsteady and viscous flow model that is derived from the momentum and mass conservation equations is described, and to enhance this physics-based model, a machine learning approach is used to determine the unknown modeling parameters. Specifically, an idealized larynx model is constructed and ten cases of three-dimensional (3D) fluid-structure interaction (FSI) simulations are performed. The flow data are then extracted to train the 1D flow model using a sparse identification approach for nonlinear dynamical systems. As a result of training, we obtain the analytical expressions for the entrance effect and pressure loss in the glottis, which are then incorporated in the flow model to conveniently handle different glottal shapes due to vocal fold vibration. We apply the enhanced 1D flow model in the FSI simulation of both idealized vocal fold geometries and subject-specific anatomical geometries reconstructed from the magnetic resonance imaging images of rabbits' larynges. The 1D flow model is evaluated in both of these setups and shown to have robust performance. Therefore, it provides a fast simulation tool that is superior to the previous 1D models.

Numerical and experimental investigations on vocal fold approximation in healthy and simulated unilateral vocal fold paralysis

We have developed a novel surgical/computational model for the investigation of unilateral vocal fold paralysis (UVFP) which will be used to inform future in silico approaches to improve surgical outcomes in type I thyroplasty. Healthy phonation (HP) was achieved using cricothyroid suture approximation on both sides of the larynx to generate symmetrical vocal fold closure. Following high-speed videoendoscopy (HSV) capture, sutures on the right side of the larynx were removed, partially releasing tension unilaterally and generating asymmetric vocal fold closure characteristic of UVFP (sUVFP condition). HSV revealed symmetric vibration in HP, while in sUVFP the sutured side demonstrated a higher frequency (10 - 11%). For the computational model, ex vivo magnetic resonance imaging (MRI) scans were captured at three configurations: non-approximated (NA), HP, and sUVFP. A finite-element method (FEM) model was built, in which cartilage displacements from the MRI images were used to prescribe the adduction and the vocal fold deformation was simulated before the eigenmode calculation. The results showed that the frequency comparison between the two sides were consistent with observations from HSV. This alignment between the surgical and computational models supports the future application of these methods for the investigation of treatment for UVFP.