電気推進入門 — イオン推進機とホール推進機の動作原理と性能比較

2010年6月13日、小惑星探査機「はやぶさ」は小惑星イトカワへの往復飛行を終え、サンプルを収めたカプセルをオーストラリアの砂漠に投下しました。地球とイトカワは平均約3億km離れていて、両者の軌道は太陽を巡る別々の楕円。この往復に必要な総 $\Delta v$ は4km/sを大きく超えます。もしこれを化学推進だけで行おうとしたら、出発質量の9割以上が推進剤になってしまい、観測機器を載せる余裕はほとんど残りません。

それでも、はやぶさは打ち上げ時の乾燥質量わずか380kgあまりで往復を成し遂げました。その心臓部にあったのがマイクロ波放電型イオンエンジン μ10です。化学推進では到底届かない比推力 $I_{sp} \approx 3000$ 秒を、4基のエンジンで継続的に絞り出すことで、薄く長く推進し続けて巨大な $\Delta v$ を稼ぐ——これが電気推進の世界です。

電気推進は今や深宇宙探査の専売特許ではありません。SpaceXのStarlink衛星はクリプトンを推進剤とするホール推進機で軌道調整と最終的な軌道離脱を行っており、コンステレーション運用全体の経済性を支えています。欧州・日本のBepiColomboは水星到達のために4基のイオンエンジンで太陽方向への巨大な減速 $\Delta v$ を稼ぎ、NASAのNEXT-Cやミシガン大のX3クラスのホールスラスタは将来の有人火星探査のキー技術として開発が進んでいます。

本記事の内容

• 電気推進が「電気エネルギーで排気速度 $v_e$ を上げる」とはどういうことか、化学推進との根本的違い
• グリッド型イオン推進機の動作原理と、Child-Langmuir則による空間電荷制限電流密度の導出
• ホール推進機の動作原理(E×Bドリフト)と、なぜ空間電荷制限を超えられるのか
• 化学推進 vs 電気推進の推力・$I_{sp}$・効率・寿命・複雑性の定量比較
• Pythonで Child-Langmuir則、E×Bドリフト、はやぶさの $\Delta v$ 予算シミュレーション

前提知識

この記事を読む前に、以下の記事を読んでおくと理解が深まります。

• ツィオルコフスキーの公式と多段ロケット
• 比推力と推進剤の物理
• 化学推進の燃焼熱力学
• 軌道マヌーバと $\Delta v$
• ガウスの法則(電場)

なぜ電気で推進できるのか — 直感

「電気で物体を加速して投げ出す」と聞くと、なんとなく非効率に思えるかもしれません。電気はあくまでエネルギーであり、運動量を直接持っているわけではないからです。しかしロケット推進の本質は「質量を後ろに投げ出した反作用で前に進む」ことであり、推進剤さえあれば、それを加速するエネルギー源は何でも構いません。化学推進では燃焼の化学エネルギーが推進剤自身を加熱・膨張させて加速する一方、電気推進では電源(太陽電池や原子炉)から得た電気エネルギーで、電場や磁場を介して推進剤を加速します。

エネルギー源を化学反応から切り離した瞬間、何が変わるでしょうか。化学推進では、燃焼室温度は推進剤分子の解離限界(おおむね3500K前後)で頭打ちになり、排気速度 $v_e$ は理論上 4.5km/s 程度が限界です。これは比推力 $I_{sp} = v_e/g_0$ で言うと約460秒。化学エネルギーが推進剤分子の中に閉じ込められているため、それを取り出した瞬間に温度の上限が決まってしまうのです。一方、電気推進では電源容量を増やせば加速電圧をいくらでも上げられ、$v_e$ を化学推進の10倍以上に押し上げることができます。実用化されているイオンエンジンの $v_e$ は 30km/s 超、$I_{sp}$ で 3000〜10000秒に達します。

ツィオルコフスキーの公式 $\Delta v = v_e \ln(m_0 / m_f)$ を思い出すと、$v_e$ が10倍になれば同じ $\Delta v$ を稼ぐのに必要な推進剤質量比が指数的に減ることがわかります。たとえば $\Delta v = 10$ km/s のミッションを考えると、化学推進($v_e = 4.5$ km/s)では $m_0/m_f \approx 9$ で乾燥質量の8倍の推進剤が必要なのに対し、イオンエンジン($v_e = 30$ km/s)では $m_0/m_f \approx 1.4$ で済みます。後者では推進剤質量は乾燥質量の4割程度。ペイロード比が劇的に改善するのが電気推進の最大の魅力です。

ただし「うまい話」には代償があります。電気推進の推力は化学推進に比べて4〜6桁小さいのが普通です。これは電源容量と排気速度の幾何学的トレードオフから直接出てくる事実で、後の節で定量的に見ます。低推力では地球から打ち上げることはできず、化学ロケットで宇宙に運んでもらう必要があるうえ、目的軌道への遷移には数か月から数年の連続加速が必要になります。電気推進ミッションは「インパルス的マヌーバを足し合わせる」古典的軌道力学とは別の、「連続加速の螺旋」として設計されます。

ここまでで電気推進の動機と全体像を押さえました。では具体的にどうやってプラズマを加速するのか。代表的な2方式であるグリッド型イオン推進機とホール推進機の動作原理を順に見ていきます。

