Astrodinámica de Lanzamiento y la Mecánica del Giro Gravitacional (Gravity Turn): De la Ecuación de Tsiolkovsky a la Inserción Orbital

1. La Paradoja Orbital: Altitud vs Velocidad Tangencial

Uno de los mayores equívocos de la divulgación generalista es asumir que un cohete gasta cientos de toneladas de combustible para "escapar" de la gravedad terrestre hacia arriba.

En la órbita de la Estación Espacial Internacional (ISS), a una altitud media de $h = 410\text{ km}$ sobre el nivel del mar, la intensidad del campo gravitatorio de la Tierra $g(h)$ apenas se reduce un $11\%$ respecto a la superficie:

$g(h) = g_0 \left( \frac{R_\oplus}{R_\oplus + h} \right)^2 = 9.81 \left( \frac{6371}{6371 + 410} \right)^2 \approx 8.68 \text{ m/s}^2$

Los astronautas en órbita no flotan porque "no haya gravedad", sino porque están en un estado permanente de caída libre continua. El cohete no necesita subir eternamente; necesita desplazarse de forma horizontal tan rápido que la curvatura de su caída balística coincida exactamente con la curvatura esférica de la superficie de la Tierra.

1.1 Velocidad Circular Orbital en LEO

Igualando la fuerza gravitatoria newtoniana con la aceleración centrípeta requerida para una órbita circular a un radio $r = R_\oplus + h$:

$\frac{G M_\oplus m}{r^2} = \frac{m v_{\text{orb}}^2}{r} \implies v_{\text{orb}} = \sqrt{\frac{G M_\oplus}{r}} = \sqrt{\frac{\mu_\oplus}{R_\oplus + h}}$

Para una órbita baja a $h = 250\text{ km}$ ($r = 6.621 \times 10^6\text{ m}$):

$v_{\text{orb}} = \sqrt{\frac{3.986004 \times 10^{14} \text{ m}^3\text{s}^{-2}}{6.621 \times 10^6\text{ m}}} \approx 7759 \text{ m/s} \approx 27.930 \text{ km/h} \quad (\sim \text{Mach } 23)$

Si un cohete volara estrictamente vertical y apagara sus motores a $400\text{ km}$ de altitud con velocidad horizontal nula ($v_x = 0$), simplemente caería de regreso hacia el océano en una trayectoria puramente balística suborbital en menos de 9 minutos.

---

2. Derivación Rigurosa de la Ecuación del Cohete de Tsiolkovsky

La propulsión por reacción espacial se rige por la conservación del momento lineal en un sistema de masa variable. Consideremos un cohete de masa instantánea $m(t)$ desplazándose a una velocidad $v(t)$ en el vacío del espacio intergaláctico. En un intervalo diferencial $dt$, el motor eyecta una masa de gas $dm > 0$ a una velocidad relativa constante de escape $v_e$ en sentido opuesto.

2.1 Balance de Momento Lineal

En el instante inicial $t$: $p(t) = m v$

En el instante $t + dt$: $p(t + dt) = (m - dm)(v + dv) + dm (v - v_e)$

Despreciando el término de segundo orden $dm \cdot dv$: $p(t + dt) = m v + m dv - v dm + v dm - v_e dm = m v + m dv - v_e dm$

Aplicando el principio de conservación del momento lineal en ausencia de fuerzas externas ($dp = 0$): $p(t + dt) - p(t) = 0 \implies m dv - v_e dm = 0$

Dado que el consumo de masa del cohete es $\Delta m = -dm$ (la masa total del cohete disminuye, $dm < 0$): $m dv = -v_e dm \implies dv = -v_e \frac{dm}{m}$

Integrando entre el estado inicial ($m_0, v_0$) y el estado final de agotamiento de combustible ($m_f, v_f$):

$\int_{v_0}^{v_f} dv = -v_e \int_{m_0}^{m_f} \frac{dm}{m} \implies \Delta v = v_f - v_0 = v_e \ln\left( \frac{m_0}{m_f} \right)$

Definiendo la velocidad de escape efectiva en términos del Impulso Específico ($I_{sp}$) a nivel del mar o en el vacío ($v_e = I_{sp} g_0$):

$\Delta v = I_{sp} g_0 \ln\left( \frac{m_0}{m_f} \right)$

