ベクトル場における発散(div)と回転(rot)の定義・意味

ベクトル解析の中でも非常に重要な,ベクトル場の発散(div)・回転(rot)について定義と物理的な意味を説明します。 微分演算子であるナブラ演算子を導入します。

ベクトル場とは

あるベクトル A\boldsymbol{A} が空間上の各点 (x,y,z)(x,y,z) で定義されているとき,A(x,y,z)\boldsymbol{A}(x,y,z)ベクトル場といいます。

たとえば,流体の速度や磁石の周りの磁力線がベクトル場です。 ベクトル場

空間微分演算子

空間微分演算子 \nabla(ナブラ)は次のように定義されます。

空間微分演算子の定義

(x,y,z)(x,y,z) 空間上で空間微分演算子 \nabla は次のように定義される。 =(x,y,z). \nabla = \left(\dfrac{\partial{}}{\partial{x}},\dfrac{\partial{}}{\partial{y}},\dfrac{\partial{}}{\partial{z}}\right).

数学的に何かに作用するものを演算子とよびます。 たとえば xx\partial_x \coloneqq \dfrac{\partial}{\partial x} も微分演算子で,関数 f(x,y,z)f(x,y,z) に作用させると xf=fx\partial_xf=\dfrac{\partial f}{\partial x} となるといった具合です(ふつう微分演算子は右へ作用します)。微分演算子に限らず,a[f(x)]af(x)a[f(x)] \coloneqq af(x) と定義したときの a[]a[\cdot] も演算子と言えるでしょう。 演算子の中には複数の対象を引数にとるものもあります(たとえば二項係数 nCr{}_n\mathrm{C}_r は二項演算子です)。

空間微分演算子 \nabla がスカラー場 f(x,y,z)f(x,y,z) に作用したものは勾配(gradient)といい, f(x,y,z)=(fx,fy,fz)\nabla f(x,y,z) = \left(\dfrac{\partial{f}}{\partial{x}},\dfrac{\partial{f}}{\partial{y}},\dfrac{\partial{f}}{\partial{z}}\right) となります。→勾配ベクトルの意味と例題

ベクトル解析の発散(div)の定義

ベクトル場の発散(divergence)の定義は以下の通りです。

発散(div)の定義

ベクトル場 A(x,y,z)=(Ax(x,y,z),Ay(x,y,z),Az(x,y,z))\boldsymbol{A}(x,y,z)=(A_x(x,y,z),A_y(x,y,z),A_z(x,y,z)) の発散は次のように定義される。 divAA=Axx+Ayy+Azz. \mathrm{div}\,\boldsymbol{A} \coloneqq \nabla\cdot\boldsymbol{A} =\dfrac{\partial{A_x}}{\partial{x}}+\dfrac{\partial{A_y}}{\partial{y}}+\dfrac{\partial{A_z}}{\partial{z}}.

ここで登場した div\mathrm{div} も演算子で,右にあるベクトル場に作用します。

たとえば,ベクトル場 A(x,y,z)=(xy,yz,zx)\boldsymbol{A}(x,y,z)=(xy,yz,zx) の発散は, divA(x,y,z)=x(xy)+y(yz)+z(zx)=y+z+x\begin{aligned} \mathrm{div}\,\boldsymbol{A}(x,y,z) &=\dfrac{\partial{}}{\partial{x}}(xy)+\dfrac{\partial{}}{\partial{y}}(yz)+\dfrac{\partial{}}{\partial{z}}(zx)\\ &=y+z+x \end{aligned} となります。

ベクトル解析の発散(div)の意味

発散の物理的な意味を考えてみましょう。 ベクトル場の例として流体の流量 V(x,y,z)\boldsymbol{V}(x,y,z) を考えてみます。点 (x,y,z)(x,y,z) から V\boldsymbol{V} の方向に単位面積当たり V=VV=|\boldsymbol{V}| の量だけ流体が流れ出ているイメージです。

ある点 (x,y,z)(x,y,z) での流体の増減を調べることにします。 点での増減を直接調べるのは難しいので,点 (x,y,z)(x,y,z) を中心とする下図のような非常に小さな直方体を考えましょう(Δx,Δy,Δz1\varDelta x,\varDelta y,\varDelta z\ll1)。後々 Δx,Δy,Δz0\varDelta x,\varDelta y,\varDelta z\to0 の極限を取ることで点 (x,y,z)(x,y,z) での増減を調べます。

