Analytische Annäherungen für die Formen dieser Diagramme von ISP vs. Düsendurchmesser? (Raptor-Motor)

Ein Kommentar zur derzeit unbeantworteten Frage Gibt es einen verbindlichen Katalog für Raketenstartkennungen? führte mich zu dieser Github -Seite, die auf r/spacexlounge verweist , auf der derzeit ein Beitrag mit einem Diagramm mit dem Titel Raptor-Effizienz für verschiedene Düsengrößen erscheint , der ebenfalls unten gezeigt wird.

  1. Gibt es für die beiden Grenzfälle der wohldefinierten Bedingungen des Umgebungsdrucks: Meereshöhe und Vakuum analytische Annäherungen an diese beiden Kurven? Die Kurven erinnern an Parabeln, wobei die Ordinate etwas ähnliches ist 1 ( 1 / D 2 A ) 2 Wo D könnte Durchmesser und sein A könnte mit dem Umgebungsdruck zusammenhängen, aber das ist nur das Putten mit geschlossenen Augen.

  2. Wie wurde (oder könnte zumindest) die Zwischenkurve für einen effektiven ISP definiert?

Hinweis: Das Erlernen der Originalquelle für die Handlung wird interessant sein, aber meine Frage bezieht sich hauptsächlich auf Raketendüsenmathematik und Physik im Allgemeinen.

Raptor-Effizienz für verschiedene Düsengrößen

HINWEIS: Wie in den Kommentaren unten angegeben, scheint das Diagramm unten eher ein individueller Beitrag als offiziell veröffentlichte Daten zu sein. Zitieren Sie nicht daraus oder verwenden Sie es als Prämisse für Folgefragen, ohne dies zu berücksichtigen.

Ich gehe davon aus, dass "effektives ISP" den zeitlichen Durchschnitt des ISP im Verlauf des Aufstiegs vom Meeresspiegel zum Vakuum darstellen soll; es könnte entweder eine Ad-hoc-Gewichtung der beiden ISP-Zahlen sein oder durch eine diskrete Stufensimulation des Aufstiegs bestimmt werden.
In den Kommentaren zum reddit-Beitrag sieht es so aus, als hätte OP einige erhebliche naive Fehler in ihrer (nicht verwandten) Turbopumpen-Effizienztabelle gemacht, also möchten Sie vielleicht, dass ein Körnchen Salz zu Ihrer Handlung passt.
@RussellBorogove danke dafür. Ich interessiere mich hauptsächlich für einfachere analytische Gleichungen, die diese Art von Verhalten zeigen, anstatt mich mit der Richtigkeit dieser bestimmten Kurven zu befassen, aber es ist wichtig zu wissen, dass diese individuell generiert und nicht offiziell sind, daher werde ich dies in der Frage vermerken.
Ich denke, dies zeigt nur die Auswirkung der Änderung des Expansionsverhältnisses, aber ohne etwas über den Halsbereich, die Durchflussraten usw. zu wissen, muss ich nachsehen, ob alles ausfällt oder zurückgesetzt werden kann. Eine Aufgabe für morgen.
@OrganicMarble Ich denke, wenn es Gleichungen gibt, die das gleiche Verhalten haben (Atm-Upside-Down-Kurvenspitzen bei endlichem Durchmesser, Vac-Kurvenspitzen bei unendlich), wäre das in Ordnung, und die Vac-Kurve wäre immer größer als die Atm-Kurve noch besser. Sie müssen überhaupt nicht exakt übereinstimmen, dafür gibt es möglicherweise zu viele motorspezifische Details.

Antworten (1)

Die einzige Gleichung, die benötigt wird, um diese Kurven zu erzeugen, ist die Schubgleichung in einigen verschiedenen Formen.

  1. F = M ˙ v e + ( P e P A T M ) A e
  2. ICH S P = F / ( G 0 M ˙ )

Sie benötigen auch isentrope Flussdiagramme oder Tabellen, um die Düsenparameter zu berechnen, die sich mit einem sich ändernden Flächenverhältnis ändern. Ein Beispiel für ein isentropisches Flussdiagramm ist hier aus dem Buch "The Dynamics and Thermodynamics of Compressible Fluid Flow" von Ascher Shapiro aus dem Jahr 1953, Seite 87, dargestellt.