Esta ecuación exponencial impone la "tiranía de la ecuación del cohete": duplicar el $\Delta v$ disponible exige elevar la fracción de masa de combustible de forma exponencial al cuadrado, lo que convierte cada metro por segundo $(\text{m/s})$ desperdiciado en pérdidas gravitatorias o aerodinámicas en decenas de toneladas adicionales de estructura en la plataforma de lanzamiento.

---

3. Ecuaciones Diferenciales de Vuelo en Coordenadas Intrínsecas

Durante la fase de ascenso atmosférico de un vehículo de lanzamiento, el cohete está sometido a tres fuerzas vectoriales simultáneas:

1. Fuerza de Empuje propulsivo ($\mathbf{T}$): Vector colineal con el eje longitudinal del vehículo.
2. Fuerza de Gravedad ($\mathbf{W} = m\mathbf{g}$): Vector dirigido hacia el centro de masa de la Tierra.
3. Fuerza Aerodinámica ($\mathbf{A} = \mathbf{L} + \mathbf{D}$): Sustentación perpendicular al vector de velocidad y arrastre opuesto a la velocidad relativa atmosférica.

Definiendo el sistema de coordenadas intrínseco formado por el vector tangente unitario $\hat{\mathbf{t}}$ (en la dirección de la velocidad instantánea $\mathbf{v}$) y el vector normal unitario $\hat{\mathbf{n}}$:

• $v$: Módulo de la velocidad de trayectoria.
• $\gamma$: Ángulo de trayectoria de vuelo respecto a la horizontal local ($\gamma = 90^\circ$ en despegue vertical, $\gamma = 0^\circ$ en vuelo horizontal).
• $\alpha$: Ángulo de ataque aerodinámico (ángulo entre el eje del empuje y el vector velocidad).
• $r = R_\oplus + h$: Distancia al centro de la Tierra.

3.1 Descomposición en los Ejes Tangencial y Normal

Proyectando la segunda ley de Newton $\sum \mathbf{F} = m \mathbf{a}$ en los ejes intrínsecos:

$\text{Eje Tangencial: } m \frac{dv}{dt} = T \cos \alpha - D - m g(r) \sin \gamma$

$\text{Eje Normal: } m v \frac{d\gamma}{dt} = T \sin \alpha + L - m g(r) \cos \gamma + \frac{m v^2}{r} \cos \gamma$

Dividiendo ambas ecuaciones por la masa instantánea $m(t)$:

$\frac{dv}{dt} = \frac{T \cos \alpha}{m} - \frac{D}{m} - g \sin \gamma$

$v \frac{d\gamma}{dt} = \frac{T \sin \alpha}{m} + \frac{L}{m} - \left( g - \frac{v^2}{r} \right) \cos \gamma$

El término $\frac{v^2}{r} \cos \gamma$ representa la aceleración centrípeta debida a la curvatura geodésica esférica de la superficie de la Tierra.

---

4. La Mecánica del Giro Gravitacional Puro (Gravity Turn)

En cohetería real, forzar al vehículo a volar con un ángulo de ataque no nulo ($\alpha \ne 0$) dentro de la atmósfera densa genera dos problemas estructurales letales:

1. Introduce momentos flectores aerodinámicos masivos (presión aerodinámica transversal $q \alpha$) que pueden partir el cohete longitudinalmente.
2. Requiere desviar las toberas mediante actuadores de guiado vectorial (TVC), introduciendo pérdidas de empuje coseno (steering losses).

La solución magistral descubierta en la ingeniería astronáutica es el Giro Gravitacional Puro (Zero-Lift Gravity Turn): diseñar la trayectoria de modo que el ángulo de ataque sea estrictamente nulo ($\alpha = 0$, $L = 0$) durante todo el vuelo a través de la atmósfera.

4.1 La Dinámica del Gravity Turn con $\alpha = 0$

Sustituyendo $\alpha = 0$ y $L = 0$ en las ecuaciones diferenciales anteriores:

$\frac{dv}{dt} = \frac{T}{m} - \frac{D}{m} - g \sin \gamma$

$\frac{d\gamma}{dt} = -\frac{g}{v} \left( 1 - \frac{v^2}{v_{\text{orb}}^2} \right) \cos \gamma$

¡Observemos con atención matemática la segunda ecuación! La tasa de cambio del ángulo de inclinación $\frac{d\gamma}{dt}$ es estrictamente negativa y puramente proporcional a la componente perpendicular de la gravedad:

$\frac{d\gamma}{dt} = -\frac{g}{v} \cos \gamma$

