Resumen Técnico: El Problema de la Masa Inicial Máxima en el Diseño de Trayectorias de Bajo Empuje

Definición del Problema

El artículo aborda un vacío en la teoría del control óptimo de bajo empuje: la falta de una formalización explícita y un tratamiento analítico del problema de la Masa Inicial Máxima (MIM, por sus siglas en inglés). Mientras que el campo está dominado por las formulaciones de mínimo propelente (máxima masa final) y de transferencia de tiempo mínimo, el problema MIM plantea una pregunta distinta: para una geometría de transferencia fija, un estado final fijo y un tiempo de transferencia fijo, ¿cuál es la mayor masa inicial desde la cual una nave espacial de bajo empuje puede realizar con éxito la transferencia?

Conceptualmente, esto invierte el paradigma estándar de "dada una masa inicial, minimizar el tiempo/propelente" por el de "dado un tiempo de vuelo, maximizar la masa inicial". Aunque esta perspectiva se ha utilizado informalmente en competiciones (por ejemplo, GTOCs) y en estudios de viabilidad de minería de asteroides, no se ha establecido rigurosamente como una formulación de control óptimo independiente con sus propias condiciones necesarias y propiedades estructurales.

Metodología

Los autores emplean el Principio de Mínimo de Pontryagin (PMP) para derivar las condiciones necesarias de primer orden para el problema MIM y compararlas estructuralmente con el problema clásico de tiempo mínimo.

1. Formulación:

   • Dinámica: El sistema se modela utilizando la dinámica estándar de bajo empuje donde el vector de estado $\mathbf{x}$ (por ejemplo, elementos equinocciales modificados) y la masa $m$ evolucionan bajo un acelerador (throttle) $u \in [0,1]$ y una dirección de empuje $\hat{\mathbf{i}}$.
   • Objetivo: Maximizar la masa inicial $m(t_0)$ sujeto a un tiempo de transferencia fijo $t_f - t_0$ y condiciones de contorno $\mathbf{x}(t_0) = \mathbf{x}_0, \mathbf{x}(t_f) = \mathbf{x}_f$. Esto se convierte en la minimización de $J = -m(t_0)$.
   • Análisis del Hamiltoniano: Los autores derivan el Hamiltoniano y su función de conmutación (switching function). Un hallazgo clave es que el acelerador óptimo es siempre $u^*=1$ (empuje total) en el caso no degenerado, y la dirección del empuje sigue una ley de vector primario ($\hat{\mathbf{i}}^* \propto -\mathbf{B}(\mathbf{x})^T\boldsymbol{\lambda}$). Esto resulta en una estructura de control continua, evitando los perfiles de "bang-bang" discontinuos que se ven a menudo en otras formulaciones, lo que facilita la convergencia rápida tanto para métodos numéricos directos como indirectos.

2. Transversalidad y Condiciones de Calibración (Gauge Conditions):

   • El artículo reexamina críticamente la condición de transversalidad para problemas de tiempo libre. Argumenta que la suposición común de que el Hamiltoniano debe ser cero en el tiempo final ($H(t_f) + \boldsymbol{\lambda}(t_f)^T\dot{\mathbf{x}}_f(t_f) = 0$) para problemas de tiempo mínimo es una elección de calibración (gauge choice) más que una condición necesaria derivada del PMP para óptimos de contorno. La verdadera condición es una desigualdad de un solo lado ($H(t_f) + \boldsymbol{\lambda}(t_f)^T\dot{\mathbf{x}}_f(t_f) \le 0$).
   • Este conocimiento permite la normalización del vector de coestado $\lambda_m(t_0) = 1$ sin violar las condiciones necesarias.

3. Correspondencia Extremal:

   • Utilizando las condiciones derivadas, los autores demuestran el Teorema 1: Existe una correspondencia biunívoca entre los extremales del problema de tiempo mínimo (con masa inicial fija) y el problema de masa inicial máxima (con tiempo de vuelo fijo). Específicamente, una trayectoria que es óptima en tiempo para una masa dada también es óptima en masa para el tiempo de vuelo resultante, y viceversa.

