ラジアルタービン
理論と物理
概要
ラジアルタービンってターボチャージャーのタービン側ですよね？
そう。ラジアルインフロータービンは外周から内径方向に流れが入り、軸方向に排出される。小型で高膨張比が得られるため、自動車ターボや小型ガスタービンに広く使われている。
全対静効率
ラジアルタービンの効率はどう定義されますか？
全対静（total-to-static）効率が一般的だ。
出口の動圧を回収しない場合（自由排気）にはこの定義が適切だ。
速度比
速度比って何ですか？
インペラ周速 $U$ と等エントロピー膨張速度 $C_s$ の比だ。
最高効率は $U/C_s \approx 0.7$ 付近で得られる。これは翼入口で旋回成分が最適になる条件に対応する。
0.7って覚えやすい数字ですね。
ラジアルタービン設計の最も基本的な指標だ。CFDで得られた効率を $U/C_s$ でプロットし、ピークが0.7付近にあるか確認するのが最初のチェックポイントになる。
ターボチャージャー用の特殊性
ターボチャージャーのラジアルタービンに特有の課題は？
排気パルスへの応答だ。エンジンの排気は間欠的で、タービン入口圧力は数msの周期で大きく変動する。定常CFDでは時間平均の性能しか予測できないから、精度を求めるなら入口圧力をパルス状に変動させる非定常計算が必要だ。
ラジアルタービンの歴史——向心式タービンの設計理論確立（1930〜60年代）
ラジアルフロータービン（向心式タービン）の理論と設計手法を確立したのは1930〜40年代で、de Lavalのインパルス蒸気タービンを基礎にGehringら（1930年代）が高効率向心型タービンの翼形設計理論を発展させた。現代のターボチャージャー用ラジアルタービンの基本的な翼形設計手法——バックスイープ角、出口翼角の最適化——は1960年代のNASAとGE Turbineの共同研究（Rohlik 1968等）で完成された。それから60年、この基礎理論の上にCFDによる精密設計が重なり、現代のVGTタービン効率は86〜90%という驚異的なレベルに達している。流体機械の設計は一朝一夕に生まれず、60年超の実験と理論の蓄積の上にCFDの精密化が重なった複合的な進化の産物だ。
各項の物理的意味
- 時間項 $\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モデルはどう組みますか？
典型的な構成は以下だ。
- ボリュート/スクロール: 静止域、非軸対称
- ノズルベーン（VGT）: 静止域、翼列
- タービンホイール: 回転域
- ディフューザ/出口管: 静止域
界面処理はボリュート-ノズル間がGGI（ピッチ差なし）、ノズル-ホイール間がFrozen RotorまたはSliding Meshだ。
VGTって可変ノズルですよね。開度を変えた解析もしますか？
する。VGTの開度を5～10段階に変えて各開度でのマップを作成する。ノズル角度の変更はBladeGenやCADでパラメトリックに行い、自動メッシュ→CFXのバッチ計算で連続処理する。
メッシュ生成
ラジアルタービンのメッシュで注意すべきことは？
領域 手法 注意点 ボリュート 非構造（テトラ+プリズム） 舌部の細かいメッシュ ノズルベーン TurboGrid or 非構造 可変開度に対応するメッシュ戦略 タービンホイール TurboGrid（構造格子） 翼間流路の曲率、スプリッタ翼対応 出口ディフューザ 構造 or 非構造 旋回流の減衰を十分に解像
タービンホイールはスプリッタ翼がありますか？
ラジアルタービンではスプリッタがない場合が多いが、一部の高性能設計ではスプリッタを設けて翼負荷を低減する。TurboGridではスプリッタ付きのトポロジにも対応している。
乱流モデル
ラジアルタービンに適した乱流モデルは？
SST k-omegaが標準だ。タービン翼面の遷移が重要な場合はGamma-Theta遷移モデルを追加する。排気脈動の非定常計算ではSASやSDESも検討対象になる。
Coffee Break よもやま話
ラジアルタービンCFDの数値設定——遷音速流と後縁衝撃波のメッシュ解像度
小型ガスタービンやターボチャージャーに使われるラジアルタービンは、出口が音速付近（Ma≈0.8〜1.0）になることが多く、遷音速の適切な数値処理が精度を決める。翼後縁から発生する斜め衝撃波（Oblique Shock）を正確に捉えるためには、後縁近傍メッシュを後縁厚さの少なくとも1/5のセルサイズに設定し、衝撃波角度に沿ったメッシュアライメントが必要だ。数値スキームはRoe Flux Differencing SchemeかHLLC（Harten-Lax-van Leer-Contact）スキームが衝撃波キャプチャ精度で推奨され、低Ma数域（ブレード圧力面）ではPreconditioning（低マッハ数前処理）が必要になる。これらの設定を怠ると、タービン効率のCFD予測が実験より3〜5%高くなる過大評価が生じることが論文で報告されている。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
ラジアルタービンのメッシュで注意すべきことは？
| 領域 | 手法 | 注意点 |
|---|---|---|
| ボリュート | 非構造（テトラ+プリズム） | 舌部の細かいメッシュ |
| ノズルベーン | TurboGrid or 非構造 | 可変開度に対応するメッシュ戦略 |
| タービンホイール | TurboGrid（構造格子） | 翼間流路の曲率、スプリッタ翼対応 |
| 出口ディフューザ | 構造 or 非構造 | 旋回流の減衰を十分に解像 |
タービンホイールはスプリッタ翼がありますか？
ラジアルタービンではスプリッタがない場合が多いが、一部の高性能設計ではスプリッタを設けて翼負荷を低減する。TurboGridではスプリッタ付きのトポロジにも対応している。
ラジアルタービンに適した乱流モデルは？
SST k-omegaが標準だ。タービン翼面の遷移が重要な場合はGamma-Theta遷移モデルを追加する。排気脈動の非定常計算ではSASやSDESも検討対象になる。
ラジアルタービンCFDの数値設定——遷音速流と後縁衝撃波のメッシュ解像度
小型ガスタービンやターボチャージャーに使われるラジアルタービンは、出口が音速付近（Ma≈0.8〜1.0）になることが多く、遷音速の適切な数値処理が精度を決める。翼後縁から発生する斜め衝撃波（Oblique Shock）を正確に捉えるためには、後縁近傍メッシュを後縁厚さの少なくとも1/5のセルサイズに設定し、衝撃波角度に沿ったメッシュアライメントが必要だ。数値スキームはRoe Flux Differencing SchemeかHLLC（Harten-Lax-van Leer-Contact）スキームが衝撃波キャプチャ精度で推奨され、低Ma数域（ブレード圧力面）ではPreconditioning（低マッハ数前処理）が必要になる。これらの設定を怠ると、タービン効率のCFD予測が実験より3〜5%高くなる過大評価が生じることが論文で報告されている。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
タービンマップの構成
タービンマップってどういう形式ですか？
横軸に膨張比（$p_{01}/p_2$）、縦軸に換算質量流量 $\dot{m}\sqrt{T_{01}}/p_{01}$ をとる。各回転数ラインを描き、効率をパラメータとして重ねる。
何点くらい計算が必要ですか？
回転数5水準 × 膨張比6～8点で30～40運転点が標準だ。VGTの場合はさらにノズル開度ごとにマップが必要になるから、全体で100点以上になることもある。
膨張比の変化方法
CFDで膨張比を変えるにはどうしますか？
入口全圧を固定して出口静圧を変えるのが一般的だ。
- 入口: 全圧 $p_{01}$（固定）、全温 $T_{01}$（固定）
- 出口: 静圧 $p_2$ を大気圧から段階的に下げる
膨張比が大きくなるとタービンがチョークに近づき、流量が飽和する。その状態を正確に予測するには翼間喉部のメッシュ解像度が重要だ。
実験との比較精度
CFDの精度はどのくらいですか？
| 指標 | 精度 |
|---|---|
| 質量流量パラメータ | ±2～3% |
| 全対静効率 | ±1～3ポイント |
| チョーク流量 | ±2% |
| 排気温度 | ±5～15K |
排気温度の誤差が大きいですね。
出口の温度分布は旋回やウェイクの影響で非常に不均一だ。測定位置とCFDの評価面を正確に合わせないと大きな乖離が出る。面積平均ではなく質量平均で比較することが重要だ。
自動車ターボチャージャーのラジアルタービン——可変ノズル（VGT）のCFD設計
可変ジオメトリターボ（VGT: Variable Geometry Turbocharger）は排気ガスの流れ方向をノズルベーン角で制御し、低回転域から高回転域まで最適なタービン入口条件を維持する。ノズルベーン角の変化に伴うタービン特性変化（効率マップ）の予測にCFDが不可欠だ。実務では5〜7段階のノズル開度（15°〜80°）でそれぞれRANS（フローズンローター）計算を行い、効率・流量係数・圧力比の3次元特性マップを生成する。BMW・Daimler・Hondaが公開した研究では、VGTのCFD特性マップと実機試験台の計測値の差は効率で±2〜3ポイント以内に収まっており、1Dエンジン性能シミュレーション（GT-SUITE）へのCFDルックアップテーブル入力として活用されている。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
VGT（可変ノズルターボ）
VGTのCFD解析で気をつけることは？
VGTはノズルベーンの開度を変えることでタービンの流量-圧力特性を制御する。CFDでは各開度でノズルの幾何形状を変更してメッシュを再生成する必要がある。
メッシュの再生成が大変ですね。
STAR-CCM+やFluent Meshingの自動メッシュなら形状変更→メッシュ再生成がスクリプト化しやすい。TurboGridでもジャーナルファイルで翼角度をパラメータとして変更できる。
排気脈動の非定常解析
エンジンの排気パルスを入れた計算はどうやりますか？
入口全圧を時間関数として与える。4気筒エンジン3000rpmなら、排気周期は10ms（100Hz）で各パルスのピーク圧力比は定常の2～3倍に達する。
Sliding Meshの非定常計算で、最低10パルス分（100ms）を計算して周期的定常に達するのを確認する。タイムステップは0.5～1度/step（回転角度基準）が推奨だ。
定常と非定常で効率は変わりますか？
一般に非定常パルス入力の時間平均効率は定常効率より2～5ポイント低くなることが報告されている。瞬間的にoff-design条件で運転するためだ。
ツイン/ツインスクロールターボ
ツインスクロールターボの解析は特別な配慮が必要ですか？
2つのスクロール入口にそれぞれ異なる排気パルスを与える必要がある。モデルは全周が必須で、セクターモデルは使えない。計算コストは単一スクロールの3～5倍だ。
| 構成 | セル数 | コア数 | 10パルス分 |
|---|---|---|---|
| シングルスクロール (Frozen Rotor) | 300万 | 32 | 4時間 |
| シングルスクロール (Sliding Mesh) | 300万 | 64 | 24時間 |
| ツインスクロール (Sliding Mesh+パルス) | 800万 | 128 | 3～5日 |
ラジアルタービンCFDツール——CFX TurboSystemとCONVERGE CFDの適用場面
ラジアルタービンCFD解析のツール選択は解析目的によって異なる。ANSYS CFX TurboSystemは翼列の性能マップ生成（P-Q-eta）に最も使いやすく、ボリュート込みのフルモデル解析も実績が豊富だ。翼形最適化をメインとするなら、アジョイント感度計算機能（CFXのAdjoint Solver）と組み合わせることで少数の反復で最適翼形を探索できる。一方CONVERGE CFD（Richards-Strauss社）は複雑な内燃機関ターボチャージャーのマルチフィジクス解析に強みがあり、排気ガスの化学反応・煤粒子輸送・タービンの熱流体を一括解析する事例がターボチャージャー開発会社から報告されている。ラジアルタービン単体の設計検証ならCFX、システムレベルのターボチャージャー連成ならCONVERGEという棲み分けが実務上の指針となっている。
選定で最も重要な3つの問い
- 「何を解くか」：ラジアルタービンに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
高膨張比タービン
ラジアルタービンで膨張比4:1以上は可能ですか？
可能だが、ノズル出口やホイール入口で超音速流れが発生する。衝撃波-境界層干渉による損失が増大するため、設計の難易度は高い。
CFDで超音速領域はどう扱いますか？
圧縮性を完全に考慮したソルバーが必須だ。CFXの結合型ソルバーは遷音速に強いが、強い衝撃波ではFluentの密度ベースソルバーのほうが安定する場合がある。
ノズルレスラジアルタービン
ノズルベーンなしのタービンもありますよね？
シンプルなターボチャージャーではボリュートが直接タービンホイールに流入するノズルレス設計が多い。この場合、ボリュートの断面積変化で流れ角を制御する。CFDでボリュート出口の流れ角分布がホイール性能に直結するから、ボリュートのメッシュ品質が重要だ。
最近の研究動向
ラジアルタービンの研究で注目のトピックは？
3Dプリンティングでタービンホイールを作れるんですか？
Inconel 718やHastelloy Xの粉末積層で製造されている。CFDで最適化された複雑な翼断面を加工制約なしで実現できるのが最大の利点だ。
ラジアルタービンの最前線——波動ロータ（Wave Rotor）との組み合わせ
超高効率のエネルギー変換を目指す「波動ロータ（Wave Rotor）」とラジアルタービンを組み合わせた複合ターボシステムが研究フロンティアで注目されている。波動ロータは圧力波を利用して圧縮・膨張を直接行う回転装置で、ラジアルコンプレッサーやタービンより高い熱力学的効率が理論上可能だ。波動ロータでガスタービンの燃焼器入口圧力を追加増圧することで、SFC（比燃料消費率）を15〜20%改善できると試算されている。CFD（Non-Linear Harmonic Method）によるロータ内の非定常圧力波の伝播と波動ロータとの干渉解析が、この概念実証の核心研究課題だ。
トラブルシューティング
チョーク流量の不一致
CFDのチョーク流量が実験と合わないのですが…
チョーク流量はノズルとホイールの喉部面積で決まる。メッシュが粗いと喉部の有効面積が変わる。以下を確認しよう。
1. ノズル喉部の実効面積をメッシュから計測し、CADの設計値と比較
2. ホイール入口の相対マッハ数分布を確認し、M=1のラインが物理的に正しい位置にあるか
3. 3水準のメッシュでチョーク流量の収束を確認
ノズル-ホイール界面
ノズルとホイールの界面でトラブルが起きやすいですか？
ラジアルタービンでは界面が環状面（cylindrical surface）になり、流れが径方向から軸方向に急激に向きを変える。この曲率の大きい領域で界面を設定すると、Frozen Rotorの位置依存性が顕著になる。
ノズルとホイールの界面でトラブルが起きやすいですか？
ラジアルタービンでは界面が環状面（cylindrical surface）になり、流れが径方向から軸方向に急激に向きを変える。この曲率の大きい領域で界面を設定すると、Frozen Rotorの位置依存性が顕著になる。
対策: 界面をホイール入口（径方向流れ）のやや上流に設定し、曲率変化の少ない位置を選ぶ。
出口ディフューザの処理
出口管で逆流が出て困っています。
ラジアルタービンの出口には強い旋回が残る。短いディフューザだと出口境界面で逆流が発生する。
対策:
- ディフューザ長さをホイール出口径の3～5倍に延長
- 出口にOpening BCを適用
- 旋回の減衰が十分になるまで計算領域を延長
高温ガスの物性
排気ガスの物性はどう設定しますか？
理想気体でγ=1.33～1.35（排気組成による）、分子量28.5～29程度を使う。温度範囲が600～1000℃と広いから、比熱の温度依存性を考慮するのが望ましい。CFXではNASA形式の多項式で温度依存比熱を設定できる。
ラジアルタービンCFDで効率が実測の5%低い——入口ボリュートの影響
ラジアルタービンのCFD解析で「効率が実験より5%程度低く計算される」という系統誤差は、計算領域にボリュート（Volute/Scroll）を含めていないことが原因のケースが多い。翼型単体（Blade Passage）だけを解析すると、ボリュートから翼車入口への流れの非一様性（流速・流入角の周方向分布）が無視され、実際より理想的な入口条件が仮定される。ボリュートがあると特定の周方向位置で局所的な失速が生じ、これが全体効率を下げる。診断手順：①ボリュート込みの全体モデルと翼型単体解析の結果を比較し、入口全圧の周方向一様性を確認する。②入口全圧周方向分布の標準偏差が平均の±2%を超える場合はボリュート込み解析が必要。ボリュートモデル化によって計算コストは3〜5倍になるが、精度向上の効果はこのコストを正当化する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ラジアルタービンの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告