バーガース渦

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