¿Aproximaciones analíticas para las formas de estos gráficos de ISP frente al diámetro de la boquilla? (Motor Raptor)

Un comentario sobre la pregunta actualmente sin respuesta . ¿Existe un catálogo autorizado para los identificadores de lanzamiento de cohetes? me llevó a esta página de Github que enlaza con r/spacexlounge en la que actualmente aparece una publicación con una trama, titulada Eficiencia de Raptor para varios tamaños de boquillas , que también se muestra a continuación.

  1. Para los dos casos límite de las condiciones bien definidas de presión ambiental: nivel del mar y vacío , ¿existen aproximaciones analíticas a estas dos curvas? Las curvas recuerdan a las parábolas, pero siendo la ordenada algo similar a 1 ( 1 / d 2 a ) 2 dónde d Puede ser diámetro y a podría estar relacionado con la presión ambiental, pero eso es solo poner con los ojos cerrados.

  2. ¿Cómo se definió (o al menos podría definirse ) la curva intermedia para un ISP efectivo ?

nota: Aprender la fuente original de la trama será interesante, pero mi pregunta es principalmente sobre las matemáticas y la física de las boquillas de los cohetes en general.

Eficiencia Raptor para varios tamaños de boquilla

NOTA: Como se señala en los comentarios a continuación, el gráfico a continuación parece ser una contribución individual en lugar de datos publicados oficialmente. No lo cites ni lo uses como premisa de preguntas de seguimiento sin tener eso en cuenta.

Supongo que "ISP efectivo" pretende representar el promedio de tiempo del ISP en el transcurso del ascenso desde el nivel del mar hasta el vacío; podría ser una ponderación ad hoc de las dos cifras de ISP, o determinada a través de una simulación de ascenso de pasos discretos.
En los comentarios en la publicación de reddit, parece que OP cometió algunos errores ingenuos sustanciales en su hoja de cálculo de eficiencia de turbobomba (no relacionada), por lo que es posible que desee un grano de sal para acompañar su trama.
@RussellBorogove gracias por eso. Estoy principalmente interesado en ecuaciones analíticas más simples que muestren este tipo de comportamiento, en lugar de abordar la veracidad de estas curvas particulares, pero es importante saber que estas son generadas individualmente, no son oficiales, así que lo anotaré en la pregunta.
Estoy pensando que esto solo muestra el efecto de cambiar la relación de expansión, pero sin saber nada como el área de la garganta, las tasas de flujo, etc., tendré que mirar y ver si todo se cae o se puede retirar. Una tarea para mañana.
@OrganicMarble Creo que si hay ecuaciones que tienen el mismo comportamiento (la curva atm invertida alcanza su punto máximo en un diámetro finito, la curva vac alcanza su punto máximo en el infinito) estaría bien, y la curva vac siempre mayor que la curva atm sería incluso mejor. No tiene que coincidir exactamente en absoluto, puede haber demasiados detalles específicos del motor para eso.

Respuestas (1)

La única ecuación necesaria para producir estas curvas es la ecuación de empuje en un par de formas diferentes.

  1. F = metro ˙ v mi + ( pag mi pag a t metro ) A mi
  2. I s pag = F / ( gramo 0 metro ˙ )

También necesita tablas o diagramas de flujo isoentrópicos para calcular los parámetros de la boquilla que cambian con el cambio de relación de área. Aquí se muestra un diagrama de flujo isoentrópico de ejemplo, del libro de 1953 "The Dynamics and Thermodynamics of Compressible Fluid Flow" de Ascher Shapiro, página 87.

ingrese la descripción de la imagen aquí

Dado que el eje X del gráfico en la pregunta es el diámetro de salida de la boquilla, que establece el área del plano de salida de la boquilla y, por lo tanto, la relación del área, reorganicé el gráfico de Shapiro para mostrar la relación de presión y la relación de velocidad frente a la relación del área. Ambos ejes son adimensionales, la relación de área está en el eje X. Empecé con una hoja de cálculo de tabla de flujo isoentrópica de aquí . La hoja de cálculo no tenía relación de velocidad, así que la calculé a partir de la relación de presión usando la ecuación 3-26 en Sutton (4.ª edición).

ingrese la descripción de la imagen aquí

Conocer la relación del área nos da estas dos relaciones que nos brindan toda la información que necesitamos para calcular los parámetros del plano de salida necesarios para la ecuación de empuje, si conocemos las propiedades de la garganta . Tendremos que hacer algunas suposiciones sobre el motor para eso.

Suposición 0 : el área de la garganta del Raptor es de 0,6 pies 2 . Esto proviene de información antigua que da un diámetro de salida de 1,7 metros y una relación de expansión de 40. No encontré ninguna referencia a la relación de expansión en la última iteración del motor, así que usé esto.

Supuesto 1 - La presión del plano de salida donde la pregunta gráfica del nivel del mar I s pag alcanza un máximo de 15 psi. Una suposición sensata porque cuando la presión del plano de salida coincide con la ambiental, las pérdidas se minimizan.

