Q准则tecplot计算（Q2不变量）

##理论部分
应变率张量：$S_{ij}=\frac{1}{2}(\frac{\partial u_i}{\partial x_j}-\frac{\partial u_j}{\partial x_i})$
涡量：${\Omega }_{ij}=\frac{1}{2}(\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i})$
$Q$ 值：$Q=\frac{1}{2}(\Vert {\Omega }^2\Vert -\Vert S^2\Vert )$
简化到三维笛卡尔坐标下的$Q$值如下：
$$Q=-\frac{1}{2}\left({\left(\frac{\partial u}{\partial x}\right)}^2+{\left(\frac{\partial v}{\partial y}\right)}^2+{\left(\frac{\partial w}{\partial z}\right)}^2\right)-\frac{\partial u}{\partial y}\frac{\partial v}{\partial x}-\frac{\partial u}{\partial z}\frac{\partial w}{\partial x}-\frac{\partial v}{\partial z}\frac{\partial w}{\partial y}$$
其中$u$，$v$，$w$分别是$x$，$y$，$z$方向上的速度。

##tecplot实现

将tecplot实现代码粘贴如下：

{Q}=-1/2*((ddx({u}))**2+(ddy({v}))**2+(ddz({w}))**2)-ddy({u})*ddx({v})-ddz({u})*ddx({w})-ddz({v})*ddy({w})