张量分析初步和矢量恒等式

此文献给今天生日的跌跤dalao，祝哈耶普哥PhD~
——————————————

问题引入

整个国庆假期基本上都是在做作业，流体力学作业从中秋节那天做到今天，让人不由自主得怒P一图：

中秋节那天做了一道有意思的问题（实际上就是第一题科科），请教了同宿舍的凝聚态物理清本dalao后，稍微掌握了一点点张量分析的技巧，强大的符号艺术工具。问题是这样的，证明：$\frac{{d{\bf{v}}}}{{dt}} = \frac{{\partial {\bf{v}}}}{{\partial t}} + \left( {\nabla \times {\bf{v}}} \right) \times {\bf{v}} + \frac{1}{2}\nabla \left( {{{\bf{v}}^2}} \right)$，实际上由全微分公式可以容易得到$\frac{{d{\bf{v}}}}{{dt}} = \frac{{\partial {\bf{v}}}}{{\partial t}} + \left( {{\bf{v}} \cdot \nabla } \right){\bf{v}}$，所以实际上我们只需要证明：

$$\left( {{\bf{v}} \cdot \nabla } \right){\bf{v}} = \left( {\nabla \times {\bf{v}}} \right) \times {\bf{v}} + \frac{1}{2}\nabla \left( {{{\bf{v}}^2}} \right),$$
这是一个人纯矢量的问题了，这种小问题怎么可能难得了我，一看 $\left( {\nabla \times {\bf{v}}} \right) \times {\bf{v}}$，嘿嘿，二话不说我们先来个二重外积公式化简一下：
$$\left( {\nabla \times {\bf{v}}} \right) \times {\bf{v}} = \left( {{\bf{v}} \cdot \nabla } \right) \cdot {\bf{v}} - \left( {{\bf{v}} \cdot {\bf{v}}} \right) \cdot \nabla ,$$
好像好厉害的样子啊，不过 $\left( {{\bf{v}} \cdot {\bf{v}}} \right) \cdot \nabla$是什么东西？【待会下面我们会讲到这里二重外积公式是不适用的】，ε=(´ο｀*)))那从 $\nabla \left( {{{\bf{v}}^2}} \right)$找找突破口吧，
$$\nabla \left( {{{\bf{v}}^2}} \right) = \sum {\frac{{\partial {{\bf{v}}^2}}}{{\partial x}}} {\bf{i}} = \sum {\frac{{\partial {{\bf{v}}^2}}}{{\partial {\bf{v}}}} \cdot \frac{{\partial {\bf{v}}}}{{\partial x}}} {\bf{i}} = 2{\bf{v}} \cdot \sum {\frac{{\partial {\bf{v}}}}{{\partial x}}{\bf{i}} = 2{\bf{v}} \cdot } \nabla {\bf{v}},$$
样子是有那么点像了，但是 $\nabla {\bf{v}}$的话，这怎么处理，矢量的梯度的话，岂不是变得更复杂了。

矢量分析初步

