| |
|
A Brief Design of Optical Resonators
Hubík, Daniel ; Nešpor, Dušan (referee) ; Kadlec, Radim (advisor)
This bachelor thesis is focused on analysis of split-ring resonators in THz region. Simulations were made by finite elements method and by finite-difference time-domain method. At first we created a resonating structure that works in GHz region. Then we were observing a dependence of movement of resonant frequency on the size of resonator. In the final chapter we assigned frequency dependent values of permitivity to such structure. As the result we simulated working resonator at frequency 500THz. All simulations have been made in program HFSS ANSYS and Lumerical FDTD Solutions.
|
|
The Selected Stochastic Programs in Engineering Design
Čajánek, Michal ; Mrázková, Eva (referee) ; Popela, Pavel (advisor)
Two-stage stochastic programming problem with PDE constraint, specially elliptic equation is formulated. The computational scheme is proposed, whereas the emphasis is put on approximation techniques. We introduce method of approximation of random variables of stochastic problem and utilize suitable numerical methods, finite difference method first, then finite element method. There is also formulated a mathematical programming problem describing a membrane deflection with random load. It is followed by determination of the acceptableness of using stochastic optimization rather than deterministic problem and assess the quality of approximations based on Monte Carlo simulation method and the theory of interval estimates. The resulting mathematical models are implemented and solved in the general algebraic modeling system GAMS. Graphical and numerical results are presented.
|
|
Analysis of Methods of Differences for Partial Differential Equations Solving
Zpěváková, Jana ; Zbořil, František (referee) ; Šátek, Václav (advisor)
In this thesis, we discuss the numerical solution of ordinary differential equation and numerical methods of solving partial differential equations. We propose and implement an application, that converts partial differential hyperbolic equation to a set of ordinary differential equations using finite difference method. After that, the system of equations is solved using the Taylor method programmed in Matlab environment. Finally, we compare the time complexity of proposed solution with parallel numerical computation.
|
|
Simulation of the Heat Diffusion with a Time-Varying Source on GPUs
Hála, Pavel ; Záň, Drahoslav (referee) ; Jaroš, Jiří (advisor)
This bachelor's thesis deals with the simulation of the heat transfer inside human tissue injected by an external time varying heat source. The proposed implemented simulation is based on a 4th order in space and 1st order in time finite-difference time domain method. First, a multithreaded CPU version was implemented. Subsequently, several GPU accelerated versions were implemented taking into account architecture aspect of the GPU. The experimental results showed that the fastest GPU kernel was the naive one using only the GPU global memory. Next, the usefulness of the Gauss-Seidel's method was investigated. The CPU implementation of the method was evaluated as usable because of being only 13% slower while saving up to 50% of memory resources. However, the GPU implementation was twice as slow as the naive version mainly due to shared memory size limits. The peak performance in terms of GFLOPS reached 32 and 135 on CPU and GPU, respectively. This corresponds to 10% and 9% of the theoretical potential of given architectures.
|
|
Rovnice vedení tepla ve fyzice planetek a meteoroidů
Pohl, Leoš ; Brož, Miroslav (advisor) ; Vokrouhlický, David (referee)
Non-gravitational forces caused by thermal emission of photons can significantly change orbits and spin states of asteroids in the long term. A solution of the Heat Conduction Equation (HCE) in an asteroid is necessary to evaluate the forces. Finite Difference Methods (FDMs) are implemented in a Fortran numerical HCE solver to calculate a temperature distribution within a system of 1-dimensional slabs which approximate the asteroid. We compare the methods w.r.t. convergence, accuracy and computational efficiency. The numerical results are compared with a simplified steady-state analytical solution. We calculate the non-gravitational accelerations and resulting semimajor axis drift from the numerical results. The implemented FDMs are shown to be convergent with denser grids and the best method has been selected. The analytical solution provides a good first-guess estimate of the temperature amplitude. The drift in semimajor axis of the tested asteroids, which is due to the non-gravitational forces, is in order-of-magnitude agreement with more accurate models and observational data.
|
|
Design optimization of packed bed for thermal energy storage
Krist, Thomas ; Charvát, Pavel (referee) ; Klimeš, Lubomír (advisor)
Tato diplomová práce se zabývá tématem výměny tepla v zásobníku tepla typu ”packed bed”. Cílem je popsat přenos tepla v zásobníku tepla obsahující kamínky malých průměrů, skrz který proudí horký vzduch. Toto je modelováno v prostředí MATLAB. Na začátku je krátký úvod do problematiky zahrnující ukládání tepla a jeho možné využití. Dále je uveden krátký přehled o základech přenosu tepla, typech přenosu tepla a termofyzikální vlastnosti systému vzduch-kámen. Ve třetí kapitole je představen zásobník tepla typu ”packed bed” a rozličné modely a dané podmínky jsou vysvětleny. Další kapitola se zabývá s numerickými metodami, převážně s metodou konečných diferencí použitou v této práci. Pátá kapitola se zaměřuje na obecnou optimalizaci daného problému přenosu tepla. Populačně založený metaheuristický optimalizační algoritmus zvaný Genetický algoritmus je popsán. Sestavení modelu je ukázáno v šesté kapitole, stejně jako prezentace výsledků získaných z programu MATLAB. V poslední kapitole je pak diskutován závěr a doporučení.
|
| |
|
Modification of Navier_Stokes equations asuming the quasi-potential flow
Navrátil, Dušan ; Pochylý, František (referee) ; Fialová, Simona (advisor)
The master's thesis deals with Navier-Stokes equations in curvilinear coordinates and their solution for quasi-potential flow. The emphasis is on detailed description of curvilinear space and its expression using Bézier curves, Bézier surfaces and Bézier bodies. Further, fundamental concepts of hydromechanics are defined, including potential and quasi-potential flow. Cauchy equations are derived as a result of the law of momentum conservation and continuity equation is derived as a result of principle of mass conservation. Navier-Stokes equations are then derived as a special case of Cauchy equations using Cauchy stress tensor of Newtonian compressible fluid. Further transformation into curvilinear coordinates is accomplished through differential operators in curvilinear coordinates and by using curvature vector of space curve. In the last section we use results from previous chapters to solve boundary value problem of quasi-potential flow, which was solved by finite difference method using Matlab environment.
|
|
Heat transfer solution of solidifying steel system with phase change with moving edge conditions
Fedorko, Tomáš ; Mauder, Tomáš (referee) ; Štětina, Josef (advisor)
Cílem diplomové práce je vytvoření 2D numerického modelu pohybujícího se řezu s proměnnými okrajovými podmínkami skutečné geometrie plynulého odlévání a chlazení předlitku v prostředí MATLAB. Model se zabývá vysoce nelineárními termofyzikálními podmínkami oceli během tuhnutí a chlazení. V práci je simulovaná nejen nelinearita termofyzikálních podmínek, ale také nelinearita při fázové změně. Fázová změna je modelovaná pomocí metody entalpie, metody zdánlivé kapacity a metody teplotního zotavení. Všechny výsledky práce jsou porovnány z více hledisek, jako např. z hlediska přesnosti, rychlosti výpočtu, nebo vhodnosti časového diskretizačního kroku pro nelineární problémy, a paralelizace.
|