Počítačová simulace a numerická analýza problémů stlačitelného proudění

Abstract

The thesis deals with the construction of an adaptive 1D and 2D mesh in the framework of the cell- centered finite volume scheme. The adaptive strategy is applied to the numerical solution of problems governed by the Euler equations, which is a hyperbolic system of PDE's. The used algorithm is applicable to nonstationary problems and consists of three independent parts, which are cyclically repeated. These steps are PDE evolution, then mesh adaptation and recovery of numerical solution from the old mesh to the newly adapted mesh. Owing to this the algorithm can be used also for other hyperbolic systems. The thesis is focused on the development of our mesh adaptation strategy, based on the anisotropic mesh adaptation, which preserves the geometric mass conservation law in each computational step. The proposed method is suitable to solve problems with moving discontinuities. Several test problems with moving discontinuity are computed to compare our algorithm with Moving Mesh algorithms.

Tato práce se zabývá konstrukcí adaptivní výpočtové sítě v 1D a ve 2D v kontextu metody konečných objemů. Adaptivní strategie je aplikována na numerické řešení Eulerových rovnic, což je hyperbolický systém parciálních diferenciálních rovnic. Použitý postup je určen pro nestacionární problémy a skládá se ze tří v podstatě nezávislých kroků, jež jsou cyklicky opakovány. Těmito kroky jsou: výpočet pomocí schématu metody konečných objemů, dále pak adaptace sítě a přepočet numerického řešení z neadaptované sítě na síť adaptovanou. Díky tomu je tento algoritmus použitelný i na jiné, nejen hyperbolické systémy. Těžiště práce spočívá v návrhu vlastní adaptační strategie, založené na anisotropní adaptivitě, která bude v každém adaptačním kroku splňovat tzv. geometrický zákon zachování. V práci je též porovnání námi navržené strategie s algoritmy typu Moving Mesh pro úlohy s pohybující se nespojitostí.