圧縮機CFD解析
圧縮機CFDの理論基礎
概要
圧縮機のCFD解析って、軸流と遠心でアプローチが違うんですか？
基本の支配方程式は同じだが、遠心圧縮機ではディフューザの役割が大きく、軸流ではブレードローディングの管理が主題になる。ただし共通するのは圧力比と断熱効率の予測精度がCFDに求められるということだ。
圧力比と断熱効率
圧力比ってどう定義されるんですか？
全圧比として定義される。
$p_0$ は全圧（よどみ圧力）。そして断熱効率は等エントロピー過程と実過程の仕事比だ。
温度から計算するんですね。CFDでは直接揚程や全圧から出せますか？
質量流量加重平均の全圧・全温を入口と出口で取得して計算する。CFX-PostやParaViewのmassFlowAve関数を使うのが標準的だ。
圧縮性効果
遠心圧縮機のインペラ先端速度って音速に近いですよね？
そう。ターボチャージャー用遠心圧縮機ではインペラ先端周速が400～500 m/sに達し、相対マッハ数は1.2を超えることもある。だから圧縮性を無視できない。
超音速流れが翼列内に出るわけですか。
インレット付近で超音速、翼間流路で衝撃波を介して減速する。衝撃波-境界層干渉による損失増大はCFDの精度を左右する重要な物理だ。
使用ソフトウェア
遠心圧縮機に強いソフトって何ですか？
Ansys CFX + TurboGridが産業界で最も実績がある。遠心インペラの子午面形状からTurboGridで構造格子を自動生成できる。NUMECA FINE/TurboのAutoGrid5も遠心に強く、スプリッタ翼のメッシュ生成が優秀だ。STAR-CCM+はポリヘドラルメッシュ+自動プリズムで手軽に始められるが、翼間流路の格子品質はTurboGridに及ばないことが多い。
圧縮機効率1ポイントの経済的インパクト
大型LNG液化プラントに使われる遠心圧縮機では、断熱効率が1ポイント上がるだけで年間数億円の電力コスト削減に直結します。圧力比が高いほど、効率改善の効果は指数関数的に大きくなる——これが圧縮機CFDに多大な開発費を投じる理由です。現代の産業用遠心圧縮機の断熱効率は85～90%に達していますが、その「最後の数%」を削り出すためにLES計算や最適化アルゴリズムが使われ続けています。
圧縮機CFDの数値計算手法
サージ予測のアプローチ
CFDでサージラインを予測するにはどうすればいいですか？
定常計算で出口背圧を徐々に上げていき、収束しなくなる点を近似的なサージ限界とする方法が最も実用的だ。ただし物理的なサージは系統全体の動的不安定だから、正確には非定常Full-Annulus計算が必要になる。
Full-Annulusって全周を計算するんですか？ それは大変ですね。
1ピッチの周期計算でサージは捉えられない。回転失速セルは周方向に伝播するから、全周（360度）を非定常で計算する必要がある。セル数は1ピッチの翼枚数倍だから、20枚翼なら計算量も20倍だ。
Harmonic Balance法
もう少し軽い方法はないんですか？
Harmonic Balance（調和バランス）法やNon-Linear Harmonic法がある。非定常変動を周波数領域で捉えるから、時間方向の計算コストを大幅に削減できる。CFXではTime Transformation法、FINE/TurboではNonlinear Harmonic法として実装されている。
どのくらいコスト削減できますか？
時間積分法の1/5～1/20で済むことが多い。ただし複数の周波数成分が干渉する強い非定常にはまだ限界がある。
遠心圧縮機のサージ
遠心圧縮機のサージは軸流と違いますか？
遠心圧縮機ではディフューザでの失速がサージの引き金になることが多い。特にベーン付きディフューザ（VD）の場合、VDのインシデンス角が大きくなると一気に失速する。ベーンレスディフューザ（VLD）のほうがサージマージンは広いが効率は劣る。
| ディフューザ型 | サージマージン | ピーク効率 | 適用 |
|---|---|---|---|
| ベーンレス（VLD） | 広い | やや低い | 自動車ターボ、可変運転 |
| ベーン付き（VD） | 狭い | 高い | 産業用、航空エンジン |
| パイプディフューザ | 中程度 | 高い | 高圧力比用途 |
自動車のターボチャージャーがベーンレスなのは運転範囲の広さが理由なんですね。
そう。エンジン回転数の広い範囲で使うから、サージマージンの確保が最優先だ。
ベーンレスディフューザが自動車ターボに多い理由
自動車用ターボチャージャーの遠心圧縮機にはベーンレスディフューザが多く採用されています。理由はシンプルで「運転範囲が広い」から。バネードディフューザは設計点では高効率ですが、off-designでサージマージンが急激に狭まる。ガソリン車のエンジンはアイドルから6000rpmまで幅広く使われるため、流量範囲の広いベーンレスが選ばれます。CFDではこの設計トレードオフを圧縮機マップ全域で確認することがエンジン開発の標準プロセスになっています。
圧縮機CFDの実務適用
解析ワークフロー
遠心圧縮機の典型的な解析フローを教えてください。
以下の手順が標準だ。
1. 1D設計: Concepts NRECのCOMPALやAxSTREAMでMean-Line設計。圧力比、流量、回転数から基本寸法を決定
2. 子午面設計: BladeGenやAxSTREAMでハブ・シュラウド曲線と翼角分布を定義
3. 3D翼形状定義: BladeGenでスプリッタ翼を含むフル3D形状を出力
4. メッシュ生成: TurboGridでH/J/L型の構造格子生成
5. CFD: CFXで定常MRF解析（設計点）→ 背圧変化で特性曲線取得
6. 最適化: optiSLangやFINE/Designで翼角・子午面形状を自動探索
BladeGenって何ですか？
Ansysのターボ翼形状定義ツールだ。ハブ/シュラウドの子午面プロファイルと各スパン位置での翼角分布（beta分布）を入力すると3D翼面を生成する。TurboGridへの直接連携が特徴だよ。
遠心特有のメッシュ注意点
遠心インペラのメッシュで軸流と違う点は？
いくつかある。
- 子午面の曲率: ハブ・シュラウドの曲率が大きく、J/L型トポロジの歪みが出やすい
- スプリッタ翼: メインブレードとスプリッタでトポロジが異なるため、TurboGridでは別ブロックとして扱う
- ディフューザ接続: インペラ出口→ベーンレスディフューザ→ボリュートの接続部でメッシュ型が変わる
ボリュートはTurboGridで作れますか？
作れない。ボリュートは非回転の非軸対称形状だから、通常のCADメッシュ（Ansys MeshingやFluent Meshing）で別途生成して、GGI（General Grid Interface）で接続するのが一般的だ。
結果評価のポイント
遠心圧縮機のCFD結果で特に確認すべき項目は？
以下を重点的に見る。
| 評価項目 | 確認方法 | 設計基準 |
|---|---|---|
| インペラ出口のジェット/ウェイク構造 | スパン断面のMach数・全圧 | ジェット/ウェイク比が過大でないか |
| ディフューザ内の圧力回復 | 子午面の静圧分布 | $C_p = 0.5 \sim 0.7$ |
| チップ漏れ渦 | シュラウド面の流線 | 渦がメインブレード前縁に達していないか |
| スプリッタ翼の入口インシデンス | スプリッタ前縁の圧力分布 | 急な吸い込み面加速がないか |
スプリッタ翼の位置が効率を決める
遠心インペラの実践設計でよく問われるのが「スプリッタ翼をどこに置くか」という問題。メインブレード前縁から何%のチャネル長の位置に入れるかによって、低流量時のサージマージンと設計点効率がトレードオフになります。経験則では50～60%位置がバランスが良いとされますが、それはあくまで出発点。実務では±5%を変えながら3～5ケースのCFDを回して最適点を探すのが定番フローです。「経験＋CFDの繰り返し」が遠心インペラ設計の本質です。
圧縮機CFDのソフトウェア比較
ターボチャージャーの特殊性
自動車用ターボチャージャーの圧縮機CFDには特有の課題がありますか？
いくつかある。
- 広い運転範囲: エンジン回転数1000～6000rpmに対応するためマップ幅が広い
- 高周速: インペラ先端周速500m/s超、先端マッハ数1.3以上
- コンパクト設計: ボリュート形状の制約が厳しい
- 過渡応答: エンジン加速時のサージ回避
先端マッハ数1.3って相当速いですね。衝撃波が出ますか？
出る。前縁付近で弓形衝撃波が形成されて、それが隣接翼の吸い込み面に入射する。衝撃波-境界層干渉による損失がチョーク側の効率低下の主因だ。
コンプレッサーマップの予測精度
CFDでターボのコンプレッサーマップはどのくらい合いますか？
典型的な精度を示そう。
| 指標 | CFD vs 実験 | 備考 |
|---|---|---|
| 質量流量（設計点） | ±2% | 良好 |
| 圧力比（設計点） | ±1～3% | 良好 |
| 断熱効率 | ±1～3ポイント | メッシュ依存性大 |
| チョーク流量 | ±3% | 喉部面積に敏感 |
| サージライン | 定性的に一致 | 定常では限界あり |
効率1～3ポイントの誤差って、実用的にはどうですか？
設計の相対比較（A案 vs B案）には十分だ。絶対値の予測にはメッシュ感度スタディとモデル校正が必要になる。
最近のトレンド
ターボチャージャーCFDの最近の動向は？
3つの大きなトレンドがある。
1. 電動ターボ（eターボ）: モータ/ジェネレータ内蔵でエアフロー制御。圧縮機-タービン間の軸トルクバランスが変わり、過渡CFDの重要性が増している
2. LES/DESによる騒音予測: NVH要求の厳格化で、圧縮機の高周波騒音をDESやLESで予測する需要が増加
3. 3Dプリンティング翼型: 積層造形で従来加工できなかった翼形状を実現。CFDでの形状最適化の自由度が大幅に拡大
eターボがCFDワークフローを変えた
電動ターボチャージャー（eターボ）の登場は、圧縮機CFDの常識を変えつつあります。従来は「タービン仕事＝圧縮機仕事」の軸トルクバランスだけを考えれば良かったのに、eターボではモータ/ジェネレータのトルクが可変で介入します。過渡時の流量条件が従来ターボより格段に複雑になり、定常マップではなく過渡応答CFDの重要性が高まっています。GT-SuiteなどのシステムシミュレーションとCFDを組み合わせた1D-3D連成解析が自動車業界の新しいスタンダードになっています。
圧縮機CFDの先端研究
RANSの限界
SST k-omegaでは捉えられない現象って何ですか？
大規模な剥離を伴う流れ場だ。ストール近傍やサージ発生時の非定常渦構造は、RANS乱流モデルでは数値拡散で平滑化されてしまう。スパイクタイプの回転失速の起源を解明するにはLESやDESが必要だ。
DES/SDESの活用
DESって具体的にどういう手法ですか？
Detached Eddy Simulation。壁近傍はRANS（SSTなど）で解き、壁から離れた領域ではLESモードに切り替える。CFXではSAS（Scale-Adaptive Simulation）やSDES（Shielded DES）が利用可能だ。
SASとSDESの違いは？
SASはvon Karmanスケールに基づいて乱流粘性を動的に調整する手法で、LESに近い解像を得られるが完全なLESではない。SDESはより明確にRANS/LES領域を分離する。圧縮機の失速解析にはSDESが推奨される。
計算コストと実行戦略
LES/DESの計算コストはどのくらいですか？
概算を示す。
| 手法 | メッシュ規模 | タイムステップ | 必要回転数 | 総コスト比（RANS=1） |
|---|---|---|---|---|
| RANS定常 | 100万/ピッチ | - | - | 1 |
| URANS | 100万/ピッチ | 翼通過の1/20 | 5～10回転 | 50～100 |
| SAS/SDES | 500万/ピッチ | 翼通過の1/50 | 10～20回転 | 500～2000 |
| Wall-Resolved LES | 5000万/ピッチ | CFL<1 | 20～50回転 | 10000～ |
LESは現実的じゃないですね…
Full-AnnulusのWall-Resolved LESは研究レベルでも一部の翼列でしか実現していない。実務ではSDESやSASが現実的な上限だ。GPUベースのソルバー（Ansys Fluent GPU, CONVERGE等）で計算コストの壁が徐々に下がりつつある。
最近の研究動向
学術的にはどんな研究が進んでますか？
注目すべきトピックを挙げる。
- スパイク型失速の起源: 翼端クリアランスの逆流が起源という説と、前縁剥離が起源という説の論争。LESで決着をつける研究が進行中
- PINN（物理インフォームドNN）: RANSの解をベースにニューラルネットワークで非定常変動を高速予測する試み
- デジタルツイン: 実機センサーデータとCFDを連携させたリアルタイム性能監視
LESで初めて見えた「スパイク型失速」の誕生瞬間
圧縮機の回転失速には「モード型」と「スパイク型」の2種類があります。スパイク型は翼端クリアランス付近の局所的な逆流から突然発生し、RANSでは捉えにくい現象です。2000年代以降にLESが普及してから、スパイク型失速の誕生メカニズムが初めて可視化されました。翼端渦が前縁を超えて上流に噴き出す「スパイキング」の瞬間をLESコンターで見たとき、研究者の間で「やっとRANSが何を見落としていたか分かった」という話が広がったといいます。
圧縮機CFDのトラブル対応
数値的チョーク問題
圧縮機のCFDで「数値的チョーク」って何ですか？
メッシュの喉部面積が実際の幾何学的喉部面積と一致しないことで、チョーク流量がずれる現象だ。特にTurboGridのH型トポロジで前縁周りの角が甘いと、喉部が実効的に狭くなってチョーク流量が過小予測される。
どう対処すればいいですか？
前縁のO-gridを十分細かくし、翼面の曲率を正確に再現する。またJ型やL型トポロジに変更することで改善することもある。チョーク流量の計算値と幾何学的喉部面積からの理論値を比較して乖離を確認しよう。
背圧制御のテクニック
特性曲線を取得するとき、背圧の上げ方にコツはありますか？
いくつかのテクニックがある。
1. 前の収束解からリスタート: 背圧を一気に上げず、前の運転点の解を初期値に使う
2. 背圧のランピング: タイムステップで背圧を徐々に上昇させる（Expert Parameter: pressure ramp）
3. 質量流量制御に切り替え: サージ近傍では出口を質量流量指定にして安定させる（ただし流量-圧力の関係は見失う）
4. スロットルモデル: 出口に仮想スロットルバルブを設定して背圧を間接制御
スロットルモデルって具体的にどういうものですか？
出口に開口面積 $A_{th}$ の仮想オリフィスを置き、流量と背圧の関係を陽に与える。CFXではOpening BCとユーザー関数の組み合わせで実装できる。系統の容積効果を模擬できるのでサージの過渡現象に近い挙動が得られる。
収束判定の基準
圧縮機CFDの収束はどう判定すればいいですか？
残差だけでなく、物理量のモニターが重要だ。
| モニター量 | 収束判定 |
|---|---|
| RMS残差 | $10^{-5}$ 以下（定常）、各タイムステップで$10^{-4}$以下（非定常） |
| 質量流量入出口差 | 0.1%以内 |
| 段圧力比 | 変動幅0.1%以内で安定 |
| 段効率 | 変動幅0.1ポイント以内で安定 |
| トルク | 最終100イテレーションで変化率0.1%以内 |
残差が下がりきらなくても、物理量が安定していれば良いんですか？
ターボ機械のMixing Plane界面では残差がある程度の水準で飽和することがある。その場合は物理量モニターの安定性で判断するのが実用的だ。
「数値チョーク」と本物のチョーク——見分けられますか？
遠心圧縮機CFDで流量を増やしていくと、ある点で計算が急に収束しなくなる——これが「数値的チョーク」です。しかし実際には物理的なチョーク（翼間のMa=1超音速遷移）と、単なる背圧設定ミスによる数値不安定が混在していることが多い。見分け方の鉄則は「ディフューザ喉部のMa分布を確認する」こと。Ma=1の等値面がインペラ出口まで達していれば本物のチョーク。そうでなければ境界条件の再設定で解決できることがほとんどです。
関連トピック
なった
詳しく
報告