Numerické řešení obyčejných diferenciálních rovnic pomocí adaptivních BDF metod
Thesis title in Czech: | Numerické řešení obyčejných diferenciálních rovnic pomocí adaptivních BDF metod |
---|---|
Thesis title in English: | Numerical solution of ODEs with the aid of adaptiv BDF methods |
Academic year of topic announcement: | 2009/2010 |
Thesis type: | Bachelor's thesis |
Thesis language: | čeština |
Department: | Department of Numerical Mathematics (32-KNM) |
Supervisor: | prof. RNDr. Vít Dolejší, Ph.D., DSc. |
Author: | hidden - assigned and confirmed by the Study Dept. |
Date of registration: | 10.11.2009 |
Date of assignment: | 10.11.2009 |
Date and time of defence: | 25.06.2010 00:00 |
Date of electronic submission: | 25.06.2010 |
Date of proceeded defence: | 25.06.2010 |
Opponents: | doc. RNDr. Josef Kofroň, CSc. |
Guidelines |
1) nastudování příslušné literatury
2) návrh odhadu lokální chyby disretizace a volby časového kroku 3) implementace a otestování metody pro jednu obyčejnou diferenciální rovnici 4) imlementace metody pro soustavu rovnic v balíku ADGFEM Pozn. bod 4) je víceméně nepovinný, případně by mohl být základem následující DP. |
References |
1) V. Dolejší, P. Kůs, Adaptive backward difference formula - discontinuous Galerkin finite element method for the solution of conservation laws, Int. J. Numer. Meth. Engng. 73(12): 1739-1766, 2008
2) E. Hairer, S. P. Norsett, G. Wanner: Solving ordinary differential equations I, Nonstiff problems, Springer Verlag, 2000. 3) E. Hairer, G. Wanner: Solving ordinary differential equations II, Stiff and differential-algebraic problems, Springer Verlag, 2002. |
Preliminary scope of work |
Efektivní numerické řešení velkých soustav obyčejných diferenciálních rovnic (ODR) představuje stále otevřenou problematiku. Pro tzv. "stiff" úlohy je vhodné použít implicitní metody (např. BDF (backward difference formula) metody), které mají velkou oblast stability. Velmi důležitým aspektem je volba časového kroku, kde příliš velký časový krok může vést ke ztrátě přesnosti a na druhou stranu příliš krátký časový krok vede ke ztrátě efektivity. Optimální volba časového kroku se provádí pomocí adaptivních metod, které odhadují lokální chybu diskretizace a navrhují nový časový krok. V DP P. Kůsa (obhájená v roce 2006) byla navržena metoda, která odhaduje lokální chybu diskretizace pomocí dvou různých metod stejného řádu přesnosti. Efektivnější přístup, který je náplní této BP, je použití pouze jedné metody, kdy k odhadu chyby poslouží diferenční formule pro aproximaci derivací řešení úlohy. |