Spektralmethode

In der numerischen Mathematik ist die Spektralmethode ein Verfahren zur Lösung von partiellen Differentialgleichungen wie den Navier-Stokes Gleichungen mittels globaler Ansatzfunktionen. Ansatzfunktionen können z. B. Fourierreihen oder Tschebyscheff-Polynome sein. Im Laufe eines numerischen Lösungsverfahrens wird die physikalische Repräsentation eines Problems in den Spektralbereich transformiert. Die Unbekannten sind dann nicht mehr physikalische Größen, wie ein diskreter Geschwindigkeits- oder Temperaturverlauf, sondern die Spektralkoeffizienten der globalen Ansatzfunktion. Daher die Bezeichnung Spektralmethode als übergeordneter Begriff. Als effizientes Verfahren zur Hin- und Rücktransformation bietet sich die schnelle Fouriertransformation (FFT) an. Auch Koeffizienten von Tschebyscheff-Polynomen, können hiermit bestimmt werden, sofern als Stützstellen die Gauß-Lobatto-Punkte verwendet werden, da sich dann die Transformation auf die Realteile einer Fouriertransformation beschränkt.

Günstige Konvergenzeigenschaften zeigen diese Verfahren bei Aufgaben, deren Lösungen ein hohes Maß an Glattheit besitzen. Darüber hinaus sollte die Ansatzfunktion dem physikalischem Problem angepasst sein. Bei periodischen Randbedingungen bieten sich Fourierreihen an. Bei festen, nicht periodischen Werten an den Rändern des Lösungsgebietes sollten Ansatzfunktionen verwendet werden, die diese Verläufe auch natürlicherweise wiedergeben können. Wenn darüber hinaus auch eine feinere Diskretisierung an den Rändern erforderlich ist, sind hier Tschebyscheff-Polynome von Vorteil (siehe Gauß-Lobatto-Punkte). Werden jedoch stattdessen Fourierreihen eingesetzt, so ist mit gibbsschen Schwingungen zu rechnen. Außerdem wird das äquidistante Gitter einer Fourierapproximation durch eine FFT der feineren Auflösung an den Rändern nicht gerecht.

Ein typischer Fall in der beide Approximationen Verwendung finden ist die dreidimensionale ebene Kanalströmung. Aufgrund der hohen Gradienten in Wandnähe und des eindeutig nicht-periodischen Verhaltens an der Wand, werden in Wandnormalen-Richtung Tschebyscheff-Polynome eingesetzt. In Hauptströmungs- und Spannweitenrichtung sind jedoch periodische Randbedingungen gefordert, um einen unendlich ausgedehnten ebenen Kanal numerisch zu modellieren.

Nachteilig ist, dass die Spektralmethode zu linearen Gleichungssystemen mit vollbesetzten und unsymmetrischen Matrizen führen kann. Zur Lösung sind dann iterative Verfahren erforderlich. Das Mehrgitterverfahren hat sich hier bewährt. Es gibt jedoch Verfahren, bei denen eine geschickte Umsortierung zu einer Matrix mit Bandstruktur führt. Hier ist eine LU-Zerlegung von Vorteil. Techniken zur Gebietszerlegung (engl. domain decomposition) sind ebenfalls von Interesse.

Literatur

• Canuto, Claudio: Spectral methods in fluid dynamics, Springer, 1993 ISBN 3-540-52205-0, 0-387-52205-0
• Heinrichs, Wilhelm: Efficient iterative solution of spectral systems for the Navier-Stokes equations, Köster, 1994 Zugl. : Düsseldorf, Univ., Habil.-Schr., 1992 . - (Wissenschaftliche Schriftenreihe Mathematik ; Bd. 2) . - ISBN 3-929937-81-6
• Javier de Frutos, Julia Novo: A Spectral Element Method for the Navier--Stokes Equations with Improved Accuracy

Weblinks

• Channelflow ist ein unter GPL gestelltes Programm basierend auf dem Algorithmus aus Kapitel 7.3 des Buches "Spectral Methods in Fluid Dynamics"