ツィオルコフスキーの公式で比推力・質量比からΔvをリアルタイム計算。LH2/LOXからケロシン・固体燃料まで推進剤を切り替えて性能を比較しよう。
推力：$F = \dot{m}\cdot v_e + (P_e - P_a) A_e$
排気速度：$v_e = I_{sp}\cdot g_0$
$g_0 = 9.80665\ \mathrm{m/s^2}$
ロケットの運動量保存則から導かれる、速度増分Δvを計算する基本公式です。ロケットが推進剤を後方に噴射することで得られる前方への速度変化を表します。
$$\Delta v = I_{sp}\cdot g_0 \cdot \ln\!\left(\frac{m_0}{m_f}\right)$$$\Delta v$: 速度増分 [m/s]
$I_{sp}$: 比推力 [s]（推進剤の効率）
$g_0$: 地球表面の重力加速度 (9.80665 m/s²)
$m_0$: 打ち上げ初期質量（機体+推進剤）
$m_f$: 燃焼終了後の質量（機体のみ）
ロケットが発生する推力の計算式です。運動量変化とノズル出口での圧力差による推力の和で表されます。
$$F = \dot{m}\cdot v_e + (P_e - P_a) A_e$$$F$: 推力 [N]
$\dot{m}$: 推進剤質量流量 [kg/s]（シミュレーターのパラメータ）
$v_e$: 排気速度 [m/s] ($v_e = I_{sp}\cdot g_0$)
$P_e$: ノズル出口圧力
$P_a$: 周囲（大気）圧力
$A_e$: ノズル出口面積
*真空中($P_a=0$)では第2項が最大となり、推力が増加します。
打ち上げロケットの段数設計：地球低軌道(LEO)投入に必要な約9.4 km/sのΔvを、1段ロケットで達成するには非現実的な質量比が必要です。そのため、燃え尽きたタンクやエンジンを切り離す「多段式ロケット」が用いられ、各段の最適な質量比と$I_{sp}$がこの式で設計されます。
深宇宙探査機の軌道設計：火星や小行星への探査機は、地球軌道離脱、目標天体への軌道投入、着陸など、複数のΔvを必要とします。ミッション全体で必要なΔvの合計から、搭載すべき推進剤の総量がツィオルコフスキーの公式で決定されます。
推進剤の選定：高$I_{sp}$の液体水素/酸素は上段エンジンに、高密度で推力が大きくタンクがコンパクトにできるケロシン/酸素は離床に必要な大推力が求められる第1段エンジンに、といった使い分けの根拠となります。
着陸機の逆噴射設計：月や火星に軟着陸するためには、重力に打ち勝ち減速するための逆噴射（Δv）が必要です。着陸機の重量とエンジン性能から、必要な燃料量と燃焼時間が計算されます。
このシミュレーターを使い始めるとき、いくつか気をつけてほしいポイントがあるよ。まず、「Δvはロケットの最終速度そのものではない」ってこと。Δvはエンジンが生み出せる「速度増分の能力」だ。実際の打ち上げでは、重力や空気抵抗でこのΔvはどんどん「消費」されちゃうんだ。例えば、地上から静止軌道まで行くには理論上約10km/sのΔvが必要だけど、重力損失や抗力損失を考慮すると、ロケットは合計で13km/s近くのΔv能力を持っていないと到達できない。シミュレーターの数字は「理想的な宇宙空間での値」と覚えておこう。
次に、比推力(Isp)と推力はトレードオフの関係にあることが多いってこと。Ispが高い液体水素エンジンは燃費はいいけど、密度が低い燃料を大量に貯める巨大なタンクが必要で、構造質量が増えがち。逆に、Ispは低くても推力が大きく密度の高いケロシンは、離床時のパンチ力に優れる。月面着陸のように精密な推力調整が必要な場面では、また別のエンジン特性が求められる。一つの数値だけを見て「優劣」を決めつけないようにね。
最後に、「推重比」は時間とともに変化するってことを忘れがち。打ち上げ直後は燃料満タンで機体が一番重いから推重比は最小。燃料が減るにつれて機体は軽くなり、推重比はどんどん上がっていく。シミュレーターで「推力」を固定値で考えがちだけど、実際の設計では「初期推重比が1.3以上あるか」と「最終段の推重比が大きくなりすぎて乗員に過負荷がかからないか」の両方をチェックするんだ。
このツールで扱う計算は、ロケット工学の根幹だけど、実は様々な工学分野と深く結びついているんだ。まず間違いなく関連するのが「材料工学」だね。高比推力を実現するには燃焼室やノズルが超高温度・高圧に耐えないといけない。ここで活躍するのが銅合金や炭素複合材などの先進材料。タンクを軽くして質量比を改善するには、アルミリチウム合金や複合材の開発が不可欠だ。
次に「熱流体力学・燃焼工学」。ノズル内での超音速流れの解析や、推進剤の混合・燃焼効率をシミュレーションするCFD（数値流体力学）は、Ispを理論値に近づけるための必須技術だ。排気速度$v_e$を上げるには、ノズルの形状設計がめちゃくちゃ重要になる。
もう一つ見落とせないのが「制御工学」との関係。推重比が時間変化する中で、ロケットの姿勢を安定させ、所定の軌道に正確に投入するには高度な制御システムが必要。例えば、複数エンジンの推力調整（スロットリング）やジンバル機構による推力方向制御は、Δvを「有効に使う」ための技術と言える。さらに、電気推進（イオンエンジン）のような非化学ロケットを学べば、電気工学・プラズマ物理学の領域にも広がっていくよ。
ツィオルコフスキーの公式に慣れてきたら、次は「ロケットの運動方程式そのものを追いかけてみよう」。このシミュレーターの基礎にあるのは、質量が減少する物体の運動を表す「変質量系の運動方程式」だ。外力（重力や抗力）を考慮すると、次のような式になる。
$$ m(t) \frac{dv}{dt} = F - m(t)g - D $$
ここで、$m(t)$は時間とともに減る質量、$g$は重力加速度、$D$は空気抵抗だ。この微分方程式を（単純化した条件下で）積分すると、あの対数関数の公式が導かれる。この過程を理解すれば、なぜΔvが「理想値」なのかが腑に落ちるはず。
学習のステップとしては、まず「パーキング軌道」「ホーマン移軌道」「重力ターン」といった軌道力学の基本概念を押さえよう。地球周回から月、他の惑星へ行くのに必要なΔvをまとめた「Δvマップ」を見ると、このシミュレーターで遊んだ各数値の実践的な意味が一気にわかる。例えば、地球低軌道から月着陸までには約6km/sの追加Δvが必要、といった感じだ。
最後に、「多段式ロケットの最適化」に挑戦してみるのがおすすめ。例えば、合計の質量が同じでも、第1段と第2段にどう質量（燃料と構造）を配分するかで、得られる総Δvは大きく変わる。これは「等分配則」として知られる最適化問題で、工学と数学が交差する面白いテーマだ。このシミュレーターのパラメータを「段」ごとに分けて考えてみると、また新しい発見があるよ。