
Finite element approximation of fluid structure interaction using TaylorHood and ScottVogelius elements
Vacek, Karel ; Sváček, P.
This paper addresses the problem of fluid flow interacting a vibrating solid cylinder described by one degree of freedom system and with fixed airfoil. The problem is described by the incompressible NavierStokes equations written in the arbitrary EulerianLagrangian (ALE) formulation. The ALE mapping is constructed with the use of a pseudoelastic approach. The flow problem is numerically approximated by the finite element method (FEM). For discretization of the fluid flow, the results obtained by both the TaylorHood (TH) element and the ScottVogelius (SV) finite element are compared. The TH element satisfies the BabuškaBrezzi infsup condition, which guarantees the stability of the scheme. In the case of the SV element the mesh, that is created as a barycentric refinement of regular triangulation, is used to satisfy the BabuškaBrezzi condition. The numerical results for two benchmark problems are shown.


Interpolation with restrictions  role of the boundary conditions and individual restrictions
Valášek, Jan ; Sváček, P.
The contribution deals with the remeshing procedure between two computational finite element meshes. The remeshing represented by the interpolation of an approximate solution onto a new mesh is needed in many applications like e.g. in aeroacoustics, here we are particularly interested in the numerical flow simulation of a gradual channel collapse connected with a~severe deterioration of the computational mesh quality. Since the classical Lagrangian projection from one mesh to another is a dissipative method not respecting conservation laws, a conservative interpolation method introducing constraints is described. The constraints have form of Lagrange multipliers enforcing conservation of desired flow quantities, like e.g. total fluid mass, flow kinetic energy or flow potential energy. Then the interpolation problem turns into an error minimization problem, such that the resulting quantities of proposed interpolation satisfy these physical properties while staying as close as possible to the results of Lagrangian interpolation in the L2 norm. The proposed interpolation scheme does not impose any restrictions on mesh generation process and it has a relatively low computational cost. The implementation details are discussed and test cases are shown.


A posteriori error estimates for numerical solution of convectiondifusion problems
Šebestová, Ivana ; Dolejší, Vít (advisor) ; Sváček, Petr (referee) ; Brandts, Jan (referee)
This thesis is concerned with several issues of a posteriori error estimates for linear problems. In its first part error estimates for the heat conduction equation discretized by the backward Euler method in time and discontinuous Galerkin method in space are derived. In the second part guaranteed and locally efficient error estimates involving algebraic error for Poisson equation discretized by the discontinuous Galerkin method are derived. The technique is based on the flux reconstruction where meshes with hanging nodes and variable polynomial degree are allowed. An adaptive strategy combining both adaptive mesh refinement and stopping criteria for iterative algebraic solvers is proposed. In the last part a numerical method for computing guaranteed lower and upper bounds of principal eigenvalues of symmetric linear elliptic differential operators is presented. 1

 

Solution of inverse problem for a flow around an airfoil
Šimák, Jan ; Feistauer, Miloslav (advisor) ; Felcman, Jiří (referee) ; Sváček, Petr (referee)
Title: Solution of inverse problem for a flow around an airfoil Author: Mgr. Jan Šimák Department: Department of Numerical Mathematics Supervisor: prof. RNDr. Miloslav Feistauer, DrSc., dr. h. c., Department of Numerical Mathematics Abstract: The method described in this thesis deals with a solution of an inverse problem for a flow around an airfoil. It can be used to design an airfoil shape according to a specified velocity or pressure distribution along the chord line. The method is based on searching for a fixed point of an operator, which combines an approximate inverse and direct operator. The approximate inverse operator, derived on the basis of the thin airfoil theory, assigns a corresponding shape to the specified distribution. The resulting shape is then constructed using the mean camber line and thickness function. The direct operator determines the pressure or velocity distribution on the airfoil surface. We can apply a fast, simplified model of potential flow solved using the Fredholm integral equation, or a slower but more accurate model of RANS equations with a komega turbulence model. The method is intended for a subsonic flow.