イオン推進機(グリッド型)の動作原理

構成と動作の流れ

グリッド型イオン推進機の構造は、放電室・引き出しグリッド系・中和器の3つに大別されます。放電室の中で中性ガス(主にキセノン Xe)をイオン化し、放電室と外部空間の電位差を利用してイオンを高速に引き出し、最後に外部で電子と再結合させて中性ビームとして放出する——これが基本の流れです。

放電室でのイオン化方法には2系統あります。ひとつはカウフマン型で、熱陰極から放出された電子を磁場と電場で閉じ込め、Xe原子に衝突させて電子をたたき出します。もうひとつがマイクロ波放電型(ECR型)で、はやぶさのμ10エンジンが採用しました。共鳴周波数(2.45GHz帯)のマイクロ波を磁場と組み合わせて電子サイクロトロン共鳴を起こし、電子を加熱してイオン化する方式で、寿命が長く電力効率も高いのが特長です。

イオン化されたXe$^+$は放電室内で約1000Vの陽電位に置かれ、グリッド系で外部空間(接地電位)との電位差により加速されます。グリッドは典型的に3枚構成で、放電室側からスクリーングリッド(高電位)・アクセラレータグリッド(負電位)・デセラレータグリッド(接地電位)の順に並びます。スクリーンとアクセラレータの間で大半の加速が行われ、デセラレータは外部電子のグリッド系内への流入(電子バックストリーミング)を防ぐ役割を担います。

引き出されたイオンビームはそのままでは正電荷の集まりであり、宇宙機本体は反作用で負に帯電して、最終的にビーム自体がブレーキをかけられてしまいます。これを防ぐのが中和器(neutralizer) で、ビーム下流に置かれた電子源からビームに合流する電子を供給して、トータルとしての電荷流出量をゼロに保ちます。中和器が止まればイオンエンジンも止まる——この事実は、電気推進が「準中性プラズマ」を扱う技術であることを象徴的に示しています。

排気速度と加速電圧

シングリッド系での加速はエネルギー保存則から簡単に書けます。質量 $m_i$、電荷 $q$ のイオンが電位差 $V$ で加速されると、出口での運動エネルギーは

$$ \frac{1}{2} m_i v_e^2 = q V $$

したがって排気速度は

$$ v_e = \sqrt{\frac{2 q V}{m_i}} $$

となります。Xe(原子量131.3, $m_i = 2.18 \times 10^{-25}$ kg)を $V = 1000$ V で加速すると $v_e \approx 38$ km/s、$I_{sp} \approx 3900$ 秒となり、化学推進の8倍以上の比推力が得られます。$V$ を上げれば $v_e$ も上がりますが、後述の通り推力密度との兼ね合いと、グリッドの絶縁・スパッタリング寿命の制約から、実用的には $V = 500\sim 3000$ V程度に収まります。

ここで重要なのは、$v_e$ がイオン質量 $m_i$ の平方根に反比例することです。なぜキセノンが選ばれるか——これは「重い原子を選ぶことで、同じ加速電圧でも単位イオン1個あたりの運動量を稼ぎ、しかも電離エネルギーが比較的低い(12.13 eV)ためイオン化効率がよい」という最適化の結果です。アルゴン(40 amu)やクリプトン(83.8 amu)、最近ではヨウ素(127 amu)も研究されており、Starlinkはコスト面からクリプトンを選択しました。

推力の表式

推力 $F$ は単位時間あたりに射出する運動量です。質量流量を $\dot{m}_i$ とすると

$$ F = \dot{m}_i v_e $$

電流 $I_b$ でイオンビームを引き出している場合、$\dot{m}_i = (m_i/q) I_b$ なので

$$ F = \frac{m_i}{q} I_b \cdot v_e = \frac{m_i}{q} I_b \sqrt{\frac{2qV}{m_i}} = I_b \sqrt{\frac{2 m_i V}{q}} $$

これを電力 $P_b = I_b V$ で表せば

$$ F = \frac{P_b}{V}\sqrt{\frac{2 m_i V}{q}} = P_b \sqrt{\frac{2 m_i}{q V}} $$

となります。同じ電力 $P_b$ を使う場合、加速電圧 $V$ を上げると推力は逆数の平方根で減るというトレードオフが明確に見えます。「電力一定なら高比推力ほど低推力」というのは電気推進の宿命です。

ここまでで「電位差で電荷を加速する」というイオン推進機の基本は理解できました。しかし「もっと電流 $I_b$ を流せば推力は上がるのでは」と思いますよね。実は、グリッドの間に流せる電流密度には物理的な上限があり、これがイオン推進機の推力密度を厳しく制約します。次節でその上限——Child-Langmuir則——を導きます。

Child-Langmuir則と空間電荷制限

直感: 自分のうしろは見えない

イオンを連続的に引き出すと、グリッド間の空間にはイオンが充満して空間電荷(space charge) を形成します。この正電荷は自分自身がつくる電場で後続のイオンの加速を妨げ、結果として「これ以上電流を増やせない」上限を生みます。蛇口を全開にすると水が排水口に追いつかず詰まるようなものです。

形式的には、ポアソン方程式

$$ \frac{d^2 \phi}{dx^2} = -\frac{\rho(x)}{\epsilon_0} $$

