🧑🎓
このシミュレーターで、マウスでぐるぐるかき混ぜると渦ができるけど、これって現実の水や空気の動きと同じ仕組みなんですか？
🎓
ざっくり言うと、その通りだよ。このツールは「渦法」という、流体を小さな渦の粒（渦要素）の集まりとして計算する方法を使っているんだ。マウスで引く軌道が渦度（うずど）という「回転の強さ」を与えて、それが周りの流体を巻き込んでいく様子を再現している。実務のCAEでも、乱流の大きな構造を解析する時に似た手法が使われることがあるね。
🧑🎓
「粘性」ってパラメータを変えると渦の広がり方が変わるみたいだけど、あれは何を表してるんですか？
🎓
粘性は、流体の「ねばつき」や「摩擦」の度合いだと思って。値が大きい（減衰率が高い）と、渦のエネルギーがすぐに周りに拡散して消えていく。逆に小さいと、きれいな渦が長く保たれる。上のスライダーで粘性をゼロに近づけてかき混ぜてみると、渦がほとんど減衰せずに複雑に絡み合うのがわかるよ。これは、水より空気の方が粘性が低いことにも対応しているんだ。
🧑🎓
プリセットの「カルマン渦列」って、円柱の後ろにできる交互の渦ですよね。これが実際に問題になるって本当ですか？
🎓
本当だよ。これが「渦励振（Vortex Induced Vibration）」という現象の原因になる。例えば、橋のケーブルや煙突、海底パイプラインに風や海流が当たると、この交互に剥離する渦が構造物に周期的な力を加える。シミュレーターで「渦核半径ε」を変えてみると、渦の大きさや間隔が変わって、発生する力の周期も変わるのが体感できる。CAEでは、この振動による疲労破壊を予測するために詳細に解析するんだ。
このシミュレーターの核となるのは「ランキン渦モデル」です。これは、渦の中心（渦核）では剛体のように回転し、外側では回転速度が距離に反比例して減衰する理想的な渦を表すモデルで、中心で速度が無限大になる特異性を避けるために使われます。
$$ v_{\theta}(r) = \begin{cases}\frac{\Gamma}{2\pi R^2}r & (r \le R) \\
\frac{\Gamma}{2\pi r}& (r > R)
\end{cases}$$
$v_{\theta}$: 中心からの距離$r$における接線方向速度、$\Gamma$: 循環（渦の強さ）、$R$: 渦核半径（パラメータ「渦核半径ε」に対応）。渦核内($r \le R$)では速度が$r$に比例して増加し、渦核外($r > R$)では$1/r$に比例して減少します。
各渦要素が他の位置に誘起する速度は、2次元渦法の基本式である「ビオ・サバールの法則」のアナロジーで計算されます。渦度$\omega$を持つ点渦が、距離$r$離れた点に誘起する速度です。
$$ \Delta \vec{v}= \frac{\vec{\omega}\times \vec{r}}{2\pi |\vec{r}|^2}$$
$\Delta \vec{v}$: 誘起される速度ベクトル、$\vec{\omega}$: 渦度ベクトル（2次元ではz成分のみ）、$\vec{r}$: 渦要素から観測点への位置ベクトル。全ての渦要素からの影響を足し合わせることで、複雑な流れ場が生成されます。
よくある誤解と注意点
このツールで遊び始めるときに、いくつか気をつけてほしいポイントがあるよ。まず、「マウスで速く動かせば強い渦ができる」と思いがちだけど、実はそう単純じゃない。このシミュレーターは「循環Γ」という渦の強さをマウスの軌跡に沿って付与している。つまり、動かす「速さ」よりも、動かした「軌跡が囲む面積」が大きく影響するんだ。ゆっくりでも大きく円を描くほうが、小刻みに速く動くよりも強い渦ができることがある。これは実務でも、境界条件の設定で「速度」と「渦度」を混同しないように、という注意に繋がるね。
次に、パラメータ「粘性」と「渦核半径ε」の関係だ。両方とも渦の広がり方に関わるから、同じような効果に見えるけど、物理的な意味は全く違う。「粘性」は流体そのものの性質で、渦のエネルギーが熱に散逸していく（減衰する）過程を表す。一方、「渦核半径ε」は計算モデル上のパラメータで、渦の中心部の構造をどのくらいの大きさで滑らかに表現するかを決めるもの。例えばεを極端に小さく（0.01など）して粘性もゼロにすると、渦が非常に鋭く、ほとんど減衰しないため、数値的不安定（渦が暴発的に強くなる）が起きやすくなる。実務のCFDでも、モデル化のためのパラメータと物性値は区別して理解することが超重要だ。
最後に、これは「2次元」のシミュレーションだという根本的な制限を忘れないで。画面は平らだよね。現実の流れは3次元なので、ここで美しく見える渦の糸も、実際には渦管として複雑に伸びたり絡まったりする。このツールでカルマン渦列を再現しても、実際の円柱後流では3次元的な乱れ（スパン方向の変動）がすぐに発生する。2次元計算は現象の本質を掴むには最高だが、定量的な値を求める実務解析では常に3次元効果を考慮する必要がある、というのが大きな落とし穴なんだ。
関連する工学分野
この「渦法」をベースにした考え方は、CAE流体解析のさまざまな分野に応用されている。まず挙げるのは航空宇宙工学だ。飛行機の翼から剥離する渦（翼端渦）は、後続機に影響を与える危険な乱気流の原因になる。このような大きな渦構造を効率的に追跡するために、ラグランジュ粒子法の一種である「渦法」が研究されてきた歴史がある。また、ヘリコプターのローター周りの流れも複雑な渦の絡み合いだ。
もう一つは自動車工学、特に空力設計だ。車体後方で発生する乱流（ウェイクフロー）は空気抵抗の主原因。ここで発生する大きな渦構造を把握するために、DES（Detached Eddy Simulation） のような高度な乱流モデルが使われる。DESは物体近くはRANS（レイノルズ平均）、剥離した渦が発達する領域はLES（ラージエディシミュレーション）で解くハイブリッド手法で、このツールで見ているような「はっきりした渦」の挙動を直接計算する思想に近い。
さらに意外なところではプラント配管設計にも関係する。ポンプから吐出された流れがパイプの曲がり部（エルボ）を通るとき、二次流れという渦が発生する。これが下流の流量計の計測精度を狂わせたり、特定部位の腐食を促進したりする。このような内部流れの渦を予測するのもCFDの重要な仕事で、その基礎となる渦度の考え方は共通しているんだ。
発展的な学習のために
このツールに慣れて「もっと知りたい」と思ったら、次のステップとして「渦度方程式」を学ぶことを強くお勧めする。ナビエ-ストークス方程式から導かれるこの方程式は、渦が「生成・輸送・拡散・消滅」するメカニズムを直接記述している。2次元の場合、粘性のある流体では次の形になるよ。
$$ \frac{\partial \omega}{\partial t} + (\vec{v} \cdot \nabla) \omega = \nu \nabla^2 \omega $$
左辺第2項が渦の輸送（このシミュレーターで粒子で可視化している部分）、右辺が粘性による渦の拡散（パラメータ「粘性」で調整している部分）に対応する。この式を見ると、3次元ではもっと複雑な「渦の伸長」の項が加わる。教科書なら「流体力学（前編）」（今井功）の渦運動の章が非常に明快だ。
学習の実践としては、このシミュレーターで遊びながら、パラメータを変えた時の挙動を仮説→検証のサイクルで記録してみよう。例えば「渦核半径εを大きくすると、二つの接近した渦の合体は早くなるのか遅くなるのか？」などだ。その感覚を養った後で、実務レベルのCFDソフト（OpenFOAMなど無料のものもあり）で「2次元キャビティ流れ」などの基礎チュートリアルを実行してみると、メッシュやソルバーといった新しい要素が出てくるが、背景にある渦の物理は変わらないことが実感できるはずだ。次のトピックとしては、乱流のエネルギーカスケードや、流体構造連成（FSI）の基礎に進むと、このツールで見ていた現象が、どのようにして実際の「構造物の振動」という問題に発展していくのかが理解できるようになるよ。