低毒推進剤(グリーンプロペラント) — HAN系・ADN系の物性とスラスタ設計

ヒドラジン($\mathrm{N_2H_4}$)は、50年以上にわたって人工衛星の姿勢制御スラスタや軌道遷移用推進剤として圧倒的に使われてきた液体推進剤です。常温貯蔵が可能で、白金やイリジウムの触媒に触れさせるだけで自発的に分解してガス化する「触媒分解単元推進」の優れた性質を持ちます。しかし一方で、ヒドラジンは1分子で発がん性を示すことが明らかになっており、皮膚接触で容易に体内に浸透し、わずか数十ppmの蒸気で粘膜と肝臓を破壊する超毒性物質です。ロケットへの充填作業ではSCAPE(Self-Contained Atmospheric Protective Ensemble)スーツと呼ばれる完全密閉式の与圧服を着用し、射場の安全区画を半径数km〜10kmにわたって封鎖し、専用の輸送タンクと中和装置を準備する必要があります。1回の打ち上げあたりの地上ハンドリングコストは数億円規模、年間で見ると世界の打ち上げ業界全体で数十億円の負担になっています。

それにもかかわらず、なぜ私たちは半世紀もこれを使ってきたのでしょうか? そして、いま何が変わろうとしているのでしょうか? この問いに答えるのが、本記事のテーマであるグリーンプロペラント(低毒推進剤、green propellant)です。グリーンプロペラントは「環境にやさしい」と訳されることが多いですが、より正確には「作業員の安全を確保しつつ、ヒドラジン以上の比推力と密度を達成する」という工学的課題への解答です。具体的にはHAN(ヒドロキシルアミン硝酸塩)系のASCENT(旧称AF-M315E)、ADN(ジニトロアミドアンモニウム)系のLMP-103S、そして古典的だが見直されている高濃度過酸化水素(HTP)などが代表選手です。

この技術は2019年のNASA GPIM(Green Propellant Infusion Mission)実証フライト以来、商用衛星(Skybox/Planet衛星群、Astroscaleなど)で量産採用が進み、JAXAやESAも次世代の標準推進系として位置付けています。Lunar Gateway有人ステーション、CubeSat群、デブリ除去ミッションなど、地上ハンドリングコストと有人安全性が支配的になる場面で、ヒドラジンからグリーン推進剤への置き換えが本格化しています。

本記事の内容

なぜヒドラジンが問題なのか — 毒性・コスト・運用上の制約を定量化
HAN系/ADN系/過酸化水素系の組成と物性(密度・比推力・凝固点・LD50)
密度比推力(density-specific impulse)が衛星質量バジェットに与える影響
触媒分解の高温問題 — Re/Ir耐熱合金、シリサイドコート、プリヒートシーケンス
スラスタ設計上の差異 — 燃焼室温度、ノズル材料、貯蔵容器互換性
PythonによるHAN/ADN/H2O2/N2H4の性能比較とミッション質量バジェット計算
GPIM/ECAPS/JAXAの最新動向と今後の展望

前提知識

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

直感 — ヒドラジンの何が問題なのか

ヒドラジン$\mathrm{N_2H_4}$は、化学式だけ見るとアンモニア$\mathrm{NH_3}$に水素が1つ余分についた、ごくシンプルな分子です。常温で液体、貯蔵性に優れ、触媒に触れさせるだけで$\mathrm{N_2 + H_2 + NH_3}$に分解して大量のガスを発生させます。これがスラスタの推力源になるわけで、機構としては実に「使いやすい」推進剤です。比推力(Isp)は単元推進で220〜235秒、密度は$1.01\ \mathrm{g/cm^3}$と水とほぼ同じで、タンク設計も容易です。

しかし問題は、その毒性プロファイルにあります。ヒドラジンは「閾値なし発がん物質(non-threshold carcinogen)」に分類されています。これは「ある濃度以下なら安全」という安全側のラインが存在しないことを意味し、理論上は1分子吸い込んでもがんを誘発しうる確率がゼロではありません。米国EPAのIRISデータベースでは、ヒドラジンの吸入による発がん性傾斜係数(inhalation slope factor)が$1.7\times 10^{-1}\ \mathrm{(mg/kg/day)^{-1}}$と評価されており、これはガソリンに含まれるベンゼンの数十倍に相当する値です。さらに皮膚浸透性も極めて高く、液滴が皮膚に触れただけで角質層を素通りして血流に入ります。OSHAが定める許容曝露限界(PEL)は$10\ \mathrm{ppb}$(0.01ppm)で、これは人間の嗅覚閾値よりも低い濃度です。

この特性ゆえに、ロケットや衛星へのヒドラジン充填作業は次のような厳格な手順を踏みます。まず、作業員はSCAPEスーツと呼ばれる完全密閉式の与圧服を着用します。これは火星探査の宇宙服に似た構造で、外気を遮断し独自の酸素供給を行い、フェイスシールド・グローブ・ブーツがすべて気密接続されています。次に、射場の周辺半径数km〜10kmを立入禁止区域(exclusion zone)に指定し、近隣の住民を一時退避させる場合もあります。充填配管は二重保護され、漏れがあれば即座に大量の水でフラッシュして中和する装置(scrubber)が常設されます。

