ウィンデージ損失
理論と物理
概要
ウィンデージ損失って何ですか？ あまり聞き慣れない用語ですが。
回転体と静止壁面の間の空気（または流体）の摩擦による動力損失だ。高速回転機械では無視できない大きさになる。タービン発電機、フライホイール、高速モータなどで特に重要だ。
回転ディスクの摩擦損失
具体的にどのくらいの損失ですか？
閉空間内の回転ディスクの摩擦損失は以下で表される。
$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年後のジェットエンジン冷却設計に直接応用されている点は特筆に値する。
各項の物理的意味
- 時間項 $\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モデルの構成
ウィンデージ損失を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で精密確認する手順が標準的だ。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
ロータ系の温度上昇
ウィンデージ損失で温度が上がるんですか？
損失は全て熱になる。密封空間では冷却が不十分だと温度が上昇し続ける。蒸気タービンの車室内やモータの空隙で問題になる。
冷却流量 $\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適用分野となっている。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ 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が選ばれる傾向が続いている。
選定で最も重要な3つの問い
- 「何を解くか」：ウィンデージ損失に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
テイラークエット流れ
回転ディスク間の流れって渦ができるんですか？
同心円筒間の流れ（テイラークエット流れ）ではテイラー渦が発生する。回転レイノルズ数が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%以内に収まれば収束だ。
よくある間違い
ウィンデージ解析でやりがちなミスは？
ウィンデージ解析の収束はどう判定しますか？
ロータ表面のトルクをモニターする。定常計算ならトルクの変動が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壁面処理を使うことが精度の基本条件だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ウィンデージ損失の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告