Budgetierung von Delta-V für atmosphärischen Widerstand, wenn die Umlaufbahn elliptisch ist?

Ich versuche, das Delta-V-Budget zu berechnen, um eine generische elliptische Umlaufbahn aufrechtzuerhalten. Wenn ich mir die atmosphärischen Luftwiderstandsverluste ansehe, kann ich nur Gleichungen finden, die den Delta-V-Verlust für kreisförmige Umlaufbahnen berechnen. SMAD beantwortet diese Frage für Kreisbahnen, scheint sie aber nicht für Ellipsen zu öffnen.

Wertz hat in seinem Artikel 'moderately low elliptical very low orbits' eine Tabelle mit dem Delta-V für einige elliptische Bahnen und sagt: 'Siehe die Referenz für die Gleichungen für eine elliptische Bahn' (Seite 6 gleich nach Gl. 1) und dann macht nicht klar, wo die Gleichungen zu finden sind. Dies impliziert, dass es allgemeine Gleichungen gibt und diese nicht numerisch gelöst wurden.

Referenz https://digitalcommons.usu.edu/cgi/viewcontent.cgi?article=1040&context=smallsat

Antworten (1)

Es ist übrigens ein interessantes Papier.

Ich denke, "die Referenz" bezieht sich auf die Referenz 25, die unmittelbar davor hochgestellt erscheint! Wertz, James R., 2011a. „Atmosphärischer Luftwiderstand und Satellitenzerfall“, Sec. 9.4.4 in Space Mission Engineering: The New SMAD, hrsg. von J. Wertz, D. Everett und J. Puschell, Hawthorne, CA: Microcosm Press, 2011.


In der ursprünglichen SMAD sieht es jedoch so aus, als ob Gl. 6.21 bis 6.23 geben bereits das, wonach Sie fragen:

a D = 1 2 ρ ( C D EIN / m ) v 2           (6-21)

Wir können die Änderungen der großen Halbachse und der Exzentrizität pro Umdrehung sowie die Lebensdauer eines Satelliten mit den folgenden Gleichungen annähern:

Δ a r e v = 2 π ( C D EIN / m ) a 2 ρ p exp ( c ) ( ich 0 + 2 e ich 1 )           (6-22)

Δ e r e v = 2 π ( C D EIN / m ) a ρ p exp ( c ) ( ich 1 + e ( ich 0 + ich 2 ) / 2 )           (6-23)

wo ρ p ist die atmosphärische Dichte am Perigäum, c = a e / H , H die Höhe der Dichteskala ist (siehe Spalte 25, Innenseite der hinteren Umschlagseite) und ich ich sind modifizierte Bessel-Funktionen der Ordnung ich und Argument c . Wir modellieren den Begriff m / ( C D EIN ) , oder ballistischer Koeffizient als Konstante für die meisten Satelliten, obwohl er je nach Ausrichtung des Satelliten um den Faktor 10 variieren kann (siehe Tabelle 8-3).


Das gibt dir Δ a r e v und Δ e r e v , die Änderung in a und e pro Umdrehung, aber es wäre immer noch an Ihnen, einen Wert für zu wählen m / ( C D EIN ) oder ballistischen Koeffizienten und zu versuchen, das variable Dichteprofil der Atmosphäre mit einer lokalen Skalenhöhen -Näherung darzustellen, die für das Funktionieren dieser analytischen Näherung erforderlich ist, und nicht mit einem geeigneten numerischen Modell, das auf einem realistischen Dichteprofil basiert).

unten: Aus Ihrem verlinkten Artikel https://en.wikipedia.org/wiki/Scale_height Wertz et al. SSC12-IV-6, 26. jährliche AIAA/USU-Konferenz über kleine Satelliten. Es gibt einige gerade Abschnitte auf diesen Log-Lin-Plots, wo eine Annäherung an die Skalenhöhe nicht allzu unvernünftig wäre.

atmosphärische Dichte vs. Höhe

Danke, uhoh, ich glaube, ich habe es in meinem Kopf zu kompliziert gemacht, als ich nur die gegebenen Gleichungen und das Vis-Viva verwenden musste. Außerdem habe ich jetzt das neue SMAD gefunden, dort gibt es ein paar zusätzliche Informationen zum Thema.
@Boomtown das ist großartig! Die Magie von SE. Ich habe hier sicherlich viel gelernt; SMAD habe ich noch nie gesehen. Ich bin daran interessiert zu sehen, dass Bessel-Funktionen jetzt in einem Störungsproblem gelandet sind. Ich werde versuchen, das aufzuspüren oder zu reproduzieren, oder noch besser, eine neue SE-Frage zu stellen. Die Möglichkeiten sind endlos!