タービンCFD解析
タービンCFDの理論基礎
概要
タービンCFDって圧縮機と何が違うんですか？
タービンは流体からエネルギーを取り出す側だ。流れが加速するため、圧縮機のような大規模剥離は起きにくい。その代わり、翼面冷却、二次流れ損失、遷音速衝撃波が主要な課題になる。
段仕事と等エントロピー効率
タービンの仕事はどう表されますか？
Euler方程式から出力と効率を定義する。
$h_{02s}$ は等エントロピー膨張後のエンタルピー。ガスタービンのHP段で $\eta_{is}=90\sim92\%$、LP段で$88\sim90\%$が現代の設計水準だ。
翼負荷係数
翼負荷の大きさはどう評価しますか？
Zweifel翼負荷係数が定番だ。
$s$:ピッチ、$c_x$:軸方向翼弦。$Z_w \approx 0.8$ が従来の最適値とされてきたが、最近の高負荷設計では$Z_w > 1.0$も研究されている。
ソフトウェア選択
タービンCFDに使われるソフトは？
タービン効率理論の礎——ランキンサイクルと熱効率（1859年）
蒸気タービンの熱効率を定量化するランキンサイクル（Rankine Cycle）の理論は、スコットランドの工学者William Rankine（1820-1872）が確立した。カルノー効率に対し、実際の蒸気タービンシステムが達成できる効率をランキンサイクル計算で求める手法は今でも発電プラント設計の基礎だ。Rankineが蒸気機関の熱力学的解析を論文にまとめた1859年から165年後の今、最新のガスタービンコンバインドサイクル（GTCC）発電では熱効率64%を達成するまでに発展した。この向上のおよそ1/3はCFDによる翼形最適化（高温域の冷却技術含む）が貢献しており、160年前の理論的枠組みが現代CFDと組み合わさって機能し続けている。
タービンCFDの数値計算手法
翼面冷却の重要性
タービン翼の冷却ってCFDでどう扱うんですか？
HPタービンの入口ガス温度は1500～1800℃に達し、翼材料の耐熱限界（Ni基超合金で約1000℃）を大幅に超える。内部冷却通路とフィルム冷却で翼面温度を下げるんだ。
冷却モデルの階層
冷却をCFDにどう取り込みますか？
精度とコストのトレードオフで複数の階層がある。
| レベル | モデル | 計算コスト | 精度 |
|---|---|---|---|
| L0 | 冷却流量なし（断熱壁） | 最低 | 冷却なしの基準評価 |
| L1 | Source Term（質量/エネルギー注入） | 低 | フィルム冷却の概算 |
| L2 | Discrete Hole（個別冷却孔BC） | 中 | フィルム効率の定量評価 |
| L3 | Resolved Cooling Holes（孔をメッシュ化） | 高 | 最高精度だが工数大 |
| L4 | CHT（流体+固体連成） | 最高 | 翼内温度分布まで予測 |
L3やL4は実用的ですか？
単一翼のL3/L4はSingleton計算として実用化されている。STAR-CCM+のCHTがこの用途で評価が高い。多段でL3/L4は現状では研究レベルだ。
フィルム冷却効率
フィルム冷却の効果はどう評価しますか？
断熱フィルム冷却効率で定義される。
$T_g$:主流ガス温度、$T_{aw}$:断熱壁面温度、$T_c$:冷却空気温度。$\eta_f = 0$ が冷却なし、$\eta_f = 1$ が完全冷却だ。CFDでは翼面の断熱壁面温度を出力して計算する。
ガスタービン翼CFDの数値設定——高温燃焼ガスの物性と熱伝達モデルの選択
ガスタービン翼のCFD解析では、作動ガス（高温燃焼ガス）の物性値の正確な設定が熱伝達予測の精度を左右する。燃焼後のガスはCO2・H2O・N2・O2の混合物で、温度依存の比熱Cp(T)・粘性mu(T)・熱伝導率k(T)を多項式近似かWSGGモデルで設定する必要がある。一般に使われる「1500K時の空気物性値を一定値として使う」近似は高温側で粘性を10〜15%過少評価し、これが翼表面の境界層厚さとHTC（熱伝達係数）の予測誤差の原因になる。NASAのCHT（Conjugate Heat Transfer）ベンチマーク実験では、物性の多項式近似を適切に設定したCFDで翼面HTC分布の実験との誤差を±8%以内に抑えた結果が示されており、物性設定の精度が高温翼CFDの基本であることが実証されている。
タービンCFDの実務適用
タービン翼列のメッシュ
タービン翼列のメッシュは圧縮機と同じですか？
基本構造は同じだが、タービン特有の注意点がある。
- 後縁の薄さ: タービン翼は後縁が非常に薄い（0.3～0.8mm）。O-gridの後縁周りに十分なセルが必要
- 冷却孔: L2/L3モデルでは冷却孔周りの局所細分化が必要
- 翼面の遷音速領域: 吸い込み面の超音速パッチと後縁衝撃波の解像
後縁が0.3mmだとメッシュも相当細かくなりますね。
後縁のO-gridは半径方向に少なくとも10セル、後縁直後のウェイク領域もメッシュを細かくする。TurboGridの trailing edge cutoff 機能で後縁形状を制御できる。
遷音速タービン翼列
タービンの流れは超音速になるんですか？
HPタービンでは翼間マッハ数1.1～1.3に達する。吸い込み面で超音速に加速した後、後縁から斜め衝撃波が発射される。この衝撃波が隣の翼に入射するTrailing Edge Shock Systemの正確な予測がCFDの精度を左右する。
衝撃波の解像にはどのくらいのメッシュが必要ですか？
衝撃波方向に直交するセルサイズが翼弦の0.5%以下、衝撃波前後に10セル以上が推奨だ。Adaptive Mesh Refinement（AMR）で衝撃波位置にメッシュを集中させるのも効果的だ。FluentやSTAR-CCM+のAMR機能が使える。
性能予測精度
タービンCFDの精度はどのくらいですか？
| 指標 | 精度 |
|---|---|
| 段効率（多段） | ±0.5～1.5ポイント |
| 翼面圧力分布 | 良好（実験と定性的に一致） |
| 翼面熱伝達係数 | ±10～20%（乱流モデル依存） |
| 後縁衝撃波位置 | 翼弦の±2% |
産業用蒸気タービンのCFD定期診断——シール摩耗とクリアランス変化の性能影響
産業用蒸気タービンは10〜20年の長期運転で翼端シールの摩耗・翼型の侵食（エロージョン）・翼端クリアランスの拡大が生じ、段効率が経年で0.5〜2ポイント低下する。この性能劣化をCFDで定量評価し、最適なオーバーホールタイミングを決定する「CFD診断」が大型火力発電所で実用化されている。解析では実機の定期点検データ（シール隙間の実測値、翼面プロファイル変化）をCFDモデルに反映し、設計値との比較で損失増加量を特定する。ある日本の電力会社の事例では、CFD診断で特定の段のシールクリアランスが設計値の3倍になっていることを予測・実機確認し、部分的なシール交換で全体オーバーホールと比べてコストを60%削減した事例が電力中央研究所の報告書で公開されている。
タービンCFDのソフトウェア比較
蒸気タービンのCFD
ガスタービンと蒸気タービンのCFDは何が違いますか？
蒸気タービンのLP段では膨張過程で蒸気が飽和線を超え、湿り蒸気（二相流）になる。湿り損失の予測が蒸気タービン特有の課題だ。
湿り蒸気ってどう扱うんですか？
Eulerian-Lagrangian法が一般的だ。蒸気（連続相）をEuler的に解き、水滴（分散相）をLagrangian粒子として追跡する。核生成による水滴発生と水滴の成長をモデル化する。
CFXの湿り蒸気モデル
CFXには蒸気タービン用のモデルがありますか？
CFXにはWet Steam Modelが搭載されている。非平衡凝縮（Spontaneous Nucleation）を考慮し、過冷却度からWilsonラインを予測できる。IAPWS-IF97の水蒸気物性テーブルも利用可能だ。
| 機能 | CFX | STAR-CCM+ | OpenFOAM |
|---|---|---|---|
| 湿り蒸気モデル | Wet Steam Model | Lagrangian＋カスタム | wetSteamFoam (研究用) |
| 物性テーブル | IAPWS-IF97 | IAPWS-IF97 | janaf or カスタム |
| 水滴追跡 | Eulerian Multiphase | Lagrangian | 限定的 |
長翼の課題
蒸気タービンのLP最終段の翼はとても長いですよね？
最新の蒸気タービンLP最終段翼は1.5m以上あり、翼端周速がマッハ1.5を超える。翼根から翼端で流れ状態が全く異なるため、3D設計が不可欠だ。翼端では超音速+湿り蒸気、翼根では亜音速+乾き蒸気という条件になる。
スパン方向でそんなに違うんですか。
翼高さ方向に50～100のスパン断面で翼型を変える3D Stackingが標準設計手法だ。CFDでスパン全体の効率分布を精密に評価する必要がある。
ガスタービンCFDツール——商用ツールと企業内製コードの産業界の棲み分け
世界のガスタービンメーカーは独自内製CFDコードと商用ツールを組み合わせる独特の生態系を持つ。Rolls-Royceの「Hydra」、GE Aerospaceの「WIND」、MTUの「TRACE」はそれぞれ50年以上の開発歴史を持つ企業内製コードで、翼型設計から認証試験まで信頼の高い内部ツールとして使われている。商用ツールとしてはANSYS CFXが航空エンジン翼型解析でのデファクト標準で、NUMECA FINE/TurboがヨーロッパのAerospace企業でよく使われる。学術・研究機関ではSU2（スタンフォード大学）やOpenFOAMも使われるが、実機設計認定の文書に使われるのは内製か商用の検証済みコードが大多数だ。ツールの信頼性は何十年もの実績による検証データの蓄積で決まる——新しいツールがいくら高機能でも信頼に足る実績がなければ採用されない現実がある。
タービンCFDの先端研究
翼列干渉
タービンの動翼と静翼の干渉はどう評価しますか？
Sliding Meshの非定常計算で、上流ノズルのウェイクが下流ロータ翼面に衝突する過程を直接捕える。ウェイクによる翼面圧力変動がHCF（高サイクル疲労）の主因だ。
クロッキング効果
クロッキングって何ですか？
多段タービンで前段のノズルウェイクが後段のノズル前縁に当たるか翼間を通過するかで、後段の性能が変わる現象だ。ノズル同士の周方向相対位置（クロッキング位置）を最適化すると段効率が0.3～1ポイント改善される。
CFDでクロッキングを評価するにはどうしますか？
ノズルの相対位置を5～10段階に変えてSliding Meshの非定常計算を行い、各位置での時間平均効率を比較する。計算コストは大きいが設計効果も大きい。
翼面遷移
タービン翼面の遷移はCFDで予測できますか？
上流ウェイクの周期的な通過によってタービン翼面では「バイパス遷移」が支配的になる。Gamma-Theta遷移モデルでは定常のバイパス遷移を予測できるが、ウェイク誘起遷移の非定常効果を正確に捕えるにはSliding Mesh + 遷移モデルの組み合わせが必要だ。
タービン冷却の最前線
翼面冷却の研究トレンドは？
タービン翼の最先端冷却——AIによるフィルム孔配置最適化とCFD検証
ガスタービン翼の冷却効率を最大化するために、AIと深層強化学習（Deep Reinforcement Learning）を使ったフィルム冷却孔の配置最適化が研究されている。数百個の孔の位置・直径・傾斜角の組み合わせは事実上無限の設計空間を持ち、従来のパラメトリック最適化では限界がある。GEとNASAの共同研究では、強化学習エージェントがCFD評価の報酬信号を使いながら10万回以上の試行で最適な孔配置パターンを探索し、冷却空気量20%削減しながら同等の冷却効果を実現した案を提示した。この手法で設計された実験翼型は3Dプリンティングで試作・検証され、CFD予測との乖離が±5%以内であることが確認されている。
タービンCFDのトラブル対応
熱伝達係数の不一致
CFDの翼面熱伝達係数が実験と合わないことが多いのですが…
翼面熱伝達予測はCFDで最も難しい項目の一つだ。典型的な不一致要因を整理しよう。
| 要因 | 影響 | 対策 |
|---|---|---|
| 乱流モデル | SSTで±15～25%の誤差 | Gamma-Theta遷移モデル追加 |
| y+の管理 | y+ > 2で熱伝達が過小 | y+ < 1を確保 |
| 入口乱流強度 | 燃焼器出口TI: 10～20% | 実験値を反映、デフォルト5%は過小 |
| フリーストリーム乱流減衰 | 翼前縁までにTIが減衰 | 乱流長さスケールも正しく設定 |
入口乱流強度が10～20%もあるんですか？
燃焼器出口では渦が残っているから、乱流強度が高い。これを5%で計算すると前縁付近の熱伝達が大幅に過小評価される。
後縁の熱伝達
後縁付近の熱伝達が特に合わないと聞きました。
後縁はウェイク領域と翼面の境界層が交差する複雑な流れ場だ。RANSでは後縁近傍の乱流構造を正確に再現できないことが多い。SASやSDESでこの領域の非定常渦を解像すると改善される。
CHT解析のTips
CHT（共役熱伝達）解析のコツを教えてください。
| Tips | 詳細 |
|---|---|
| 固体メッシュの整合 | 流体-固体界面で節点を一致させると精度向上 |
| 固体の熱伝導率 | Ni基超合金: 11～25 W/(m・K)、温度依存を考慮 |
| TBCの扱い | Thin Wall BCで薄い断熱コーティングを模擬 |
| 内部冷却通路 | 1D流れ網モデルで簡略化可能（CFXのBoundary Source Term） |
| 収束判定 | 翼面温度の変動が±1K以内で安定 |
タービンCFDで翼面熱伝達が実験の2倍——フィルム冷却孔周辺のメッシュ不足
ガスタービン翼のCFD解析で「翼面の熱伝達係数（HTC）が実験より2倍程度高い」という問題は、フィルム冷却孔周辺の局所メッシュ解像度不足が原因のケースが多い。冷却孔径D（典型的に0.5〜1mm）に対して、孔出口から下流10D以内は孔径の少なくとも1/20のセルサイズが必要で、これを満たさない粗いメッシュでは孔出口の速度・温度分布が著しく誤算される。さらに孔周辺に生成される腎臓渦対（Counter-Rotating Vortex Pair）がAdiabaticフィルム冷却効率の主要制御因子だが、渦コア径がD/2程度と小さいため解像度が不足すると渦が消えてしまう。実用的な対策は孔周辺のメッシュにAMR（適応細分化）を適用し、熱流束収束確認を孔間隔とノズルベーン配列の両方で実施する手順が推奨される。
関連トピック
なった
詳しく
報告