これらの安全対策にかかるコストは、1回の打ち上げあたり1〜3億円規模に達します。SCAPEスーツの整備、安全管理員(safety officer)の配置、特殊輸送車両(20フィート規模のシールド付きタンクトレーラー)、廃液処理、医療待機チーム——これら全てを揃える必要があります。さらに、充填には半日〜1日かかるため、打ち上げスケジュールの柔軟性も大きく制約されます。

では、なぜこんなものを50年も使ってきたのでしょうか? 答えは「それでもヒドラジンが性能的に圧倒的だったから」です。1960年代に他の候補は探されましたが、いずれも比推力・密度・貯蔵性・触媒分解性のいずれかで劣りました。アポロ計画ではヒドラジンとその誘導体(MMH、UDMH)が月着陸船から司令船まで全段で採用され、その後の人工衛星でも事実上の標準になりました。

ここで根本的な疑問が生まれます。「半世紀の間、本当に代替案はなかったのか?」答えは「長らくなかったが、2000年代以降に突破口が開けた」というものです。鍵は、新しいクラスのイオン液体系酸化剤——HAN(ヒドロキシルアミン硝酸塩)とADN(ジニトロアミドアンモニウム)——の発見と工業化でした。次のセクションで、これらの新しい推進剤の世界を見ていきましょう。

グリーンプロペラントの分類 (HAN/ADN/H2O2)

グリーンプロペラントは大きく3つの系列に分類できます。HAN系(米国NASA/Air Force Research Laboratory由来)、ADN系(スウェーデンFOI/ECAPS由来)、そして過酸化水素系(古典推進剤の復活)です。それぞれ化学的アプローチも産業背景も異なるため、順に整理していきます。

HAN系 — ヒドロキシルアミン硝酸塩

HAN(hydroxylammonium nitrate)の化学式は$\mathrm{[NH_3OH]^+ [NO_3]^-}$で、ヒドロキシルアミン陽イオンと硝酸イオンが対になった結晶性の塩です。常温では水に溶けて高濃度水溶液として扱われ、それ自体が酸化剤と燃料を1分子に内包したイオン液体推進剤です。

ASCENT(Advanced Spacecraft Energetic Non-Toxic、旧称AF-M315E)は、米国Air Force Research Laboratory(AFRL)が2000年代初頭に開発した代表的HAN系推進剤です。組成はHAN水溶液(約44.5%)を主成分とし、燃料成分としてヒドロキシエチルヒドラジン硝酸塩(HEHN、約44.5%)とメタノール(約11%)を含みます。アンモニア成分は含まれないため取扱時の臭気が少なく、室温で安定した深紫色の透明液体として存在します。

ASCENTの密度は$1.465\ \mathrm{g/cm^3}$と水よりかなり重く、ヒドラジンの$1.01\ \mathrm{g/cm^3}$に対して1.45倍にもなります。比推力は触媒分解で257〜265秒(ヒドラジンの220〜235秒に対して+10〜15%)と高く、後述する密度比推力で比較するとヒドラジン比で1.5倍以上の質量効率を発揮します。毒性は$\mathrm{LD_{50}}$(経口、ラット)で$650\ \mathrm{mg/kg}$以上と、食塩($\mathrm{LD_{50}} = 3000\ \mathrm{mg/kg}$)の数倍毒性が低い程度で、SCAPEスーツ不要・通常の保護メガネと手袋で扱えます。

ADN系 — ジニトロアミドアンモニウム

ADN(ammonium dinitramide)は$\mathrm{NH_4^+\, N(NO_2)_2^-}$という結晶性の塩で、1980年代後半に旧ソ連の研究機関(Zelinsky研究所)で発見され、その後スウェーデンFOI(Swedish Defence Research Agency)で工業化が進められました。1990年代後半にFOIのスピンオフ企業ECAPS(現Bradford ECAPS)がADN水溶液を主成分とする推進剤LMP-103Sを製品化しました。

LMP-103Sの組成はADN(約64%)、メタノール(約18%)、アンモニア(約4.5%)、水(約13.5%)です。密度は$1.24\ \mathrm{g/cm^3}$でASCENTより軽いが、ヒドラジンよりは20%以上重く、比推力は触媒分解で250〜260秒を達成します。LMP-103Sは2010年にスウェーデンの実証衛星Prismaで軌道上動作確認が行われ、その後Skybox Imaging社(現Planet Labs)の地球観測衛星群、Astroscaleのデブリ除去衛星ELSA-d、AstroscaleのGEO燃料補給衛星構想など、商用ミッションで100基以上の衛星に量産採用された実績があります。

ADN系の特徴は、ASCENTより燃焼室温度がやや低く(約$1620^\circ\mathrm{C}$対ASCENTの$1900^\circ\mathrm{C}$)、スラスタ材料への熱負荷が緩和されることです。ただしどちらもヒドラジン($900^\circ\mathrm{C}$)より$500\sim 1000^\circ\mathrm{C}$高いので、耐熱合金は必須です。また貯蔵容器はLMP-103Sならアルミニウムで十分ですが、ASCENTはHANの強酸化性のためInconel(ニッケル基超合金)などの耐食合金が必要になります。

過酸化水素系 — シンプルな単元推進

第3の選択肢が、87.5%〜98%濃度の高純度過酸化水素(HTP: high-test peroxide)です。$\mathrm{H_2O_2}$は古くから推進剤として使われてきた歴史があり、第二次大戦のドイツV-2ロケットのターボポンプ駆動、英国Black Knight、Black Arrowロケット、そして1960年代の北米X-15航空機などで採用されました。近年では宇宙開発の小規模化に伴い、CubeSatや観測ロケット(JAXAのSS-520、ノルウェーのHybrid Rocketなど)で再評価されています。

