Delta-v-Budgets und modifizierte Schwerkraft

Ich baue eine Soft-Sci-Fi-Welt, in der die Schwerkraft etwas anders funktioniert: Die Beschleunigung aufgrund der Schwerkraft fällt mit einer Rate ab, die viel größer ist als unser umgekehrtes Quadrat; und die Beschleunigung aufgrund der Schwerkraft "kippt" und oszilliert in regelmäßigen Abständen und wird eher abstoßend als anziehend und umgekehrt (ähnlich wie MOND , aber in kleineren planetarischen Maßstäben).

Die Beschleunigung G Die Erfahrung eines Objekts mit vernachlässigbarer Masse kann modelliert werden mit:

G S e λ D u cos ( D u + ϕ )

Wo

  • S - Gravitationsbeschleunigung auf einer Planetenoberfläche 1 (in M / S 2 ),
  • λ - dimensionslose Zerfallskonstante (beschreibt, wie allmählich der Kraftabfall ist),
  • D - Abstand von der Oberfläche in Metern,
  • u - Skalarkonstante 2 (wandelt Meter in Bogenmaß um),
  • ϕ - Phasenwinkel, fast immer Null und kann ignoriert werden.

1 Wenn Sie der Mathematik folgen, oszillieren die Gravitationskräfte unter der Oberfläche mit exponentiell größeren Magnituden, was einen Planeten offensichtlich in Stücke reißen würde. Die Formel ist also nur gut für die Schwerkraft über der Oberfläche. Nachfolgend gelten andere (und unwichtige) Regeln.

2 u = 0,5 π / D 0 , Wo D 0 ist der Abstand zwischen zwei Null-g-Punkten oder 1 / 2 die Kosinuswellenlänge.

e λ D u Und cos ( D u ) sind dimensionslos, sodass das Ergebnis in Einheiten von bleibt M / S / S . Das Multiplizieren des Ergebnisses mit der Masse des Objekts ergibt eine auf das Objekt ausgeübte Kraft, die positiv oder negativ sein kann.

Da es sich hier um Worldbuilding handelt, sind viele dieser Parameter umschaltbar und entstehen (scheinbar) nicht aus grundlegenden Interaktionen. Das ist okay, denke ich. Ich muss nicht vollständig ausarbeiten und begründen, warum diese Konstanten so sind, wie ich sie gewählt habe (oder warum die Schwerkraft anders ist), ich kann sie einfach in eine schwarze Ecke stecken und mit ihrer Ausgabe arbeiten.


Es gibt einige seltsame Eigenschaften bei dieser Form der Gravitation. Hier ist eine Grafik:

Gravitationsdiagramm2

λ = 0,5 , S = 1

Wenn sich ein Objekt über die Oberfläche erhebt, sinkt die Gravitationsbeschleunigung bis in einiger Entfernung π / u , die Beschleunigung wird Null und kehrt sich um und wird abstoßend. Die Kraftvektoren zeigen von der Planetenoberfläche weg und das Objekt wird wegbeschleunigt, bis es den wiederkehrenden Null-g-Punkt kreuzt und mit der Verzögerung beginnt.
Da jedoch die Entfernung, über die Beschleunigung/Verzögerung auftritt, konstant ist und die Größe der Beschleunigung/Verzögerung immer abnimmt, wird ein nicht antreibendes Objekt immer vom Planeten wegbeschleunigen – für immer – sobald es den anfänglichen Null-G-Punkt überschritten hat über der Planetenoberfläche. Um sich dem Planeten überhaupt zu nähern, muss ein Objekt eine gewisse minimale „Eintrittsgeschwindigkeit“ haben.


Wie sehen Umlaufbahnen aus?

Schwerkraftdiagramm 1

Die Schwerkraft oszilliert zwischen positiv und negativ, anziehend und abstoßend. (In der Grafik sind attraktive Regionen solide.) Stabile Umlaufbahnen scheinen in zwei Varianten vorzukommen: Umlaufbahnen, die ausschließlich innerhalb attraktiver Regionen verbleiben (weitgehend kreisförmig) und Umlaufbahnen, die zwischen anziehenden und abstoßenden Regionen oszillieren (weitgehend sinusförmig). Hochelliptische Bahnen sind unmöglich; Wenn ein Objekt genug Orbitalenergie hat, um mehr als einen Null-G-Punkt zu überqueren, wird es aus dem System ausgeworfen. Flugbahnen, die fast einen zweiten Null-G-Punkt überqueren, sind „kritisch“.

kritischer Weg 1 kritischer Weg 2

