Numerical Analysis of the Sulcus Vocalis Disorder on the Function of the Vocal Folds

2016 ◽  
Vol 33 (4) ◽  
pp. 513-520 ◽  
Author(s):  
A. Vazifehdoostsaleh ◽  
N. Fatouraee ◽  
M. Navidbakhsh ◽  
F. Izadi

AbstractSpeaking is a very complex process resulting from the interaction between the air flow along the larynx and the vibrating structure of the vocal folds. Sulcus is a disease missing layers in the vocal folds result in cracks resulting in some disorders in producing sounds. Sulcus and its effects on the vocal cord vibrations are numerically studied for the first time in this paper. An ideal model of healthy vocal folds and Sulcus vocalis has been two-dimensionally defined and the finite element model of vocal folds is solved in a fully coupled form. The proposed calculative model was used in a fluid range of the computational fluid dynamics, arbitrary Lagrangian-Eulerian (ALE), incompressible continuity and Navier-Stokes equations and in a structure range of a three-layer elastic linear model. Self-excited oscillations were presented for vocal folds among type II patients and compared with healthy models. Responses were qualitatively and quantitatively studied. The healthy model was compared with numerical and empirical results. In addition, the effects of the disease on the flow parameters and the vibration frequency of the vocal folds were studied. According to the simulated model, the oscillation frequency decreased 25% and the average and instantaneous volume flux significantly increased compared to healthy samples. Results may help present a guideline for surgery and subsequently evaluate patients’ improvement.

2018 ◽  
Vol 34 (6) ◽  
pp. 791-800 ◽  
Author(s):  
A. Vazifehdoostsaleh ◽  
N. Fatouraee ◽  
M. Navidbakhsh ◽  
F. Izadi

AbstractThe effect of sulcus vocalis on vocal folds function is investigated. A type II sulcus vocalis is defined, parameterized and incorporated into a three-dimensional, fully coupled finite element model of vocal folds and laryngeal airway. The proposed Fluid-Structure Interaction (FSI) model is utilized in computational fluid dynamics, Arbitrary Lagrangian-Eulerian (ALE), incompressible continuity and Navier-Stokes equations and in a structure range of a three-layer elastic linear model. Flow parameters, vibration behavior and glottal jet aerodynamics of healthy and patient vocal folds models are compared with each other. Flow visualization is utilized to characterize Coanda effect and three dimensionality of flow patterns. The vibration frequency of vocal folds having sulcus vocalis decreases in comparison with that of healthy ones. Upon increasing the volume flux in the sulcus vocalis model, the non-periodic and disordered behavior of it is visible for patient vocal folds. Underlying mechanisms for the observed changes, possible implications for treatments of sulcus vocalis and human perfect voice production are also discussed.


2021 ◽  
Vol 11 (11) ◽  
pp. 4748
Author(s):  
Monika Balázsová ◽  
Miloslav Feistauer ◽  
Jaromír Horáček ◽  
Adam Kosík

This study deals with the development of an accurate, efficient and robust method for the numerical solution of the interaction of compressible flow and nonlinear dynamic elasticity. This problem requires the reliable solution of flow in time-dependent domains and the solution of deformations of elastic bodies formed by several materials with complicated geometry depending on time. In this paper, the fluid–structure interaction (FSI) problem is solved numerically by the space-time discontinuous Galerkin method (STDGM). In the case of compressible flow, we use the compressible Navier–Stokes equations formulated by the arbitrary Lagrangian–Eulerian (ALE) method. The elasticity problem uses the non-stationary formulation of the dynamic system using the St. Venant–Kirchhoff and neo-Hookean models. The STDGM for the nonlinear elasticity is tested on the Hron–Turek benchmark. The main novelty of the study is the numerical simulation of the nonlinear vocal fold vibrations excited by the compressible airflow coming from the trachea to the simplified model of the vocal tract. The computations show that the nonlinear elasticity model of the vocal folds is needed in order to obtain substantially higher accuracy of the computed vocal folds deformation than for the linear elasticity model. Moreover, the numerical simulations showed that the differences between the two considered nonlinear material models are very small.


2003 ◽  
Vol 3 ◽  
pp. 195-207
Author(s):  
A.M. Ilyasov ◽  
V.N. Kireev ◽  
S.F. Urmancheev ◽  
I.Sh. Akhatov

The work is devoted to the analysis of the flow of immiscible liquid in a flat channel and the creation of calculation schemes for determining the flow parameters. A critical analysis of the well-known Two Fluids Model was carried out and a new scheme for the determination of wall and interfacial friction, called the hydraulic approximation in the theory of stratified flows, was proposed. Verification of the proposed approximate model was carried out on the basis of a direct numerical solution of the Navier–Stokes equations for each fluid by a finite-difference method with phase-boundary tracking by the VOF (Volume of Fluid) method. The graphical dependencies illustrating the change in the interfase boundaries of liquids and the averaged over the occupied area of the phase velocities along the flat channel are presented. The results of comparative calculations for two-fluid models are also given, according to the developed model in the hydraulic approximation and direct modeling. It is shown that the calculations in accordance with the hydraulic approximation are more consistent with the simulation results. Thus, the model of hydraulic approximation is the most preferred method for calculating stratified flows, especially in cases of variable volumetric content of liquids.