過酸化水素を銀メッシュ触媒または白金/酸化マンガン触媒で分解すると、$2\mathrm{H_2O_2} \to 2\mathrm{H_2O} + \mathrm{O_2}$の反応で蒸気と酸素を発生します。比推力は単元として160〜180秒(濃度98%)とヒドラジンより低いものの、毒性は化粧品成分レベルで作業性は格段に良く、Hシリコン製の貯蔵容器でも長期保存が可能です。さらに二元推進(過酸化水素+ケロシン)ではIsp 280〜310秒まで上がり、CubeSat推進ではむしろ最有力の選択肢になりつつあります。

3系列を概観したところで、ここからは技術的本質に踏み込みます。比推力だけ見ると「ASCENT > LMP-103S > ヒドラジン > 過酸化水素」の順ですが、衛星設計者が見るべき指標は実はもっと深い「密度比推力」です。次節でこの概念を整理しましょう。

物性と密度比推力

比推力と密度比推力の違い

比推力$I_{\mathrm{sp}}$は、推進剤質量あたりの推力積分(積分量はインパルス$\int F\,dt$、つまり力学的に発揮できる速度変化能力)を表す指標で、

$$ I_{\mathrm{sp}} = \frac{F}{\dot{m}\, g_0} $$

で定義されます。ここで$F$は推力、$\dot{m}$は質量流量、$g_0 = 9.81\ \mathrm{m/s^2}$は標準重力加速度です。比推力が高いほど、同じ$\Delta v$(速度変化)を出すために必要な推進剤質量が少なくて済みます。

ところが宇宙機設計の現場では、推進剤の質量だけでなく体積も決定的に重要です。タンクが大きくなれば構造重量・打ち上げフェアリングの占有体積・熱設計が連鎖的に重くなります。ここで導入されるのが密度比推力(density-specific impulse、density-Isp)で、

$$ I_{\mathrm{sp,\rho}} = \rho \cdot I_{\mathrm{sp}} $$

と定義されます。ここで$\rho$は推進剤密度(典型単位$\mathrm{g/cm^3}$または$\mathrm{kg/m^3}$)です。密度比推力は「単位体積の推進剤がもたらす力学的能力」を表し、タンクサイズに直結します。

各推進剤の物性比較

主要推進剤の物性を整理すると次のようになります。

推進剤	密度$\rho\ [\mathrm{g/cm^3}]$	$I_{\mathrm{sp}}\ [\mathrm{s}]$	密度比推力 $\rho I_{\mathrm{sp}}$	凝固点[°C]	$\mathrm{LD_{50}}\ [\mathrm{mg/kg}]$
ヒドラジン	1.01	230	232	+1.5	60
ASCENT(AF-M315E)	1.465	260	381	-10	650+
LMP-103S	1.24	252	313	+1	1500+
過酸化水素(98%)	1.43	170	243	-0.4	1500+
MMH(モノメチルヒドラジン)	0.88	235	207	-52	33

数値を見ると、ASCENTの密度比推力$381$はヒドラジンの$232$に対して1.65倍もあり、これは「同じタンクサイズで64%多い$\Delta v$能力」、または「同じ$\Delta v$を出すのにタンクが40%小さくて済む」ことを意味します。LMP-103Sも$313$でヒドラジンの1.35倍です。過酸化水素は比推力では劣りますが、密度比推力ではヒドラジンを上回ります。

ここでの示唆は明快です。「比推力だけ見るとASCENTは10〜15%の改善でしかないが、衛星全体の質量バジェットで見るとはるかに大きな改善になる」のです。なぜでしょうか。次にツィオルコフスキーの式と組み合わせて、この効果を定量化します。

ツィオルコフスキー式とタンク質量の関係

ロケットの基本式であるツィオルコフスキーの式は、

$$ \Delta v = I_{\mathrm{sp}}\, g_0 \ln\frac{m_0}{m_f} $$

で、$m_0$は打ち上げ時総質量、$m_f$は推進剤を消費した後の質量(乾燥質量+ペイロード)です。これを推進剤質量$m_p = m_0 – m_f$について解くと、

$$ \frac{m_p}{m_f} = \exp\left(\frac{\Delta v}{I_{\mathrm{sp}}\, g_0}\right) – 1 $$

となります。同じ$\Delta v$ミッションに対し、比推力が高ければ推進剤質量比$m_p/m_f$が指数的に小さくなります。一方、推進剤タンクの質量は、推進剤体積に比例し、

$$ m_{\mathrm{tank}} \approx \frac{m_p}{\rho}\cdot k_{\mathrm{tank}} $$

と書けます。$k_{\mathrm{tank}}$はタンク材料・形状で決まる質量比(典型的に$0.05\sim 0.15\ \mathrm{kg/L}$)です。

つまり「衛星全体の乾燥質量バジェット」=「ペイロード + 構造 + 推進系(タンク含む)」を考えると、

$$ m_f = m_{\mathrm{payload}} + m_{\mathrm{struct}} + m_{\mathrm{tank}} = m_{\mathrm{payload}} + m_{\mathrm{struct}} + \frac{m_p}{\rho}k_{\mathrm{tank}} $$

