マッハ数・圧縮性流れ計算機 戻る
空気力学 / 圧縮性流体

マッハ数・圧縮性流れ計算機

マッハ数と比熱比 γ を設定して等エントロピー関係・垂直衝撃波特性を即時計算。亜音速から極超音速まで対応し、流れ場を可視化します。

入力パラメータ

超音速
T/T₀
P/P₀
ρ/ρ₀
A/A*
q/P₀ (動圧比)
V/a₀

垂直衝撃波(M₁ → M₂)

M₂ (衝撃波後)
P₂/P₁
T₂/T₁
ρ₂/ρ₁
等エントロピー
$\frac{T}{T_0}=\left(1+\frac{\gamma-1}{2}M^2\right)^{-1}$
$\frac{P}{P_0}=\left(\frac{T}{T_0}\right)^{\gamma/(\gamma-1)}$
垂直衝撃波
$M_2^2=\frac{M_1^2+\frac{2}{\gamma-1}}{\frac{2\gamma}{\gamma-1}M_1^2-1}$
流れ場の可視化
等エントロピー特性曲線(M 0.1〜5)

マッハ数・圧縮性流れとは

🧑‍🎓
マッハ数って、飛行機の速さを表す「マッハ2」みたいなやつですよね?具体的に何を計算しているんですか?
🎓
ざっくり言うと、流れの速さと音速の比だね。式は $M = V / a$ で、$V$が流速、$a$がその場の音速だ。このツールの左のスライダーでマッハ数を0から5まで変えられるよ。例えばM=0.8は亜音速の旅客機、M=2.0は超音速戦闘機くらいのイメージだ。
🧑‍🎓
「圧縮性」って何ですか?マッハ数が変わると、空気の性質が変わるってこと?
🎓
その通り!空気は速く流れると密度が変わっちゃうんだ。車の窓から手を出すくらいの低速なら密度はほぼ一定だけど、マッハ0.3を超えるとこの「圧縮性」の影響を無視できなくなる。このツールでマッハ数を上げていくと、温度比$T/T_0$や圧力比$P/P_0$がどう変わるか、グラフで一目瞭然だよ。
🧑‍🎓
「垂直衝撃波」って、超音速機についてるあの三角の雲みたいなやつですか?ツールの右のスライダー「衝撃波上流マッハ数」で何がわかるんですか?
🎓
鋭いね!あの雲は衝撃波の前後で空気が急冷縮して発生する現象だ。垂直衝撃波は、流れに対して真っ直ぐ立つ最もシンプルな衝撃波だよ。右のスライダーで上流のマッハ数$M_1$を変えると、衝撃波を通過した後の下流のマッハ数$M_2$がどう急激に低下するか、圧力や温度がどうジャンプするかが計算できる。ロケットノズルの出口や超音速インテークの設計で必須の知識だ。

物理モデルと主要な数式

等エントロピー流れの関係式
摩擦や熱の出入りがない断熱的な流れ(等エントロピー流れ)では、流れが止まった仮想的な「よどみ点」の状態($T_0, P_0$)と、実際に流れている場所の静状態($T, P$)の間に以下の関係が成り立ちます。これが航空機やエンジンの基本的な性能計算の根幹です。

$$\frac{T}{T_0}=\left(1+\frac{\gamma-1}{2}M^2\right)^{-1}$$

$$\frac{P}{P_0}=\left(\frac{T}{T_0}\right)^{\gamma/(\gamma-1)}$$ $T$: 静温度, $T_0$: よどみ温度, $P$: 静圧, $P_0$: よどみ圧, $M$: マッハ数, $\gamma$: 比熱比(空気では約1.4)

垂直衝撃波の関係式
超音速流れ中に発生する垂直衝撃波を通過すると、流れは不連続に変化し、超音速から亜音速になります。衝撃波の前後で質量・運動量・エネルギーの保存則を適用すると、以下の関係が導かれます。

$$M_2^2=\frac{M_1^2+\frac{2}{\gamma-1}}{\frac{2\gamma}{\gamma-1}M_1^2-1}$$

$M_1$: 衝撃波上流のマッハ数(>1), $M_2$: 衝撃波下流のマッハ数(<1)
この式から、上流マッハ数$M_1$が大きいほど、下流マッハ数$M_2$は1/√2(約0.7)に近づき、圧力と温度の上昇が極めて大きくなることがわかります。

実世界での応用

航空機・宇宙機の設計:機体表面の圧力分布、翼まわりの流れ、エンジンインテークからコンプレッサー入口への空気の流入状態を計算するために必須です。特に遷音速(M=0.8〜1.2)での挙動は設計の難所です。

