収縮拡大ノズル設計計算機 戻る
超音速流体力学

収縮拡大ノズル(ド・ラバルノズル)設計計算機

出口マッハ数・面積比・総圧を設定すると、等エントロピー流れの理論からマッハ数・圧力・温度の軸方向分布を即時計算。垂直衝撃波位置も可視化できます。

パラメータ設定
出口マッハ数 Me
スロート面積 A* [cm²]
cm²
総圧 P₀ [bar]
bar
総温 T₀ [K]
K
比熱比 γ
プリセット
出口性能
計算結果
出口速度 [m/s]
推力 [N]
比推力 Isp [s]
質量流量 [kg/s]
面積比 Ae/A*
Pe/P₀
ノズル
理論・主要公式

等エントロピー流れの支配方程式:

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

収縮拡大ノズル(ド・ラバルノズル)とは

🙋
ド・ラバルノズルって、ただの管がくびれてまた広がってるだけに見えますけど、何がすごいんですか?
🎓
大まかに言うと、この特殊な形が「亜音速の流れを超音速に変える魔法の管」なんです。例えばロケットエンジンのノズルはこれ。上のシミュレーターで「出口マッハ数」のスライダーを2.0にしてみて。くびれた部分(スロート)でマッハ数がちょうど1(音速)になって、その先でどんどん加速してるのがわかるよね。
🙋
え、そうなんですか!確かにくびれのところで1になってます。でも、なぜそこで音速になるんですか?
🎓
良いところに気づいたね。実は、流路が狭まると流れは加速する(マッハ数が増える)けど、音速(M=1)に達すると逆転現象が起きるんだ。そこから先は、流路を広げないとさらに加速できない。だから、超音速にするには「収縮してから拡大する」この形が必須なんだよ。試しに「比熱比γ」を1.67(アルゴンガスみたい)に変えてみると、同じマッハ数でも必要な断面積の変化が変わってくるのがわかる。
🙋
なるほど!でも、FAQに「衝撃波」って書いてあります。あの「衝撃波位置」のスライダーを動かすと、グラフが大きく変わりますけど、あれは何が起きてるんですか?
🎓
まさに現場で起こるトラブルを再現してるんだ。例えばロケットが設計より高い高度(低い外圧)で飛ぶと、ノズル内に「垂直衝撃波」が発生する。衝撃波を境に、一瞬で超音速から亜音速に戻り、圧力が跳ね上がる。シミュレーターで衝撃波をノズルの中ほどに移動させると、その前後でマッハ数と圧力が不連続に変化するのが見えるでしょ? これは推力損失の原因だから、設計では避けたい現象なんだ。

よくある質問

用途によります。設計したいノズルの出口速度が決まっている場合は出口マッハ数を先に入力し、逆にノズル形状(面積比)が先に決まっている場合は面積比を入力してください。どちらか一方を入力すれば、等エントロピー関係式からもう一方が自動計算されます。
ノズル背圧(出口外の圧力)が設計条件より高い場合、ノズル内部に垂直衝撃波が発生します。本ツールでは、背圧条件に応じて衝撃波前後のマッハ数・圧力比から位置を計算し、グラフ上に赤い破線で表示します。衝撃波通過後の圧力回復も確認できます。
空気の場合は1.4、燃焼ガス(ロケットエンジン)では約1.2〜1.3、ヘリウムでは1.67が標準です。気体の種類が不明な場合は、一般的な二原子分子ガスとして1.4を推奨します。比熱比が変わるとマッハ数と面積比の関係が変化するため、正確な値を入力してください。
ノズル内部で流れが剥離したり、垂直衝撃波が発生した場合、等エントロピー流れの仮定が成立しなくなるためです。本ツールでは衝撃波位置以降の等エントロピー計算を停止します。背圧を下げるか、ノズル面積比を大きくすることで、より広範囲の計算が可能になります。

実世界での応用