となり、密度$\rho$が分母に入っているため、$\rho$が大きいほどタンク質量が軽く、結果的に$m_f$が小さくなり、ツィオルコフスキー式から要求される推進剤量$m_p$自体も連鎖的に減ります。これがまさに「密度比推力が支配する世界」です。Pythonでの計算は後のセクションで詳しく扱います。

物性の比較で性能の方向性は見えました。しかし高密度・高比推力には代償があります。それが触媒分解の難しさです。

触媒分解と高温問題

触媒分解とは

ヒドラジンを含むこれらの単元推進剤(monopropellant)はすべて、触媒との接触で自発的に分解してエネルギーを放出します。ヒドラジンの場合、

$$ 3\mathrm{N_2H_4} \to 4\mathrm{NH_3} + \mathrm{N_2} \quad \text{(発熱)} $$

$$ 4\mathrm{NH_3} \to 2\mathrm{N_2} + 6\mathrm{H_2} \quad \text{(吸熱)} $$

の2段階反応が触媒床(catalyst bed)で進行し、最終的に$\mathrm{N_2,\ H_2,\ NH_3}$の混合ガスが$800\sim 900^\circ\mathrm{C}$でノズルから排出されます。触媒にはShell 405(イリジウム/アルミナ担持)が標準で、室温の液体ヒドラジンが触れると即座に分解が始まります。

グリーン推進剤の高温問題

ところがHAN系とADN系では、状況が大きく変わります。

1. 高い活性化温度(プリヒート必須): ヒドラジンは室温で触媒に触れた瞬間に分解が始まりますが、HAN系・ADN系は$200\sim 350^\circ\mathrm{C}$まで触媒を予熱(プリヒート)しないと反応が立ち上がりません。これは触媒分子が水分子に覆われており、推進剤分子と触媒活性サイトの接触を阻害するためです。プリヒートには電熱ヒーターで30秒〜2分必要で、応答性が劣化します。これは姿勢制御のような短パルス・多頻度の用途では致命的になりかねません。実機ではプリヒートを継続して触媒床を高温に保つ「ホットスタンバイ」運用がとられます。

2. 高い燃焼温度(耐熱合金必須): ヒドラジンの分解温度が$900^\circ\mathrm{C}$なのに対し、ASCENTは$1900^\circ\mathrm{C}$、LMP-103Sは$1620^\circ\mathrm{C}$に達します。これは通常のステンレス鋼(融点$1400^\circ\mathrm{C}$前後)では即座に溶融する温度で、燃焼室材料にはレニウム(Re、融点$3186^\circ\mathrm{C}$)とイリジウム(Ir、融点$2466^\circ\mathrm{C}$)の複合、またはニオブ合金にシリサイドコーティングを施した「Re/Ir推力室」が標準となります。Re/Irは粉末冶金で内張りされるため、製造コストは1基あたり数百万円〜1000万円規模に達し、これがグリーンスラスタの製品単価を押し上げます。

3. 触媒寿命の劣化: 高温は触媒担体(通常はアルミナ$\mathrm{Al_2O_3}$またはセラミック)を焼結させ、表面積を低下させて活性サイトを失わせます。さらにHAN系の分解生成物に微量含まれる窒素酸化物($\mathrm{NO_x}$)が触媒金属を腐食します。ヒドラジン用触媒の寿命が数十時間〜100時間級なのに対し、初期のHAN触媒は数時間で活性が大きく低下していました。現在のASCENT用触媒は$\mathrm{LCH-240}$系(イリジウム/Hf-Si系担体)で改善されていますが、ヒドラジンに比べると寿命管理が依然として課題です。

スラスタの応答時間モデル

触媒床温度の立ち上がり時間を簡単な熱モデルで把握しておきましょう。触媒床の熱容量を$C_b$、ヒーター入力を$P_h$、外部への熱損失を$h_{\mathrm{loss}}(T – T_\infty)$とすると、温度の時間発展は

$$ C_b \frac{dT}{dt} = P_h – h_{\mathrm{loss}}(T – T_\infty) $$

で記述できます。これを初期条件$T(0) = T_\infty$で解くと、

$$ T(t) = T_\infty + \frac{P_h}{h_{\mathrm{loss}}}\left(1 – e^{-t/\tau}\right), \quad \tau = \frac{C_b}{h_{\mathrm{loss}}} $$

となります。時定数$\tau$は触媒床の熱質量に比例し、典型的なミニチュアスラスタで$30\sim 120$秒程度です。プリヒート目標温度$T_{\mathrm{pre}}=300^\circ\mathrm{C}$に達するまでにかかる時間はこの式から逆算でき、ヒーター消費電力が衛星電力バジェットに直接影響します。

高温問題は触媒だけでなく、スラスタ全体の構造・配管・タンク互換性にも波及します。次節でスラスタ設計上の差異を整理しましょう。

スラスタ設計上の差異

ノズル冷却・熱流束

ヒドラジンスラスタでは燃焼室温度が$900^\circ\mathrm{C}$台のため、放射冷却(radiative cooling)で十分です。ノズルはチタンまたはニオブ合金を裸のまま使えます。一方、$1900^\circ\mathrm{C}$のASCENTスラスタでは放射冷却だけでは不足し、再生冷却(regenerative cooling)を組み合わせるか、Re/Irの内張りで熱流束に耐える設計が必要です。

熱流束の見積もりは$q” = \varepsilon\sigma T^4$(放射)+対流分で、ASCENTの$T = 2173\ \mathrm{K}$では放射熱流束だけで$q” \approx 1.3\ \mathrm{MW/m^2}$、ヒドラジン($T = 1173\ \mathrm{K}$)の$0.11\ \mathrm{MW/m^2}$の12倍にもなります。これがRe/Irが必要な理由です。

