ベクトル解析の公式

POINT

  • ベクトル解析の公式と,その導出方法の一覧.
  • 行列計算も統一的に理解できる.

ベクトル解析の公式と,その導出方法を一覧にまとめました.

力学・電磁気学・流体力学などを学ぶ上で,これらの計算はとても重要です.計算練習をして,すぐに公式を導出できるようになりましょう!

【関連記事】

前提知識と記法

微分演算子の略号を使って計算すると,負担を減らすことができます.

微分演算子

この記事では簡単のため,微分演算子を
\begin{aligned}
\partial_i:=\frac{\partial}{\partial x^i}
\end{aligned}
と書くことにします($x_1=x$, $x_2=y$, $x_3=z$).

※ 一般的に使われる記法です.分数を使わず,一行で書くことが出来るので便利です.



例えば,

  • $\partial_i\phi$は,「$\boldsymbol{\nabla\,}\phi$の$i$成分」という意味
  • $\partial_2 f$は,$\dfrac{\partial f}{\partial y}$という意味 ($\partial_y f$とかくこともあります)

です.

ベクトル演算の成分表記

以下では,ダブっている添字については和を取ることにします(Einsteinの規約).すると,ベクトル演算の成分表記が次のように簡単に表せます.

名称 ベクトル表記 成分表記
内積 $\boldsymbol{A}\cdot\boldsymbol{B}$ $A_i B_i$
外積 $\boldsymbol{A}\times\boldsymbol{B}$ $\epsilon_{ijk}A_j B_k$
勾配 / Gradient $\mathrm{grad\,}\phi=\boldsymbol{\nabla}\phi$ $\partial_i\phi$
回転 / Rotation / Curl
\begin{aligned}\mathrm{rot\,}\boldsymbol{A}&=\mathrm{curl\,}\boldsymbol{A}\\&=\boldsymbol{\nabla}\times\boldsymbol{A}\end{aligned}
$\epsilon_{ijk}\partial_jA_k$
発散 / Divergence $\mathrm{div\,}\boldsymbol{A}=\boldsymbol{\nabla}\cdot\boldsymbol{A}$ $\partial_iA_i$
ラプラシアン / Lapracian
\begin{aligned} \Delta \phi&=\mathrm{div\,grad\,}\phi\\&=\boldsymbol{\nabla\,}\cdot\boldsymbol{\nabla\,}\phi \end{aligned}
$\partial_i\partial_i \phi$

行列の成分表記

一般的には,ベクトル演算は行列やテンソルの演算に含まれます.そのような場合には,計算結果が行列式の形にまとめられる場合があります(例:積分の座標変換でヤコビアンを計算する場合).そこで,行列式がどのように表されるのか,知っておく必要があります:
$n\times n$行列$A=(a_{i,j})$の行列式は,
\begin{aligned}
\mathrm{det\,}A
&=\epsilon_{i_1i_2\cdots i_n} a_{1i_1}a_{2i_2}\cdots a_{ni_n}\\
&=\epsilon_{i_1i_2\cdots i_n} a_{i_11}a_{i_22}\cdots a_{i_nn}
\end{aligned}
と書くことができます.

この式は,$\epsilon_{i_1i_2\cdots i_n}$の反対称性から

行列式
\begin{aligned}
\epsilon_{j_1j_2\cdots j_n}\mathrm{det\,}A
&=\epsilon_{i_1i_2\cdots i_n} a_{j_1 i_1}a_{j_2 i_2}\cdots a_{j_n i_n}\\
&=\epsilon_{i_1i_2\cdots i_n} a_{i_1 j_1}a_{i_2 j_2}\cdots a_{i_n j_n}
\end{aligned}
と一般化することができます.

公式の導出手順

ベクトル解析のどんな公式も,以下の手順で導けます!
手順

【一覧】公式と導出

上の手順に従い,公式を導出します.

積の微分:$\boldsymbol{\nabla}\cdot \left( \psi \boldsymbol{A}\right) $$= \boldsymbol{\nabla}\psi \cdot \boldsymbol{A} +\psi\left( \boldsymbol{\nabla}\cdot\boldsymbol{A}\right)$

積の微分のベクトル解析版です.成分表記すれば,通常の積の微分の公式が使えます.
\begin{aligned}
\boldsymbol{\nabla}\cdot \left( \psi \boldsymbol{A}\right)
&=\partial_i \left( \psi A_i\right) \\
&=\left(\partial_i \psi \right) A_i + \psi \left(\partial_i A_i \right) \\
&= \boldsymbol{\nabla}\psi \cdot \boldsymbol{A}
+\psi\left( \boldsymbol{\nabla}\cdot\boldsymbol{A}\right)
\end{aligned}

