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

球の後流遷移は Johnson & Patel (1999) 以降、精力的に研究されてきた。軸対称 → 非軸対称の遷移は regular bifurcation（ピッチフォーク型）で、Re $\approx 210$ で起こる。このとき後流は1つの対称面のみを持つ「片翼」構造になる。

さらに Re $\approx 270$ で Hopf 分岐が起き、周期的なヘアピン渦の放出が始まる。渦放出の軸方位は Re = 210 で選ばれた対称面に拘束される。

球の場合は片側からヘアピン渦が連続的に放出される。これは後流の対称性が円柱（平面対称）と球（軸対称）で異なるためだ。

抗力危機は境界層の層流→乱流遷移が剥離点を後方に移動させることで起こる。これをCFDで捉えるには遷移モデルが必須だ。

• $\gamma$-$Re_\theta$ 遷移モデル (Langtry & Menter): RANS で遷移を予測。Fluent, CFX, STAR-CCM+ で利用可能
• LES: 格子が十分細かければ遷移を自然に捉えるが、コストが膨大
• DNS: Re $\sim 10^4$ までは実現されている（Rodriguezら 2011）

Achenbach (1972) の実験では、滑らかな球の抗力危機は Re $\approx 3.7 \times 10^5$ で起こる。フリーストリーム乱流強度が $0.45\%$ の場合だ。

回転球には Magnus 力が作用する。スピンパラメータ $\alpha = \omega a / U_\infty$ で無次元化すると、横力係数 $C_L$ は、

サッカーボールの無回転キック（ナックルボール効果）は、回転なし($\alpha \approx 0$)の球で剥離点の位置が不規則に変動し、横力が不安定に振動する現象だ。これは抗力危機のRe域で顕著になる。

IBM は固定直交格子上で物体をペナルティ力や内挿で表現する手法だ。球の沈降、衝突、多体問題に適している。

• Direct Forcing IBM: Uhlmann (2005) の手法。粒子系の大規模DNS に使われる
• Volume Penalization: 物体内部の速度をゼロに強制。実装が簡単
• Cut-Cell法: 物体と格子の交差部で保存的な離散化を行う

その通り。Particle-Resolved DNS (PR-DNS) では IBM が事実上の標準になっている。Breugem (2012) のコードや PeleLM が有名だ。