Geben Sie hier die Bildbeschreibung ein

Da die X-Achse des Diagramms in der Frage der Düsenaustrittsdurchmesser ist, der die Fläche der Düsenaustrittsebene und damit das Flächenverhältnis festlegt, habe ich das Shapiro-Diagramm neu angeordnet, um das Druckverhältnis und das Geschwindigkeitsverhältnis anzuzeigen, die gegen das Flächenverhältnis aufgetragen sind. Beide Achsen sind dimensionslos, das Flächenverhältnis liegt auf der X-Achse. Ich begann mit einer Tabellenkalkulationstabelle für isentropische Strömungen von here . Die Tabelle hatte kein Geschwindigkeitsverhältnis, also habe ich es aus dem Druckverhältnis mit Gleichung 3-26 in Sutton (4. Ausgabe) berechnet.

Geben Sie hier die Bildbeschreibung ein

Wenn wir das Flächenverhältnis kennen, erhalten wir diese beiden Verhältnisse, die uns alle Informationen liefern, die wir benötigen, um die Parameter der Austrittsebene zu berechnen, die für die Schubgleichung benötigt werden – wenn wir die Halseigenschaften kennen . Dazu müssen wir einige Annahmen über den Motor treffen.

Annahme 0 – Der Halsbereich des Raptors beträgt 0,6 ft 2 . Dies stammt aus alten Informationen , die den Austrittsdurchmesser mit 1,7 Metern und das Expansionsverhältnis mit 40 angeben. Ich habe keinen Hinweis auf das Expansionsverhältnis bei der neuesten Iteration des Motors gefunden, also habe ich diesen verwendet.

Annahme 1 – Der Austrittsdruck der Ebene, wo der Fragegraph des Meeresspiegels war ICH S P erreicht maximal 15 psi. Eine vernünftige Annahme, denn wenn der Druck der Austrittsebene dem Umgebungsdruck entspricht, werden die Verluste minimiert.

Annahme 2 – Der Meeresspiegelschub des Raptor beträgt 380.000 lbf ( Wikipedia )

Mit diesen Annahmen können wir den Massendurchfluss mit Gleichung 2 berechnen. Ich lese den maximalen Meeresspiegel ab ICH S P aus dem Diagramm als 334 ergibt dies eine Massendurchflussrate von 35,3 Schwall/s.

Unter Verwendung dieses Massendurchflusses können wir Gleichung 1 verwenden, um die Abgasgeschwindigkeit zu berechnen. Der Delta-Druckterm verschwindet, wenn der Austrittsebenendruck 15 beträgt, sodass die Abgasgeschwindigkeit 10.750 ft/s beträgt.

Mit diesen Annahmen kennen wir das Flächenverhältnis, die Abgasgeschwindigkeit und den Austrittsdruck an der Spitze der Meeresspiegelkurve. Dies ermöglicht uns die Berechnung der Halseigenschaften unter Verwendung der Isentropentabellen und ergibt

v T = 4898 Fuß/s

P T = 640.900 lbf/ ft² .

Wir können jetzt den spezifischen Impuls für jeden Düsenaustrittsdurchmesser durch den folgenden Prozess berechnen:

  1. Berechnen Sie die Fläche der Austrittsebene basierend auf dem Durchmesser; Teilen Sie dies durch die Halsfläche, um das Flächenverhältnis zu erhalten.
  2. Schlagen Sie die Druck- und Geschwindigkeitsverhältnisse in den Isentropentabellen basierend auf dem Flächenverhältnis nach.
  3. Berechnen Sie den Austrittsebenendruck unter Verwendung des Druckverhältnisses und des Halsdrucks.
  4. Berechnen Sie die Abgasgeschwindigkeit unter Verwendung des Geschwindigkeitsverhältnisses und der Halsgeschwindigkeit.
  5. Verwenden Sie mit Massendurchfluss, Abgasgeschwindigkeit, Austrittsebenendruck und Austrittsebenenfläche Gleichung 1, um den Schub zu berechnen.
  6. Verwenden Sie Gleichung 2, um den spezifischen Impuls zu berechnen.

