No label defined (Q4380358)

Aus dem Text: In der vorliegenden Arbeit wird eine neue a posteriori adaptive Lösungsstrategie für das konvektionsdominante Modellproblem \(\varepsilon\Delta u+\vec p\cdot\nabla u+ qu=f\) in \(\Omega\subset\mathbb{R}^2\), \(u|_{\partial\Omega}= 0\) vorgestellt. Die Gittergenerierung beruht auf speziellen Informationen einer Wavelet-Dekomposition der diskreten Lösung. Erstmalig wurde es durch den Einsatz dieser Technik möglich, ausgehend von einer uniformen Triangulierung, anisotrope Gitter zu erzeugen. Zu diesem Zweck konstruieren wir Spline-Waveletbasen, die den homogenen Randbedingungen unseres Modellproblems angepaßt sind. Hierbei gehen wir vom eindimensionalen Fall aus, übertragen diese Basen mit einem Tensorproduktansatz auf den zweidimensionalen Fall und konstruieren spezielle Elemente auf Friedrichs-Keller-Triangulierungen. Dabei werden ausführlich die Eigenschaften dieser Basisfunktionen, insbesondere die Approximationseigenschaften, untersucht. Die numerische Lösung erfolgt mit einem Mittelpunkt-upwind-Verfahren, das ein gleichmäßig stabiles Verfahren bzgl. des Parameters \(\varepsilon\) ist. Im Gegensatz zum einfachen upwind-Verfahren besitzt das Mittelpunkt-upwind-Verfahren eine höhere formale Konsistenzordnung. Wir untersuchen verschiedene Methoden adaptiver (a posteriori) Strategien und versuchen, die hier angewandte adaptive Methode einzuordnen. Es wurden zwei MATLAB-Programme erstellt, die das Modellproblem auf dem Einheitsintervall bzw. Einheitsquadrat mit der a posteriori adaptiven Strategie numerisch lösen. Eine Kurzdokumentation dieser Programme befindet sich im Anhang dieser Arbeit. Fehlerabschätzungen konnten bisher mit der angegebenen a posteriori Methode nicht bewiesen werden. Aus diesem Grund vergleichen wir unsere numerischen Ergebnisse mit Ergebnissen verschiedener Programmsysteme, so z.B. Femlab und UG3.