(Attraktive Regionen in Schwarz, abstoßende in Dunkelgrau. Pfade, die mit einem nicht rotierenden Körper in der Mitte verfolgt werden.)
Hier sehen wir ein Objekt (zweites Bild) auf einer kritisch oszillierenden Flugbahn, das auf den Planeten fällt, obwohl es nicht genug kinetische Energie hat, um es zu überwinden die abstoßende Region und wird wieder heraus beschleunigt – obwohl nicht stark genug, um in die zweite abstoßende Region einzudringen – ad infinitum. Die Tangentialgeschwindigkeit scheint fast unabhängig von der Radialgeschwindigkeit zu sein, was bedeutet, dass zwei Objekte eine Orbitalregion teilen und dennoch sehr unterschiedliche Geschwindigkeiten im Vergleich zueinander haben können. Daher das dickere Band nachgezeichneter Umlaufbahnen (die eine schnellere Bewegung darstellen) innerhalb der größeren, langsameren, spärlicheren Flugbahn.

Stabile Umlaufbahnen gibt es nur in attraktiven Regionen, in denen die Kraftvektoren auf den Planeten zeigen. Kreisbahnen in diesen Regionen können aufgrund des großen Gravitationsgradienten sehr unterschiedliche Umlaufgeschwindigkeiten haben. An der Spitze des Bereichs, wo der Kraftvektor am stärksten ist, ist die Orbitalgeschwindigkeit am höchsten, wohingegen zu beiden Seiten der Spitze die Geschwindigkeiten abnehmen und sich Null nähern. Geosynchrone Trajektorien sind nicht auf den Raum der Null-g-Regionen beschränkt, sondern können auch zwischen anziehenden und abstoßenden Regionen oszillieren. Auf der Planetenoberfläche würde dies wie ein Satellit aussehen, der im Laufe des Tages näher und weiter wächst (in derselben Position oder entlang einer Achterbahn).

wenige geostationäre 1 wenige geosynchrone 2

In der Nähe der höheren Umlaufbahnen befinden sich einige synchrone Flugbahnen, die fast "auf der Stelle springen".


Womit ich Probleme habe

Ich denke, mein Problem ist hauptsächlich ein mathematisches (ich weiß über Mathematik so viel wie ein fortgeschrittener Programmierer sollte, aber nicht genug, um die Gleichungen zu analysieren und abzuleiten, die ich brauche, um meine eigene Frage zu beantworten, denke ich).
Ich möchte, dass meine Reisenden solche Objekte erkunden, und obwohl die Eigenschaften dieses Gravitationssystems magisch sind, sind die Reisenden selbst durch die Tsiolkovsky-Raketengleichung eingeschränkt .
Was ich brauche, ist eine Methode zur Berechnung von Delta-V-Budgets unter Berücksichtigung der modifizierten Schwerkraft.

Angenommen, meine Reisenden wollen vom Planeten in eine Höhe über der Oberfläche abheben (angenommen, es gibt keine Atmosphäre). Sie bräuchten eine Methode, um herauszufinden, wie viel Delta-v erforderlich ist, um dorthin zu gelangen (was wahrscheinlich eine Verlangsamung an "dort" beinhaltet, da die abstoßenden Regionen daran arbeiten, sie aus dem System zu beschleunigen) sowie wie viel Delta -v müssten sie zurückkommen . Denken Sie daran, dass die abstoßenden Regionen ihnen dabei helfen werden, den Planeten zu verlassen (und vielleicht sogar ein wenig bei der Zirkularisierung ihrer Umlaufbahn), aber sie werden auf dem Weg zurück nach unten gegen dieselben Regionen kämpfen.
Vielleicht wäre das einfachste Problem, herauszufinden, wie viel Delta-v erforderlich ist, um sich radial von einer Region – abstoßend oder anziehend – zu einer anderen zu bewegen.
Der Einfachheit halber machen wir uns nur Sorgen um koplanare Änderungen der Flugbahn.

Ich denke, Sie müssen die Vis-Viva-Gleichungen (das ist der Kompromiss zwischen potenzieller und kinetischer Energie) neu formulieren. Ist Ke in deiner Welt immer noch halb m V zum Quadrat?
@JamesMcLellan Ja, ich denke schon. Meine Absicht war nur, das Verhalten der Schwerkraft zu ändern. (Eigentlich ist das falsch, aber ich hielt die Unterscheidung nicht für zu wichtig, um sie in die Frage aufzunehmen. In meiner fiktiven Welt gibt es immer noch reguläre Schwerkraft. Diese modifizierte Schwerkraft ist entweder die wahre Einsteinsche Schwerkraft, die sich (irgendwie) auf diese andere Weise verhält . , oder ein zusätzliches, überlagerndes Feld, das ähnlich wie die Schwerkraft wirkt, nur eben auf diese modifizierte Weise.)
Schlagen Sie das vor, weil die Erhaltung der mechanischen Energie verletzt wird?
Wollte nur bestätigen
sieht aus wie Elektronenwellen ;D
Mit welcher App zeichnest du?
Nur um sicherzugehen, wir haben es mit einem reinen 1-Körper-Problem zu tun, richtig? Kein Zentralstern oder andere Planeten, um die man sich kümmern muss?
@Alice stimmt. Keine anderen störenden Systeme.

