円柱周りの流れ — 先端技術と研究動向

Williamson (1996) の分類が有名だ。2Dカルマン渦列が3D化する際に2つの不安定モードが現れる。

• Mode A ($Re \approx 190$): スパン方向波長 $\lambda_A \approx 3\text{--}4D$。楕円型不安定性に起因。渦管のうねりとして観察される
• Mode B ($Re \approx 260$): $\lambda_B \approx 1D$。ブレード間のストリーク構造。双曲型不安定性に関連

そうだ。基底状態（2D周期解）に微小な3D擾乱を重ねて Floquet 乗数を求める。$|\mu| > 1$ となるスパン方向波数でモードが不安定化する。Barkley & Henderson (1996) のFloquet安定性解析が決定的な仕事だ。

wall-resolved LES は格子点数が $Re^{13/7}$ でスケールする。高Re数では wall-modeled LES (WMLES) が必須だ。壁面近傍をRANS（例: Spalart-Allmaras）で解き、離壁後にLESに切り替える。

DES (Detached-Eddy Simulation) はまさにこの発想だ。Spalart-Allmaras ベースの DES では、壁面距離 $d$ を $\tilde{d} = \min(d, C_{DES} \Delta)$ で置き換える。DDES はさらに遅延関数を加えて、境界層内でのLES化（Grid Induced Separation）を防ぐ。

Nek5000 や Nektar++ では高次スペクトル要素法（$p$次 $= 7\text{--}15$）を使う。指数的収束により、同じ精度を少ない自由度で達成できる。Karniadakis & Sherwin のグループが円柱DNSで多くの成果を出している。

実用的にも重要なテーマだ。

• スプリッタープレート: 円柱背面に薄板を設置。渦の相互作用を抑制
• 吹出し/吸込み: 壁面からの周期的吹出しで渦放出を制御
• 回転制御: 円柱自体を回転させる。Magnus 効果と渦抑制の組み合わせ
• フィードバック制御: 圧力センサで渦位相を検出し、アクチュエータで対向渦度を注入

最近は強化学習（Deep Reinforcement Learning）を使った最適制御が注目されている。Rabault et al. (2019) は2D円柱の渦放出を強化学習エージェントが完全に抑制する結果を報告している。

3つのアプローチがある。

• ALE法 (Arbitrary Lagrangian-Eulerian): メッシュ自体が変形。小振幅向き
• Overset（Chimera）メッシュ: 円柱周りのメッシュが背景メッシュ上をスライド。大変位対応
• Immersed Boundary法: 固定格子上に仮想力で物体を表現。メッシュ再生成不要

そうだ。OpenFOAM v2306 以降では overset の安定性がかなり改善されている。Fluent でも overset mesh が使えるし、STAR-CCM+ は overset の実績が豊富だ。