渦度

200px|right|thumb|北半球における高気圧 (H) ・低気圧 (L) の回転方向
**渦度**（うずど、、）は、流れの回転するありさまを表現する量である。渦度はベクトル量（さらに言えば擬ベクトル）であり、流れの速度ベクトルのなすベクトル場の回転である。

渦度ベクトル  は流速ベクトル  により、以下のように表される。
:\begin{align} \boldsymbol{\Omega} &= \operatorname{rot} \boldsymbol{v} = \nabla \times \boldsymbol{v} \\
&= \begin{pmatrix}\dfrac{\partial v_z}{\partial y} - \dfrac{\partial v_y}{\partial z},& \dfrac{\partial v_x}{\partial z} - \dfrac{\partial v_z}{\partial x},& \dfrac{\partial v_y}{\partial x} - \dfrac{\partial v_x}{\partial y} \end{pmatrix}\end{align}

渦度ベクトルを流線のようにつなげた曲線を**渦線**という。流れの中のある閉曲線上の各点を通る渦線によってできる曲面を**渦管**という。断面積を無限小にした渦管を**渦糸**という。渦糸が閉曲線になっている場合、これを**渦輪**という。

## 性質
* 流れの中に微小領域をとったとき、微小時間で考えれば、その領域は渦度ベクトルの方向を軸に角速度  で剛体的に回転する。
* 1本の渦管を考えたとき、渦管の表面を一周する任意の閉曲線  に沿った循環（渦度  の面積分）
*::\Gamma=\oint_{\partial S} \boldsymbol{v}\cdot \mathrm d\boldsymbol{s}=\int_{S} \boldsymbol{\omega}\cdot \mathrm d\boldsymbol{S}
*:は、渦管に固有の量となり、**渦管の強さ**と呼ばれる。ここで  は閉曲線  の線要素、 は  で囲まれる曲面の面要素である。
* 同様に、渦糸に対してその断面積  と渦度  との積  は渦糸に固有の量であり、**渦糸の強さ**と呼ばれる。このことから、渦糸が細いところほど渦度が大きくなることが分かる。
* 渦管や渦糸は流体内部の点で途切れることはなく、流体領域の境界（無限遠を含む）まで伸びているか、閉曲線となり渦輪を作るかのどちらかである。

### テンソルとしての性質
渦度は、速度勾配テンソルの反対称成分
:
\big( \partial_{[i} v_{j]} \big) =
\frac{1}{2}
\begin{pmatrix}
0 &
\dfrac{\partial v_1}{\partial x_2}-\dfrac{\partial v_2}{\partial x_1} &
\dfrac{\partial v_1}{\partial x_3}-\dfrac{\partial v_3}{\partial x_1}\\\\
\dfrac{\partial v_2}{\partial x_1}-\dfrac{\partial v_1}{\partial x_2} &
0 & 
\dfrac{\partial v_2}{\partial x_3}-\dfrac{\partial v_3}{\partial x_2}\\\\
\dfrac{\partial v_3}{\partial x_1}-\dfrac{\partial v_1}{\partial x_3} &
\dfrac{\partial v_3}{\partial x_2}-\dfrac{\partial v_2}{\partial x_3} &
0
\end{pmatrix}

を、関係式
:\frac{1}{2}\left(\frac{\partial v_i}{\partial x_j}-\frac{\partial v_j}{\partial x_i}\right)=\epsilon_{jik}\omega_k
を用いてベクトルとして表したものである。ここで  はエディントンのイプシロンである。

## ビオ・サバールの法則
無限に広がる流体領域の渦度分布が与えられており、かつ無限遠で速度が0であるとき、速度分布はビオ・サバールの法則により求められる。
:\begin{align}\boldsymbol{v}&=\frac{1}{4\pi}\iiint\frac{\boldsymbol{\omega}\times\boldsymbol{r}}{r^3} dV \\
&=\operatorname{rot}\iiint\frac{\boldsymbol{\omega}}{4\pi r}dV\end{align}

## 渦度方程式
流体の運動方程式の回転を取ることによって、次の**渦度方程式**が導かれる。
:\frac{\partial\boldsymbol{\omega}}{\partial t}=\operatorname{rot}(\boldsymbol{v}\times\boldsymbol{\omega})+\nu\nabla^2\boldsymbol{\omega}
ここで&nu;は流体の動粘性係数である。

## 大気力学における渦度
大気力学においては、渦度ベクトルの中で3方向の成分の中で特に重要なのは鉛直方向の成分であり、これは**鉛直渦度**と呼ばれ、単に渦度といった場合、暗黙の了解として鉛直渦度を指すことが多い。鉛直渦度は、東西南北方向の長方形の単位領域に対して、長方形の辺に沿う方向の風速を考え、
:（上辺の西向きの風速－下辺の西向きの風速）／縦方向の距離  - （右辺の南向きの風速－左辺の南向きの風速）／横方向の距離
により与えられる。北半球においては、低気圧の鉛直渦度は正になり、高気圧の鉛直渦度は負になる。南半球においては逆符号になる。

地表からは静止しているように見える空気も地球の自転に伴って運動しているので、慣性系から見れば渦度を持っている。この自転に伴う渦度を**惑星渦度**という。惑星渦度  は （ は自転の角速度、 は緯度）である。惑星渦度と地表から見たときの風による渦度（**相対渦度**  という）と足し合わせた渦度、つまり慣性系から見たときの風による渦度を**絶対渦度**という。

中緯度の総観気象規模における渦度方程式は次の式で表される。

\frac{D_h (\zeta+f)}{Dt} = -(\zeta+f)\biggl(\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}\biggr)

ただし、\frac{D_h}{Dt} \equiv \frac{\partial}{\partial t} + u\frac{\partial}{\partial x} + v\frac{\partial}{\partial y} である。

風に発散や収束が無い場合、上の式の右辺がゼロになる。すなわち、空気が別の場所へ移動してもその空気の絶対渦度は保存される。大気中では  の面がこの状況に近い。そのため将来の渦度の分布を予測することが可能であり数値予報の重要な要素となっている。

## 脚注

### 注釈

### 出典

## 関連項目
* 循環 (流体力学) - ストークスの定理によって、渦度と関連付けられる。
* 空気砲 (科学教材) - 渦輪を可視化する工作物。