|
|
On numerical approximation of fluid-structure interactions of air flow with a model of vocal folds
Valášek, J. ; Horáček, Jaromír ; Sváček, P.
This paper deals with flow driven vibration of an elastic body. Our goal is to develop and mathematically describe a simplified model of the human vocal fold. The developed numerical schemes for viscous incompressible fluid flow in ALE formulation and the elastic body are implemented by two solvers, specific for each domain. The studied problem is coupled by Dirichlet-Neumann boundary conditions. Both solvers are based on the finite element method. Particularly, for the fluid model the crossgrid elements are used. Numerical results focus on the verification of the developed program.
|
|
|
Large eddy simulation of airflow in human vocal folds
Šidlof, Petr
Human phonation is a complex physiological process involving flow-induced oscillations of the vocal folds and aeroacoustic sound generation. The flow fields encountered in phonation are highly unsteady, feature massive flow separation and recirculation in the supraglottal spaces and generation of coherent vortex structures from the shear layer of the jet. For the sake of computational aeroacoustic modeling of human voice generation, an accurate resolution of the airflow through the vocal folds is essential. The Reynolds-averaged Navier-Stokes turbulence modeling is inappropriate, since it provides only the averaged flow field. The paper presents the first results obtained with a large-eddy simulation of flow through a model of human vocal folds using a second-order finite volume discretization of incompressible Navier-Stokes equations. In the first step, the flow field was resolved on a fine 2D mesh covering a short subglottal region, the glottis and a part of the supraglottal channel. The simulation was parallelized using domain decomposition method and run in parallel on a shared-memory supercomputer. The results compare two large eddy simulations using the algebraic Smagorinsky and one-equation sub-grid scale models against a simulation without turbulence model.
|
|
|
On application of finite element method for approximation of 3D flow problems
Sváček, P. ; Horáček, Jaromír
This paper is interested to the interactions of the incompressible flow with a flexibly supported airfoil. The bending and the torsion modes are considered. The problem is mathematically described. The numerical method is based on the finite element method. A combination of the streamline-upwind/Petrov-Galerkin and pressure stabilizing/Petrov-Galerkin method is used for the stabilization of the finite element method. The numerical results for a three-dimensional problem of flow over an airfoil are shown.
|
|
|
Finite element simulation of aeroelasticity problems
Horáček, Jaromír ; Sváček, P. ; Feistauer, M.
The paper presents results achieved by the authors in development of in-house codes based the finite element (FE) method and applied to solution of fluid-structure problems in aeroelasticity of airfoils. We consider flexibly supported airfoil with two- or three degrees of freedom (2- or3-DOF) in two-dimensional (2D) incompressible viscous flow. The airfoil vibration is described by nonlinear ordinary differential equations of motion for large vibration amplitudes. The flow is medeled by the Navier-Stokes equations for laminar flow or by the Reynolds averaged Navier-Stokes (RANS) equations for turbulent model.
|
|
|
Numerical Solution of 3D Airflow in Channel Representing a Vocal Tract
Pořízková, P. ; Kozel, Karel ; Horáček, Jaromír
This study deals with the numerical solution of a 3D compressible flow of a viscous fluid in a channel for low inlet airfloe velocity. the channel is a simplified model of the glottal space in the human vocal tract. The system of Navier-Stokes equations has been used as mathematical model of laminar flow of the compressible viscous fluid in a domain. The numerical solution is implemented using the finite volume method (FVM)and the predictor-corrector MacCormack scheme with artifical viscosity using a grid of quadrilateral cells.
|
|
|
Numerical solution of compressible subsonic flows in 3D channel
Pořízková, P. ; Kozel, K. ; Horáček, Jaromír
The channels shape is a simplified geometry of the glottal space in the human vocal tract. Goal is numerical simulation of flow in the channels which involves attributes of real flow causing acoustic perturbations. The system of Navier-Stokes equations closed with static pressure expression for ideal gas describes the unsteady laminar flow of compressible viscous fluid. The numerical solution is implemented using the finite volume method and the predictor-corrector MacCormack scheme with artificial viscosity using a grid of quadrilateral cells. The unsteady grid of quadrilateral cells is considered in the form of conservation laws using Arbitrary Lagrangian-Eulerian method.
|
|
|
Analytical and numerical approach to the airfoil stability computation
Chládek, Štěpán ; Horáček, Jaromír
The fluid structure interaction (FSI) represents an important task in many applications. This paper deals with interaction of airfoil and fluid flow and compares two possible ways of computation of the stability boundaries of the system. The airfoil has two degrees of freedom represented by translation h(t) and rotation a(t). The flow acts on the airfoil with aerodynamic forces depending on the flow velocity. Two methods for calculation of the critical flow velocity, when the airfoil loses aeroleastic stability, are described.
|
|
|
Spektrální analýza kontinuálního a syntetizovaného proudu
Broučková, Z. ; Trávníček, Zdeněk ; Šafařík, P.
Tato experimentální práce se zabývá charakterizací tří druhů proudů vzduchu (kontinuálního, syntetizovaného a kontinuálního řízeného syntetizovaným proudem) pomocí spektrální výkonové hustoty (PSD). Tyto charakteristiky byly vyhodnocovány z bodových měření rychlosti, prováděných pomocí anemometru se žhaveným drátkem (CTA). Z naměřeného signálu, popř. z časového průběhu rychlosti, byla spektrální výkonová hustota vyhodnocována pomocí programu MATLAB. Výsledné grafy ukázaly vliv řízení na hlavní proud. U samotného syntetizovaného proudu byla navíc nalezena výrazná kvalitativní změna v okolí hranice jeho existence. Dále bylo prokázáno, že takto zásadní změna proudového pole může být velmi dobře kvantifikována.
|
|
|
On mathematical modelling of gust response using the finite element method
Sváček, P. ; Horáček, Jaromír
In this paper the numerical approximation of aeroelastic response to sudden gust is presented. The fully coupled formulation of two dimensional incompressible viscous fluid flow over a flexibly supported structure is used. The flow is modelled with the system of Navier-Stokes equations written in Arbitrary Lagrangian-Eulerian form and coupled with system of ordinary differential equations describing the airfoil vibrations with two degrees of freedom. The Navier-Stokes equations are spatially discretized by the fully stabilized finite element method. The numerical results are shown.
|
|
|
Parallel numerical simulation of airflow past an oscillating NACA0015 airfoil
Řidký, Václav ; Šidlof, Petr ; Vlček, Václav
This paper focuses on 3D and 2D parallel computation of pressure and velocity fields around an elastically supported airfoil self-oscillating due to interaction with the airflow The results of numerical simulations are compared with data measured in a wind tunnel, where physical model of a NACA0015 airfoil was mounted and tuned to exhibit the flutter instability. The experimental results were obtained previously in the Institute of Thermomechanics by interferometric measurements in a subsonic wind tunnel in Nový Knín. For the numerical solution is implemented in OpenFOAM, an open-source software package based on finite volume method. In the numerical solution is prescribed displacement of the airfoil, which corresponds to the experiment.
|
host ::
RSS.