球周りの流れ — トラブルシューティングガイド
より充実した内容を flow-around-sphere.html でご覧いただけます。
よくあるトラブル
球周りの流れの計算でよくある失敗を教えてください。
1. 軸対称計算で後流が非対称にならない
原因: Re > 210 でも軸対称解が出てしまう場合、メッシュの対称性が厳密すぎるか、初期条件が完全に軸対称だ。
対策: フル3D計算に切り替え、微小な非対称擾乱を初期条件に加える。もしくはwedge（2D軸対称）メッシュで計算している場合は、フル360度メッシュに変更する。
2. Stokes 解との不一致（低Re数）
Re = 0.1 で計算したのに $C_D = 24/Re$ と合わないんですが。
確認ポイント:
- 計算領域のサイズ: 低Re数では Stokes 流れの影響範囲が広い。外部境界が球から $100D$ 以上離れていないと、境界の影響で $C_D$ がずれる
- 入口・側面条件: 一様流入口が近すぎると球の上流にまで影響が及ぶ（Stokes流の影響は $1/r$ で減衰するため極めて遠くまで到達する）
- Oseen 補正: Re = 0.1 でもStokes解からの偏差は数%ある。$C_D = (24/Re)(1 + 3Re/16) = 240.45$ であり、$C_D = 240$ とは $0.2\%$ 異なる
3. 抗力危機が捉えられない
RANS の SST $k$-$\omega$ では、乱流モデルの仮定により境界層は最初から乱流として計算される。そのため、層流剥離→遷移→乱流再付着というプロセスが再現されず、抗力危機が起こらない。
対策: $\gamma$-$Re_\theta$ 遷移モデルを有効にする。Fluent なら Models → Viscous → Transition SST (4-equation) を選択。
4. 後流の渦構造が崩れる
ヘアピン渦がうまく出ないんです。
確認ポイント:
| 症状 | 原因 | 対策 |
|---|---|---|
| 渦が下流で消える | 数値拡散 | 2次精度以上のスキーム使用。中心差分の混合率を上げる |
| 渦の形状が不自然 | メッシュ異方性 | 後流域のアスペクト比を3以下に。等方的なメッシュが理想 |
| 渦放出が不規則 | 統計収集時間不足 | 少なくとも50周期以上の統計を取る |
| 3Dモードが出ない | メッシュ解像度不足 | 方位角方向に最低64分割（5.6度間隔）以上 |
5. 計算が発散する（高Re数LES）
高Re数の LES で中心差分スキームを使うと、格子Reynolds数が大きい領域で振動・発散することがある。
対策:
- 混合スキーム（中心差分 + 少量の風上）を使用。Fluent の Bounded Central Differencing がこの方式
- SGS モデルの拡散を利用。Dynamic Smagorinsky は局所的に拡散を適応的に調整する
- メッシュを細かくして格子Reynolds数を下げる（根本的だが高コスト）
結局、問題が起きたらまずメッシュと境界条件を疑えということですね。
その通り。そして低Reの検証済みケースから始めて、段階的にRe数を上げていくのが最も確実だ。Re=100 で $C_D = 1.09$、$L_w/D = 1.7$ が出ることを確認してから進むべきだ。
「CDが実験と3割違う」——境界層遷移モデルを疑う前に確認すること
球の抗力係数CDをCFDで計算すると実験値より30%ほど高く出ることがあります。よくある原因は「ドメインが狭すぎる閉塞効果」と「表面メッシュが粗くて曲面の解像が不十分」の2つです。ドメイン閉塞率（球断面積/ドメイン断面積）が5%を超えると壁の影響で流速が加速して抗力が増加します。これは実験でも同様で、風洞試験では閉塞補正（blockage correction）が標準的に行われます。CFDではドメイン断面積を十分大きくして閉塞率を1%以下にするのが基本。それでも合わない場合に初めて乱流モデルや遷移モデルの精度を疑う順番で調査することが重要です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——球周りの流れの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告