Aerodynamic transfer of energy to vibrating vocal folds for different driving mechanisms
Valášek, J. ; Sváček, P. ; Horáček, Jaromír
This paper studies the mutual energy transfer between the fluid flow, described by incompressible NavierStokes equations, and the elastic body represented by vocal folds. The aerodynamic energy transfer function describes the amount and more importantly the sign of the energy exchange. It determines if the vocal fold vibrations are selfexcited or prescribed.The energy transfer function is studied for three different driving mechanisms introduced by different inlet boundary conditions (BC). The most frequently used inlet BCs for incompressible model of fluid flow approximated by the finite element method are either Dirichlet BC giving the inlet velocity or donothing type of BC prescribing the pressure difference between the inlet and the outlet. Since the numerical simulations with both aforementioned BCs do not provide results observed experimentally the newly introduced BC based on the penalization approach seems as remedy. The numerical model consists of strongly coupled partitioned scheme based on the stabilized finite element method.\n\n


The Gibbs phenomenon in the discontinuous Galerkin method
Stará, Lenka ; Kučera, Václav (advisor) ; Sváček, Petr (referee)
The solution of the Burgers' equation computed by the standard finite element method is degraded by oscillations, which are the manifestation of the Gibbs phenomenon. In this work we study the following numerical me thods: Discontinuous Galerkin method, stable low order schemes and the flux corrected technique method in order to prevent the undesired Gibbs phenomenon. The focus is on the reduction of severe overshoots and under shoots and the preservation of the smoothness of the solution. We consider a simple 1D problem on the interval (0, 1) with different initial conditions to demonstrate the properties of the presented methods. The numerical results of individual methods are provided. 1


The Influence of Different Geometries of Human Vocal Tract Model on Resonant Frequencies
Valášek, Jan ; Sváček, Petr ; Horáček, Jaromír
This paper presents the transfer function approach in order to determine the acoustic resonant frequencies of a human vocal tract model. The transfer function is introduced here as an acoustic pressure ratio between input amplitude at glottis position and output amplitude at mouth opening given by the solution of Helmholtz equation. This equation is numerically approximated by finite element method. The influence of different boundary conditions are studied and also different locations of excitation and microphone. Four different vocal tract geometries motivated by vocal tract geometry for vowel [u:] are investigated. Its acoustic resonance frequencies in range 100  2500 Hz are computed and compared with published results. Further, the transient acoustic computation with different acoustic analogies are performed. The frequency spectra of Lighthill analogy, acoustic wave equation and perturbed convective wave equation are compared, where the vocal tract model with best frequency agreement with published results was chosen. The dominant frequencies correspond with predicted frequencies of transfer function approach.\n


On Finite Element Approximation of Flow Induced Vibration of Elastic Structure
Valášek, J. ; Sváček, P. ; Horáček, Jaromír
In this paper the fluidstructure interaction problem is studied on a simplified model of the human vocal fold. The problem is mathematically described and the arbitrary LagrangianEulerian method is applied in order to treat the time dependent computational domain. The viscous incompressible fluid flow and linear elasticity models are considered. The fluid flow and the motion of elastic body is approximated with the aid of fininite element method. An attention is paid to the applied stabilization technique. The whole algorithm is implemented in an inhouse developed solver. Numerical results are presented and the influence of different inlet boundary conditions is discused.


Incompressible and compressible viscous flow with low Mach numbers
Balázsová, M. ; Feistauer, M. ; Sváček, Petr ; Horáček, Jaromír
In this paper we compare incompressible flow and low Mach number compressible viscous flow. Incompressible NavierStokes equations were treated with the aid of discontinuous Galerkin method in space and backward difference method in time. We present numerical results for a flow in a channel which represents a simplified model of the human vocal tract. Presented numerical results give a good correspondence between the incompressible flow and the compressible flow with low Mach numbers.