において、空間電荷密度 $\rho(x)$ がイオン流自身から決まり、電位 $\phi(x)$ が加速電位を歪めるという自己無撞着問題を解くことで上限電流密度を求めます。境界条件は「陰極側(x=0)で $\phi = 0$、陽極側($x = d$)で $\phi = -V$」、加えて「陰極側で電場ゼロ($d\phi/dx|_0 = 0$、これが電流が限界に達した条件)」とします。

導出

定常流では電流密度 $j$ は一定で、速度 $v(x)$ とイオン密度 $n(x)$ で $j = q n v$ と書けます。エネルギー保存から $\frac{1}{2} m_i v^2 = -q \phi$ なので($\phi$ は陽極側で負と定義)、$v = \sqrt{-2q\phi/m_i}$。これを用いてイオン密度を電位で表すと

$$ n(x) = \frac{j}{q v(x)} = \frac{j}{q}\sqrt{\frac{m_i}{-2 q \phi(x)}} $$

ポアソン方程式に空間電荷 $\rho = q n$ を代入すると

$$ \frac{d^2\phi}{dx^2} = -\frac{q n}{\epsilon_0} = -\frac{j}{\epsilon_0}\sqrt{\frac{m_i}{-2 q \phi}} $$

ここで $u = -\phi$($u \geq 0$)とおいて符号を整理すると、

$$ \frac{d^2 u}{dx^2} = \frac{j}{\epsilon_0}\sqrt{\frac{m_i}{2 q u}} $$

両辺に $du/dx$ をかけて $x$ で積分すると、左辺は $\frac{1}{2}(du/dx)^2$ となり、右辺は

$$ \int \frac{j}{\epsilon_0}\sqrt{\frac{m_i}{2 q u}} \, du = \frac{j}{\epsilon_0}\sqrt{\frac{m_i}{2 q}} \cdot 2 \sqrt{u} $$

両辺を整理して

$$ \frac{1}{2}\left(\frac{du}{dx}\right)^2 = \frac{2 j}{\epsilon_0}\sqrt{\frac{m_i u}{2 q}} + C $$

陰極での電場ゼロ条件 $du/dx|_0 = 0$ かつ $u(0) = 0$ から積分定数 $C = 0$。両辺の平方根を取って $du/dx > 0$ の枝を選ぶと

$$ \frac{du}{dx} = 2\sqrt{\frac{j}{\epsilon_0}}\left(\frac{m_i}{2 q}\right)^{1/4} u^{1/4} $$

これを変数分離して $u^{-1/4} du = \cdots dx$ の形にし、$x$ で積分すると

$$ \frac{4}{3} u^{3/4} = 2 \sqrt{\frac{j}{\epsilon_0}}\left(\frac{m_i}{2q}\right)^{1/4} x $$

陽極側 $x = d$ で $u = V$ を代入し、$j$ について解くと

$$ \boxed{\;j_{\rm CL} = \frac{4 \epsilon_0}{9} \sqrt{\frac{2 q}{m_i}} \cdot \frac{V^{3/2}}{d^2}\;} $$

これがChild-Langmuir則(空間電荷制限電流則) です。電流密度は加速電圧の $3/2$ 乗に比例し、ギャップ距離の2乗に反比例するという有名な式が得られました。

推力密度の上限

電流密度 $j_{\rm CL}$ に対応する推力密度(単位面積あたりの推力)$F/A$ は、$F = \dot m v_e = (m_i/q) j A \cdot v_e$ から

$$ \frac{F}{A} = \frac{m_i}{q} j_{\rm CL} v_e = \frac{m_i}{q} \cdot \frac{4\epsilon_0}{9}\sqrt{\frac{2q}{m_i}}\frac{V^{3/2}}{d^2} \cdot \sqrt{\frac{2qV}{m_i}} $$

ここで $\sqrt{2q/m_i}\cdot\sqrt{2qV/m_i} = (2q/m_i)\sqrt{V}$ なので、

$$ \frac{F}{A} = \frac{8\epsilon_0}{9} \frac{V^2}{d^2} $$

推力密度は加速電圧の2乗に比例し、ギャップ距離の2乗に反比例します。実用的な値を入れてみると、$V = 1500$ V、$d = 2$ mm で $F/A \approx 4$ N/m$^2$。これがどれほど小さいか感覚をつかむには、化学推進(液体ロケットの燃焼室では1MN/m$^2$オーダー)と比べると一目瞭然で、5桁以上の差があります。

推力密度を上げるには $V$ を上げるか $d$ を縮めるしかなく、$V$ は絶縁とイオン速度(高すぎるとグリッドへのスパッタダメージが急増)で頭打ち、$d$ はグリッド製作精度と熱変形で頭打ち。結果として、グリッド型イオンエンジンの推力は1基あたり数十〜数百mNが限界です。はやぶさのμ10は1基あたり 8 mN(! 1円玉約1枚の重さ)で、4基並列にして 32 mN を発生していました。

この空間電荷制限を打ち破る発想として、「正電荷のイオンだけでなく、電子も一緒に混ぜて準中性プラズマ全体を加速すれば、空間電荷効果はキャンセルされるのでは」というアイデアが出てきます。これを実現するのがホール推進機です。

ホール推進機の動作原理 (E×Bドリフト)

設計思想

