ニューホライズンが木星をスイングバイする時の双曲線軌道

初期条件はこれだ!

(1)ロケットの影響圏境界侵入速度(vR = 19.068 km/s (68,644.8 km/h))
(2)侵入角度(半径(R)の円軌道との接線に対する角度)(θR = 14.06 °)
(3)中心天体の基準線に対する角度(θ = 94.96748 °)

(4)木星の質量(M = 1.89813 e27kg)、ニューホライズンの質量(暫定m = 1000 kg)

これだけで、双曲線軌道を算出するのか〜
その通り。です!


求める値はこれだ!

(1)AR:ロケットの動的エネルギ=
(2)A VR:ロケットの中心天体方向の動的エネルギー
(3)AHR:ロケットの影響圏境界の円軌道方向の動的エネルギー
(4)a:基準軌道半径(半交軸と同値)
(5)R:影響圏境界までの距離
(6)f:振幅の最小値(焦点距離と同値)
(7)r:近点距離
(8)vr:近点速度

では、一つづつ行きます。
(1)AR:ロケットの動的エネルギー
   AR = m x vR2 = 1000 x 68,644.82 = 4.71211 e12 je=kg・(km/h)2

(2)A VR:ロケットの中心天体方向の動的エネルギー
   A VR = AR x cos2θR = 4.71211 e12 x cos214.06° = 4.43401E+12 je

(3)AHR:ロケットの影響圏境界の円軌道方向の動的エネルギー
   A VR = AR x sin2θR = 4.71211 e12 x sin214.06° = 2.78103E+11 je

(4)a:基準軌道半径(半交軸と同値)
   基準軌道の静的エネルギー(Sa) = AVR = 4.43401 e12 je
  
   ・面積エネルギー(EA)
    EA = Em x ac = 1.64144 e16 je・km
      Em:質量エネルギー Em =mc2 = 1000 x 1.07925 e9 = 1.16479 e21 je
      ac:光速時久順軌道半径 ac = U(M + m)= 7.42426 e-31 km

    a = EA / Sa = 370.193.8 km