スカラー三重積:$\boldsymbol{a}\cdot (\boldsymbol{b}\times \boldsymbol{c})$$=(\boldsymbol{a}\times \boldsymbol{b})\cdot \boldsymbol{c}$

ベクトルの外積が
\begin{aligned}
(\boldsymbol{a}\times\boldsymbol{b})_i=\epsilon_{ijk}a_j b_k
\end{aligned}
で表わされることから,
\begin{aligned}
\boldsymbol{a}\cdot (\boldsymbol{b}\times \boldsymbol{c})
&=a_i (\epsilon_{ijk} b_j c_k)\\
&=\color{red}{\epsilon_{ijk} a_i b_j c_k}\\
&=(\epsilon_{kij} a_i b_j ) c_k\\
&=(\boldsymbol{a}\times \boldsymbol{b})\cdot \boldsymbol{c}
\end{aligned}

ベクトル三重積:$\boldsymbol{a}\times (\boldsymbol{b}\times \boldsymbol{c})$$=(\boldsymbol{a}\cdot \boldsymbol{c})\boldsymbol{b} -(\boldsymbol{a}\cdot \boldsymbol{b})\boldsymbol{c}$

縮約公式を使うと,
\begin{aligned}
[\boldsymbol{a}\times (\boldsymbol{b}\times \boldsymbol{c})]_i
&=\epsilon_{ijk}a_j(\epsilon_{klm}b_l c_m) \\
&=(\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl}) a_j b_l c_m \\
&= (\boldsymbol{a}\cdot \boldsymbol{c})b_i -(\boldsymbol{a}\cdot \boldsymbol{b})c_i
\end{aligned}

$\mathrm{rot\,}\mathrm{grad\,}=0$, $\mathrm{div\,}\mathrm{rot\,}=0$

微分の順序を交換できるとき*1 以下の重要な性質が成り立ちます:
\begin{aligned}
\epsilon_{ijk}\partial_j\partial_k=0
\end{aligned}

この式は,反対称性から簡単に導けます:
【証明】

\begin{aligned}
\epsilon_{ijk}\partial_j\partial_k=\epsilon_{ijk}\partial_k\partial_j.
\end{aligned}
ここで$j\leftrightarrow k$とし,$\epsilon_{ijk}$の反対称性を用いると
\begin{aligned}
\epsilon_{ijk}\partial_j\partial_k
=\epsilon_{ikj}\partial_j\partial_k
=-\epsilon_{ijk}\partial_j\partial_k.
\end{aligned}
であるから,上の式が導かれた.//

これを利用すると,以下の公式が導けます:

\begin{aligned}
&(\mathrm{rot\,}\mathrm{grad\,}\phi)_i
=\epsilon_{ijk}\partial_j\partial_k\phi
=0\\
&\mathrm{div\,}\mathrm{rot\,}\boldsymbol{A}
=\partial_i\epsilon_{ijk}\partial_j A_k
=0
\end{aligned}

$\mathrm{rot\,}\mathrm{rot\,}\boldsymbol{A}$$=\mathrm{grad\,}\mathrm{div\,}\boldsymbol{A}-\boldsymbol{\Delta\,A}$

\begin{aligned}
(\mathrm{rot\,}\mathrm{rot\,}\boldsymbol{A})_i
&=\epsilon_{ijk}\partial_j(\epsilon_{klm}\partial_l A_m) \\
&=(\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl})\partial_j\partial_l A_m \\
&=\partial_i\partial_j A_j-\partial_j\partial_jA_i \\
&=[\mathrm{grad\,}\mathrm{div\,}\boldsymbol{A}-\boldsymbol{\Delta\,A}]_i
\end{aligned}

ここで,右辺の$\boldsymbol{\Delta\,A}$は「ベクトルラプラシアン」です.

$\mathrm{div\,}(\boldsymbol{A}\times\boldsymbol{B})$$=\boldsymbol{B}\cdot\mathrm{rot\,}\boldsymbol{A}-\boldsymbol{A}\cdot\mathrm{rot\,}\boldsymbol{B}$

