Improved Directional State Transition Tensors for Accurate Aerocapture Performance Analysis

Resumen Técnico: Tensores de Transición de Estado Direccional Mejorados para un Análisis Preciso del Rendimiento de la Captura Aerodinámica

Enunciado del Problema
La captura aerodinámica es una tecnología crítica para las misiones a planetas exteriores, ya que ofrece ahorros en la masa de propelente en comparación con la inserción orbital propulsiva. Sin embargo, la dinámica de la captura aerodinámica está dominada por fuerzas aerodinámicas no conservativas que varían significativamente a lo largo de la trayectoria, generando no linealidades inherentes. Los métodos tradicionales de propagación de incertidumbre lineal, como el Análisis de Covarianza Lineal (LCA), no logran capturar con precisión estos efectos no lineales, especialmente cuando las incertidumbres de entrada son grandes o se utilizan modelos atmosféricos complejos. Aunque existen métodos de orden superior como la Expansión del Caos Polinomial (PCE) y la Transformada Unscented (UT), estos adolecen de la maldición de la dimensionalidad o de altos costos computacionales que dificultan su implementación a bordo.

Los Tensores de Transición de Estado (STT), que representan expansiones de series de Taylor de orden superior del flujo de la solución, ofrecen una alternativa semi-analítica computacionalmente eficiente para escenarios orbitales. No obstante, su aplicación a la captura aerodinámica ha estado limitada por el crecimiento exponencial en el número de ecuaciones diferenciales requeridas para calcular STT de orden superior a medida que aumenta la dimensión del estado. Los Tensores de Transición de Estado Direccional (DSTT) se desarrollaron para mitigar este costo proyectando el estado en una base de dimensión reducida a lo largo de una dirección dinámicamente importante. Trabajos anteriores establecieron que, para la dinámica orbital, la dirección del estiramiento lineal máximo (el vector propio dominante del Tensor de Cauchy-Green de segundo orden, o CGT) es relativamente constante y se alinea con las direcciones de estiramiento de orden superior. Este artículo investiga si estas hipótesis se mantienen para la captura aerodinámica y busca mejorar el rendimiento de los DSTT identificando bases de reducción más efectivas.

Metodología
Los autores analizan el problema de la captura aerodinámica utilizando un vector de estado de 7 grados de libertad que incluye la distancia radial, la latitud, la longitud, la velocidad, el ángulo de trayectoria de vuelo, el ángulo de rumbo y el logaritmo natural de la densidad atmosférica. El estudio emplea dos escenarios basados en una misión propuesta de sonda orbital de Urano: una entrada de referencia y una entrada de alta velocidad.

La metodología central implica tres pasos principales:

1. Descomposición de la Dinámica: Los autores descomponen la dinámica de la captura aerodinámica en componentes conservativos (gravitatorios) y disipativos (aerodinámicos). Al analizar los CGT de estos componentes por separado, demuestran que la dirección del estiramiento lineal máximo varía con el tiempo en la captura aerodinámica, a diferencia de la dinámica orbital, debido al dominio de las fuerzas aerodinámicas durante el pico de presión dinámica.
2. Identificación de una Base Novel: Reconociendo que la dirección del estiramiento lineal (CGT de segundo orden) no se alinea con las direcciones de estiramiento de orden superior en la captura aerodinámica, los autores proponen el uso de Tensores de Cauchy-Green de Orden Superior (HOCGT). Introducen dos formulaciones aumentadas:
   • CGT Selectivos (sCGT): Estos identifican la dirección del estiramiento máximo para un subconjunto específico de estados (por ejemplo, posición, velocidad, ángulo de trayectoria de vuelo).
   • CGT de Magnitud de Interés (qCGT): Estos identifican la dirección del estiramiento máximo para una función no lineal específica del estado, como la energía orbital específica o el radio del apoapsis.
     Los vectores propios dominantes de estos tensores se calculan utilizando el Método de Potencia de Orden Superior Simétrico Desplazado (SS-HOPM).
