Mecánica Celeste de Tres Cuerpos y Aerotermodinámica de Reentrada: El Perfil de Misión de Retorno Libre de Artemis 2

1. Arquitectura y Perfil Operativo de la Misión Artemis 2

La misión tripulada Artemis 2 está configurada como una expedición circunlunar de 10 días de duración. A diferencia de las misiones Apolo posteriores, que encendían el motor SPS (Service Propulsion System) para frenar e insertarse en órbita lunar circular baja (LOI, Lunar Orbit Insertion), Artemis 2 opera bajo una restricción de redundancia de seguridad máxima de soporte vital: la nave no puede quedar atrapada en el campo lunar si el motor del Módulo de Servicio Europeo (ESM) falla.

La misión utiliza una Trayectoria de Retorno Libre Híbrida (Free Return Trajectory):

1. Despegue con el cohete SLS (Space Launch System) Bloque 1 desde la plataforma 39B en el Centro Espacial Kennedy.
2. Inserción preliminar en una Órbita Terrestre Baja (LEO) de $185 \times 185\text{ km}$.
3. Encendido de elevación de perigeo con la etapa criogénica ICPS (Interim Cryogenic Propulsion Stage) hacia una Órbita Terrestre Alta Elíptica (HEO) de $42 \text{ horas}$ ($185 \times 74.000\text{ km}$) para verificar todos los subsistemas de control ambiental y soporte de vida (ECLSS) antes de comprometer a la tripulación hacia el espacio profundo.
4. Maniobra de Inyección Trans-Lunar (TLI, Trans-Lunar Injection): un incremento vectorial de velocidad $\Delta v_{\text{TLI}} \approx 3.150\text{ m/s}$ que estira el apogeo hasta los $400.000\text{ km}$.

---

2. El Problema Circular Restringido de Tres Cuerpos (CR3BP)

En la región de influencia compartida entre la Tierra y la Luna, el modelo kepleriano de dos cuerpos con cónicas simples de Newton-Kepler no es válido. La nave espacial Orion (de masa $m \ll M_{\text{Luna}} \ll M_\oplus$) se desplaza en el campo gravitatorio combinado de dos cuerpos primarios masivos que orbitan en torno a su baricentro común.

2.1 Ecuaciones de Movimiento en el Sistema Sinódico Adimensional

Definimos un sistema de coordenadas no inercial baricéntrico rotante $(x, y, z)$ donde el origen es el centro de masas del sistema Tierra-Luna, el eje $x$ apunta constantemente de la Tierra hacia la Luna, el plano $xy$ coincide con el plano orbital lunar, y el sistema rota con la velocidad angular del sistema selenocéntrico $\boldsymbol{\omega} = \omega \hat{\mathbf{z}}$.

Normalizamos las unidades de forma canónica:

• Distancia Tierra-Luna media: $D = R_\oplus + D_L = 1$.
• Masa total del sistema: $M_\oplus + M_L = 1$.
• Constante de gravitación universal y velocidad angular orbital: $G = 1$, $\omega = 1$.
• Período orbital normalizado: $T = 2\pi$.

El parámetro de masa reducido del sistema Tierra-Luna se define como:

$\mu = \frac{M_L}{M_\oplus + M_L} \approx \frac{7.342 \times 10^{22}}{5.972 \times 10^{24} + 7.342 \times 10^{22}} \approx 0.0121505856$

Las posiciones baricéntricas adimensionales de la Tierra ($P_1$) y la Luna ($P_2$) son: $P_1 = (-\mu, 0, 0), \quad P_2 = (1 - \mu, 0, 0)$

Las distancias de la nave espacial a cada uno de los primarios son: $r_1 = \sqrt{(x + \mu)^2 + y^2 + z^2}, \quad r_2 = \sqrt{(x - (1 - \mu))^2 + y^2 + z^2}$

2.2 Derivación de las Ecuaciones de Lagrange del CR3BP

El potencial gravitatorio más el pseudopotencial centrífugo rotante define el Potencial Efectivo de Jacobi $\Omega(x,y,z)$:

$\Omega(x,y,z) = \frac{1}{2}(x^2 + y^2) + \frac{1 - \mu}{r_1} + \frac{\mu}{r_2}$

Las ecuaciones diferenciales de movimiento de la nave Orion en el sistema rotante son:

