Stokessche Stromfunktion
Die Stokes’sche Stromfunktion (Formelzeichen , Dimension L³ T−1) von George Gabriel Stokes ist in der Strömungsmechanik ein analytisches Hilfsmittel zur Lösung der Euler-Gleichungen in drei dimensionalen, axialsymmetrischen, stationären Strömungen inkompressibler, reibungsfreier Fluide. Die Stokes’sche Stromfunktion ist also die Anwendung des Konzepts der Stromfunktion auf axialsymmetrische Strömungen, die dann auch analoge Eigenschaften besitzen. Aus Ableitungen der Stokes’schen Stromfunktion ergibt sich das Geschwindigkeitsfeld, das automatisch divergenzfrei und die Strömung mithin volumenerhaltend und dichtebeständig ist. Die Höhenlinien der Stokes’schen Stromfunktion stellen wie im ebenen Fall Stromlinien dar, die hier wegen der Axialsymmetrie Stromröhren beranden. Wie im ebenen Fall ist der Volumenstrom zwischen zwei Stromlinien – im von ihnen berandeten Stromröhrenring – überall gleich.
Definition
Betrachtet wird eine dichtebeständige und stationäre Strömung mit einem ortsabhängigen aber nicht zeitabhängigen weil stationärem Geschwindigkeitsfeld Der Ortsvektor kann bei axialsymmetrischer Strömung vorteilhaft mit Zylinder- oder Kugelkoordinaten parametrisiert werden.
Stokes’sche Stromfunktion in Zylinderkoordinaten
![](../I/Cylindrical_with_grid.svg.png.webp)
Das Zylinderkoordinatensystem wie im Bild wird so ausgerichtet, dass die -Richtung die Richtung ist, um die die Strömung axialsymmetrisch ist. Den Abstand eines Punktes von -Achse gibt die Koordinate an, die hier mit einem großen bezeichnet wird, um eine Verwechselung mit der Dichte zu vermeiden. Der Winkel zählt in Umfangsrichtung senkrecht zur -Achse. Die Geschwindigkeit darf nicht von abhängen und auch keine Komponente in tangentialer -Richtung besitzen. Die Geschwindigkeiten in - und -Richtung ergeben sich dann durch folgende Ableitungen der Stromfunktion :
Der Operator rot bildet die Rotation, grad den Gradient und das Kreuzprodukt.
Stokes’sche Stromfunktion in Kugelkoordinaten
![](../I/Spherical_with_grid.svg.png.webp)
In Kugelkoordinaten ist die Achse mit die Richtung, um die die Strömung axialsymmetrisch ist. Den Abstand eines Punktes vom Ursprung gibt der Radius an und der Winkel zählt – wie in Zylinderkoordinaten – in Umfangsrichtung senkrecht zur -Achse. Wiederum darf die Geschwindigkeit nicht von abhängen und auch keine Komponente in tangentialer -Richtung besitzen. Die Geschwindigkeiten in - und -Richtung berechnen sich dann durch folgende Ableitungen der Stromfunktion :
Der Zusammenhang mit den Zylinderkoordinaten ist durch und bzw. und gegeben.
Eigenschaften von mit Stokes’schen Stromfunktionen beschriebenen Strömungen
Stromlinien
Der Gradient der Stromfunktion in Zylinderkoordinaten ist wegen
senkrecht zur Geschwindigkeit und in Kugelkoordinaten gilt dasselbe:
Die Geschwindigkeit ist per definitionem überall tangential zur Stromlinie, auf der der Wert der Stromfunktion also konstant ist. In der hier vorausgesetzten Axialsymmetrie repräsentiert die Stromlinie eine Stromröhre.
Dichtebeständigkeit
Wenn das Geschwindigkeitsfeld einer axialsymmetrischen Strömung durch eine Stoke’sche Stromfunktion gegeben ist, dann gilt in Zylinderkoordinaten
wie in Kugelkoordinaten
weil die Divergenz von Rotationsfeldern immer null ist. In einer divergenzfreien Strömung verschwindet auf Grund der Massenbilanz überall die substantielle Zeitableitung der Dichte, die daher zeitlich konstant ist.
Eine divergenzfreie Strömung enthält weder Quellen noch Senken, so dass unter den gegebenen Voraussetzungen Stromlinien im Inneren der Flüssigkeit weder beginnen noch enden können. Die Stromlinien sind also entweder torusförmig geschlossen, sind buchstäblich unendlich oder enden auf dem Rand des Strömungsgebiets.
Rotation der Strömung
Die Rotation des Geschwindigkeitsfeldes ist die Wirbelstärke, die in Zylinderkoordinaten wegen
nur eine Komponente ω in tangentialer Umfangsrichtung hat, weswegen die Wirbelstärke als Skalarfeld behandelt werden kann. In Kugelkoordinaten ist das auch so:
Anders als in ebenen Strömungen steht hier auf der rechten Seite des Gleichheitszeichens nicht der Laplace-Operator.
Volumenstrom zwischen Stromlinien
![](../I/Stokesvolstrom.png.webp)
Der Volumenstrom zwischen zwei Stromlinien ist überall gleich. Dies wird anhand zweier Stromlinien gezeigt, auf denen die Stromfunktion die Werte bzw. annimmt, siehe Bild. Um den Volumenstrom zu berechnen, der zwischen diesen beiden Stromlinien hindurchtritt, wird eine Linie mit der Bogenlänge und definiert, die also auf der einen Stromlinie beginnt und auf der anderen Stromlinie endet. Die Parametrisierung mit der Bogenlänge bewirkt, dass die Länge der Linie ist und der Tangenteneinheitsvektor gleich der Ableitung des Ortsvektors ist. Auf Grund der Axialsymmetrie definiert diese Linie eine Fläche , auf der der übertretende Volumenstrom zu bestimmen ist. Der Volumenstrom , der über diese Fläche tritt, berechnet sich mit einem Kurvenintegral und der Normale an die Kurve zu
Indem ersetzt wird, ergibt sich dasselbe Ergebnis in Kugelkoordinaten. Daher gilt hier dasselbe wie bei der Stromfunktion in der Ebene: Unabhängig vom speziellen Kurvenverlauf ist der Volumenstrom zwischen zwei Stromlinien überall gleich. Wenn die Linie auf derselben Stromlinie startet und endet, dann verschwindet der über sie hinweglaufende Volumenstrom. Wenn die gewählte Linie ein Stück einer Stromlinie ist, dann zeigt sich, dass an keiner Stelle einer Stromlinie Fluid über sie hinwegströmt. Eine Stromlinie wirkt auch hier wie eine undurchdringliche Wand.
Bestimmungsgleichungen für die Stromfunktion
Nicht jede Stromfunktion repräsentiert eine physikalisch realistische Strömung. Damit die Stromfunktion im Einklang mit den physikalischen Gesetzen ist, muss sie den eulerschen Gleichungen gehorchen, aus denen sich – wie sich zeigt – die Stromfunktion unabhängig vom Druck berechnen lässt. In einem konservativen Schwerefeld gestaltet sich die Suche nach der Stromfunktion besonders einfach. Anders als im ebenen Fall ergeben sich Differentialgleichungen mit variablen Koeffizienten, was die Lösung erschwert.
Eulersche Gleichungen
Die Euler-Gleichungen liefern nach Bildung der Rotation
Bestimmungsgleichungen für die Stromfunktion aus der Tabelle:
Koordinatensystem | Bestimmungsgleichung |
---|---|
Zylinderkoordinaten | |
Kugelkoordinaten |
Darin ist die -Komponente der Wirbelstärke, siehe oben. Diese Gleichungen muss die Stromfunktion erfüllen, damit sie eine physikalisch realistische Strömung beschreibt.
Beweis |
Ausnutzung der Grassmann-Entwicklung zeigt bei der Bildung der Rotation in den Euler-Gleichungen: denn Gradientenfelder sind immer rotationsfrei. Mit der Wirbelstärke : und der Stromfunktion ergibt sich in
|
Konservatives Beschleunigungsfeld
In einem konservativen Beschleunigungsfeld – wie es die Schwerkraft eines ist – verschwinden die rechten Seiten der Bestimmungsgleichungen wegen der Rotationsfreiheit solcher Felder. Dann kann – wie im ebenen Fall – argumentiert werden: die im obigen Beweis als Zwischenergebnis angefallene Bestimmungsgleichung
wird mit
und einer beliebigen Funktion f immer erfüllt:
In Kugelkoordinaten gilt Analoges mit dem Endergebnis:
Zylinderkoordinaten | |
---|---|
Kugelkoordinaten |
Insbesondere ist erlaubt.
Randbedingungen
Ein Strömungsfeld kann nur bei festen Wänden stationär sein. Die Randbedingungen werden entlang von Meridiankurven vorgegeben, die mit und der Bogenlänge s definiert werden. Dann lautet der Tangenteneinheitsvektor und die Normale der Kurve in radialer Richtung . Fließt nirgends Fluid über die Linie, dann ist sie ein Teil einer Stromlinie und die Linie stellt gleichzeitig eine Wand dar.
Die Dirichlet-Randbedingungen geben den Wert der Stromfunktion entlang einer solchen Linie vor und in Zylinderkoordinaten folgt:
weswegen mit Dirichlet-Randbedingungen die radiale Geschwindigkeit senkrecht zu Linien festgelegt wird. Ist der Wert der Stromfunktion auf der Linie konstant, dann ist die Linie ein Teil einer Stromlinie und die Normalkomponente der Geschwindigkeit verschwindet entlang der Linie.
Die Neumann-Randbedingungen geben die Ableitungen der Stromfunktion senkrecht zu Linien vor:
Durch die Neumann-Randbedingungen wird also die Geschwindigkeitskomponente tangential zur Linie vorgegeben. In Kugelkoordinaten ergibt sich Gleiches mit
Beispiel
In Zylinderkoordinaten gilt in einer ebenen Strömung und die Geschwindigkeit hat nur eine radiale Komponente. Mit der obigen Bestimmungsgleichung ergibt sich mit dann:
Also verschwindet die zweite Ableitung der Stokes’schen Stromfunktion nach der z-Koordinate und die erste Ableitung ist mithin eine Konstante . Dann lautet die Radialgeschwindigkeit
was die Geschwindigkeitsverteilung der ebenen Quelle/Senke ist.
Eine in Kugelkoordinaten nur vom Radius abhängige Strömung ist die drei-dimensionale Quelle/Senke. Mit und ergibt sich aus der obigen Bestimmungsgleichung
woraus das Geschwindigkeitsfeld einer drei dimensionalen Quelle/Senke folgt:
Hier nimmt die Geschwindigkeit also mit dem Quadrat des Abstands zur Quelle ab.
Literatur
- G. K. Batchelor: An Introduction to Fluid Dynamics. Cambridge University Press, 1967, ISBN 0-521-66396-2.
- P. K. Kundu: Fluid Mechanics. Academic Press, 2015, ISBN 978-0-12-405935-1.