円柱周りの流れ — トラブルシューティングガイド
より充実した内容を cylinder-flow.html でご覧いただけます。
よくある問題と対策
円柱周り流れの計算でよくハマるポイントを教えてください。
実務で頻出するトラブルを整理しよう。
1. 渦が出ない（定常解に収束してしまう）
原因と対策:
- 1次精度風上スキームを使っている → 数値拡散が渦を消す。最低2次精度に変更
- 計算領域が狭すぎる → ブロッケージ効果で実効Reが変わる。側面を $10D$ 以上離す
- 初期条件が対称 → 対称性が破れない。微小な非対称擾乱を初期条件に加える
- 定常ソルバーを使っている → 非定常ソルバー（PISO/PIMPLE）に切り替える
- 時間刻みが大きすぎる → CFL < 1 を確認。1渦放出周期あたり200ステップ以上
初期条件に擾乱を加えるって、物理的に正しいんですか？
実際の流れには必ず微小な擾乱がある。数値計算ではそれがないため、不安定でも分岐が起きないことがある。$y$方向速度に $O(10^{-3})$ 程度のランダムノイズを加えれば十分だ。
2. Strouhal数がずれる
確認ポイント:
- 時間平均の統計が十分か → 初期過渡（最低20周期）を捨ててから統計を取る
- メッシュ解像度 → 後流域のメッシュが粗いと渦が早く散逸して周波数が変わる
- ブロッケージ効果 → 計算領域が狭いとStが高くなる。Behr et al. の補正式がある
- FFTの窓関数 → 十分長い時間データでFFTを取る。ゼロパディングで周波数分解能を確保
3. $C_D$ が実験値と合わない
$C_D$ が $1.5$ とか出ちゃうんですけど。
Re=100 で $C_D > 1.4$ なら以下を確認だ。
- 2D計算で3D実験値と比較していないか: 2D計算の $C_D$ は3D実験値より $5\text{--}10\%$ 高い。これは物理的に正しい（3D効果が抗力を下げる）
- 壁面メッシュの解像度: 壁面圧力分布が正しく解像できているか確認
- 出口境界の影響: 出口が近すぎると背圧が変わり $C_D$ に影響
4. 計算が発散する
症状別の対策:
| 症状 | 原因 | 対策 |
|---|---|---|
| Courant数が爆発 | 時間刻みが大きすぎる | $\Delta t$ を半分に。adjustableRunTime を使用 |
| 圧力が振動 | チェッカーボード現象 | Rhie-Chow 補間を確認。非構造格子ならセル重心配置確認 |
| 速度が壁面近傍で発散 | $y^+$ がモデルの適用範囲外 | 壁関数なら $30 < y^+ < 300$、壁面解像なら $y^+ < 1$ に調整 |
| LES で非物理的な渦 | SGSモデル不適切 | Smagorinsky ($C_s = 0.1$)は壁面で散逸過大。WALE に変更 |
5. 3D計算でスパン方向の周期性が出ない
3D計算でスパン方向に均一な渦しか出ないんですが。
スパン方向の計算長さが不足しているか、スパン方向のメッシュが粗すぎる可能性がある。Mode A を捉えるには $L_z > 4D$、Mode B を捉えるには $\Delta z < 0.5D$ が必要だ。また、スパン方向端面は必ず周期境界条件にすること。壁面条件にすると端部効果が支配的になる。
なるほど。2Dで検証してから3Dに進むのが安全ですね。
その通り。まず2D・低Reで $C_D$, St を検証 → 3D低Reで Mode A/B を確認 → 目標Reに段階的に上げる、というのが最も確実なアプローチだ。
「渦が左右対称に出てくる」——過渡解析の初期化問題
円柱周り非定常計算で「カルマン渦が発達しない、渦が左右対称のままで流れてしまう」というトラブルは初学者が踏みやすい罠です。原因は「初期場が完全に対称で、数値的な対称性を破る擾乱がない」こと。Navier-Stokes方程式は左右対称な境界条件と初期条件なら対称解を出し続けます。渦放出を引き起こすには微小な非対称擾乱が必要です。対策は「円柱のわずかに非対称な位置に格子節点を置く」か「最初の数タイムステップだけ微小な横方向の速度擾乱を加える」かのどちらか。「物理的に起きるはずのことがCFDで起きない」のは、初期条件の対称性が原因のことが多い好例です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——円柱周りの流れの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告