Ecuación del cohete de Tsiolkovski

La ecuación del cohete de Tsiolkovski considera el principio del cohete: un aparato que puede acelerarse a sí mismo (empuje) expulsando parte de su masa a alta velocidad en el sentido opuesto a la aceleración obtenida debido a la conservación de la cantidad de movimiento.

La ecuación lleva el nombre del científico ruso Konstantín Tsiolkovsky que, de forma independiente, la derivó y publicó en su obra de 1903.^[1]

Introducción

Un volumen de control que pierde masa

Considerando un dispositivo, un cohete, por ejemplo, que en un momento dado tiene una masa $m$ y que se desplaza hacia delante con una velocidad $v$ , en ese mismo instante, el dispositivo expele la cantidad de masa $m_{e}$ con una velocidad de flujo $v_{e}$ . El volumen de control incluirá tanto la masa $m$ del dispositivo como la masa expelida $m_{e}$ .^[2]

Durante el tiempo dt, su velocidad se incrementa de v a v+ dv puesto que una cantidad de masa $dm_{e}$ ha sido expulsada y por tanto se incrementó el escape. Este incremento de la velocidad hacia delante, sin embargo, no cambia la velocidad $v_{e}$ de la masa expelida, como lo vería un espectador fijo, puesto que la masa se mueve a una velocidad constante una vez que ha sido expulsada. Los impulsos son creados por ( $\textstyle \sum F_{cv}$ ) la cual representa la resultante de todas las fuerzas externas, como resistencia al avance y peso, que actúan en el volumen de control en la dirección del movimiento.^[2]

Esta resultante de fuerzas no incluye la fuerza que impulsa al volumen de control hacia delante, puesto que esta fuerza(llamada empuje) es interna al volumen de control; es decir, el empuje actúa con magnitud igual pero dirección opuesta en la masa m del dispositivo y la masa expelida $m_{e}$ .

Al aplicar el principio de impulso y cantidad de movimiento al volumen de control tenemos:^[2]

\sum F_{cv}=-v{\frac {\partial {m_{e}}}{\partial {t}}}+m{\frac {\partial {v}}{\partial {t}}}-{\frac {\partial {m_{e}}}{\partial {t}}}dv-{\frac {\partial {v_{e}}}{\partial {t}}}dm_{e}

La velocidad del dispositivo puede escribirse como:^[3]

\sum F_{cv}=m{\dfrac {dv}{dt}}-v_{D/e}{\dfrac {dm_{e}}{dt}}

Aquí el término $dm_{e}/dt$ representa la tasa a la cual se expulsara la masa.

Expresión de Tsiolkovski

La expresión de Tsiolkovski expresa que para cualquier maniobra o viaje que incluya maniobras:

$\Delta v=v_{e}\ln {\frac {m_{0}}{m_{1}}}$

o equivalentemente:

$m_{1}=m_{0}e^{-\Delta v\ /v_{e}},\qquad m_{0}=m_{1}e^{\Delta v\ /v_{e}}$

donde:

$m_{0}\,$ es la masa total inicial.
$m_{1}\,$ la masa total final
$v_{e}\,$ la velocidad de los gases de salida con respecto al cohete (impulso específico).

Por otro lado el término:

$1-{\frac {m_{1}}{m_{0}}}=1-e^{-\Delta v\ /v_{e}}$

es la fracción de masa (la parte de la masa total inicial que se utiliza para propulsar el cohete).

$\Delta v$ (delta-v) es el resultado de integrar en el tiempo la aceleración producida por el uso del motor del cohete (no la aceleración debida a otras fuentes como rozamiento o gravedad). En el caso típico de aceleración en el sentido de la velocidad, es el incremento de la velocidad. En el caso de aceleración en el sentido contrario (desaceleración) es el decremento de la velocidad. La gravedad y el rozamiento cambian también la velocidad pero no forman parte de delta-v. Por ello, delta-v no es simplemente el cambio en la velocidad. Sin embargo, el empuje se aplica en corto tiempo, y durante ese periodo las otras fuentes de aceleración pueden ser despreciables, así que la delta-v de un momento determinado puede aproximarse al cambio de velocidad. La delta-v total puede ser simplemente añadida, aunque entre momentos de propulsión la magnitud y cantidad de velocidad cambia debido a la gravedad, como por ejemplo en una órbita elíptica.

La ecuación se obtiene integrando la ecuación de conservación del momento lineal.

$mdv=v_{e}dm$

para un cohete simple que emite masa a velocidad constante (la masa que se emite es $dm$ ).

Aunque es una simplificación extrema, la ecuación del cohete muestra lo esencial de la física del vuelo del cohete en una única y corta ecuación. La magnitud delta-v es una de las cantidades más importantes en mecánica orbital que cuantifica lo difícil que es cambiar de una trayectoria a otra.