Author(s):  
Elisabeth Longatte

This work is concerned with the modelling of the interaction of a fluid with a rigid or a flexible elastic cylinder in the presence of axial or cross-flow. A partitioned procedure is involved to perform the computation of the fully-coupled fluid solid system. The fluid flow is governed by the incompressible Navier-Stokes equations and modeled by using a fractional step scheme combined with a co-located finite volume method for space discretisation. The motion of the fluid domain is accounted for by a moving mesh strategy through an Arbitrary Lagrangian-Eulerian (ALE) formulation. Solid dyncamics is modeled by a finite element method in the linear elasticity framework and a fixed point method is used for the fluid solid system computation. In the present work two examples are presented to show the method robustness and efficiency.


2005 ◽  
Vol 9 (1) ◽  
pp. 67-78 ◽  
Author(s):  
J. Socolowsky

iscous two‐fluid channel flows arise in different kinds of coating technologies. The corresponding mathematical models represent two‐dimensional free boundary value problems for the Navier‐Stokes equations. In this paper the solvability of the related stationary problems is discussed and computational results are presented. Furthermore, it is shown that depending on the flow parameters like viscosity or density ratios and on the fluxes there can happen nonexistence of steady‐state solutions. For other parameter sets the solution is even unique. Dvieju, tekančiu kanale, klampiu skysčiu srauto uždavinys iškyla taikant ivairias skirtingu rušiu paviršiu padengimo technologijas. Atitinkamas matematinis modelis išreiškiamas dvimačiu kraštiniu uždaviniu su laisvu paviršiumi Navje-Stokso lygtims. Straipsnyje nagrinejamas santykinai stacionaraus uždavinio išsprendžiamumas ir pateikiami skaičiavimo rezultatai. Be to parodoma, kad priklausomai nuo sroves parametru kaip ir nuo klampumo ir tankio santykio stacionarus sprendiniai gali neegzistuoti. Su kitais parametrais egzistuoja tiksliai vienas sprendinys.


Author(s):  
Dongho Shin ◽  
John C. Strikwerda

AbstractWe consider several methods for solving the linear equations arising from finite difference discretizations of the Stokes equations. The two best methods, one presented here for the first time, apparently, and a second, presented by Bramble and Pasciak, are shown to have computational effort that grows slowly with the number of grid points. The methods work with second-order accurate discretizations. Computational results are shown for both the Stokes equations and incompressible Navier-Stokes equations at low Reynolds number.


Author(s):  
Jaromi´r Hora´cˇek ◽  
Miloslav Feistauer ◽  
Petr Sva´cˇek

The contribution deals with the numerical simulation of the flutter of an airfoil with three degrees of freedom (3-DOF) for rotation around an elastic axis, oscillation in the vertical direction and rotation of a flap. The finite element (FE) solution of two-dimensional (2-D) incompressible Navier-Stokes equations is coupled with a system of nonlinear ordinary differential equations describing the airfoil vibrations with large amplitudes taking into account the nonlinear mass matrix. The time-dependent computational domain and a moving grid are treated by the Arbitrary Lagrangian-Eulerian (ALE) method and a suitable stabilization of the FE discretization is applied. The developed method was successfully tested by the classical flutter computation of the critical flutter velocity using NASTRAN program considering the linear model of vibrations and the double-lattice aerodynamic theory. The method was applied to the numerical simulations of the post flutter regime in time domain showing Limit Cycle Oscillations (LCO) due to nonlinearities of the flow model and vibrations with large amplitudes. Numerical experiments were performed for the airfoil NACA 0012 respecting the effect of the air space between the flap and the main airfoil.


1994 ◽  
Vol 116 (4) ◽  
pp. 870-876 ◽  
Author(s):  
Ko-Foa Tchon ◽  
Ion Paraschivoiu

In order to study the dynamic stall phenomenon on a Darrieus wind turbine, the incompressible flow field around a moving airfoil is simulated using a noninertial stream function-vorticity formulation of the two-dimensional unsteady Navier-Stokes equations. Spatial discretization is achieved by the streamline upwind Petrov-Galerkin finite element method on a hybrid mesh composed of a structured region of quadrilateral elements in the vicinity of solid boundaries, an unstructured region of triangular elements elsewhere, and a layer of infinite elements surrounding the domain and projecting the external boundary to infinity. Temporal discretization is achieved by an implicit second order finite difference scheme. At each time step, a nonlinear algebraic system is solved by a Newton method. To accelerate computations, the generalized minimum residual method with an incomplete triangular factorization preconditioning is used to solve the linearized Newton systems. The solver is applied to simulate the flow around a NACA 0015 airfoil in Darrieus motion and the results are compared to experimental observations. To the authors’ knowledge, it is the first time that the simulation of such a motion has been performed using the Navier-Stokes equations.


Sign in / Sign up

Export Citation Format

Share Document