アポロ・スペースシャトル・スターシップのプリセットで突入加熱を比較。速度・角度・弾道係数を変えて最大G・熱流束・壁面温度を探ろう。
$\rho(h) = \rho_0 e^{-h/H}$（H=7000 m）
$\dot{q}\approx C_h \rho^{0.5}V^3 R_n^{-0.5}$
平衡壁温: $\varepsilon\sigma T_w^4 = \dot{q}$
大気密度は高度とともに指数関数的に減少します。ここでは簡便な「指数大気モデル」を使っています。
$$\rho(h) = \rho_0 e^{-h/H}$$$\rho$: 高度hにおける大気密度 [kg/m³], $\rho_0$: 地表での基準密度, $H$: スケールハイト (約7000 m)。高度がH上がると密度は1/e倍になります。
機体表面（特にノーズ）に加わる熱流束は、密度の平方根と速度の3乗に比例します。これは実用的な近似式です。
$$\dot{q}\approx C_h \rho^{0.5}V^3 R_n^{-0.5}$$$\dot{q}$: ステグナント熱流束 [W/m²], $C_h$: 熱伝達係数, $V$: 速度 [m/s], $R_n$: ノーズ先端半径 [m]。速度が2倍になると熱は8倍になるので、超高速の恐ろしさがわかります。
表面が十分な時間熱にさらされると、入ってくる熱と放射で宇宙空間に逃がす熱が平衡状態に達します。その時の壁面温度は以下の式で求められます。
$$\varepsilon\sigma T_w^4 = \dot{q}$$$\varepsilon$: 表面の放射率, $\sigma$: ステファン・ボルツマン定数, $T_w$: 平衡壁面温度 [K]。放射率εを上げると（黒くすると）、同じ熱流束でも温度を下げられることがシミュレーターで確認できます。
有人宇宙船の帰還：アポロ計画やソユーズ、クルードラゴンなど、地球帰還時の熱防護システム（TPS）設計の基本計算に用いられます。どの角度で突入するか（再突入回廊）は、この熱流束と加速度（最大G）の計算から決定されます。
再利用型宇宙機の設計：スペースシャトルやスペースXのスターシップのように、何度も再突入を繰り返す機体では、熱防護タイルや耐熱材の寿命予測に空力加熱の解析が不可欠です。シミュレーターの「スターシップ」プリセットはその一例です。
惑星探査機の着陸：火星や金星など、大気を持つ惑星への探査機突入でも同じ物理が働きます。大気組成や密度が異なるため、パラメータを調整した計算が必要になります。
極超音速ミサイルの開発：極超音速滑空体（HGV）など、長時間大気圏内を極超音速で飛行する兵器では、先端部の継続的な空力加熱と材料の耐熱性が最大の課題の一つとなっています。
このツールを使い始めるときに、特に気をつけてほしいポイントがいくつかあるよ。まず「熱流束が大きい場所＝機体が一番熱い場所」とは限らないってこと。確かにノーズ先端には猛烈な熱流束がかかる。でも、熱が流れ込む「強さ」は熱流束だけど、部品が実際にどの温度まで上がるかは、その材料の熱容量や放熱のしやすさ（熱伝導率）で決まるんだ。例えば、熱流束が高いけど熱を素早く内部に伝えたり、後ろに捨てたりできる構造なら、表面温度は意外と低く抑えられる。逆に、断熱性の高い素材だと、熱がこもってじわじわ温度が上がり、内部機器がダメになることもある。シミュレーターの「平衡壁面温度」は、あくまで「十分な時間が経って、入ってくる熱と放射で出ていく熱が釣り合った状態」の理論値だからね。
次に、パラメータ設定の落とし穴。「速度」と「高度」は独立した変数ではないんだ。実際の再突入では、高度が下がるにつれて空気抵抗で速度がガクンと落ちる。このツールではシンプルにするために、ある瞬間の「条件」を入力して計算しているけど、実務では「軌道」全体を通した時間変化を追う「軌道計算」が必須だ。例えば、高度70kmで秒速7kmの状態と、高度40kmで秒速7kmの状態では、大気密度が100倍以上違うから、熱流束も全然違う値になる。ツールで遊ぶときは、高度と速度の組み合わせを現実的に考えるクセをつけよう。例えば「アポロ」プリセットの初期値（高度120km、速度11km/s）からスタートするのがおすすめだ。
最後に、この計算は「局所」の評価だという点。機体全体の温度分布や、熱が構造内部にどう伝わっていくか（熱伝導解析）は、もっと複雑なCAEソフト（熱流体連成解析：CHT）の領域だ。このツールは、TPS（熱防護システム）を設計する第一歩として「機体のどこが一番熱的に厳しいか」を見極める「スクリーニング」に使うイメージだよ。
この空力加熱計算の考え方は、実は宇宙機の再突入だけじゃなくて、いろんな高速飛翔体の設計に応用されているんだ。まず真っ先に挙がるのが極超音速機（マッハ5以上）の開発だ。例えば、スクラムジェットエンジンを搭載した実験機や、将来の極超音速旅客機。大気圏内を極超音速で巡航するときも、機首や翼の前縁には再突入に匹敵する空力加熱が発生する。ここで学んだ熱流束の式 $\dot{q} \propto V^3$ は、速度の怖さを如実に物語っているよね。
もう一つは、ロケットの上昇時の空力加熱だ。再突入とは逆に、ロケットが打ち上げで大気圏を抜けるときも、最大動圧（Max Q）付近で機体表面に熱が加わる。特に先端のフェアリングや翼の部分だ。また、隕石やスペースデブリの大気圏突入を研究する「天体物理学」や「スペースデブリ対策」の分野でも、同じ物理モデルが基礎として使われている。隕石が燃え尽きずに地上に落ちてくるかどうかは、まさにこのツールで扱っている「弾道係数」と「突入角度」で大きく変わるんだ。
さらに視野を広げると、プラズマ物理学にも繋がる。速度が非常に高いと、圧縮された空気の分子が分解・電離してプラズマ状態（イオン化したガス）になる。これが電波を遮断する「通信ブラックアウト」現象の原因だ。このツールで計算する高温環境は、そのプラズマ層の発生条件を考える入り口にもなるんだ。
このツールに慣れて「もっと詳しく知りたい！」と思ったら、次のステップに進んでみよう。まずは「運動方程式」を時間発展で解くことを理解するのが大事だ。このツールの計算は瞬間風速的だけど、実現象は連続的だ。再突入体の運動は、重力、空気抵抗（抗力）、そして揚力（もしあれば）のバランスで決まる。これを数値積分（例えばオイラー法）で解けば、高度や速度の時間変化、つまり「軌道」が計算できる。ExcelやPythonで、重力加速度 $g$ と抗力 $D = 0.5 \times \rho \times V^2 \times C_d \times A$ を使った簡単なプログラムを組んでみるのが最高の練習だ。
次に、大気モデルを精緻化しよう。このツールで使っている指数関数的なモデルは簡便だが、実際の大気は高度によって温度や組成が変わる。より現実的な「USスタンダード大気」などのデータをプログラムに組み込むと、計算精度が一段上がる。また、熱流束の式も、より複雑だが精度の高い「フェイ-リデル式」などを学ぶと、実務の解析手法に近づける。
最終的には、専門書や論文で「熱防護システム（TPS）の材料科学」を学ぶことをおすすめする。アブレータ（焼け落ちて熱を逃がす材料、アポロで使用）と、繰り返し使える放熱型タイル（スペースシャトルで使用）の根本的な違いは何か。スターシップで採用されているステンレス鋼の外皮は、どのように熱を放射しているのか。空力加熱の「量」を計算できるようになったら、次はその「熱」とどう戦うかという、もう一つの深い技術の世界が待っているんだ。