ホール推進機(Hall thruster)は、グリッドを使わずに準中性プラズマ全体に電場をかけて加速する装置です。正電荷の空間電荷をキャンセルするためには、加速領域に電子を共存させる必要があります。しかし電子は質量がイオンの10万倍以上小さいので、同じ電場をかければ電子の方が圧倒的に速く加速されてしまい、加速電場をショートカットしてしまいます。

ここで使う巧妙な仕掛けが、半径方向の磁場による電子の閉じ込めです。電子は軸方向の電場 $E$ と半径方向の磁場 $B$ を同時に感じると、ローレンツ力 $\bm F = -e(\bm E + \bm v \times \bm B)$ により方位方向の $\bm E \times \bm B$ ドリフト運動に拘束されます。ドリフト速度は

$$ \bm v_{E\times B} = \frac{\bm E \times \bm B}{B^2} $$

このドリフトは方位方向(円周方向)の閉じた経路を描くため、電子は陰極から陽極まで「直線的に」走ることができず、軸方向の電流に対しては高い実効抵抗として働きます。一方、イオンはラーモア半径が大きすぎて(イオンは磁化されない)、電場をそのまま受けて軸方向に加速されます。

つまりホール推進機は「電子は閉じ込め、イオンは出す」という選択的な振る舞いを磁場で実現する装置です。加速領域には電子とイオンがほぼ同密度で共存し、準中性が保たれるため、空間電荷制限を回避できます。

ホール電流

方位方向ドリフトする電子の総量は、その経路にそってホール電流 $I_H$ として流れます。「ホール推進機」の名前はこの方位電流に由来します。半径方向磁場 $B_r$ と軸方向電場 $E_z$ に対するホール電流の電流密度は

$$ j_\theta = -e n_e v_{E\times B} = -e n_e \frac{E_z}{B_r} $$

電子の磁化を表す指標がホールパラメータ $\Omega_e = \omega_{ce}\tau_{ei}$ で、サイクロトロン角周波数 $\omega_{ce} = eB/m_e$ と電子-イオン衝突周期 $\tau_{ei}$ の積です。$\Omega_e \gg 1$ で電子は十分に磁化されており、$\Omega_e \sim 100$ オーダーが典型的なホール推進機の設計条件です。

SPTとTAL

ホール推進機の2大方式は、SPT(Stationary Plasma Thruster, 定常プラズマ推進機) とTAL(Thruster with Anode Layer, 陽極層型推進機) です。SPTはロシアで発展した方式で、加速チャネルが長くセラミック壁(BN-SiO$_2$)で囲まれており、壁との電子衝突がプラズマ生成を助けます。Starlinkや多くの商用衛星で使われているのはこの系統です。一方TALは加速領域が陽極のごく近傍に集中し、金属壁を用いる方式で、より高い比推力に向きます。

排気速度と効率

イオンは陽極から陰極(中和器)までの電位差 $V_d$ で加速されるため、グリッド型と同じく

$$ v_e = \sqrt{\frac{2 e V_d}{m_i}} $$

ただし実機では電圧降下が完全には加速エネルギーに転換されず、効率 $\eta_v \approx 0.8\sim 0.9$ が掛かります。典型的なホール推進機の動作電圧は $V_d = 200\sim 400$ V、$I_{sp} = 1500\sim 2500$ 秒で、グリッド型イオンエンジンより低めです。これは「電子閉じ込めに磁場を使う」アーキテクチャ上、現実的な磁束密度で安定動作させるには加速電圧があまり高くできないためです。

中和器とプラズマ準中性

ホール推進機にも中和器が必要です。中和器カソードから供給される電子の一部はイオンビームに合流して下流の準中性を維持し、もう一部は方位ドリフトしながら陽極へ流れ込んで放電電流を担います。「電子はホール電流として方位方向にぐるぐる回り、ゆっくり軸方向にも染み込んで陽極に達する」という二重の役割を果たすのが、ホール推進機の動作の妙味です。

ここまでで2大方式の動作原理がそろいました。実機の性能はそれぞれどう違い、化学推進と比べてどんな位置付けになるのか、定量比較していきます。

化学推進との性能比較

性能数値の俯瞰

横軸を $I_{sp}$、縦軸を推力で対数プロットすると、化学推進と電気推進は明確に分離した2つのクラスタを形成します。化学推進は「右上(高推力・中$I_{sp}$)」、電気推進は「左下(低推力・高$I_{sp}$)」。間を埋めるアークジェットや抵抗加熱は「電熱推進」として中間に位置します。

電力と質量のトレードオフ

電気推進の本質的制約として、ジェット出力 $P_{\rm jet} = \frac{1}{2}\dot m v_e^2$ は電源出力に比例することがあります。$\dot m = F/v_e$ を代入すると

$$ P_{\rm jet} = \frac{1}{2} F v_e $$

すなわち推力 $F$ と排気速度 $v_e$ の積の半分がジェット出力。電源効率を $\eta$ とすれば必要電源容量は $P_e = P_{\rm jet}/\eta$ となり、$P_e$ は宇宙機質量(太陽電池や核電源の質量)を制約します。電源比質量 $\alpha = m_{\rm power}/P_e$ [kg/W] を導入すると、電気推進ミッションは

$$ m_{\rm tot} = m_{\rm dry} + m_p + \alpha P_e $$

