パラメータ
$\dfrac{T}{T_0}=\left(1+\dfrac{\gamma-1}{2}M^2\right)^{-1}$
$\dfrac{P}{P_0}=\left(\dfrac{T}{T_0}\right)^{\gamma/(\gamma-1)}$
垂直衝撃波後マッハ数:
$M_2=\sqrt{\dfrac{M_1^2+\frac{2}{\gamma-1}}{\frac{2\gamma}{\gamma-1}M_1^2-1}}$
マッハ数と比熱比γから、圧力比・温度比・密度比・面積比を即座に算出。垂直衝撃波後の状態量も計算。ジェットエンジン・超音速ノズル設計の基礎。
等エントロピー流れでは、流れの全状態(滞止状態)は保存され、静的な状態(温度T, 圧力P, 密度ρ)はマッハ数Mと比熱比γの関数として以下のように表されます。
$$ \frac{T}{T_0}= \left(1 + \frac{\gamma - 1}{2}M^2 \right)^{-1}$$ $$ \frac{P}{P_0}= \left(1 + \frac{\gamma - 1}{2}M^2 \right)^{-\gamma/(\gamma-1)}$$ $$ \frac{\rho}{\rho_0}= \left(1 + \frac{\gamma - 1}{2}M^2 \right)^{-1/(\gamma-1)}$$ここで、$T_0, P_0, \rho_0$は流速がゼロの時の滞止温度、滞止圧力、滞止密度です。$M$はマッハ数、$\gamma$は比熱比(空気では約1.4)です。流れが加速すると(Mが増えると)、静圧と静温は低下することが式から読み取れます。
流路面積とマッハ数の関係は、連続の式と等エントロピー関係から導かれ、ラバルノズルの設計に直接使われる重要な式です。また、垂直衝撃波を通過する際の状態変化はランキン・ユゴニオの関係式で記述されます。
$$ \frac{A}{A^*}= \frac{1}{M}\left[ \left( \frac{2}{\gamma+1}\right) \left( 1 + \frac{\gamma-1}{2}M^2 \right) \right]^{(\gamma+1)/(2(\gamma-1))}$$ $$ M_2^2 = \frac{M_1^2 + \frac{2}{\gamma-1}}{\frac{2\gamma}{\gamma-1} M_1^2 - 1} $$$A$は任意断面の面積、$A^*$はマッハ数が1となる臨界断面(スロート)の面積です。$M_1$は衝撃波上流のマッハ数、$M_2$は下流のマッハ数です。衝撃波ではエントロピーが増加するため、全圧は減少します。
航空宇宙エンジンのノズル設計:ロケットエンジンやジェットエンジンの排気ノズルは、燃焼ガスを効率よく加速・膨張させ推力を得るために設計されます。ノズル形状(特にA/A*)は目標の排気マッハ数を達成するために、等エントロピー流れの関係式を用いて決定されます。
超音速機・再突入カプセルの空力設計:機体周りの流れ場には衝撃波が発生します。特にエンジンへの空気取り入れ口(インテーク)では、衝撃波の位置と強度を制御してエンジンが効率よく働くように設計されます。本ツールの衝撃波計算はその基礎理解に役立ちます。
CAE(数値流体力学)解析の前処理・検証:圧縮性流体を扱うCFD解析を行う際、計算結果の妥当性をチェックするための参照解として、等エントロピー関係や垂直衝撃波の理論解が使われます。複雑な形状の解析を始める前に、本ツールで単純な流路の挙動を確認するのは有効です。
風洞実験データの解釈:風洞実験で計測された静圧・静温のデータから、流れのマッハ数や滞止状態を推定する際に、等エントロピー関係式が逆に用いられます。実験と理論を結びつける基本的なツールとなっています。
この計算機を使い始めるとき、いくつか気をつけてほしいポイントがあるよ。まず、「等エントロピー」は「等温」や「等圧」ではないってこと。摩擦や熱の出入りがない「理想的な」流れを仮定しているんだ。例えば、実際のノズル壁面には摩擦があるから、計算結果は「理論上の上限値」と考えるのが実務的だね。次に、比熱比γの値は流体で変わることを忘れがち。空気(1.4)で計算した結果を、そのまま燃焼ガス(γ=1.2〜1.3程度)に適用すると、圧力や温度の見積もりが大きく外れるから注意してね。
もう一点、面積比「A/A*」について。これは同じ流路の中での、ある断面とスロート断面の比だ。別々のノズル同士の面積を単純比較する値じゃないんだ。例えば、A/A*=2の流路とA/A*=4の流路があったら、後者の方が同じ比熱比ならより高いマッハ数になる、という関係を表している。最後に、衝撃波計算で「上流マッハ数M1」を入力するとき、M1は必ず1より大きい値にしないと意味がないよ。衝撃波は超音速流れで発生する現象だからね。M1=2で計算すると、M2は約0.58になるはずだ。これが1以下にならない場合は、入力を見直してみよう。
この等エントロピー流れと衝撃波の計算は、CAEの世界を超えて、さまざまな先端工学の根幹を支えているんだ。まず航空宇宙工学では、ロケット・ノズルの形状最適化や、超音速旅客機・再使用型宇宙往還機のインテーク設計に直結する。インテークでは、外部で発生する衝撃波の位置を制御して、エンジンへと流れを効率よく送り込むことが求められるよ。
次にターボ機械工学だ。ガスタービンの静翼・動翼間の流路は、ミニチュアのノズルや拡大管の連続だと考えられる。ここでの流れの加速・減速と損失の評価に、等エントロピー変化が基準として使われる。また、自動車工学でも、ターボチャージャーのタービンやコンプレッサー、さらにはF1マシンのエンジンインテークの設計では、圧縮性流体の考え方が欠かせない。最近では超音速燃焼(スクラムジェット)エンジンの研究でも、超音速流路内での燃料混合や燃焼を理解するための第一歩として、この基礎計算が活用されているんだ。
このツールで遊んで感覚をつかんだら、次は「なぜそうなるのか」の理論的背景を学ぶと、一気に視界が開けるよ。おすすめのステップはまず、熱力学の第一法則と第二法則を復習すること。等エントロピーの仮定が、いかに強力な制約(断熱かつ可逆)を意味するかがわかるはずだ。その上で、ベルヌーイの式の圧縮性流体版と言える、エネルギー保存式($$h + \frac{V^2}{2} = const.$$)を導出してみよう。ここでhはエンタルピーだ。
数学的には、ツールで使われている数式は、連続の式・運動量保存式・エネルギー保存式・状態方程式を連立させて導出されている。特に面積比の式は少し複雑だが、マッハ数Mを独立変数として、他の変数をMで微分していく過程で現れるものだ。この「微分形式」を追うことで、流路面積が変化することがなぜ加速・減速を引き起こすのか、その本質が理解できる。次のトピックとしては、斜め衝撃波や膨張波(プラントル・マイヤー流れ)に進むのがいいね。これらを学べば、超音速機周りの複雑な流れ場を、衝撃波と膨張波の組み合わせとして読み解く、という面白い世界が見えてくるんだ。