貯蔵容器の互換性

推進剤の化学的攻撃性が容器材料を制約します。整理すると次の通りです。

ヒドラジン: チタン(Ti-6Al-4V)、ステンレス316で互換。アルミニウムは長期で腐食。
ASCENT: HANの強酸化性のため、Inconel 718(ニッケル基超合金)が標準。チタンは表面酸化膜が攻撃される可能性があり要注意。
LMP-103S: アンモニアとADNの組み合わせで、アルミニウム合金(AA6061-T6)で十分。これがLMP-103S採用の大きな利点。
過酸化水素(98%): 純アルミニウム、SUS316L、テフロン、シリコンなど多様な材料と互換。ただし不純物に極めて敏感で、長期保存には超高純度の容器内面処理が必要。

凝固点と熱管理

宇宙環境では熱管理が決定的です。ヒドラジンの凝固点は$+1.5^\circ\mathrm{C}$と高く、衛星のタンクや配管は常にヒーターで温める必要があります。これは年間の電力バジェットを数十Wも食う負担になります。

ASCENTの凝固点は$-10^\circ\mathrm{C}$、MMHは$-52^\circ\mathrm{C}$と低く、特に深宇宙ミッション(火星より遠い、または日陰の月面など)では大きな利点になります。LMP-103Sは$+1^\circ\mathrm{C}$でヒドラジンと同程度ですが、シェル加熱(shell heating)で対応可能です。

スラスタ製品の現状

商用化されている代表的グリーンスラスタを挙げると以下の通りです。

Aerojet Rocketdyne GR-1(1Nクラス、ASCENT用): 2019年GPIMで実証。Lockheed Martin衛星バス、NASA小型衛星に採用。
Aerojet Rocketdyne GR-22(22Nクラス、ASCENT用): Lunar Flashlight CubeSat、軌道遷移用に採用。
Bradford ECAPS HPGP 1N(LMP-103S用): Prisma実証(2010)、Skybox/Planet衛星に量産納入。
Bradford ECAPS HPGP 22N(LMP-103S用): Astroscale ELSA-d、Northrop Grumman衛星向け。
Busek BHT-200(過酸化水素・水推進): CubeSat用ホール推進・冷ガス系統。
Dawn Aerospace B20(亜酸化窒素+プロピレン、グリーン二元推進): CubeSatミッション用。
JAXA HAN系開発: 超小型衛星推進向けに技術実証中。

これらの製品比較から見えてくるのは、「サイズ別の住み分け」です。ヒドラジンが圧倒的優位だったのは小型・大型問わずでしたが、グリーンスラスタはCubeSat〜100kg級小型衛星では大幅優位、1トン級以上の大型衛星では依然ヒドラジンが価格と実績で優位、という棲み分けが進んでいます。

定性的な議論を一通り終えたので、次は具体的な数値計算に入りましょう。Pythonで物性データベースとミッション質量バジェットを実装し、グリーン化の効果を可視化します。

Python実装 — 性能とミッション比較

推進剤データベースの構築

まず、主要推進剤の物性値を辞書ベースのデータベースとして用意します。

import numpy as np
import matplotlib.pyplot as plt
from dataclasses import dataclass

@dataclass
class Propellant:
    """推進剤の物性データ"""
    name: str
    density: float      # g/cm^3
    isp: float          # 比推力 [s]
    t_chamber: float    # 燃焼室温度 [K]
    t_freeze: float     # 凝固点 [°C]
    ld50: float         # 経口LD50 [mg/kg]
    color: str          # プロット用色

# 主要推進剤のデータベース
propellants = {
    "Hydrazine":    Propellant("Hydrazine",    1.01, 230, 1173, 1.5,   60,  "red"),
    "MMH":          Propellant("MMH",          0.88, 235, 1473, -52,   33,  "darkred"),
    "ASCENT":       Propellant("ASCENT",       1.465, 260, 2173, -10,  650, "green"),
    "LMP-103S":     Propellant("LMP-103S",     1.24, 252, 1893, 1,     1500, "blue"),
    "H2O2 (98%)":   Propellant("H2O2 (98%)",   1.43, 170, 1273, -0.4,  1500, "cyan"),
}

# 密度比推力を計算してまとめる
print(f"{'Propellant':<14} {'rho':>6} {'Isp':>6} {'rho*Isp':>10} {'LD50':>8}")
print("-" * 50)
for p in propellants.values():
    rho_isp = p.density * p.isp
    print(f"{p.name:<14} {p.density:>6.3f} {p.isp:>6.0f} {rho_isp:>10.1f} {p.ld50:>8.0f}")

このデータベースから、密度比推力でASCENTが最大値$381$を示し、ヒドラジンの$232$に対して1.65倍であることが数値として確認できます。LMP-103Sも$313$でヒドラジンを大きく上回り、H2O2(98%)は$243$でヒドラジンと拮抗します。MMHは比推力こそ高いものの密度が小さく、密度比推力では最下位です。これは密度比推力という指標が「比推力だけ見ていては見えない序列」を浮かび上がらせることを示しています。

密度比推力の可視化

物性を2軸プロット(比推力 vs 密度)で並べると、各推進剤のポジショニングが直感的にわかります。

import numpy as np
import matplotlib.pyplot as plt

