Computational Investigations of the Fluid-Structure Interaction During Phonation: The Role of Vocal Fold Elasticity and Glottal Flow Unsteadiness

Human voice production arises from the biomechanical interaction between vocal fold vibrations and airflow dynamics. Changes in vocal fold stiffness can lead to changes in vocal fold vibration patterns and further changes in voice outcomes. A good knowledge of the cause-and-effect relationship between vocal fold stiffness and voice production can not only deepen the understanding of voice production mechanisms but also benefit the treatment of voice disorders associated with vocal fold stiffness changes. This constitutes the first objective of this dissertation. The second objective of this dissertation is to further examine the range of validity of the quasi-steady assumption of glottal flow during phonation. The assumption is of vital importance for phonation modeling since it enables to eliminate the unsteady aspects of glottal flow, which greatly simplifies the flow modeling. A three-dimensional flow-structure interaction model of voice production is employed to investigate the effects of vocal fold stiffness parameters on voice production. The vocal fold is modeled as the cover-ligament-body structure with a transversely isotropic constitutive relation. Stiffness parameters in both the transverse plane and the longitudinal direction of each layer of the vocal fold are systematically varied. The results show that varying the stiffness parameters has obvious monotonic effects on the fundamental frequency, glottal flow rate and glottal opening, but has non-monotonic effects on the glottal divergent angle, open quotient and closing velocity. Compared to the transverse stiffness parameters, the longitudinal stiffness parameters generally have more significant impacts on glottal flows and vocal fold vibrations. Additionally, the sensitivity analysis reveals that the stiffness parameters of the ligament layer have the largest effect on most output measures. Next, flow-structure interaction simulations are carried out to study the effect of fiber orientation in the conus elasticus on voice production. Two continuum vocal fold models with different fiber orientations in the conus elasticus are constructed. The more realistic fiber orientation (caudal-cranial) in the conus elasticus is found to yield smaller structural stiffness and larger deflection at the junction of the conus elasticus and ligament than the anterior-posterior fiber orientation, which facilitates vocal fold vibrations and eventually causes a larger peak flow rate and higher speed quotient. The generated voice is also found to have a lower fundamental frequency and smaller spectral slope. Finally, the validity of the quasi-steady assumption for glottal flow is systematically examined by considering the voice frequency range, complexity of glottal shapes and air inertia in the vocal tract. The results show that at the normal speech frequency (~ 100 Hz), the dynamics of the quasi-steady flow greatly resembles that of a dynamic flow, and the glottal flow and glottal pressure predicted by the quasi-steady approximation have very small errors. However, the assumption produces huge errors at high frequencies (~ 500 Hz). In addition, air inertia in the vocal tract can undermine the validity of the assumption via the nonlinear interaction with the unsteady glottal flow. The role of glottal shapes in the validation is found to be insignificant

Constricted channel flow with different cross-section shapes

Pressure driven steady flow through a uniform circular channel containing a constricted portion is a common problem considering physiological flows such as underlying human speech sound production. The influence of the constriction’s cross-section shape (circle, ellipse, circular sector) on the flow within and downstream from the constriction is experimentally quantified. An analytical boundary layer flow model is proposed which takes into account the hydraulic diameter of the cross-section shape. Comparison of the model outcome with experimental and three-dimensional numerically simulated flow data shows that the pressure distribution within the constriction can be modeled accurately so that the model is of interest for analytical models of fluid–structure interaction without the assumption of two-dimensional flow

Measurement of flow separation in a human vocal folds model

International audienceThe 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

Validation of a flow-structure-interaction computation model of phonation

Computational models of vocal fold (VF) vibration are becoming increasingly sophisticated, their utility currently transiting from exploratory research to predictive research. However, validation of such models has remained largely qualitative, raising questions over their applicability to interpret clinical situations. In this paper, a computational model with a segregated implementation is detailed. The model is used to predict the fluid–structure interaction(FSI) observed in a physical replica of the VFs when it is excited by airflow. Detailed quantitative comparisons are provided between the computational model and the corresponding experiment. First, the flow model is separately validated in the absence of VF motion. Then, in the presence of flow-induced VF motion, comparisons are made of the flow pressure on the VF walls and of the resulting VF displacements. Self-similarity of spatial distributions of flow pressure and VF displacements is highlighted. The self-similarity leads to normalized pressure and displacement profiles. It is shown that by using linear superposition of average and fluctuation components of normalized computed displacements, it is possible to determine displacements in the physical VF replica over a range of VF vibration conditions. Mechanical stresses in the VF interior are related to the VF displacements, thereby the computational model can also determine VF stresses over a range of phonation conditions

High Fidelity Computational Modeling and Analysis of Voice Production