Esto significa que la propia gravedad de la Tierra es la que inclina y curva el cohete suavemente hacia la horizontal, sin necesidad de gastar energía de combustible en maniobras de giro forzado.

4.2 La Maniobra de Cabeceo Inicial (Pitch Kick)

Sin embargo, hay una singularidad matemática al inicio: en el momento del despegue sobre la plataforma, el vuelo es vertical ($\gamma = 90^\circ$). En ese punto:

$\cos(90^\circ) = 0 \implies \frac{d\gamma}{dt} = 0$

Si el cohete subiera perfectamente vertical a $\gamma = 90^\circ$, la gravedad actuaría exactamente a lo largo del eje longitudinal y jamás produciría torsión lateral.

Por ello, todos los cohetes (como el Falcon 9, el Saturn V o el SLS) ejecutan una maniobra crucial denominada Pitch Kick (o cabeceo inicial) pocos segundos después de abandonar la torre de lanzamiento (típicamente entre $100 \text{ m}$ y $1.000 \text{ m}$ de altitud, cuando la velocidad supera los $50 - 100\text{ m/s}$):

• El sistema de guiado inclina el vehículo apenas $1.5^\circ - 3.0^\circ$ hacia el Este.
• A partir de ese milisegundo, el ordenador de a bordo alinea el vector de empuje estrictamente paralelo al flujo de aire ($\alpha = 0$).
• La gravedad hace el resto del trabajo durante los siguientes 3 minutos: a medida que el cohete gana velocidad, el vector de trayectoria $\gamma$ se aplana de forma natural y elegante hasta rozar los $0^\circ$ en el borde del vacío espacial.

---

5. Presupuesto de Delta-V y Cuantificación de Pérdidas de Ascenso

Integrando la ecuación de aceleración tangencial a lo largo del tiempo total de combustión $t_b$:

$\int_0^{t_b} \frac{dv}{dt} dt = \int_0^{t_b} \frac{T \cos \alpha}{m} dt - \int_0^{t_b} \frac{D}{m} dt - \int_0^{t_b} g \sin \gamma dt$

El incremento neto de velocidad del vehículo $v_{\text{final}}$ se desglosa en el balance energético:

$\Delta v_{\text{ideal}} = v_{\text{final}} + \Delta v_{\text{grav}} + \Delta v_{\text{drag}} + \Delta v_{\text{steer}} + \Delta v_{\text{press}}$

Para alcanzar una órbita baja estable a $250\text{ km}$, se requiere una velocidad final de $v_{\text{orb}} \approx 7.800\text{ m/s}$. Sin embargo, el cohete debe proporcionar entre $9.200\text{ m/s}$ y $9.600\text{ m/s}$ de $\Delta v_{\text{ideal}}$.

5.1 Pérdidas por Gravedad (Gravity Losses)

Las pérdidas por gravedad provienen de quemar combustible oponiéndose al peso estático del cohete:

$\Delta v_{\text{grav}} = \int_0^{t_b} g(t) \sin \gamma(t) dt$

• Si el cohete vuela vertical ($\gamma = 90^\circ$), $\sin \gamma = 1$: el cohete sufre la máxima penalización gravitatoria posible ($9.81\text{ m/s}$ perdidos por cada segundo de encendido de motor). Para un encendido de 180 segundos, esto representaría una pérdida catastrófica de más de $1.765\text{ m/s}$.
• A medida que el cohete se curva horizontalmente hacia $\gamma \to 0^\circ$, $\sin \gamma \to 0$: las pérdidas por gravedad caen a cero absoluto.

5.2 Pérdidas por Arrastre y la Región Max-Q

El arrastre aerodinámico viene dado por:

$D = \frac{1}{2} \rho(h) v^2 S C_D$

donde la densidad atmosférica $\rho(h)$ decae exponencialmente con la altura barométrica $\rho(h) = \rho_0 e^{-h/H}$ (con altura de escala $H \approx 7.5\text{ km}$):

$\Delta v_{\text{drag}} = \int_0^{t_b} \frac{\frac{1}{2} \rho(h) v^2 S C_D}{m(t)} dt \approx 100 - 150 \text{ m/s}$

Aproximadamente a los 70–90 segundos de vuelo ($h \approx 12 - 14\text{ km}$, $M \approx 1.5 - 2.0$), el producto $\frac{1}{2}\rho v^2$ alcanza su máximo absoluto conocido como Max-Q (Máxima Presión Dinámica). En este punto, los motores reducen ligeramente su potencia (throttle down) para no superar los límites de carga estructural del fuselaje antes de atravesar las capas enrarecidas de la mesosfera.