3. Construcción y Validación de DSTT: Los autores construyen DSTT de tercer orden utilizando estas nuevas direcciones (hoDSTT, sDSTT, qDSTT) y los comparan con los DSTT tradicionales construidos utilizando los vectores propios del CGT de segundo orden (1-DSTT, 3-DSTT, 6-DSTT) y los STT completos. El rendimiento se evalúa utilizando el error de la norma de Frobenius normalizada, la precisión de la propagación de perturbaciones y la precisión de la predicción para el radio del apoapsis terminal y la energía específica.

Contribuciones Clave

• Identificación de Estiramiento Variable en el Tiempo: El artículo demuestra que, para la captura aerodinámica, la dirección del estiramiento lineal máximo no es constante y no se alinea con las direcciones de estiramiento de orden superior, invalidando las hipótesis previas derivadas de la dinámica orbital.
• Formulaciones Aumentadas de HOCGT: El desarrollo de sCGT y qCGT permite la identificación de direcciones de sensibilidad adaptadas a subconjuntos específicos del estado o a magnitudes de interés específicas, en lugar de solo a la norma del estado completo.
• Mayor Precisión de los DSTT: El estudio muestra que construir DSTT a lo largo de direcciones derivadas de HOCGT de tercer y cuarto orden (y sus variantes aumentadas) reduce significativamente el error de la norma de Frobenius en comparación con los DSTT tradicionales, particularmente durante la región de presión dinámica máxima donde las no linealidades son más fuertes.
• Análisis de Compensaciones: Los autores proporcionan un análisis detallado que muestra que el uso de una única dimensión latente a lo largo de una dirección sensible de orden superior puede ofrecer una mejor precisión que el uso de múltiples dimensiones (por ejemplo, 6-DSTT) basadas únicamente en el estiramiento lineal.

Resultados

• Precisión de la Aproximación: Los DSTT construidos utilizando las nuevas direcciones de orden superior (hoDSTT, sDSTT, qDSTT) exhiben errores de norma de Frobenius normalizada más bajos que los DSTT tradicionales durante la fase crítica de presión dinámica máxima.
• Propagación de Perturbaciones: Para perturbaciones alineadas con la dirección de estiramiento no lineal máximo, los nuevos DSTT mantienen una alta precisión. Incluso para perturbaciones ortogonales a esta dirección, el hoDSTT de una sola dimensión funciona de manera comparable o mejor que la Matriz de Transición de Estado (STM).
• Predicción de Magnitudes de Interés:
  • Energía Específica: El uso de sCGT (seleccionando estados constituyentes) o qCGT (utilizando la función de energía) para construir DSTT redujo el error mediano de predicción de energía en aproximadamente medio orden de magnitud en comparación con los métodos tradicionales.
  • Radio del Apoapsis: El enfoque qDSTT que utiliza directamente la función del radio del apoapsis funcionó mal debido a problemas numéricos derivados de la discontinuidad en el semieje mayor cerca del límite de captura. Sin embargo, el enfoque sDSTT, que selecciona los estados constituyentes ($r, V, \gamma$) utilizados en el cálculo del apoapsis, superó a otros métodos DSTT en la predicción del radio del apoapsis.
• Eficiencia Computacional: Los nuevos DSTT requieren significativamente menos términos para integrar y propagar que los STT completos (por ejemplo, 14 términos para un DSTT de una sola dimensión de tercer orden frente a 343 términos para un STT de segundo orden), ofreciendo una alternativa computacionalmente eficiente para la guía y navegación a bordo.

Significado
El artículo establece que los STT y DSTT son herramientas prácticas para el análisis del rendimiento de la captura aerodinámica, siempre que la base de reducción se elija cuidadosamente. Al ir más allá de la hipótesis de que las direcciones de estiramiento lineal son suficientes, los autores demuestran que aprovechar la dinámica de orden superior y adaptar la base a magnitudes de interés específicas conduce a una propagación de incertidumbre más robusta y precisa. Estos resultados permiten el uso de métodos semi-analíticos eficientes para la guía y navegación estocástica a bordo en entornos altamente no lineales y no conservativos, abordando una brecha crítica en las capacidades actuales de diseño de misiones de captura aerodinámica. La metodología propuesta es generalizable a cualquier sistema dinámico no lineal y no conservativo.