スイングバイ宇宙航法（惑星重力アシストと運動量保存）

Q: スイングバイ宇宙航法（惑星重力アシストと運動量保存）とは何ですか？

太陽に対する惑星の公転速度を \mathbf{V}_p、接近前の探査機の速度を \mathbf{v}_1、離脱後の速度を \mathbf{v}_2 とする。 惑星を基準とする相対座標系では、エネルギー保存則より接近時の相対速度 u_1 と離脱時の相対速度 u_2 の大きさは等しい。 \mathbf{u}_1 = \mathbf{v}_1 - \mathbf{V}_p, \quad \mathbf{u}_2 = \mathbf{v}_2 - \mathbf{V}_p, \quad |\mathbf{u}_1| = |\mathbf{u}_2| = u 惑星の重力による軌道の偏向角を \delta とすると、太陽静止系における離脱速度 \mathbf{v}_2 はベクトル和として求まる。 \mathbf{v}_2 = \mathbf{u}_2 + \mathbf{V}_p 探査機が惑星の後方を通過して180度反転した場合（完全正面衝突的弾性衝突）、最大離脱速度は次式となる。 v_2^{\text{max}} = v_1 + 2 V_p

Q: スイングバイ宇宙航法（惑星重力アシストと運動量保存）は社会やプログラミングでどのように使われていますか？

深宇宙探査機の軌道計画ソフトウェア（SPICE / GMAT）: NASAやJAXAの軌道設計エンジニアは、ランベルト問題の解法と多段階スイングバイ探索アルゴリズムを用いて、消費デルタV（推進剤）が最小となる打ち上げウィンドウと天体フライバイ順序を最適化している。

燃料を一切使わずに、惑星の公転運動量を盗んで探査機を超加速させる天体力学の幾何学

この単元の発見をシェア

3行要約（Key Takeaways / AIOサマリー）

1: スイングバイ（重力アシスト）は、惑星の背後を通過することで惑星の公転速度ベクトルを機体に加算し、燃料ゼロで劇的に加速・方向転換する宇宙航法
2: 惑星から見た探査機の相対速度の「大きさ」は接近前と離脱後で変わらないが、惑星が太陽系内を高速移動しているため、太陽静止系で見ると最大で「惑星公転速度の2倍」の速度増加が得られる
3: ボイジャー1号・2号（太陽系脱出）、カッシーニ（土星探査）、はやぶさ（地球帰還）など、すべての深宇宙探査の命綱となる数理技術

なぜこれを学ぶのか？ — 現実社会とテクノロジーでの価値

Purpose & Motivation

なぜ生まれたのか？

地球から積んでいけるロケット燃料の重量限界を打ち破り、太陽系の遥か彼方（木星・土星・冥王星）まで探査機を飛ばすために考案された天体力学の裏ワザ。

身につく武器・思考力

天体の公転運動量保存則を味方につけ、「燃料消費ゼロ」で機体を脱出速度まで一気に加速・方向転換させるダイナミックなベクトル幾何学が手に入る。

実社会・ITでの活躍現場

•ボイジャー1号・2号による太陽系外（星間空間）への脱出航法
•小惑星探査機「はやぶさ」「はやぶさ2」の地球スイングバイによる軌道精密調整
•水星探査機「ベピ・コロンボ」の多段階減速スイングバイ（金星・地球利用）
•スペースデブリ（宇宙ゴミ）の効率的捕獲・軌道離脱ミッション設計

Layer 1 — 直感・シミュレーター

探査機の進入角度 $\theta$ や突入速度をスライダーで変えてみよう。公転する惑星の後ろ側を絶妙な距離で回り込んだ瞬間、探査機がグイッと引き込まれ、太陽系脱出速度へと急激に加速されて飛び去る様子を体感できる。

進入角 $\theta$（後方通過で加速）45°

探査機の突入初速 $v_0$2.8 km/s

惑星重力（木星・地球質量）1500

重力アシストによる加速倍率: 1.00倍に急上昇！

公転する惑星の後ろ側をかすめることで、惑星の公転運動量を「燃料消費ゼロ」で奪って脱出速度に達します

オレンジ＝公転中の木星｜青色＝探査機（ボイジャー/はやぶさ）の双曲線フライバイ軌道

直感チェック！ミニクイズ

探査機が太陽系基準で「加速スイングバイ」を行うためには、公転している惑星のどちら側を通過する必要があるか？

Layer 2 — 定義・受験・証明

定義

太陽に対する惑星の公転速度を $\mathbf{V}_p$ 、接近前の探査機の速度を $\mathbf{v}_1$ 、離脱後の速度を $\mathbf{v}_2$ とする。惑星を基準とする相対座標系では、エネルギー保存則より接近時の相対速度 $u_1$ と離脱時の相対速度 $u_2$ の大きさは等しい。

\mathbf{u}_1 = \mathbf{v}_1 - \mathbf{V}_p, \quad \mathbf{u}_2 = \mathbf{v}_2 - \mathbf{V}_p, \quad |\mathbf{u}_1| = |\mathbf{u}_2| = u

惑星の重力による軌道の偏向角を $\delta$ とすると、太陽静止系における離脱速度 $\mathbf{v}_2$ はベクトル和として求まる。