---

6. Aporte Inercial de la Rotación Terrestre (Efecto Tirachinas)

La Tierra rota sobre su eje polar con una velocidad angular de:

$\omega_\oplus = \frac{2\pi}{86164.1\text{ s}} \approx 7.2921 \times 10^{-5} \text{ rad/s}$

Cualquier punto sobre la superficie a una latitud geográfica $\phi$ posee una velocidad lineal tangencial intrínseca hacia el Este debida a la rotación planetaria:

$v_{\text{rot}}(\phi) = \omega_\oplus R_\oplus \cos \phi$

• En el Ecuador terrestre ($phi = 0^circ$, p. ej. Puerto Espacial de Kourou en la Guayana Francesa): $v_{\text{rot}}(0^circ) = (7.2921 \times 10^{-5}) \times (6.378 \times 10^6) \approx 465.1 \text{ m/s} \approx 1674 \text{ km/h}$
• En Cabo Cañaveral ($phi = 28.5^circ$, Florida, EE.UU.): $v_{\text{rot}}(28.5^circ) = 465.1 \times \cos(28.5^circ) \approx 408.7 \text{ m/s} \approx 1471 \text{ km/h}$

Al despegar y curvar la trayectoria hacia el Este (acimut de lanzamiento $\beta \approx 90^\circ$), el cohete suma gratuitamente esos más de $400\text{ m/s}$ a su balance de velocidad orbital. Si un cohete intentara despegar hacia el Oeste en dirección retrógrada, tendría que quemar suficiente propelente para cancelar primero esa velocidad y luego generar la contraria, penalizando el presupuesto en más de $800\text{ m/s}$ y reduciendo la carga útil a una fracción insignificante.

---

7. Optimización de Masas y Etapado de Cohetes Multi-fase

Dado que el cociente de masa $m_0/m_f$ de una sola etapa está limitado por la metalurgia de materiales ligeros a valores prácticos de $\lambda = m_{\text{estructura}}/m_0 \sim 0.05 - 0.08$, ningún vehículo monietapa (SSTO, Single-Stage-To-Orbit) que queme propelentes químicos puede alcanzar la velocidad orbital de inserción en la Tierra con una carga útil económicamente viable.

Para un cohete de $N$ etapas en serie, el $\Delta v$ total es la suma acumulada de cada etapa independiente:

$\Delta v_{\text{total}} = \sum_{i=1}^N I_{sp,i} g_0 \ln\left( \frac{m_{0,i}}{m_{f,i}} \right)$

Al desprender la primera etapa vacía (cuyo peso muerto en tanques gigantescos de aleación aluminio-litio o fibra de carbono ya no aporta energía), la segunda etapa comienza su encendido en el vacío con una masa inicial drásticamente reducida. Además, las toberas de la segunda etapa están optimizadas con relaciones de expansión geométrica extremas ($epsilon = A_e / A_t > 150 - 250$), elevando el impulso específico de vacío a más de $348\text{ s}$ en motores kerolox y más de $450\text{ s}$ en motores criogénicos hidrolox.

---

8. Referencias Bibliográficas y Textos de Referencia

1. Tsiolkovsky, K. E. (1903). The Exploration of Cosmic Space by Means of Reaction Devices. The Science Survey (Nauchnoe Obozrenie), San Petersburgo.
2. Bate, R. R., Mueller, D. D., & White, J. E. (1971). Fundamentals of Astrodynamics. Dover Publications, New York. ISBN: 978-0486600611.
3. Curtis, H. D. (2020). Orbital Mechanics for Engineering Students (4th Edition). Butterworth-Heinemann / Elsevier Aerospace Engineering Series. ISBN: 978-0081021330.
4. Humble, R. W., Henry, G. N., & Larson, W. J. (1995). Space Propulsion Analysis and Design. McGraw-Hill Science/Engineering/Math. ISBN: 978-0070313200.
5. NASA SP-8008. (1965). Prelaunch Inflight and Ascent Flight Guidance and Trajectory Optimization Analysis. NASA Space Vehicle Design Criteria (Guidance and Control).
6. Griffin, M. D., & French, J. R. (2004). Space Vehicle Design (2nd Edition). American Institute of Aeronautics and Astronautics (AIAA) Education Series. DOI: 10.2514/4.862426.