Antworten (2)

Ich denke, Sie müssen die vis-viva-Gleichungen neu formulieren. Probieren wir das aus und sehen, was passiert:

Kreisbahn :

  • G = G M R 2
  • v 2 = G M R v 2 = G R
  • G S e λ D u cos ( D u + ϕ )

Seit D Ist R R , u ist ein Frequenzterm, S Ist G M R 2 Und ϕ kann ignoriert werden, ich werde einfach ersetzen:

  • G G M R 2 e λ ( R R ) u cos ( ( R R ) u + 0 )
  • v 2 = [ G M R 2 e λ ( R R ) u cos ( ( R R ) u + 0 ) ] R

Probieren wir es aus! Verwenden Sie die Erde als Referenz (M = 5.9 × 10 24 , R=6378km), versuche es mal mit den Sinuseffekten weit weg u = π R , und r an der Oberfläche (6378 km ~ 6.378.000 m); Ich bekomme 9,67 m/s/s für g, das ist derselbe Wert, den ich beim Rechnen erhalte G M R 2 . Ich erhalte eine Umlaufgeschwindigkeit von 7.855 m/s; was passt G M R

Wenn ich mich bei r = 1,5 R, g ~ 0 und v ~ 0 in den sinusförmigen Effekt bewege, ist das genau das, was ich erwarten würde.

Versuchen wir also mit etwas mehr Selbstvertrauen, dies anzuwenden auf ...

Elliptische Umlaufbahn:

  • v 2 = G M ( 2 R 1 A ) v 2 = ( 2 G M R G M A ) v 2 = ( 2 G R G R 2 A )
  • G G M R 2 e λ ( R R ) u cos ( ( R R ) u )
  • v 2 = ( 2 G M R 2 e λ ( R R ) u cos ( ( R R ) u ) R G M R 2 e λ ( R R ) u cos ( ( R R ) u ) R 2 A )

Nochmals testen, immer noch die Erde verwenden und eine elliptische Umlaufbahn versuchen, bei der die Hauptachse gleich R ist (eine kreisförmige Umlaufbahn): 7.855 m / s ... genau das, was ich von einer kreisförmigen Umlaufbahn erwarten würde.

Wie wäre es mit dem Sinuseffekt bei r = 1,5 R? v ~ 0 m/s; was ich erwarten würde.

Überweisungen:

Boosten aus einer niedrigeren Umlaufbahn R 1 zu einem Höheren R 2

  • Δ v 1 = μ R 1 ( 2 R 2 R 1 + R 2 1 ) Δ v 1 = G R 1 ( 2 R 2 R 1 + R 2 1 )
  • Δ v 1 = ( G M R 2 e λ ( R R ) u cos ( ( R R ) u ) ) R 1 ( 2 R 2 R 1 + R 2 1 )

Und die neue Umlaufbahn kreisförmig machen

  • Δ v 2 = μ R 2 ( 1 2 R 1 R 1 + R 2 ) Δ v 2 = G R 2 ( 1 2 R 1 R 1 + R 2 )
    Δ v 2 = ( G M R 2 e λ ( R R ) u cos ( ( R R ) u ) ) R 2 ( 1 2 R 1 R 1 + R 2 )

Fluchtgeschwindigkeit:

Kann hier nicht einfach das neue g ersetzen.

  • 1 2 M v e 2 = G δ R

Kreditaufnahme BMFs ausgezeichnete Antwort. G = S e λ D u ( Sünde ( D u + ϕ ) λ cos ( D u + ϕ ) ) ( λ 2 + 1 ) u

Ersetzen von d und S und Weglassen von Phi

  • G δ R = G M R 2 e λ ( R R ) u ( Sünde ( ( R R ) u ) λ cos ( ( R R ) u ) ) ( λ 2 + 1 ) u

Dann wieder in die Fluchtgeschwindigkeitsgleichung einsetzen

  • 1 2 M v e 2 = G M R 2 e λ ( R R ) u ( Sünde ( ( R R ) u ) λ cos ( ( R R ) u ) ) ( λ 2 + 1 ) u

Fallenlassen von m und Auflösen nach v e

  • v e 2 = 2 [ G M R 2 e λ ( R R ) u ( Sünde ( ( R R ) u ) λ cos ( ( R R ) u ) ) ( λ 2 + 1 ) u ]