Sie können die Annahmen sicherlich ändern und andere Werte erhalten. Dies zeigt jedoch, dass Sie Kurven erzeugen können, die wie die in Frage stehenden geformt sind, indem Sie diese Parameter konstant halten und die isentropischen Strömungseigenschaften ausschließlich auf der Grundlage des Flächenverhältnisses berechnen.

Dies stimmt nicht genau mit dem Diagramm in der Frage überein; Ich nehme an, der Hersteller dieses Diagramms hatte bessere / andere Informationen über den Motor.

Geben Sie hier die Bildbeschreibung ein

Betreff: der "effektive ISP"

Es soll der durchschnittliche ISP für einen Booster sein, da er vom Meeresspiegel bis zum Abschalten des Hauptmotors reicht. Ich weiß nicht genug über das Flugprofil, um den effektiven ISP wirklich zu berechnen, also sind es nur der Meeresspiegel und der Vakuum-ISP, gemittelt, zusammen mit dem Vakuum, das als doppelt so wichtig gewichtet wird.

Source-Reddit-Thread in der Frage

Auf meiner Startsim-Todo-Liste ist (in mehrfacher Hinsicht) die ISP-gegen-Zeit-Kurve aus einem der Apollo-Trägerraketenberichte zu integrieren. Es ist ein ziemlich schönes Sigmoid und symmetrisch genug, dass ich denke, dass die Gewichtungsfunktion 1: 1 sein sollte, anstatt in Richtung Vakuum zu verzerren. en.wikipedia.org/wiki/Sigmoid_function
Vielleicht ist es 2: 1, wenn es über die Geschwindigkeit integriert wird (statt über die Zeit)? Ich wäre nicht überrascht, da ISP in deltaV verwendet wird, was keinen Bezug zur Zeit hat :) [nicht sicher, ob es Sinn macht, nur ein zufälliger Gedanke]
@BlueCoder Guter Gedanke!
blöde Frage; was ist M die horizontale Achse des Diagramms? Ist es auch irgendwie normalisierter Austrittsdurchmesser oder Austrittsdruck oder irgendein Verhältnis? Ein Großteil der Antwort befasst sich mit der Handlung in der Frage (und ich bin im Allgemeinen leicht verwirrt / abgelenkt).
Hey, die sind wirklich klasse! :D Während ich in Excel noch ziemlich eingerostet bin, werde ich mein Bestes geben und versuchen, zu sehen, ob ich die analytischen Gleichungen im verlinkten Dokument finden kann.
@uhoh Ich habe dies komplett neu geschrieben, um das Diagramm in der Frage besser mit einigen angenommenen Raptor-Eigenschaften abzugleichen. Ich habe auch versucht, den Berechnungsprozess explizit durchzugehen. Dies ist wahrscheinlich mein letzter Versuch, diese Frage zu beantworten; Ich habe es wirklich genossen, darin zu graben, danke!
Dies ist perfekt! Jetzt wird Kaffee gebrüht ... (Mir hat es aus rein ästhetischen Gründen besser gefallen, als der Plot auf 2,5 hinausging)
Ja, die von mir heruntergeladene komprimierbare Durchflusstabelle ging nur bis Mach 5. Als ich also die realistischeren Annahmen eingab, rannte ich am Ende der Tabelle vorbei.
@uhoh Ich habe darauf geantwortet, aber kurz gesagt, der Begriff des ankommenden Luftstroms ist nur für Düsentriebwerke wichtig, die diesen ankommenden Luftstrom beschleunigen. Der Druckschubterm sollte nicht weggelassen werden, er wird nur zu Null ausgewertet, wenn die Drücke übereinstimmen.
@uhoh danke für die Annahme, aber hauptsächlich für die Frage. Ich hatte mehr Spaß daran, diese Frage zu beantworten, als jede Stack-Tausch-Frage seit langem.
Ich war überrascht zu sehen, dass dieses Zeug verständlich ist, wenn man die richtige Antwort gibt. Danke!