という重量配分で、$v_e$ を上げると $m_p$ は減る一方で $P_e \propto v_e$ で増えるため $m_{\rm power}$ が増えます。両者の和を最小化する最適排気速度が、ミッション $\Delta v$ と $\alpha$ から決まります。これは「化学推進ではあり得ない」電気推進固有の最適化問題で、$\alpha = 0.01$ kg/W、$\Delta v = 10$ km/s なら最適 $I_{sp} \approx 3000\sim 5000$ 秒となり、ちょうどグリッド型イオンエンジンの動作領域に対応します。

ミッション設計の違い: インパルス vs 連続加速

化学推進ミッションは、エンジン点火時間(分〜十分オーダー)に対して軌道周期(時間〜日)が長いため、推力をデルタ関数(インパルス)で近似するインパルシブマヌーバで設計します。ホーマン遷移やリンドフラー遷移のような閉解が使えるのはこの近似のおかげです。

一方、電気推進では推力が小さく連続的に何か月も加速し続けるため、軌道は太陽中心の螺旋(spiral) を描きます。地球周回でも、低軌道から静止軌道への電気推進遷移は数か月かかり、軌道は何百周もの螺旋になります。設計には軌道平均化(orbital averaging)による低推力軌道最適化や、間接法・直接コロケーション法による最適制御が用いられ、化学推進のインパルス近似とは別世界の数学が必要です。

寿命とスケーリング

電気推進固有の寿命要因として、スパッタリングがあります。グリッド型ではイオンがアクセラレータグリッドに当たって金属を削り、寿命の主因(典型10000時間以上、はやぶさμ10で18000時間を達成)になります。ホール推進機ではチャネル壁(SPT)や陽極(TAL)のスパッタリングが寿命を決め、現代の磁気シールディング技術(壁面と磁力線をほぼ平行にしてイオン衝突を抑える)により50000時間級の寿命が見えてきました。

ここまでで「いつ・どんなとき」に電気推進が選ばれるのか、設計判断の根拠が固まりました。実際にこれらの方程式をPythonで動かして、はやぶさのような深宇宙ミッションを定量的に再現してみましょう。

Python実装 — 性能計算とミッション設計

Child-Langmuir則と最大電流密度

まず空間電荷制限電流密度を計算します。実用上のグリッド型イオンエンジン(キセノン、$V = 500\sim 3000$ V、$d = 0.5\sim 3$ mm)について、$j_{\rm CL}$ と推力密度がどう変わるかを見ます。

import numpy as np
import matplotlib.pyplot as plt

# 物理定数
eps0 = 8.854e-12       # 真空誘電率 [F/m]
qe   = 1.602e-19       # 素電荷 [C]
amu  = 1.661e-27       # 原子質量単位 [kg]

# キセノンイオン
m_xe = 131.3 * amu     # Xe+ の質量 [kg]

def j_child_langmuir(V, d, m_i=m_xe, q=qe):
    """Child-Langmuir則による空間電荷制限電流密度 [A/m^2]"""
    return (4 * eps0 / 9) * np.sqrt(2 * q / m_i) * V**1.5 / d**2

def thrust_density(V, d):
    """空間電荷制限下での推力密度 [N/m^2]"""
    return (8 * eps0 / 9) * V**2 / d**2

# 加速電圧と推力密度のスキャン
Vs = np.linspace(300, 3000, 100)
ds = [0.5e-3, 1.0e-3, 2.0e-3]                # 各ギャップ [m]

plt.figure(figsize=(10, 5))
for d in ds:
    plt.plot(Vs, thrust_density(Vs, d), label=f'd={d*1000:.1f} mm')
plt.xlabel('Accelerator voltage V [V]')
plt.ylabel(r'Thrust density $F/A$ [N/m$^2$]')
plt.title('Space-charge limited thrust density (Xe ion thruster)')
plt.legend()
plt.grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig('cl_thrust_density.png', dpi=150, bbox_inches='tight')
plt.show()

# はやぶさμ10相当の数値を点検
V_mu10 = 1500       # [V]
d_mu10 = 1.0e-3     # [m]
A_mu10 = np.pi * (0.05)**2 * 0.6  # グリッド開口面積(直径10cm × 開口率0.6)
F_mu10 = thrust_density(V_mu10, d_mu10) * A_mu10
print(f"V={V_mu10} V, d={d_mu10*1000} mm:")
print(f"  j_CL = {j_child_langmuir(V_mu10, d_mu10):.1f} A/m^2")
print(f"  F/A  = {thrust_density(V_mu10, d_mu10):.2f} N/m^2")
print(f"  推定推力(μ10想定) = {F_mu10*1000:.1f} mN")

このグラフは、推力密度が加速電圧の2乗とギャップ距離の逆2乗で急激にスケールすることを示します。$d = 0.5$ mm まで詰めれば $V = 3000$ V で約 90 N/m$^2$ に達しますが、実際にはこの精度でグリッドの平行度を保ち、熱変形・スパッタ寿命に耐えさせるのが工学的難題です。$d = 1$ mm、$V = 1500$ V のμ10相当条件で計算した推力 8 mN 前後は、はやぶさの実機(8 mN/基)とよく一致しています。電気推進が「グリッド面積を稼げないと推力が出ない」のがこの式から定量的に分かります。

ホール推進機の E×B ドリフト速度

電子のホールドリフトを可視化します。径方向磁場 $B_r$ と軸方向電場 $E_z$ のもとで、電子が方位方向にどれくらいの速度で回るか。

