円柱周りの流れ

見た目はシンプルだが、Re数によって流れの構造が劇的に変わる。定常対称剥離からカルマン渦列、さらには乱流遷移まで、流体力学の教科書に必ず出てくる古典的かつ本質的な問題だよ。

Re が約47を超えると、後流の対称性が絶対不安定になる。片側の剥離せん断層が反対側の渦度を巻き込むことでフィードバックループが生じ、上下交互に渦が放出される。これが Benard-von Karman 渦列だ。

整理するとこうなる。

そうだ。この閾値は Hopf 分岐に対応している。線形安定性解析で臨界 Re を精密に求めると $Re_{cr} \approx 46.7$ という値が得られる。

非圧縮性 Navier-Stokes 方程式だ。

ここで $\nu = \mu / \rho$ は動粘度、$\mathbf{u}$ は速度ベクトル、$p$ は圧力だ。無次元化すると Reynolds 数 $Re = U_\infty D / \nu$ が唯一の支配パラメータになる。

Strouhal 数で整理する。

亜臨界域（$300 < Re < 3 \times 10^5$）では $\mathrm{St} \approx 0.2$ がほぼ一定値を取る。これは Roshko の実験で確認された有名な結果だ。低 Re 側では Re の関数として次の経験式がよく使われる。

単位スパン長さあたりの定義はこうだ。

定常流の $C_D$ はRe によって大きく変わる。Stokes 領域では $C_D \propto Re^{-1}$、亜臨界域では $C_D \approx 1.0 \text{--} 1.2$ 、臨界域で $C_D \approx 0.3$ まで低下し、超臨界域で再び $C_D \approx 0.6 \text{--} 0.7$ まで回復する。

揚力の変動周波数は渦放出周波数 $f$ そのものだ。一方、抗力の変動周波数は $2f$ になる。上下交互に渦が出るたびに抗力は増減するが、揚力は片側の渦放出で正負が入れ替わるからだ。

Re数によって最適な手法が変わる。整理するとこうだ。

高Re数のDNSは格子点数が $Re^{9/4}$ でスケールする。Re=$10^6$ なら $10^{13}$ 点以上必要で、現状のスパコンでも非現実的だ。

非圧縮Navier-Stokesでは圧力のPoisson方程式を解く必要がある。代表的な手法は以下の通りだ。

その通り。OpenFOAM の pimpleFoam では PISO ループと SIMPLE の外部反復を組み合わせた PIMPLE 法が使えるから、大きめの時間刻みでも安定して計算できる。

対流項の扱いが最も重要だ。

カルマン渦を解像したいなら、最低でも2次精度が必要だ。1次風上では数値拡散で渦が消滅する。

渦放出を正しく捉えるには時間刻みが重要だ。CFL 条件 $\mathrm{CFL} = u \Delta t / \Delta x < 1$ を満たすのが基本で、非定常計算では2次精度の陰解法（例: 後退2次差分 BDF2）がよく使われる。Strouhal 数 $\mathrm{St} \approx 0.2$ から渦放出周期 $T = D / (\mathrm{St} \cdot U_\infty)$ が分かるので、1周期あたり少なくとも200ステップ以上は取りたい。

円柱周り流れのメッシュ設計のポイントはこうだ。

$y^+ = y u_\tau / \nu$ で、$u_\tau = \sqrt{\tau_w / \rho}$ は摩擦速度だ。亜臨界域の $C_f$ から見積もると、Re=$10^4$ で第1層厚さは $10^{-3}D$ 程度になる。これが高Re円柱LESのコスト要因だ。

以下のフローが標準的だ。

各境界の設定はこうだ。

良い質問だ。出口が近すぎると渦が境界に到達して反射波が戻ってくる。だから出口を $20D$ 以上離すか、対流条件（convective outflow: $\partial \phi / \partial t + U_c \partial \phi / \partial x = 0$）を使うのが望ましい。

以下のベンチマークデータと比較するんだ。

そうだ。まず Re=100 の2D計算で $C_D$ と St を検証するのが定石だ。ここで合わなければメッシュか時間刻みに問題がある。

少なくとも3水準のメッシュで Grid Convergence Index (GCI) を計算する。Richardson 外挿を適用して、離散化誤差のオーダーと格子収束解を推定するんだ。

Re=100 の2D計算なら数万セルで十分だ。Re=$3900$ の3D LES だと $500$ 万〜$2000$ 万セル、Re=$10^5$ のDES だと $1000$ 万〜$5000$ 万セル程度が目安になる。

各ツールの特徴を円柱流れの観点で比較しよう。

pimpleFoam を使う場合のディレクトリ構成と主要設定はこうだ。

O型メッシュなら blockMeshDict で作れる。円柱をブロックで囲み、arc エッジで円弧を定義する。ただし高Re数で境界層を精密に解像する場合は、snappyHexMesh や外部メッシャ（Pointwise, ICEM CFD, Gmsh）を使うほうが効率的だ。

Fluent では以下が推奨だ。

Fluent の TUI コマンドで渦放出のモニタリングを設定する場合、円柱表面のリフト係数を1タイムステップごとに出力して FFT にかけると St 数が得られる。

判断基準はこうだ。

その通り。渦励振はStrouhal周波数と構造の固有振動数が一致するとロックイン現象が起こる。$U_r = U_\infty / (f_n D) = 1/\mathrm{St} \approx 5$ 付近で最も危険だ。

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

そうだ。基底状態（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で多くの成果を出している。

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

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

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

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

実務で頻出するトラブルを整理しよう。

原因と対策:

実際の流れには必ず微小な擾乱がある。数値計算ではそれがないため、不安定でも分岐が起きないことがある。$y$方向速度に $O(10^{-3})$ 程度のランダムノイズを加えれば十分だ。

確認ポイント:

Re=100 で $C_D > 1.4$ なら以下を確認だ。

症状別の対策:

スパン方向の計算長さが不足しているか、スパン方向のメッシュが粗すぎる可能性がある。Mode A を捉えるには $L_z > 4D$、Mode B を捉えるには $\Delta z < 0.5D$ が必要だ。また、スパン方向端面は必ず周期境界条件にすること。壁面条件にすると端部効果が支配的になる。

その通り。まず2D・低Reで $C_D$, St を検証 → 3D低Reで Mode A/B を確認 → 目標Reに段階的に上げる、というのが最も確実なアプローチだ。