Nablaの演算法則メモ
まず簡単な場合を考えたいため、基底を正規直交直線基底$ \mathsf{E}:=\{\pmb{e}_i\}_{i\in A\subseteq\N}に固定しておく
$ \pmb{\nabla}:=\sum_i\pmb{e}_i\frac{\partial}{\partial x_i}とする
$ \pmb{\nabla}\phi=\sum_i\pmb{e}_i\frac{\partial \phi}{\partial x_i}=\sum_i\frac{\partial \phi}{\partial x_i}\pmb{e}_i
$ \pmb{\nabla}(\phi+\psi)=\pmb{\nabla}\phi+\pmb{\nabla}\psi
$ \pmb{\nabla}(\phi\psi)=\sum\frac{\partial (\phi\psi)}{\partial x_i}\pmb{e}_i=\psi\pmb{\nabla}\phi+\phi\pmb{\nabla}\psi
$ \pmb{\nabla}\cdot\pmb{a}=\sum_i\pmb{e}_i\cdot\frac{\partial \pmb{a}}{\partial x_i}=\sum_i\frac{\partial [\pmb{a}]_i^\mathsf{E}}{\partial x_i}
$ \pmb{\nabla}\cdot(\pmb{a}+\pmb{b})=\pmb{\nabla}\cdot\pmb{a}+\pmb{\nabla}\cdot\pmb{b}
$ \pmb{\nabla}\pmb{a}=\sum\pmb{e}_i\frac{\partial \pmb{a}}{\partial x_i}=\sum\frac{\partial [\pmb{a}]_j^\mathsf{E}}{\partial x_i}\pmb{e}_i\pmb{e}_j
$ \mathrm{tr}(\pmb{\nabla}\pmb{a})=(\pmb{\nabla}\pmb{a}):\pmb{I}=\pmb{\nabla}\cdot\pmb{a}
$ \pmb{\nabla}(\pmb{a}+\pmb{b})=\sum\left(\pmb{e}_i\frac{\partial (\pmb{a}+\pmb{b})}{\partial x_i}\right)
$ =\pmb{\nabla} \pmb{a}+\pmb{\nabla} \pmb{b}
$ \pmb{\nabla} (\pmb{\nabla}\phi)=\sum\pmb{e}_i \frac{\partial \pmb{\nabla}\phi}{\partial x_i}
$ =\sum_{i,j}\frac{\partial^2 \phi}{\partial x_i\partial x_j}\pmb{e}_i \pmb{e}_j
曲線座標の場合、この展開が怪しくなるので注意
内側のnablaの基底に、外側のnablaの偏微分演算子がかかってくる
$ \mathrm{tr}(\pmb{\nabla} (\pmb{\nabla}\phi))=\pmb{\nabla}\cdot(\pmb{\nabla}\phi)となる
$ \because \mathrm{tr}(\pmb{\nabla} \pmb{a})=\pmb{\nabla}\cdot\pmb{a}
$ \pmb{\nabla} \pmb{\nabla}=\sum_{i,j}\pmb{e}_i \pmb{e}_j\frac{\partial^2}{\partial x_i\partial x_j}とするなら、
$ \pmb{\nabla} (\pmb{\nabla}\phi)=(\pmb{\nabla} \pmb{\nabla})\phi
と書ける(結合法則)
また$ \mathrm{tr}(\pmb{a} \pmb{b})=\pmb{a}\cdot\pmb{b}の類推も効く
$ \mathrm{tr}(\pmb{\nabla} (\pmb{\nabla}\phi))=\mathrm{tr}(\pmb{\nabla} \pmb{\nabla})\phi=\pmb{\nabla}^2\phi
$ \pmb{\nabla}\cdot(\pmb{\nabla}\phi)=\sum\frac{\partial^2 \phi}{{\partial x_i}^2}
$ \pmb{\nabla}\cdot\pmb{\nabla}=\sum_i\left|\pmb{e}_i\frac{\partial}{\partial x_i}\right|^2=\sum_i\frac{\partial^2}{{\partial x_i}^2}とするなら、
$ \pmb{\nabla}\cdot(\pmb{\nabla}\phi)= (\pmb{\nabla}\cdot\pmb{\nabla})\phi
$ = |\pmb{\nabla}|^2\phi
$ =: \Delta\phi
$ \pmb{\nabla} \pmb{T}=\sum\pmb{e}_i \frac{\partial \pmb{T}}{\partial x_i}
$ =\sum\frac{\partial [\pmb{T}]_{jk}^\mathsf{EE}}{\partial x_i}\pmb{e}_i \pmb{e}_j \pmb{e}_k
$ \pmb{\nabla} (\pmb{\nabla} \pmb{a})=\sum\frac{\partial^2 [\pmb{a}]_k^\mathcal{S}}{\partial x_i\partial x_j}\pmb{e}_i \pmb{e}_j \pmb{e}_k
$ =(\pmb{\nabla} \pmb{\nabla}) \pmb{a}
$ \pmb{\nabla}(\pmb{\nabla}\cdot\pmb{a})=\sum\frac{\partial^2 [\pmb{a}]_j^\mathsf{E}}{\partial x_i\partial x_j}\pmb{e}_i
$ \pmb{\nabla}\wedge\pmb{a}=\pmb{\nabla} \pmb{a}-\pmb{a} \pmb{\nabla}
$ \pmb{a} \pmb{\nabla}:=(\pmb{\nabla} \pmb{a})^\topとすると、
$ =\sum_{i\neq j}\left(\frac{\partial [\pmb{a}]_j^\mathsf{E}}{\partial x_i}-\frac{\partial [\pmb{a}]_i^\mathsf{E}}{\partial x_j}\right)\pmb{e}_i \pmb{e}_j
(左発散)$ \mathrm{div}(\pmb{T}):=\sum_j\frac{\partial \pmb{T}}{\partial x_j}\pmb{e}_j=\sum_{i,j}\frac{\partial [\pmb{T}]_{ij}^\mathsf{SS}}{\partial x_j}\pmb{e}_i
(右発散)$ \mathrm{div}_R(\pmb{T}):=\sum_i\frac{\partial \pmb{T}^\top}{\partial x_i}\pmb{e}_i=\sum_{i,j}\frac{\partial [\pmb{T}]_{ij}^\mathsf{SS}}{\partial x_i}\pmb{e}_j
定義より$ \mathrm{div}_R(\pmb{T})=\mathrm{div}(\pmb{T}^\top)
もし$ \pmb{a}\cdot\pmb{T}:=\pmb{T}^\top\pmb{a}とするなら、$ \mathrm{div}_R(\pmb{T})=\pmb{\nabla}\cdot\pmb{T}とできる
$ \mathrm{div}(\pmb{T})=(\pmb{\nabla}\otimes\pmb{T}):\pmb{I}
$ =\sum\frac{\partial [\pmb{T}]_{jj}^\mathsf{EE}}{\partial x_i}\pmb{e}_i
$ =\pmb{\nabla}\mathrm{tr}(\pmb{T})
これはおかしいtakker.icon
本来は$ \sum_{i,j}\frac{\partial [\pmb{T}]_{ij}^\mathsf{EE}}{\partial x_j}\pmb{e}_i (左発散)になるはず わかった。↑での勾配の定義が違ったんだ
このメモでは$ \pmb{\nabla}\otimes\bulletをtensorの勾配として用いている
これなら$ (\pmb{\nabla}\pmb{T}):\pmb{I}=\sum_{i,j}\frac{\partial [\pmb{T}]_{ij}^\mathsf{EE}}{\partial x_j}\pmb{e}_i となる
左発散の場合
$ \mathrm{div}(\pmb{a} \pmb{b})=\sum_{i,j}\frac{\partial [\pmb{a}]_i^\mathsf{E}[\pmb{b}]_j^\mathsf{E}}{\partial x_j}\pmb{e}_i
$ =\sum_{i,j}[\pmb{a}]_i^\mathsf{E}\frac{\partial[\pmb{b}]_j^\mathsf{E}}{\partial x_j}\pmb{e}_i+\sum_{i,j}[\pmb{b}]_j^\mathsf{E}\frac{\partial[\pmb{a}]_i^\mathsf{E}}{\partial x_j}\pmb{e}_i
$ =\pmb a(\pmb\nabla\cdot\pmb b)+\pmb b\cdot\pmb\nabla\pmb a
$ =(\pmb{\nabla}\cdot\pmb{b})\pmb{a}+(\pmb{b}\cdot\pmb{\nabla})\pmb{a}
別解
$ \mathrm{div}(\pmb{a} \pmb{b})=\sum_i(\pmb{\nabla}\cdot([\pmb{a}]^\mathsf{E}_i\pmb{b}))\pmb{e}_i
$ =\sum_i\left((\pmb{\nabla}[\pmb{a}]^\mathsf{E}_i)\cdot\pmb{b}+[\pmb{a}]^\mathsf{E}_i\pmb{\nabla}\cdot\pmb{b}\right)\pmb{e}_i
$ =\sum_i(\pmb{e}_i \pmb{\nabla}[\pmb{a}]^\mathsf{E}_i)\pmb{b}+(\pmb{\nabla}\cdot\pmb{b})\pmb{a}
$ =(\pmb{\nabla}\cdot\pmb{b})\pmb{a}+(\pmb{b}\cdot\pmb{\nabla})\pmb{a}
別解
$ \mathrm{div}(\pmb{a} \pmb{b})=\sum_i\frac{\partial \pmb{a} \pmb{b}}{\partial x_i}\pmb{e}_i
$ =\pmb{a} \sum_i\frac{\partial\pmb{b}}{\partial x_i}\pmb{e}_i+\sum_i\left(\frac{\partial\pmb{a}}{\partial x_i} \pmb{b}\right)\pmb{e}_i
$ =\left(\sum_i\pmb{e}_i\cdot\frac{\partial\pmb{b}}{\partial x_i}\right)\pmb{a}+\sum_i(\pmb{b}\cdot\pmb{e}_i)\frac{\partial\pmb{a}}{\partial x_i}
$ =(\pmb{\nabla}\cdot\pmb{b})\pmb{a}+(\pmb{b}\cdot\pmb{\nabla})\pmb{a}
右発散の場合
$ \mathrm{div}_R(\pmb{a}\otimes\pmb{b})=\mathrm{div}(\pmb{b}\otimes\pmb{a})
$ =(\pmb{\nabla}\cdot\pmb{a})\pmb{b}+(\pmb{a}\cdot\pmb{\nabla})\pmb{b}
$ \pmb{a}\cdot(\pmb{b}\otimes\pmb{c})=(\pmb{a}\cdot\pmb{b})\pmb{c}=(\pmb{b}\cdot\pmb{a})\pmb{c}より、$ \pmb{\nabla}\cdot(\pmb{a}\otimes\pmb{b})=(\pmb{\nabla}\cdot\pmb{a})\pmb{b}+(\pmb{a}\cdot\pmb{\nabla})\pmb{b}ともできる
(左発散)$ \mathrm{div}(\pmb{\nabla}\otimes\pmb{a})=\sum_{i,j}\left(\frac{\partial}{\partial x_j}\left(\pmb{e}_i\otimes\frac{\partial \pmb{a}}{\partial x_i}\right)\right)\pmb{e}_j
$ =\sum_{i,j}\frac{\partial}{\partial x_j}\left(\frac{\partial\pmb{a}}{\partial x_i}\cdot\pmb{e}_j\right)\pmb{e}_i
$ =\sum_{i,j}\frac{\partial^2[\pmb{a}]^\mathsf{E}_j}{\partial x_j\partial x_i}\pmb{e}_i
もし分母の偏微分を交換できるなら(普通は交換できる)、
$ =\pmb{\nabla}(\pmb{\nabla}\cdot\pmb{a})
$ \pmb{a}\cdot(\pmb{b}\otimes\pmb{c})^\top=(\pmb{b}\otimes\pmb{c})\pmb{a}=\pmb{b}(\pmb{a}\cdot\pmb{c})のanalogyからも直感的に導ける
(左発散)$ \mathrm{div}((\pmb{\nabla}\otimes\pmb{a})^\top)=\sum_{i,j}\frac{\partial}{\partial x_j}\frac{\partial (\pmb{a}\otimes\pmb{e}_i)\pmb{e}_j}{\partial x_i}
$ =\sum_{i,j}\frac{\partial}{\partial x_j}\frac{\partial (\pmb{e}_i\cdot\pmb{e}_j)\pmb{a}}{\partial x_i}
$ =\sum_i\frac{\partial^2\pmb{a}}{{\partial x_j}^2}
$ =\pmb{\nabla}^2\pmb{a}
$ \pmb{a}\cdot(\pmb{b}\otimes\pmb{c})=(\pmb{a}\cdot\pmb{b})\pmb{c}のanalogyで、$ \pmb{\nabla}\cdot(\pmb{\nabla}\otimes\pmb{a})=\pmb{\nabla}^2\pmb{a}と直ちに導出することもできる
$ \bm\nabla*(\bm a*\bm b)系
$ \bm\nabla*(\bm a*\bm b)=\sum_i\bm e_i*\frac{\partial (\bm a*\bm b)}{\partial x_i}
$ = (\bm\nabla*\bm a)*\bm b+\sum_i\bm e_i*\bm a*\frac{\partial\bm b}{\partial x_i}
$ *に交換法則が成立する演算子などを持ってくると、もう少しきれいになる
$ \bm{\nabla}\cdot(\phi\bm{a})=(\bm{\nabla}\phi)\cdot\bm{a}+\phi\bm{\nabla}\cdot\bm{a}
一般式から求まらない形
一般式が一般式じゃなかったtakker.icon
$ \bm{\nabla} (\phi\bm{a})=(\bm{\nabla}\phi) \bm{a}+\phi(\bm{\nabla} \bm{a})
$ \bm\nabla(\bm a\cdot\bm b)=(\bm\nabla \bm a)\cdot\bm b+(\bm\nabla \bm b)\cdot\bm a
$ =(\bm\nabla \bm a)\cdot\bm b+\bm a\cdot(\bm b\overleftarrow{\bm\nabla})
$ \bm{\nabla} (\bm{a} \bm{b})=(\bm{\nabla} \bm{a}) \bm{b}+\sum_i\bm e_i\bm a\frac{\partial\bm b}{\partial x_i}
第2項はこれ以上簡略化できない
$ \bm{\nabla}\cdot(\bm{T}\bm{a})=(\bm{\nabla}\cdot\bm{T})\bm{a}+\bm{T}^\top\cdot(\bm{\nabla}\bm{a})
$ \bm{\nabla}\cdot(\bm{T}\cdot\bm{A})=(\bm{\nabla}\cdot\bm{T})\cdot\bm{A}+\bm{T}:\bm\nabla\bm A
$ \bm Tが2階tensorのときしか成立しない
$ \bm{\nabla}\cdot(\bm a\cdot\bm T)=\sum_i\bm e_i\cdot\frac{\partial (\bm a\cdot\bm T)}{\partial x_i}
$ = \bm a\overleftarrow{\bm\nabla}:\bm T+\bm a\cdot\bm T\cdot \overleftarrow{\bm\nabla}
これも一般式では表現できない
一般化すると
References