import numpy as np
import matplotlib.pyplot as plt

# 電子
me = 9.109e-31         # 電子質量 [kg]

def vExB(E, B):
    """E×Bドリフト速度 [m/s] (大きさのみ)"""
    return E / B

def cyclotron_freq(B):
    """電子サイクロトロン角周波数 [rad/s]"""
    return qe * B / me

# ホール推進機の典型条件
V_d = 300              # 放電電圧 [V]
L_acc = 0.025          # 加速領域長 [m]
E_z = V_d / L_acc      # 軸方向電場 [V/m]

Bs = np.linspace(0.005, 0.05, 100)   # 5 ... 50 mT
v_drift = vExB(E_z, Bs)
omega_ce = cyclotron_freq(Bs)
# 電子-中性衝突周波数の典型値(プラズマ密度 1e18 /m^3, Xe密度 1e19 /m^3)
nu_en = 1e7                          # [Hz]
Omega_e = omega_ce / (2*np.pi*nu_en) # ホールパラメータ

fig, ax1 = plt.subplots(figsize=(10, 5))
color = 'tab:blue'
ax1.plot(Bs*1000, v_drift/1000, color=color, lw=2, label='vExB drift [km/s]')
ax1.set_xlabel('Radial magnetic field $B_r$ [mT]')
ax1.set_ylabel('E×B drift velocity [km/s]', color=color)
ax1.tick_params(axis='y', labelcolor=color)
ax1.grid(True, alpha=0.3)

ax2 = ax1.twinx()
color = 'tab:red'
ax2.plot(Bs*1000, Omega_e, color=color, lw=2, ls='--', label='Hall parameter')
ax2.set_ylabel('Hall parameter $\\Omega_e$', color=color)
ax2.tick_params(axis='y', labelcolor=color)
ax2.axhline(1, color='red', alpha=0.3, ls=':')

plt.title(f'Hall thruster: E×B drift & Hall parameter ($V_d$={V_d}V, $L$={L_acc*1000}mm)')
fig.tight_layout()
plt.savefig('hall_drift.png', dpi=150, bbox_inches='tight')
plt.show()

# SPT-100相当の動作点
B_typ = 0.020         # 20 mT
print(f"B={B_typ*1000:.0f} mT, E={E_z:.0f} V/m:")
print(f"  E×B drift = {vExB(E_z, B_typ)/1000:.1f} km/s")
print(f"  Hall parameter Omega_e = {cyclotron_freq(B_typ)/(2*np.pi*nu_en):.1f}")

このプロットから、典型的なホール推進機の動作領域(20 mT付近)で電子のドリフト速度は約 600 km/s に達することがわかります。これは加速されるイオン速度(20 km/s 程度)の30倍以上であり、電子は加速方向(軸方向)にはほとんど動かず、方位方向にぐるぐる回るだけで「直接的に陽極までショートカットしない」ことが直感的に確認できます。ホールパラメータ $\Omega_e$ は20 mTで100オーダーに達しており、電子が十分に磁化されている設計条件を満たします。$B_r$ を小さくすれば $\Omega_e$ が1に近づき、電子が衝突で散乱されて加速電場を埋めてしまう「電子バックストリーミング」モードに入ります。

化学推進 vs 電気推進のミッション質量比較

ツィオルコフスキー公式を使って、同じ $\Delta v$ ミッションで両者の質量配分を比較します。

import numpy as np
import matplotlib.pyplot as plt

g0 = 9.81

def mass_ratio(dv, Isp):
    """質量比 m0/mf = exp(dv / (Isp * g0))"""
    return np.exp(dv / (Isp * g0))

def prop_fraction(dv, Isp):
    """推進剤質量比 mp/m0"""
    r = mass_ratio(dv, Isp)
    return 1 - 1/r

dvs = np.linspace(1000, 15000, 200)    # 1 ... 15 km/s
Isps = {
    'Solid (Isp=280s)':       280,
    'LOX/RP-1 (Isp=320s)':    320,
    'LOX/LH2 (Isp=450s)':     450,
    'Hall (Isp=1800s)':      1800,
    'Ion (Isp=3500s)':       3500,
    'Ion (Isp=8000s)':       8000,
}

plt.figure(figsize=(10, 6))
for name, isp in Isps.items():
    plt.plot(dvs/1000, prop_fraction(dvs, isp)*100, lw=2, label=name)
plt.xlabel(r'Mission $\Delta v$ [km/s]')
plt.ylabel('Propellant mass fraction [%]')
plt.title('Propellant fraction vs delta-v for various propulsion systems')
plt.legend(loc='lower right')
plt.grid(True, alpha=0.3)
plt.axhline(90, color='gray', ls=':', alpha=0.7)
plt.text(1.5, 91, '90% propellant — practical upper limit', fontsize=9, color='gray')
plt.tight_layout()
plt.savefig('prop_fraction_comparison.png', dpi=150, bbox_inches='tight')
plt.show()

# はやぶさ相当: dv=4 km/s
dv_hayabusa = 4000
print(f"Δv = {dv_hayabusa/1000} km/s でのpropellant fraction:")
for name, isp in Isps.items():
    print(f"  {name:30s}: {prop_fraction(dv_hayabusa, isp)*100:5.1f}%")

