logo

磁場の発散と回転 📂電磁気学

磁場の発散と回転

定理

磁場発散回転は以下の通りだ。

B= 0×B= μ0J \begin{align*} \nabla \cdot \mathbf{B} =&\ 0 \\ \nabla \times \mathbf{B} =&\ \mu_{0} \mathbf{J} \end{align*}

説明

電場が常に0\mathbf{0}の回転を持つ特別なベクトル関数だったように、磁場もまたそうだ。体積電流に対するビオ・サバールの法則を使って発散と回転を計算しよう。

B(r)=μ04πJ(r)×2dτ \mathbf{B} (\mathbf{r}) = \dfrac{\mu_{0}}{4 \pi} \int \dfrac{\mathbf{J} (\mathbf{r}^{\prime}) \times \crH }{\cR^2} d \tau

ここでは、各関数がどの座標に依存しているかを明確に区別することが重要だ。

B\mathbf{B}観測点の座標(x,y,z)(x, y, z)に関する関数で、J\mathbf{J}源点の座標(x,y,z)(x^{\prime}, y^{\prime}, z^{\prime})に関する関数だ。分離ベクトルと体積要素はそれぞれ次の通り。

=(xx)x+(yy)y+(zz)z \bcR = (x-x^{\prime})\mathbf{x} + (y-y^{\prime})\mathbf{y} +(z-z^{\prime})\mathbf{z}

dτ=dxdydz d\tau^{\prime}=dx^{\prime}dy^{\prime}dz^{\prime}

積分は源点(プライム^{\prime}あり)に関する操作であり、発散と回転は観測点(プライム^{\prime}なし)に関する操作だ。

証明

磁場の発散

ビオ・サバールの法則により、磁場の発散は以下の通り。

B=μ04π(J×2)dτ \nabla \cdot \mathbf{B}=\dfrac{\mu_{0}}{4\pi} \int \nabla \cdot \left( \mathbf{J} \times \dfrac{\crH}{\cR^2}\right)d\tau^{\prime}

デル演算子が含まれる乗算規則

(A×B)=B(×A)A(×B) \nabla \cdot ( \mathbf{A} \times \mathbf{B} ) = \mathbf{B} \cdot (\nabla \times \mathbf{A}) - \mathbf{A} \cdot (\nabla \times \mathbf{B} )

上記の乗算規則を被積分関数に適用すると

(J×2)=2(×J)J(×2) \nabla \cdot \left( \mathbf{J} \times \dfrac{\crH}{\cR^2}\right) = \dfrac{\crH}{\cR^2}\cdot (\nabla \times \mathbf{J} ) -\mathbf{J} \cdot \left( \nabla \times \dfrac{\crH}{\cR^2} \right)

第一項で、\nablax,y,zx, y, zに対する微分操作であり、J\mathbf{J}(x,y,z)(x^{\prime}, y^{\prime}, z^{\prime})に関する関数だから×J=0\nabla \times \mathbf{J} = \mathbf{0}となる。さらに、分離ベクトルの回転が0\mathbf{0}だから、二番目の括弧の中も同様に0\mathbf{0}となる。したがって、磁場の発散は

B=μ04π0dτ=0 \nabla \cdot \mathbf{B} = \dfrac{\mu_{0}}{4\pi} \int 0 d\tau^{\prime} = 0

磁場の回転

ビオ・サバールの法則により、磁場の回転は以下の通り。

×B=μ04π×(J×2)dτ \nabla \times \mathbf{B} = \dfrac{\mu_{0}}{4\pi} \int \nabla \times \left( \mathbf{J} \times \dfrac{\crH}{\cR^2}\right)d\tau^{\prime}

デル演算子が含まれる乗算規則

×(A×B)=(B)A(A)B+A(B)B(A) \nabla \times (\mathbf{A} \times \mathbf{B}) = (\mathbf{B} \cdot \nabla)\mathbf{A} - (\mathbf{A} \cdot \nabla)\mathbf{B} + \mathbf{A} (\nabla \cdot \mathbf{B}) - \mathbf{B} (\nabla \cdot \mathbf{A})

上記の乗算規則を被積分関数に適用すると次のようになる。

×(J×2)=(2)J(J)2+J(2)2(J) \nabla \times \left( \mathbf{J} \times \dfrac{\crH}{\cR^2} \right)= \left( \dfrac{\crH}{\cR^2} \cdot \nabla \right) \mathbf{J} - \left( \mathbf{J} \cdot \nabla \right) \dfrac{\crH}{\cR^2} + \mathbf{J} \left( \nabla \cdot \dfrac{\crH}{\cR^2} \right) - \dfrac{\crH}{\cR^2} \left( \nabla \cdot \mathbf{J} \right)

ここで、第一項と第四項は\nablaJ\mathbf{J}に適用されるから、先の発散の計算で述べた理由と同様、00となる。第二項は直接00ではないが、積分を行うと00となる。この過程は下で別に説明する。最後に、第三項は分離ベクトルの発散なのでディラックのデルタ関数と同じだ。