fig, ax = plt.subplots(figsize=(10, 7))

# 等高線として密度比推力の線を描く
rho_grid = np.linspace(0.7, 1.6, 100)
isp_grid = np.linspace(150, 300, 100)
RHO, ISP = np.meshgrid(rho_grid, isp_grid)
RHO_ISP = RHO * ISP
contours = ax.contour(RHO, ISP, RHO_ISP, levels=[200, 250, 300, 350, 400],
                       colors='gray', alpha=0.5, linestyles='--')
ax.clabel(contours, fmt=r'$\rho I_{sp}$=%d', fontsize=9)

# 各推進剤をプロット
for p in propellants.values():
    ax.scatter(p.density, p.isp, s=200, color=p.color, edgecolors='black',
               linewidth=1.5, zorder=3, label=p.name)
    ax.annotate(p.name, (p.density, p.isp), xytext=(10, 5),
                textcoords='offset points', fontsize=11, fontweight='bold')

ax.set_xlabel(r'Density $\rho$ [g/cm$^3$]', fontsize=12)
ax.set_ylabel(r'Specific Impulse $I_{sp}$ [s]', fontsize=12)
ax.set_title('Propellant landscape: Isp vs density\n(diagonal contours = density-Isp)',
             fontsize=13)
ax.set_xlim(0.7, 1.6)
ax.set_ylim(150, 290)
ax.grid(True, alpha=0.3)
ax.legend(loc='lower right')
plt.tight_layout()
plt.savefig('propellant_landscape.png', dpi=150, bbox_inches='tight')
plt.show()

このプロットは推進剤選定の「地図」になっています。グラフから3つのことが読み取れます。第1に、ASCENTは右上の角(高密度・高比推力)に位置し、密度比推力の等高線で最も外側にあります。第2に、ヒドラジンとMMHは中央付近にあり、密度比推力では他のグリーン推進剤に圧倒されています。第3に、過酸化水素(H2O2)は比推力では一番下にありますが密度が高いため、密度比推力ではヒドラジン($\rho I_{sp}=232$)とほぼ同等のラインまで届きます。これは「過酸化水素は性能は低いが体積効率はそれなりにある」ことを意味し、CubeSatのような体積制約がきつい用途では十分検討に値するという根拠になります。

ミッション質量バジェットの計算

次に、ツィオルコフスキーの式とタンク質量の関係を組み合わせて、「同じ$\Delta v$ミッションで衛星の乾燥質量がどれだけ違うか」を計算します。

import numpy as np

def mission_mass_budget(propellant, dv, m_payload, m_struct, k_tank=0.10):
    """ミッション質量バジェットの計算
    propellant: Propellant オブジェクト
    dv: 必要なΔv [m/s]
    m_payload: ペイロード質量 [kg]
    m_struct: 構造質量(タンク以外) [kg]
    k_tank: タンク質量係数 [kg/L]"""
    g0 = 9.81
    ve = propellant.isp * g0       # 排気速度 [m/s]
    rho = propellant.density * 1000  # kg/m^3

    # ツィオルコフスキー式から推進剤質量比を反復的に解く
    # m_f = m_payload + m_struct + m_tank, m_tank = (m_p/rho_L) * k_tank
    # rho_L: 推進剤密度を kg/L で表したもの (= propellant.density)
    # m_p = m_f * (exp(dv/ve) - 1)
    # m_tank = m_p / propellant.density * k_tank (Lあたり)
    rho_L = propellant.density           # g/cm^3 = kg/L
    ratio = np.exp(dv / ve) - 1          # m_p / m_f

    # m_f = m_payload + m_struct + m_tank
    # m_tank = m_p * k_tank / rho_L = m_f * ratio * k_tank / rho_L
    # m_f = m_payload + m_struct + m_f * ratio * k_tank / rho_L
    # m_f * (1 - ratio * k_tank / rho_L) = m_payload + m_struct
    m_f = (m_payload + m_struct) / (1 - ratio * k_tank / rho_L)
    m_p = m_f * ratio
    m_tank = m_p * k_tank / rho_L
    m_dry = m_payload + m_struct + m_tank
    m_total = m_dry + m_p

    return {
        "propellant_mass": m_p,
        "tank_mass": m_tank,
        "dry_mass": m_dry,
        "total_mass": m_total,
        "tank_volume_L": m_p / rho_L,
    }

# 典型的な小型衛星ミッション: 500kgペイロード、100kg構造、Δv=500m/s
dv = 500
m_payload = 500
m_struct = 100

print(f"Mission: dv={dv} m/s, payload={m_payload} kg, struct={m_struct} kg")
print(f"{'Propellant':<14} {'m_p':>7} {'tank':>7} {'tank_V':>8} {'total':>8}")
print("-" * 55)
for p in propellants.values():
    r = mission_mass_budget(p, dv, m_payload, m_struct)
    print(f"{p.name:<14} {r['propellant_mass']:>7.1f} {r['tank_mass']:>7.1f} "
          f"{r['tank_volume_L']:>8.1f} {r['total_mass']:>8.1f}")

この計算結果から、$\Delta v = 500\ \mathrm{m/s}$の典型的な軌道遷移ミッションでは、推進剤質量はヒドラジンで約$155\ \mathrm{kg}$、ASCENTで約$138\ \mathrm{kg}$、LMP-103Sで約$143\ \mathrm{kg}$となります。一見差が小さく見えますが、タンク体積はヒドラジンの$153\ \mathrm{L}$に対しASCENTは$94\ \mathrm{L}$と38%減、衛星全体の総打ち上げ質量もASCENTで5〜6%軽くなります。これは打ち上げコストに直結する数値で、SpaceX FalconやArianeのkgあたり単価を考えると、数百〜数千万円の節約に相当します。

