Kowalewskaja-Kreisel

Kreiseltyp mit analytisch lösbaren Bewegungsgleichungen

Der Kowalewskaja-Kreisel ist in der Kreiseltheorie einer der drei Kreiseltypen, deren Bewegungsgleichungen bei beliebigen Anfangsbedingungen analytisch gelöst werden können. Er ist ein schwerer symmetrischer Kreisel, dessen drei Hauptträgheitsmomente A, B und C die Bedingung A = B = 2C erfüllen, und dessen Massenmittelpunkt in der Ebene liegt, die von den zu A und B gehörenden Hauptachsen aufgespannt wird.

Sofja Wassiljewna Kowalewskaja gab die Entdeckung des Kreisels 1888 in einer preisgekrönten Arbeit über die Rotation eines Starrkörpers um einen Fixpunkt[1] und gleich auch die Lösung der Euler-Poisson-Gleichungen bekannt. Für die Darstellung der Lösungsfunktionen benötigte Kowalewskaja etwa 50 Seiten, was die Komplexität der Bewegungsmöglichkeiten des Kreisels unterstreicht. Die Vereinfachung der Lösung und die Klassifizierung der Bewegungsmodi ist so aufwändig, dass sie die Wissenschaft noch bis ins 21. Jahrhundert beschäftigten[2].

Die Kowalewskaja-Konstante K = k² verschwindet für Staude-Drehungen, siehe die Animationen in Abb. 3 und 4. Das Bifurkationsdiagramm des Kowalewskaja-Kreisels ist durchsetzt mit Separatrizen, die Phasen, wo gleichartige, topologisch verwandte Bewegungstypen auftreten, voneinander trennen, siehe Abb. 5. Eine typische Bewegung des Kowalewskaja-Kreisels weist im körperfesten Raum der Drehimpulse Trajektorien auf, die einen ringförmigen Schlauch umlaufen und ihn dabei dicht ausfüllen, siehe Abb. 6a-d. Die Separatrizen sind die Phasengrenzen, in deren Zentren eine instabile periodische Bahn liegt. In fast allen Phasen durchläuft der Kreisel Gebiete stabilen und instabilen Verhaltens. Diese Tatsache kann mit Fomenko-Graphen skizziert werden, siehe Abb. 7. Poincaré-Schnitte sind eine weitere Möglichkeit das Verhalten zu veranschaulichen, indem wie in einer Landkarte die stabilen Bereiche eingefärbt werden, an deren Grenzen die Separatrizen liegen, siehe Abb. 8.

Der Kowalewskaja-Kreisel ist von kreiseltheoretischer Bedeutung und kann zur Modellierung gewisser Phänomene in stationären Strömungen inkompressibler Flüssigkeiten verwendet werden; eine technische Anwendung gibt es jedoch Anfang des 21. Jahrhunderts noch nicht.[3]

Realisierungen

Bearbeiten

Das Trägheitsellipsoid des Kowalewskaja-Kreisels ist rotationssymmetrisch aber der Massenmittelpunkt befindet sich nicht auf der Symmetrieachse. Das ist nur bei inhomogener Massenverteilung möglich[4]. Auf Hermann Schwarz geht eine Konstruktion mit zwei Zylindern zurück[5] und P. H. Richter benutzte die Form mit einem Kreisring und drei Massenpunkten[6], die auch in den Animationen unten verwendet wurde. In den Abbildungen 1 und 2 bezeichnet s den Massenmittelpunkt und das Koordinatenkreuz den Stützpunkt.

Phänomenologie

Bearbeiten

Karussell-Bewegungen

Bearbeiten

Die Karussell-Bewegungen sind Staude-Drehungen des Kowalewskaja-Kreisels, bei denen der Kreisel um eine lotrechte Achse rotiert, die mit keiner der Hauptträgheitsachsen übereinstimmt. In Abb. 3 tut er das in stabiler Weise und in Abb. 4 ist der Schwerpunkt zu hoch und die Karussell-Bewegung wird instabil. In beiden Animationen wurde der Neigungswinkel zwischen der 1-Achse und der Lotlinie um 0,01 rad gegenüber dem anfänglichen Sollwert erhöht, was in Abb. 3 unbemerkt bleibt, in Abb. 4 jedoch zum periodischen Ausbrechen des Kreisels führt.

Karussell-Bewegungen entstehen unter den Bedingungen

n2 = q = p² + c n1 = k = 0.

siehe #Gleichgewichtslösungen.

Bifurkationsdiagramm des Kowalewskaja-Kreisels

Bearbeiten

 
Abb. 5: Bifurkationsdiagramm des Kowalewskaja-Kreisels (schematisch)

