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.     

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.     

An investigation of bimodal jet trajectory in flow through scaled models of the human vocal tract

Pulsatile two-dimensional flow through static divergent models of the human vocal folds is investigated. Although the motivation for this study is speech production, the results are generally applicable to a variety of engineering flows involving pulsatile flow through diffusers. Model glottal divergence angles of 10, 20, and 40° represent various geometries encountered in one phonation cycle. Frequency and amplitude of the flow oscillations are scaled with physiological Reynolds and Strouhal numbers typical of human phonation. Glottal velocity trajectories are measured along the anterior–posterior midline by using phase-averaged particle image velocimetry to acquire 1,000 realizations at ten discrete instances in the phonation cycle. The angular deflection of the glottal jet from the streamwise direction (symmetric configuration) is quantified for each realization. A bimodal flow configuration is observed for divergence angles of 10 and 20°, with the flow eventually skewing and attaching to the vocal fold walls. The deflection of the flow toward the vocal fold walls occurs when the forcing function reaches maximum velocity and zero acceleration. For a divergence angle of 40°, the flow never attaches to the vocal fold walls; however, there is increased variability in the glottal jet after the forcing function reaches maximum velocity and zero acceleration. The variation in the jet trajectory as a function of divergence angle is explained by performance maps of diffuser flow regimes. The smaller angle cases are in the unstable transitory stall regime while the 40° divergent case is in the fully developed two-dimensional stall regime. Very small geometric variations in model size and surface finish significantly affect the flow behavior. The bimodal, or flip-flopping, glottal jet behavior is expected to influence the dipole contribution to sound production.     

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.     

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.     

Solving for unsteady airflow in a glottal model with immersed moving boundaries

This work builds on previous efforts to characterize the dynamic development of the airflow in the glottis from a fluid mechanical point of view. A multigrid finite-difference method with immersed boundaries is implemented to solve the Navier–Stokes equations in a channel constricted by a vibrating rigid structure with a shape conforming to the human vocal folds. For the dynamically evolving boundaries we apply a forced oscillation glottal model. The large scale deformations of the boundaries are handled without regridding and tracheal input velocity is either set to a constant value or synchronized with wall motion. Particular attention is paid to the mobility of the point where the airflow detaches from the flapping walls. Results illustrate the relevance and the diversity of flow separation dynamics within the constriction standing for the glottis, while flow instabilities past the constriction are not found to affect flow behavior between the moving walls significantly. A comparison between static and dynamic numerical experiments show that mobility of the flow separation point is nontrivial in general and only rarely quasi-static.     

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.     

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.     

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.     

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.     

Existence of Global Strong Solutions to a Beam–Fluid Interaction System

We study an unsteady nonlinear fluid–structure interaction problem which is a simplified model to describe blood flow through viscoelastic arteries. We consider a Newtonian incompressible two-dimensional flow described by the Navier–Stokes equations set in an unknown domain depending on the displacement of a structure, which itself satisfies a linear viscoelastic beam equation. The fluid and the structure are fully coupled via interface conditions prescribing the continuity of the velocities at the fluid–structure interface and the action–reaction principle. We prove that strong solutions to this problem are global-in-time. We obtain, in particular that contact between the viscoelastic wall and the bottom of the fluid cavity does not occur in finite time. To our knowledge, this is the first occurrence of a no-contact result, and of the existence of strong solutions globally in time, in the frame of interactions between a viscous fluid and a deformable structure.     

Analysis of flow-structure coupling in a mechanical model of the vocal folds and the subglottal system

An analysis is made of the nonlinear interactions between flow in the subglottal vocal tract and glottis, sound waves in the subglottal system and a mechanical model of the vocal folds. The mean flow through the system is produced by a nominally steady contraction of the lungs, and mechanical experiments frequently involve a ‘lung cavity’ coupled to an experimental subglottal tube of arbitrary or ill-defined effective length L, on the basis that the actual value of L has little or no influence on excitation of the vocal folds. A simple, self-exciting single-mass mathematical model of the vocal folds is used to investigate the sound generated within the subglottal domain and the unsteady volume flux from the glottis for experiments where it is required to suppress feedback of sound from the supraglottal vocal tract. In experiments where the assumed absorption of sound within the sponge-like interior of the lungs is small, the influence of changes in L can be very significant: when the subglottal tube behaves as an open-ended resonator (when L is as large as half the acoustic wavelength) there is predicted to be a mild increase in volume flux magnitude and a small change in waveform. However, the strong appearance of second harmonics of the acoustic field is predicted at intermediate lengths, when L is roughly one quarter of the acoustic wavelength. In cases of large lung damping, however, only modest changes in the volume flux are predicted to occur with variations in L.     

Finite element computations for unsteady fluid and elastic membrane interaction problems

The interaction between the hydrodynamic forces of a flow field and the elastic forces of adjacent deformable boundaries is described by elastohydrodynamics, a coupled fluid–elastic membrane problem. Direct numerical solution of the unsteady, highly non-linear equations requires that the dynamic evolution of both the flow field and the domain shape be determined as part of the solution, since neither is known a priori. This paper describes a numerical algorithm based on the deformable spatial domain space–time (DSD/ST) finite element method for the unsteady motion of an incompressible, viscous fluid with elastic membrane interaction. The unsteady Navier–Stoke and elastic membrane equations are solved separately using an iterative procedure by the GMRES technique with an incomplete lower-upper (ILU) decomposition at every time instant. One-dimensional, two-dimensional and deformable domain model problems are used to demonstrate the capabilities and accuracy of the present algorithm. Both steady state and transient problems are studied. © 1997 John Wiley & Sons, Ltd.     