Claramente, para conseguir un Δv elevado, debe ser $m_{0}$ elevada (crece exponencialmente con delta-v), o $m_{1}$ debe ser pequeña, o $v$ debe ser elevada, o una combinación de estos.

En la práctica, esto se consigue con cohetes muy grandes (aumentando $m_{0}$ ), con varias etapas (decrementando $m_{1}$ ), y cohetes con combustibles con velocidades de escape muy elevadas. Los cohetes Saturno V utilizados en el Proyecto Apolo y los motores de iones usados en sondas no tripuladas de larga distancia son un buen ejemplo de esto.

La ecuación del cohete muestra un «decaimiento exponencial» de masa, pero no como función del tiempo, sino conforme a mientras se produce la Δv. La Δv que corresponde a la «vida media» es $v_{e}\ln 2\approx 0,693v_{e}$ .

Esta ecuación había sido derivada antes por el matemático británico William Moore en 1813.

Etapas

En el caso de cohetes de varias etapas, la ecuación se aplica a cada etapa, y en cada etapa la masa inicial del cohete es la masa total del cohete después de dejar la etapa anterior y la masa final es la del cohete justo antes de dejar la etapa que se está calculando. El impulso específico para cada etapa puede ser diferente. Por ejemplo, si el 80% de la masa es el combustible de la primera etapa y el 10% es masa en vacío de la primera etapa y el 10% es el resto del cohete, entonces:

$\Delta v=v_{e}\ln 5=1,61v_{e}$

Con tres etapas similares más pequeñas, se tiene:

$\Delta v=3v_{e}\ln 5=4,83v_{e}$

y la carga útil es un 0,1% de la masa inicial.

Un cohete de una etapa a órbita, también con un 0,1% de carga útil puede tener una masa del 11% para depósitos y motores y el 88,9% de combustible. Esto da

$\Delta v=v_{e}\ln(100/11,1)=2,20v_{e}$

Si el motor de una nueva etapa se enciende antes de que la etapa anterior haya caído y los motores que trabajan simultáneamente tienen un impulso específico diferente (como es muchas veces el caso en cohetes de combustible sólido y etapas líquidas), la situación es más complicada.

Energía

En el caso ideal $m_{1}$ es la carga útil y $m_{0}-m_{1}$ es la masa que reacciona (que corresponde a depósitos vacíos sin masa, etc.). La energía necesaria es:

${\frac {1}{2}}(m_{0}-m_{1})v_{e}^{2}$ .

Esta es la energía cinética de la masa de reacción y no la energía cinética requerida por la carga, pero si $v_{e}$ =10 km/s y la velocidad del cohete es 3 km/s, entonces la velocidad de la masa de reacción solo cambia desde 3 a 7 km/s; La energía «ahorrada» corresponde al incremento de la energía cinética específica (energía cinética por kg) para el cohete. En general:

$d\left({\frac {1}{2}}v^{2}\right)=vdv=vv_{e}dm/m={\frac {1}{2}}\left(v_{e}^{2}-(v-v_{e})^{2}+v^{2}\right)dm/m$

Se tiene:

$\Delta \epsilon =\int v\,d(\Delta v)$

donde $\epsilon$ es la energía específica del cohete y $\Delta v$ es una variable separada, no solo el cambio en $v$ . En el caso de usar el cohete parar decelerar, es decir, expeler masa de reacción en la dirección de la velocidad, $v$ es negativa.

La fórmula es para el caso ideal sin pérdidas de energía por calor, etc. Esta última causa una reducción del empuje, así que es una desventaja aun cuando el objetivo es perder energía (decelererar).

Si la energía se produce por la masa misma, como en un cohete químico, el valor del combustible tiene que ser: ${\frac {1}{2}}v_{e}^{2}$ , donde para el valor del combustible se tiene que tomar también la masa del oxidante. Un valor típico es v_e=4,5 km/s, correspondiente a 10,1 MJ/kg. La valor real es más alto pero parte de la energía se pierde en forma de calor que sale como radiación. La energía necesaria es:

$E={\frac {1}{2}}m_{1}\left(e^{\Delta v\ /v_{e}}-1\right)v_{e}^{2}$

Conclusiones:

Para $\Delta v\ll v_{e}$ se tiene $E\approx {\frac {1}{2}}m_{1}v_{e}\Delta v$
Para una $\Delta v$ dada, la energía mínima se necesita si $v_{e}=0,6275\Delta v$ , requiriendo una energía de

E=0,772m_{1}(\Delta v)^{2}

.

