Keplersche Gesetze zur Bestimmung des Kreisbahnradius

„Im nichtrelativistischen Grenzfall der Allgemeinen Relativitätstheorie gibt es eine Korrektur der Newtonschen Gravitationspotentialenergie H / R 3 mit H = a L 2 / ( M C ) 2 , Wo C ist die Lichtgeschwindigkeit, a = G M M Und L ist der Drehimpuls"

Mit diesem Wissen soll ich den Radius von Kreisbahnen für eine gegebene Zeit finden M Und L und entscheiden, welche davon stabil ist.

Meine Frage hat damit zu tun, wie ich die Radien tatsächlich bestimmen kann, aber ertragen Sie mich, wie ich meinen bisherigen Prozess zeige:

Ich schickte meinem Professor eine E-Mail und ihm wurde gesagt, dass ich diesen Korrekturfaktor von der potenziellen Gravitationsenergie SUBTRAHIEREN muss, was mir Folgendes gibt:

v ( R ) = G M M R G M M L 2 R 3 ( M C ) 2

Ich kann das effektive Potenzial finden, um zu sein:

v eff ( R ) = L 2 2 M R 2 G M M R G M M L 2 R 3 ( M C ) 2

Basierend auf den Informationen in meinem Lehrbuch würde ich mir vorstellen, dass ich das effektive Potenzial und die gerade Linie meiner konstanten Energie grafisch darstellen muss E , und die beiden Radien sind die Punkte, an denen die Energielinie die schneidet v eff ( R ) Kurve. Mein Problem tritt auf, wenn ich versuche, es grafisch darzustellen.

Wenn ich einen qualitativen Graphen zeichne, einfach indem ich benutze 1 R 2 1 R 1 R 3 Ich erhalte eine Kurve ohne offensichtliche Extrema, die sich 0 nähert j Achse als R geht ins Unendliche. Mit anderen Grafiken, die ich gezeichnet habe v eff , ich habe Kurven bekommen, die Sinn machen - ich sehe ein lokales Maximum oder Minimum, und ich gehe davon aus, dass der Planet zwischen den beiden "Wänden" des Minimums "gefangen" bleiben könnte. In diesem Fall sehe ich davon natürlich nichts.

Meine Frage ist (hoffentlich) viel allgemeiner als nur dieses Beispiel: "Wie finde ich die potenziellen Radien des umlaufenden Planeten mit Keplers Gesetzen?" Wenn ich in meinem Prozess keinen Fehler gemacht habe, glaube ich nicht, dass ich sie mit dieser Methode finden kann. Ich könnte mir vorstellen, dass ich sie mit viel Kalkül und Neuordnung finden könnte, aber ich bin sicher, dass es einen einfacheren Weg geben muss.

Was ist das 1 / R 2 Begriff für?
v e F F ( R ) = L 2 2 M R 2 + v ( R ) , also die 1 / R 2 Begriff kommt von L 2 2 M R 2 !
Wenn wir von der "nicht relativistischen Grenze" der allgemeinen Relativitätstheorie sprechen, impliziert dies, dass die 1 / R 3 viel kleiner ist als die beiden anderen Terme, damit diese Annäherung gültig ist
Versuchen Sie es mit Plotten R 2 R 1 0,01 R 3 stattdessen. Die Korrektur ist gering. C 2 ist klein.
1 / R 2 : Aha. Ihr OP hatte einen Tippfehler ...
Ah ich sehe! Ich habe einige dieser Werte falsch berechnet und etwas um E92 herum bekommen. Das führt mich zum zweiten Teil meiner Frage – ich habe versucht, es anders anzugehen und festgestellt, dass seit einer kreisförmigen Umlaufbahn dies der Fall ist R M ich N = R M A X , wären die beiden möglichen Radien einfach an den Extrema, richtig?
Exakt; Suchen Sie nach den Extrema des effektiven Potentials.

Antworten (2)

Deine Gleichung hat die Form

v eff ( R ) = a R 2 β R γ R 3
Wenn Sie einstellen a = β = γ = 1 , dann überschätzen Sie die R 3 Begriff, der eine kleine Korrektur sein soll. Sie werden nur dann zwei Extrema finden, wenn die Ableitung zwei Nullstellen hat:
v eff ' ( R ) = 2 a R 3 + β R 2 + 3 γ R 4 = 0
was impliziert
β R 2 2 a R + 3 γ = 0
Diese Gleichung hat die Diskriminante
Δ 2 = 4 ( a 2 3 β γ )
So v ( R ) hat zwei Extrema if
a 2 3 β γ
was wahr ist, wenn γ ausreichend klein ist.

Vielen Dank, dass du mir hier geholfen hast! Das bringt mich zu meiner zweiten Frage, wie man die Radien findet. Da die Umlaufbahnen kreisförmig sind, R M ich N = R M A X Sie wären also nur am Extrem, oder? Gibt es einen einfacheren Weg, dies zu lösen? Ich hatte grafisch kein Glück, da die Werte so groß sind und die Formel zu kompliziert ist, um eine schöne Ableitung gleich Null zu ergeben.
@Alex Ja, die Radien sind an den Extrema. Einfach lösen v eff ' ( R ) = 0 .

Das Finden der möglichen Radien ist eigentlich ganz einfach. Ich habe schon:

v e F F ( R ) = L 2 2 M R 2 G M M R G M M L 2 R 3 ( M C ) 2

Ich habe irrtümlich grafisch dargestellt 1 / R 2 1 / R 1 / R 3 , obwohl es in Wirklichkeit sinnvoller ist, zu nehmen 1 / R 2 1 / R 0,1 / R 3 , seit 1 / R 3 ausreichend klein ist.

Da die Umlaufbahnen kreisförmig sind, liegen die potentiellen Radien an den Extrema der Funktion v e F F ( R ) , also muss ich einfach nehmen:

v e F F ' ( R ) = 0
L 2 M R 3 + G M M R 2 + 3 G M M L 2 R 4 ( M C ) 2 = 0
Multiplizieren Sie jede Seite mit R 4 , kann ich sehen, dass ich eine ähnliche Gleichung habe wie:
A R 2 + B R + C = 0
Also kann ich die quadratische Formel verwenden, um zu erhalten:
R = | L 2 / M + L 4 C 2 12 G 2 M 2 M 2 L 2 / M C 2 G M M |
In Bezug auf mein Diagramm kann ich sehen, dass der erste dieser Radien instabil sein wird, weil die Energie größer ist als er, also ins Unendliche geht, während für den zweiten (größeren) Radius die Gesamtenergie größer als das lokale Minimum ist also Der Planet / das Teilchen wird in der Potentialmulde gefangen sein, also ist es stabil.

Wenn ich akzeptierte Werte für G, M, m, c und L in die Gleichung einsetze, die ein größeres r ergibt (wobei ich am +-Zeichen subtrahiere), erhalte ich den tatsächlichen Radius von der Erde zur Sonne, ~ 1.5 10 8 k M !