渦度方程式 — 実践ガイドとベストプラクティス

渦度場の可視化は流れ場の理解に極めて有用だ。主な活用方法を紹介しよう。

• 渦度の等値面（3D）: 翼端渦やKarman渦列の構造を可視化。Ansys FluentのContour機能でVorticity Magnitudeを選択
• 渦度の等値線（2D）: 渦の発達・剥離の追跡。ParaViewではCompute Derivativesフィルターで渦度を計算可能
• Q基準・λ2基準: 渦度だけでは純粋せん断と渦の区別ができないため、渦構造の同定には $Q = \frac{1}{2}(\|\Omega\|^2 - \|S\|^2) > 0$ が使われる

Fluent では Custom Field Function で 0.5*(Vorticity_Magnitude^2 - Strain_Rate_Magnitude^2) を定義すれば計算できる。STAR-CCM+ では Field Function で同様に定義する。OpenFOAM では Q フィールドを postProcess -func Q で直接出力可能だよ。

Lid-driven cavity（蓋駆動キャビティ）は最も広く使われるベンチマークだ。正方形領域の上壁が一定速度 $U$ で移動し、他の3壁は静止。Ghia et al.（1982）の結果が比較データとして有名だよ。

Re が大きくなるほど主渦の中心がキャビティ中央に移動し、コーナー部の二次渦が発達する。Re=5000 では下部コーナーに明確な二次・三次渦が現れる。メッシュが粗いとこれらが捕捉できないから、収束性確認の良い指標になるんだ。

もちろん。Re=100 の円柱周りでカルマン渦が形成される問題は、渦度場の非定常解析の好例だ。ストルーハル数 $St = fD/U \approx 0.164$（Re=100のとき）が検証の目安になる。

実装上のポイントは以下だ。

• 計算領域: 円柱から上流に $10D$、下流に $30D$ 以上確保
• メッシュ: 円柱表面に少なくとも100点以上配置。$y^+<1$ なら壁面第一層厚 $\Delta y \approx D/Re^{0.5}$
• 時間刻み: 渦放出周期を1周期あたり200ステップ以上で解像
• 初期擾乱: 完全対称な初期条件では渦放出が始まらないため、微小擾乱を加える

Fluent では Transient 解析を選択し、Pressure-Based ソルバーでSIMPLE法を使う。空間離散化は Second Order Upwind、時間離散化は Second Order Implicit だ。渦放出の周期は Monitor で Lift coefficient の時間変化を記録すると把握できる。

渦度場は速度勾配の微分だから、速度場より1階高い空間解像度が必要だ。具体的には

• 渦核の解像: 渦核直径に最低10セル以上。渦が移流する領域全体で均一に細かくする必要がある
• 数値拡散の抑制: 1次精度風上スキームは渦度を過剰に減衰させる。必ず2次精度以上を使うこと
• 非構造格子での注意: ポリヘドラルメッシュは六面体に比べて渦度の数値拡散が大きい。渦の軌跡が予測できる場合はその領域を六面体メッシュで解像するとよい

そこがまさに渦度法（vortex method）の出番だ。ラグランジュ的に渦要素を追跡すれば、渦がない領域に格子が不要になる。ただし壁面を含む問題では境界処理が複雑になるので、一長一短だね。