Empezando desde velocidad cero es el 54,4% más que la energía cinética de la carga útil. Empezando desde una velocidad que no es cero, la energía requerida puede ser "menos" que el incremento de energía cinética de la carga. Este puede ser el caso cuando la masa de reacción tiene una velocidad menor después de ser expelida que antes. Por ejemplo, desde una órbita baja de 300 km de altitud a una órbita de escape es un incremento de 29,8 MJ/kg, lo cual, usando un impulso específico de 4,5 km/s, tiene un coste neto de 20,6 MJ/kg (

\Delta v

= 3,20 km/s; las energías son por kg de carga útil).

Esta optimización no tiene en cuenta las masa de los diferentes tipos de cohetes.

Además, para un objetivo determinado, como por ejemplo cambiar de una órbita a otra, la $\Delta v$ requerida dependa mucho de la velocidad a la que el motor produce $\Delta v$ y determinadas maniobras pueden ser imposibles si esta es muy baja. Por ejemplo, un lanzamiento a LEO requiere normalmente una $\Delta v$ de alrededor de 9,5 km/s (mayormente para conseguir la velocidad), pero si el motor pudiese producir $\Delta v$ a una velocidad solo algo más elevada que g, sería un lanzamiento lento y requeriría una $\Delta v$ mucho más elevada (costaría una $\Delta v$ de 9,8 m/s cada segundo). Si la aceleración posible es $g$ o menor, no es posible ir a órbita con ese motor.

La potencia se obtiene de

$P={\frac {1}{2}}mav_{e}={\frac {1}{2}}Fv_{e}$

donde $F$ es el empuje y $a$ es la aceleración debida a ella. Por ello, el empuje teórico posible por unidad de potencia es 2 dividido por el impulso específico en m/s. La eficiencia de empuje es el empuje real entre empuje teórico.

Si se usa energía solar se restringe $a$ ; en el caso de $v_{e}$ elevadas, la aceleración posible es inversamente proporcional a la velocidad de escape, así que el tiempo necesario para conseguir una Δv es proporcional a $v_{e}$ ; con el 100% de eficiencia:

para $\Delta v\ll v_{e}$ tenemos que $t\approx {\frac {mv_{e}\Delta v}{2P}}$

Ejemplos:

potencia 1000 W, masa 100 kg, $\Delta v$ = 5 km/s, $v_{e}$ = 16 km/s, lleva 1,5 meses.
potencia 1000 W, masa 100 kg, $\Delta v$ = 5 km/s, $v_{e}$ = 50 km/s, lleva 5 meses.

Por ello, la $v_{e}$ no puede ser demasiado alta.

Ejemplos

Se asume un impulso específico de 4,5 km/s y una $\Delta v$ de 9,7 km/s (Tierra a LEO).

Un cohete de una etapa a órbita: $1-e^{-9,7/4,5}$ = 0,884, por ello el 88,4% de la masa total inicial será propelente. El restante 11,6% es para los motores, el tanque y la carga.

Un cohete de dos etapas a órbita: se supone que la primera etapa da una $\Delta v$ de 5,0 km/s; $1-e^{-5,0/4,5}$ = 0,671, por ello, el 67,1%. El restante es el 32,9%. Después de dejar la primera etapa, la masa será este 32,9% menos el tanque y el motor de la primera etapa. Si se asume que esto es el 8% de la masa total inicial, queda el 24,9%. La segunda etapa da una $\Delta v$ de 4,7 km/s; $1-e^{-4,7/4,5}$ = 0,648, por ello, el 64,8% de la masa restante debe ser propelente, que es el 16,2%, y el 8,7% el tanque, el motor y la carga de la segunda etapa, Así que hay disponible el 16,7% para motores, tanques y carga útil.

Véase también

Referencias

↑ К. Э. Циолковский, Исследование мировых пространств реактивными приборами, 1903. It is available online here Archivado el 15 de agosto de 2011 en Wayback Machine. in a RARed PDF
↑ ^a ^b ^c Hibbeler, R. C. (2010). Cruz Castillo, Luis Miguel, ed. Dinámica. México: Pearson. p. 282. ISBN 978607442560-4.
↑ Hibbeler, R. C. (2010). Cruz Castillo, Luis Miguel, ed. Dinámica. México: Pearson. p. 283. ISBN 978607442560-4.

Datos: Q834199

[1] К. Э. Циолковский, Исследование мировых пространств реактивными приборами, 1903. It is available online here Archivado el 15 de agosto de 2011 en Wayback Machine. in a RARed PDF

[ReferenceA-2] Hibbeler, R. C. (2010). Cruz Castillo, Luis Miguel, ed. Dinámica. México: Pearson. p. 282. ISBN 978607442560-4.

[3] Hibbeler, R. C. (2010). Cruz Castillo, Luis Miguel, ed. Dinámica. México: Pearson. p. 283. ISBN 978607442560-4.

[1]

[2]

[3]