技術要約：低推力軌道設計における最大初期質量問題

問題の定義

本論文は、低推力最適制御理論における空白、すなわち**最大初期質量（Maximum Initial Mass: MIM）**問題の明示的な定式化および解析的扱いの欠如に対処している。この分野は、最小推進剤（最大最終質量）や最小時間転移の定式化が支配的であるが、MIM問題は異なる問いを投げかける。すなわち、固定された転移幾何学、固定された最終状態、および固定された転移時間に対して、低推力宇宙機がその転移を成功させるために可能な最大の初期質量はいくらであるか、という問いである。

概念的には、これは標準的な「与えられた初期質量に対して、時間または推進剤を最小化する」というパラダイムを、「与えられた飛行時間に対して、初期質量を最大化する」というものへと反転させている。この視点は、コンペティション（例：GTOC）や小惑星採掘の実現可能性調査において非公式に利用されてきたが、独自の必要条件や構造的特性を持つ独立した最適制御定式化として厳密に確立されてはいない。

手法

著者らは、**ポントリャーギンの最小原理（PMP）**を用いて、MIM問題の一次必要条件を導出し、それを古典的な最小時間問題と構造的に比較している。

1. 定式化:

   • 動力学: システムは、スロットル $u \in [0,1]$ と推力方向 $\hat{\mathbf{i}}$ の下で進化する状態ベクトル $\mathbf{x}$（例：修正正準要素）および質量 $m$ を用いた標準的な低推力動力学を用いてモデル化される。
   • 目的関数: 固定された転移時間 $t_f - t_0$ および境界条件 $\mathbf{x}(t_0) = \mathbf{x}_0, \mathbf{x}(t_f) = \mathbf{x}_f$ の下で、初期質量 $m(t_0)$ を最大化する。これは $J = -m(t_0)$ の最小化に変換される。
   • ハミルトニアン解析: 著者らはハミルトニアンと関連するスイッチング関数を導出している。主要な知見は、非退化的なケースにおいて最適なスロットルは常に $u^*=1$（全推力）であり、推力方向はプライマー・ベクトル則（$\hat{\mathbf{i}}^* \propto -\mathbf{B}(\mathbf{x})^T\boldsymbol{\lambda}$）に従うことである。これにより、連続的な制御構造が得られ、他の定式化で見られる不連続なバンバン制御のプロファイルが回避されるため、直接法および間接法の両方において急速な収束が促進される。

2. 横断条件とゲージ条件:

   • 著者らは、自由時間の問題における横断条件を批判的に再検討している。最小時間問題において、ハミルトニアンが最終時刻でゼロになる必要がある（$H(t_f) + \boldsymbol{\lambda}(t_f)^T\dot{\mathbf{x}}_f(t_f) = 0$）という一般的な仮定は、PMPから導かれる必要条件ではなく、一つのゲージ選択であると論じている。真の条件は、片側不等式（$H(t_f) + \boldsymbol{\lambda}(t_f)^T\dot{\mathbf{x}}_f(t_f) \le 0$）である。
   • この洞察により、必要条件を破ることなく余剰コストの正規化 $\lambda_m(t_0) = 1$ を行うことが可能となる。

3. 極値の対応関係:

   • 導出された条件を用いて、著者らは定理1を証明している。すなわち、最小時間問題の極値（固定された初期質量を伴う）と、最大初期質量問題（固定された転移時間を伴う）の間には一対一の対応関係が存在する。具体的には、ある質量に対して時間最適である軌道は、結果として得られる飛行時間に対して質量最適でもある、という相互関係である。

4. ホモトピーの適用:

   • 著者らは、MIM定式化を多回転最小時間問題に適用している。彼らは、回転数（$N_{rev}$）を変化させながら最大初期質量の曲面 $M(t_f, N_{rev})$ を追跡することで、ホモトピー経路を構築している。このアプローチは、多回転シナリオにおいてしばしば不良条件（ill-conditioning）に陥る従来の推力継続ホモトピーと対比されている。

主要な結果

• 数値的特性: MIM問題は、7つの方程式と7つの未知数（固定端点の場合）を持つ、良好に定義された境界値問題をもたらし、間接射撃法（indirect shooting methods）で一般的なスケーリングの不定性の問題を回避する。最適制御は厳密に全推力であり、滑らかな数値的収束を保証する。
• 理論的明確化: 本論文は、終端ハミルトニアンの等式条件が時間最適性のための必要条件ではなく、正規化の選択であることを明確にしている。これにより、文献における「共通の誤解」を解消し、数値的安定性を向上させるための代替ゲージ条件の使用を正当化している。
• 解のランドスケープの可視化: 地球静止遷移軌道（GTO）から静止軌道（GEO）へのベンチマーク転移において、MIMホモトピーは解のランドスケープの写像に成功した。この手法は以下のことを明らかにした。
  • 局所的な最小値および等質量曲線に沿った飛行時間の極大値を含む、複数の極値軌道を特定した。
  • 以前に報告された解（ベンチマークの解Eなど）が、最小時間ではなく、飛行時間の局所的な極大値に対応していることを示した。
• ホモトピーの効率性: MIMフレームワークにおける継続パス（$N_{rev}$ を変化させて等質量曲線を追跡する手法）は、滑らかで良好な条件を持つことが観察された。これは、複雑な多回転ランドスケープをナビゲートするための、より堅牢な代替手段を提供している。

意義と主張

本論文は、最大初期質量問題が、低推力軌道設計における「第一級」の最適制御定式化へと昇格されるべきであると主張している。その意義は、主に以下の3点にある：

1. 補完的な設計プリミティブ: 構造的には最小時間問題と関連しているものの、MIM定式化は、時間を固定して質量を解放するという異なる視点を提供する。これは、大量の転移を評価する必要がある予備解析において、ミッション設計のための補完的なツールとなる。
2. 解析的厳密性: 低推力最適化の必要条件を正式に導出し、時間最適問題における横断条件の理解を修正することにより、本論文は低推力最適化に対してより厳密な理論的基礎を提供している。
3. アルゴリズム的有用性: MIM定式化は、従来の推力継続法で見られる不良条件に陥りにくい、滑らかなホモトピー経路の構築を可能にする。これにより、複雑な解空間（多回転転移など）の体系的な探索や、サロゲートモデリングおよび到達可能性解析のための大規模な最適軌道データベースの生成が可能となる。

著者らは、MIM定式化を標準的な設計プリミティブとして採用することが、特にデータ駆動型のアプローチや予備的なミッション解析において、低推力ミッション設計の効率性と堅牢性を向上させることができると結論付けている。