(^2)=4πδ2() \nabla \cdot \left( \dfrac{ \hat { \bcR} }{\cR^2} \right) = 4\pi \delta^2(\bcR )

したがって、磁場の回転は以下の通り。

×B=μ04πJ(r)4πδ3(rr)dτ \nabla \times \mathbf{B} = \dfrac{\mu_{0}}{4\pi} \int \mathbf{J} (\mathbf{r}^{\prime} ) 4\pi \delta^3(\mathbf{r} - \mathbf{r}^{\prime} ) d\tau^{\prime}

ディラックのデルタ関数の定義によれば

×B=μ0J(r) \nabla \times \mathbf{B} = \mu_{0} \mathbf{J}(\mathbf{r})

今、第二項の積分がなぜ00となるかを確認することが残っている。この部分が難しい場合は、数学的にではなく、物理的な意味によって00と理解しても良い。分離ベクトルが=(xx)x^+(yy)y^+(zz)z^\bcR = (x-x^{\prime})\hat{\mathbf{x}} + (y-y^{\prime})\hat{\mathbf{y}} + (z-z^{\prime})\hat{\mathbf{z}}となるので、以下の式が成り立つ。

x=x^=xandy=y^=yandz=z^=z -\dfrac{\partial \bcR}{\partial x} = - \hat{\mathbf{x}} = \dfrac{\partial \bcR}{\partial x^{\prime}} \quad \text{and} \quad -\dfrac{\partial \bcR}{\partial y} = - \hat{\mathbf{y}} = \dfrac{\partial \bcR}{\partial y^{\prime}} \quad \text{and} \quad -\dfrac{\partial \bcR}{\partial z} = - \hat{\mathbf{z}} = \dfrac{\partial \bcR}{\partial z^{\prime}}

したがって、以下の等式が成り立つ。

(J)2= (Jxx+Jyy+Jzz)2= (Jxx+Jyy+Jzz)2=(J)2 \begin{align*} -( \mathbf{J} \cdot \nabla) \dfrac{\crH}{\cR^2} =&\ -\left( J_{x^{\prime}} \dfrac{\partial }{\partial x} + J_{y^{\prime}} \dfrac{\partial }{\partial y} + J_{z^{\prime}} \dfrac{\partial }{\partial z}\right) \dfrac{\crH}{\cR^2} \\ =&\ \left( J_{x^{\prime}} \dfrac{\partial }{\partial x^{\prime}} + J_{y^{\prime}} \dfrac{\partial }{\partial y^{\prime}} + J_{z^{\prime}} \dfrac{\partial }{\partial z^{\prime}}\right) \dfrac{\crH}{\cR^2} \\ =&( \mathbf{J} \cdot \nabla^{\prime} ) \dfrac{\crH}{\cR^2} \end{align*}

このベクトルのx\mathbf{x}成分は以下の通り。

[(J)2]x=(J)(xx3) \left[ ( \mathbf{J} \cdot \nabla^{\prime} ) \dfrac{\crH}{\cR^2} \right]_{x}=( \mathbf{J} \cdot \nabla^{\prime} ) \left( \dfrac{ x-x^{\prime}}{\cR^3} \right)

これを乗算規則 (fA)=f(A)+A(f)\nabla \cdot (f\mathbf{A}) = f(\nabla \cdot \mathbf{A}) + \mathbf{A} \cdot (\nabla f)で解くと以下の通り。

(J)(xx3)=[xx3J](xx3)(J) ( \mathbf{J} \cdot \nabla^{\prime} ) \left( \dfrac{ x-x^{\prime}}{\cR^3} \right) = \nabla^{\prime} \cdot \left[ \dfrac{x-x^{\prime}}{\cR^3} \mathbf{J} \right] - \left( \dfrac{x-x^{\prime}} {\cR^3} \right) (\nabla^{\prime} \cdot \mathbf{J} )

さらに、定常電流ならばJ\mathbf{J}の発散が00であるから、第二項は00となる。残りの項を積分すると、発散定理により

V[xx3J]dτ=Sxx3Jda \int_\mathcal{V} \nabla^{\prime} \cdot \left[ \dfrac{x-x^{\prime}}{\cR^3}\mathbf{J} \right]d\tau=\oint_\mathcal{S} \dfrac{x-x^{\prime}}{\cR^3}\mathbf{J} \cdot d\mathbf{a}^{\prime}

上で、\nabla\nabla^{\prime}に変えた理由はこの積分のためだ。左辺の積分範囲は全ての電流を含めるために十分に大きく設定しなければならず、この場合、領域の境界面でJ=0\mathbf{J}=0である。だから、右辺の面積積分で領域全体でJ=0\mathbf{J}=0であるため、積分結果は00となる。同じ論理で、y\mathbf{y}成分とz\mathbf{z}成分も00であるので、第二項の積分は00となる。