|
|
A posteriori error estimates for numerical solution of convection-difusion 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 k-omega turbulence model. The method is intended for a subsonic flow.
|
|
|
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
|
|
|
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 Navier-Stokes 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.
|
|
|
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 k-omega turbulence model. The method is intended for a subsonic flow.
|
|
|
A posteriori error estimates for numerical solution of convection-difusion 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
|
|
|
Discontinuous Galerkin method for the solution of compressible viscous flow
Česenek, Jan ; Feistauer, Miloslav (advisor) ; Najzar, Karel (referee) ; Sváček, Petr (referee)
Title: Discontinuous Galerkin method for solving compressible viscous flow Author: Jan Česenek Department: Department of Numerical Mathematics Supervisor: prof. RNDr. Miloslav Feistauer, DrSc., dr.h.c., Department of Numerical Mathematics Abstract: The subject of this PhD thesis is the numerical simulation of the interaction of two-dimensional compressible viscous flow and a vibrating airfoil. We consider a solid airfoil with two degrees of freedom which can rotate around the elastic axis and oscillate in the vertical direction. The numerical simulation of this problem consist of the dis- continuous Galerkin finite element method solving Navier-Stokes equations coupled with a system of nonlinear ordinary differential equations describing the airfoil motion. The time-dependent domain is taken into account with the aid of the Arbitrary Lagrangian- Eulerian(ALE) formulation. Theoretical part of this paper is concerned with error esti- mates of the space-time discontinuous Galerkin method for scalar nonstationary equations with nonlinear convection and nonlinear diffusion. Keywords: convection-diffusion problems, discontinuous Galerkin method, interaction of a fluid with a vibrating airfoil, ALE method
|
|
|
Numerical Solution of a Fredholm Integral Equation of the Second Kind Related to Induction Heating
Rak, Josef ; Kofroň, Josef (advisor) ; Feistauer, Miloslav (referee) ; Sváček, Petr (referee)
This thesis deals with numerical solution of an integral equation of the second kind with special singular kernel function related to induction heating. The numerical solution is based on collocation and Nyström methods. The idea of collocation methods is to choose a finite-dimensional space of candidate solutions (usually polynomials up to a certain degree). The Nyström methods are based on approximation of the integral in equation by numerical integration rule. This thesis describes and gives error estimates of both methods. Error estimates are compared to the exact solutions in simple cases.
|
guest ::
RSS feed.