キャビティ流れ（蓋駆動） — トラブルシューティングガイド
より充実した内容を lid-driven-cavity.html でご覧いただけます。
よくある問題と対策
キャビティ流れの計算で起きるトラブルを教えてください。
1. Ghia データと合わない
確認ポイント:
- メッシュ密度: Re = 1000 なら最低 $64 \times 64$ 以上。$32 \times 32$ だと $5\text{--}10\%$ の誤差が出る
- 空間離散化スキーム: 1次風上（First Order Upwind）だと数値拡散で渦が弱まる。最低2次精度を使う
- 参照点のずれ: Ghia のデータは $129 \times 129$ の格子点上の値。自分のメッシュとの内挿誤差に注意
- Re数の定義: $Re = UL/\nu$ であることを確認。$L$ をキャビティの一辺の長さとしているか
2. 圧力が振動する
圧力のコンター図がチェッカーボード模様になるんですが。
原因: collocated配置（速度と圧力を同じ格子点に配置）で Rhie-Chow 補間が不適切な場合に起こる。
対策:
- Rhie-Chow momentum interpolation が有効になっているか確認（OpenFOAM では標準で有効）
- staggered grid（速度と圧力を異なる格子点に配置）を使う手法もある（古典的FDMコード）
- 圧力の under-relaxation を $0.2\text{--}0.3$ に設定
3. 高Re数で収束しない
Re > 5000 では非線形性が強く、SIMPLE法の収束が遅くなる。
対策:
- 擬似時間進行法: 非定常ソルバー（PIMPLE）で十分長い時間走らせて定常解に収束させる
- continuation法: 低Re（例: Re=100）の収束解を初期条件にして、Re を段階的に上げる
- マルチグリッド: AMG前処理を圧力方程式に使用。OpenFOAMでは
GAMGソルバー
Re を段階的に上げるのは良い方法ですね。
いわゆる parameter continuation だ。Re=100 → 400 → 1000 → 3200 → 5000 → 10000 と段階的に上げる。各段階で前段の解を初期条件にする。
4. コーナー渦が解像できない
Moffatt 渦は幾何級数的に小さくなるので、均一格子では解像が困難だ。
対策: コーナーに向かって格子を幾何級数的に細かくする。OpenFOAM の blockMesh なら simpleGrading で壁面への grading ratio を $10\text{--}20$ に設定。
5. 3D計算で2Dと異なる結果が出る
3Dで計算したら2Dの結果と違うんですが。
Re > 800 程度では3D不安定性（Taylor-Goertler渦）が発現するので、3Dの結果が2Dと異なるのは物理的に正しい。
確認事項:
- スパン方向に周期境界を使っているか（端壁面がある場合は端壁効果が加わる）
- スパン方向の解像度が十分か（不安定波長の $1/10$ 以下のメッシュ幅）
- 2D計算との差が不安定性によるものか、数値誤差によるものかを区別するため、低Re（Re=100）で2Dと3Dが一致することを先に確認
段階的に検証を進める姿勢が大事なんですね。
その通りだ。cavity flow は「簡単な問題」と思われがちだが、高Re数での分岐現象、コーナー特異性、3D不安定性など、数値流体力学の本質的な問題を含んでいる。だからこそ40年以上にわたってベンチマーク問題として使われ続けているんだ。
「中央渦の中心位置がずれる」のはメッシュのせいか数値拡散か
キャビティ流れで「Ghia論文の渦中心位置と自分の計算が合わない」場合、真っ先に疑うのはメッシュ解像度ですが、実はそれ以外にも原因があります。移流項の離散化スキームが低精度（1次精度の上流差分）の場合、数値拡散で渦が「実際より下方向にぼやける」効果が出ます。少なくとも2次精度以上のスキームを使わないと中心位置はずれがちです。また、高Re数（Re=10000）ではキャビティが定常解に落ち着かず、弱い振動解になることがある。収束判定を「残差が10⁻⁶」だけで行うと「振動しながら残差が下がっている」状態を正常終了と誤認することがあります。速度の時系列も確認することが大切です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——キャビティ流れ（蓋駆動）の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告