このグラフから、$\Delta v$ が大きくなるほど化学推進と電気推進の差が指数的に拡大することが鮮明に読み取れます。$\Delta v = 4$ km/sのはやぶさ相当ミッションでは、化学推進(Isp=320s)で約72%の推進剤質量を要するのに対し、イオン推進(Isp=3500s)では約11%で済みます。$\Delta v = 10$ km/s級の惑星間ミッションになると、化学推進は推進剤質量比が95%超(乾燥質量の19倍の推進剤)となって事実上実現不可能になる一方、Isp=3500sのイオンエンジンなら25%、Isp=8000sのVASIMR級なら12%と、極めて省推進剤に飛べることがわかります。これが「深宇宙=電気推進」が定着した理由です。

はやぶさ ミッションの $\Delta v$ 予算シミュレーション

最後に、はやぶさが地球出発からイトカワランデブーまで連続加速で稼いだ $\Delta v$ を、scipy.integrate を使って簡略化したロケット方程式を時間積分して再現します。

import numpy as np
from scipy.integrate import solve_ivp
import matplotlib.pyplot as plt

g0 = 9.81

# はやぶさ μ10 諸元(公開資料に基づく代表値)
F_thrust = 8e-3 * 3        # 3基並列運転時の総推力 [N]   ≈ 24 mN
Isp_mu10 = 3000            # 比推力 [s]
ve_mu10  = Isp_mu10 * g0   # 排気速度 [m/s]
m_dry    = 380             # 乾燥質量 [kg]
m_prop0  = 60              # 搭載キセノン [kg]
m0       = m_dry + m_prop0 # 出発時総質量 [kg]
mdot     = F_thrust / ve_mu10   # 質量流量 [kg/s]

print(f"推力 F = {F_thrust*1000:.1f} mN, ve = {ve_mu10/1000:.1f} km/s")
print(f"質量流量 = {mdot*1e6:.2f} mg/s = {mdot*86400*1000:.2f} g/day")

# 連続加速での Δv を時間積分
def dynamics(t, y):
    """y = [v, m].  推進剤がある間だけ推力作用"""
    v, m = y
    if m <= m_dry:
        return [0.0, 0.0]
    dvdt = F_thrust / m
    dmdt = -mdot
    return [dvdt, dmdt]

t_end = m_prop0 / mdot           # 推進剤枯渇時刻 [s]
print(f"推進剤枯渇まで: {t_end/86400:.0f} 日 = {t_end/(86400*365.25):.2f} 年")

sol = solve_ivp(dynamics, (0, t_end), [0.0, m0], dense_output=True,
                rtol=1e-9, atol=1e-12, max_step=86400*30)

t = np.linspace(0, t_end, 500)
v, m = sol.sol(t)

fig, axes = plt.subplots(1, 2, figsize=(12, 5))
axes[0].plot(t/(86400*365.25), v/1000, lw=2)
axes[0].set_xlabel('Time [years]')
axes[0].set_ylabel('Cumulative $\\Delta v$ [km/s]')
axes[0].set_title('Hayabusa-class continuous EP acceleration')
axes[0].grid(True, alpha=0.3)

axes[1].plot(t/(86400*365.25), m, lw=2, color='tab:red')
axes[1].axhline(m_dry, color='gray', ls=':', label=f'dry mass {m_dry} kg')
axes[1].set_xlabel('Time [years]')
axes[1].set_ylabel('Spacecraft mass [kg]')
axes[1].set_title('Mass depletion vs time')
axes[1].legend()
axes[1].grid(True, alpha=0.3)

plt.tight_layout()
plt.savefig('hayabusa_dv.png', dpi=150, bbox_inches='tight')
plt.show()

# ツィオルコフスキー公式との比較
dv_tsiolkovsky = ve_mu10 * np.log(m0 / m_dry)
print(f"最終 Δv (積分結果)        = {v[-1]/1000:.2f} km/s")
print(f"Tsiolkovsky 理論値        = {dv_tsiolkovsky/1000:.2f} km/s")

このシミュレーションから、はやぶさが3基並列のμ10を連続運転して稼げる $\Delta v$ は約 4.4 km/s(2.4年間の連続点火)で、ツィオルコフスキー公式の理論値とぴったり一致することが確認できます。実機ではエンジン故障や運用上の停止があり、実際の往復ミッションでは約 2 年弱の累積点火時間で同等の $\Delta v$ を達成しました。重要なのは「2.4 年も加速し続ける」という時間スケールで、化学推進では数分で終わるマヌーバが電気推進では年単位になることが、低推力ミッション設計の本質を表しています。質量グラフの右側を見ると、最終時点で乾燥質量に到達し、推進剤を使い切った瞬間に加速が止まる様子も再現されています。

ここで実機性能と理論モデルが整合することが定量的に確認できたので、最後にこれらの電気推進が実際の宇宙ミッションでどう活躍しているか、応用事例を整理して締めましょう。

応用 — はやぶさ・Starlink・深宇宙

電気推進の実用例は近年劇的に増えています。代表例を方式別に整理します。

グリッド型イオンエンジンの代表ミッション

