円柱周りの流れ — 商用ツール比較と選定ガイド

各ツールの特徴を円柱流れの観点で比較しよう。

pimpleFoam を使う場合のディレクトリ構成と主要設定はこうだ。

```

constant/

transportProperties: nu = 1e-4 (Re=100ならU=1, D=0.01)

turbulenceProperties: simulationType laminar (Re=100の場合)

system/

controlDict: deltaT=0.005, endTime=300

fvSchemes: ddt(backward), div(phi,U)(Gauss linearUpwind grad(U))

fvSolution: PIMPLE{nOuterCorrectors 2; nCorrectors 1}

```

O型メッシュなら blockMeshDict で作れる。円柱をブロックで囲み、arc エッジで円弧を定義する。ただし高Re数で境界層を精密に解像する場合は、snappyHexMesh や外部メッシャ（Pointwise, ICEM CFD, Gmsh）を使うほうが効率的だ。

Fluent では以下が推奨だ。

• ソルバー: Pressure-Based, Transient
• 圧力-速度連成: Coupled（収束が速い）またはPISO
• 空間離散化: Second Order Upwind（RANS）、Bounded Central Differencing（LES）
• 時間離散化: Second Order Implicit
• 乱流モデル: SST $k$-$\omega$（RANS）、WALE SGS（LES）、SST-DDES（ハイブリッド）

Fluent の TUI コマンドで渦放出のモニタリングを設定する場合、円柱表面のリフト係数を1タイムステップごとに出力して FFT にかけると St 数が得られる。

判断基準はこうだ。

• 学習・研究目的で低Re: OpenFOAM（無償）または COMSOL（GUIが直感的）
• 工学設計でRANS: Fluent か STAR-CCM+（自動メッシュとワークフロー統合）
• 高精度LES/DNS: OpenFOAM（自由度が高い）またはNek5000（スペクトル要素法）
• VIV（渦励振）解析: Fluent/STAR-CCM+のFSI機能、またはOpenFOAMのoverset mesh

その通り。渦励振はStrouhal周波数と構造の固有振動数が一致するとロックイン現象が起こる。$U_r = U_\infty / (f_n D) = 1/\mathrm{St} \approx 5$ 付近で最も危険だ。