ウィンデージ損失
ウィンデージ損失の理論基礎
概要
ウィンデージ損失って何ですか？ あまり聞き慣れない用語ですが。
回転体と静止壁面の間の空気（または流体）の摩擦による動力損失だ。高速回転機械では無視できない大きさになる。タービン発電機、フライホイール、高速モータなどで特に重要だ。
回転ディスクの摩擦損失
具体的にどのくらいの損失ですか？
閉空間内の回転ディスクの摩擦損失は以下で表される。
$C_M$はモーメント係数で、回転レイノルズ数 $Re_\theta = \rho \omega r^2 / \mu$ と隙間比 $s/r$ の関数だ。
$\omega^3$ に比例するんですか。回転数の影響が大きいですね。
そう。回転数を2倍にするとウィンデージ損失は8倍になる。だから高速回転機械では支配的な損失源になりうるんだ。
流れレジーム
隙間の流れにも層流や乱流があるんですか？
ある。Daily & Nece (1960) の分類が広く使われる。
| レジーム | $Re_\theta$ | 隙間比 | 特徴 |
|---|---|---|---|
| I: 層流・統合BL | 低 | 小 | 両壁面BLが結合 |
| II: 層流・分離BL | 低 | 大 | BLが独立、コア回転 |
| III: 乱流・統合BL | 高 | 小 | 最も一般的な工業条件 |
| IV: 乱流・分離BL | 高 | 大 | コア回転あり |
産業用ターボ機械はほとんどがレジームIIIまたはIVだ。
$C_M$ の経験式
モーメント係数の経験式はありますか？
レジームIIIの場合、Dailyの式が代表的だ。
CFDの結果をこの経験式と比較することで、計算の妥当性を簡便に検証できる。
回転円板の流体力学——von Karmanのディスク解析（1921年）とストークスの先駆
静止した流体中で回転する円板の流れを解析したのはTheodore von Karman（1921年）だ。彼はロータのカルマン渦だけでなく、回転ディスクの境界層（von Karman境界層）の厳密解析解を導き、ディスクのトルクと渦粘性の関係を定式化した。Cochran(1934)がこの解をより高精度に計算し、Daily & Nece(1960)が実験でモデルを拡張して今日の相関式の基礎を作った。現代のターボ機械ディスクキャビティCFDは100年前のvon Karmanの解析を出発点とし、多段ステージ・多体回転・ホットガス侵入といった複雑な三次元効果を数値的に追加した発展形だ。数学的アプローチで革命的業績を残した科学者の仕事が、100年後のジェットエンジン冷却設計に直接応用されている点は特筆に値する。
ウィンデージ損失の数値計算手法
CFDモデルの構成
ウィンデージ損失をCFDで計算するにはどうモデル化しますか？
回転ディスクと静止壁の隙間空間を流体領域としてモデル化する。
- 回転壁面: ディスク表面（No-Slip、回転速度を指定）
- 静止壁面: ケーシング内面（No-Slip、速度ゼロ）
- 入口/出口: 隙間へのリーク流れがある場合は設定
軸対称であれば2D軸対称モデルで十分だ。3D効果（ボルト頭、冷却孔など）がある場合はセクターモデルを使う。
乱流モデルの選択
どの乱流モデルが適していますか？
SST k-omegaが推奨だ。回転壁面近傍のBLが重要だから、Low-Re解法でy+ < 1を確保する。k-epsilonは壁近傍の回転効果を正確に捕えにくいため推奨しない。
| 乱流モデル | ウィンデージ予測精度 | 備考 |
|---|---|---|
| SST k-omega (Low-Re) | ±5～10% | 推奨 |
| k-epsilon + 壁関数 | ±15～25% | 壁近傍精度不足 |
| RSM (Reynolds Stress) | ±3～8% | 旋回効果を捕える、コスト大 |
| LES | ±2～5% | 最高精度だがコスト大 |
RSMが良いんですか？
回転ディスク流れでは旋回方向と径方向でレイノルズ応力の異方性が強い。RSM（BSL-RSM等）はこの異方性を直接モデル化するため、SST より精度が上がる場合がある。ただし計算コストが1.5～2倍になる。
メッシュ要件
隙間の中にどのくらいのメッシュが必要ですか？
隙間方向（軸方向）に40～60セル、径方向に100～200セルが目安だ。両壁面にプリズム層を配置し、y+ < 1を確保する。2D軸対称なら数万セルで十分だ。
ウィンデージ損失の計算——Daily-Nece相関式とCFDの精度比較
ローテーティングディスクのウィンデージ損失トルクMは、Daily & Nece(1960)の実験相関式から算出できる。4つの流れ体制（G-Reomegaの組み合わせ）それぞれに異なる係数式が適用され、乱流域（体制III,IV）ではM はロータ回転数の2乗・半径の5乗に比例しReの1/5乗に反比例する形になる。この相関式はエンジン・タービン設計では今でも初期見積もりに広く使われているが、適用範囲（G=0.02〜0.2、Reomega=10^4〜10^7）を逸脱した条件では誤差が30%以上になることがある。CFDとの比較研究では、ディスク間隙G<0.02（狭い間隙）やディスクに窪み・穴がある実形状では、Daily-Nece式がCFDより10〜20%低い損失を予測する傾向があると複数の論文が報告している。形状が複雑なエンジンディスクでは、相関式で初期見積もり後に必ずCFDで精密確認する手順が標準的だ。
ウィンデージ損失の実務適用
ロータ系の温度上昇
ウィンデージ損失で温度が上がるんですか？
損失は全て熱になる。密封空間では冷却が不十分だと温度が上昇し続ける。蒸気タービンの車室内やモータの空隙で問題になる。
冷却流量 $\dot{m}_{cooling}$ が少ないと温度上昇が大きくなる。
CFDで温度分布も予測できますか？
エネルギー方程式を有効にしてCHTで解けば、ロータ表面温度分布も予測できる。ただし対流だけでなく輻射も考慮する必要がある場合がある（高温タービン）。
シール流れとの連成
ラビリンスシールのリーク流れもウィンデージに影響しますか？
大きく影響する。ラビリンスシールを通過するリーク流れがディスク間の隙間に流入すると、コア回転率が変化してウィンデージ損失も変わる。CFDではシールと隙間を一体でモデル化するのが望ましい。
実験との検証
ウィンデージ損失のCFD結果はどう検証しますか？
| 検証方法 | 精度 |
|---|---|
| Daily & Nece経験式 | ±10%（単純ディスク） |
| トルクメータ実測 | 最も正確だが実験コスト高 |
| 温度上昇からの逆算 | ±15%（熱損失推定の不確実性） |
単純ディスクならCFDの精度は±5～10%で経験式とよく一致する。突起物やボルト頭がある場合は形状を正確にモデル化すれば±10～15%の精度が得られる。
航空エンジンのディスクキャビティ——二次流れ冷却システムのCFD熱設計
ジェットエンジンの高圧タービンでは、ホットガスが回転ディスクとケーシングの間の「ディスクキャビティ（Disk Cavity）」に侵入するとディスクが過熱・破損する危険がある。これを防ぐために「二次流れ冷却システム（Secondary Air System: SAS）」が圧縮機から空気を抽気してキャビティに注入し、ホットガス侵入を防ぐパージ流れを形成する。このパージ流量が多すぎるとエンジン効率が低下し、少なすぎるとディスクが危険な温度に達する。CFD（Rotating Disk + Compressible RANS）でキャビティ内の圧力・温度・流量の最適設計を行うことで、パージ流量を最小化しながら熱安全率を確保するトレードオフ設計が航空エンジンメーカーの重要なCFD適用分野となっている。
ウィンデージ損失のソフトウェア比較
損失低減の方法
ウィンデージ損失を減らすにはどうすればいいですか？
以下の方法がある。
| 方法 | 原理 | 効果 |
|---|---|---|
| 隙間縮小 | 流体体積の減少 | レジームIII→I移行で$C_M$低下 |
| 低密度ガス充填 | $P_w \propto \rho$ | Heで空気の1/7に |
| 真空化 | 密度低下の極限 | 大幅低減だがシール要 |
| 表面粗さ低減 | 摩擦係数低下 | 研磨面で5～15%低減 |
| ロータ表面の突起除去 | 圧力抗力の低減 | ボルト頭カバーで20～30%低減 |
フライホイールは真空化されてますよね？
そう。フライホイールエネルギー貯蔵装置では10Pa以下の真空中で運転し、ウィンデージ損失をほぼゼロにしている。ただし真空シールの設計が難しい。
CFDによるパラメトリック評価
CFDで損失低減策を評価する手順は？
2D軸対称モデルなら1ケース数分で計算できるから、パラメトリックスタディに最適だ。隙間比、回転数、ガス種を変えて$C_M$のマップを作成し、最適な設計条件を探索する。
ソルバーの選択
どのソルバーがウィンデージ解析に向いていますか？
ディスクキャビティCFDの専用ツール——国際的なSASコンソーシアムとCFX実績
航空エンジンの二次流れシステム（SAS）解析では、産業界と学術界が連携した国際研究プログラムが進められている。欧州のRFICAS（Internal Cooling Flows and Turbine Airfoils Study）、米国のNASA Turbofan SAS研究などが産業界標準解析手法の確立に貢献してきた。ツール面ではANSYS CFXがロールスロイス・GE・プラット&ホイットニーでの実績が多く、ターボ機械分野でのデファクト標準に近い。特にCFXのSST乱流モデルはディスクキャビティの旋回流（Swirl）予測で他モデルより一貫して高い精度を示しており、ディスクキャビティ内の旋回係数（Swirl Ratio）予測で±5%の精度が実績データとして蓄積されている。OpenFOAMのrotating Frameworkも研究用途で使われるが、商用ライセンス認証を要する設計業務ではCFXが選ばれる傾向が続いている。
ウィンデージ損失の先端研究
テイラークエット流れ
回転ディスク間の流れって渦ができるんですか？
同心円筒間の流れ（テイラークエット流れ）ではテイラー渦が発生する。回転レイノルズ数がTaylor数の臨界値を超えると、軸対称のトロイダル渦が自発的に形成される。
テイラー渦はウィンデージ損失にどう影響しますか？
渦による運動量輸送が増大し、$C_M$が急増する。テイラー渦の遷移をCFDで正確に予測するには、メッシュ解像度と乱流モデルの選択が重要だ。LESやDNSでは渦構造を直接解像できるが、RANSでは渦の発生は捉えられるものの構造の詳細は失われる。
軸方向流れのあるテイラー渦
シールにリーク流れがあるとテイラー渦はどう変わりますか？
軸方向流れが加わるとスパイラルテイラー渦になる。軸方向レイノルズ数と回転レイノルズ数の比で流れレジームが決まる。CFDではこの遷移を非定常計算で追うことで、シール内部の流れ構造と損失の関係を明らかにできる。
最近の研究
ウィンデージ損失研究の最前線は？
ウィンデージ損失の最前線——超高速フライホイールと真空ケーシング設計
フライホイールエネルギー貯蔵（FES: Flywheel Energy Storage）システムでは、回転体が毎分40,000〜100,000回転で駆動するため、周速度が500〜1000m/sに達しウィンデージ損失が全損失の主要部分になる。真空ケーシング（0.1〜1 Pa）内での回転体の空力損失は、稀薄気体（分子流）の特殊な摩擦モデルが必要で、通常の連続体CFD（Knudsen数Kn>0.1）が成立しない。DSMCシミュレーション（Direct Simulation Monte Carlo）による分子論的アプローチか、格子ボルツマン法（LBM）が適用される研究フロンティアだ。フライホイールの損失最小化で超高速回転体形状（スポーク配置・ハブ形状）のCFD最適化が行われており、ウィンデージ損失を50%削減した報告が超電導磁気軸受フライホイールの研究論文に掲載されている。
ウィンデージ損失のトラブル対応
2D vs 3Dモデル
ウィンデージ解析は2Dで十分ですか？
場合による。
| 条件 | 推奨 |
|---|---|
| 軸対称ロータ（滑らかな面） | 2D軸対称で十分 |
| ボルト頭や冷却孔がある | 3Dセクターモデル |
| ラビリンスシール付き | 2D軸対称（シール歯が軸対称なら） |
| 非軸対称形状（翼付きディスク） | 3Dフルモデル or セクター |
2Dと3Dで結果はどのくらい違いますか？
滑らかなディスクなら差は5%以内だ。ボルト頭が突出している場合は3Dでないとドラッグが過小評価される。ボルト頭1つで$C_M$が10～30%増加することがある。
収束判定
ウィンデージ解析の収束はどう判定しますか？
ロータ表面のトルクをモニターする。定常計算ならトルクの変動が0.1%以内に収まれば収束だ。
よくある間違い
ウィンデージ解析でやりがちなミスは？
2. 層流計算: 産業機械のRe_thetaは$10^6$以上で確実に乱流。層流設定だと$C_M$が桁違いに低くなる
3. コア回転率の無視: 隙間内の流体はディスクと同じ速度では回転しない。コア回転率（典型的に0.3～0.5）を正しく予測するには十分なメッシュ解像が必要
4. 突起物の省略: ボルト、配線ガード等の突起を省くと損失を大幅に過小評価する
コア回転率って何ですか？
ディスク間の中央（コア）の流体がディスク角速度の何割で回転しているかを示す指標だ。隙間が大きいとコア回転率は0.3～0.4程度で、ディスクの半分以下の速度で回転する。この値がウィンデージ損失の大きさを左右する。
ウィンデージ損失CFDが実測の2倍——ディスク間隙の流れモデルの選択ミス
ターボ機械のディスク・ハウジング間の「ウィンデージ損失（Windage Loss）」のCFD予測が実測の2倍になる場合、ロータリングの流れ体制の判断ミスが多い。ディスク間の流れはGap Ratio G = s/R（s:間隙、R:半径）と回転Reynolds数Reomega = omega*R^2/v で4つの体制（I〜IV）に分類される（Daily & Nece 1960）。G<0.01の狭い間隙でReomega>10^5の乱流域（体制IV）では境界層が合流し、流速分布が大きく変わる。この体制判別を誤るとウィンデージトルクの予測が2〜3倍ズレる。CFDでは間隙をGeometry上で正確に再現し（簡略化・省略しない）、y+<1のLow-Re壁面処理を使うことが精度の基本条件だ。
関連トピック
なった
詳しく
報告