辺の長さが dx, dy, dz の直方体 x 方向に垂直な面の流量変化

まず,xx 方向に垂直な面の流入・流出を考えます。 点 (x±Δx/2,y,z)(x\pm\varDelta x/2,y,z) における xx 方向の(単位面積あたりの)流量は Vx(x±Δx/2,y,z)V_x(x\pm\varDelta x/2,y,z) なので, この面を通過する流量を Vx(x±Δx/2,y,z)ΔyΔzV_x(x\pm\varDelta x/2,y,z)\varDelta y\varDelta z と近似してやります。 yΔy/2y+Δy/2zΔz/2z+Δz/2Vx(x±Δx2,y,z)dydzVx(x±Δx2,y,z)ΔyΔz \begin{aligned} &\quad\int_{y-\varDelta y/2}^{y+\varDelta y/2}\int_{z-\varDelta z/2}^{z+\varDelta z/2} V_x\left(x\pm\dfrac{\varDelta x}{2},y,z\right)\,dy\,dz\\ &\simeq V_x\left(x\pm\dfrac{\varDelta x}{2},y,z\right)\varDelta y\varDelta z \end{aligned} Δy,Δz\varDelta y,\varDelta z は極めて小さいので,面上での VxV_x を一点の値として近似しています。

これを Taylor 展開して Δx\varDelta x の1次まで近似すると, Vx(x±Δx2,y,z)=Vx(x,y,z)±Vx(x,y,z)xΔx2+122Vx(x,y,z)x2(Δx2)2+Vx(x,y,z)±Vx(x,y,z)xΔx2\begin{aligned} &\quad V_x\left(x\pm\dfrac{\varDelta x}{2},y,z\right) \\ &=V_x(x,y,z)\pm\dfrac{\partial{V_x(x,y,z)}}{\partial{x}}\frac{\varDelta x}{2}+\dfrac{1}{2}\dfrac{\partial^2{V_x(x,y,z)}}{\partial{x}^2}\left(\frac{\varDelta x}{2}\right)^2+\cdots\\ &\simeq V_x(x,y,z)\pm\dfrac{\partial{V_x(x,y,z)}}{\partial{x}}\frac{\varDelta x}{2} \end{aligned} となります(Taylor 展開についてはテイラーの定理の例と証明を参照)。 正味の流量変化は, [Vx(x+Δx2,y,z)Vx(xΔx2,y,z)]ΔyΔz=[(Vx+VxxΔx2)(VxVxxΔx2)]ΔyΔz=VxxΔxΔyΔz \begin{aligned} &\quad\left[ V_x \left(x+ \dfrac{\varDelta x}{2},y,z\right) -V_x\left(x-\dfrac{\varDelta x}{2},y,z\right)\right] \varDelta y\varDelta z \\ &=\left[\left(V_x+\frac{\partial V_x}{\partial x}\frac{\varDelta x}{2}\right)-\left(V_x-\frac{\partial V_x}{\partial x}\frac{\varDelta x}{2}\right)\right]\varDelta y\varDelta z\\ &=\dfrac{\partial{V_x}}{\partial{x}}\varDelta x\varDelta y\varDelta z \end{aligned} となります。同様の計算を yy 方向に垂直な2面,zz 方向に垂直な2面にも適応すると,正味の流量変化はそれぞれ VyyΔxΔyΔz,VzzΔxΔyΔz\dfrac{\partial{V_y}}{\partial{y}}\varDelta x\varDelta y\varDelta z, \dfrac{\partial{V_z}}{\partial{z}}\varDelta x\varDelta y\varDelta z となるので,3方向の正味の流量変化は, (Vxx+Vyy+Vzz)ΔxΔyΔz \left( \dfrac{\partial{V_x}}{\partial{x}}+\dfrac{\partial{V_y}}{\partial{y}}+\dfrac{\partial{V_z}}{\partial{z}} \right)\varDelta x\varDelta y\varDelta z となります。V\boldsymbol{V} の発散の因子 V\nabla\cdot\boldsymbol{V} が出てきたことに気づいたと思います。V\nabla\cdot\boldsymbol{V} を用いて正味の流入変化を書き換えると VΔxΔyΔz\nabla\cdot\boldsymbol{V}\varDelta x\varDelta y\varDelta z と表せます。

