スライディングメッシュ法
理論と物理
概要
Sliding Meshって非定常のターボ機械解析でよく聞きますが、何がMRFと違うんですか？
MRFはロータ位置を固定して疑似定常解を求めるが、Sliding Meshは回転域のメッシュを時間ステップごとに物理的に回転させる。非一致格子界面（Non-Conformal Interface）で情報を補間するんだ。
界面でメッシュがスライドするからSliding Meshなんですね。
そう。これにより翼-翼干渉（ウェイクの通過、ポテンシャル干渉）を物理的に捉えられる。圧力脈動や非定常翼面荷重の予測には不可欠な手法だ。
界面処理
界面の補間はどうやるんですか？
各タイムステップで回転側と静止側の界面メッシュの重なりを計算し、フラックスを保存的に補間する。CFXのTransient Rotor-Stator はGGIベースの重み付き補間を使う。FluentのSliding Meshは面対面の交差判定による補間だ。
適用範囲
どういう場合にSliding Meshが必須ですか？
| 解析目的 | Sliding Mesh必須度 |
|---|---|
| 設計点効率の予測 | 不要（MRF/Mixing Planeで十分） |
| 翼通過による圧力脈動 | 必須 |
| ウェイク-翼干渉 | 必須 |
| 非定常翼面荷重（振動評価） | 必須 |
| 騒音予測（FW-H入力） | 必須 |
| サージ/ストール | 必須（全周） |
スライディングメッシュ法の誕生——非定常ロータ-ステータ干渉CFDの夜明け
スライディングメッシュ（Sliding Mesh）法の工業的普及は1990年代後半に始まった。それ以前のターボ機械CFDは定常のフローズンローターかミキシングプレーン法に限られ、ロータとステータの非定常干渉（ウェイク干渉、ポテンシャル干渉）は捉えられなかった。転換点となったのはANSYSのFluentがバージョン5（1998年）でSliding Interface機能を商用CFDとして初めて本格実装したことだ。これにより渦のロータからステータへの輸送、BPF（Blade Passing Frequency）の脈動力計算、音響的な圧力変動の予測が可能になった。現代のターボ機械開発（ジェットエンジン圧縮機の失速解析、ポンプの水力脈動低減）では非定常スライディングメッシュが不可欠なツールとなっており、1998年の実装から25年でCFDの常識を変えた手法だ。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
タイムステップの決定
Sliding Meshの時間刻みはどう決めますか？
翼1枚通過に対して20～50タイムステップが一般的な目安だ。
例: 3000rpm、翼12枚、30steps/翼通過 → $\Delta t = 60/(3000 \times 12 \times 30) = 55.6 \mu s$
もっと細かくする必要がある場合は？
騒音予測やDES/LESではBPFの10次高調波まで解像する必要があり、100～200 steps/翼通過が要求される。
周期的定常の判定
何回転分計算すれば十分ですか？
圧力脈動のモニタリングで判断する。
1. 適当な点（カットオフ付近等）に圧力モニターを配置
2. 連続する2回転の圧力波形を重ねて比較
3. 振幅の変化が2%以内に収まれば周期的定常と判断
一般に5～15回転で安定する。初期値にMRFやFrozen Rotorの定常解を使うと3～5回転に短縮できる。
ピッチ比の処理
ロータとステータの翼枚数比が整数でない場合はどうしますか？
理想的には全周モデルだが、計算コストが膨大になる。翼枚数の最大公約数でセクターモデルを作る方法がある。例えばロータ7枚、ステータ12枚なら最大公約数1で全周が必要だが、ロータ6枚、ステータ12枚なら1/6セクターで計算できる。
CFXのTime Transformation法やFINE/TurboのNLH法は、ピッチ比が非整数でも1ピッチ計算で非定常干渉を近似できる手法だ。
スライディングメッシュ法の数値実装——AMIとフラックス修正の役割
スライディングメッシュ（Sliding Mesh）法は、回転するロータ領域と静止するステータ領域の境界を「任意メッシュインターフェース（AMI: Arbitrary Mesh Interface）」で結合し、毎タイムステップで境界のメッシュ相対位置を更新して非定常ロータ-ステータ干渉を計算する。AMIでの変数補間は双線形補間か重み付き最小二乗（WLS）で行われるが、界面のセルが非一致（Non-Conformal）の場合、フラックス保存のために「フラックス修正（Flux Correction）」が必要だ。OpenFOAMではAMI面のフラックス補正が自動化されているが、2回以上のローター回転後に質量保存誤差が蓄積する問題が知られており、定期的にpressure correctionを追加で実施する設定が推奨される。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
圧力脈動の抽出
Sliding Meshの結果から圧力脈動をどう評価しますか？
以下の手順だ。
1. 関心点（ボリュート壁面、配管接続部等）に圧力モニターを配置
2. 周期的定常に達した後の2～5回転分のデータをサンプリング
3. FFT（高速フーリエ変換）でスペクトル解析
4. BPFとその高調波のピーク振幅を評価
BPFの振幅はどのくらいが正常ですか？
機種によるが、遠心ポンプの場合BPFの圧力振幅は平均揚程の1～5%が典型的。これを超えると配管振動や構造共振のリスクがある。
翼面の非定常荷重
翼面の非定常力も評価できますか？
各翼面に作用する力のx, y, z成分を時刻歴で出力し、FFTでスペクトル解析する。翼の固有振動数とBPF高調波が一致する場合は共振（フラッタ）のリスクがある。
後処理の注意点
Sliding Meshの後処理で気をつけることは？
洗濯機ドラムのCFD——スライディングメッシュで衣類への水分配を最適化
家電製品の意外なCFD応用として洗濯機ドラムの水流解析がある。ドラム回転による洗濯液（界面活性剤含有水）の動きをVOF法＋スライディングメッシュで解析し、衣類との接触効率を向上させる設計に役立てている。ドラムに設けられたリフター（バッフル）の形状・数・配置が洗濯液の掬い上げと衣類全体への均一分配を決める。パナソニックの事例では、CFD（スライディングメッシュ+VOF+粒子追跡）で衣類汚れの分散モデルとの連成解析を実施し、洗浄効率の高い新しいリフター形状を特定することで従来比30%の節水を実現した洗濯機の開発に貢献した。製造業のCFD応用が高度化製品の差別化技術として機能している好例だ。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
ソルバー間比較
Sliding Meshの実装でソルバーに違いはありますか？
| 機能 | Ansys CFX | Ansys Fluent | STAR-CCM+ | OpenFOAM |
|---|---|---|---|---|
| 名称 | Transient Rotor-Stator | Sliding Mesh | Rigid Body Motion | cyclicAMI |
| 界面補間 | GGI (フラックス保存) | Face intersection | 内挿/外挿 | AMI (Arbitrary Mesh Interface) |
| ピッチ比 | Profile Transformation | 整数比推奨 | Profile Scaling | 整数比推奨 |
| 並列効率 | 良好 | 良好 | 良好 | やや劣る(AMI通信) |
| GPU対応 | なし | 2024R1以降対応 | 部分対応 | 一部(AmgX) |
GPU対応は進んでいますか？
Fluent 2024R1でGPUネイティブソルバーがリリースされ、Sliding Meshの非定常計算でもGPU加速が可能になった。CPUの5～10倍の高速化が報告されている。これにより従来1週間かかった計算が1日で終わるようになりつつある。
計算コストの目安
Sliding Meshの計算コストはどのくらいですか？
| ケース | セル数 | コア数 | 回転数 | 所要時間 |
|---|---|---|---|---|
| 単段ポンプ(1ピッチ) | 200万 | 32 | 10回転 | 12～24時間 |
| 単段ファン(全周) | 1000万 | 128 | 10回転 | 2～4日 |
| 多段タービン(1ピッチ) | 500万 | 64 | 15回転 | 3～5日 |
| Full-Annulus圧縮機 | 5000万 | 512 | 20回転 | 1～2週間 |
Full-Annulusは相当な計算資源が要りますね。
HPCクラスタやクラウド計算（AWS, Azure等）の利用が現実的だ。Ansys Cloud やSimcenter Cloud Engineeringで必要な時だけ大規模計算資源を確保できる。
スライディングメッシュCFDツール比較——ANSYS FluentとOpenFOAMのAMI設定の差
スライディングメッシュの実装をANSYS FluentとOpenFOAMで比較する。Fluentでは「Sliding Mesh」機能をDynamic Meshパネルで有効化し、ロータ・ステータ間のInterface Zonesを指定するだけでAMIが自動設定される。回転数の入力もCell Zone ConditionsのFrame Motionで直接設定できる簡単さが強みだ。OpenFOAMではdynamicMeshDictとfvConstraintsでAMI境界条件を手動設定し、0ディレクトリのpatchFieldにcyclicAMIを指定するなど複数の設定ファイルを整合させる必要がある。ただしOpenFOAMはParallelメッシュ分割とAMI界面の一貫性に関するバグが過去のバージョンで報告されており、v2112以降でほぼ解消された。最新バージョン確認と上流/下流のパッチペア設定の整合確認が安定運用の鍵だ。
選定で最も重要な3つの問い
- 「何を解くか」：スライディングメッシュ法に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
Sliding Mesh + DES
Sliding MeshとDESを組み合わせるのはどういう場合ですか？
非定常の渦構造と翼-翼干渉の両方を解像したい場合だ。典型的にはファン騒音予測やタービンの冷却孔-主流干渉の解析だ。
メッシュ要件は URANS と違いますか？
大幅に厳しくなる。DES/SDESでは壁から離れた領域でLESモードに切り替わるため、翼間流路の主流域でもLESに適したイソトロピックなセルが必要だ。
Sliding MeshとDESを組み合わせるのはどういう場合ですか？
非定常の渦構造と翼-翼干渉の両方を解像したい場合だ。典型的にはファン騒音予測やタービンの冷却孔-主流干渉の解析だ。
メッシュ要件は URANS と違いますか？
大幅に厳しくなる。DES/SDESでは壁から離れた領域でLESモードに切り替わるため、翼間流路の主流域でもLESに適したイソトロピックなセルが必要だ。
| 項目 | URANS | SDES/SAS |
|---|---|---|
| 翼面y+ | < 2 | < 1 |
| 翼間セルサイズ | 翼弦の2～5% | 翼弦の0.5～1% |
| セル数/ピッチ | 50～100万 | 300～1000万 |
| タイムステップ | 20～50/翼通過 | 100～500/翼通過 |
音響解析との連携
Sliding Mesh + DES + FW-Hで騒音予測する場合の手順は？
1. Sliding Mesh + SDES/LESで翼面と透過面（permeable surface）の圧力時刻歴を取得
2. 最低10BPF周期分のデータを蓄積
3. FW-Hソルバーに入力して遠方場の音圧スペクトルを算出
4. A特性重み付けで全体騒音レベル（dBA）を評価
最近の動向
Sliding Meshの最前線のトピックは？
スライディングメッシュ法の最前線——数億セルの大規模非定常LES解析
スライディングメッシュ法とLES（大渦シミュレーション）を組み合わせた大規模ターボ機械の非定常解析が、スパコンの普及で現実的な研究手段になりつつある。特にターボポンプの配管系誘起振動（FIV）やジャイロコプターロータの音響予測では、RANSでは捉えられない非定常渦構造の詳細解析が必要だ。ECNのPROEL3D（欧州水力タービン研究）では2億セル規模のカプラン水車のLES解析が実施され、ドラフトチューブの非定常圧力脈動（Vortex Rope）のRMS値を実験値±10%以内で予測することに成功した。計算には3000コアのスパコン100時間を要したが、翼疲労寿命評価への活用で得られる経済的価値はこのコストを大きく上回ると評価されている。
トラブルシューティング
典型的なトラブル
Sliding Mesh計算でよくあるトラブルは？
1. 初期発散
いきなりSliding Meshで非定常計算を始めると、初期の流れ場が未発達で発散しやすい。MRFまたはFrozen Rotorの定常解を初期値にすることで回避できる。
2. タイムステップが大きすぎる
CFLが大きいと発散しますか？
CFXの結合型ソルバーは陰的なのでCFL制限は緩いが、Sliding Mesh界面の補間精度はタイムステップに依存する。界面でのクーラン数が10を超えると補間精度が急激に劣化する。翼通過1ステップでの回転角度が5度以内になるようにする。
3. 界面のメッシュ密度不足
回転側と静止側で界面付近のメッシュ密度が大きく異なると、補間誤差が非定常計算で蓄積して振動やドリフトが発生する。界面前後のセルサイズは1:2以内に揃える。
結果検証のチェックリスト
Sliding Meshの結果をどう検証すればいいですか？
時間平均効率がMRFと大きくずれたら何が問題ですか？
2ポイント以上ずれる場合は、タイムステップ不足、メッシュ品質不良、または初期過渡の切り捨て不十分が考えられる。時間平均の開始点を十分に遅らせて（5回転以降から平均開始等）再評価すること。
スライディングメッシュCFDで非定常力が振動する——タイムステップとAMI補間の関係
スライディングメッシュCFDで「ロータへの力がタイムステップごとに大きく振動し、収束判定ができない」問題は、タイムステップとAMI補間の組み合わせに起因することが多い。ロータが1タイムステップで大きく回転すると、AMI界面のセル相対位置が大きく変わり補間エラーが発生する。ガイドライン：翼通過角度あたりのタイムステップ数を最低20ステップ以上確保すること（翼枚数Zと回転数Nから Dt_max = 1/(20*N*Z) を計算）。またロータ-ステータ力を「翼通過周期（BPP: Blade Passing Period）」で平均化することで短周期の数値振動を除去する。さらにAMI補間のorder（1次 vs 2次）を上げると振動が低減するが計算コストが増す——収束性の確認後に選択することが実務的な手順だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——スライディングメッシュ法の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告