[論文レビュー] Emulating satellite drag from large simulation experiments

Furong Sun, Robert B. Gramacy+3 more|arXiv (Cornell University)|Dec 1, 2017

Advanced Multi-Objective Optimization Algorithms参考文献 28被引用数 4

ひとこと要約

本稿では、大規模なシミュレーションからの衛星の空気抵抗の正確でリアルタイムなエミュレーションを可能にする、マルチレベルの局所的近似ガウス過程（laGP）フレームワークを提案する。マルチスケールのグローバル／ローカルモデリングとセットワイズの局所サブセット選択を統合することで、70,000コアアワーのシミュレーションを用いて、ハッブル宇宙望遠鏡およびGRACE衛星のデータにおいて1%未塔の相対誤差を達成し、標準的なGPおよびlaGP手法を上回る性能を発揮した。

ABSTRACT

Obtaining accurate estimates of satellite drag coefficients in low Earth orbit is a crucial component in positioning and collision avoidance. Simulators can produce accurate estimates, but their computational expense is much too large for real-time application. A pilot study showed that Gaussian process (GP) surrogate models could accurately emulate simulations. However, cubic runtime for training GPs means that they could only be applied to a narrow range of input configurations to achieve the desired level of accuracy. In this paper we show how extensions to the local approximate Gaussian Process (laGP) method allow accurate full-scale emulation. The new methodological contributions, which involve a multi-level global/local modeling approach, and a set-wise approach to local subset selection, are shown to perform well in benchmark and synthetic data settings. We conclude by demonstrating that our method achieves the desired level of accuracy, besting simpler viable (i.e., computationally tractable) global and local modeling approaches, when trained on seventy thousand core hours of drag simulations for two real-world satellites: the Hubble space telescope (HST) and the gravity recovery and climate experiment (GRACE).

研究の動機と目的

• ガウス過程（GP）エミュレーションの計算上の非現実性を解消する。これは、立方体時間計算量と二次的メモリスケーリングに起因する。
• 入力空間全体にわたる抵抗係数表面の非定常性が、標準的なGPエミュレーションの精度を損なうという課題を克服する。
• 低地球軌道における衝突回避および軌道位置決めのためのリアルタイムかつ高精度な抵抗係数予測を実現する。
• 実際の衛星形状（HSTおよびGRACE）から得られる数百万回のシミュレーション実行を処理できる、スケーラブルで並列処理可能なエミュレーションフレームワークを開発する。
• マルチレベルモデリングと最適化された局所サブセット選択を導入することで、先行するlaGPおよびグローバルGP手法を改善し、より高い予測精度を実現する。

提案手法

• サブサンプリングされたデータからのグローバルGPの長さスケール推定（Liu & Hung, 2015）と局所的近似ガウス過程（laGP）のフィッティングを組み合わせるマルチレベルのグローバル／ローカルモデリング戦略を採用し、非定常予測の安定化を図る。
• 高次元入力空間における近傍選択の効率性と予測精度を向上させるために、セットワイズの局所サブセット選択手法を導入する。
• 非定常な抵抗表面へのより良い局所フィットを可能にするために、laGPに非等方的局所共分散構造を拡張する。
• 活性学習を用い、ALC（平均対数予測密度）基準により、情報量の多い局所トレーニングポイントを動的に選択することで、局所モデルの忠実度を向上させる。
• HSTおよびGRACEのシミュレーションにわたり、数百万回の実行を対象として、スーパーコンputingアーキテクチャを活用した大規模並列処理により、70,000コアアワーのシミュレーションデータを生成する。
• 軌道に基づく局所エミュレーションを統合し、衛星軌道に沿った点ごとの予測を可能にすることで、リアルタイムナビゲーションに不可欠なサポートを提供する。

実験結果

リサーチクエスチョン

• RQ1 大規模なシミュレーションアンサンブルから、正確に衛星の空気抵抗係数を予測できるスケーラブルで非定常なエミュレーションフレームワークを開発できるか？
• RQ2 局所的近似ガウス過程（laGP）手法をどのように改善すれば、軌道抵抗シミュレーションに一般的な高次元で非定常な入力空間においても高い精度を維持できるか？
• RQ3 マルチスケールのグローバル／ローカルモデリング戦略は、標準的なグローバルまたはローカルGP手法と比較して、どの程度エミュレーション精度を向上させるか？
• RQ4 セットワイズの局所サブセット選択と非等方的局所フィッティングは、大規模な抵抗シミュレーションエミュレーションにおいて予測誤差を顕著に低減できるか？
• RQ5 提案手法は、低地球軌道における衝突回避などのリアルタイム応用に十分に計算的に実行可能で、精度も十分に高いと評価できるか？

主な発見

• 提案されたマルチレベルlaGPフレームワークは、実世界のハッブル宇宙望遠鏡およびGRACE衛星のシミュレーションにおいて、1%未塔の平均二乗平方誤差（RMSPE）を達成し、精度の基準基準を満たした。
• 非定常性および局所構造のより良い取り扱いのおかげで、数百万回のシミュレーション実行でトレーニングされたにもかかわらず、標準的なグローバルGPおよびベースラインlaGP手法を上回る性能を発揮した。
• セットワイズの局所サブセット選択と非等方的フィッティングの拡張により、従来のlaGPにおける近傍選択と比較して、予測精度が顕著に向上した。
• マルチスケールのグローバル／ローカルモデリング戦略により、laGPの非定常挙動が安定化され、多様な入力設定において一貫した性能を発揮できるようになった。
• 固定された局所近傍サイズ（n=50）を用いることで、大規模なデータセット（過去のGP手法の10倍以上）に対しても、計算コストがほぼ一定（壁時計時間）を維持するため、スケーラビリティが確保された。
• 本手法により、衛星軌道に沿った点ごとの抵抗係数予測がリアルタイムで可能となり、衝突確率評価およびナビゲーションシステムにとって不可欠な要件を満たした。