Einen Eindruck von der Komplexität des Verhaltens des Kowalewskaja-Kreisels vermittelt das Bifurkationsdiagramm in Abb. 5. Bifurkation (Verzweigung) ist das Phänomen, dass ein System an einem kritischen Punkt in verschiedene Zustände verzweigen kann. So kann ein lotrecht stehender Stab bei kleiner Störung nach links oder rechts umfallen sofern er irgendwie an eine Ebene gebunden ist. Die Trennlinie zwischen den beiden Pfaden ist eine Separatrix, von der auch beim unsymmetrischen Euler-Kreisel eine existiert. Die Bewegungen entlang der Separatrizen sind oft instabil.

Die kritischen Punkte beim Kowalewskaja-Kreisel stellen sich bei bestimmten Werten der Gesamtenergie h und des vertikalen Drehimpulses l ein. Das Bifurkationsdiagramm wird durch die #Gleichgewichtslösungen in fünf Flächen geteilt, wobei in der weißen Fläche O unphysikalische Parameterkombinationen auftreten, die dem Kreisel verboten sind. Die eingefärbten Bereiche sind dem Kreisel zugänglich und er zeigt in den roten, gelben, blauen und grauen Gebieten jeweils topologisch ähnliches Verhalten.

Der linke Rand a–b definiert den lotrecht hängenden Kreisel, dessen Bewegung immer stabil ist, und die Linie d–e den lotrecht aufrechten Kreisel, dessen Bewegung immer instabil ist. Charakteristisch für den Kowalewskaja-Kreisel sind die Karussell-Bewegungen, die auf der Linie u-v-w stattfinden und die auf der Linie u–v stabil und im Bereich v–w instabil sind[7], siehe Abb. 3 und 4. Entlang der zusätzlich eingezeichneten Linien verzweigen die periodischen Orbits, die der Kreisel beim minimalen oder maximalen Wert der Kowalewskaja-Konstante verfolgt, siehe Abb. 7.

An den im Bifurkationsdiagramm eingezeichneten Linien finden Bifurkationen statt:

In den Gebieten dazwischen sind die Energieflächen einheitlich geblättert.

Analytische Form des Bifurkationsdiagramms
Die Tabelle gibt die analytischen Darstellungen der Phasengrenzen in Abb. 5 für c = 1.
Bifurkation Darstellung Kowalewskaja-
Konstante
Definitionsbereich
a–b l = 2p, h = p2 - 1 k2 = ( p2 + 1 )2 p > 0
d–e l = -2p, h = p2 + 1 k2 = ( p2 - 1 )2 p > 0
u–v l = -p3 - p–1, h = 32 p2 + 12 p–2 k2 = 0 0 < p < 3–0,25
v–w l = -p3 - p–1, h = 32 p2 + 12 p–2 k2 = 0 3–0,25 < p < 1
B–C, B–J l2 = 1627 h3 k2 = 1 - h2/3 0 < h < √3
A–B, F–G l2 = 2h k2 = 0 h > 0
E-D, E–H l2 = h ± √(h2 - 2) k2 = l–4 h > √2

Familien gleicher Tori

Bearbeiten

Bei einer gegebenen Bewegung des Kowalewskaja-Kreisels sind die Gesamtenergie, der Drehimpuls in Lotrichtung und die Kowalewskaja-Konstante unveränderlich. Weil der Kowalewskaja-Kreisel ein integrables System darstellt, bewegt sich der Drehimpuls im körperfesten Hauptachsensystem auf einem ringförmigen Schlauch oder Torus und die Trajektorien füllen diesen dicht aus[8], siehe Abb. 6a. Der Torus kann sich dabei selbst durchdringen, die Außenseite ist heller als die Innenseite getönt. Bei gleicher Gesamtenergie und gleichem Drehimpuls in Lotrichtung verändern sich die Tori stetig mit der Kowalewskaja-Konstante und bilden eine Familie, was die Farbgebung andeuten soll. An einem Verzweigungspunkt, gehen die Tori in andere Familien über. Der Verzweigungspunkt ist eine Separatrix, die diesen angrenzenden Torifamilien gemeinsam ist, Abb. 6b und c. Im gezeigten Fall aus dem Gebiet B treffen in der Separatrix ein roter und zwei zueinander symmetrische grüne Tori zusammen. Im Zentrum der Separatrix liegt ein instabiler periodischer Orbit, den der Kreisel bei geringster Störung entlang der Separatrix verlässt, wobei er in einen der benachbarten Tori verzweigt. Nimmt die Kowalewskaja-Konstante zu oder ab, kann der Kreisel auf einen weiteren Verzweigungspunkt treffen, oder bei einem extremalen Wert der Konstanten in einem stabilen periodischen Orbit enden, wie dem grauen in Abb. 6d. Mit ansteigender Kowalewskaja-Konstante durchläuft der Kreisel in den Gebieten A bis J immer in gleicher Weise stabile Orbits, zwischen denen Separatrizen liegen. Das lässt sich mit Fomenko-Graphen versinnbildlichen.