\begin{aligned}
\mathrm{div\,}(\boldsymbol{A}\times\boldsymbol{B})
&=\partial_i(\epsilon_{ijk}A_jB_k) \\
&=\epsilon_{ijk}(\partial_iA_j)B_k+\epsilon_{ijk}A_j(\partial_iB_k) \\
&=B_k\epsilon_{kij}\partial_iA_j - A_j\epsilon_{jik}\partial_iB_k \\
&=\boldsymbol{B}\cdot\mathrm{rot\,}\boldsymbol{A}-\boldsymbol{A}\cdot\mathrm{rot\,}\boldsymbol{B}
\end{aligned}

$\boldsymbol{\nabla} (\boldsymbol{A}\cdot\boldsymbol{B})$$=(\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{B}+(\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A}$$+\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B}+\boldsymbol{B}\times\mathrm{rot\,}\boldsymbol{A}$

\begin{aligned}
(\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B})_i
&=\epsilon_{ijk} A_j (\mathrm{rot\,}\boldsymbol{B})_k \\
&=\epsilon_{ijk} A_j \cdot \epsilon_{klm}\partial_l B_m \\
&=\epsilon_{kij}\epsilon_{klm} A_j \partial_l B_m \\
&=(\delta_{il}\delta_{jm} - \delta_{im}\delta_{jl}) A_j \partial_l B_m \\
&=A_j \partial_i B_j - A_j \partial_j B_i
\end{aligned}
に着目する.$\boldsymbol{A}$と$\boldsymbol{B}$を入れ替えたものを足し合わせると
\begin{aligned}
&(\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B}+\boldsymbol{B}\times\mathrm{rot\,}\boldsymbol{A})_i \\
&=(A_j \partial_i B_j+B_j \partial_i A_j) - (A_j \partial_j B_i+B_j \partial_j A_i) \\
&=\partial_i(A_j B_j) - (A_j \partial_j B_i) - (B_j \partial_j A_i) \\
&=[\boldsymbol{\nabla} (\boldsymbol{A}\cdot\boldsymbol{B}) - (\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{B} - (\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A} ]_i
\end{aligned}

$\mathrm{rot\,}(\boldsymbol{A}\times \boldsymbol{B})$$=\boldsymbol{A}\mathrm{\,div\,}\boldsymbol{B}-\boldsymbol{B}\mathrm{\,div\,}\boldsymbol{A}$$+(\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A}-(\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{B}$

\begin{aligned}
&[\mathrm{rot\,}(\boldsymbol{A}\times \boldsymbol{B})]_i \\
&=\epsilon_{ijk} \partial_j (\boldsymbol{A}\times \boldsymbol{B})_k \\
&=\epsilon_{ijk} \partial_j (\epsilon_{klm} A_l B_m) \\
&=(\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl}) [(\partial_j A_l) B_m + A_l (\partial_j B_m)]\\
&=[(\partial_j A_i) B_j + A_i (\partial_j B_j)] - [(\partial_j A_j) B_i + A_j (\partial_j B_i)] \\
&=[(\boldsymbol{B}\cdot \boldsymbol{\nabla})\boldsymbol{A} + \boldsymbol{A}\mathrm{\,div\,}\boldsymbol{B}
- \boldsymbol{B}\mathrm{\,div\,}\boldsymbol{A} - (\boldsymbol{A}\cdot \boldsymbol{\nabla})\boldsymbol{B}]_i
\end{aligned}

$(\boldsymbol{B}\cdot\boldsymbol{\nabla})\boldsymbol{A}$$=\dfrac{1}{2}\bigl[\boldsymbol{\nabla} (\boldsymbol{A}\cdot\boldsymbol{B}) +\mathrm{rot\,}(\boldsymbol{A}\times \boldsymbol{B}) $$-\boldsymbol{A}\mathrm{\,div\,}\boldsymbol{B}+\boldsymbol{B}\mathrm{\,div\,}\boldsymbol{A}$$-\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{B}-\boldsymbol{B}\times\mathrm{rot\,}\boldsymbol{A}\bigr]$

上で求めた2つの関係式を辺々加えることで示される.

$(\boldsymbol{A}\cdot\boldsymbol{\nabla})\boldsymbol{A}$$=\boldsymbol{\nabla} (|\boldsymbol{A}|^2/2)$$-\boldsymbol{A}\times\mathrm{rot\,}\boldsymbol{A}$

上で$\boldsymbol{A}=\boldsymbol{B}$とすれば導かれる.

この関係式は,流体力学でよく使われます:物質微分の意味と関係式(流体力学) - Notes_JP

*1:フツーできるときしか考えないので,あまり気にしないように.