等エントロピー流れ計算機 戻る
圧縮性流体

等エントロピー流れ計算機

マッハ数と比熱比γから、圧力比・温度比・密度比・面積比を即座に算出。垂直衝撃波後の状態量も計算。ジェットエンジン・超音速ノズル設計の基礎。

パラメータ

T/T₀
温度比
P/P₀
圧力比
A/A*
面積比
M₂ (衝撃波後)
下流マッハ数
等エントロピー関係式
$\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}}$

等エントロピー流れ計算機とは

🧑‍🎓
「等エントロピー流れ」って何ですか?教科書で見たけど、具体的にどんな時に使うんですか?
🎓
ざっくり言うと、摩擦や熱の出入りが無視できる理想的な圧縮性流体の流れだよ。例えば、ロケットのノズルの中や、ジェットエンジンの吸気口の設計の基礎になる考え方なんだ。このシミュレーターで、上の「マッハ数」スライダーを動かしてみて。流れが加速・減速するにつれて、温度や圧力がどう変わるかがリアルタイムでわかるよ。
🧑‍🎓
え、そうなんですか!「A/A*」って表示されてますけど、これって何の面積比なんですか?
🎓
良いところに気が付いたね。A/A*は、ある断面の流路面積と、流速がちょうど音速(マッハ数=1)になる「のど」(スロート)の面積の比だ。超音速ノズル(ラバルノズル)の設計では、この値がマッハ数を決める超重要なパラメータになるんだ。右の模式図を見ながら、マッハ数を変えるとA/A*がどう変わるか確かめてみよう。
🧑‍🎓
なるほど!じゃあ「衝撃波タイプ」のボタンは?衝撃波を入れると、何が変わるんですか?
🎓
衝撃波は、流れが急に圧縮される超薄い領域だ。例えば、超音速機が飛行する時、機首や翼の前によく発生するよ。「垂直衝撃波あり」に切り替えて、衝撃波前後のマッハ数M1, M2を見てごらん。衝撃波を通過すると、流れは必ず亜音速になるのがわかる。実務では、エンジンインテークの設計で、どこに衝撃波を立てるかを考えるんだ。

物理モデルと主要な数式

等エントロピー流れでは、流れの全状態(滞止状態)は保存され、静的な状態(温度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で微分していく過程で現れるものだ。この「微分形式」を追うことで、流路面積が変化することがなぜ加速・減速を引き起こすのか、その本質が理解できる。次のトピックとしては、斜め衝撃波膨張波(プラントル・マイヤー流れ)に進むのがいいね。これらを学べば、超音速機周りの複雑な流れ場を、衝撃波と膨張波の組み合わせとして読み解く、という面白い世界が見えてくるんだ。