凍結ロータ法
理論と物理
概要
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年後の現代でも設計探索フェーズの標準手法として生き続けている。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
位置依存性の問題
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）だ。また界面の流れ変数の補間精度が精度に影響するため、スライディングメッシュ（より高精度）との比較検証を重要な部位では実施することが推奨される。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
遠心ポンプのモデル構成
遠心ポンプの全体モデルはどう組むんですか？
典型的な構成は以下だ。
- 吸込管: 静止域、非回転
- インペラ: 回転域、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週間に短縮した実績がある。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
ソルバー別の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を手動で設定するため、設定ファイルへの理解が必要でターボ専用ウィザードはない。ただしスクリプト化されれば多数のパラメータ変更も自動化でき、最適化ループには向く。
選定で最も重要な3つの問い
- 「何を解くか」：凍結ロータ法に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
水力タービンでの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を超える場合は補間誤差の主因になるため、必ず品質チェックを実施する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——凍結ロータ法の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告