\mathbf{v}_2 = \mathbf{u}_2 + \mathbf{V}_p

探査機が惑星の後方を通過して180度反転した場合（完全正面衝突的弾性衝突）、最大離脱速度は次式となる。

v_2^{\text{max}} = v_1 + 2 V_p

公式・数理モデル

\mathbf{v}_2 = (\mathbf{v}_1 - \mathbf{V}_p)_{\text{rotated}} + \mathbf{V}_p, \quad v_2^{\text{max}} = v_1 + 2 V_p

証明・導出

完全弾性衝突の1次元アナロジーによる最大加速の導出：

質量 $M$ の巨大な惑星が速度 $V_p$ で進んでおり、質量 $m$ （ $m \ll M$ ）の探査機が逆向きから速度 $v_1$ で正面接近し、重力散乱によって180度Uターンして戻る極限を考える。運動量保存則および力学的エネルギー保存則は、質量比 $m/M \to 0$ の弾性衝突と完全に等価である。

惑星静止系で見ると、探査機の接近相対速度は $u_1 = v_1 + V_p$ である。惑星の質量が無限大と見なせるため、接近と離脱で相対スピードは変わらず、離脱相対速度は $u_2 = -(v_1 + V_p)$ （逆向き）となる。これを太陽静止系に戻すため、惑星の速度 $V_p$ を足し戻す。

v_2 = u_2 + V_p = (v_1 + V_p) + V_p = v_1 + 2 V_p

すなわち、探査機は惑星の公転速度の「2倍」の速度をまるまる手に入れることができる。惑星側の公転速度は運動量保存により $\Delta V_p = -\frac{m}{M} \Delta v$ だけ減速するが、 $M \approx 10^{27} \text{kg}$ （木星）に対し $m \approx 10^3 \text{kg}$ であるため、惑星の速度低下は $10^{-24} \text{m/s}$ （検出不可能な微量）にとどまり、実質的に無尽蔵のエネルギー源となる。

押さえておきたいポイント

ボイジャー2号は木星・土星・天王星・海王星の4惑星を連続スイングバイする「グランドツアー」により、数十年の航路短縮と打ち上げ燃料の大幅削減を達成した
水星探査機「ベピ・コロンボ」や太陽探査機「パーカー・ソーラー・プローブ」は、逆に惑星の前方をかすめる「減速スイングバイ」を何度も繰り返し、太陽の強大な引力圏の内側へと安全に降下している
軌道計算は二体問題のパッチド・コニック法（円錐曲線接続法）によって高精度に事前設計される

参照基準・学習指導要領との対応

本トピックの内容・数理モデルは、文部科学省学習指導要領（算数・数学・情報）および国際標準カリキュラムに準拠し、エンジニアリング現場の実装事例と照らし合わせてMath World Lab編集部により検証・作成されています。

Layer 3 — 工学応用・現場実装

深宇宙探査機の軌道計画ソフトウェア（SPICE / GMAT）

NASAやJAXAの軌道設計エンジニアは、ランベルト問題の解法と多段階スイングバイ探索アルゴリズムを用いて、消費デルタV（推進剤）が最小となる打ち上げウィンドウと天体フライバイ順序を最適化している。

1// 2次元平面におけるスイングバイ速度計算
2interface Vector2D {
3  x: number;
4  y: number;
5}
6
7function calculateGravityAssist(
8  vProbeSun: Vector2D,   // 太陽系基準の探査機初速度 (km/s)
9  vPlanetSun: Vector2D,  // 太陽系基準の惑星公転速度 (km/s)
10  deflectionDeg: number  // 重力偏向角 delta (度)
11) {
12  // 1. 惑星基準の相対速度ベクトル
13  const uX = vProbeSun.x - vPlanetSun.x;
14  const uY = vProbeSun.y - vPlanetSun.y;
15  const uSpeed = Math.sqrt(uX * uX + uY * uY);
16
17  // 2. 相対速度ベクトルの回転
18  const rad = (deflectionDeg * Math.PI) / 180;
19  const uRotatedX = uX * Math.cos(rad) - uY * Math.sin(rad);
20  const uRotatedY = uX * Math.sin(rad) + uY * Math.cos(rad);
21
22  // 3. 太陽系基準の離脱速度ベクトル
23  const vOutX = uRotatedX + vPlanetSun.x;
24  const vOutY = uRotatedY + vPlanetSun.y;
25  const vOutSpeed = Math.sqrt(vOutX * vOutX + vOutY * vOutY);
26
27  const initialSpeed = Math.sqrt(vProbeSun.x * vProbeSun.x + vProbeSun.y * vProbeSun.y);
28
29  return {
30    initialSpeed: initialSpeed.toFixed(2),
31    finalSpeed: vOutSpeed.toFixed(2),
32    speedGain: (vOutSpeed - initialSpeed).toFixed(2),
33  };
34}
35
36// 例: 木星（公転速度 13 km/s）によるスイングバイ
37const result = calculateGravityAssist(
38  { x: 10, y: 5 },
39  { x: -13, y: 0 },
40  60
41);
42console.log(`接近前速度: ${result.initialSpeed} km/s ➔ 離脱後速度: ${result.finalSpeed} km/s (加増: +${result.speedGain} km/s)`);