This research aims to improve the fundamental understanding of the multiphysics nature of voice production, particularly, the dynamic couplings among glottal flow, vocal fold vibration and airway acoustics through high-fidelity computational modeling and simulations. Built upon in-house numerical solvers, including an immersed-boundary-method based incompressible flow solver, a finite element method based solid mechanics solver and a hydrodynamic/aerodynamic splitting method based acoustics solver, a fully coupled, continuum mechanics based fluid-structure-acoustics interaction model was developed to simulate the flow-induced vocal fold vibrations and sound production in birds and mammals. Extensive validations of the model were conducted by comparing to excised syringeal and laryngeal experiments. The results showed that, driven by realistic representations of physiology and experimental conditions, including the geometries, material properties and boundary conditions, the model had an excellent agreement with the experiments on the vocal fold vibration patterns, acoustics and intraglottal flow dynamics, demonstrating that the model is able to reproduce realistic phonatory dynamics during voice production. The model was then utilized to investigate the effect of vocal fold inner structures on voice production. Assuming the human vocal fold to be a three-layer structure, this research focused on the effect of longitudinal variation of layer thickness as well as the cover-body thickness ratio on vocal fold vibrations. The results showed that the longitudinal variation of the cover and ligament layers thicknesses had little effect on the flow rate, vocal fold vibration amplitude and pattern but affected the glottal angle in different coronal planes, which also influenced the energy transfer between glottal flow and the vocal fold. The cover-body thickness ratio had a complex nonlinear effect on the vocal fold vibration and voice production. Increasing the cover-body thickness ratio promoted the excitation of the wave-type modes of the vocal fold, which were also higher-eigenfrequency modes, driving the vibrations to higher frequencies. This has created complex nonlinear bifurcations. The results from the research has important clinical implications on voice disorder diagnosis and treatment as voice disorders are often associated with mechanical status changes of the vocal fold tissues and their treatment often focus on restoring the mechanical status of the vocal folds

Numerical Study of Laryngeal Control of Phonation using Realistic Finite Element Models of a Canine Larynx

While many may take it for granted, the human voice is an incredible feat. An average person can produce a great variety of voices and change voice characteristics agilely even without formal training. Last several decades of research has established that the production of voice is largely a mechanical process: i.e., the sustained vibration of the vocal folds driven by the glottal air flow. Since one only has a single pair of vocal folds, the versatility comes with the ability to change the mechanical status of the vocal folds, including vocal fold length and thickness, tension, and level of adduction, through activation of the laryngeal muscles. However, the relationship between laryngeal muscle activity and the characteristics of voice is not well understood due to limitations in experimental observation and simplifications in modelling and simulations. The science is still far behind the art. The current research aims to investigate first the relationship between laryngeal muscle activation and the posture of the vocal folds and second the relationship between voice source characteristics and vocal fold mechanical status using more comprehensive numerical models and simulations, thus improving the understanding of the roles of each laryngeal muscle in voice control. To do so, (1) the mechanics involved in vocal fold posturing and vibration, especially muscle contraction; (2) the realistic anatomical structure of the larynx must be considered properly. To achieve this goal, a numerical model of the larynx as realistic as possible was built. The geometry of the laryngeal components was reconstructed from high resolution MRI (Magnetic Resonance Imaging) data of an excised canine larynx, which makes more accurate the representation of the muscles and their sub-compartments, cartilages, and other important anatomical features of the larynx. A previously proposed muscle activation model was implemented in a 3D finite element package and applied to the larynx model to simulate the action of laryngeal muscles. After validation of the numerical model against experimental data, extensive parametric studies involving different combination of muscle activations were conducted to investigate how the voice source is controlled with laryngeal muscles. In the course of this study, some work was done to couple the same finite element tool with a Genetic Algorithm program to inversely determine model parameters in biomechanical models. The method was applied in a collaborated study on shape changes of a fish fin during swimming. This study is presented as a separate chapter at the end of this thesis. The method has potential application in determining parameters in vocal fold models and optimizing clinical vocal fold procedures. This thesis is essentially an assembly of the papers published by the author during the doctoral study, with the addition of an introductory chapter. Chapter 1 reviews the overall principles of voice production, the biomechanical basis of voice control, and past studies on voice control with a focus on the fundamental frequency. Chapter 2 describes the major numerical methods employed in this research with an emphasis on the finite element method. The muscle activation model is also described in this chapter. Chapter 3 describes the building of the larynx model from MRI data and its partial validation. Chapter 4 presents the application of the larynx model to posturing studies, including parametric activation of muscle groups and specific topics related to vocal fold posturing. Chapter 5 describes the change of vocal fold vibration dynamics under the influence of the interaction of the cricothyroid muscle and the thyroarytenoid muscle. The Flow-structure interaction simulations was realized by coupling the larynx model to a simple Bernoulli flow model and a two-stage simulation technique. Chapter 6 concludes the current thesis study. Suggestions for future studies are proposed. Chapter 7 is an independent study that is not related to voice control. It describes a numerical framework that inversely determines and validates model parameters of biomechanical models. The application of the proposed framework to a finite element model of a fish fin is presented

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