渦度方程式

渦度 $\boldsymbol{\omega}$ は流体の局所的な回転を表すベクトル量で、速度場の回転（curl）として定義される。

その通り。例えば直線的なせん断流 $u = ky$、$v = 0$ でも $\omega = -k \neq 0$ だから渦度がある。つまり渦度は「渦」の存在とは別の概念なんだ。

Navier-Stokes方程式の両辺に $\nabla \times$ を作用させると渦度方程式が得られる。非圧縮性流体の場合こうなる。

第1項 $(\boldsymbol{\omega} \cdot \nabla)\mathbf{u}$ は渦管の伸長・傾斜項（vortex stretching）だ。3次元流れで渦管が引き伸ばされると断面が細くなり、角運動量保存から渦度が増大する。竜巻やバスタブの排水渦が細くなると回転が速くなるのはこの効果だよ。

第2項 $\nu \nabla^2 \boldsymbol{\omega}$ は粘性拡散項で、渦度が分子粘性により拡散する効果を表す。動粘性係数 $\nu$ が大きいほど渦度は速やかに拡散して消散する。

よく気づいたね。2次元では $\boldsymbol{\omega}$ は $z$ 方向のみで $(\boldsymbol{\omega} \cdot \nabla)\mathbf{u} = 0$ になるから、渦度方程式は単純な移流拡散方程式になる。

非圧縮・バロトロピック流体では、渦度は流体内部では生成されない。渦度の主な発生源は以下だ。

壁面近傍のメッシュ解像度が決定的に重要だ。壁面第一セルの $y^+$ を1以下にするか、壁関数を使って壁面せん断応力から渦度フラックスを近似するかの選択になる。

非粘性のバロトロピック流体でポテンシャルな体積力のみが作用する場合、閉曲線に沿った循環 $\Gamma = \oint \mathbf{u} \cdot d\mathbf{l}$ は保存される。Stokesの定理により $\Gamma = \int_S \boldsymbol{\omega} \cdot d\mathbf{S}$ なので、これは渦度方程式で粘性項とバロクリニック項がゼロのときに対応する。

その理解で正しい。渦度方程式は粘性・非バロトロピック効果を全て含んだ、より一般的な記述なんだ。

2次元非圧縮流れでは渦度-流れ関数法（vorticity-stream function method）が古典的な手法だ。渦度方程式と流れ関数のPoisson方程式を連立させる。

2次元では未知数が $\omega$ と $\psi$ の2つだけで、圧力-速度連成の問題（SIMPLE法などの反復が不要）を回避できる。ただし3次元への拡張は渦度が3成分ベクトルになるため複雑で、実用的にはあまり使われない。

有限差分法（FDM）での典型的な離散化を見てみよう。均一格子間隔 $h$ の中心差分を用いると、Poisson方程式は

FVMではセル界面での渦度フラックスを評価する。移流フラックスには2次精度風上スキーム（Second Order Upwind）やTVDスキーム、拡散フラックスには中心差分を使う。OpenFOAMのscalarTransportFoamソルバーが渦度の移流拡散をFVMで解く実装の参考になるよ。

代表的な時間積分スキームをまとめよう。

もちろん。移流のCFL条件は $\Delta t < h/|u_{max}|$ だ。拡散の安定性条件は $\Delta t < h^2/(4\nu)$ で、高Reynolds数では移流のCFLが支配的、低Reynolds数では拡散の制約が厳しくなる。

壁面での渦度は直接与えられないから、流れ関数から間接的に導出する。代表的な手法はThomの公式だ。壁面（$j=0$）で滑りなし条件 $u=v=0$ と $\psi=0$ を使うと

Woods の公式は $\omega_{i,0} = -\frac{3\psi_{i,1}}{h^2} - \frac{\omega_{i,1}}{2}$ で3次精度になる。またBriley の公式は壁面近傍の渦度分布をTaylor展開して4次精度を達成する。精度が上がるほど収束が遅くなる傾向があるので、実用上はThomの公式で十分な場合が多い。

その通りで、ここがボトルネックになる。代表的な解法は以下だ。

格子数が $100 \times 100$ 程度ならSOR法で十分だ。$1000 \times 1000$ 以上ならマルチグリッド法が圧倒的に速い。OpenFOAMのfvMatrixではGAMGソルバーがマルチグリッド法に対応している。

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

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のとき）が検証の目安になる。

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

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

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

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