つまり,発散 V\nabla\cdot\boldsymbol{V} はその点における単位体積当たりの流量変化を表しているとわかります。

V>0\nabla\cdot\boldsymbol{V}>0 の点では流体が湧き出しています(湧き出し)。 逆に V<0\nabla\cdot\boldsymbol{V}<0 の点では流体を吸い込んでいます(吸い込み)。 V=0\nabla\cdot\boldsymbol{V}=0 の点では湧き出しも吸い込みもないことになります。

ベクトル解析の回転(rot)の定義

ベクトル場の回転(rotation)の定義は以下の通りです。

回転(rot)の定義

ベクトル場 A(x,y,z)=(Ax(x,y,z),Ay(x,y,z),Az(x,y,z))\boldsymbol{A}(x,y,z)=(A_x(x,y,z),A_y(x,y,z),A_z(x,y,z)) の回転は次のように定義される。 rotA×A=(AzyAyz,AxzAzx,AyxAxy) \mathrm{rot}\,\boldsymbol{A} \coloneqq \nabla\times\boldsymbol{A} =\left(\dfrac{\partial{A_z}}{\partial{y}}-\dfrac{\partial{A_y}}{\partial{z}},\dfrac{\partial{A_x}}{\partial{z}}-\dfrac{\partial{A_z}}{\partial{x}},\dfrac{\partial{A_y}}{\partial{x}}-\dfrac{\partial{A_x}}{\partial{y}}\right)

たとえば,発散の例と同じベクトル場 A(x,y,z)=(xy,yz,zx)\boldsymbol{A}(x,y,z)=(xy,yz,zx) の回転は, rotA=(y(zx)z(yz),z(xy)x(zx),x(yz)y(xy))=(0y,0z,0x)=(y,z,x)\begin{aligned} \mathrm{rot}\,\boldsymbol{A} &=\left(\dfrac{\partial{}}{\partial{y}}(zx)-\dfrac{\partial{}}{\partial{z}}(yz),\dfrac{\partial{}}{\partial{z}}(xy)-\dfrac{\partial{}}{\partial{x}}(zx),\dfrac{\partial{}}{\partial{x}}(yz)-\dfrac{\partial{}}{\partial{y}}(xy)\right)\\ &=(0-y,0-z,0-x)\\ &=(-y,-z,-x) \end{aligned} となります。

ベクトル解析の回転(rot)の意味

回転の物理的な意味を考えてみましょう。 発散の場合と同様に,ベクトル場の例として流体の流れ V(x,y,z)\boldsymbol{V}(x,y,z) を考えます。

ある点 (x,y,z)(x,y,z) を中心として流体の回転の度合いを計算してみます。 下図のように,点 (x,y,z)(x,y,z) を中心とする xyxy 平面に平行で非常に小さな長方形を考えます(Δx,Δy1\varDelta x,\varDelta y\ll1)。

一辺の長さが dx, dy の長方形 x 方向に沿った流れ

まず xx 方向に沿った流れを計算します。 回転を考えたいので,y=y+Δyy=y+\varDelta y では xx 正の方向を負の回転,y=yΔyy=y-\varDelta y では xx 正の方向を正の回転と考えることにします。 点 (x,y±Δy/2,z)(x,y\pm\varDelta y/2,z) における xx 方向に沿った流れは Vx(x,y±Δy/2,z)V_x(x,y\pm\varDelta y/2,z) なので, y=y±Δyy=y\pm\varDelta y の線に沿った回転量を Vx(x,y±Δy/2,z)Δx\mp V_x(x,y\pm\varDelta y/2,z)\varDelta x に近似します。 xΔx/2x+Δx/2Vx(x,y±Δy2,z)dxVx(x,y±Δy2,z)Δx \int_{x-\varDelta x/2}^{x+\varDelta x/2} V_x\left(x,y\pm\dfrac{\varDelta y}{2},z\right)\,dx \simeq V_x\left(x,y\pm\dfrac{\varDelta y}{2},z\right)\varDelta x Δx\varDelta x は極めて小さいので,線上での VxV_x を一点の値として近似しています。