$\ddot{x} - 2\dot{y} = \frac{\partial \Omega}{\partial x}$

$\ddot{y} + 2\dot{x} = \frac{\partial \Omega}{\partial y}$

$\ddot{z} = \frac{\partial \Omega}{\partial z}$

donde los términos $2\dot{y}$ y $-2\dot{x}$ son las aceleraciones aparentes de Coriolis debidas a la rotación del marco sinódico.

---

3. La Integral de Jacobi y Superficies de Velocidad Cero

Multiplicando las ecuaciones de movimiento por $\dot{x}, \dot{y}, \dot{z}$ respectivamente y sumándolas:

$\dot{x}\ddot{x} + \dot{y}\ddot{y} + \dot{z}\ddot{z} = \frac{\partial \Omega}{\partial x}\dot{x} + \frac{\partial \Omega}{\partial y}\dot{y} + \frac{\partial \Omega}{\partial z}\dot{z} = \frac{d\Omega}{dt}$

Integrando directamente con respecto al tiempo:

$\frac{1}{2} (\dot{x}^2 + \dot{y}^2 + \dot{z}^2) = \Omega(x,y,z) - \frac{1}{2} C_J$

donde $C_J$ es la única constante escalar conservada del movimiento en el CR3BP, conocida como la Integral de Jacobi:

$C_J = 2\Omega(x,y,z) - v^2 = (x^2 + y^2) + 2\left( \frac{1-\mu}{r_1} + \frac{\mu}{r_2} \right) - (\dot{x}^2 + \dot{y}^2 + \dot{z}^2)$

3.1 Superficies de Velocidad Cero y Regiones Prohibidas de Hill

Haciendo la velocidad relativa nula ($v = 0$):

$2\Omega(x,y,z) = C_J$

Esta ecuación define superficies cerradas tridimensionales en el espacio. Dado que la energía cinética real no puede ser negativa ($v^2 \ge 0$), las regiones donde $2\Omega(x,y,z) < C_J$ son mecánicamente inaccesibles para la nave espacial.

• Para valores muy altos de $C_J$, las superficies de Hill encierran a la Tierra y a la Luna en lóbulos cerrados separados.
• Para el valor crítico de $C_J$ correspondiente al punto de Lagrange $L_1$, se abre el cuello de botella que conecta energéticamente el dominio terrestre con el dominio lunar, permitiendo la transferencia orbital de Artemis 2.

---

4. Astrodinámica de la Trayectoria de Retorno Libre (Free-Return)

Una trayectoria de retorno libre circunlunar aprovecha la perturbación gravitatoria lunar para curvar la elipse heliocéntrica y reconvertirla en una trayectoria de regreso sin consumo de combustible.

4.1 Geometría Circumlunar en Forma de Ocho (Figure-8)

En el sistema sinódico, la trayectoria describe una curva asimétrica en forma de ocho:

1. Paso por delante del limbo de avance lunar: La Luna se mueve en su órbita con una velocidad orbital media de $v_L \approx 1.022\text{ km/s}$ en sentido antihorario.
2. Al ingresar en la Esfera de Influencia de Laplace de la Luna ($R_{\text{SOI}} \approx D (M_L/M_\oplus)^{2/5} \approx 66.100\text{ km}$), la nave Orion pasa por delante de la trayectoria orbital de la Luna.
3. El campo gravitatorio lunar ejerce una atracción que frena la energía geocéntrica de la nave respecto a la Tierra, curvando fuertemente el vector de velocidad hiperbólico alrededor del hemisferio no visible (lado oculto) a una altitud de periapsis selenocéntrico de:

$h_p \approx 1.400 - 1.800 \text{ km}$

4. Al girar alrededor de la Luna por el lado opuesto, el vector de salida resultante redirige el apogeo directamente hacia el perigeo atmosférico en la Tierra ($h_{\text{perigeo}} \approx 50 - 65\text{ km}$), garantizando la reentrada pasiva aunque fallen todos los ordenadores y motores de propulsión de servicio.

---

5. Vector de Asistencia Gravitatoria Selenocéntrica

En el marco de referencia inercial centrado en la Luna, el paso de Orion es una hipérbola kepleriana pura con energía específica orbital positiva:

$\mathcal{E} = \frac{v_\infty^2}{2} = \frac{v^2}{2} - \frac{\mu_L}{r} > 0$