Fomenko-Graphen

Bearbeiten

 
Abb. 7: Fomenko-Graphen

Die Abfolge der Bifurkationen bei zunehmender Kowalewskaja-Konstante lassen sich in Fomenko-Graphen darstellen, die in den Gebieten A bis J jeweils gleich geartet sind. In den freien Enden der Äste der Graphen finden stabile Bewegungen statt. Wo Äste verschiedener Farben aufeinander treffen findet eine Bifurkation entlang von Separatrizen statt. Die verschiedenen Farben gruppieren Bewegungstypen. Entlang der lila Äste beispielsweise finden karussellartige Bewegungen statt[9].

Poincaré-Schnitte

Bearbeiten

 
Abb. 8: Poincare-Schnitte entlang eines Pfades D-C-B-A im Bifurkationsdiagramm

Einen anschaulichen Eindruck vermitteln die Poincaré-Schnitte, siehe Abb. 8. Eine Trajektorie füllt im Phasenraum der Drehimpulse einen ringförmigen Schlauch, einen Torus, dicht aus. Die Tori werden durch die Kowalewskaja-Konstante k charakterisiert. Bei einem kritischen Wert von k entarten die Tori zu einer Separatrix, die einer instabilen periodischen Bahn entspricht.

Die Trajektorien des Kreisels werden mit einer Kugel verschnitten, die Oberfläche dieser Kugel gemäß der Familie der geschnittenen Tori eingefärbt und in die Ebene projiziert. Die Farben entsprechen den Familien der Tori und die Helligkeit wächst mit k. Die stabilen Orbits sind die hellen und dunklen Zentren der Flächen. Die Flächen werden durch Separatrizen getrennt, an denen deutliche Farbsprünge auftreten. Die Schnittpunkte von Separatrizen sind instabile periodische Orbits.

Zwischen D und C findet eine Bifurkation der Energiefläche statt und der gelbe Bereich verschwindet. Zwischen C und B verschwindet blau in einer Tangenten-Bifurkation. Zwischen B und A verschwindet grün in einer Pitchfork-Bifurkation[10].

Lösung der Bewegungsgleichungen

Bearbeiten

Euler-Poisson-Gleichungen

Bearbeiten

Aus Symmetriegründen kann die zu A gehörende erste Hauptträgheitsachse so gewählt werden, dass der Massenmittelpunkt des Kreisels auf ihr liegt. Damit spezialisieren sich die Euler-Poisson-Gleichungen beim Kowalewskaja-Kreisel zu[11]

 

Darin stehen p, q, r = ω1,2,3 für die Winkelgeschwindigkeiten, der Vektor   für die Richtung der Gewichtskraft des Kreisels, jeweils ausgedrückt im mitrotierenden Hauptachsensystem, auf das sich die Indizes 1,2,3 beziehen. Der Parameter

 

setzt sich aus der Masse m, der Schwerebeschleunigung g, dem Abstand s1 des Massenmittelpunkts vom Stützpunkt und dem Hauptträgheitsmoment A um die 1-Achse oder C um die 3-Achse zusammen. Die Konstante c quantifiziert die Wirkung der Gewichtskraft auf den Kreisel und mit c = 0 entsteht ein symmetrischer kräftefreier Euler-Kreisel.

Durch Skalierung der Zeit mit   und der Winkelgeschwindigkeiten mit   entstehen Bewegungsgleichungen mit c = 1, sodass sie keinen freien Parameter mehr besitzen. Mathematisch reicht es aus, nur diesen einen Fall zu untersuchen[12]. Kowalewskaja skalierte die Längen so, dass A = 2C = 2 wird[13].

Integrale der Bewegung

Bearbeiten
 
Abb. 9: Kurven mit k2 = 0,36 und n2 = 0 in Abhängigkeit von c n1 in der 1-2-Ebene

Wie bei jedem schweren Kreisel ist die Norm des Richtungsvektors der Gewichtskraft, der Drehimpuls in Lotrichtung und die Gesamtenergie E konstant:

 