4. Aplicación de Homotopía:

   • Los autores aplican la formulación MIM al problema de tiempo mínimo de múltiples revoluciones. Construyen un camino de homotopía variando el número de revoluciones ($N_{rev}$) mientras rastrean la superficie de masa inicial máxima $M(t_f, N_{rev})$. Este enfoque se contrasta con las homotopías tradicionales de continuación de empuje, que a menudo sufren de mal condicionamiento en escenarios de múltiples revoluciones.

Resultados Clave

• Propiedades Numéricas: El problema MIM produce un problema de valor de contorno bien planteado con siete ecuaciones y siete incógnitas (para extremos fijos), evitando los problemas de indeterminación de escala comunes en los métodos de disparo indirecto (indirect shooting) de tiempo mínimo. El control óptimo es estrictamente de empuje total, lo que asegura una convergencia numérica suave.
• Clarificación Teórica: El artículo aclifica que la condición de igualdad del Hamiltoniano terminal no es una condición necesaria para la optimalidad de tiempo, sino una elección de calibración. Esto resuelve un "malentendido común" en la literatura y valida el uso de condiciones de calibración alternativas para mejorar la estabilidad numérica.
• Visualización del Paisaje de Soluciones: En una transferencia de referencia de Órbita de Transferencia Geoestacionaria (GTO) a Órbita Geoestacionaria (GEO), la homotopía MIM logró mapear el paisaje de soluciones. El método reveló una superficie suave y diferenciable de masa inicial máxima.
  • Identificó múltiples trayectorias extremales, incluyendo mínimos y máximos locales del tiempo de vuelo a lo largo de curvas de masa constante (iso-mass).
  • Localizó con éxito la solución de tiempo óptimo global (1500 kg de masa inicial) y demostró que algunas soluciones reportadas anteriormente (por ejemplo, la solución E del benchmark) corresponden a un máximo local de tiempo de vuelo, no a un mínimo.
• Eficiencia de la Homotopía: Se observó que el camino de continuación en el marco de MIM (variando $N_{rev}$ para rastrear curvas iso-masa) es suave y bien condicionado, ofreciendo una alternativa más robusta a los métodos de continuación de empuje para navegar paisajes complejos de múltiples revoluciones.

Significado y Reivindicaciones

El artículo argumenta que el problema de la Masa Inicial Máxima merece ser elevado a una formulación de control óptimo de "primera clase" en el diseño de trayectorias de bajo empuje. Su importancia radica en tres áreas principales:

1. Primitiva de Diseño Complementaria: Aunque está estructuralmente relacionada con el problema de tiempo mínimo, la formulación MIM ofrece una perspectiva distinta que fija el tiempo y libera la masa. Esto proporciona una herramienta complementaria para el diseño de misiones, particularmente para el análisis preliminar donde deben evaluarse grandes conjuntos de transferencias.
2. Rigor Analítico: Al derivar formalmente las condiciones necesarias y corregir la comprensión de las condiciones de transversalidad en problemas de tiempo óptimo, el artículo proporciona una base teórica más rigurosa para la optimización de bajo empuje.
3. Utilidad Algorítmica: La formulación MIM permite la construcción de caminos de homotopía suaves que son menos propensos al mal condicionamiento encontrado en los métodos tradicionales de continuación de empuje. Esto permite la exploración sistemática de espacios de soluciones complejos (como las transferencias de múltiples revoluciones) y la generación de grandes bases de datos de trayectorias óptimas para el modelado de sustitutos (surrogate modeling) y el análisis de alcanzabilidad.

Los autores concluyen que la adopción de la formulación MIM como una primitiva de diseño estándar puede mejorar la eficiencia y la robustez del diseño de misiones de bajo empuje, particularmente en enfoques basados en datos y análisis de misiones preliminares.