ケルビンの渦定理

非粘性バロトロピック流体の保存外力下での流れにおいて、流体とともに動く閉曲線に沿う循環は時間的に不変である[1]。

数式では

${\displaystyle {\frac {\mathrm {D} {\mathit {\Gamma }}}{\mathrm {D} t}}=0}$

と表現される。

ここで、物質微分 ${\displaystyle \mathrm {D} /\mathrm {D} t}$ は流体と一緒に動く観測者から見た時間変化率、循環 ${\displaystyle {\it {\Gamma }}}$ は流体要素から成る(流体と一緒に動く)閉曲線 ${\displaystyle C(t)}$ 上の流体の速度 ${\displaystyle {\boldsymbol {v}}}$ の線積分

${\displaystyle {\mathit {\Gamma }}(t)=\oint _{C(t)}{\boldsymbol {v}}\cdot \mathrm {d} {\boldsymbol {l}}}$

を表す。

保存外力のもとでの非粘性バロトロピック流体の支配方程式はオイラー方程式 (流体力学) ${\displaystyle {\frac {\mathrm {D} {\boldsymbol {v}}}{\mathrm {D} t}}=-\nabla \int {\mathrm {d} p \over \rho }-\nabla \Omega }$ で表される。ここで、 ${\displaystyle \rho }$ は密度、 ${\displaystyle p}$ は圧力、 ${\displaystyle \Omega }$ は外力のポテンシャルを表す。なお、バロトロピック性( ${\displaystyle \rho =\rho (p)}$ )より ${\displaystyle -{1 \over \rho }\nabla p=-\nabla \int {\mathrm {d} p \over \rho }}$ となることを用いた。 循環の定義式の物質微分をとると ${\displaystyle {\frac {\mathrm {D} {\it {\Gamma }}}{\mathrm {D} t}}=\oint _{C(t)}{\frac {\mathrm {D} {\boldsymbol {v}}}{\mathrm {D} t}}\cdot \mathrm {d} {\boldsymbol {l}}+\oint _{C(t)}{\boldsymbol {v}}\cdot {\frac {\mathrm {D} \mathrm {d} {\boldsymbol {l}}}{\mathrm {D} t}}}$ となる。 第1項に支配方程式を代入すると、 ${\displaystyle {\begin{aligned}\oint _{C(t)}{\frac {\mathrm {D} {\boldsymbol {v}}}{\mathrm {D} t}}\cdot \mathrm {d} {\boldsymbol {l}}&=-\oint _{C(t)}\nabla \left(\int {\mathrm {d} p \over \rho }+\Omega \right)\cdot \mathrm {d} {\boldsymbol {l}}\\&=-\left[\int {\mathrm {d} p \over \rho }+\Omega \right]_{C(t)}\end{aligned}}}$ また、第2項は ${\displaystyle {\begin{aligned}\oint _{C(t)}{\boldsymbol {v}}\cdot {\frac {\mathrm {D} \mathrm {d} {\boldsymbol {l}}}{\mathrm {D} t}}&=\oint _{C(t)}{\boldsymbol {v}}\cdot \mathrm {d} {\boldsymbol {v}}\\&=\oint _{C(t)}\mathrm {d} \left({|{\boldsymbol {v}}|^{2} \over 2}\right)\\&=\left[{|{\boldsymbol {v}}|^{2} \over 2}\right]_{C(t)}\end{aligned}}}$ となるので、 ${\displaystyle {\begin{aligned}{\frac {\mathrm {D} {\it {\Gamma }}}{\mathrm {D} t}}&=\left[{|{\boldsymbol {v}}|^{2} \over 2}-\int {\mathrm {d} p \over \rho }-\Omega \right]_{C(t)}\end{aligned}}}$ となる。ただし、${\displaystyle [f]_{C}}$ は閉曲線 ${\displaystyle C}$ を一周したときの ${\displaystyle f}$ の差を表す。また、 ${\displaystyle {\begin{aligned}\left.{\frac {\mathrm {D} \mathrm {d} {\boldsymbol {l}}}{\mathrm {D} t}}\right|_{\boldsymbol {r}}&={\frac {\mathrm {D} {\boldsymbol {r}}'}{\mathrm {D} t}}-{\frac {\mathrm {D} {\boldsymbol {r}}}{\mathrm {D} t}}&({\boldsymbol {r}}'-{\boldsymbol {r}}=\left.\mathrm {d} {\boldsymbol {l}}\right|_{\boldsymbol {r}})\\&={\boldsymbol {v}}({\boldsymbol {r}}')-{\boldsymbol {v}}({\boldsymbol {r}})\\&=\left.\mathrm {d} {\boldsymbol {v}}\right|_{\boldsymbol {r}}\end{aligned}}}$ であることを用いた。 速度、圧力、密度、ポテンシャルは座標の連続一価関数であるので、 ${\displaystyle {\frac {\mathrm {D} {\it {\Gamma }}}{\mathrm {D} t}}=0}$ が証明された。

1. 巽友正『流体力学』（1982年 4月15日初版発行）培風館。ISBN 456302421X。

• 循環 (流体力学)
• ヘルムホルツの渦定理