Lagrangeの渦定理

from 渦度方程式

完全流体かつbarotropic流体な流体に保存力しか作用していない時、渦の不生不滅が成立する

すなわち、

ある時刻で$ \bm\Omega=\bm 0なら、ずっと$ \bm\Omega=\bm 0である

ある時刻で$ \bm\Omega\neq\bm 0なら、ずっと$ \bm\Omega\neq\bm0である

$ \bm\Omega：渦度

粘性のすくないサラサラの液体中で渦を作ったとき、渦がなかなか消えないのはこの定理に基づく

浴槽の水などで試してみるといいtakker.icon

導出

物質表示された完全流体の運動方程式$ \left(\tilde{\bm f}-\tilde\rho\frac{\partial^2\bm\phi}{{\partial t}^2}\right)\cdot\bm F=\bm\nabla\tilde Pから出発する

$ -\tilde\bm f+\bm\nabla\tilde P=\tilde\rho\bm\nabla Qとなる$ Qを定義する

Lagrangeの渦定理の前提条件から存在が保証される

$ \left(\tilde{\bm f}-\tilde\rho\frac{\partial^2\bm\phi}{{\partial t}^2}\right)\cdot\bm F=\bm\nabla\tilde P

$ \implies\ddot\bm\phi\cdot\bm F=\bm\nabla Q

$ \implies(\ddot\bm\phi\cdot\bm F)\times\overleftarrow\bm\nabla=\bm 0

$ \because∇⨯∇φ=0

scalar potentialはrotationを取ると消えてしまうということが、この定理の重要なポイント

$ \iff(\bm\nabla\ddot\bm\phi\cdot\bm F)^\top:{\Large\bm\epsilon}+\ddot\bm\phi\cdot(\bm F\times\overleftarrow\bm\nabla)=\bm 0

$ {\Large\bm\epsilon}：3階完全反対称tensor

$ \iff(\bm\nabla\ddot\bm\phi\cdot\bm F)^\top:{\Large\bm\epsilon}=\bm 0

$ \because∇⨯∇φ=0だから$ \bm F\times\overleftarrow\bm\nabla=\bm0

$ \iff\ddot\bm F^\top\cdot\bm F=\left(\ddot\bm F^\top\cdot\bm F\right)^\top

$ \because\bm A:{\Large\bm\epsilon}=\bm 0\iff\bm A=\bm A^\top

$ \iff\frac{\partial}{\partial t}(\dot\bm F^\top\cdot\bm F)-\dot\bm F^\top\cdot\dot\bm F=\frac{\partial}{\partial t}(\bm F^\top\cdot\dot\bm F)-\dot\bm F^\top\cdot\dot\bm F

$ \iff\frac{\partial}{\partial t}(\dot\bm F^\top\cdot\bm F-\bm F^\top\cdot\dot\bm F)=\bm0

$ \iff\frac{\partial}{\partial t}(\bm F^\top\cdot\left.\bm l\right|_{\bm x=\bm\phi(\bm X,t)}^\top\cdot\bm F-\bm F^\top\cdot\left.\bm l\right|_{\bm x=\bm\phi(\bm X,t)}\cdot\bm F)=\bm0

$ \bm l：速度勾配tensor

$ \iff\frac{\partial}{\partial t}(-2\bm F^\top\cdot\left.\bm w\right|_{\bm x=\bm\phi(\bm X,t)}^\top\cdot\bm F)=\bm0

$ \bm w：spin tensor

$ \iff\frac{\partial}{\partial t}(\bm F^\top\cdot\left.\bm w\right|_{\bm x=\bm\phi(\bm X,t)}\cdot\bm F)=\bm0

Lie時間微分と似ている？takker.icon

$ \implies\bm F^\top\cdot\left.\bm w\right|_{\bm x=\bm\phi(\bm X,t)}\cdot\bm F=\left.\left(\bm F^\top\cdot\left.\bm w\right|_{\bm x=\bm\phi(\bm X,t)}\cdot\bm F\right)\right|_{t=0}

$ =\left.\bm w\right|_{\bm x=\bm X,t=0}

$ \because\left.\bm F\right|_{t=0}=\bm I,\bm\phi(\bm X,0)=\bm X

$ \implies{\Large\bm\epsilon}:\left(\bm F^\top\cdot\left.\bm w\right|_{\bm x=\bm\phi(\bm X,t)}\cdot\bm F\right)=\left.\bm\Omega\right|_{\bm x=\bm X,t=0}

$ \bm\Omega={\Large\bm\epsilon}:\bm w=\bm\nabla\times\bm v：渦度

$ \implies{\Large\bm\epsilon}\vdots\left(\bm F^\top\cdot\left.\bm w\right|_{\bm x=\bm\phi(\bm X,t)}\cdot\bm F\right)\bm F^\top=\left.\bm\Omega\right|_{\bm x=\bm X,t=0}\cdot\bm F^\top

$ \iff \left.\bm\Omega\right|_{\bm x=\bm\phi(\bm X,t)}=\frac1J\bm F\cdot\left.\bm\Omega\right|_{\bm x=\bm X,t=0}

$ \because{\Large\bm\epsilon}\vdots(\bm A^\top\cdot\bm B\cdot\bm A)\bm A^\top=\epsilon_{ijk}A_{li}B_{lm}A_{mj}A_{nk}

$ = \epsilon_{ijk}\epsilon_{lmn}A_{0i}A_{1j}A_{2k}B_{lm}

$ \becauseLevi-Civita記号#6610c0c21280f00000e72ec2

$ = (\det\bm A){\Large\bm\epsilon}:\bm B

係数が面積変化率と酷似している点がめっちゃ気になるtakker.icon

以上より、以下が示された

ある時刻で$ \bm\Omega=\bm 0なら、ずっと$ \bm\Omega=\bm 0である

ある時刻で$ \bm\Omega\neq\bm 0なら、ずっと$ \bm\Omega\neq\bm0である

計算メモ

参照先の式b をこのページの記号で大ざっぱに書き直して、式展開がずれていないか確かめる

$ \dot F_{xx}F_{xy}-\dot F_{xy}F_{xx}+\dot F_{yx}F_{yy}-\dot F_{yy}F_{yx}=\dot F^\top_{xi}F_{iy}-F^\top_{xi}\dot F_{iy}