フィルム冷却
フィルム冷却の理論基礎
フィルム冷却の原理
先生、フィルム冷却って具体的にどういう仕組みで翼を冷やすんですか？
タービン翼の内部から冷却空気を翼表面に吹き出して、高温の主流ガスとの間に低温の「膜（フィルム）」を形成するんだ。この膜が断熱層として機能し、翼表面への熱流束を低減する。ガスタービンの入口温度が1500度Cを超える現代では不可欠な技術だよ。
冷却効率を表す指標はあるんですか？
断熱冷却効率（adiabatic film cooling effectiveness） $\eta$ が基本指標だ。
ここで $T_{\infty}$ は主流温度、$T_{aw}$ は断熱壁温度、$T_c$ は冷却空気温度。$\eta = 1$ なら壁が完全に冷却空気温度、$\eta = 0$ なら冷却効果なしということだ。
ブローイング比と運動量比
冷却空気をたくさん吹けば吹くほどいいんですか？
そう単純ではない。ブローイング比 $M$ が重要なパラメータだ。
$M$ が小さすぎると冷却空気が主流に呑まれて効果が薄い。$M$ が大きすぎると冷却ジェットが壁面から剥離（lift-off）して、かえって冷却効率が低下する。円孔の場合、$M \approx 0.5$〜$1.0$ が最適域と言われている。
それって孔の形状にも依存しますよね？
もちろん。Shaped hole（扇形孔、レイドバック孔）では出口が拡大するため冷却空気の運動量が下がり、lift-offが起きにくい。同じ $M$ でも円孔より $\eta$ が30〜50%向上するというデータがある。GE、Rolls-Royceなどのエンジンメーカーは独自の孔形状を特許化しているほどだよ。
フィルム冷却の誕生——1940年代のジェットエンジン黎明期
フィルム冷却の概念は1940年代後半、英国のRolls-RoyceがDerwentエンジンの燃焼器ライナーに多孔板から冷却空気を噴出させた試みに遡る。当初は「空気膜冷却（Air Film Cooling）」と呼ばれ、メカニズムは経験則頼みだった。理論的基盤はGoldstein(1971)の冷却効率式（η = (T∞ - Taw)/(T∞ - Tc)）の整備で確立し、以後50年でLES(大渦シミュレーション)による冷却孔周辺渦構造の解明まで進化した。現代のジェットエンジンタービン入口温度は1700℃を超え、ニッケル超合金融点（約1350℃）を上回るためフィルム冷却なしには一秒も動かせない。
フィルム冷却の数値計算手法
RANS vs LES の選択
フィルム冷却のCFDって、RANSで十分ですか？
主流のアプローチはRANSだけど限界がある。標準k-εは冷却ジェットの横方向拡散を過大評価する傾向があり、$\eta$ を過大予測しがちだ。Realizable k-εやSST k-ωモデルが推奨されるけど、それでもジェットのlift-off後の再付着を正確に捉えるのは難しい。
LESならうまくいくんですか？
LESは冷却ジェットと主流の混合を物理的に正しく再現できるので、$\eta$ の分布がRANSより実験に近い。ただし計算コストがRANSの100〜1000倍になる。実務では設計段階でRANS、最終検証でLESという使い分けが現実的だ。最近ではHybrid RANS/LES（DESやSBESモデル）が妥協案として使われている。
メッシュ戦略
冷却孔まわりのメッシュってどうすればいいですか？
冷却孔の直径 $D$ を基準にする。孔の周方向に最低20分割、孔内部の長さ方向にも20分割程度が目安。壁面第一層は $y^+ < 1$ を目指す。孔出口の下流域は $\Delta x / D \approx 0.1$ のメッシュ密度で少なくとも $x/D = 30$ まで確保するのが理想だ。
全翼面に何百もの孔がある場合、全部メッシュを切るのは非現実的ですよね？
そのとおり。実務では3つのアプローチがある。(1) 代表的な孔のみ詳細にメッシュを切り、他は質量流量境界条件で模擬する。(2) Source termモデルを使い、孔を個別にメッシュ化せずに冷却効果を導入する。(3) Ansys FluentのFilm Cooling Modelのような専用機能を使う。STAR-CCM+にもCooling Film機能がある。
OpenFOAMでは？
OpenFOAMには標準でフィルム冷却専用機能はないけど、codedFixedValueを使って孔出口に速度・温度プロファイルを処方する方法がある。あるいはAMI（Arbitrary Mesh Interface）で孔周辺だけ細かいメッシュの別領域を接続するのも手だよ。
フィルム冷却孔のメッシュ——「解像度不足」が最大の敵
フィルム冷却CFDで最も計算コストが嵩む原因は「冷却孔内部と主流の界面を十分に解像する」ことの難しさだ。孔径は通常0.5〜2mm程度と非常に小さいが、孔内部の乱流境界層や孔出口での噴流と主流の混合域は、少なくとも孔径に対して10〜20分割以上の解像度が必要とされる。実際の翼では数百〜数千個の冷却孔があるため、全孔を高解像度でモデル化するとメッシュ数が1億を超えることも珍しくない。そのため実務では「代表的な数個の孔を詳細解析して局所特性を把握し、翼全体はモデル化された境界条件で解く」という2スケール解析アプローチが採用される。LESは孔出口付近の混合精度が格段に高いが、計算コストはRANSの数十倍になる。
フィルム冷却の実務適用
実験データとの照合
フィルム冷却CFDの検証にはどんな実験データを使えばいいですか？
定番はGoldsteinらの単一円孔フィルム冷却実験（1968年）と、最近ではミシガン大学やオハイオ州立大学のshaped hole実験データだ。感温塗料（PSP/TSP）やIRサーモグラフィで取得した2次元 $\eta$ 分布と比較するのが標準的な検証手順だよ。
CFDと実験のずれってどれくらいが許容範囲ですか？
RANSで中心線上の $\eta$ が実験値と20%以内、横方向平均 $\eta$ が15%以内なら「妥当」とされることが多い。LESなら10%以内を目指せる。ただしlift-off領域やジェットの端部では乖離が大きくなる傾向がある。
設計パラメータと感度
設計で重要なパラメータって何ですか？
主要パラメータとその影響を整理しよう。
| パラメータ | 典型的範囲 | $\eta$ への影響 |
|---|---|---|
| ブローイング比 $M$ | 0.3〜2.0 | $M \approx 0.5$〜$1.0$ で最大 |
| 孔間隔比 $P/D$ | 3〜6 | 小さいほど横方向カバレッジ向上 |
| 噴射角 $\alpha$ | 20〜45度 | 小さいほど壁面密着性向上 |
| 孔形状 | 円孔、扇形、layback | Shaped holeで大幅改善 |
| 密度比 $DR$ | 1.0〜2.0 | 実機では$DR \approx 2$、実験では$CO_2$で模擬 |
噴射角20度ってかなり浅いですよね。製造できるんですか？
従来の放電加工（EDM）では30度が限界だったけど、最近のレーザー加工や3Dプリンティング（付加製造）で20度以下の浅い角度も実現可能になってきた。設計の自由度が大幅に広がっているよ。
多列孔の干渉効果
実機では複数列の冷却孔がありますよね。列間の干渉はどう扱うんですか？
上流列からのフィルムが下流列の吹き出しに影響を与える。一般に上流のフィルムが壁面近傍の温度を下げるので、下流列の $\eta$ は単列の場合より向上する。ただし上流ジェットがlift-offしていると、下流のジェットも巻き上げられやすくなる。このような多列干渉の正確な予測はLESでないと難しい場合が多いよ。
ジェットエンジン翼のフィルム冷却——0.1%の効率差が命取り
航空エンジンのタービン翼は1700℃を超えるガス中で動作するが、翼材料の耐熱温度は1100℃程度が限界だ。その差を埋めるのがフィルム冷却で、翼面に開けた無数の孔から冷却空気を吹き出して薄い保護膜を形成する。現代の高バイパス比エンジンでは、冷却効率ηが0.01（1%）上がるだけで燃料消費が数%改善できるとされる。一方でフィルム冷却空気はエンジン出力に寄与しない「ロス」でもある——冷却しすぎると効率が落ちる。この微妙なトレードオフを最適化するためにCFDが使われており、実機試験前にフィルム冷却孔の配置・角度・直径を何十パターンも計算して絞り込む工程は、エンジン開発の中で最も計算時間を消費するフェーズのひとつだ。
フィルム冷却のソフトウェア比較
主要ソルバーのフィルム冷却対応
各ソフトでのフィルム冷却モデリングの違いを教えてください。
主要ソルバーの対応状況はこうだ。
| ソルバー | 専用機能 | 推奨乱流モデル | 備考 |
|---|---|---|---|
| Ansys Fluent | Film Cooling Model (source term) | Realizable k-ε, SST k-ω | CHT連成可能。DESもサポート |
| Ansys CFX | 冷却孔のパラメトリックモデル | SST k-ω (デフォルト) | 回転系フィルム冷却に実績あり |
| STAR-CCM+ | Cooling Film機能 | SST k-ω, IDDES | ポリヘドラルメッシュで孔周りの品質確保が容易 |
| OpenFOAM | 専用機能なし | kOmegaSST | codedFixedValueやAMIで対応 |
| COMSOL | 専用機能なし | k-ε, k-ω | 小規模研究用途。大規模翼面計算には不向き |
FluentのFilm Cooling Modelってどういう仕組みですか？
冷却孔を個別にメッシュ化せず、壁面セルにmass/momentum/energy sourceを付加することでフィルム冷却効果を模擬する。孔の位置、直径、噴射角、ブローイング比をパラメトリックに指定するだけで済むので、設計初期段階の多案比較に向いている。ただし孔出口近傍の詳細な渦構造は再現できないから、最終的にはresolved（孔をメッシュ化した）計算で検証する必要がある。
大規模計算ではどのソルバーが有利ですか？
産業用途の大規模フィルム冷却CFDでは、STAR-CCM+のポリヘドラルメッシュとtrimmerメッシュの組合せが冷却孔周辺の自動メッシュ生成で優位性がある。FluentはMosaic meshingで同等の品質が得られるようになった。並列計算性能はどちらも十分高い。
乱流Prandtl数の影響
乱流Prandtl数 $Pr_t$ の設定って重要ですか？
非常に重要だ。RANSでは乱流熱流束を $q_t = -\rho c_p \nu_t / Pr_t \cdot \nabla T$ で評価するので、$Pr_t$ の値が $\eta$ に直結する。デフォルトの $Pr_t = 0.85$ では横方向拡散を過大評価しがちで、$Pr_t = 1.0$〜$1.2$ に上げると実験値に近づくことが多い。ただしこれは問題依存なので、感度分析を推奨するよ。
フィルム冷却CFDツール対決——ANSYS Fluent vs CFX vs StarCCM+
ガスタービンのフィルム冷却孔シミュレーションでは、各商用ツールの予測精度に顕著な差が生じる。ANSYSのFluent（圧力ベースソルバ）はBreaker-out孔の非定常剥離を捉えやすい一方、CFXの連成ソルバはFSI（流固連成）との統合が強み。StarCCM+はPolymesherによる複雑形状メッシュ生成で優位性を発揮するが、乱流モデル（SST vs v2-f）の選択がAdiabatic効率予測に±5%の誤差をもたらすことが論文では報告されている。ツール選択より物理モデルの精査が先決だ。
フィルム冷却の先端研究
LESによるフィルム冷却解析
LESでフィルム冷却を解くと何がわかるんですか？
冷却ジェットと主流の界面に発生するKelvin-Helmholtz渦やCounter-Rotating Vortex Pair（CRVP）の非定常挙動を直接捉えられる。CRVPはジェットを壁面から引き剥がす主因で、この渦の強度と寿命がフィルム冷却効率を決定するんだ。
LESのメッシュ要件はどれくらいですか？
単一孔で500万〜2000万セル、壁面解像を含めると $\Delta x^+ \approx 20$〜$50$、$\Delta y^+ < 1$、$\Delta z^+ \approx 10$〜$30$ が目安だ。時間刻みは$\Delta t \cdot u_{\infty} / D \approx 0.01$ 程度。統計的に定常な結果を得るには少なくとも50〜100 flow-through timeの計算が必要で、計算コストは膨大になる。
付加製造（AM）との融合
3Dプリンティングでフィルム冷却がどう変わるんですか？
従来の鋳造+放電加工では実現不可能だった複雑な冷却構造が作れるようになった。格子状（lattice）の内部冷却構造、マイクロチャネル、可変断面の冷却孔など、CFDで最適化した形状をそのまま製造できる。GEのLEAPエンジンやSiemens Energyのガスタービンで既に実用化されている。
CFDによる形状最適化も進んでいますか？
アジョイント法を使った冷却孔形状の最適化が活発に研究されている。Fluent Adjointは目的関数（例：翼表面の平均温度を最小化）に対して壁面形状の感度を効率的に計算できる。STAR-CCM+のDesign Managerと組合せてDOE（実験計画法）やKriging代替モデルを使った多目的最適化も実務で使われ始めているよ。
機械学習の活用
機械学習も使えますか？
使える。LESデータを教師データとして、RANSの乱流Prandtl数マップを機械学習で補正するアプローチが注目されている。また、CFDの代替モデル（サロゲート）としてCNNやGraph Neural Networkを使い、冷却孔配置と $\eta$ 分布の関係を高速に予測する研究も進んでいるよ。
アディティブ製造（AM）で実現する次世代フィルム冷却孔形状
従来のフィルム冷却孔は機械加工の制約から円形断面が主流だったが、金属3Dプリンティング（SLM/DMLS）の普及で卵形・ファン型・ダブル孔など複雑断面が実用段階に入った。CFDと実験の比較では、ファン型孔（Laid-back Fan-shaped）が円形に比べて冷却効率（ηad）を最大50%向上させるとNASAの研究で示されている。ただし積層造形特有の表面粗さ（Ra=10〜30μm）が境界層に与える影響も無視できず、粗さ修正壁関数を組み込んだCFD検証が最先端の研究課題となっている。
フィルム冷却のトラブル対応
冷却効率が実験値と大きくずれる
CFDで計算した $\eta$ が実験値の1.5倍もあるんですけど…
RANSでの典型的な症状だ。まず乱流モデルを確認しよう。標準k-εは冷却ジェットの横方向拡散を過大評価するから、SST k-ωかRealizable k-εに変更する。それでも過大なら乱流Prandtl数 $Pr_t$ を0.85から1.0〜1.2に上げてみよう。
逆にCFDの $\eta$ が実験値より低い場合は？
メッシュが粗すぎてジェットの壁面への再付着が解像できていない可能性がある。孔の下流域で $\Delta x / D > 0.5$ になっていたらメッシュを細かくすべきだ。また、孔のプレナム（供給室）を省略して一様流入条件を置いている場合、孔内の速度分布が実際と異なり $\eta$ が過小評価されることがある。
冷却孔内の逆流が発生する
孔の入口付近で逆流が起きてしまいます。
ブローイング比が小さい場合（$M < 0.3$）や、主流の圧力変動が大きい場合に起こりうる。まず境界条件の設定を確認しよう。冷却孔の入口にmass flow inletを使っている場合はpressure inletに変更してみる。逆流が物理的に妥当な場合もあるので、孔のL/D比や供給プレナムの圧力分布も含めて判断する必要がある。
非定常計算が安定しない
DESで計算しようとしたら発散してしまいました。
DESやLESは数値スキームの選択が重要だ。対流項にFirst Orderを使っていると数値拡散でLESモードがうまく機能しないし、Central Differenceだと安定しない。推奨はBounded Central DifferenceかBlended Second Order/Centralだ。Fluentなら「Bounded Second Order Implicit」、STAR-CCM+なら「Hybrid 2nd/Central」を選択しよう。
時間刻みの目安は？
CFL数を1以下に保つのが基本だ。孔の直径 $D$ と主流速度 $u_{\infty}$ を使って $\Delta t < D / (20 \cdot u_{\infty})$ 程度を初期値にして、CFL数の分布を確認しながら調整するのがよい。STAR-CCM+のConvective Courant Number場やFluentのCFL contourで可視化できるよ。
初期過渡をどれくらい捨てればいいですか？
冷却ジェットの主流方向のflow-through time $\tau = L_{domain} / u_{\infty}$ を基準に、最低10〜20 $\tau$ を初期過渡として捨て、その後30〜50 $\tau$ で時間平均を取るのが標準的なやり方だ。
フィルム冷却CFDで「効率が負になる」謎のバグ
フィルム冷却の解析で「冷却効率η（フィルム冷却効率）が負の値になる」という不思議なトラブルが報告されることがある。ηの定義はη = (T∞ - Taw)/(T∞ - Tc) なので、通常は0〜1の範囲になるはずだ。負になる原因の多くは「熱い主流ガスが冷却穴に逆流している」か「温度の参照値を間違えている」かのどちらかだ。逆流はブローイング比Mが大きすぎる場合（M ＞ 2程度）や主流の圧力勾配が強い場合に起きる。CFDでこれが発生したら、まず穴出口付近の速度ベクトルを可視化して逆流の有無を確認するのが最初の一手だ。実験では捉えにくい局所的な逆流をCFDで発見できた、というのがフィルム冷却解析の大きな価値のひとつでもある。
関連トピック
なった
詳しく
報告