Diese Konstanten werden in der Kreiseltheorie Integrale genannt, die ersten beiden auch Casimir-Invarianten. Die Gesamtenergie wird in der analytischen Mechanik auch als Hamilton-Funktion bezeichnet, was ihre Benennung mit h begründet. Beim Kowalewskaja-Kreisel gibt es noch ein viertes rationales Integral, die Kowalewskaja-Konstante

 

die also nie negativ ist und nur von den Winkelgeschwindigkeiten und Komponenten der Gewichtskraft in der 1-2-Ebene abhängt. Im pqr-Raum der Winkelgeschwindigkeiten um die Hauptachsen definiert diese Konstante momentan also eine zylindrische Fläche mit in r-Richtung konstantem Querschnitt. Wegen K(p,q) = K(-p,-q) ist die Fläche zudem symmetrisch, siehe Abb. 9. Bei negativem c n1 tauschen p und q die Rollen und es entstehen an einer Diagonalen gespiegelte Figuren. Kowalewskaja benutzte k² als Bezeichnung, es ist aber auch K gebräuchlich[14].

Gleichgewichtslösungen

Bearbeiten

Relative Gleichgewichte sind Fixpunkte der Kreiselgleichungen, die bei  , also bei Konstanz dieser Größen, auftreten[15]:

0 = q r
0 = -r p – c n3
0 = c n2
0 = r n2q n3
0 = p n3r n1
0 = q n1p n2

Somit folgt sofort n2 = 0 und bei q r = 0 sind zwei Fälle zu unterscheiden:

Fall 1, r=0
Aus r = 0 folgt weiter n3 = n2 = 0, also n1 = ±1, und daher schließlich auch q = 0. Dieser Fall entspricht dem lotrecht hängenden oder aufrechten Kreisel, was sich im Bifurkationsdiagramm in Abb. 5 auf den Linien a–b bzw. d–e wieder findet. Erstere Bewegung ist stabil, letztere instabil. Die Kowalewskaja-Konstante lautet hier k² = (p² ± c)², das Drehimpulsquadrat l² = 4p² und je nachdem l² < 4c oder l² > 4c unterscheidet sich der Typ der Instabilität. Die Bereiche werden in Abb. 5 durch den Punkt w getrennt, in dem eine Karussell-Bewegung mit k = 0 stattfindet.
Fall 2, q=0
Im zweiten Fall, q = 0, folgt r p = -c n3 und aus der vorletzten Euler-Poisson-Gleichung p² + c n1 = 0. In diesem Fall verschwindet daher die Kowalewskaja-Konstante k und es muss n1 < 0 sein, sodass der Kreisel ein aufrechter ist mit dem Massenmittelpunkt über dem Stützpunkt. Diese Bewegungen definieren die sogenannten Karussell-Bewegungen des Kreisels und finden sich im Bifurkationsdiagramm in Abb. 5 auf der Strecke u-v-w.

Lagrange- und Hamilton-Funktion des Kowalewskaja-Kreisels

Bearbeiten

Die theoretische Behandlung des Kowalewskaja-Kreisels geschieht mit der Hamilton-Funktion, die gleich der Gesamtenergie des Kreisels ist. Eine mögliche Formulierung benutzt die lokalen Drehimpulskomponenten L1,2,3 und den Richtungsvektor der Lotlinie[16]

 

mit

 

Bei Benutzung der Euler-Winkel ψ, ϑ und φ schreibt sich das

 

siehe Bezugssysteme und Euler-Winkel in der Kreiseltheorie. Darin sind

 

die konjugierten Impulse zu den Winkeln und

 

die Lagrange-Funktion des Kreisels.

Kowalewskajas Ausarbeitung

Bearbeiten

Sofia Kowalewskaja[1] gelang die Darstellung aller nötigen Variablen außer dem Präzessionswinkel ψ als hyperelliptische Funktionen der ersten Art, deren beide Argumente linear von der Zeit abhängen. Dazu nutzte sie Variablen s1,2, die zu einem gewissen Grad die Winkelgeschwindigkeiten p und q ersetzen[17]. Sie ging dabei wie folgt vor.

Die Winkelgeschwindigkeiten p und q sowie die Richtungskomponenten n1,2 werden mit der imaginären Einheit i jeweils zu einer komplexen Größe zusammengefasst:

 

Mit den #Euler-Poisson-Gleichungen können die Zeitableitungen

 

gebildet werden, aus denen n3 eliminiert werden kann:

 

Darin steht ln für den natürlichen Logarithmus. Das gleiche Verfahren führt mit den konjugiert komplexen Größen   zu   und die Summe liefert das reelle, nicht negative Integral

 

Mit den Variablen x1,2 konnte Kowalewskaja auch r, n3 sowie ξ1,2 ausdrücken.

