凍結ロータ法
凍結ロータ法の理論基礎
概要
Frozen Rotor法って名前がカッコいいですけど、何をしてるんですか？
回転体と静止体の相対位置を固定したまま定常解析する手法だ。動翼の座標を回転系に置き、静翼の座標を静止系に置くが、界面で周方向の平均をとらずに「そのまま」情報を渡す。
Mixing Planeとの違いは？
| 手法 | 界面処理 | 翼枚数制約 | 計算コスト | 精度 |
|---|---|---|---|---|
| Frozen Rotor | 位置固定、直接補間 | 1ピッチが異なるとピッチ比補正必要 | 低い | 位置依存性あり |
| Mixing Plane | 周方向質量平均 | ピッチ比任意 | 低い | 周方向変動を平滑化 |
| Sliding Mesh | 時間進行で回転 | ピッチ比は整数比が望ましい | 高い | 最も正確 |
いつ使うか
Frozen Rotorはどういう場面で使うんですか？
以下のケースに適している。
- 初期設計段階での高速スクリーニング: 多数の形状案を比較する際にMixing Planeより速い
- ボリュート付き遠心機: ボリュートは非軸対称だからMixing Planeが適用しにくい。Frozen Rotorで翼-ボリュート干渉を概算する
- 水力タービン: ドラフトチューブとの干渉評価
ボリュートのようにMixing Planeが使えない場合に重宝するわけですね。
そう。ただし結果は翼とボリュートの相対位置に依存する。正確な評価には複数の相対位置で計算して平均するか、最終的にSliding Meshに進むべきだ。
CFXでのFrozen Rotor設定
CFXでFrozen Rotorを設定する手順を教えてください。
ドメインインターフェースのタイプを「Frozen Rotor」に設定するだけだ。ピッチ比が1:1でない場合は「Pitch Change」オプションで自動スケーリングされる。注意点としてGGI面のメッシュピッチは概ね合わせておくこと。大きくずれると補間誤差が増える。
ターボ機械CFDの歴史——1970年代のスリム・スロット法から3D非粘性解析へ
ターボ機械の内部流れをCFDで解析する試みは1970年代に始まった。当初は計算機の制限から翼列を2D解析が主流で、3Dのブレード形状は考慮できなかった。転換点は1970年代後半のANSYS Fluent前身コードの開発と、Denton(1982)がIBMメインフレームで実現した3D非粘性（オイラー方程式）ターボ機械解析だ。さらに1990年代にHarvey & Denton、Arnone等がロータ-ステータの非定常干渉を含む粘性解析を実現し、現代のターボCFD（RANS+スライディングメッシュ）の原型が確立した。フローズンローター法はこの歴史の中で「定常近似の合理的な最初の選択」として生まれ、50年後の現代でも設計探索フェーズの標準手法として生き続けている。
凍結ロータ法の数値計算手法
位置依存性の問題
Frozen Rotorの結果が翼の位置で変わるって本当ですか？
本当だ。例えば遠心ポンプで翼がボリュートの舌部（カットオフ）の正面にある場合とずれている場合で、揚程が5～10%変わることがある。これはFrozen Rotorの本質的な限界だ。
じゃあ信頼できないんですか？
単一位置の結果をそのまま性能値として報告するのは危険だ。翼ピッチの1/3～1/2間隔で3～5位相を計算して平均するのが推奨される。
ピッチ比補正
ピッチ比って何ですか？
ロータとステータの1ピッチの角度の比だ。例えばロータ7枚（ピッチ51.4度）、ステータ12枚（ピッチ30度）ならピッチ比は51.4/30=1.71。Frozen Rotorでは界面の両側でこのピッチ差を何らかの形で処理する必要がある。
CFXでは界面にPitch Ratioを設定すると、周方向にスケーリング補間が行われる。ただしピッチ比が2を超えると精度が大きく劣化する。
Mixing PlaneとFrozen Rotorの使い分け
実務ではどう使い分けるんですか？
| 状況 | 推奨手法 |
|---|---|
| 軸対称ディフューザ/ボリュートなし | Mixing Plane |
| ボリュート付き遠心ポンプ | Frozen Rotor（複数位相）→ Sliding Mesh |
| 多段軸流 | Mixing Plane |
| 予備設計のパラメトリック | Frozen Rotor（高速） |
| 圧力脈動・騒音評価 | Sliding Mesh（必須） |
最終的にはSliding Meshで検証するのが安全ですね。
その通り。Frozen Rotorは「速い概算」、Mixing Planeは「安定した定常近似」、Sliding Meshは「物理的に正しい非定常解析」という位置づけだ。
フローズンローター法の数値設定——インターフェース処理と収束のコツ
フローズンローター（Frozen Rotor）法はロータ-ステータ間の周方向に一様な仮定のもとでMRFを固定した定常解析法だ。実装では「ロータ-ステータインターフェース」の処理が精度を決める。このインターフェースでは座標系変換（回転座標系から静止座標系へ）を行うが、周方向の流れが不均一な場合（ウェイク干渉が強い）は仮定が成立せず、異なる周方向位置で解を評価すると結果が変化する「周方向依存性」が生じる。対策は複数の周方向位置でフローズンローター計算を実施し、周方向平均を取ること（Pitch Average）だ。また界面の流れ変数の補間精度が精度に影響するため、スライディングメッシュ（より高精度）との比較検証を重要な部位では実施することが推奨される。
凍結ロータ法の実務適用
遠心ポンプのモデル構成
遠心ポンプの全体モデルはどう組むんですか？
典型的な構成は以下だ。
- 吸込管: 静止域、非回転
- インペラ: 回転域、MRFまたはSliding Mesh
- ボリュート: 静止域、非軸対称
- インペラ-ボリュート界面: Frozen Rotor（定常）or Sliding Mesh（非定常）
ボリュートが非軸対称だからMixing Planeは適用できない。これがFrozen Rotorが重宝される理由だ。
ボリュートのメッシュはどうやって作りますか？
ボリュートは断面形状がスパイラル状に変化するから、TurboGridでは作れない。Ansys MeshingやFluent Meshingの非構造メッシュ、またはSTAR-CCM+の自動メッシュを使う。断面形状をスイープ方向に並べてヘキサドミナントメッシュにするのが品質面で有利だ。
ボリュート舌部の処理
カットオフ（舌部）付近は何が難しいんですか？
舌部はインペラ出口流と再循環流が衝突する領域で、圧力勾配が急峻だ。メッシュを特に細かくする必要がある。また、Frozen Rotorでは翼と舌部の相対位置で流れ場が大きく変わるため、圧力脈動評価にはSliding Meshが不可欠だ。
性能マップの作成
Frozen Rotorで遠心ポンプのH-Q曲線は作れますか？
作れるが、各流量点で複数位相の平均をとることを推奨する。手順は以下だ。
1. 設計流量でFrozen Rotor計算を収束
2. 翼を10度間隔で3～5位相回転させて再計算
3. 各位相の揚程を算術平均
4. 流量を変えて2-3を繰り返し
結構手間ですね…
だから実務ではFrozen Rotorで設計点付近の1～2点を素早く評価し、最終的な性能マップはSliding Meshの非定常計算で確定させるケースが多い。
水力タービン（カプラン型）のフローズンローター解析——ガイドベーン角最適化
可動ガイドベーン（入口案内羽根）を持つカプラン水車のCFD設計では、フローズンローター法が多数のガイドベーン開度（典型的に10段階）と回転数の組み合わせで性能マップを高速生成するのに活用される。全流量・水頭の運転条件範囲でη-Q曲線を計算し、最高効率点（BEP: Best Efficiency Point）付近の設計指針を確立する作業は設計初期に必須だ。各ガイドベーン角度でのフローズンローター計算（定常RANS）は1ケース当たり数時間で完了し、10ケースの性能マップが1〜2日で完成する——スライディングメッシュなら同じ計算に数週間かかる。スウェーデンの水力発電設備メーカーが公開した事例では、フローズンローター+最適化で設計期間を3ヶ月から3週間に短縮した実績がある。
凍結ロータ法のソフトウェア比較
ソルバー別のFrozen Rotor実装
ソルバーごとにFrozen Rotorの実装は違いますか？
名称と設定方法が異なる。
| ソルバー | 名称/設定 | ピッチ比補正 | 備考 |
|---|---|---|---|
| Ansys CFX | Domain Interface → Frozen Rotor | Pitch Change自動 | 最も使い慣らされた実装 |
| Ansys Fluent | MRF + Interface | なし（1:1が前提） | Cell Zone Motionで回転を設定 |
| STAR-CCM+ | Frozen Rotor Interface | Profile Scaling | Rigid Body Motion指定 |
| OpenFOAM | MRF (steadyState) | なし | cyclicAMIで界面処理 |
Fluentにはピッチ比補正がないんですか？
FluentのMRFでは界面の両側が同じ角度範囲を持つことが前提だ。ピッチ比が異なる場合は全周モデルか、整数比になるセクターモデルを作る必要がある。この点はCFXのほうが柔軟だ。
注意すべき制限事項
Frozen Rotorの落とし穴を教えてください。
1. ウェイクの凍結: 翼後流が固定された位置に留まるため、下流への影響が過大/過小になりうる
2. 非軸対称入口条件: ボリュート下流のディフューザで非軸対称な入口流れが発生するが、Frozen Rotorでは位置依存
3. トルク変動の見落とし: 実機では翼とカットオフの相対位置で変動するトルクが、Frozen Rotorでは1つの値に固定される
4. 全周モデルの必要性: 非軸対称システムではセクターモデルが使えない場合が多い
つまりFrozen Rotorは「定常の近似」であることを常に意識しないといけませんね。
そう。Frozen Rotorの結果を最終設計根拠にするのは危険だ。必ずSliding Meshで検証するか、少なくとも複数位相平均で不確実性を見積もること。
ターボ機械向けCFDツール——NUMECA FINE/Turbo vs ANSYS CFX vs OpenFOAMの設定難易度
ターボ機械の定常解析（フローズンローター・MRF）で主要3ツールを設定の難易度で比較する。NUMECA FINE/Turboはターボ専用のGUIで、回転方向・回転数・ロータ-ステータ界面の設定がウィザード形式で完結し、初期設定が最も容易だ。翼列解析に特化したメッシュ品質指標とエラーチェックが充実しており、初期収束エラーを自動診断する機能もある。ANSYS CFXはCFX-Preのウィザードでほぼ同等の便利さを持ち、Ansysライセンスを持つ企業への普及率が高い。OpenFOAMはrotatingWallVelocity境界条件とMRFfvPatchFieldを手動で設定するため、設定ファイルへの理解が必要でターボ専用ウィザードはない。ただしスクリプト化されれば多数のパラメータ変更も自動化でき、最適化ループには向く。
凍結ロータ法の先端研究
水力タービンでのFrozen Rotor
水力タービンのCFDでもFrozen Rotorが使われるんですか？
フランシス水車やカプラン水車の初期設計でよく使われる。ステイベーン→ガイドベーン→ランナー→ドラフトチューブの4ドメイン構成で、ランナー前後にFrozen RotorまたはMixing Planeを設定する。
水力タービンだとドラフトチューブが問題になるって聞きました。
ドラフトチューブ内にはランナー出口の旋回残りによる渦ロープ（vortex rope）が発生する。これは非定常現象だからFrozen Rotorでは捉えられない。部分負荷での渦ロープによる圧力脈動評価にはSliding Meshが必須だ。
IEC 62006ガイドライン
水力タービンのCFDには規格がありますか？
水車のヒル図作成
ヒル図って何ですか？
単位流量 $Q_{11}$ vs 単位回転数 $N_{11}$ の平面上に効率をコンター表示したものだ。
ガイドベーン開度、回転数、流量を変えて100点以上の運転点をCFDで計算し、効率マップを生成する。Frozen Rotorの高速性がここで活きるんだ。
100点以上！ それはFrozen Rotorじゃないと回せませんね。
そう。1点あたり数時間で済むFrozen Rotorなら並列で数日で完了するが、Sliding Meshだと1点に1日かかるから100点は非現実的だ。
フローズンローター法の精度限界——非線形干渉が強いケースの見極め方
フローズンローター（Frozen Rotor）法は定常MRF解析で計算が速い便利な手法だが、適用できる物理条件には限界がある。ロータ-ステータ間の干渉が強い条件——①ロータ-ステータ間距離が短い（s/R<0.05）、②翼の負荷が高く強いウェイクが発生している（Zweiffel数>0.8）、③再循環域がインターフェースを跨ぐ低流量・高負荷運転——では周方向依存性が大きくなり、解が「インターフェース相対位置」に敏感になる。定量的な判断基準は「周方向3〜5点でフローズンローター計算を行い、性能（効率・全圧）の変動が±2%以上なら非定常スライディングメッシュが必要」という実務ガイドラインがある。まずフローズンローターで感度チェックを実施し、必要なら非定常に切り替えるという2段階アプローチが合理的な選択だ。
凍結ロータ法のトラブル対応
Frozen Rotor界面の収束問題
Frozen Rotor計算で界面付近の残差が下がりません。何が原因ですか？
いくつかの原因がある。
1. ピッチ比が大きすぎる
ピッチ比が2以上だと補間による数値的なオーバーシュートが発生しやすい。翼枚数の最大公約数でセクターモデルを組めないか検討しよう。
2. 界面前後の大きな速度差
回転域と静止域で速度が大きく違うと問題になりますか？
当然だ。回転系の相対速度から静止系の絶対速度に変換する際に大きな不連続が生じると収束が悪化する。初期条件として前の運転点の解を使うか、回転数を徐々に上げるランピングが有効だ。
3. GGI面のメッシュ不整合
界面両側のメッシュサイズが大きく異なると補間精度が劣化する。回転側と静止側で界面近傍のセルサイズを概ね揃える（比1:2以内）こと。
回転数ランピング
回転数ランピングって具体的にはどうするんですか？
CFXではExpert Parameterで回転速度をステップごとに増加させる設定が可能だ。例えば目標3000rpmに対して、最初300rpmで100イテレーション、次に1000rpm、最後に3000rpmと段階的に上げる。FluentではUDF内でZone Motionの回転速度を時間関数として記述する。
Frozen Rotor → Sliding Meshへの移行
Frozen Rotorの結果をSliding Meshの初期値に使えますか？
使える。CFXではFrozen Rotorの結果ファイル（.res）をSliding Meshの初期条件としてリスタートできる。これで非定常計算の立ち上がりを大幅に短縮できるんだ。
それは便利ですね。
設計のワークフローとしては、Frozen Rotorで概算→Sliding Meshで精緻化という2段階アプローチが最も効率的だ。
フローズンローターCFDで圧力比が安定しない——界面インターポレーションのミスマッチ
フローズンローター解析で「ロータ-ステータ界面を挟んで圧力に急なジャンプが生じる」問題は、インターフェースの流量保存精度に起因する。特に非一致界面（Non-Conformal Interface）では、ロータ側とステータ側のセル数・サイズが異なり、面積加重平均によるフラックス交換で質量が保存されないケースがある。診断手順：①界面の両側で積分した質量流量を比較し、1%以上の差があれば設定ミス。②インターフェースに混在した設定（ミキシングプレーンとフローズンローターの同時指定）がないか確認。③ロータ側とステータ側のメッシュが界面付近で大きさが揃っていない場合はリファイン。界面近傍のセルアスペクト比が10を超える場合は補間誤差の主因になるため、必ず品質チェックを実施する。
関連トピック
なった
詳しく
報告