二重楕円遷移のソースを表示

[[File:Bi-elliptic transfer.svg|thumb|高度の低い円軌道(青色)から高度の高い円軌道(赤色)に移る二重楕円遷移]]
[[File:Hohmann transfer orbit.svg|thumb|ホーマン遷移軌道(比較用)]]
{{Astrodynamics}}
[[宇宙工学]]と[[航空宇宙工学]]で'''二重楕円遷移'''（にじゅうだえんせんい、{{Lang-en|Bi-elliptic transfer}}) とは、宇宙船を元の[[軌道]]から異なる軌道へ移す[[軌道マヌーバ]]の一つであり、特定の状況では[[ホーマン遷移軌道]]より推進剤の消費が少ない。

二重楕円遷移での宇宙船の経路は、近点の異なる二つの[[楕円軌道]]を接合したものである。最初の燃焼では宇宙船が{{仮リンク|中心天体|en|Primary body}}から任意の長さ<math>r_b</math>を遠点をもつ楕円軌道へ投入される。次の燃焼はこの楕円軌道の遠点で行われ、宇宙船は目的とする軌道に[[近点]]をもつ別の楕円軌道に投入される。新しい楕円軌道の近点で最後の燃焼が行われ、宇宙船は目的の高度の軌道に乗る<ref name="Curtis">{{Cite book | last = Curtis | first = Howard | title = Orbital Mechanics for Engineering Students | page = 264 | publisher = [[Elsevier]] | year = 2005 | isbn = 0-7506-6169-0 | url = https://books.google.com/books?id=6aO9aGNBAgIC}}</ref>。

二重楕円遷移はホーマン遷移に比べ、必要な燃焼回数が一回分増えて時間もかかるという欠点があるが、中間の楕円軌道の選び方によって遷移する軌道の[[軌道長半径]]の比率が11.94以上のとき、必要な総ΔV(速度変化量)が小さいという利点がある<ref name="Vallado">{{Cite book | last = Vallado | first = David Anthony | title = Fundamentals of Astrodynamics and Applications | page = 318 | publisher = Springer | year = 2001 | isbn = 0-7923-6903-3 | url = https://books.google.com/books?id=PJLlWzMBKjkC}}</ref>。

二重楕円遷移の概念は{{仮リンク|Ary Sternfeld|en|Ary Sternfeld}}によって1934年に初めて提案された{{citation needed|date=January 2014}}<!-- 一次ソースであるSternfeldの論文では、彼が「最初」だとは書かれていない(a primary source of Sternfeld's paper does not establish that he was the "first") --><ref name=sternfeld1934>{{Citation |last=Sternfeld |first=Ary J. |title=Sur les trajectoires permettant d'approcher d'un corps attractif central à partir d'une orbite keplérienne donnée |language=French |trans-title=On the allowed trajectories for approaching a central attractive body from a given Keplerian orbit |url=http://gallica.bnf.fr/ark:/12148/bpt6k31506/f711.image.langEN |journal=Comptes rendus de l'Académie des sciences |location=Paris |volume=198 |number=1 |date=1934-02-12 |pages=711–713}}.</ref>。

==計算==
===ΔV===
三回の燃焼でのΔVは{{仮リンク|Vis-vivaの式|en|vis-viva equation}}から直接得られる。

<math display="block">v^2 = \mu \left( \frac{2}{r} - \frac{1}{a} \right),</math>

ただし
* <math>v</math> は宇宙船の軌道速度
* <math>\mu = GM</math> は中心天体の重力係数([[万有引力定数]]と質量の積)i
* <math>r</math> は宇宙船と中心天体との距離(円軌道の場合は半径)
* <math>a</math> は宇宙船の[[軌道長半径]]

である。

以下では
* <math>r_1</math> は最初の円軌道の半径
* <math>r_2</math> は最後の円軌道の半径
* <math>r_b</math> は遷移途中の中間軌道の遠点でマヌーバで任意に決めることができる変数
* <math>a_1</math> と <math>a_2</math> は二つの楕円軌道の軌道長半径（それぞれ
<math display="block">a_1 = \frac{r_1 + r_b}{2},</math> と <math display="block">a_2 = \frac{r_2 + r_b}{2}.</math>
によって得られる。）
である。

[[File:Bi-parabolicTransfer.png|thumb|青い低軌道から赤い高軌道へ遷移する二重楕円遷移の概念図]]

半径が<math>r_1</math>の円軌道(図中青)から始める。[[順行]]噴射(図中1)で宇宙船は最初の楕円軌道(図中青緑)に移る。必要なΔVの大きさは
<math display="block">\Delta v_1 = \sqrt{ \frac{2 \mu}{r_1} - \frac{\mu}{a_1}} - \sqrt{\frac{\mu}{r_1}}. </math>
である。

宇宙船が中心天体から距離<math>r_b</math>のところにある最初の楕円軌道の遠点に達すると、二回目の順行噴射(図中2)が軌道近点を目標軌道と接するまで上昇させ、宇宙船を二つ目の楕円軌道(図中オレンジ)に投入する。二回目の噴射で要求されるΔVの大きさは

<math display="block">\Delta v_2 = \sqrt{\frac{2 \mu}{r_b} - \frac{\mu}{a_2}} - \sqrt{\frac{2 \mu}{r_b} - \frac{\mu}{a_1}}.</math>
である。

最後に、近点で逆行噴射(図中3)を行い、半径が<math>r_2</math>目標の円軌道(図中赤)に到達する。要求されるΔVの大きさは

<math display="block">\Delta v_3 = \sqrt{\frac{2 \mu}{r_2} - \frac{\mu}{a_2}} - \sqrt{\frac{\mu}{r_2}}.</math>

である。

もし、<math>r_b = r_2</math>なら、ホーマン遷移と同じである(その場合、<math>\Delta v_3</math>が0という扱いになる)。こうして、ホーマン遷移は噴射を二回しか必要としない特殊な二重楕円遷移とすることで、軌道遷移をより一般的な枠組みで扱えるようになる。

二重楕円遷移による、ホーマン遷移に対する燃料の最大節約量は、<math>r_b = \infty</math>を仮定して計算される。この場合では<math>\Delta v</math>の節約量は合計で <math display="inline">\sqrt{\mu/r_1} \left(\sqrt{2} - 1\right) \left(1 + \sqrt{r_1/r_2}\right)</math>となる。この場合、遷移軌道はもはや楕円ではなく放物線であるため、遷移軌道の経路は''二重放物線''と呼ぶべきものになり、遷移にかかる時間は無限大になる。

===二重楕円遷移に必要な時間===
ホーマン遷移のように、二重楕円遷移に使われる二つの遷移軌道は楕円軌道を半分だけ取り出したものである。つまり、それぞれの遷移軌道に滞在する時間は、完全な楕円軌道の周期の半分である。[[軌道周期]]の式と上で定義した文字を使って

<math display="block">T = 2 \pi \sqrt{\frac{a^3}{\mu}}.</math>

遷移にかかる合計時間<math>t</math>はそれぞれの半分の軌道の合計だから



<math display="block">t_1 = \pi \sqrt{\frac{a_1^3}{\mu}} \quad \text{and} \quad t_2 = \pi \sqrt{\frac{a_2^3}{\mu}},</math>
そして最後に

<math display="block">t = t_1 + t_2.</math>

== ホーマン遷移との比較 ==
=== ΔＶ ===
[[File:Bielliptic transfers comparison.svg|thumb|upright=1.5|遷移する円軌道の半径比を変えた時の、ホーマン遷移に必要なΔＶ(太い黒)と二重楕円遷移に必要なΔＶ(色付き)の変化]]

この図は、軌道半径が<math>r_1</math>の円軌道から<math>r_2</math>の円軌道に移るときの総<math>\Delta v</math>を示している。ここでの<math>\Delta v</math>は出発する軌道の軌道速度<math>v_1</math>で規格化されており、出発軌道と到達軌道の半径比の関数<math>R \equiv r_2 / r_1</math>でプロットし、一般の比較に使えるようにしてある(つまり<math>r_1</math>や<math>r_2</math>の実際の値ではなく、その比率にのみ従う)。<ref name="Vallado" />

太い黒の曲線はホーマン遷移での<math>\Delta v</math>を表し、薄い色付きの曲線はパラメーター<math>\alpha \equiv r_b / r_1</math>を変えたさまざまな二重楕円遷移での<math>\Delta v</math>を表す。 <math>r_b</math>は遷移途中の楕円軌道の軌道長半径であり、それを出発軌道の軌道半径で規格化して、曲線の横に示している。差し込み図は、二重楕円遷移の効率がホーマン遷移を初めて上回る箇所を拡大したものである。

この図からは、半径比<math>R</math>が11.94を下回るとき常に二重楕円遷移よりホーマン遷移が有利であり、一方でもし半径比<math>R</math>が15.58を上回れば常に二重楕円遷移の方が有利であることがわかる。この関係は軌道長半径によらない(最初の軌道の半径を超えている限りは)。半径比<math>R</math>が11.94から15.58の間にあるときは、遷移に用いる楕円軌道の遠点距離<math>r_b</math>によってどちらの遷移方法が有利かは変わり、それを超えると二重楕円遷移が有利となるような遠点距離<math>r_b</math>の値が存在する。、
下の図はいくつかの条件での二重楕円遷移が有利になる<math>\alpha \equiv r_b / r_1</math>を並べたものである<ref>{{cite journal
| last1      = Gobetz
| first1     = F. W.
| last2      = Doll
| first2     = J. R.
| date       = May 1969
| title      = A Survey of Impulsive Trajectories
| journal    = AIAA Journal
| publisher  = [[American Institute of Aeronautics and Astronautics]]
| volume     = 7
| issue      = 5
| pages      = 801–834
| doi        = 10.2514/3.5231 | bibcode = 1969AIAAJ...7..801D }}</ref>

{| class="wikitable"
|+ 二重楕円遷移が最小の<math>\Delta v</math>を必要とする最小の <math>\alpha \equiv r_b/r_1</math><ref>{{Cite book
 | first = Pedro R.
 | last = Escobal
 | title = Methods of Astrodynamics
 | location = New York
 | publisher = [[John Wiley & Sons]]
 | year = 1968
 | isbn = 978-0-471-24528-5
}}</ref>
! 半径比 <math>\frac{r_2}{r_1}</math>
!最小 <math>\alpha \equiv \frac{r_b}{r_1}</math>
! 備考
|-
| <11.94 || {{n/a}} || ホーマン遷移が常に有利
|-
|  11.94 || <math>\infty</math> ||二重楕円遷移
|-
|  12    || 815.81 ||
|-
|  13    || 48.90 ||
|-
|  14    || 26.10 ||
|-
|  15    || 18.19 ||
|-
|  15.58 || 15.58 ||
|-
| >15.58 || <math>> \frac{r_2}{r_1}</math> || 常に二重楕円遷移が有利
|}

=== 遷移時間 ===
二重楕円遷移での、最長で無限に達する長い遷移時間は、二重楕円遷移の主要な欠点である。二重楕円遷移の遷移時間は次の式で表される

<math display="block">t = \pi\sqrt{\frac{a_1^3}{\mu}} + \pi\sqrt{\frac{a_2^3}{\mu}},</math>

ホーマン遷移軌道は一つの楕円軌道しか遷移に用いないため、二重楕円遷移の半分以下の遷移時間である。具体的には


<math display="block">t = \pi\sqrt{\frac{a^3}{\mu}}.</math>

=== 他の軌道マヌーバとの組み合わせ ===
二重楕円遷移は円軌道間の遷移でΔＶの点で厳格にホーマン遷移に対して優位なパラメーターの範囲が小さいにもかかわらず、推進剤の節約効果はまったく小さく、他のマヌーバと組み合わせた時に真価を発揮する。

二重楕円遷移の途中の遠点で宇宙船の速度は遅くなり、小さなΔＶで大きな近点の移動を行える。このとき、同時に軌道面の変更を行うと、通常より圧倒的に少ないΔＶしか要しない。

同じように、大気のエアロブレーキによる減速は推進剤を消費しないので、遠点で減速マヌーバを行い近点を大気圏まで下げるのは、燃焼の回数が増えるが単純に軌道の高度を下げるよりΔＶの節約になることがある。

==例==
半径{{nowrap|''r''<sub>0</sub> {{=}} 6700&nbsp;km}}の地球低軌道から半径{{nowrap|''r''<sub>1</sub> {{=}} 93 800&nbsp;km}}の軌道へ遷移する場合、ホーマン遷移では{{nowrap|2825.02 + 1308.70 {{=}} 4133.72&nbsp;m/s}}のΔＶが必要であるが、半径比が{{nowrap|''r''<sub>1</sub> {{=}} 14''r''<sub>0</sub> &gt; 11.94''r''<sub>0</sub>}}の関係にあるため、二重楕円遷移の方が推進剤を節約して遷移できる。まず、宇宙船は、3061.04&nbsp;m/sだけ加速して遠点が{{nowrap|''r''<sub>2</sub> {{=}} 40''r''<sub>0</sub> {{=}} 268 000&nbsp;km}}の楕円軌道に入り、 遠点でだけ608.825&nbsp;m/sまで加速して近点が{{nowrap|''r''<sub>1</sub> {{=}} 93 800&nbsp;km}}の楕円軌道に入る。最後に近点で447.662&nbsp;m/sだけ減速し目標軌道に入る。ΔＶは合計で4117.53&nbsp;m/sであり、ホーマン遷移と比べて16.19&nbsp;m/s (0.4%)の節約になる。

ΔVの節約は、遷移時間の増加を許容して遷移軌道の遠点を伸ばすことで改善する。例えば、遠点を{{nowrap|75.8''r''<sub>0</sub> {{=}} 507 688&nbsp;km}}(月軌道の1.3倍)にするとホーマン遷移に対して1%のΔVの節約になる。現実的ではないが、遠点を{{nowrap|1757''r''<sub>0</sub> {{=}} 11 770 000&nbsp;km}}(月軌道の30倍)にすると2%のΔVの節約になる。しかし遷移には4.5年かかる(し地球の重力圏を離脱してしまう)。一方ホーマン遷移は15時間34分しかかからない。

{| class='wikitable'
|+ いろいろな遷移でのΔ''v''
! colspan=2 | 遷移の種類
!             ホーマン遷移
! colspan=4 | 二重楕円遷移
|-
! colspan=2 | 遠点 (km)
|  93 800            || 268 000            || 507 688            || 11 770 000          || ∞ <!-- The infinite limit -->
|-
! rowspan=3 | 燃焼<br/>(m/s)
! 1
| {{increase}} 2825.02 || {{increase}} 3061.04  || {{increase}} 3123.62  || {{increase}} 3191.79 || {{increase}} 3194.89
|-
! 2
| {{increase}} 1308.70 || {{increase}} 608.825  || {{increase}} 351.836  || {{increase}} 16.9336 || {{steady}} 0
|-
! 3
| {{steady}} 0         || {{decrease}} 447.662  || {{decrease}} 616.926  || {{decrease}} 842.322 || {{decrease}} 853.870
|-
! colspan=2 |合計 (m/s)
| 4133.72           || 4117.53           || 4092.38           || 4051.04           || 4048.76
|-
!  colspan=2 | ホーマン遷移に対する節約
| 100%              || 99.6%             || 99.0%             || 98.0%             || 97.94%
|}
{{ubl
| {{increase}} Δ''v'' applied ''prograde''
| {{decrease}} Δ''v'' applied ''retrograde''
}}

明らかに、二重楕円遷移は最初にΔVを多く費やす
これにより、{{仮リンク|比軌道エネルギー|en|specific orbital energy}} への寄与が大きくなり、[[オーベルト効果]]により、必要なデルタ v が正味減少する原因となる。

==脚注==
{{Reflist}}

==関連項目==
{{Portal|宇宙開発}}
* {{仮リンク|ΔV予算|en|Delta-v budget}}
* [[ホーマン遷移軌道]]
* {{仮リンク|オーベルス効果|en|Oberth effect}}

{{軌道}}
{{Astro-stub}}
{{DEFAULTSORT:にしゆうたえんせんい}}
[[Category:人工衛星の軌道]]