これを Taylor 展開して Δy\varDelta y の1次まで近似すると, Vx(x,y±Δy2,z)=Vx(x,y,z)±Vx(x,y,z)yΔy2+122Vx(x,y,z)y2(Δy2)2+Vx(x,y,z)±Vx(x,y,z)yΔy2\begin{aligned} &\quad V_x\left(x,y\pm\dfrac{\varDelta y}{2},z\right)\\ &=V_x(x,y,z) \pm\dfrac{\partial{V_x(x,y,z)}}{\partial{y}}\frac{\varDelta y}{2} +\dfrac{1}{2}\dfrac{\partial^2{V_x(x,y,z)}}{\partial{y}^2}\left(\frac{\varDelta y}{2}\right)^2 +\cdots\\ &\simeq V_x(x,y,z) \pm\dfrac{\partial{V_x(x,y,z)}}{\partial{y}}\frac{\varDelta y}{2} \end{aligned} となるので,正味の回転量は, [Vx(x,y+Δy2,z)+Vx(x,yΔy2,z)]Δx=[(Vx+VxyΔy2)(VxVxyΔy2)]Δx=VxyΔxΔy\begin{aligned} &\quad \left[ -V_x\left(x,y+\dfrac{\varDelta y}{2},z\right) +V_x\left(x,y-\dfrac{\varDelta y}{2},z\right)\right] \varDelta x \\ &=-\left[ \left(V_x+\frac{\partial V_x}{\partial y}\frac{\varDelta y}{2}\right) -\left(V_x-\frac{\partial V_x}{\partial y}\frac{\varDelta y}{2} \right) \right]\varDelta x\\ &=-\dfrac{\partial{V_x}}{\partial{y}}\varDelta x\varDelta y \end{aligned} となります。同様の計算を yy 方向に沿った2本の経路にも適応すると,正味の回転量は VyxΔxΔy\dfrac{\partial{V_y}}{\partial{x}}\varDelta x\varDelta y となります。yy に沿った方向では x=x+Δx/2x=x+\varDelta x/2yy 正の方向を正の回転,x=xΔx/2x=x-\varDelta x/2yy 正の方向を負の回転とするので,回転量の符号が変わることに気をつけてください。したがって4辺の正味の回転量は, (VyxVxy)ΔxΔy \left( \dfrac{\partial{V_y}}{\partial{x}}-\dfrac{\partial{V_x}}{\partial{y}} \right)\varDelta x\varDelta y となります。今度は V\boldsymbol{V} の回転の zz 成分 (×V)z(\nabla\times\boldsymbol{V})_z が出てきました。(×V)z(\nabla\times\boldsymbol{V})_z を用いて正味の回転量を書き換えると (×V)zΔxΔy(\nabla\times\boldsymbol{V})_z\varDelta x\varDelta y と表せます。 yzyz 平面,zxzx 平面上の回転量も同様に (×V)xΔyΔz,(×V)yΔzΔx(\nabla\times\boldsymbol{V})_x\varDelta y\varDelta z, (\nabla\times\boldsymbol{V})_y\varDelta z\varDelta x となります。

つまり,回転 V\nabla\cdot\boldsymbol{V} はその点まわりの(単位面積あたりの)物理的な回転を表していることがわかります。

回転 ×V\nabla\times\boldsymbol{V} はベクトルなので,スカラーである回転量と対応しないのではという疑問を抱いた方もいるかもしれません。 例えば,ベクトル n\boldsymbol{n} に垂直で面積が ΔS\varDelta S の平面を考えると,その平面上での回転量は (×V)ΔSn(\nabla\times\boldsymbol{V})\cdot \varDelta S\boldsymbol{n} と表されます。先の導出では n\boldsymbol{n}ez\boldsymbol{e}_z に取ることに対応します(実際に xyxy 平面を考えていました)。n\boldsymbol{n} を変化させたときに回転量 (×V)ΔSn(\nabla\times\boldsymbol{V})\cdot \varDelta S\boldsymbol{n} が最大になるのは n\boldsymbol{n}×V\nabla\times\boldsymbol{V} と平行なときです。

発散・回転の計算の具体例

例1:回転が0になる例

ベクトル場として V(x,y,z)=(x,y,z)\boldsymbol{V}(x,y,z)=(x,y,z) を考えます。 発散と回転を計算すると次の通りです。 V=xx+yy+zz=3×V=(yzzy,zxxz,xyyx)=(0,0,0)\begin{aligned} \nabla\cdot\boldsymbol{V} &=\dfrac{\partial{}}{\partial{x}}x+\dfrac{\partial{}}{\partial{y}}y+\dfrac{\partial{}}{\partial{z}}z =3\\ \nabla\times\boldsymbol{V} &=\left( \dfrac{\partial{}}{\partial{y}}z-\dfrac{\partial{}}{\partial{z}}y, \dfrac{\partial{}}{\partial{z}}x-\dfrac{\partial{}}{\partial{x}}z, \dfrac{\partial{}}{\partial{x}}y-\dfrac{\partial{}}{\partial{y}}x \right) =(0,0,0) \end{aligned}

