ラジアルタービン
ラジアルタービンの理論基礎
概要
ラジアルタービンってターボチャージャーのタービン側ですよね？
そう。ラジアルインフロータービンは外周から内径方向に流れが入り、軸方向に排出される。小型で高膨張比が得られるため、自動車ターボや小型ガスタービンに広く使われている。
全対静効率
ラジアルタービンの効率はどう定義されますか？
全対静（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の精密化が重なった複合的な進化の産物だ。
ラジアルタービンの数値計算手法
モデル構成
ラジアルタービンの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も検討対象になる。
ラジアルタービン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%高くなる過大評価が生じることが論文で報告されている。
ラジアルタービンの実務適用
タービンマップの構成
タービンマップってどういう形式ですか？
横軸に膨張比（$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ルックアップテーブル入力として活用されている。
ラジアルタービンのソフトウェア比較
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という棲み分けが実務上の指針となっている。
ラジアルタービンの先端研究
高膨張比タービン
ラジアルタービンで膨張比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の位置依存性が顕著になる。
対策: 界面をホイール入口（径方向流れ）のやや上流に設定し、曲率変化の少ない位置を選ぶ。
出口ディフューザの処理
出口管で逆流が出て困っています。
ラジアルタービンの出口には強い旋回が残る。短いディフューザだと出口境界面で逆流が発生する。
対策:
- ディフューザ長さをホイール出口径の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倍になるが、精度向上の効果はこのコストを正当化する。
関連トピック
なった
詳しく
報告