Euler方程的涡度方程

Euler方程的涡度方程

${\mathbb{R}}^3$上的Euler方程如下

$$\{\begin{array}{rl}& {\partial }_{\mathrm{t}}u+{\mathrm{d}\mathrm{i}\mathrm{v}}_{\mathrm{x}}(\rho u\otimes u)+{\nabla }_{\mathrm{x}}p(\rho )=0\\ & {\partial }_{\mathrm{t}}\rho +{\mathrm{d}\mathrm{i}\mathrm{v}}_{\mathrm{x}}(\rho u)=0\end{array}$$

其中$p=\frac{1}{\gamma }{\rho }^{\gamma }$, $\gamma >1$. 定义涡度为$\omega ={\mathrm{c}\mathrm{u}\mathrm{r}\mathrm{l}}_{\mathrm{x}}u$, 以及定义单位质量的涡度为$\zeta =\omega /\rho$. 现在我们来导出$\zeta$满足的方程. 在动量守恒方程两边取旋度, 我们有

0 = ε i j k ∂ j ( ∂ t u k + u l ∂ l u k + ∂ k ( ρ γ − 1 / γ ) ) = ∂ t ω i + ε i j k ∂ j u l ∂ l u k + ε i j k u l ∂ l ∂ j u k = ∂ t ( ρ ζ i ) + u l