物理的な意味を考えましょう。ベクトル場は下の図のようになっています。 原点から広がるようなベクトル場になっており,離れれば離れるほど V=VV=|\boldsymbol{V}| は大きくなります。 発散は V=3\nabla\cdot\boldsymbol{V}=3 なので,空間上すべての点で等しい湧き出しが生じています。 一方で ×V=0\nabla\times\boldsymbol{V}=\boldsymbol{0} なので,回転がどこにも存在していないということになります。 例1

例2:発散が0になる例

ベクトル場として V(x,y,z)=(y,z,x)V(x,y,z)=(-y,-z,-x) を考えます。 発散と回転は次のように計算できます。 V=x(y)+y(z)+z(x)=0×V=(y(x)z(z),z(y)x(x),x(z)y(y))=(1,1,1)\begin{aligned} \nabla\cdot\boldsymbol{V} &=\dfrac{\partial{}}{\partial{x}}(-y)+\dfrac{\partial{}}{\partial{y}}(-z)+\dfrac{\partial{}}{\partial{z}}(-x) =0\\ \nabla\times\boldsymbol{V} &=\left( \dfrac{\partial{}}{\partial{y}}(-x)-\dfrac{\partial{}}{\partial{z}}(-z), \dfrac{\partial{}}{\partial{z}}(-y)-\dfrac{\partial{}}{\partial{x}}(-x), \dfrac{\partial{}}{\partial{x}}(-z)-\dfrac{\partial{}}{\partial{y}}(-y) \right)\\ &=(1,1,1) \end{aligned} 例1とは反対に,発散が V=0\nabla\cdot\boldsymbol{V}=0 なので,どの点においても湧き出しが生じていないということになります。回転は ×V=(1,1,1)\nabla\times\boldsymbol{V}=(1,1,1) となり,例えば zz 平面に平行な面を考えると,(×V)ez=1(\nabla\times\boldsymbol{V})\cdot{\boldsymbol{e}}_z=1 なので正の方向(反時計回り)に回転が生じていることがわかります。 例2

2次元空間上の発散と回転

いままでは3次元の xyzxyz 空間を考えてきましたが,2次元の (x,y)(x,y) 平面でも発散と回転を定義することができます(3次元の場合で z=0,Az(x,y,0)=0z=0, A_z(x,y,0)=0 とした場合に対応します)。

2次元空間上の発散と回転

ベクトル場 A(x,y)=(Ax(x,y),Ay(x,y))\boldsymbol{A}(x,y)=(A_x(x,y),A_y(x,y)) に対して,2次元上の発散と回転は次のように定義される。 divA=Axx+AyyrotA=(AyxAxy)ez\begin{aligned} \mathrm{div}\,\boldsymbol{A} &=\dfrac{\partial{A_x}}{\partial{x}}+\dfrac{\partial{A_y}}{\partial{y}}\\ \mathrm{rot}\,\boldsymbol{A} &=\left(\dfrac{\partial{A_y}}{\partial{x}}-\dfrac{\partial{A_x}}{\partial{y}}\right){\boldsymbol{e}}_z \end{aligned}

具体例 V=(y,x)\boldsymbol{V}=(-y,x) を用いて発散と回転を計算してみましょう。 divV=x(y)+yx=0rotV=(xxy(y))ez=2ez\begin{aligned} \mathrm{div}\,\boldsymbol{V} &=\dfrac{\partial{}}{\partial{x}}(-y)+\dfrac{\partial{}}{\partial{y}}x =0\\ \mathrm{rot}\,\boldsymbol{V} &=\left( \dfrac{\partial{}}{\partial{x}}x-\dfrac{\partial{}}{\partial{y}}(-y)\right){\boldsymbol{e}}_z =2{\boldsymbol{e}}_z \end{aligned} 図を見ると回転になっていることが容易にわかると思います。

例3

流体を見つけたら落ち葉を浮かべましょう。