• La velocidad asintótica de aproximación al entrar en la esfera de influencia es $\mathbf{v}_\infty^-$.
• La velocidad asintótica de salida es $\mathbf{v}_\infty^+$.
• Por conservación de energía en un campo estático conservativo selenocéntrico:

$|\mathbf{v}_\infty^-| = |\mathbf{v}_\infty^+| = v_\infty$

5.1 Ángulo de Deflexión Hiperbólica $\delta$

El cambio de dirección del vector de velocidad viene dado exactamente por la relación hiperbólica:

$\sin\left( \frac{\delta}{2} \right) = \frac{1}{e} = \frac{1}{1 + \frac{r_p v_\infty^2}{\mu_L}}$

donde $r_p = R_L + h_p$ es el radio del periapsis selenocéntrico y $e$ es la excentricidad hiperbólica.

El cambio neto vectorial en la velocidad geocéntrica inercial de la nave $\Delta \mathbf{V}$ tras abandonar el encuentro con la Luna es:

$\Delta \mathbf{V} = \mathbf{v}_\infty^+ - \mathbf{v}_\infty^- \implies |\Delta \mathbf{V}| = 2 v_\infty \sin\left( \frac{\delta}{2} \right) = \frac{2 v_\infty}{1 + \frac{r_p v_\infty^2}{\mu_L}}$

Este $\Delta \mathbf{V}$ natural de aproximadamente $1.100\text{ m/s}$ es proporcionado íntegramente por la masa de la Luna a coste cero de propelente químico.

---

6. Aerotermodinámica de Reentrada Hipersónica (Mach 32)

Al retornar de la Luna en una órbita hiperbólica/elíptica de energía cercana a cero, la nave espacial alcanza la interfaz de entrada atmosférica (definida por la NASA en $h_{\text{EI}} = 121.920\text{ m}$ o $400.000\text{ pies}$) con una velocidad aterradora:

$v_{\text{reentrada}} = \sqrt{\frac{2GM_\oplus}{R_\oplus + h_{\text{EI}}}} \approx 11.02 \text{ km/s} \approx 39.670 \text{ km/h} \quad (M \approx 32.4)$

A esta velocidad colosal, la energía cinética específica por unidad de masa que debe disiparse térmicamente es astronómica:

$e_k = \frac{1}{2} v^2 = \frac{1}{2} (11.020)^2 \approx 60.72 \text{ MJ/kg}$

Esta densidad de energía supera con creces el calor de combustión del TNT ($4.184\text{ MJ/kg}$). Si incluso una fracción del $1\%$ de esta energía penetrara hacia la cabina presurizada de la tripulación, la estructura de aleación de aluminio-litio se vaporizaría instantáneamente.

6.1 El Estrecho Corredor de Reentrada Atmosférica

El ángulo de trayectoria de vuelo en la interfaz de entrada ($\gamma_{\text{EI}}$) debe ser controlado con una precisión de centésimas de grado:

$\gamma_{\text{EI}} \in [-7.4^\circ, -6.0^\circ]$

1. Límite de Entrada Empinada ($\gamma < -7.4^\circ$): El vehículo penetra demasiado rápido en las capas atmosféricas densas. Las cargas de deceleración superan los $15 - 20 g$, aplastando fisiológicamente a los astronautas, y el flujo térmico supera el límite de falla térmica del escudo de ablación.
2. Límite de Rebote Somero (Skip-Out, $\gamma > -6.0^\circ$): La atmósfera enrarecida genera sustentación aerodinámica nula o ascendente insuficiente para capturar la cápsula. La nave rebota como una piedra en el agua hacia una órbita altamente elíptica con un período de días, donde el soporte de oxígeno y baterías se agotaría antes de un segundo reintento.

---

7. Física del Choque Disociativo y Ablación de Avcoat

Delante de la superficie convexa del escudo térmico de Orion (un cuerpo romo con radio de nariz $R_n \approx 4.6\text{ m}$), se forma una onda de choque desprendida (bow shock wave) hiperbólica.

A través del choque normal, la conservación de masa, momento y energía (relaciones de Rankine-Hugoniot) desacelera el gas de Mach 32 a régimen subsónico en la capa de choque, convirtiendo la energía cinética macroscópica en energía interna térmica:

$T_{\text{post-choque}} \approx \frac{2(\gamma - 1)}{(\gamma + 1)^2} M^2 T_\infty \sim 11.000 - 13.000 \text{ Kelvin}$

A estas temperaturas extremas (superiores a la fotosfera del Sol), el aire ya no se comporta como un gas ideal diatómico: las moléculas de nitrógeno ($N_2$) y oxígeno ($O_2$) sufren disociación molecular endotérmica completa e ionización parcial:

$N_2 + \Delta E_1 \to 2N, \quad O_2 + \Delta E_2 \to 2O, \quad N + O \to NO^+ + e^-$

7.1 Flujo de Calor Radiativo y Convectivo (Ecuación de Fay-Riddell)

El flujo de calor en el punto de estancamiento $q_s$ se compone de dos mecanismos acoplados:

$q_s = q_{\text{convectivo}} + q_{\text{radiativo}}$

El flujo convectivo laminar se modela mediante la correlación clásica de Fay-Riddell:

$q_{\text{conv}} = 0.763 \, Pr^{-0.6} (\rho_w \mu_w)^{0.1} (\rho_s \mu_s)^{0.4} \left( \frac{du_e}{dx} \right)_s^{0.5} (h_s - h_w) \left[ 1 + (Le^{0.63} - 1) \frac{h_D}{h_s} \right]$

donde $Le = \frac{\rho D_{12} c_p}{k}$ es el número de Lewis que cuantifica la recombinación química en la superficie de la pared. En primera aproximación fenomenológica:

$q_{\text{conv}} \propto \sqrt{\frac{\rho_\infty}{R_n}} v_\infty^3$

A velocidades lunares ($> 10.5\text{ km/s}$), el flujo radiativo emitido por el plasma brillante (transiciones electrónicas libres-ligadas y ligadas-ligadas) se vuelve dominante, escalando con una potencia extremadamente agresiva de la velocidad:

$q_{\text{rad}} \propto R_n^a \rho_\infty^b v_\infty^8$

En el punto álgido de reentrada de Artemis 2, el flujo de calor integrado alcanza $1.200 - 1.400 \text{ W/cm}^2$.

7.2 Termodinámica del Escudo de Ablación Monolítico Avcoat

Para absorber semejante infierno térmico, el escudo utiliza Avcoat (resina epoxi rellena de microesferas de fibra de vidrio y resina fenólica):

1. Pirólisis Endotérmica: La resina fenólica absorbe energía descomponiéndose químicamente en carbón y gases orgánicos volátiles ligeros ($H_2, CH_4, CO$).
2. Rechazo por Soplado de Capa Límite (Blowing Effect): La inyección masiva de gases fríos desde los poros del escudo hacia la capa límite empuja el choque caliente lejos de la nave, reduciendo el gradiente térmico conductivo en un $60\%$.
3. Sublimación y Desprendimiento: La capa carbonizada se desgasta milímetro a milímetro de forma controlada, garantizando que la estructura de titanio interior no supere jamás los $150^\circ\text{C}$.

---

8. Referencias Bibliográficas y Documentos Técnicos NASA

1. Szebehely, V. (1967). Theory of Orbits: The Restricted Problem of Three Bodies. Academic Press, New York. ISBN: 978-0126830507.
2. NASA Technical Report. (2020). Artemis II Mission Concept of Operations and Trajectory Profile. National Aeronautics and Space Administration, Lyndon B. Johnson Space Center. NASA-TM-2020-220451.
3. Fay, J. A., & Riddell, F. R. (1958). Theory of Stagnation Point Heat Transfer in Dissociated Air. Journal of the Aeronautical Sciences, 25(2), 73-85. DOI: 10.2514/8.7517.
4. Tauber, M. E., Bowles, J. V., & Yang, L. (1991). Atmospheric Trajectories and Aerodynamic Heating for Lunar Return Vehicles. Journal of Spacecraft and Rockets, 28(6), 640-647. DOI: 10.2514/3.26297.
5. Anderson, J. D. (2006). Hypersonic and High-Temperature Gas Dynamics (2nd Edition). American Institute of Aeronautics and Astronautics (AIAA). ISBN: 978-1563477805.
6. Battin, R. H. (1999). An Introduction to the Mathematics and Methods of Astrodynamics. AIAA Education Series. ISBN: 978-1563473425.