imported>Bithisarea: /* growthexperiments-addlink-summary-summary:2|0|0 */

2025-01-25T22:30:55Z

growthexperiments-addlink-summary-summary:2|0|0

Neue Seite

In der [[Numerische Mathematik|numerischen Mathematik]] ist die '''Spektralmethode''' ein Verfahren zur Lösung von [[Partielle Differentialgleichung|partiellen Differentialgleichungen]], wie den [[Gleichungen von Navier-Stokes|Navier-Stokes Gleichungen]], mittels globaler [[Ansatzfunktion]]en. Ansatzfunktionen können z. B. [[Fourierreihe]]n oder [[Tschebyschow-Polynom|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 dient die Bezeichnung Spektralmethode als übergeordneter Begriff. Besonders bieten sich hier orthogonale Funktionen an, z. B. [[Orthogonale Polynome]]. Denn diese haben die Eigenschaft sich nicht gegeneinander zu beeinflussen. Dadurch ist bei jedem neuen Polynom sichergestellt, sich der Lösung besser annähern zu können, ohne das erreichte zu verschlechtern.
Als effizientes Verfahren zur Hin- und Rücktransformation bietet sich die [[Schnelle Fourier-Transformation|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 physikalischen 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 [[Gibbssches Phänomen|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 [[Lineares Gleichungssystem|linearen Gleichungssystemen]] mit vollbesetzten und unsymmetrischen [[Matrix (Mathematik)|Matrizen]] führen kann. Zur Lösung sind dann [[Iteratives Verfahren|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.

== Verbindung zur Spektral Elemente Methode ==
Die Spektralmethode weist besonders gute Konvergenzeigenschaften auf, je glatter die Lösung ist. Um dennoch gute Konvergenzen, aber auch gleichzeitig gute geometrische Abbildungen zu bekommen, wurde die ''Spektral Elemente Methode'' eingeführt. Diese nutzt die Herangehensweise der Spektralmethode für einzelne Elemente und verknüpft diese dann. Das Ergebnis ist damit mit einer ''Finite Elemente Methode'' (FEM), von sehr hoher Ordnung vergleichbar.<ref>{{Literatur |Autor=David Kopriva |Titel=Implementing Spectral Methods for Partial Differential Equations: Algorithms for Scientists and Engineers |Verlag=Springer Netherlands |Datum=2009 |Reihe=Scientific Computation |ISBN=978-90-481-2260-8 |Online=https://www.springer.com/de/book/9789048122608 |Abruf=2020-10-10}}</ref>

== Literatur ==
* Claudio Canuto: ''Spectral methods in fluid dynamics.'' Springer, 1993, ISBN 3-540-52205-0.
* Wilhelm Heinrichs: ''Efficient iterative solution of spectral systems for the Navier-Stokes equations.'' Köster, 1994, ISBN 3-929937-81-6. (Zugl.: Düsseldorf, Univ., Habil.-Schr., 1992) (Wissenschaftliche Schriftenreihe Mathematik Bd. 2)
* Javier de Frutos, Julia Novo: ''A Spectral Element Method for the Navier--Stokes Equations with Improved Accuracy.''] In: ''SIAM J. Numer. Anal.'' 38(3), S. 799–819.

== Weblinks ==
* [http://channelflow.org/ Channelflow] ist ein unter GPL gestelltes Programm basierend auf dem Algorithmus aus Kapitel 7.3 des Buches "Spectral Methods in Fluid Dynamics"

== Einzelnachweise ==
<references />
[[Kategorie:Numerische Mathematik]]

Spektralmethode - Versionsgeschichte

imported>Bithisarea: /* growthexperiments-addlink-summary-summary:2|0|0 */