タービン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に使われるソフトは？
ベッツ限界59.3%の謎——なぜ「100%」は不可能なのか
1919年にアルベルト・ベッツが導出したベッツ限界（16/27 ≈ 59.3%）は、風力タービンが風から取り出せるエネルギーの理論上の上限だ。直感的には「全部のエネルギーを取ればよい」と思えるが、そうすると下流の風速がゼロになり流れが止まって新しい風が入ってこなくなる。風速を完全に落とさず、適度に「通過させる」ことが最大出力の鍵だ。現代の大型風力タービンの実効率は45〜50%で、機械・電気損失を除けばほぼベッツ限界の80〜85%に達している。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でどう扱うんですか？
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では翼面の断熱壁面温度を出力して計算する。
アクチュエータディスク法の100年——BEM理論からLESへ
風力タービンの流体解析の最初の数学的モデルは、1920年代のBetz（ベッツ）の運動量理論に遡る。ロータを「推力を生成する無限薄の円盤」として扱うアクチュエータディスク法は、現在でも風力発電所全体のウェイク解析（ファームCFD）で広く使われる。2000年代以降はアクチュエータラインモデル（ALM）が登場し、個々の翼の揚力・抗力を体積力としてLES流れ場に印加することで、翼端渦の発生と崩壊まで再現できるようになった。BEMから始まった理論が100年でLESと融合した進化のストーリーは、CFDの発展史の縮図だ。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
タービン翼列のメッシュ
タービン翼列のメッシュは圧縮機と同じですか？
基本構造は同じだが、タービン特有の注意点がある。
- 後縁の薄さ: タービン翼は後縁が非常に薄い（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% |
ウィンドファームの「陰」——ウェイク干渉が隠す損失
洋上ウィンドファームでは後流（ウェイク）干渉による出力損失が全体の10〜20%に達することがある。2011年のHorns Rev洋上風力発電所（デンマーク）の衛星雲画像には、80基の風車が生み出すウェイクが数km下流まで続く様子が鮮明に映っており、CFDコミュニティに衝撃を与えた。その後のCFD研究でウェイク回復距離は大気安定度に強く依存することが判明し、中立・安定・不安定の大気境界層条件を適切に設定しない解析は信頼性が著しく低下する。実務ではWind Atlasデータとの整合を必ずチェックすることが業界標準となっている。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ 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ツール選定において他のターボ機械分野と異なる独特の生態系を持つ。デンマーク工科大学（DTU）発のEllipSysやRISØコードが欧州の主要メーカーで長年使われてきた一方、商用分野ではSIMPACK（MBD連成）やStar-CCM+（ファームCFD）が勢力を伸ばす。さらにWAsP（ウィンドアトラス解析）やFLAWS（疲労荷重解析）など風力専用ツールとのワークフロー統合が重要で、汎用CFDツール単独では対応できない業界固有の要件が多い。2020年代からはNREL（米国国立再生可能エネルギー研究所）のOpenFAST/OpenFOAM連成が業界標準の地位を確立しつつある。
選定で最も重要な3つの問い
- 「何を解くか」：タービンCFD解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
翼列干渉
タービンの動翼と静翼の干渉はどう評価しますか？
Sliding Meshの非定常計算で、上流ノズルのウェイクが下流ロータ翼面に衝突する過程を直接捕える。ウェイクによる翼面圧力変動がHCF（高サイクル疲労）の主因だ。
クロッキング効果
クロッキングって何ですか？
多段タービンで前段のノズルウェイクが後段のノズル前縁に当たるか翼間を通過するかで、後段の性能が変わる現象だ。ノズル同士の周方向相対位置（クロッキング位置）を最適化すると段効率が0.3～1ポイント改善される。
CFDでクロッキングを評価するにはどうしますか？
ノズルの相対位置を5～10段階に変えてSliding Meshの非定常計算を行い、各位置での時間平均効率を比較する。計算コストは大きいが設計効果も大きい。
翼面遷移
タービン翼面の遷移はCFDで予測できますか？
上流ウェイクの周期的な通過によってタービン翼面では「バイパス遷移」が支配的になる。Gamma-Theta遷移モデルでは定常のバイパス遷移を予測できるが、ウェイク誘起遷移の非定常効果を正確に捕えるにはSliding Mesh + 遷移モデルの組み合わせが必要だ。
タービン冷却の最前線
翼面冷却の研究トレンドは？
浮体式洋上風力のCFD——6自由度の嵐の中で
浮体式洋上風力発電機（FOWT）のCFD解析は、風力タービン界で最も計算負荷が高い問題のひとつだ。波・風・タービン回転・係留力学の連成が必要で、ロータが6自由度で動くプラットフォーム上に載っている。2020年代初頭の研究では、このような連成CFDのフル解析に数百のCPUコアで数週間という計算時間が報告されている。実用的には「アクチュエータラインモデル（ALM）」でロータを簡略化し、構造・係留ソルバーと弱連成させる手法が標準化されつつある。日本の洋上風力30GW目標に向けて、この分野の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が収束しない理由——大気境界層の「壁」
風力タービンCFDで初学者がはまる典型的な問題が「入口境界条件の設定ミス」だ。大気境界層では風速が高度のべき乗則（v ∝ z^α）に従い、地上粗度によってαが0.1〜0.4まで変わる。フラットな流入境界条件を設定すると、タービン前で境界層が計算領域内で発達してしまい、タービン位置の流入条件が設計とかけ離れる。正しい手順は「発達した大気境界層を先に別途計算し、その結果を入口条件として与える」か、解析的プロファイルと乱流量の整合した組み合わせを使うことだ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——タービンCFD解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告