ターボ機械・エンジン:ジェットエンジンの圧縮機やタービン、ロケットエンジンのノズル内部は超音速流れが発生する領域です。等エントロピー関係を用いて効率的なエネルギー変換を設計します。

風洞実験データの解析:風洞実験で計測されたモデル表面の静圧データを、等エントロピー関係式を用いてよどみ点条件に関連付け、実際の飛行条件との対応を取ります。

CAE(数値流体力学:CFD)解析:超音速・極超音速流れの数値シミュレーションを行う際、衝撃波を正しく捉えることは最大の課題の一つです。本ツールで学ぶ垂直衝撃波の理論は、CFD結果の検証や境界条件設定の基礎知識として活用されます。

よくある誤解と注意点

このツールを使い始める際、特にCAE初心者がハマりがちなポイントがいくつかあります。まず「よどみ状態は止まっている場所の値ではない」という点。$T_0$や$P_0$は、流れを「断熱的かつ等エントロピーで減速させて止めたらどうなるか」という仮想的な基準値です。例えば、飛行機の機首のよどみ点では実際に流速はゼロですが、エンジンインテーク内部など流れている場所で使う「よどみ圧」は、センサーで測る全圧のことで、流れが止まっているわけではありません。

次に、比熱比 $\gamma$ の扱い。ツールでは空気の1.4で固定していますが、実務ではこれが大きな落とし穴に。燃焼ガスが関わるロケットノズルやジェットエンジンのタービン後流では、$\gamma$ は1.3や1.2程度まで低下します。この値を間違えると、温度比や圧力比の計算が大幅にずれ、性能予測を誤ります。常に「どんな流体か?」を第一に確認しましょう。

最後に、「垂直衝撃波は現実ではレア」という認識。ツールで学ぶ垂直衝撃波は最もシンプルなモデルです。実際の超音速機周りで発生するのは、斜め衝撃波や膨張波の複雑な組み合わせ。垂直衝撃波に近い現象は、超音速流れが直壁に正面衝突するような極端なケースで主に見られます。基礎を学ぶには最適ですが、その限界も理解しておくことが実務への第一歩です。

関連する工学分野

このツールの計算ロジックは、航空宇宙以外の様々な先端分野で応用されています。例えば「内燃機関」。自動車エンジンの過給機(ターボチャージャー)では、コンプレッサー出口からインテークマニホールドへの流れは亜音速〜遷音速。ここでの圧力回復の設計には、等エントロピー関係が基礎として使われています。特に、スロットル弁での絞りは一種のノズル流れとみなせます。

もう一つは「化学プラントの安全設計」。高圧ガスが配管の破裂部から大気中に噴出する際、出口でマッハ1になるチャークフローが発生します。この時の質量流量を正確に見積もるには、$A/A^*$(臨界断面積比)の計算が不可欠。ガス漏洩量を求め、避難区域を設定する安全工学(セーフティエンジニアリング)の根幹です。

さらに「MEMS(微小電気機械システム)」の分野でも応用が。例えば、微細なガスチャネルを使った分析デバイスでは、流路寸法と分子平均自由行程の関係から、連続体の圧縮性流れ理論が適用できる限界を見極める必要があります。マッハ数は低くても、微小スケール特有の現象を理解するための基礎パラメータとして重要な役割を果たしています。

発展的な学習のために

このツールで基礎に慣れたら、次のステップとして「1次元流れの体系的理解」に進むことをお勧めします。具体的には、等エントロピー流れに「摩擦」の影響を加えたファンノー流れ(Fanno flow)と、「加熱・冷却」の影響を加えたレイリー流れ(Rayleigh flow)を学びましょう。これらは合わせて「1次元圧縮性流れの三大基礎」と呼ばれ、例えばジェットエンジンの燃焼器(加熱)や長い排気管(摩擦)の設計に直結します。

数学的背景としては、ツールで使われている関係式のほとんどが「保存則の代数からの導出」です。次の学びでは、質量・運動量・エネルギーの各保存式を微分形で書き、それを簡単な仮定(断熱、摩擦なしなど)で積分して、今のツールの式がどう出てくるかを自分で追ってみてください。このプロセスを経験すると、仮定を変えたとき(例えば比熱が温度変化するとき)にどこをどう修正すればいいか、応用力が格段に上がります。

最終的には、これらの1次元モデルを「流れ場の一部」として捉え直すのが目標です。CAEシミュレーション(CFD)では、これらはノズル入口・出口の境界条件を決めたり、衝撃波の位置を初期推定したりする「知識」として生かされます。ツールでパラメータ感覚を磨いたら、簡単な2次元ノズルのCFD解析を実行し、計算結果の中に等エントロピー関係や衝撃波条件がどこに現れているかを探検してみると、理論と実践が結びつく最高の体験になるでしょう。