Wow, das sieht vielversprechend aus. Ich werde später Zeit haben, es weiter zu analysieren, aber die ersten beiden Grafiken von G Und v 2 schau mal nach rechts
Wie sind Sie in den Hohmann-Transfergleichungen zum Substituieren gekommen? μ R N mit G R N ? (Auch mit welchem ​​Radius 1 oder 2 ist G berechnet?)
Ich habe die Erweiterung von "g"= ... nicht gemacht. Ich dachte, es wäre in den kreisförmigen und elliptischen Beispielen zweimal gemacht worden. Ich hatte das Gefühl, dass es chaotisch war. Ich werde es hinzufügen.
Ich denke, was mich verwirrt, ist, dass g einen Radiusparameter mit kleinem r annimmt, und mit den Übertragungsgleichungen haben wir 2 kleine r
@BMF darf ich mir dein Integral ausleihen?
Ja, natürlich habe ich es geschrieben, falls es bei einer Antwort helfen könnte :)
Innerhalb der Hohmann-Transfergleichungen G könnte man sagen, eine Funktion von zu sein R , Und R 1 Und R 2 und die unteren und höheren Umlaufbahnen, machen Δ v das momentane delta-v, das benötigt wird, um es zu erreichen R bei einer Flugbahnänderung aus R 1 Zu R 2 ?

Ich glaube, ich könnte einer Antwort nahe sein, aber ich bin mir nicht ganz sicher. Zumindest kann dies jemandem zu einer vollständigeren Antwort verhelfen.
Ursprünglich dachte ich daran, meine modifizierte Erdbeschleunigungsformel zu integrieren

G S e λ D u cos ( D u + ϕ )

über die Dimension D , das ist die Entfernung in Metern. Ich konnte die Rückkehr davon nicht verstehen (ich dachte, es könnte nützlich sein, um die Geschwindigkeit zu finden), also ging ich mit dieser Frage weiter, aber jetzt glaube ich, dass ich ein bisschen besser verstehe, was dieses Ergebnis haben könnte gewesen.

G S e λ D u cos ( D u + ϕ ) D D = S e λ D u ( Sünde ( D u + ϕ ) λ cos ( D u + ϕ ) ) ( λ 2 + 1 ) u + v ich

v ich ist die "Anfangsgeschwindigkeit". Wenn es auf das Minimum des Diagramms entlang der positiven x-Achse eingestellt ist, liegt das gesamte Geschwindigkeitsdiagramm über der x-Achse, was meiner Meinung nach bedeutet, dass der Mindestwert die "Eintrittsgeschwindigkeit" ist, die ein Objekt ab unendlich erreichen muss D = 0 .

Beschleunigung über Geschwindigkeit

Die Gleichung für die Erdbeschleunigung in Rot (für die Sichtbarkeit um den Faktor 20.000 aufgeblähte y-Werte), die Integration dieser Gleichung in Blau.

Der Graph der integrierten Formel hat einen Minimalwert entlang der positiven x-Achse, den ich als die minimale Geschwindigkeit interpretiert habe, die erforderlich ist, um sich dem Planeten zu nähern (vorausgesetzt, die integrierte Formel stellt die Momentangeschwindigkeit dar). Dieser Mindestwert existiert in der Grafik bei D = D 0 , Wo D 0 ist der Abstand in Metern zwischen zwei Zero-G-Punkten (ich glaube, sie werden Wendepunkte genannt).
Bewerten Sie also die integrierte Formel für D = D 0 und wenn wir den Graphen um diesen Betrag auf der y-Achse nach oben verschieben, sollten wir ein Diagramm der Geschwindigkeit eines Objekts finden, das aus dem Unendlichen auf die Oberfläche des Planeten fällt, die oben erwähnte "Eintrittsgeschwindigkeit".

Wenn das Objekt weniger als die Eintrittsgeschwindigkeit hat, dann ist die Geschwindigkeit des Objekts am letzten Null-g-Punkt negativ, was ich so interpretiert habe, dass das Objekt nicht genug Energie hat, um es zu übertreffen, und zurückgestoßen wird, woher es kam. Wenn wir der Handlung rückwärts folgen, sollten wir die Mindestgeschwindigkeit für ein Objekt bei jedem Schritt finden, die erforderlich ist, um den Planeten vollständig zu verlassen. Diese Energie scheint dem Objekt unmittelbar nach dem Überschreiten des ersten Null-G-Punktes zugeführt zu werden.

Welches Programm verwendest du da?
@JamesMcLellan Desmos . Insbesondere das Grafiktool.
G kann aus dem Gravitationspotential ermittelt werden v von G = v , was Sie also gefunden haben, ist (mit einem Minuszeichen) das Gravitationspotential.
@ HDE226868 Ja, ich denke, du hast Recht. P = G D , was sehr ähnlich aussieht P E G = M G H Sie müssten irgendwie eine Gleichung für die Beschleunigung über die Zeit extrapolieren (anstelle über die Höhe über dem Boden) und diese integrieren. Scheint kein einfacher Weg zu sein, das zu tun.