(5)R:影響圏境界までの距離
   R = 2 x EA / SR = 11,804,557.74 km
   影響圏境界の静的エネルギー(SR
   
   SR = EA / R = 1.39051 e11 je = 1.39051 e11

(6)f:振幅の最小値(焦点距離と同値)
   Rx = R x cosθ=7.06145 e11
   f = (-Rx + sqrt(Rx2 + 4a(a+R)))/2
   = (-7.06145 + sqrt((-7.06145)2 + 4 x 370,193.8 x (370,193.8 + 11,804,557.74)))
   = 2,694,705.8 km

(7)r:近点距離
   r = f – a = 2,694,705.8 – 370.193.8 = 2,324,512 km

(8)vr:近点速度
   近点(r)の静的エネルギー(Sr)は、
   
   S r = EA / r = 1.64144 e16 / 2324512 = 7.06145 e11 
   
   ロケットの動的エネルギーは。
   
   Ar = 2Sr + Sa =2 x 7.06145 e11 + 4.43401 e12 =5.84629 e12 je
   
   ロケットの近点速度(vr) は、  

   vr = sqrt(Ar / m) =sqrt(5.8429 e12 / 1000)
    
     = 76,461.1 km/h(=21.24 km / s)

   となります。

ちなみに、「惑星探査機の軌道計算入門」での計算結果は、
基準軌道半径(半交軸)a = 3.7019 e5 km(当方:3.7019e5 km)
近点軌道速度 vr = 21.244 km/s(当方:21.24 km/s)
と、当ラボの結果とほぼ同じ値でした。

しかし、影響圏境界までの距離(Rは)
R=4.820e7 km(当方 1.1805 e7 km)となり、
約4倍の差が出てしまった。

どちらが正しいか、検証しなければわからないが、
概ね、この考え方には、間違いがないのではないかという
確証を得ることができた。


スイングバイ少しわかった

宇宙の構造を考えているうちに、説明の最後の方だと思うが、スイングバイがうまく説明できず、ここ3年くらい考えていたのだった。しかし、ここにきてなんとなく見通しがよくなった。

参考にしていた書物は、半楊稔雄著「惑星探査機の軌道計算入門 宇宙飛翔力学への誘い」です。この中で、惑星探査機ニューホライズンの軌道計算についての記述がある。興味のある方は、購入してみてください。

その中で、自分の理論と相違点というか、わからなかったのは、スイングバイの軌道が双曲線軌道になるというところだ。

双曲線軌道にならないというのが、自論だった。その理由は、動的エネルギーがある位置の2倍の静的エネルギーになると、基準軌道が無限大となってしまうので、親が移動するという理論だからだ。それなのに、スイングバイは双曲線軌道なのだ。

ついに自分の理論崩壊か?

そこで、再度見直してみた。

エネルギーバランス理論(調和理論)からスイングバイ

スイングバイは大きく3つのフェーズに分けられる。

(1)影響圏内に入る前

(2)影響圏の中

(3)影響圏から出た後

今回は(2)について語ってみることにします。実は、(1)(3)は、まだ十分に考えていないので。

とりあえず、双曲線軌道のエネルギー構造を調べてみた。

双曲線軌道のエネルギー構造

where

  a:半交軸

  r:最近点

  R:影響圏

Sr:位置rの静的エネルギー

SR:位置Rの静的エネルギー

AR:宇宙船の動的エネルギー

Sa:半交軸の頂点の静的エネルギー

になります。なんだこれって思う人も多数いると思いますが、中心天体と距離によって、静的エネルギがどのように変化するかを表しています。用語は、半揚先生の本の言葉を使っていますが、実際の自分の理論の意味とは異なります。

このグラフも1年くらい考えて完成しました。難しそうに見えませんが、実はーSαに辿り着くまでに、試行錯誤しました。

このグラフの見方ですが、宇宙船は右から左へ動いて、中心天体に近づいて来ます。例を使った方がわかりやすいので、中心天体を木星、その親を太陽とします。地球から打ち上げた宇宙船が、木星でスイングバイするという設定です。本当は地球のパーキング軌道から木星に向かうところの軌道も考察しなくてはいけないのですが、今回は省いて、木星の近くにきてからの軌道について考えてみます。

そのまま左に宇宙船が進むとグラフのRの位置に来ます。Rは影響圏です。半揚先生の著書では、影響圏の距離は、木星から4.820e7kmになっています。しかし、自分の理論では、その10分の1くらいになります。Rの位置の静的エネルギー(SR)の2倍に動的エネルギー(AR)になります。その位置で、宇宙船の親は木星に代わります。

という感じなるはずでしたが、宇宙船は、基準軌道Rの円軌道に対して、角度を持って侵入します。つまり宇宙船の動的エネルギー(AR)は、宇宙船全体の動的エネルギーの基準軌道方向のエネルギーになります。ということは、そのほかの動的エネルギーは、どうなっちゃうのってことです。ここで気づいたのが、木星方向の動的エネルギー(Aa)です。これが、半交軸を決める要因であるはずではないかと。

スウィングバイの軌道のエネルギー関係

木星の影響圏に入った宇宙船は、グラフの2Sの2倍静的エネルギー曲線(赤色グラフ)上を移動します。

そして、rの位置まで近づいて、反対側から出て行きます。

ここで、最近点rは、影響圏に入った時の距離の半分になります。ここが重要です。

計算してみます

宇宙エネルギー定数 U = 7.42426E-31 km/kg

木星質量 m = 1.89813e+27 Kg

宇宙船質量 mx = 1.0e3 kg (質量エネルギーEmx = 1.16479e+21)

光速時基準軌道 ac = 0.00141 km

光速 c = 1.07925e+9 km/h

宇宙船の動的エネルギー AX = 4.45007E+12 kg(km/h)2=je (宇宙船、木星影響圏突入速度 VX = 66,708.828 Km/h)

宇宙船突入角度 θ=80.9553度

とすると

基準軌道方向の動的エネルギーは、AR = AX × cos( θ ) よりAR = 6.99573e+11jeになります

AR=2SRだから、影響圏の半径Rは

より、R = 4,692,696.69445 km

になります。有効桁数がめちゃくちゃですいません、

半揚先生の影響圏 R’=4.820e7 km と比較すると、10分の1くらいの大きさになります。半揚先生の影響圏をどうやって算出したのかは、本に載っていなかったので、わからないのですが、とりあえず、自分で算出した、影響圏を使うことにします。

宇宙船の木星最近点距離は、Rの2分の1なので、r = 2,346,348.34722 kmになります。実際は、ニューホライズンは、最近距離は2.25302e6 kmでした。入射角と入射速度は推定なので、若干違うかもしれませんので、ほぼ同値といっても良いのではないかと思います。

木星の影響圏での宇宙船の動的エネルギーは、木星からの距離x kmでは、

Ax = 2Sx – ( -Sa )

になります。

この式を使用すれば、木星上の宇宙船の速度がもとまります。

ということで、スイングバイ影響圏の軌道の説明でした。

では。