圧縮機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に及ばないことが多い。
ターボチャージャーの熱力学——なぜ排気で空気を圧縮できるのか
ターボチャージャーは「捨てる排気エネルギー」でコンプレッサーを駆動するという、エネルギー的に非常に巧みな仕組みだ。典型的な乗用車用ターボでは、コンプレッサー側が空気を圧力比2〜3倍に圧縮し、タービン側が排気の膨張エネルギーを機械仕事に変換してコンプレッサーを回す。CFD解析では理想的な断熱効率がコンプレッサー70〜80%、タービン70〜75%程度が設計目標とされる。この効率が1%下がると燃費が約0.5%悪化するため、ターボチャージャーメーカーは数値解析で翼型を0.1mm単位まで最適化している。
圧縮機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） | 狭い | 高い | 産業用、航空エンジン |
| パイプディフューザ | 中程度 | 高い | 高圧力比用途 |
自動車のターボチャージャーがベーンレスなのは運転範囲の広さが理由なんですね。
そう。エンジン回転数の広い範囲で使うから、サージマージンの確保が最優先だ。
高速回転体のメッシュ——ターボチャージャーCFDの空間的課題
ターボチャージャーのCFD解析で最も頭を悩ませるのがメッシュ生成だ。コンプレッサーインペラの後退翼（後ろ向きに曲がった翼）では、Hub側とShroud側で流路高さが大きく変わり、均質な品質のメッシュを作るのに専門的なスキルが必要になる。タービン側では高温排気（900℃超）による熱膨張で実際のクリアランスが設計値から変わるため、熱変形後の形状でメッシュを切り直す手順が必要なケースがある。業界の経験則では「クリアランスの部分に最低5〜8要素」が精度確保の目安とされており、全体メッシュ数が1000万要素を超えることも珍しくない。
圧縮機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$ |
| チップ漏れ渦 | シュラウド面の流線 | 渦がメインブレード前縁に達していないか |
| スプリッタ翼の入口インシデンス | スプリッタ前縁の圧力分布 | 急な吸い込み面加速がないか |
サージングマージンとCFD——設計者の命綱
ターボチャージャーのコンプレッサーがサージング（流れの周期的逆流）を起こすと、激しい振動と騒音が発生して最悪の場合インペラが破損する。CFDでサージングマージンを評価するには、コンプレッサーマップ（流量-圧力比特性）全体をスキャンし、サージラインを特定する必要がある。実務では1設計点あたり数時間のCFDを流量を変えて10〜15点計算し、マップを構築する。自動車メーカーはこのプロセスをスクリプト化し、一晩でコンプレッサーマップを自動生成するパイプラインを構築している。CFDで見つけたサージラインは実機試験と通常±5%以内で一致することが業界の経験則だ。
圧縮機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での形状最適化の自由度が大幅に拡大
ターボチャージャーCFDの商用ツール——NUMECAとANSYSの覇権争い
ターボチャージャーCFD市場では長年ANSYS CFX（ターボ機械特化機能）とStar-CCM+（自動車産業の強さ）が2強として競ってきた。2010年代後半から台頭したのが独NUMECA社のFINE/Turboで、BladeGenとの完全統合による翼型設計〜メッシュ生成〜CFDのワンストップ環境が中小メーカーに受け入れられた。NUMECAはその後Cadenceに買収され、大規模電子設計ツールとの統合が進んでいる。ターボチャージャーメーカー（BorgWarner, Honeywell Turbo, IHI）の多くは複数ツールのベンチマーク結果を蓄積しており、ツール依存のバイアスを排除するため意図的に複数ソルバーを使い分けている。
圧縮機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を連携させたリアルタイム性能監視
電動ターボチャージャーの挑戦——CFDが追いかける10万rpmの世界
電動ターボチャージャー（e-Turbo）は従来のターボチャージャーにモーター・ジェネレーターを組み込み、過渡応答性（ターボラグ）の改善と排気エネルギー回収を同時に実現する技術だ。最高回転数は10万〜15万rpmに達し、インペラの先端速度が550m/sを超えるものもある。この領域では遠心力による応力とCFD流体荷重の連成（流体構造連成）が不可欠で、単純な流体解析では動翼の変形によるクリアランス変化を見落とす。Formula 1では2022年規定からe-Turboが解禁され、各チームがCFD + FEA連成解析を競って開発している。
圧縮機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で追う
ターボチャージャーの不具合で最も多く報告されるのが特定の回転数で発生する異音だ。このNVH（騒音・振動・ハーシュネス）問題をCFDで追跡するには、定常解析では不十分で非定常（時間変化）解析が必要になる。翼通過周波数（ブレードパス周波数 = 回転数 × 枚数）付近の圧力変動を解析すると、スプリッタ翼の配置が不適切な場合に特定の翼枚数比で共鳴が起きるパターンが見えてくる。あるメーカーでは主翼6枚・スプリッタ6枚という構成を7枚・7枚に変更することで問題の音が消えた——CFDの非定常圧力スペクトルが設計変更の決め手となった。
なった
詳しく
報告