渦度そのものを直接解く専用ソルバーは少ないが、N-S方程式を解いた結果から渦度場を算出する機能はどのツールにもある。渦度解析に関連する各ツールの特徴を比較しよう。

OpenFOAMのソルバーフレームワークを使えば、渦度輸送方程式のソルバーを自作できる。scalarTransportFoamをベースに渦度のスカラー輸送方程式を実装し、別途Poisson方程式で流れ関数を求めるアプローチだ。研究用途では実際にそういう実装が公開されているよ。

後処理での渦度解析機能を比較しよう。

STAR-CCM+のVortex Core機能は速度勾配テンソルの固有値から渦核線を自動抽出できる。航空機の翼端渦やターボ機械の渦構造解析で重宝するよ。ParaView にもVTKベースの渦核抽出フィルタがあるが、STAR-CCM+ のほうが使いやすい印象だ。

渦度解析が特に重要になるのはLES（Large Eddy Simulation）やDES（Detached Eddy Simulation）だ。渦構造を時間的に追跡する必要があるからね。

研究目的でカスタマイズしたい場合はOpenFOAM一択だ。産業用途で渦構造の可視化と定量評価を効率的に行いたい場合はSTAR-CCM+が強い。Fluentは汎用性が高く、UDFで拡張すればほぼ何でもできるが、渦度に特化した機能はやや少ないね。

渦度解析そのものは追加ライセンスが不要なケースが多い。ただしLES/DESの時間解像解析は大量の計算リソースが必要だ。

コスト面ではそうだが、セットアップの手間とサポートの有無も考慮する必要がある。初めてLESをやるなら、まずFluentやSTAR-CCM+で小規模なケースを回して感覚を掴み、大規模化するときにOpenFOAMに移行するのも一つの戦略だよ。

渦法（Vortex Method）はラグランジュ的に渦要素を追跡する手法だ。渦度場を離散的な渦要素（点渦、渦ブロブ）の集合で表現する。

Biot-Savartの法則を使って渦度分布から速度を復元する。2次元では

渦が局在する外部流れ問題に強い。翼の後流、噴流混合、渦輪の衝突などだ。一方、壁面を含む内部流れや乱流の壁面層の解像は苦手で、格子法とのハイブリッド手法が研究されている。

乱流のエネルギーカスケードは渦度の観点から理解できる。大スケールの渦がvortex stretchingにより小スケールの渦を生成し、最終的に粘性でエンストロフィー（渦度の二乗の積分）が散逸する。

これはミレニアム懸賞問題の一つだ。3次元N-S方程式の滑らかな解が有限時間で特異性を持つかどうかは未解決。数値的には超高解像度DNS（$4096^3$格子など）で調べられているが、決定的な結論は出ていない。

通常のDNSでは原始変数（速度・圧力）を解いて渦度は後処理で計算する。ただし渦度ベースの定式化を使ったDNSも研究されていて、Kim & Moin（1985）が先駆的だ。

注目すべき研究動向は以下だ。

実験データが速度場の一部しか得られない場合でも、渦度方程式という物理法則を制約として入れることで、観測されていない領域の渦度場を推定できる。PIV（粒子画像流速計）データの補完・超解像に応用されているよ。

渦度解析特有のトラブルとその対策を整理しよう。

これは最もよくある問題だ。原因は数値拡散（numerical diffusion）で、以下をチェックしよう。

Fluentの場合、Spatial Discretization > Momentum を Second Order Upwind 以上にすること。理想的には Bounded Central Differencing（LES用）を使うと数値拡散が最小限になる。

考えられる原因と対策はこうだ。

壁面の渦度は $\omega_{wall} = -\frac{\partial u_{tangential}}{\partial n}$ だから、壁面近傍の速度勾配が正確に解けていないと振動する。対策は

非粘性（$\nu=0$）の2次元流れではエンストロフィーは保存されるはずだ。増大する場合は数値スキームの散逸・分散特性に問題がある。

代表的なトラブルと対策をまとめよう。

NaN発生の主な原因は発散だ。渦度解析で発散しやすいケースと対策を挙げよう。

その通り。渦度は速度の1階微分、ひずみ速度テンソルも1階微分だ。だから速度場の2次精度は渦度の1次精度にしかならない。渦度の精度を上げたいなら、速度場を3次精度以上で解くか、十分なメッシュ解像度を確保する必要がある。