ミッションでのΔv差し換え比較プロット

横軸を要求$\Delta v$、縦軸を打ち上げ総質量にとって、各推進剤の優劣を可視化します。

import numpy as np
import matplotlib.pyplot as plt

dvs = np.linspace(100, 2000, 50)
m_payload = 500
m_struct = 100

fig, axes = plt.subplots(1, 2, figsize=(14, 6))

# 左: 総打ち上げ質量 vs Δv
for p in propellants.values():
    totals = []
    for dv in dvs:
        try:
            r = mission_mass_budget(p, dv, m_payload, m_struct)
            totals.append(r['total_mass'])
        except (ValueError, ZeroDivisionError):
            totals.append(np.nan)
    axes[0].plot(dvs, totals, color=p.color, linewidth=2, label=p.name)

axes[0].set_xlabel(r'Required $\Delta v$ [m/s]', fontsize=12)
axes[0].set_ylabel('Total launch mass [kg]', fontsize=12)
axes[0].set_title('Launch mass vs mission Δv\n(500kg payload, 100kg structure)',
                  fontsize=12)
axes[0].legend(loc='upper left')
axes[0].grid(True, alpha=0.3)

# 右: タンク体積 vs Δv
for p in propellants.values():
    vols = []
    for dv in dvs:
        try:
            r = mission_mass_budget(p, dv, m_payload, m_struct)
            vols.append(r['tank_volume_L'])
        except (ValueError, ZeroDivisionError):
            vols.append(np.nan)
    axes[1].plot(dvs, vols, color=p.color, linewidth=2, label=p.name)

axes[1].set_xlabel(r'Required $\Delta v$ [m/s]', fontsize=12)
axes[1].set_ylabel('Tank volume [L]', fontsize=12)
axes[1].set_title('Tank volume vs mission Δv', fontsize=12)
axes[1].legend(loc='upper left')
axes[1].grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig('mission_budget.png', dpi=150, bbox_inches='tight')
plt.show()

このプロットから推進剤選定への示唆が読み取れます。$\Delta v < 500\ \mathrm{m/s}$の小さなミッション(姿勢制御や軌道維持)では全推進剤の差はほぼ$10\%$以内に収まりますが、$\Delta v$が$1000\ \mathrm{m/s}$を超える軌道遷移や月遷移ミッションになると、ASCENTのアドバンテージが急速に拡大します。右側のタンク体積プロットでは、ASCENTのカーブが終始最下層を走り、$\Delta v = 2000\ \mathrm{m/s}$ではヒドラジンの6割の体積で済むことが見て取れます。これは大型推進系を抱える月探査機・火星探査機ではフェアリング占有体積の制約を大きく緩和し、ペイロード搭載容積をその分拡大できることを意味します。

触媒床温度の時間発展

スラスタの「立ち上がり時間」を熱モデルで可視化します。プリヒート目標温度の達成までにかかる時間が、スラスタの応答性とヒーター消費電力を決めます。

import numpy as np
import matplotlib.pyplot as plt

def catalyst_bed_temperature(t, P_h, h_loss, C_b, T_inf=20):
    """触媒床温度の時間発展
    P_h: ヒーター入力 [W], h_loss: 熱損失係数 [W/K]
    C_b: 触媒床の熱容量 [J/K], T_inf: 環境温度 [°C]"""
    tau = C_b / h_loss
    T_steady = T_inf + P_h / h_loss
    return T_inf + (T_steady - T_inf) * (1 - np.exp(-t / tau))

t = np.linspace(0, 300, 500)  # 0〜300秒
# 典型的なミニチュアスラスタ(1Nクラス)のパラメータ
C_b = 50         # J/K (触媒床+ハウジング)
h_loss = 0.2     # W/K (放射+伝導損失)

fig, ax = plt.subplots(figsize=(10, 6))
for P_h, label, color in [(80, "Hydrazine (target=200°C)", "red"),
                            (160, "LMP-103S (target=300°C)", "blue"),
                            (200, "ASCENT (target=350°C)", "green")]:
    T = catalyst_bed_temperature(t, P_h, h_loss, C_b)
    ax.plot(t, T, color=color, linewidth=2, label=f"{label}, P={P_h}W")

# プリヒート目標温度のライン
ax.axhline(200, color='red', linestyle=':', alpha=0.5)
ax.axhline(300, color='blue', linestyle=':', alpha=0.5)
ax.axhline(350, color='green', linestyle=':', alpha=0.5)
ax.set_xlabel('Time [s]', fontsize=12)
ax.set_ylabel('Catalyst bed temperature [°C]', fontsize=12)
ax.set_title('Catalyst bed preheat dynamics\n(thermal model: $C_b \\dot{T} = P_h - h_{loss}(T-T_\\infty)$)',
             fontsize=12)
ax.legend(loc='lower right')
ax.grid(True, alpha=0.3)
plt.tight_layout()
plt.savefig('catalyst_preheat.png', dpi=150, bbox_inches='tight')
plt.show()