Die Zeitableitungen stellte Kowalewskaja mit neuen Variablen s1,2 dar, die Nullstellen der Funktion 2. Grades

 

in s sind, worin die Koeffizienten

 

Polynome in x1,2 sind und l1 = h/3 ist. Für die Nullstellen s1,2 berechnen sich die Zeitableitungen

 

oder gleichbedeutend[18]

 

Darin ist

 

ein Polynom 5. Grades in s und

 

sind Konstanten der Bewegung. Der erste Faktor im Polynom R5 hat einen direkten Bezug zur Weierstraß’schen ℘-Funktion. Der Zusammenhang der Differentiale ds1,2 impliziert, dass sich s1,2 als Quotienten von Theta-Funktionen, deren Argumente lineare Funktionen der Zeit sind, ausdrücken lassen[19].

Damit ist das Problem im Prinzip gelöst. Kowalewskaja zeigte auf weiteren 50 Seiten, wie alle Variablen mit hyperelliptischen Funktionen dargestellt und mit ihnen beliebige Anfangsbedingungen berücksichtigt werden können. Die Materie ist derart kompliziert, dass es die Wissenschaft noch über ein Jahrhundert beschäftigte, diese Lösungen zu untersuchen und ihre Eigenschaften herauszuarbeiten[20].

Siehe auch

Bearbeiten

Die anderen beiden, bei beliebigen Anfangsbedingungen immer integrablen Kreiseltypen

oder den

Fußnoten

Bearbeiten
  1. a b Die Artikel
    Sophie Kowalevski: Sur le probleme de la rotation d’un corps solide autour d’un point fixe. In: Acta Mathematica. Volume 12, 1889, S. 177–232. (projecteuclid.org [abgerufen am 21. März 2018]).
    und
    Sophie Kowalevski: Sur une propriété du système d’équations différentielles qui définit la rotation d’un corps solide autour d’un point fixe. In: Acta Mathematica. Volume 14, 1890, S. 81–93 (projecteuclid.org [abgerufen am 21. März 2018]).
    enthalten die Ausarbeitungen, die am 25. Dezember 1888 mit dem Bordin-Preis ausgezeichnet wurden.
  2. siehe A. M. Permelov: Kovalevskaya Top: An Elementary Approach. In: Journal of Applied Mathematics and Physics. Vol. 131, Issue 2, 2002, S. 612–620, arxiv:math-ph/0111025v1 (Permelov vereinfacht die Bewegungsfunktionen soweit, dass sie nur wenige Seiten beanspruchen.).
    und
    Fawzy Mohamed Fahmy El-Sabaa, Alshimaa Abdelbasit Mohamed, Salma Khalel Zakria: The Qualitative and Quantitative Methods of Kovalevskys Case. In: Journal of Applied Mathematics and Physics. Vol. 5, Issue 9, September 2017, ISSN 2327-4379, doi:10.4236/jamp.2017.59155 (scirp.org [abgerufen am 23. März 2018] Studium des Phasen-Portraits, der singulären Punkte und des Bifurkationsdiagramms.).
  3. Tuschmann und Hawig (1993), S. 124.
  4. Magnus (1971), S. 126.
  5. Audin (2008), S. 99.
  6. Richter, Dullin und Wittek (1997), S. 36f.
  7. Richter, Dullin und Wittek (1997), S. 11ff.
  8. Richter, Dullin und Wittek (1997), S. 56f.
  9. Richter, Dullin und Wittek (1997), S. 25f.
  10. Richter, Dullin und Wittek (1997), S. 22f.
  11. Magnus (1971), S. 126, Richter, Dullin und Wittek (1997), S. 12.
  12. Richter, Dullin und Wittek (1997), S. 36.
  13. Kowalewskaja (1888), S. 184.
  14. Richter, Dullin und Wittek (1997), S. 10.
  15. Richter, Dullin und Wittek (1997), S. 11f.
  16. Richter, Dullin und Wittek (1997), S. 36. Dort sind die Bezeichnungen der Winkel ψ und φ jedoch gegenüber der hiesigen vertauscht.
  17. Leimanis (1965), S. 39.
  18. Kowalewskaja benutzte die Bezeichnung R1(s) statt R5(s). Für die Wurzeln ist nicht immer das positive Vorzeichen richtig. Wird das Vorzeichen entsprechend angepasst, stimmen die Gleichungen.
  19. Tuschmann und Hawig (1993), S. 129.
  20. siehe beispielsweise Permelov (2002) und Richter, Dullin und Wittek (1997).

Literatur

Bearbeiten