円柱まわりの流れ — CAE用語解説

いくつかの理由がある。①形状がシンプルで解析解や実験データが豊富——Williamson（1996年）やPapaioannou（2006年）など権威ある実験・数値データが公開されている。②レイノルズ数ReによってStokes流れ（Re<1）→定常はく離（Re=10〜40）→カルマン渦列（Re=40〜200）→乱流（Re>1000）と多様な流れ体制をカバーする。③カルマン渦列のStrouhal数（St = f*D/U）が実験値St≈0.2と比較できて、渦の発生周波数を正しく捕捉できるかというソルバー検証に使える。CFDの動的精度を見るのに最適な問題だよ。

円柱後方に周期的に交互に渦が放出される現象で、1911年にTheodore von Karmanが理論解析したことから名前がついた。Re ≈ 40を超えると円柱後ろに形成された対称なはく離泡が不安定になって、上下交互に渦が放出され始める。これが円柱に周期的な揚力変動を生んで——この振動が構造物の固有振動数と一致すると共振する「渦励振（VIV: Vortex Induced Vibration）」が起きる。橋梁や煙突が風でうなる音の原因で、タコマナローズ橋崩落（1940年）もこの一種だ。

Re=100のカルマン渦を正確に捕捉するには——①2D解析ならブロック六面体メッシュで円柱表面をy+ < 1に保つ（乱流モデルが不要なLaminar解析でも境界層解像が重要）、②時間刻みをCo < 0.5程度に保つ（CFL条件）、③計算時間が十分長く（定常な周期振動が確立するまで）解く——が基本だ。検証指標は抗力係数Cd、揚力係数Clの振幅、そしてStrouhal数St。St = 0.166〜0.195（Re依存）という実験値と±5%以内に一致できれば合格ラインだ。

Re > 200程度から3D効果が出始める。2D解析では渦構造がz方向に無限に伸びた理想的な状態だが、実際の3D流れではz方向にも不安定性が発展して渦構造がねじれ・分裂する——これをMode A/Bと呼ぶ。Re = 190〜260でMode A（波長≈3〜4D）が発現、Re = 250〜300でMode B（波長≈0.8〜1D）が現れる。3D LES解析はこれを正確に捕捉できるが計算コストが大幅に増える。工学的目的なら「渦放出周波数を求めたい」だけなら2D非定常RANSで十分なことが多い。

煙突・橋梁・海洋ライザー管の渦励振設計が典型だ。円柱モデルで渦放出周波数 f = St * U / D を計算して、構造物の固有振動数と比較する。もし一致するとVIVが発生するから——①形状を変える（断面を楕円やD型にして渦放出を弱める）、②構造物に制振ダンパーを追加、③螺旋状フィン（ストレークス）を取り付けて渦の規則性を壊す——という対策を打つ。石油プラットフォームの海中ライザー管は潮流によるVIVが疲労設計の支配要因で、CFDとFEM疲労の連成解析が実務で行われている。