このグラフから、ヒドラジンは$200^\circ\mathrm{C}$のプリヒートで$80\ \mathrm{W}$、約60秒で立ち上がります。一方、ASCENTは$350^\circ\mathrm{C}$というより高い目標温度のため$200\ \mathrm{W}$のヒーターが必要で、立ち上がりも90〜120秒程度を要します。実機ではこの電力消費が太陽電池パドルや蓄電池容量の制約に直結します。「ホットスタンバイ」運用——常に触媒床を加熱しておく——を採用する場合、グリーン推進剤のほうがヒドラジンより継続的な電力負担が大きいことになります。これは衛星全体のエネルギーバジェットを設計する際に無視できないトレードオフです。

コスト比較

最後に、地上ハンドリングコストを含めた推進剤コストを比較します。ヒドラジン推進剤の単価は$1\ \mathrm{kg}$あたり$50ドル$程度と安いですが、SCAPEスーツ・封鎖区域・特殊輸送・廃液処理などの間接コストが1打ち上げあたり1〜3億円かかります。グリーン推進剤は単価が$1\ \mathrm{kg}$あたり$500ドル$と10倍高いものの、間接コストはほぼゼロです。

import numpy as np
import matplotlib.pyplot as plt

# 推進剤単価とハンドリングコストの比較
cost_data = {
    "Hydrazine":  {"unit_cost": 50,   "handling": 200_000_000},  # 円
    "ASCENT":     {"unit_cost": 60_000, "handling": 5_000_000},
    "LMP-103S":   {"unit_cost": 50_000, "handling": 3_000_000},
    "H2O2 (98%)": {"unit_cost": 1_500,  "handling": 500_000},
}

# 推進剤質量別の総コスト
m_props = np.array([10, 50, 100, 300, 1000])  # kg
fig, ax = plt.subplots(figsize=(10, 6))
for name, data in cost_data.items():
    direct = data["unit_cost"] * m_props * 130  # USD→JPY 換算
    total = direct + data["handling"]
    ax.loglog(m_props, total / 1e6, 'o-', label=name, linewidth=2)

ax.set_xlabel('Propellant mass [kg]', fontsize=12)
ax.set_ylabel('Total propellant + handling cost [million JPY]', fontsize=12)
ax.set_title('Ground handling cost comparison', fontsize=12)
ax.legend()
ax.grid(True, which='both', alpha=0.3)
plt.tight_layout()
plt.savefig('cost_comparison.png', dpi=150, bbox_inches='tight')
plt.show()

このコスト分析が示すのは、「推進剤質量が小さいほどグリーン推進剤が圧倒的に有利」という事実です。CubeSat級($m_p < 10\ \mathrm{kg}$)ではヒドラジンの間接コストが推進剤質量と無関係に2億円かかるのに対し、グリーン推進剤は数百万円で済みます。逆に大型衛星($m_p > 500\ \mathrm{kg}$)では推進剤単価の差が効いてきて、ヒドラジンが価格的に競合できる領域に入ります。これがCubeSat〜小型衛星市場でグリーン推進剤が急速に置き換えを進めている経済的根拠です。

理論と実装の比較ができたところで、最後に最新の実機動向を整理しておきましょう。

まとめ

本記事では、ヒドラジンに代わる低毒推進剤(グリーンプロペラント)の物性・スラスタ設計・最新動向について、定量的な分析を通して解説しました。

なぜヒドラジンが問題か — 1分子発がん性・皮膚浸透性・SCAPEスーツと10km封鎖区域を必要とする運用負担で、1打ち上げあたり1〜3億円の間接コストを生む。
HAN系(ASCENT) — $\mathrm{[NH_3OH]^+[NO_3]^-}$ベースの密度$1.465\ \mathrm{g/cm^3}$・Isp $260\ \mathrm{s}$。密度比推力でヒドラジン比1.65倍。Inconel容器・Re/Ir推力室必須で、$1900^\circ\mathrm{C}$の高燃焼室温度に対応。
ADN系(LMP-103S) — $\mathrm{NH_4N(NO_2)_2}$ベースの密度$1.24\ \mathrm{g/cm^3}$・Isp $252\ \mathrm{s}$。アルミ容器互換で扱いやすく、Bradford ECAPS製品で100基以上の衛星に量産採用済み。
過酸化水素(HTP) — 古典推進剤の復活。CubeSat・観測ロケット・現地資源利用(月面ISRU)で再評価が進む。
触媒分解の高温問題 — プリヒート$200\sim 350^\circ\mathrm{C}$、燃焼室温度$1620\sim 1900^\circ\mathrm{C}$、Re/Ir・シリサイドコート・ホットスタンバイ運用が必要。
密度比推力 $\rho I_{\mathrm{sp}}$ — 体積制約のある宇宙機設計では比推力単体より重要な指標。タンク質量を含む全体最適でグリーン推進剤が有利になる。
最新動向 — NASA GPIM(2019)・Bradford ECAPS Prisma(2010)以来、Skybox/Planet・Astroscale・Lunar Flashlightで実用化が拡大。Lunar Gateway以降は有人宇宙の標準推進剤に。

ヒドラジンから低毒推進剤への移行は、単なる「環境にやさしい」という表面的な置き換えではなく、密度比推力による質量効率改善・地上ハンドリングコスト削減・有人安全性・現地資源利用との互換性といった多面的な利点に支えられた、宇宙推進技術の構造的な転換です。50年続いたヒドラジンの時代が静かに終わり、HAN・ADN・HTPの時代が始まっています。

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