圧縮機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単位まで最適化している。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
サージ予測のアプローチ
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万要素を超えることも珍しくない。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
解析ワークフロー
遠心圧縮機の典型的な解析フローを教えてください。
以下の手順が標準だ。
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の解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でも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での形状最適化の自由度が大幅に拡大
ターボチャージャーCFDの商用ツール——NUMECAとANSYSの覇権争い
ターボチャージャーCFD市場では長年ANSYS CFX（ターボ機械特化機能）とStar-CCM+（自動車産業の強さ）が2強として競ってきた。2010年代後半から台頭したのが独NUMECA社のFINE/Turboで、BladeGenとの完全統合による翼型設計〜メッシュ生成〜CFDのワンストップ環境が中小メーカーに受け入れられた。NUMECAはその後Cadenceに買収され、大規模電子設計ツールとの統合が進んでいる。ターボチャージャーメーカー（BorgWarner, Honeywell Turbo, IHI）の多くは複数ツールのベンチマーク結果を蓄積しており、ツール依存のバイアスを排除するため意図的に複数ソルバーを使い分けている。
選定で最も重要な3つの問い
- 「何を解くか」：圧縮機CFD解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
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の計算コストはどのくらいですか？
概算を示す。
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で「数値的チョーク」って何ですか？
メッシュの喉部面積が実際の幾何学的喉部面積と一致しないことで、チョーク流量がずれる現象だ。特に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の非定常圧力スペクトルが設計変更の決め手となった。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——圧縮機CFD解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告