Suposición 2 : el empuje a nivel del mar del Raptor es de 380 000 lbf ( Wikipedia )

Con estas suposiciones podemos calcular el caudal másico usando la ecuación 2. Leí el nivel máximo del mar I s pag del gráfico como 334, esto da un caudal másico de 35,3 slugs/seg.

Usando ese caudal másico, podemos usar la ecuación 1 para calcular la velocidad de escape. El término de presión delta desaparece cuando la presión del plano de salida es 15, por lo que la velocidad de escape es 10 750 ft/s.

Con estas suposiciones, conocemos la relación del área, la velocidad de escape y la presión de salida en el pico de la curva del nivel del mar. Esto nos permite calcular las propiedades de la garganta usando las tablas isoentrópicas y da

V t = 4898 pies/s

PAG t = 640 900 lbf/pie 2 .

Ahora podemos calcular el impulso específico para cualquier diámetro de salida de la boquilla mediante el siguiente proceso:

  1. Calcule el área del plano de salida en función del diámetro; divida esto por el área de la garganta para obtener la relación del área.
  2. Busque las relaciones de presión y velocidad en las tablas isoentrópicas según la relación de área.
  3. Calcule la presión del plano de salida usando la relación de presión y la presión de garganta.
  4. Calcule la velocidad de escape usando la relación de velocidad y la velocidad de garganta.
  5. Con el caudal másico, la velocidad de escape, la presión del plano de salida y el área del plano de salida, utilice la Ecuación 1 para calcular el empuje.
  6. Use la Ecuación 2 para calcular el impulso específico.

Ciertamente puede cambiar las suposiciones y obtener valores diferentes. Sin embargo, esto demuestra que al mantener esos parámetros constantes y calcular las propiedades de flujo isoentrópico únicamente en función de la relación de área, puede generar curvas que tienen la forma de las de la pregunta.

Esto no coincide exactamente con el gráfico de la pregunta; Supongo que el creador de ese gráfico tenía información mejor/diferente sobre el motor.

ingrese la descripción de la imagen aquí

Re: el "ISP efectivo"

Se supone que es el ISP promedio para un propulsor a medida que va desde el nivel del mar hasta el corte del motor principal. No sé lo suficiente sobre el perfil de vuelo para calcular realmente el ISP efectivo, por lo que es solo el nivel del mar y el ISP de vacío promediados junto con el vacío ponderado como el doble de importante.

Fuente: hilo de Reddit en la pregunta.

En mi lista de tareas pendientes del simulador de lanzamiento está integrar (en múltiples sentidos) la curva ISP versus tiempo de uno de los informes del vehículo de lanzamiento Apollo. Es un sigmoide bastante agradable, y lo suficientemente simétrico como para pensar que la función de ponderación debería ser 1:1 en lugar de inclinarse hacia el vacío. en.wikipedia.org/wiki/Sigmoid_function
¿Quizás es 2: 1 si se integra sobre la velocidad (en lugar del tiempo)? No me sorprendería ya que ISP se usa en deltaV, que no hace referencia al tiempo :) [no estoy seguro si tiene sentido, solo un pensamiento aleatorio]
@BlueCoder ¡Buen pensamiento!
pregunta tonta; qué es METRO el eje horizontal de la trama? ¿También está normalizado el diámetro de salida de alguna manera, o la presión de salida, o alguna relación? Gran parte de la respuesta analiza la trama en la pregunta (y me confundo/distraigo fácilmente en general).
Hey estos son realmente excelentes! :D Aunque todavía estoy bastante oxidado en Excel , haré lo mejor que pueda e intentaré ver si puedo encontrar las ecuaciones analíticas en el documento vinculado.
@uhoh, reescribí esto por completo para que coincida mejor con el gráfico en la pregunta usando algunas propiedades asumidas de Raptor; También traté de pasar explícitamente por el proceso de cálculo. Este es probablemente mi último intento en esta pregunta; Realmente disfruté profundizar en él, ¡gracias!
¡Esto es perfecto! Ya se está haciendo café... (me gustó más cuando la trama se bajó a 2.5 por cuestiones puramente estéticas)
Sí, la tabla de flujo comprimible que descargué solo salió a Mach 5. Entonces, cuando puse las suposiciones más realistas, me salí del final de la tabla.
@uhoh lo respondí, pero brevemente, el término de flujo de aire entrante solo es importante para los motores a reacción, que aceleran ese flujo de aire entrante. El término de empuje de presión no debe descartarse, simplemente se evalúa como cero si las presiones coinciden.
@uhoh gracias por aceptar, pero sobre todo por la pregunta. Me divertí más respondiendo esta que cualquier pregunta de intercambio de pila en mucho tiempo.
Me sorprendió ver que estas cosas son comprensibles dada la respuesta correcta. ¡Gracias!