ロケット・ミサイルエンジン:燃焼室で生成された高温高圧ガスを超音速まで加速し、大きな推力(反作用力)を得るために使用されます。宇宙空間のような低圧環境でも効率的に推力を発生できる形状は、この理論に基づいて設計されています。

超音速・極超音速風洞:実験部で超音速流れを発生させるためのノズルとして使われます。目標マッハ数に応じて精密に設計されたド・ラバルノズルが、試験モデルに均一な超音速流を供給します。

蒸気タービン:特に発電所用の大型タービンの最終段などで、蒸気を超音速まで膨張させて効率よく仕事を取り出すために採用されることがあります。比熱比γは空気とは異なる値で計算します。

産業用スプレー・サンドブラスト:高速の流体ジェットを生成するノズルとして応用されます。流体の加速により微粒化や洗浄・加工の効率を高めます。衝撃波の発生はジェットの乱れの原因となるため、設計が重要です。

よくある誤解と注意点

このツールを使い始めるとき、いくつかつまずきやすいポイントがあるから気をつけてね。まず「入口条件は無限大のタンクから来る」という仮定を理解しておこう。シミュレータは入口で流速がほぼゼロ(マッハ数≒0)の状態から計算している。でも実機の燃焼室は有限の大きさで、そこにも流れがあるんだ。入口マッハ数が0.1以上あると、このツールの単純な面積比の式だけでは正確な形状が決まらないから注意が必要だよ。

次に、「比熱比γ」の設定ミス。空気(γ=1.4)でロケットノズルを設計してしまうのはよくある勘違い。ロケットの燃焼ガスは水蒸気や二酸化炭素が多く、γは1.2前後になることもある。例えばγを1.4から1.25に変えるだけで、同じ出口マッハ数5.0を達成するのに必要な面積比は約40も大きくなるんだ。使う気体を絶対に確認しよう。

最後は「理想と現実のギャップ」。この計算は壁面摩擦や熱損失、二次元流れを一切無視した「一次元等エントロピー流れ」が前提だ。実際のノズル、特に小型のものでは境界層の影響で有効な流路面積が狭まり、設計通りの性能が出ないことがある。シミュレーションで「最適設計」できても、そこからが本当の設計の始まりだと覚えておいて。

使い方ガイド

  1. 出口マッハ数(Me)と喉部面積A*を入力します。ロケットエンジンの場合、典型的なMeは3~4.5、液体酸素・ケロシン燃焼ガスではγ=1.2程度を想定
  2. 全温T₀(燃焼室温度、通常3000~3500K)と全圧P₀(燃焼室圧力、1~20 MPa)を設定
  3. 計算ボタンを押すと、等エントロピー関係式により喉部から出口までの軸方向のマッハ数・圧力・温度分布、衝撃波位置、推力と比推力Ispが自動計算される

具体的な計算例

液体ロケットエンジン(LOX/ケロシン)の設計例:Me=4.0、T₀=3200K、P₀=15 MPa、A*=0.005 m²、作業流体γ=1.2、分子量28 g/molの場合、出口マッハ数4.0での圧力比Pe/P₀≈0.027、出口温度Te≈1200K、出口速度Ve≈4100 m/s、質量流量約45 kg/sとなります。喉部からの面積比Ae/A*≈4.8で、衝撃波が拡大部に生じる条件では推力が500 kN、比推力Isp≈360秒程度に低下する可能性があります。

実務での注意点

  1. 燃焼ガスの比熱比γは燃料組成により変動します。LOX/ケロシンはγ=1.2~1.25、固体推進薬(アルミニウム/アンモニウムパークロレート)ではγ=1.15程度となり、Isp計算に±10秒の誤差が生じます
  2. 喉部での流れは必ず音速(M=1)であることが前提です。過膨張(Pe < 環境圧力)になる場合、拡大部で衝撃波が発生し、推力喪失が生じるため出口面積を再設計が必要です
  3. 高温ガスの実際の粘度・熱伝導率変化は無視されるため、15 MPa以上の高圧条件では燃焼室損失補正係数0.95~0.98を乗じて実現推力を見積もります