Resumo Técnico: O Problema da Massa Inicial Máxima no Design de Trajetórias de Baixo Empuxo

Definição do Problema

O artigo aborda uma lacuna na teoria de controle ótimo de baixo empuxo: a falta de uma formalização explícita e tratamento analítico do problema da Massa Inicial Máxima (MIM - Maximum Initial Mass). Enquanto o campo é dominado por formulações de transferência de propelente mínimo (massa final máxima) e de tempo mínimo, o problema MIM faz uma pergunta distinta: para uma geometria de transferência fixa, estado final fixo e tempo de transferência fixo, qual é a maior massa inicial da qual uma espaçonave de baixo empuxo pode realizar a transferência com sucesso?

Conceitualmente, isso inverte o paradigma padrão de "dada a massa inicial, minimizar tempo/propelente" para "dado o tempo de voo, maximizar a massa inicial". Embora essa perspectiva tenha sido utilizada informalmente em competições (ex: GTOCs) e para estudos de viabilidade de mineração de asteroides, ela não foi rigorosamente estabelecida como uma formulação de controle ótimo independente, com suas próprias condições necessárias e propriedades estruturais.

Metodologia

Os autores empregam o Princípio do Mínimo de Pontryagin (PMP) para derivar as condições necessárias de primeira ordem para o problema MIM e compará-las estruturalmente com o problema clássico de tempo mínimo.

1. Formulação:

   • Dinâmica: O sistema é modelado usando a dinâmica padrão de baixo empuxo, onde o vetor de estado $\mathbf{x}$ (ex: elementos equinociais modificados) e a massa $m$ evoluem sob um acelerador (throttle) $u \in [0,1]$ e direção de empuxo $\hat{\mathbf{i}}$.
   • Objetivo: Maximizar a massa inicial $m(t_0)$ sujeito a um tempo de transferência fixo $t_f - t_0$ e condições de contorno $\mathbf{x}(t_0) = \mathbf{x}_0, \mathbf{x}(t_f) = \mathbf{x}_f$. Isso é convertido em uma minimização de $J = -m(t_0)$.
   • Análise do Hamiltoniano: Os autores derivam o Hamiltoniano e a função de comutação associada. Uma descoberta fundamental é que o acelerador ótimo é sempre $u^*=1$ (empuxo total) no caso não degenerado, e a direção do empuxo segue uma lei de vetor primário ($\hat{\mathbf{i}}^* \propto -\mathbf{B}(\mathbf{x})^T\boldsymbol{\lambda}$). Isso resulta em uma estrutura de controle contínua, evitando os perfis bang-bang descontínuos frequentemente vistos em outras formulações, o que facilita a convergência rápida para métodos numéricos diretos e indiretos.

2. Transversalidade e Condições de Gauge:

   • O artigo reexamina criticamente a condição de transversalidade para problemas de tempo livre. Argumenta-se que a suposição comum de que o Hamiltoniano deve ser zero no tempo final ($H(t_f) + \boldsymbol{\lambda}(t_f)^T\dot{\mathbf{x}}_f(t_f) = 0$) para problemas de tempo mínimo é uma escolha de gauge e não uma condição necessária derivada do PMP para ótimos de contorno. A verdadeira condição é uma desigualdade de um lado só ($H(t_f) + \boldsymbol{\lambda}(t_f)^T\dot{\mathbf{x}}_f(t_f) \le 0$).
   • Este insight permite a normalização do vetor de coestado $\lambda_m(t_0) = 1$ sem violar as condições necessárias.

3. Correspondência Extremal:

   • Usando as condições derivadas, os autores provam o Teorema 1: Existe uma correspondência biunívoca entre os extremais do problema de tempo mínimo (com massa inicial fixa) e o problema da massa inicial máxima (com tempo de transferência fixo). Especificamente, uma trajetória que é ótima em tempo para uma determinada massa também é ótima em massa para o tempo de voo resultante, e vice-versa.

4. Aplicação de Homotopia:

   • Os autores aplicam a formulação MIM ao problema de tempo mínimo de múltiplas revoluções. Eles constroem um caminho de homotopia variando o número de revoluções ($N_{rev}$) enquanto rastreiam a superfície de massa inicial máxima $M(t_f, N_{rev})$. Esta abordagem é contrastada com as tradicionais homotopias de continuação de empuxo, que frequentemente sofrem de má condicionamento em cenários de múltiplas revoluções.

Principais Resultados

• Propriedades Numéricas: O problema MIM gera um problema de valor de contorno bem posto com sete equações e sete incógnitas (para endpoints fixos), evitando os problemas de indeterminação de escala comuns em métodos de tiro indireto de tempo mínimo. O controle ótimo é estritamente de empuxo total, garantindo uma convergência numérica suave.
• Clarificação Teórica: O artigo esclarece que a condição de igualdade do Hamiltoniano terminal não é uma condição necessária para a otimalidade de tempo, mas sim uma escolha de normalização. Isso resolve um "equívoco comum" na literatura e valida o uso de condições de gauge alternativas para melhorar a estabilidade numérica.
• Visualização do Cenário de Solução: Em uma transferência de referência de Órbita de Transferência Geoestacionária (GTO) para Órbita Geoestacionária (GEO), a homotopia MIM mapeou com sucesso o cenário de solução. O método revelou uma superfície suave e diferenciável de massa inicial máxima.
  • Identificou múltiplas trajetórias extremais, incluindo mínimos e máximos locais de tempo de voo ao longo de curvas de massa iso.
  • Localizou com sucesso a solução global de tempo ótimo (1500 kg de massa inicial) e demonstrou que algumas soluções previamente relatadas (ex: solução E do benchmark) correspondem a máximos locais de tempo de voo, não mínimos.
• Eficiência de Homotopia: O caminho de continuação na estrutura MIM (variando $N_{rev}$ para rastrear curvas de iso-massa) foi observado como sendo suave e bem condicionado, ofereando uma alternativa mais robusta aos métodos de continuação de empuxo para navegar em cenários complexos de múltiplas revoluções.

Significância e Alegações

O artigo argumenta que o problema da Massa Inicial Máxima merece ser elevado a uma formulação de controle ótimo de "primeira classe" no design de trajetórias de baixo empuxo. Sua significância reside em três áreas principais:

1. Primitiva de Design Complementar: Embora estruturalmente relacionada ao problema de tempo mínimo, a formulação MIM oferece uma perspectiva distinta que fixa o tempo e libera a massa. Isso fornece uma ferramenta complementar para o design de missões, particularmente para análise preliminar onde grandes conjuntos de transferências devem ser avaliados.
2. Rigor Analítico: Ao derivar formalmente as condições necessárias e corrigir a compreensão das condições de transversalidade em problemas de tempo ótimo, o artigo fornece uma base teórica mais rigorosa para a otimização de baixo empuxo.
3. Utilidade Algorítmica: A formulação MIM permite a construção de caminhos de homotopia suaves que são menos propensos ao mau condicionamento encontrado nos métodos tradicionais de continuação de empuxo. Isso permite a exploração sistemática de espaços de solução complexos (como transferências de múltiplas revoluções) e a geração de grandes bancos de dados de trajetórias ótimas para modelagem de substitutos (surrogate modeling) e análise de alcançabilidade.

Os autores concluem que a adoção da formulação MIM como uma primitiva de design padrão pode melhorar a eficiência e a robustez do design de missões de baixo empuxo, particularmente em abordagens baseadas em dados e análise de missão preliminar.