No label defined (Q4254902)

Wir untersuchen hauptsächlich dreidimensionale Konvektionszustände bei positiver Soretkopplung. Das Problem wird durch nichtlineare partielle Differentialgleichungen beschrieben. Das vorrangige Ziel dieser Arbeit war es, einen effizienten numerischen Code zu erstellen, zu implementieren und zu testen. Er sollte in der Lage sein, mit einem ausgewogenen Verhältnis aus Aufwand und Ergebnis dreidimensionale Konvektionsdynamik zu simulieren. Während Teil II dieser Arbeit den physikalischen Aspekten vorbehalten bleibt, beschreibt Teil I den Algorithmus zur Integration der dreidimensionalen Oberbeck-Boussinesq-Gleichungen. Grundlage ist ein explizites Finite-Differenzen-Verfahren vom MAC (Marker-And-Cell) Typ. Die Zeitdiskretisierung erfolgt mit einem DuFort-Frankel-Verfahren. Zur Bestimmung des Druckes wurde eine Poisson-Gleichung abgeleitet, für deren Lösung ein schnelles Mehrgitter-Verfahren entwickelt wurde. Es nutzt eine spezielle Eigenschaft von klassischen Iterationsverfahren (z.B. SOR-Verfahren): Abhängig von der Feinheit der Diskretisierung, dämpfen sie hochfrequente Fehleranteile recht stark, die glatten Beiträge werden jedoch nur langsam vermindert. Nach Durchführung von ein bis zwei Schritten eines klassischen Iterationsverfahrens, nun als Glättungsschritte bzw. Glättungen bezeichnet, ist die Abweichung von der exakten Lösung, der Defekt, nicht klein, aber glatt. Es ist deshalb möglich, ihn auf einem gröberen Gitter darzustellen, indem man etwa im 1D-Fall jeden zweiten Gitterpunkt wegfallen lässt. Diesen Prozeß nennt man Restriktion. Die exakte Lösung der niedrigdimensionalen Defektgleichung auf dem gröberen Gitter wird anschließend wieder geeignet auf das feine Gitter interpoliert; im 1D-Fall kann man z.B. durch je zwei Werte Geraden berechnen und die entsprechenden Feingitterwerte ergänzen. Dieser Prozeß heißt Prolongation. Auf die Darstellung des MAC-Mehrgitter-Verfahrens wurde in Teil I besonders viel Wert gelegt, denn ähnliche Versuche ein Mehrgitter-Verfahren in den MAC-Code einzubinden, sind nur spärlich dokumentiert. Die Teil I und II können unabhängig voneinander gelesen werden.