バーガース渦
バーガース渦またはバーガース・ロット渦(バーガースうず、バーガース・ロットうず、Burgers vortex, Burgers–Rott vortex)とは、粘性を考慮したナビエ–ストークス方程式の厳密解であり、定常で自己相似な流れによる渦である。ヤン・バーガースとニコラス・ロットの名前に由来する。
対称軸中心向きの放射状の流れは、軸周りに渦度を集中させる傾向にある。また同時に粘性拡散は渦度を拡散する傾向にある。この2つの効果が釣り合うときに定常なバーガース渦が発生する。
バーガース渦は竜巻のような連続的な対流駆動などによる渦の伸長によって渦度が与えられる流れを説明することができる。
バーガース渦の流れは円筒座標系 (r,\theta,z) で表現される。軸対称(\theta に依存しない)を仮定し、軸に対称な滞留点のある流れ場を考える。
v_r= -\alpha r,
v_z=2\alpha z,
v_\theta=\frac{\Gamma}{2\pi r}g(r).
ここで \alpha > 0 は伸長度、\Gamma > 0 は循環である。最初の二式 v_r, v_z により、流れ場は連続の式を満たす。
ナビエ–ストークス方程式の圧力項と外力項を無視した方位角方向の運動方程式は次のように表現される。
r \frac{\mathrm{d}^2 g}{\mathrm{d} r^2} + \left(\frac{\alpha r^2}{\nu} -1 \right) \frac{\mathrm{d}g}{\mathrm{d} r} = 0.
ここで \nu は動粘性係数である。
無限遠ではポテンシャル渦のような振る舞いをし、有限の場所では回転する流れとなるように、方程式は g(\infty)=1 を満たすように積分される。この仮定により、軸上ではg(0)=0 すなわち v_\theta=0が保証される。このとき解は
g= 1 - \exp \left(-\frac{\alpha r^2}{2\nu}\right).
したがって渦度は非自明な z 軸方向の成分のみ与えられ、以下のように表現される。
\omega_z = \frac{\alpha\Gamma}{2\pi\nu} \exp \left(-\frac{\alpha r^2}{2\nu}\right).
\omega_z のための渦度方程式の3つの項によって直感的に流れ場を理解することができる。軸方向の速度 v_z は渦の伸長によって軸方向の渦中心の渦度を強める。強まった渦度は動径方向に拡散しようとするが、v_r による動径方向の渦対流により拡散は阻止される。この3つの流れのバランスにより、定常的な渦場が得られる。
バーガース渦レイヤーまたはバーガース渦シートは、バーガース渦を2次元的に近似したシアー層である。これもナビエ–ストークス方程式の厳密解であり、1951年にA.A.Townsendによって最初に記述された。直交座標系で表される速度場 (v_x,v_y,v_z) は
v_x= -\alpha x,
v_z=\alpha z,
v_y= U\mathrm{erf}\left(
\frac{\sqrt\alpha x}{\sqrt{2\nu}\right).
ここで \alpha > 0 は伸長度、v_y(+\infty) = U、v_y(-\infty) = -U である。2U の値は渦レイヤーの強度と解釈できる。渦度は非自明な z 軸方向の成分のみ与えられ、以下のように表現される。
\omega_z = \frac{2U}{\sqrt{2\pi}\sqrt{\frac{\alpha}{\nu} \exp \left(-\frac{\alpha x^2}{2\nu}\right).
非軸対称な流れ場の中では、非軸対称なバーガース渦が現れる。渦のレイノルズ数 Re=\Gamma/(2\pi \nu) が小さい場合の非軸対称なバーガース渦の理論は1984年に A. C. Robinson と Philip Saffman によって構築されたが、1994年にはKeith Moffatt、S. Kida、K. Ohkitani によってレイノルズ数が大きい場合 Re\gg 1 の理論を構築した。
渦のレイノルズ数を任意に設定した場合の非軸対称なバーガース渦は数値積分によって考えることができ、速度場は次のような形になる。
v_x = -\alpha x + u(x,y),
v_y = -\beta y + v(x,y),
v_z = \gamma z.
ここで \gamma=\alpha+\beta である。一般性を損なわないように \alpha > 0、\gamma > 0 を仮定する。渦の断面は xy 平面上にあり、z 方向に0ではない渦度成分が存在し、以下のように表現される。
\omega_z = \frac{\partial u}{\partial y}-\frac{\partial v}{\partial x}.
軸対称なバーガース渦は \alpha=\beta=\gamma/2 を満たす場合であり、バーガース渦レイヤーは \alpha=\gamma and \beta=0 を満たす場合である。
伸長した円筒座標の淀み面上におけるバーガース渦についてのナビエ–ストークス方程式の陽解法が2021年に議論された。
円筒座標系において方程式は以下の速度場に従う。
v_r= -\alpha \left(r - \frac{r_s^2}{r}\right),
v_z=2\alpha z,
v_\theta= \frac{\Gamma}{2\pi r}P\left(1+ \frac{\alpha r_s^2}{2\nu},\frac{\alpha r^2}{2\nu}\right),
ここで \alpha > 0 は伸長度、r_s \geq 0 は円筒座標上の淀み面の場所である。\Gamma > 0 は循環。P は不完全ガンマ関数である。
この解はポテンシャル流による線状の強制 Q=2\pi \alpha r_s^2 が存在する場合のバーガース渦の表現に他ならない。
渦度は非自明な z 軸方向の成分のみ与えられ、以下のように表現される。
\omega_z = \frac{\alpha\Gamma}{2\pi\nu \tilde\Gamma(1+\alpha r_s^2/2\nu)}\left(\frac{\alpha r^2}{2\nu}\right)^{\alpha r_s^2/2\nu} \exp \left(-\frac{\alpha r^2}{2\nu}\right)
ここで \tilde\Gamma はガンマ関数である。r_s\rightarrow 0 のとき関数はバーガース渦へ近づき、r_s\rightarrow \infty のときバーガース渦レイヤーの解に近づく。
円筒座標上の淀み面におけるサリバン渦の陽な解も存在する。
1959年に Roger D. Sullivan が以下のような速度場を考えてバーガース渦を拡張した。
v_r=- \alpha r + \frac{1}{r} f(\eta),
v_z=2\alpha z\left[1-\frac{f'(\eta)}{2\nu}\right],
v_\theta=\frac{\Gamma}{2\pi r}\frac{g(\eta)}{g(\infty)}.
ここで \eta = \alpha r^2/(2\nu) である。関数 f(\eta) と g(\eta) は次のように表現される。
f(\eta) = 6\nu (1-e^{-\eta}),
g(\eta)=\frac{\nu}{\alpha} \int_0^\eta \exp \left[-t + 3 \int_0^t \frac{1-e^{-s}{s}\mathrm{d}s\right]\mathrm{d} t.
バーガース渦において v_r < 0 であり、v_\theta > 0 と v_z/z>0 は常に正である。サリバンの結果は \eta \leq 2.821 のとき v_r > 0 となり、\eta \leq 1.099 のとき v_z/z < 0 となる。
したがってサリバン渦は \eta > 2.821 のときバーガース渦と似た構造をもつが、v_z/z の符号が軸付近で異なる2セル構造が発達する。