PIV measurements combined with the motion tracking technique to analyze flow around a moving porous structure

To gain a better understanding of the fluid–structure interaction and especially when dealing with a flow around an arbitrarily moving body, it is essential to develop measurement tools enabling the instantaneous detection of moving deformable interface during the flow measurements. A particularly useful application is the determination of unsteady turbulent flow velocity field around a moving porous fishing net structure which is of great interest for selectivity and also for the numerical code validation which needs a realistic database. To do this, a representative piece of fishing net structure is used to investigate both the Turbulent Boundary Layer (TBL) developing over the horizontal porous moving fishing net structure and the turbulent flow passing through the moving porous structure. For such an investigation, Time Resolved PIV measurements are carried out and combined with a motion tracking technique allowing the measurement of the instantaneous motion of the deformable fishing net during PIV measurements. Once the two-dimensional motion of the porous structure is accessed, PIV velocity measurements are analyzed in connection with the detected motion. Finally, the TBL is characterized and the effect of the structure motion on the volumetric flow rate passing though the moving porous structure is clearly demonstrated.     

A computational framework for fluid–porous structure interaction with large structural deformation

We study the effect of poroelasticity on fluid–structure interaction. More precisely, we analyze the role of fluid flow through a deformable porous matrix in the energy dissipation behavior of a poroelastic structure. For this purpose, we develop and use a nonlinear poroelastic computational model and apply it to the fluid–structure interaction simulations. We discretize the problem by means of the finite element method for the spatial approximation and using finite differences in time. The numerical discretization leads to a system of non-linear equations that are solved by Newton’s method. We adopt a moving mesh algorithm, based on the Arbitrary Lagrangian–Eulerian method to handle large deformations of the structure. To reduce the computational cost, the coupled problem of free fluid, porous media flow and solid mechanics is split among its components and solved using a partitioned approach. Numerical results show that the flow through the porous matrix is responsible for generating a hysteresis loop in the stress versus displacement diagrams of the poroelastic structure. The sensitivity of this effect with respect to the parameters of the problem is also analyzed.

     

Simulation of forced deformable bodies interacting with two-dimensional incompressible flows: Application to fish-like swimming

We present an efficient algorithm for simulation of deformable bodies interacting with two-dimensional incompressible fluid flows. The temporal and spatial discretizations of the Navier–Stokes equations in vorticity stream-function formulation are based on classical fourth-order Runge–Kutta scheme and compact finite differences, respectively. Using a uniform Cartesian grid we benefit from the advantage of a new fourth-order direct solver for the Poisson equation to ensure the incompressibility constraint down to machine zero over an optimal grid. For introducing a deformable body in fluid flow, the volume penalization method is used. A Lagrangian structured grid with prescribed motion covers the deformable body which is interacting with the surrounding fluid due to the hydrodynamic forces and the torque calculated on the Eulerian reference grid. An efficient law for controlling the curvature of an anguilliform fish, swimming toward a prescribed goal, is proposed which is based on the geometrically exact theory of nonlinear beams and quaternions. Validation of the developed method shows the efficiency and expected accuracy of the algorithm for fish-like swimming and also for a variety of fluid/solid interaction problems.     

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.     

Motions of elastic solids in fluids under vibration

Motion of a rigid or deformable solid in a viscous incompressible fluid and corresponding fluid–solid interactions are considered. Different cases of applying high frequency vibrations to the solid or to the surrounding fluid are treated. Simple formulas for the mean velocity of the solid are derived, under the assumption that the regime of the fluid flow induced by its motion is turbulent and the fluid resistance force is nonlinearly dependent on its velocity. It is shown that vibrations of a fluid’s volume slow down the motion of a submerged solid. This effect is much pronounced in the case of a deformable solid (i.e., gas bubble) exposed to near-resonant excitation. The results are relevant to the theory of gravitational enrichment of raw materials, and also contribute to the theory of controlled locomotion of a body with an internal oscillator in continuous deformable (solid or fluid) media.     

Shear deformable composite plates with nonlinear curvatures: modeling and nonlinear vibrations of symmetric laminates

Moving from a general plate theory, a modified general shear deformable laminated plate theory (MGFSDT) exhibiting nonlinear curvatures but still allowing for some worth features of linear curvature models (von Karman) is formulated. Starting from MGFSDT partial differential equations, a minimal discretized model (Duffing equation) for symmetric cross-ply laminates nonlinear vibrations, whose coefficients account for shear deformability and nonlinear curvatures, is obtained via the Galerkin procedure. The variable features of such a Duffing model as obtainable via alternative kinematic assumptions at the continuum level are highlighted. Through the comparison of a number of underlying models in different technical situations, information on the influence of shear deformability on system nonlinear response and on the influence of nonlinear curvatures are obtained. Frequency–response curves through a multiple scale analysis are presented for different continuous models, kinds of material, mode numbers and boundary conditions.     

Numerical simulation of a flow around an impulsively started radially deforming circular cylinder

The development of a two‐dimensional viscous incompressible flow generated by a deformable circular cylinder impulsively started into rectilinear motion is studied numerically for the Reynolds numbers equal to 550 and 3000. The vorticity transport equation is solved by a second‐order finite difference method in both directions of the domains. The Poisson equation for the streamfunction is solved by a Fourier–Galerkin method in the direction of the flow that is assumed to remain symmetrical and a second‐order finite difference for the radial direction. The advance in time is achieved by a second‐order Adams–Bashforth scheme. The computed results are compared qualitatively with experimental and numerical results done before in the particular non‐deformable case. The comparison is found to be satisfactory. The influence of the deformation of the cylinder on the flow structure and the drag coefficient is then analyzed. Copyright © 1999 John Wiley & Sons, Ltd.