混合面法
混合面法の理論基礎
概要
Mixing Planeって多段ターボ機械の定常計算で必ず出てきますよね。正確にはどういう処理ですか？
動翼と静翼の界面で流量を周方向に質量平均し、下流への境界条件として渡す手法だ。これにより異なるピッチ数の翼列を定常的に接続できる。
$\phi$ は全圧、全温、流れ角、乱流量などの保存変数だ。
質量平均するとウェイク（後流）の情報は消えますよね？
その通り。Mixing Planeの最大の限界は、上流翼列の翼後流や二次流れの周方向変動が下流に伝わらないことだ。それでもスパン方向の変動は保持されるから、段性能の予測には十分な精度が得られる。
保存性の確保
Mixing Planeで質量やエネルギーの保存は担保されますか？
質量流量、運動量、全エンタルピーの保存は数値的に強制される。CFXではフラックス保存型の実装がされていて、界面の上流と下流で質量流量が完全に一致する。ただし混合エントロピーの増加（混合損失）は数値的に避けられない。
混合損失ってどのくらいですか？
段効率に0.1～0.5ポイント程度の影響がある。実機でも動翼-静翼間で混合は起きるから、ある程度は物理的に正当化されるが、Mixing Planeの位置や手法によって変わるため注意が必要だ。
ミキシングプレーン法の誕生——Denton & Dawes（1988年）とターボ定常解析の標準化
ミキシングプレーン（Mixing Plane）法は1988年にJohn Denton（ケンブリッジ大）とW.N. Dawesが発表した論文で確立された。それ以前の多段ターボ機械解析は各段を独立に解析して出力・入力を繋ぐしかなく、段間の流れ場の整合的な予測が困難だった。DentonとDawesは「ステージ境界で周方向流れ変数を平均化して交換する」というMixing Planeの考え方を提案し、複数段の連続した定常解析を可能にした。この概念は3〜5年後には主要なCFDコードに実装され、1990年代のガスタービン設計の標準ツールとなった。Dentonは後にターボ機械CFDの開発貢献を称えられてAIME Melchior Jackman賞を受賞しており、ミキシングプレーン法はターボCFD史上最も影響力のある発明の一つとして認識されている。
混合面法の数値計算手法
実装のバリエーション
Mixing Planeの実装はソルバーごとに違いますか？
大きく2つのアプローチがある。
| 手法 | 説明 | ソルバー |
|---|---|---|
| Band-averaged | スパン方向にバンド分割し各バンドで周方向平均 | CFX (Stage Interface) |
| Profile-transfer | スパン方向のプロファイルを連続関数として転写 | FINE/Turbo, STAR-CCM+ |
Band-averagedのバンド数は多いほうがいいですか？
バンドが少なすぎるとスパン方向の変動が平滑化される。CFXのデフォルトは概ね十分だが、端壁付近の急な変動を捉えたい場合はバンド数を増やすか、ユーザー定義のバンド分布を設定する。
多段計算のMixing Plane配置
多段軸流タービンでMixing Planeはどこに置きますか？
各動翼-静翼間に1つずつ置くのが基本だ。Mixing Planeの位置は翼後縁と次の翼前縁の中間付近が推奨される。前縁に近すぎるとポテンシャル干渉が大きくなり、後縁に近すぎるとウェイクの混合が不十分になる。
NUMECA FINE/TurboではどうMixing Planeを設定しますか？
AutoGrid5で翼列境界を定義すると、Row Interface として自動的にMixing Planeが配置される。Non-Reflectingオプションを選ぶと、界面での圧力波の反射を抑制してサージ近傍の計算安定性が向上する。
Non-Reflecting Mixing Plane
Non-Reflectingってどういう意味ですか？
通常のMixing Planeでは界面での圧力変動が反射して数値的な振動を起こすことがある。Non-Reflecting処理はGiles特性条件を適用して、界面での波動を透過させる手法だ。高負荷の圧縮機段やサージ近傍で特に効果がある。
ミキシングプレーン法の数値的実装——周方向平均の取り方とラジアル分布の保存
ミキシングプレーン（Mixing Plane）法では、ロータ-ステータ境界で周方向平均化した流れ変数を交換することで、定常解析でロータとステータを結合する。数値実装上の重要点は「どの変数を平均するか」だ。圧力・速度の単純算術平均では、高エンタルピー域でエネルギー保存が崩れる——代わりに質量流量加重平均（Mass-Flux-Weighted Average）を使うことが精度上必要だ。またラジアル（径方向）の分布は保持し、周方向のみ平均するラジアル分布保存型の設定が高精度の標準だ。CFXでは界面の適切な平均化が自動実装されているが、OpenFOAMのMRFfvPatchFieldは設定によって単純平均になるケースがあり、実装仕様の確認が必要だ。
混合面法の実務適用
多段軸流タービンの解析
ガスタービンの多段タービンをCFDで解析する際のMixing Plane活用法を教えてください。
HP（高圧）タービンの1～2段＋LP（低圧）タービンの3～5段で合計8～12翼列になる。各翼列間にMixing Planeを配置して全段定常計算を行う。
セル数はどのくらいになりますか？
1翼列50～100万セルとして、10翼列で500～1000万セル。128コアで12～24時間程度だ。Mixing Planeの恩恵で1ピッチ計算ができるから、全周比で翼枚数分の1に圧縮される。
結果評価のポイント
多段タービンのCFD結果で何を確認すべきですか？
| 評価項目 | 確認方法 | 基準 |
|---|---|---|
| 段間の質量流量保存 | Mixing Plane前後の流量差 | 0.01%以内 |
| 各段の圧力比 | Mixing Plane面での質量平均全圧比 | 1D設計値と±2%以内 |
| スパン方向の効率分布 | スパン断面での断熱効率 | ハブ/チップで低下が自然 |
| 翼面マッハ数 | 翼面の等マッハ数コンター | 衝撃波位置の確認 |
Mixing Planeの混合損失は結果にどう影響しますか？
段数が多いほどMixing Plane面の数が増え、累積する混合損失も大きくなる。10段で累積0.5～1ポイントの効率過小評価が生じうる。非定常計算と比較して補正係数を持っておくのが実務的だ。
クーリング流のモデル化
タービン翼の冷却孔からの噴流はどう扱いますか？
フィルム冷却の全孔をメッシュ化するのは現実的でない。以下のアプローチが使われる。
- Source Term Model: 冷却流量を翼面近くのセルに質量・運動量・エネルギーのソースとして注入
- Discrete Hole Model: 各冷却孔を小さなInletとしてモデル化（CFXのInlet BC）
- Conjugate Heat Transfer: 翼内部の冷却通路もメッシュ化して流体-固体連成
多段軸流コンプレッサーのCFD——ミキシングプレーンで6段の定常特性を計算
航空エンジンや産業用ガスタービンの多段軸流コンプレッサーは5〜15段以上の翼列が直列に並ぶ複雑なターボ機械だ。全段を非定常スライディングメッシュで計算するのは研究用途でも困難で、実際の設計・最適化ではミキシングプレーン法による定常多段解析が主流だ。6段コンプレッサーの1/5扇形周期性モデルでミキシングプレーンによる定常RANS計算を行うと、全ステージのインタラクションを含む圧力比・効率を実験±3%以内で予測した実績がGE Aerospaceの内部報告で示されている。特に段間の再循環域（Surge寸前の条件）では定常仮定が破綻するため、設計マージン（Stall Margin）の評価には非定常SlidingMesh解析を組み合わせることが安全な設計フローだ。
混合面法のソフトウェア比較
CFX vs FINE/Turbo vs STAR-CCM+
Mixing Planeの機能でソルバーを選ぶとしたら、どこが違いますか？
NUMECAが多段に一番強い印象ですね。
FINE/Turboは元々多段ターボ機械のために開発されたソルバーだから、多段のセットアップとMixing Planeの扱いは最も洗練されている。一方CFXはユーザーベースが最大で、サポート情報が豊富だ。
OpenFOAMでのMixing Plane
OpenFOAMにMixing Planeはありますか？
標準にはない。cyclicAMIで回転界面を設定できるが、これはSliding Meshに相当する。Mixing Planeを実装するにはカスタムBCの開発が必要で、研究グループがいくつか公開しているが品質はまちまちだ。
ワークフローの自動化
多段の特性マップ作成は自動化できますか？
CFXではCFX-Pre のJournal（Pythonスクリプト）で背圧を変えた連続計算を自動実行できる。FINE/Turboには「AutoRun」機能があり、特性曲線の自動取得がGUIから設定可能だ。STAR-CCM+ではJavaマクロで同様の自動化が可能だ。
ターボCFDツール比較——NUMECA FINE/TurboとANSYS CFXのミキシングプレーン実装
ターボ機械向けCFDの2大商用ツール、NUMECA FINE/TurboとANSYS CFXを比較する。FINE/Turboはターボ専用で開発されており、ミキシングプレーン・スライディングメッシュ・非線形調和法（NLH）をシームレスに切り替えられる設計が強みだ。翼型の形状最適化機能（AutoDesign）と連携したCFD最適化ループが容易で、航空エンジン・コンプレッサの研究開発機関での採用率が高い。ANSYS CFXはターボ機械向けインターフェースの成熟度が高く、産業用ターボポンプ・風力発電機など幅広い応用実績がある。FINE/TurboがCFXよりやや低い傾向があるが、Ansys全ソリューションとの連携（構造・電磁場・システム）が必要な場合はANSYSエコシステムが優位になる。
混合面法の先端研究
Mixing Planeの限界を超える手法
Mixing Planeの限界を克服する手法にはどんなものがありますか？
Mixing Planeの主な限界は「非定常干渉を捕えない」ことだ。これを解決する手法を精度・コスト順に並べよう。
| 手法 | 非定常情報 | コスト（vs Mixing Plane） | 適用ソルバー |
|---|---|---|---|
| Mixing Plane | なし（周方向平均） | 1x | 全主要ソルバー |
| NLH (Non-Linear Harmonic) | 基本周波数＋高調波 | 3～5x | FINE/Turbo |
| Time Transformation | 位相ずれの周期解 | 5～10x | CFX |
| Sliding Mesh (URANS) | フル非定常 | 20～100x | 全主要ソルバー |
NLHってどういう手法ですか？
非定常変動を時間平均成分＋基本周波数成分＋高調波に分解して周波数領域で解く。時間方向のステップが不要だから、Sliding Meshの1/5～1/20のコストで翼列干渉を捕えられる。
CFXのTime Transformationは？
Pitch Ratioが整数比でない翼列の非定常干渉を、位相変換（time-shift）で周期的に処理する手法だ。1ピッチ計算で全周の非定常干渉をある程度再現できる。
どの手法を選ぶべきか
実務での使い分けの指針はありますか？
非線形調和解析法（NLH）——ミキシングプレーンを超えるターボ非定常解析
ミキシングプレーン法はロータ-ステータ間の周方向流れを平均化して定常解を得るが、非定常干渉（ウェイク、衝撃波）の情報が失われる。この限界を超えるために「非線形調和解析法（NLH: Non-Linear Harmonic Method）」が開発された。NLHは非定常流れ場をFourier級数で展開し、各調波成分の結合方程式を定常的に解く手法で、完全非定常解析の1/10〜1/5の計算コストでロータ-ステータ干渉を再現できる。NUMECA社のFINE/Turboに実装されたNLHは、ガスタービンの多段解析でMixing Planeより20〜30%精度が高い効率予測を実現しており、高精度ターボCFDの次世代標準として注目されている。
混合面法のトラブル対応
Mixing Plane界面のトラブル
Mixing Planeで計算が不安定になる原因は何ですか？
代表的な問題を挙げよう。
1. 界面での逆流
低流量や部分負荷でMixing Plane面の一部に逆流が発生すると、境界条件の設定が破綻して発散する。Non-Reflecting Mixing Planeを使うか、界面位置を翼から離すことで改善する場合がある。
2. スパン方向メッシュの不整合
上流と下流でスパンのメッシュ分布が違うと問題ですか？
Mixing Planeでのスパン方向補間精度が落ちる。特にハブ・シュラウド付近のプリズム層がMixing Plane面に集中している場合、補間が不安定になる。界面から少し離した位置でスパンメッシュを統一しておくのが対策だ。
3. 翼列間のBleeding（抽気/注入）
段間で冷却空気の注入やブリーディングがある場合、Mixing Plane面を通過する質量流量が上下流で異なる。この不整合を境界条件で正しく処理しないと質量保存が崩れる。
収束改善テクニック
Mixing Plane計算の収束を改善するコツは？
| テクニック | 詳細 |
|---|---|
| 初期値の工夫 | 1D設計の速度・圧力分布を初期値にする |
| 背圧ランピング | 出口背圧を徐々に目標値まで上げる |
| 擬似タイムステップ | Auto Timescale Factorを0.5～0.75に下げる |
| 乱流モデルの2段階 | まずk-epsilonで粗く収束→SSTに切り替え |
| MPI並列のドメイン分割 | Mixing Plane界面をまたぐパーティションを避ける |
k-epsilonで初期収束させてからSSTに切り替えるのは面白いですね。
k-epsilonはSSTより収束性が良い傾向があるから、まず全体の流れ場を大まかに確立させる。その後SSTに切り替えることで、壁近傍の精度を上げつつ安定した収束が得られる。
ミキシングプレーンCFDで効率が急変する——界面流れの異常診断
ミキシングプレーンCFDで「ロータ出口とステータ入口の全圧差が数%ある」「全圧効率の計算値が95%超という非現実的な値になる」という問題は、界面の変数定義の混乱に起因することが多い。特によく起きるのは静圧（Static Pressure）と全圧（Total Pressure）の混同——インターフェースで静圧を渡すべき所に全圧を設定するか、またはその逆の設定ミスだ。流れ変数のどちらを界面で渡すか受け取るかの設定はツールごとに異なり、CFXではInlet/Outlet BCのType設定で区別する必要がある。診断は界面の両側で積分した総エンタルピーを比較し、保存されているかを確認する基本的なエネルギーバランスチェックから始める。
関連トピック
なった
詳しく
報告