• はやぶさ・はやぶさ2(JAXA, 2003/2014) — マイクロ波放電型 μ10/μ20、推力 8 mN/基、$I_{sp}$ ≈ 3000 s。小惑星イトカワ・リュウグウへの往復飛行を実現。サンプルリターンの成否は電気推進の信頼性にかかっていた。
• Deep Space 1(NASA, 1998) — NSTAR(キセノン、推力 92 mN、$I_{sp}$ ≈ 3100 s)。電気推進を主推進系とした初の惑星間ミッション。
• Dawn(NASA, 2007) — NSTAR 3基。小惑星帯のVestaとCeresの両方をランデブーし、軌道間遷移に電気推進の真価を発揮。
• BepiColombo(ESA/JAXA, 2018) — T6 イオンエンジン4基(推力145 mN/基、$I_{sp}$ 4200 s)、太陽方向への巨大減速 $\Delta v$ を稼いで2025年に水星到着予定。
• GOES 静止気象衛星 — 軌道遷移用に電気推進(XR-5A, XIPS-25等)を採用、ロケット上段の代替により打ち上げ質量を大幅削減。

ホール推進機の代表ミッション

• Starlink v2 衛星(SpaceX, 2022〜) — Argon/Krypton ホール推進機を自社開発。コンステレーション軌道調整・最終的なデオービットに使用。低コスト推進剤の選択がコンステ経済性を支える。
• SMART-1(ESA, 2003) — PPS-1350 SPT、推力 70 mN、地球から月への電気推進遷移を世界で初めて実証。
• OneWeb / Galileo / GOES — 多くの商用・準商用衛星がSPT-100/PPS-1350系を南北制御や軌道遷移に採用。

開発中・次世代

• NEXT-C(NASA) — Glenn研で開発中の次世代グリッド型(推力 235 mN、$I_{sp}$ 4170 s)、Psyche ミッション(2023年打ち上げ、小惑星Psycheへ)に搭載済み。
• X3(ミシガン大/NASA) — 100kW級ネステッドチャネルホール推進機。将来の有人火星ミッションのキー技術。
• AEPS(Advanced Electric Propulsion System) — Gateway(月軌道有人ステーション)用12.5kW級ホール推進機。

応用分野の本質

電気推進が選ばれる場面は、共通して「$\Delta v$ が大きいか、ミッション期間が長いか、その両方」です。GEO衛星の南北制御(年間50 m/s × 15年 = 750 m/s)、深宇宙ミッション(数 km/s 以上)、LEOコンステの軌道維持・離脱(累積数十〜数百 m/s)。いずれも「化学推進だと推進剤質量に圧迫されて、ミッションが成立しない/利益が出ない」領域です。逆に、打ち上げや緊急マヌーバのように「秒〜分で大きな $\Delta v$ を出す」場面では化学推進が圧倒的に有利で、両者は補完的に共存しています。

将来的には、超大電力(MW級)を伴う有人火星ミッション、月-火星間の貨物輸送、デブリ除去のような商業用途、さらにはVASIMRや核電気推進(NEP)による太陽系外縁部探査まで、電気推進の活躍領域は広がり続けます。

まとめ

本記事では、宇宙の電気推進技術について、その動作原理から性能比較・実装まで一通り解説しました。

• 電気推進の本質: 電気エネルギーで排気速度 $v_e$ を化学推進の10倍($v_e = 30\sim 100$ km/s, $I_{sp} = 3000\sim 10000$ 秒)に引き上げ、ツィオルコフスキー公式が示す指数的な質量比改善を獲得する。
• グリッド型イオン推進機: 放電室でガスをイオン化し、グリッド間電位差 $V$ で $v_e = \sqrt{2qV/m_i}$ に加速、中和器で電子を供給。空間電荷制限により電流密度は Child-Langmuir則 $j_{\rm CL} = \frac{4\epsilon_0}{9}\sqrt{2q/m_i}\, V^{3/2}/d^2$ で頭打ちになり、推力密度 $F/A = (8\epsilon_0/9) V^2/d^2$ は数N/m$^2$ オーダーに制限される。
• ホール推進機: 半径方向磁場で電子を $\bm E \times \bm B$ ドリフトに閉じ込めて方位ホール電流を作り、軸方向電場でイオンだけ加速。準中性プラズマ全体を加速するため空間電荷制限を回避でき、グリッド型より高推力密度を実現。
• 性能比較: 化学推進は高推力・低$I_{sp}$、電気推進は低推力・高$I_{sp}$。ジェット出力 $P_{\rm jet} = Fv_e/2$ は電源容量に比例するため、電源比質量 $\alpha$ とミッション $\Delta v$ から最適 $I_{sp}$ が決まる固有の設計問題が生じる。
• ミッション設計: 化学推進のインパルシブマヌーバとは異なり、電気推進では連続加速の螺旋軌道を低推力軌道最適化で設計する。はやぶさは2年級の連続点火で4 km/s以上の $\Delta v$ を稼いだ。

電気推進は、もはや「特殊技術」ではなく、深宇宙探査・GEO衛星・LEOコンステの商用運用すべてを支える基盤技術です。化学推進と電気推進は競合関係ではなく、それぞれの得意領域で共存しながら、人類の宇宙活動の射程を着実に広げています。

次のステップとして、以下の記事も参考にしてください。

• ツィオルコフスキーの公式と多段ロケット
• 比推力と推進剤の物理
• 化学推進の燃焼熱力学
• 軌道マヌーバと $\Delta v$
• ガウスの法則(電場)