ノズル流れ
ノズル流れの理論基礎
概要
先生、ラバールノズルって収束-発散形状で超音速を出すやつですよね？ なんであの形だと音速を超えられるんですか？
いい質問だ。圧縮性流体力学の基本定理で、断面積変化と速度変化の関係は
で表される。亜音速（M<1）では面積を絞ると加速し、超音速（M>1）では面積を広げると加速する。だからM=1を達成するスロート（最小断面積）を経由して超音速に遷移できるんだ。
スロートでちょうどM=1になるんですね。これがチョーク条件ですか？
その通り。チョーク（choking）とは、スロートでマッハ数が1に達して質量流量が最大になる状態だ。この臨界質量流量は
で決まる。背圧をいくら下げてもスロートを通過する質量流量はこれ以上増えない。
面積-マッハ数関係
ノズルの断面積からマッハ数を求めるにはどうするんですか？
等エントロピー流れの仮定のもと、面積比とマッハ数の関係は
となる。これは与えられた $A/A^*$ に対して亜音速解と超音速解の2つの解を持つ。どちらが実現するかは背圧条件で決まるんだ。
非線形方程式だから解析的には解けないですよね？
Newton-Raphson法で簡単に解ける。初期値を亜音速側（M<1）か超音速側（M>1）のどちらかに設定すれば、対応する解に収束する。等エントロピー関係も合わせて
で温度比・圧力比が求まる。空気（$\gamma=1.4$）でM=2なら $T_0/T=1.8$, $p_0/p=7.824$ だ。
ノズル内の衝撃波
背圧が設計値と違うとどうなるんですか？
背圧が高すぎると、発散部内に垂直衝撃波が立つ。衝撃波前後の圧力比はRankine-Hugoniot関係
で決まる。衝撃波後は亜音速になるから、ノズル出口まで減速する。背圧をさらに下げると衝撃波がノズル出口に移動し、最終的にはノズル外の膨張波-圧縮波パターンになる。
ロケットエンジンの排気で見えるダイヤモンドパターンがそれですか？
そう。過膨張（出口圧<背圧）では斜め衝撃波、不足膨張（出口圧>背圧）ではPrandtl-Meyer膨張波が形成される。あのマッハディスクの位置は背圧比で決まるんだ。
ド・ラバルノズルの発明——蒸気タービンを救った「くびれ」の発見
収縮拡大ノズル（de Laval nozzle）を考案したのは、スウェーデンの技術者Gustav de Lavalです。1880年代、蒸気タービンの効率改善に取り組んでいた彼は、蒸気を超音速まで加速するには喉部で一度音速に達した後に拡大させる必要があることを実験的に発見しました。当時は気体の圧縮性理論が十分に整備されていなかったにもかかわらず、試行錯誤でその形状にたどり着いたのです。現代のロケットエンジン・ジェットエンジンのノズルはすべてこの原理で動いています。理論より先に実験で「答え」を見つけた好例といえます。
ノズル流れの数値計算手法
準1次元ノズルの数値解法
準1次元ノズル流れって、CFDの入門として最適なんですよね？
そう。Euler方程式の1次元版を面積変化項付きで解く問題で、圧縮性CFDの基本スキームを全て試せる。保存形で書くと
となる。右辺の $p \, dA/dx$ がソース項として効く。
これをどのスキームで解くのが一般的ですか？
教育的にはMacCormack法（予測子-修正子の2段階陽解法）が古典だ。Anderson教科書の定番問題だね。実用ではRoe法やAUSMでフラックスを計算し、MUSCL再構築で2次精度にする。
2D/3Dノズルのメッシュ戦略
2次元や3次元のノズル解析ではメッシュをどう作りますか？
軸対称ノズルなら2D軸対称メッシュで十分だ。構造格子が望ましく、スロート部は格子を密にする。
| 領域 | 格子方針 | 理由 |
|---|---|---|
| 収束部 | 壁面法線方向に20層以上 | 境界層の解像 |
| スロート近傍 | 軸方向に高密度 | 音速条件の正確な捕獲 |
| 発散部 | 衝撃波位置に適応的に密 | 衝撃波の鮮明な捕獲 |
| 中心軸 | 特異性回避（wedgeセル等） | 軸対称計算の数値安定性 |
3Dだと推力偏向ノズル（TVC）とかもありますよね？
推力偏向ノズルや多ノズルクラスターでは3D計算が必須だ。ノズル壁面の曲率が大きい部分ではプリズム層（inflation layer）を入れて境界層を解像し、コア部はポリヘドラルか六面体で埋める。
境界条件の設定
ノズル計算の境界条件はどう設定するんですか？
これがノズルCFDの核心部分だ。
- 入口: 全温 $T_0$ と全圧 $p_0$ を固定（燃焼室条件）。流れ方向は軸方向
- 出口: 超音速流出の場合は全変数を外挿（情報が上流に伝わらない）。亜音速域がある場合は背圧 $p_b$ を指定
- 壁面: 断熱壁（adiabatic）or 等温壁。No-slip条件
- 対称軸: 軸対称条件 or 対称面条件
Fluentでは「Pressure Inlet」で $p_0, T_0$ を、「Pressure Outlet」で背圧を設定する。超音速出口なら「Pressure Outlet」の背圧値は実質的に無視される。
亜音速-超音速遷移の数値的扱い
スロートでのM=1遷移は数値的に問題にならないですか？
良い着眼点だ。Euler方程式の双曲型性質がM=1で変化するため、数値スキームによっては遷移点近傍で問題が生じうる。Roe法のエントロピー修正（entropy fix）やHarten-Hymanの修正が必要になることがある。AUSM系は元々この問題に強い設計になっている。
圧力ベースソルバーでは超音速は解けないんですか？
最新のPressure-Based Coupled Solver（例：Fluentのcoupled scheme）は全速度域に対応しているが、衝撃波を含むノズル流れでは密度ベースソルバーの方が安定で高精度だ。特にノズル内衝撃波の位置と強度を正確に予測するなら密度ベース一択だね。
ノズル数値解法で「密度基準か圧力基準か」で悩む理由
ラバルノズルの超音速領域を数値計算するとき、初心者がよく迷うのが「密度基準ソルバーと圧力基準ソルバーのどちらを使うべきか」だ。亜音速から超音速まで連続して解くノズル流れでは、マッハ数が1を跨ぐため圧力基準ソルバーのPresto!や隣接補正が上手く機能しないことがある。一般的には「超音速・圧縮性が支配的」なら密度基準が安定しやすく、特に時間前進（density-based explicit）を使うと衝撃波の捕捉が鮮明になる。ただし壁面近傍のY+管理が緩いと境界層が暴れるため、メッシュと解法のセットで考えるのが正解だ。
ノズル流れの実務適用
ロケットノズルの設計解析
ロケットノズルの設計CFDはどんな流れで進めるんですか？
実務的なフローを見てみよう。
1. 等エントロピー設計: MOC（特性曲線法）で最小長さノズル or ベル型ノズルの壁面輪郭を決定
2. 冷態流CFD: まず理想気体で流れ場を確認。衝撃波やスロート付近の流れパターンを検証
3. 高温気体CFD: 燃焼ガスの実在気体効果（解離・再結合）を考慮した計算
4. 壁面冷却: 再生冷却チャネルやフィルム冷却の熱伝達解析
5. 構造連成: 熱荷重と圧力荷重による構造応力解析
MOC（特性曲線法）でノズル形状を決めるんですね。
Mach数の設計点（例えばM=3.0）とスロート半径を指定すれば、特性曲線の交点を追跡して一様平行流を出口で実現するノズル壁面形状が得られる。Anderson「Modern Compressible Flow」に詳しいアルゴリズムが載っている。
背圧変動のパラメトリック解析
背圧を変えたときのノズル内流れはどう変わりますか？
背圧比 $p_b/p_0$ によって以下のレジームが現れる。
| $p_b/p_0$ | 状態 | 流れの特徴 |
|---|---|---|
| $= 1$ | 流れなし | 等圧 |
| $> p^*/p_0$ | 全域亜音速 | ベンチュリ効果のみ |
| $= p^*/p_0$ | スロートでM=1 | チョーク開始 |
| 中間値 | ノズル内垂直衝撃波 | 衝撃波位置は背圧に依存 |
| 設計値 | 完全膨張 | 出口で設計マッハ数 |
| 設計値より低い | 過膨張/不足膨張 | ノズル外の波パターン |
ここで $p^*/p_0 = (2/(\gamma+1))^{\gamma/(\gamma-1)}$ で、空気では約0.528だ。
CFDでパラメトリックに背圧を振ると、衝撃波がノズル内を移動する様子が見えるんですね。
そう。Fluentなら背圧をパラメータにしたDesign Pointsを作成し、バッチ計算で一連の背圧条件を走らせる。結果をParaViewで動画にすると教育的にも素晴らしい可視化ができる。
3次元効果と非対称性
実際のノズルは軸対称じゃない場合もありますよね？
過膨張ノズルでは壁面から衝撃波が剥離する「フロー・セパレーション」が起きる。これが非対称になると横力（side load）が発生し、ノズルの構造的な問題になる。SSME（スペースシャトルメインエンジン）の設計でも大きな課題だった。
非対称な剥離を計算するには3Dが必要ですか？
検証のポイント
ノズルCFDの結果検証はどうするんですか？
以下を確認するのが基本だ。
- 質量流量: 理論値 $\dot{m}_{max}$ と1%以内で一致するか
- スロート圧力: $p^*/p_0 \approx 0.528$（空気）を確認
- 出口マッハ数: 面積比から計算される理論値と一致するか
- 壁面圧力分布: 実験データまたは等エントロピー解との比較
- 推力: $F = \dot{m} V_e + (p_e - p_a) A_e$ の推力方程式で確認
ノズル形状を「手計算」で最適化していた時代
現代ではノズル輪郭設計にCFDと最適化ループを組み合わせるのが当たり前だが、1950〜60年代の設計者は「特性曲線法（Method of Characteristics）」を手計算で解いてノズル形状を決めていた。流れ場の特性線（マッハ線）を一本一本グラフ用紙に引いて、等エントロピー膨張を実現する壁面輪郭を求める作業だ。NASA初期のロケットエンジンノズルはこうして設計された。今も教科書に載っているこの手法を理解すると、CFDの計算結果の「なぜこの形が最適か」が直感的に見えてくる。
ノズル流れのソフトウェア比較
ノズル流れ解析ツール
ノズル流れの解析にはどのソフトが適していますか？
ノズル流れは圧縮性CFDの典型問題だから、主要な商用・OSSツールは全て対応している。ただし機能の充実度には差がある。
| ツール | 密度ベース | 実在気体 | 化学反応 | MOC設計 | 向いている用途 |
|---|---|---|---|---|---|
| Ansys Fluent | Density-Based | NIST RGP | Species Transport | なし | 汎用ノズルCFD |
| Ansys CFX | Coupled | 実在気体EOS | 反応流 | なし | ターボ機械ノズル |
| STAR-CCM+ | Coupled Flow | Tables/EOS | 詳細化学 | なし | 燃焼ノズル |
| OpenFOAM (rhoCentralFoam) | 密度ベース | 限定的 | reactingFoam | なし | 研究・教育 |
| NPSS (NASA) | 0D/1D | 各種 | CEA連携 | あり | エンジンサイクル解析 |
| NASA CEA | 0D | 多種 | 化学平衡 | 出力のみ | 燃焼ガス組成計算 |
NASA CEAって何ですか？
Chemical Equilibrium with Applications。燃焼室条件（燃料・酸化剤の種類と混合比）を入力すると、化学平衡状態のガス組成と熱力学特性を計算してくれる。ロケットノズルの入口条件を決める際の必須ツールで、Webインターフェースが無料公開されている。
Ansys Fluentでのノズル設定
Fluentでノズル計算する際の具体的な設定を教えてください。
基本設定はこうなる。
- ソルバー: Density-Based, Steady
- 空間離散化: 2nd Order Upwind（初期収束が難しければ1st Orderでスタート）
- フラックスタイプ: Roe-FDS（衝撃波をシャープに捕獲）
- 気体モデル: Ideal Gas（冷態流）/ Real Gas Property Table（高温ガス）
- 乱流: SST k-omegaモデル（壁面近傍の精度重視）
- 入口BC: Pressure-Inlet（$p_0$, $T_0$ 指定）
- 出口BC: Pressure-Outlet（背圧指定。超音速出口ではSupersonic/Initial Gauge Pressureを設定）
超音速出口の場合、背圧の値は意味がないんですよね？
そう。全出口面が超音速なら全変数が内部から外挿されるので背圧は無視される。ただし起動過程の過渡計算や部分的に亜音速になる場合は背圧が効くから注意が必要だ。
OpenFOAMでのノズル計算
OpenFOAMだとどのソルバーを使いますか？
rhoCentralFoamが最適だ。KNP（Kurganov-Noelle-Petrova）スキームベースの密度ベースソルバーで、衝撃波を含むノズル流れを安定に計算できる。設定例を示そう。
thermophysicalProperties:hePsiThermo+perfectGasfvSchemes:localEulerddt（定常）、vanLeerdiv（2次精度TVD）fvSolution: smoothSolver + GaussSeidel0/: p, T, U の初期条件と境界条件
BoundaryConditions:
- inlet:
totalPressure+totalTemperature - outlet:
waveTransmissive（非反射境界） - walls:
zeroGradient（断熱壁）
waveTransmissive境界って何ですか？
特性線に基づく非反射境界条件で、波が計算領域から出ていくのを妨げない。ノズル出口に人工的な反射波を作らないために重要だ。
ノズル設計で各ツールが最初に検証するベンチマーク
ノズル流れの計算ツールを評価するとき、CFD業界が長年使い続けているベンチマークが「1D等エントロピーノズル流れ」です。解析解が存在するため、計算値と厳密解を直接比較できます。圧力比を変えてノズル内に垂直衝撃波（正常衝撃波）が形成される条件でも解析解があり、衝撃波の位置と強さをkPa単位で検証できます。FLUENTでもOpenFOAMでも、商用コードの開発チームは新しいスキームを追加するたびにこのノズルケースで回帰テストを実施します。ツール導入時にまずノズル流れで検証する習慣は、エンジニアとして鉄板の作法です。
ノズル流れの先端研究
デュアルベルノズルとアダプティブノズル
ノズル設計の最新トレンドは何ですか？
大気圧が変化する上昇中にノズル膨張比を最適化する「高度補償ノズル」が注目されている。代表的な概念は3つだ。
- デュアルベルノズル: 内側ベルと外側ベルの2段構成。低高度では内側ベルの端で流れが剥離し、高高度では外側ベルまで完全に膨張する
- プラグノズル（エアロスパイクノズル）: 中心プラグの周りに膨張。外側は大気圧で自動的に調整される
- 伸展-偏向ノズル（E-Dノズル）: 壁面の一部を省略して高度補償効果を得る
SpaceXのRaptorエンジンもノズル最適化してますよね？
二相流ノズル
固体ロケットのノズルには固体粒子も混じってますよね？
そう。固体推進薬の燃焼ガスにはAl₂O₃粒子（直径1-10μm）が含まれ、比推力の低下やノズル侵食の原因になる。この二相流ノズルの計算にはEulerian-Lagrangianアプローチが使われる。
- 連続相（ガス）: Euler方程式
- 分散相（粒子）: Lagrangian追跡。粒子-ガス間の抗力・熱伝達のカップリング
FluentのDPM（Discrete Phase Model）やSTAR-CCM+のLagrangian Multiphaseで対応可能だ。粒子がスロートを通過するときの速度遅れ（velocity lag）が推力損失に寄与する。
推力偏向制御（TVC）のCFD
推力偏向ってどうやってCFDで予測するんですか？
推力偏向の方式は複数あり、CFDの扱いも異なる。
| TVC方式 | CFD手法 | 課題 |
|---|---|---|
| ジンバリング | メッシュ変形/オーバーセット | 可動部のメッシュ |
| 二次流体噴射（SITVC） | 横噴射のジェットインタラクション | 衝撃波-境界層干渉 |
| フレキシブルノズル | FSI連成 | 構造変形のカップリング |
| ピントルノズル | 可動ピンドルのメッシュ変形 | スロート面積変化 |
二次流体噴射は最も活発に研究されている。主流の超音速膨張流に横方向から亜音速or超音速ジェットを噴射し、衝撃波を誘起して横力を発生させる。この衝撃波-境界層干渉を正確に予測するにはDES以上の乱流モデルが望ましい。
形状最適化
CFDとノズル形状最適化を組み合わせた研究はありますか？
アジョイント法による感度解析が強力だ。推力やIsp（比推力）を目的関数として、ノズル壁面形状パラメータ（ベジエ曲線の制御点等）の勾配を効率的に計算できる。
- Ansys Fluent: adjoint solver搭載。壁面形状のメッシュモーフィングと連携
- SU2: オープンソースCFD。アジョイント設計最適化に強い。Stanford大学で開発
- DAFoam: OpenFOAMベースの離散アジョイント最適化フレームワーク
SU2は航空宇宙分野で人気ですよね。
そう。SU2はEuler/RANSのアジョイント法が洗練されていて、翼型やノズルの形状最適化で多数の論文がある。無料で使えるから大学の研究にも最適だ。
「エアロスパイクノズル」——次世代ロケットが再注目する理由
通常のベルノズルはある特定の高度（設計高度）でしか最適効率にならない。大気圧が変わると過膨張や不足膨張が生じて推力が落ちる。これを解決するのがエアロスパイクノズルで、スパイク状のプラグ周りに燃焼ガスが流れる構造だ。理論上は全高度で自動的に最適膨張に近づく「高度補償機能」を持つ。NASAは1970年代に研究を進めたが、冷却設計の難しさから実用化が遅れた。現在は3Dプリンティングで複雑な冷却流路が製造できるようになり、民間宇宙企業が改めて開発を進めている。CFDでの形状最適化が再び脚光を浴びている分野だ。
ノズル流れのトラブル対応
よくあるトラブルと対処法
ノズルのCFDで初心者がハマりやすいポイントを教えてください。
ノズルCFD特有のトラブルパターンを見ていこう。
1. スロートでチョークしない
症状: スロートのマッハ数が1に達しない
原因と対策:
- 背圧が高すぎる → 背圧を下げて $p_b/p_0 < 0.528$（空気）にする
- 入口にPressure-Inletではなくmass-flow-inletを使っている → チョーク以上の質量流量を指定するとソルバーが全圧を勝手に上げてしまう。全温全圧で入口を設定すべき
- メッシュが粗い → スロート近傍の軸方向メッシュを細かくする
2. 衝撃波が振動して収束しない
残差が下がらなくて衝撃波が行ったり来たりするんですが…
原因: ノズル内衝撃波は背圧に敏感で、数値的な振動が衝撃波位置を揺さぶることがある。
対策:
- CFL数を下げる（Fluentなら0.5-1.0程度）
- 1次精度で計算して定常解を得てから、2次精度に切り替える
- Fluentの場合、FMG（Full Multigrid）初期化を使う
- 残差が完全に定常にならない場合はURANS計算に切り替え（物理的に非定常な可能性）
3. 質量流量が理論値と合わない
理論式で計算した質量流量とCFDの値が数%ずれるのは正常ですか？
2-3%程度のずれは壁面境界層による有効断面積の減少（displacement thickness effect）で説明できる。ただし5%以上ずれるなら以下をチェック。
- スロート断面積の定義が正しいか（CAD形状のスロート径を確認）
- 気体定数 $R$ と比熱比 $\gamma$ が入力ミスしていないか
- 2D軸対称計算でwedge角度やaxis条件が正しいか
- 入口の全温全圧が正しく設定されているか
4. 壁面圧力にギザギザが出る
壁面圧力分布をプロットするとノコギリ状になってしまいます。
原因: メッシュ品質の問題。特に壁面法線方向に非均等なセル配列がある場合。
対策:
- 壁面に沿った方向のセルサイズ変化率（growth rate）を1.2以下にする
- スロートの曲率が大きい箇所で格子を十分に密にする
- 非構造格子の場合、壁面にプリズム層を入れる
- 結果の後処理でlocal averagingをかける（物理を変えずに可視化を改善）
5. 過膨張ノズルの壁面剥離
過膨張条件でノズル壁面からの剥離が予測できないんですが。
RANSでは過膨張ノズルの壁面剥離を正確に予測するのは困難だ。以下の経験則で剥離位置を推定してCFDと比較するといい。
Schmucker-Summerfieldの相関式:
ここで $p_{sep}$ は剥離点の壁面圧力、$p_a$ は周囲大気圧。SST k-omegaモデルで壁面圧力がこの値まで回復する位置を剥離点と推定できる。
RANSで限界がある場合はどうしますか？
ノズル流れが「詰まった」と思ったら——過膨張と不足膨張の見分け方
ロケットノズルのCFDで「背圧を変えているのに質量流量が変わらない」という現象に遭遇した場合、多くは「チョーク（詰まり）」状態に入っている。ノズルスロート部でマッハ1に達すると、下流の情報が上流に伝わらなくなるため、背圧をどう変えても流量は変わらない——これは正常動作だ。問題は出口での過膨張と不足膨張の判断。出口圧力が背圧より低ければ過膨張で斜め衝撃波が外側で発生し、高ければ不足膨張で膨張扇が出る。CFDの可視化で出口付近の圧力分布を見るだけでどちらの状態かわかるので、まずそこを確認するのが鉄則だ。
関連トピック
なった
詳しく
報告