• （符号艺术1）爱因斯坦求和约定：对于一个表达式中的某项，如果一个指标重复出现两次，那么就将该项在该指标下遍历求和，该指标成为哑变量。给一个例子，如$\sum_{i=1}^3a_ib_i=a_ib_i$，这是一种简化的写法。
• （符号艺术2）克罗内克符号：用${\delta _{ij}}$表示，如果$i=j$，那么${\delta _{ij}}=1$，否则为0，克罗内克符号起到了一种神奇的换下标的功能，比如：${\delta _{ij}}{a_{jk}} = {a_{ik}}$，从它的爱因斯坦求和约定展开后轻松看出。【简单地说就是如果和$\delta$相乘的量有相同的，比如这里的$j$，那么将重复的量$j$换为$\delta$的另一个下标】
• （符号艺术3）置换符号：在行列式学习的时候，我们知道有：
  $$\det = \sum\limits_{{j_1},{j_2},...,{j_n}} {{{\left( { - 1} \right)}^{\tau \left( {{j_1},{j_2},...,{j_n}} \right)}}{a_{1{j_1}}}{a_{2{j_2}...}}{a_{n{j_n}}}} ,$$
  其中的${{{\left( { - 1} \right)}^{\tau \left( {{j_1},{j_2},...,{j_n}} \right)}}}$我们就称为置换符号，由于矢量运算一般是三维的，我们常用$e_{rst}$来表示，我们可以用$e_{rst}=\frac{1}{2}\cdot (r-s)(s-t)(t-r)$来定义，也可以如下标准定义：
  $${e_{rst}} = \left\{ \begin{array}{l} 1,逆序对为偶数个\\ - 1,逆序对为奇数个\\ 0,有相同的下标（为了方便爱因斯坦求和约定） \end{array} \right.，$$
  有了这三个表述之后，我们可以稍微搞点事情，考虑三阶行列式：
  $$\left| {\begin{array}{*{20}{c}} {{a_{11}}}&{{a_{12}}}&{{a_{13}}}\\ {{a_{21}}}&{{a_{22}}}&{{a_{23}}}\\ {{a_{31}}}&{{a_{32}}}&{{a_{33}}} \end{array}} \right| = {e_{ijk}}{a_{i1}}{a_{j2}}{a_{k3}} = {e_{ijk}}{a_{1i}}{a_{2j}}{a_{3k}},$$
  简单地可以证明（列换序）有：
  $$\left| {\begin{array}{*{20}{c}} {{a_{1r}}}&{{a_{1s}}}&{{a_{1t}}}\\ {{a_{2r}}}&{{a_{2s}}}&{{a_{2t}}}\\ {{a_{3r}}}&{{a_{3s}}}&{{a_{3t}}} \end{array}} \right| = {e_{rst}}{e_{ijk}}{a_{i1}}{a_{j2}}{a_{k3}},$$
  进一步地，行换序得到一般形式有：
  $$\left| {\begin{array}{*{20}{c}} {{a_{or}}}&{{a_{os}}}&{{a_{ot}}}\\ {{a_{pr}}}&{{a_{ps}}}&{{a_{pt}}}\\ {{a_{qr}}}&{{a_{qs}}}&{{a_{qt}}} \end{array}} \right| = {e_{opq}}{e_{rst}}{e_{ijk}}{a_{i1}}{a_{j2}}{a_{k3}},$$
  现在我们根据上式导出一个关键公式：
  $$\left| {\begin{array}{*{20}{c}} {{\delta _{or}}}&{{\delta _{os}}}&{{\delta _{ot}}}\\ {{\delta _{pr}}}&{{\delta _{ps}}}&{{\delta _{pt}}}\\ {{\delta _{qr}}}&{{\delta _{qs}}}&{{\delta _{qt}}} \end{array}} \right| = {e_{opq}}{e_{rst}}{e_{ijk}}{\delta _{i1}}{\delta _{j2}}{\delta _{k3}} = {e_{opq}}{e_{rst}}{e_{123}} = {e_{opq}}{e_{rst}},$$
  取$o=r$立马有如下关键公式（联系克罗内克符号和置换符号）：
  $${e_{rpq}}{e_{rst}} = \left| {\begin{array}{*{20}{c}} {{\delta _{rr}}}&{{\delta _{rs}}}&{{\delta _{rt}}}\\ {{\delta _{pr}}}&{{\delta _{ps}}}&{{\delta _{pt}}}\\ {{\delta _{qr}}}&{{\delta _{qs}}}&{{\delta _{qt}}} \end{array}} \right| = \left| {\begin{array}{*{20}{c}} 1&0&0\\ {{\delta _{pr}}}&{{\delta _{ps}}}&{{\delta _{pt}}}\\ {{\delta _{qr}}}&{{\delta _{qs}}}&{{\delta _{qt}}} \end{array}} \right| = {\delta _{ps}}{\delta _{qt}} - {\delta _{qs}}{\delta _{pt}},$$

对于这些初步只是做个总结：

• 克罗内克符号和置换符号的关系：${e_{rpq}}{e_{rst}} = {\delta _{ps}}{\delta _{qt}} - {\delta _{qs}}{\delta _{pt}}$
• 克罗内克符号的换下标作用： ${\delta _{ij}}{a_{jk}} = {a_{ik}}$

矢量恒等式

下面我们回到$\left( {{\bf{v}} \cdot \nabla } \right){\bf{v}} = \left( {\nabla \times {\bf{v}}} \right) \times {\bf{v}} + \frac{1}{2}\nabla \left( {{{\bf{v}}^2}} \right)$的证明，

$$\begin{array}{l} \left( {\nabla \times {\bf{v}}} \right) \times {\bf{v}}\\ = \left( {{e_{klm}} \cdot {{\bf{i}}_k} \cdot {\partial _{2l}}{v_{3m}}} \right) \times {\bf{v}}\\ = {e_{klm}} \cdot {{\bf{i}}_k} \cdot {\partial _{2l}}{v_{3m}} \cdot {e_{abc}} \cdot {{\bf{i}}_a} \cdot {v_{2b}}【只有{\bf{i}_a=\bf{i}_k}时，点积才为1，否则为0】\\ = {e_{klm}}{e_{kbc}} \cdot {v_{2b}} \cdot {\partial _{2l}}{v_{3m}}【由上解释有e_{abc}只剩下e_{kbc}】\\ = \left( {{\delta _{lb}}{\delta _{mc}} - {\delta _{mb}}{\delta _{lc}}} \right){v_{2b}} \cdot {\partial _{2l}}{v_{3m}}\\ = {\delta _{lb}}{\delta _{mc}}{v_{2b}} \cdot {\partial _{2l}}{v_{3m}} - {\delta _{mb}}{\delta _{lc}}{v_{2b}} \cdot {\partial _{2l}}{v_{3m}}【运用了克罗内克符号的换下标功能】\\ = {v_{2b}} \cdot {\partial _{2b}}{v_{3c}} - {v_{2b}} \cdot {\partial _{2c}}{v_{3b}}\\ = \left( {{\bf{v}} \cdot \nabla } \right){\bf{v}} - \frac{1}{2}{\partial _c}v_b^2\\ = \left( {{\bf{v}} \cdot \nabla } \right){\bf{v}} - \frac{1}{2}\nabla \left( {{{\bf{v}}^2}} \right) \end{array}$$

于是乎命题得证，有人可能还不过瘾，那么我们顺便来秒杀二重外积公式吧，即证明：$\left( {{\bf{A}} \times {\bf{B}}} \right) \times {\bf{C}} = \left( {{\bf{A}} \cdot {\bf{C}}} \right) \cdot {\bf{B}} - \left( {{\bf{B}} \cdot {\bf{C}}} \right) \cdot {\bf{A}}$,同样利用张量分析：

$$\begin{array}{l} \left( {{\bf{A}} \times {\bf{B}}} \right) \times {\bf{C}}\\ = \left( {{e_{klm}} \cdot {{\bf{i}}_k} \cdot {A_l}{B_m}} \right){e_{tij}} \cdot {{\bf{i}}_t} \cdot {C_i}\\ = \left( {{\delta _{il}}{\delta _{jm}} - {\delta _{jl}}{\delta _{im}}} \right){A_l}{B_m}{C_i}.........①\\ = {A_i}{B_j}{C_i} - {A_j}{B_i}{C_i}................②\\ = \left( {{\bf{A}} \cdot {\bf{C}}} \right) \cdot {\bf{B}} - \left( {{\bf{B}} \cdot {\bf{C}}} \right) \cdot {\bf{A}} \end{array}$$

上面推导中①到②是需要$A,B,C$之间可交换的，而最开始$\left( {\nabla \times {\bf{v}}} \right) \times {\bf{v}} = \left( {{\bf{v}} \cdot \nabla } \right) \cdot {\bf{v}} - \left( {{\bf{v}} \cdot {\bf{v}}} \right) \cdot \nabla$这个操作并不科学，这是因为$\frac{\partial }{{\partial x}}a \cdot b \ne \frac{\partial }{{\partial x}}b \cdot a$，其中$\frac{\partial }{{\partial x}}a \equiv \frac{{\partial a}}{{\partial x}}$。
这告诉我们公式不是乱套用的，尤其是整天给我逼逼CNN的大佬们，不懂其数学原理是后患无穷的。