はじめに
宇宙船が効率よく加速できるスイングバイという航法がある。天体を利用してスピードアップを図るという技術だが,皆さんは疑問に思ったことはないだろうか?
天体に近づいて,そして再び天体から離れるだけでどうして加速できるのだろうか?そもそも天体から見れば,宇宙船が近づいてくると重力によって加速されるが,離れるときは逆に重力によって減速される。エネルギー保存則から考えて元の速度に戻るだけではないだろうか?
つまり進入前の速度 $v$ と脱出後の速度 $v'$ の大きさは等しく,$|v| = |v'|$ になるのではないか?という疑問だ。
高校レベルの数学で考える
簡単のため一次元の運動で考えよう。質量 $m$ の宇宙船が質量 $M$ の天体と完全弾性衝突して跳ね返るものとする。衝突前の宇宙船の速さ $v$ とし,衝突後 $v'$ とする。一方,衝突前の天体の速さ $V$,衝突後 $V'$ とおく。なお座標軸の方向は左向きとする。衝突前の宇宙船は右向きに進んでいるので $v < 0$ である。
なお,いくらなんでも完全弾性衝突はおかしいと思う人は図4に示すように重力を利用して180度ターンしたものとイメージされたい。
いずれにせよ宇宙船の速度 $v$ および $v'$,天体の速度 $V$ および $V'$ は,お互い十分離れた距離における速度とする。近づくと重力による加減速が生じるが,エネルギー損失は発生しないことからその過程をブラックボックスとみなして,初期状態と最終状態だけを考えようという訳だ。
まず,エネルギー保存則より
\frac{1}{2}mv^2 + \frac{1}{2}MV^2 = \frac{1}{2}m(v')^2 + \frac{1}{2}M(V')^2 \tag{1}
となる。一方,運動量保存則より
mv + MV = mv' + MV' \tag{2}
が得られる。式 $(1)$ の両辺を2倍したものと式 $(2)$ の両辺をそれぞれ $M$ で割る。宇宙船と天体の質量比 $a = m/M$ とおくと
av^2 + V^2 = a(v')^2 + (V')^2 \tag{3}
av + V = av' + V' \tag{4}
が得られる。式 $(4)$ より,$V' = a(v - v') + V$ を式 $(3)$ に代入して $V'$ を消去する。
av^2 + V^2 = a(v')^2 + \left( a(v - v') + V \right)^2 \tag{5}
両辺から $V^2$ の項が消せる。
\require{cancel}
av^2 + \xcancel{V^2} = a(v')^2 + a^2(v - v')^2 + 2a(v - v')V + \xcancel{V^2} \tag{6}
式 $(6)$ の両辺を $a$ で割る。
v^2 = (v')^2 + a(v - v')^2 + 2(v - v')V \tag{7}
式 $(7)$ の右辺を左辺に移項する。
v^2 - (v')^2 - a(v - v')^2 - 2(v - v')V = 0 \tag{8}
うっすらと共通項 $v - v'$ が見えるので因数分解してまとめると
(v - v')(v + v') - a(v - v')^2 - 2(v - v')V = 0 \tag{9}
(v - v')\left\lbrace (v + v') - a(v - v') - 2V \right\rbrace = 0 \tag{10}
(v - v')\left\lbrace (1 - a) v + (1 + a)v' - 2V \right\rbrace = 0 \tag{11}
となる。これより
\left\lbrace\begin{align}
v - v' = 0\\
(1 - a) v + (1 + a)v' - 2V = 0
\end{align} \right. \tag{12}
が得られる。$v'$ について求めると,有意な解は($v' = v$ を除いて)
v' = \frac{2}{1 + a}V - \frac{1 - a}{1 + a}v \tag{13}
となる。宇宙船と天体の質量比 $a = m/M \ll 1$ だから
v' \approx 2V - v \tag{14}
と近似できる。仮に天体が静止していた場合 $V = 0$ を考えると $v' \approx -v$ となる。これは壁に衝突して跳ね返ったようなもので,速度の大きさを維持したまま向きが正反対になっている。
衝突後の天体の速度 $V'$ も求めてみよう。
V' = a(v - v') + V \tag{15}
式 $(15)$ に式 $(13)$ を代入する。
V' = a \left[v - \left( \frac{2}{1 + a}V - \frac{1 - a}{1 + a}v \right) \right] + V \tag{16}
V' = a \left(v - \frac{2}{1 + a}V + \frac{1 - a}{1 + a}v \right) + V \tag{17}
V' = a \left(\frac{2}{1 + a}v - \frac{2}{1 + a}V \right) + V \tag{18}
V' = \frac{1 - a}{1 + a}V + \frac{2a}{1 + a}v \tag{19}
これも $a \ll 1$ とすれば
V' \approx V \tag{20}
と近似できる。これは予想通りであり,天体レベルの質量になれば宇宙船がぶつかってもビクともしないのだ。
式 $(14)$ から式 $(20)$ を引いて天体から見た衝突後の宇宙船の相対速度 $v' - V$ の近似値を求めてみると,
v' - V = (2V - v) - V = V - v \tag{21}
となり,衝突前の相対速度 $v - V$ と大きさは同じで向きが正反対になる。念のため,厳密に求めてみよう。式 $(13)$ から 式 $(19)$ を引くと
v' - V' = \left( \frac{2}{1 + a} V - \frac{1 - a}{1 + a} v \right)
- \left( \frac{1 - a}{1 + a} V + \frac{2a}{1 + a} v\right) \tag{22}
v' - V' = \left( \frac{2}{1 + a} - \frac{1 - a}{1 + a} \right) V
- \left( \frac{1 - a}{1 + a} + \frac{2a}{1 + a} \right) v \tag{23}
v' - V' = V - v \tag{24}
となり,厳密に計算しても同じ結果になることが示された。
まとめ
衝突前後の速度を表1にまとめる。静止座標系から見た絶対速度としては,宇宙船は元の速度とは正反対の方向に跳ね返り,さらに天体の速度 $V$ の2倍だけ加速されている。 ただし,天体から見た宇宙船の相対速度としては衝突前後で大きさは変わらず,向きが正反対になっただけである。
記号 | 衝突前 | 衝突後 | |
---|---|---|---|
宇宙船の絶対速度 | ① | $v$ | $\approx 2V - v$ |
天体の絶対速度 | ② | $V$ | $\approx V$ |
天体から見た宇宙船の相対速度 | ①-② | $v - V$ | $V - v$ |
最後に
冒頭で完全弾性衝突はさすがに非現実的過ぎるので,重力ターンをイメージして欲しいと述べたが,たぶん重力だけで180度ターンするのは無理だ。というのもスイングバイの単純な軌道は放物線または双曲線になるので,来た方向に戻ることはできないからだ。
高校レベルの数学シリーズ
筆者の過去シリーズである。
- 高校レベルの数学でケプラーの法則を導く - Qiita
- 高校レベルの数学で歳差周期2万6000年を導く - Qiita
- 高校レベルの数学で考える~スペースコロニーの中でジャンプしたら無重力状態になるの? - Qiita