静磁場




静磁場（せいじば、英語: static magnetic field）とは、時間的に変動しない磁場のことである。
[1][2][3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28][29]
本記事では、静磁気学 (Magnetostatics) の視点から、静磁場について述べる。




目次






  • 1 定常な電流が作り出す静磁場の一般論


  • 2 時間的に定常な磁化が作り出す静磁場の一般論


    • 2.1 時間的に定常な磁化が作り出す磁気ベクトルポテンシャル


    • 2.2 時間的に定常な磁化が作り出す磁束密度


    • 2.3 磁荷密度の導入


    • 2.4 磁化電流の導入




  • 3 強制電流と磁化の両方が既知のとき(一般論)


    • 3.1 全系の磁束密度


    • 3.2 B=μ0(H+M)について


    • 3.3 透磁率の導入


    • 3.4 静磁場の境界条件


    • 3.5 静磁場のエネルギーと停留値関数




  • 4 ベクトル解析に関する補足


    • 4.1 ベクトル場の代数演算と微分作用素


      • 4.1.1 (1) <F|∇>について


      • 4.1.2 (2)rot とスカラー倍、ベクトル積


      • 4.1.3 (3)div とスカラー倍、ベクトル積




    • 4.2 ディラックのデルタと、重積分


    • 4.3 発散微分とディラックのデルタ


    • 4.4 分数関数の偏微分


    • 4.5 面積分に関する補足


      • 4.5.1 (1)普通の面積分


      • 4.5.2 (2)絶対値による面積分


      • 4.5.3 (3)外積面積分






  • 5 脚注


    • 5.1 注釈


    • 5.2 参考文献




  • 6 関連項目


  • 7 外部リンク





定常な電流が作り出す静磁場の一般論


本節では、真空中に定常な（つまり時刻tに依存しない）電流密度が作り出す磁束密度について、一般に成り立つ事柄について述べる。ただし、時間的な変動の影響はもちろんのこと、これ以外にも、電場や強制電荷、分極電荷の影響は排除されているものとする。
本記事では、専ら体積電流密度を中心に扱い、線電流近似については、例えば[30][31][32]等に委ねることとする。




真空中に定常な（つまり時刻tに依存しない）電流密度
i(r){displaystyle {boldsymbol {i}}(mathbf {r} )}{displaystyle {boldsymbol {i}}(mathbf {r} )}が与えられた
とする。このとき、i(r){displaystyle {boldsymbol {i}}(mathbf {r} )}{displaystyle {boldsymbol {i}}(mathbf {r} )}
は、以下の磁気ベクトルポテンシャル
Ai(r){displaystyle mathbf {A} _{boldsymbol {i}}(mathbf {r} )}{displaystyle mathbf {A} _{boldsymbol {i}}(mathbf {r} )}を空間内に作り出す。



Ai(r)=μ04π∫s∈R3(i(s)|r−s|) d3s{displaystyle mathbf {A} _{boldsymbol {i}}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in mathbb {R} ^{3}}left({frac {{boldsymbol {i}}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right) d^{3}mathbf {s} }{displaystyle mathbf {A} _{boldsymbol {i}}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in mathbb {R} ^{3}}left({frac {{boldsymbol {i}}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right) d^{3}mathbf {s} } （1-1）

となる。


Bi=rot⁡[Ai]{displaystyle mathbf {B} _{boldsymbol {i}}=operatorname {rot} [mathbf {A} _{boldsymbol {i}}]}{displaystyle mathbf {B} _{boldsymbol {i}}=operatorname {rot} [mathbf {A} _{boldsymbol {i}}]}を考え併せると、
i(r){displaystyle {boldsymbol {i}}(mathbf {r} )}{displaystyle {boldsymbol {i}}(mathbf {r} )}が直接的に作り出す磁束密度Bi{displaystyle mathbf {B} _{boldsymbol {i}}}{displaystyle mathbf {B} _{boldsymbol {i}}}は、



Bi(r)=μ04π∫s∈R3(i(s)×(r−s)|r−s|3)d3s{displaystyle mathbf {B} _{boldsymbol {i}}(mathbf {r} )={frac {mu _{0}}{4pi }}int _{mathbf {s} in mathbb {R} ^{3}}left({frac {{boldsymbol {i}}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right)d^{3}mathbf {s} }{displaystyle mathbf {B} _{boldsymbol {i}}(mathbf {r} )={frac {mu _{0}}{4pi }}int _{mathbf {s} in mathbb {R} ^{3}}left({frac {{boldsymbol {i}}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right)d^{3}mathbf {s} } （1-2）

となる。これは、すなわち、ビオ・サバールの法則である。


上記のBi(r){displaystyle mathbf {B} _{boldsymbol {i}}(mathbf {r} )}{displaystyle mathbf {B} _{boldsymbol {i}}(mathbf {r} )}に対し、
新たな場Hi{displaystyle {mathbf {H} }_{boldsymbol {i}}}{displaystyle {mathbf {H} }_{boldsymbol {i}}}を、



Hi(r):=1μ0Bi(r){displaystyle mathbf {H} _{boldsymbol {i}}(mathbf {r} ):={frac {1}{{mu }_{0}}}mathbf {B} _{boldsymbol {i}}(mathbf {r} )}{displaystyle mathbf {H} _{boldsymbol {i}}(mathbf {r} ):={frac {1}{{mu }_{0}}}mathbf {B} _{boldsymbol {i}}(mathbf {r} )} (1-3)

と定義する。この場のことを、「電流密度iが作り出す磁場」と呼ぶ。ここでμ0は、真空の透磁率である。
尚、定義の上では、「電流密度iが作り出す磁場」Hi{displaystyle {mathbf {H} }_{boldsymbol {i}}}{displaystyle {mathbf {H} }_{boldsymbol {i}}}は、
透磁率がμの場所でも、:Hi(r):=1μ0Bi(r){displaystyle mathbf {H} _{boldsymbol {i}}(mathbf {r} ):={frac {1}{{mu }_{0}}}mathbf {B} _{boldsymbol {i}}(mathbf {r} )}{displaystyle mathbf {H} _{boldsymbol {i}}(mathbf {r} ):={frac {1}{{mu }_{0}}}mathbf {B} _{boldsymbol {i}}(mathbf {r} )}であることに特に注意されたい。


式(1-2),式(1-3)より、



Hi(r)=14π∫s∈R3(i(s)×(r−s)|r−s|3)d3s{displaystyle mathbf {H} _{boldsymbol {i}}(mathbf {r} )={frac {1}{4pi }}int _{mathbf {s} in mathbb {R} ^{3}}left({frac {{boldsymbol {i}}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right)d^{3}mathbf {s} }{displaystyle mathbf {H} _{boldsymbol {i}}(mathbf {r} )={frac {1}{4pi }}int _{mathbf {s} in mathbb {R} ^{3}}left({frac {{boldsymbol {i}}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right)d^{3}mathbf {s} } (1-4)

である。これに、回転微分を作用させると、



rotr⁡[Hi]=i{displaystyle operatorname {rot} _{mathbf {r} }[mathbf {H} _{boldsymbol {i}}]={boldsymbol {i}}}{displaystyle operatorname {rot} _{mathbf {r} }[mathbf {H} _{boldsymbol {i}}]={boldsymbol {i}}} (1-5)

が得られる。実際、ベクトル解析の公式より、



rotr⁡[(i(s)×(r−s)|r−s|3)]=4π(divr⁡[(r−s)|r−s|3])i(s)=4πδ3(r−s)i(s)(1-6){displaystyle {begin{aligned}operatorname {rot} _{mathbf {r} }left[left({frac {{boldsymbol {i}}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right)right]&=4pi left(operatorname {div} _{mathbf {r} }left[{frac {(mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right]right){boldsymbol {i}}(mathbf {s} )\&=4pi {delta }^{3}(mathbf {r} -mathbf {s} ){boldsymbol {i}}(mathbf {s} )qquad {text{(1-6)}}end{aligned}}}{displaystyle {begin{aligned}operatorname {rot} _{mathbf {r} }left[left({frac {{boldsymbol {i}}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right)right]&=4pi left(operatorname {div} _{mathbf {r} }left[{frac {(mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right]right){boldsymbol {i}}(mathbf {s} )\&=4pi {delta }^{3}(mathbf {r} -mathbf {s} ){boldsymbol {i}}(mathbf {s} )qquad {text{(1-6)}}end{aligned}}}

従って、


rotr⁡[Hi(r)]=14π∫s∈R3rotr⁡[(i(s)×(r−s)|r−s|3)]d3s=∫s∈R3δ3(r−s)i(s)d3s=i(r)(1-7){displaystyle {begin{aligned}operatorname {rot} _{mathbf {r} }[mathbf {H} _{boldsymbol {i}}(mathbf {r} )]&={frac {1}{4pi }}int _{mathbf {s} in mathbb {R} ^{3}}operatorname {rot} _{mathbf {r} }left[left({frac {{boldsymbol {i}}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right)right]d^{3}mathbf {s} \&={int }_{mathbf {s} in mathbb {R} ^{3}}{delta }^{3}(mathbf {r} -mathbf {s} ){boldsymbol {i}}(mathbf {s} )d^{3}mathbf {s} \&={boldsymbol {i}}(mathbf {r} )qquad {text{(1-7)}}end{aligned}}}{displaystyle {begin{aligned}operatorname {rot} _{mathbf {r} }[mathbf {H} _{boldsymbol {i}}(mathbf {r} )]&={frac {1}{4pi }}int _{mathbf {s} in mathbb {R} ^{3}}operatorname {rot} _{mathbf {r} }left[left({frac {{boldsymbol {i}}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right)right]d^{3}mathbf {s} \&={int }_{mathbf {s} in mathbb {R} ^{3}}{delta }^{3}(mathbf {r} -mathbf {s} ){boldsymbol {i}}(mathbf {s} )d^{3}mathbf {s} \&={boldsymbol {i}}(mathbf {r} )qquad {text{(1-7)}}end{aligned}}}

が得られる。ここでδ3は、3変数のδ関数（ディラックのデルタ）を意味する。




時間的に定常な磁化が作り出す静磁場の一般論


本節では、定常な（つまり時刻tに依存しない）磁化が作り出す磁束密度について、一般に成り立つ事柄について述べる。ただし、時間的な変動の影響はもちろんのこと、これ以外にも、電場や強制電荷、分極電荷の影響は排除されているものとする。





時間的に定常な磁化が作り出す磁気ベクトルポテンシャル


空間内の領域Ω{displaystyle Omega }Omega に物質が置かれ、前記物質が、定常な（つまり時刻tに依存しない）磁化M(r){displaystyle mathbf {M} (mathbf {r} )}{displaystyle mathbf {M} (mathbf {r} )}を帯びているとする。このとき、磁化ベクトルは、「単位体積当たりの磁気モーメントの密度」を表すものであるため、磁化の定義より、



物質内の各点s∈Ω{displaystyle mathbf {s} in Omega }{displaystyle mathbf {s} in Omega }それぞれに、それぞれ

M(s)d3s{displaystyle mathbf {M} (mathbf {s} ){d}^{3}mathbf {s} }{displaystyle mathbf {M} (mathbf {s} ){d}^{3}mathbf {s} } （2-1-1）


で与えられる磁気モーメントが配置されている


と考えることができる[2]。


まず、時間的に定常な磁化が作り出す磁気ベクトルポテンシャルについて考えよう。



(1)原点に置かれた磁気モーメントmは、空間上の位置rに作り出す磁気ベクトルポテンシャルは、



a0(r)=μ04π(m×r|r|){displaystyle mathbf {a} _{0}(mathbf {r} )={frac {mu _{0}}{4pi }}left({frac {mathbf {m} times mathbf {r} }{|mathbf {r} |}}right)}{displaystyle mathbf {a} _{0}(mathbf {r} )={frac {mu _{0}}{4pi }}left({frac {mathbf {m} times mathbf {r} }{|mathbf {r} |}}right)} （2-1-2）
である。



(2)従って、"(1)"を平行移動すれば、位置s∈R3{displaystyle mathbf {s} in mathbb {R} ^{3}}{displaystyle mathbf {s} in mathbb {R} ^{3}}に置かれた磁気モーメントmは、空間上の位置rに作り出す磁気ベクトルポテンシャルは、



as(r)=μ04π(m×(r−s)|r−s|){displaystyle mathbf {a} _{s}(mathbf {r} )={frac {mu _{0}}{4pi }}left({frac {mathbf {m} times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right)}{displaystyle mathbf {a} _{s}(mathbf {r} )={frac {mu _{0}}{4pi }}left({frac {mathbf {m} times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right)} （2-1-3）
である。



従って、上記の磁気モーメントM(s)ds{displaystyle mathbf {M} (mathbf {s} )dmathbf {s} }{displaystyle mathbf {M} (mathbf {s} )dmathbf {s} }それぞれは、磁気ベクトルポテンシャル



μ04π((M(s)ds)×(r−s)|r−s|){displaystyle {frac {mu _{0}}{4pi }}left({frac {(mathbf {M} (mathbf {s} )dmathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right)}{displaystyle {frac {mu _{0}}{4pi }}left({frac {(mathbf {M} (mathbf {s} )dmathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right)} （2-1-4）

を作り出す。


上記の磁気ベクトルポテンシャルそれぞれを、全てのs∈Ω{displaystyle mathbf {s} in Omega }{displaystyle mathbf {s} in Omega }に渡って足し合わせると、



AM(r)=μ04π∫s∈Ω((M(s))×(r−s)|r−s|) d3s{displaystyle mathbf {A} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }left({frac {(mathbf {M} (mathbf {s} ))times (mathbf {r} -{mathbf {s} })}{|mathbf {r} -mathbf {s} |}}right) {d}^{3}mathbf {s} }{displaystyle mathbf {A} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }left({frac {(mathbf {M} (mathbf {s} ))times (mathbf {r} -{mathbf {s} })}{|mathbf {r} -mathbf {s} |}}right) {d}^{3}mathbf {s} } （2-1-5）

を得る。即ち、物質の磁化M(r){displaystyle mathbf {M} (mathbf {r} )}{displaystyle mathbf {M} (mathbf {r} )}
は、上記の磁気ベクトルポテンシャル
AM{displaystyle mathbf {A} _{M}}{displaystyle mathbf {A} _{M}}を空間内に作り出す[2][注釈 1][注釈 2]



時間的に定常な磁化が作り出す磁束密度


次に、時間的に定常な磁化が作り出す磁束密度について考えよう。磁気ベクトルポテンシャルの回転微分をとれば、磁束密度が得られる。



(1)原点に置かれた磁気モーメントmは、空間上の位置rに作り出す磁気ベクトルポテンシャルは、前述の通り、




a0(r)=μ04π(m×r|r|){displaystyle mathbf {a} _{0}(mathbf {r} )={frac {mu _{0}}{4pi }}left({frac {mathbf {m} times mathbf {r} }{|mathbf {r} |}}right)}{displaystyle mathbf {a} _{0}(mathbf {r} )={frac {mu _{0}}{4pi }}left({frac {mathbf {m} times mathbf {r} }{|mathbf {r} |}}right)} （2-2-1）

である。これの回転微分をとることで、原点に置かれた磁気モーメントmが、空間上の位置rに作り出す磁束密度は、

b0(r)=rotr⁡[μ04π(m×r|r|)]=μ04πgradr⁡[<m|r>|r|3]+μ04πdivr⁡[r|r|]m{displaystyle mathbf {b} _{0}(mathbf {r} )=operatorname {rot} _{mathbf {r} }left[{frac {mu _{0}}{4pi }}left({frac {mathbf {m} times mathbf {r} }{|mathbf {r} |}}right)right]={frac {{mu }_{0}}{4pi }}operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {m} |mathbf {r} >}{{|mathbf {r} |}^{3}}}right]+{frac {{mu }_{0}}{4pi }}operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} }{|mathbf {r} |}}right]mathbf {m} }{displaystyle mathbf {b} _{0}(mathbf {r} )=operatorname {rot} _{mathbf {r} }left[{frac {mu _{0}}{4pi }}left({frac {mathbf {m} times mathbf {r} }{|mathbf {r} |}}right)right]={frac {{mu }_{0}}{4pi }}operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {m} |mathbf {r} >}{{|mathbf {r} |}^{3}}}right]+{frac {{mu }_{0}}{4pi }}operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} }{|mathbf {r} |}}right]mathbf {m} } （2-2-2）


であることが判る。ここで、< | >{displaystyle < | >}{displaystyle < | >}は、内積を表す。



(2)従って、"(1)"を平行移動すれば、位置s∈R3{displaystyle mathbf {s} in mathbb {R} ^{3}}{displaystyle mathbf {s} in mathbb {R} ^{3}}に置かれた磁気モーメントmが、空間上の位置r
に作り出す磁束密度が、




bs(r)=rotr⁡[μ04π(m×(r−s)|r−s|)]=μ04π gradr⁡[<m|(r−s)>|r−s|3]+μ04π divr⁡[r−s|r−s|]m{displaystyle mathbf {b} _{s}(mathbf {r} )=operatorname {rot} _{mathbf {r} }left[{frac {mu _{0}}{4pi }}left({frac {mathbf {m} times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right)right]={frac {{mu }_{0}}{4pi }} operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {m} |(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]+{frac {{mu }_{0}}{4pi }} operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {m} }{displaystyle mathbf {b} _{s}(mathbf {r} )=operatorname {rot} _{mathbf {r} }left[{frac {mu _{0}}{4pi }}left({frac {mathbf {m} times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right)right]={frac {{mu }_{0}}{4pi }} operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {m} |(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]+{frac {{mu }_{0}}{4pi }} operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {m} } （2-2-3）

であることが判る。




従って、上記の磁気モーメントM(s)ds{displaystyle mathbf {M} (mathbf {s} )dmathbf {s} }{displaystyle mathbf {M} (mathbf {s} )dmathbf {s} }それぞれは、磁束密度



μ04πgradr⁡[<M(s)d3s|(r−s)>|r−s|3]+μ04πdivr⁡[r−s|r−s|]M(s)d3s{displaystyle {frac {{mu }_{0}}{4pi }}operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} ){d}^{3}mathbf {s} |(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]+{frac {{mu }_{0}}{4pi }}operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {M} (mathbf {s} ){d}^{3}mathbf {s} }{displaystyle {frac {{mu }_{0}}{4pi }}operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} ){d}^{3}mathbf {s} |(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]+{frac {{mu }_{0}}{4pi }}operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {M} (mathbf {s} ){d}^{3}mathbf {s} } （2-2-4）

を作り出す。


上記の磁気ベクトルポテンシャルそれぞれを、全てのs∈Ω{displaystyle mathbf {s} in Omega }{displaystyle mathbf {s} in Omega }に渡って足し合わせると、



BM(r)=μ04π∫s∈Ωgradr⁡[<M(s)|(r−s)>|r−s|3]+divr⁡[r−s|r−s|]M(s) d3s{displaystyle mathbf {B} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]+operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {M} (mathbf {s} ) {d}^{3}mathbf {s} }{displaystyle mathbf {B} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]+operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {M} (mathbf {s} ) {d}^{3}mathbf {s} } （2-2-5）

を得る。即ち、物質の磁化M(r){displaystyle mathbf {M} (mathbf {r} )}{displaystyle mathbf {M} (mathbf {r} )}
は、上記の磁気ベクトルポテンシャル
BM{displaystyle mathbf {B} _{M}}{displaystyle mathbf {B} _{M}}を空間内に作り出す
[2][注釈 3][注釈 4]


さて、ベクトル解析の公式から、



M(s) = ∫s∈Ω divr⁡[r−s|r−s|]M(s) d3s{displaystyle mathbf {M} (mathbf {s} ) = {int }_{mathbf {s} in Omega } operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {M} (mathbf {s} ) {d}^{3}s}{displaystyle mathbf {M} (mathbf {s} ) = {int }_{mathbf {s} in Omega } operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {M} (mathbf {s} ) {d}^{3}s} （2-2-6）

が判る。[注釈 5]


従って、上記の磁気モーメントM(s)ds{displaystyle mathbf {M} (mathbf {s} )dmathbf {s} }{displaystyle mathbf {M} (mathbf {s} )dmathbf {s} }それぞれが作り出す磁束密度は、



μ04πgradr⁡[<M(s)d3s|(r−s)>|r−s|3]+μ0δ(r−s)M(s)d3s{displaystyle {frac {{mu }_{0}}{4pi }}operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} ){d}^{3}mathbf {s} |(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]+{mu }_{0}delta (mathbf {r} -mathbf {s} )mathbf {M} (mathbf {s} ){d}^{3}mathbf {s} }{displaystyle {frac {{mu }_{0}}{4pi }}operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} ){d}^{3}mathbf {s} |(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]+{mu }_{0}delta (mathbf {r} -mathbf {s} )mathbf {M} (mathbf {s} ){d}^{3}mathbf {s} } （2-2-7）

と書けることが判る。従って、



BM(r)=μ04π∫s∈Ωgradr⁡[<M(s)|(r−s)>|r−s|3] d3s+μ0M(r){displaystyle mathbf {B} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right] {d}^{3}{s}+mu _{0}mathbf {M} (mathbf {r} )}{displaystyle mathbf {B} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right] {d}^{3}{s}+mu _{0}mathbf {M} (mathbf {r} )} （2-2-8）

を得る。[注釈 6]


今、新たな場HM{displaystyle mathbf {H} _{M}}{displaystyle mathbf {H} _{M}}を、



HM(r):=14π∫s∈Ωgradr⁡[<M(s)|(r−s)>|r−s|3] d3s{displaystyle mathbf {H} _{M}(mathbf {r} ):={frac {1}{4pi }}{int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right] {d}^{3}{s}}{displaystyle mathbf {H} _{M}(mathbf {r} ):={frac {1}{4pi }}{int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right] {d}^{3}{s}} （2-2-9）

と定義すると
[注釈 7]、



BM(r)=μ0(HM(r)+M(r)){displaystyle mathbf {B} _{M}(mathbf {r} )={mu }_{0}(mathbf {H} _{M}(mathbf {r} )+mathbf {M} (mathbf {r} ))}{displaystyle mathbf {B} _{M}(mathbf {r} )={mu }_{0}(mathbf {H} _{M}(mathbf {r} )+mathbf {M} (mathbf {r} ))} （2-2-10）

が得られる。


さて、以上の議論から「物質の磁化Mが既知である場合に限れば、その磁化Mが作り出す磁束密度BMを計算する術が得られた」
ことになる。然しながら、物質の磁化Mが既知でない場合には、上述の関係式のみからは、BMもMも、判らない。上述の関係式は、
一つの拘束条件を与えているに過ぎないのである。（だが、極めて重要な関係式ではある）。


例えば、磁化Miniを帯びた鉄心が空間におかれていたとき、外部からの磁束密度Bext(r){displaystyle mathbf {B} _{text{ext}}(mathbf {r} )}{displaystyle mathbf {B} _{text{ext}}(mathbf {r} )}
が与えられたとき、鉄心の磁化は、元々の磁化Miniと、外部からの磁束密度Bext(r){displaystyle mathbf {B} _{text{ext}}(mathbf {r} )}{displaystyle mathbf {B} _{text{ext}}(mathbf {r} )}
の影響で、元々の磁化Miniとは異なる新たな磁化Mconを得ることになる。仮に、このMcon
が判れば、全系の磁束密度が計算できるのだが、


（難所）元々磁化Miniを帯びている物質に、外部から磁束密度Bext(r){displaystyle mathbf {B} _{text{ext}}(mathbf {r} )}{displaystyle mathbf {B} _{text{ext}}(mathbf {r} )}を印加したとき、最終的に、物質がどのようなMconを得るか？

が、実のところは難しい。そこで、一般には。B-H曲線等の実測結果と、上記の拘束条件を考え合わせ、数値計算によって磁束密度や磁化が計算されるのである。但し、線形物質に関して言えば、上の難所は比較的簡単である。このような特殊な物質に関する問題については、次章で述べることにする。



磁荷密度の導入


空間内の領域Ω{displaystyle Omega }Omega に物質が置かれ、前記物質が、定常な（つまり時刻tに依存しない）磁化M(r){displaystyle mathbf {M} (mathbf {r} )}{displaystyle mathbf {M} (mathbf {r} )}を帯びているとする。


先に定義した、（磁化による）磁場



HM(r)=14π∫s∈Ωgradr⁡[<M(s)|(r−s)>|r−s|3] d3s{displaystyle mathbf {H} _{M}(mathbf {r} )={frac {1}{4pi }}{int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right] {d}^{3}{s}}{displaystyle mathbf {H} _{M}(mathbf {r} )={frac {1}{4pi }}{int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right] {d}^{3}{s}} （2-3-1）

の原因が、磁荷であると考えた場合に、それは、どのようなものであるのかを検討しよう。


ベクトル解析の公式より、


divs⁡[M(s)|r−s|]=⟨grads⁡[1|r−s|]|M(s)⟩+divs⁡[M]|r−s|=⟨r−s|r−s|3|M(s)⟩+divs⁡[M]|r−s|=<M(s)|(r−s)>|r−s|3+divs⁡[M]|r−s|(2-3-2){displaystyle {begin{aligned}operatorname {div} _{s}left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right]&=leftlangle operatorname {grad} _{mathbf {s} }left[{frac {1}{|mathbf {r} -mathbf {s} |}}right]left|right.mathbf {M} (mathbf {s} )rightrangle +{frac {operatorname {div} _{mathbf {s} }[mathbf {M} ]}{|mathbf {r} -mathbf {s} |}}\&=leftlangle {frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |^{3}}}left|right.mathbf {M} (mathbf {s} )rightrangle +{frac {operatorname {div} _{mathbf {s} }[mathbf {M} ]}{|mathbf {r} -mathbf {s} |}}\&={frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}+{frac {operatorname {div} _{mathbf {s} }[mathbf {M} ]}{|mathbf {r} -mathbf {s} |}}qquad {text{(2-3-2)}}end{aligned}}}{displaystyle {begin{aligned}operatorname {div} _{s}left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right]&=leftlangle operatorname {grad} _{mathbf {s} }left[{frac {1}{|mathbf {r} -mathbf {s} |}}right]left|right.mathbf {M} (mathbf {s} )rightrangle +{frac {operatorname {div} _{mathbf {s} }[mathbf {M} ]}{|mathbf {r} -mathbf {s} |}}\&=leftlangle {frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |^{3}}}left|right.mathbf {M} (mathbf {s} )rightrangle +{frac {operatorname {div} _{mathbf {s} }[mathbf {M} ]}{|mathbf {r} -mathbf {s} |}}\&={frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}+{frac {operatorname {div} _{mathbf {s} }[mathbf {M} ]}{|mathbf {r} -mathbf {s} |}}qquad {text{(2-3-2)}}end{aligned}}}

従って、



<M(s)|(r−s)>|r−s|3=divs⁡[M(s)|r−s|]−divs⁡[M]|r−s|{displaystyle {frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}=operatorname {div} _{s}left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right]-{frac {operatorname {div} _{mathbf {s} }[mathbf {M} ]}{|mathbf {r} -mathbf {s} |}}}{displaystyle {frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}=operatorname {div} _{s}left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right]-{frac {operatorname {div} _{mathbf {s} }[mathbf {M} ]}{|mathbf {r} -mathbf {s} |}}} （2-3-3）

今、スカラー値関数ϕM{displaystyle {phi }_{M}}{displaystyle {phi }_{M}}を、



ϕM(r):=−∫s∈Ω<M(s)|(r−s)>|r−s|3 d3Ω{displaystyle {phi }_{M}(mathbf {r} ):=-{int }_{mathbf {s} in Omega }{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}} {d}^{3}Omega }{displaystyle {phi }_{M}(mathbf {r} ):=-{int }_{mathbf {s} in Omega }{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}} {d}^{3}Omega } （2-3-4）

と定める（磁化が定める磁位、あるいは、磁化が定めるスカラーポテンシャルと言う）と、



ϕM(r)=−∫s∈Ωdivs⁡[M(s)|r−s|] d3Ω+∫s∈Ωdivs⁡[M]|r−s| d3Ω{displaystyle {phi }_{M}(mathbf {r} )=-{int }_{mathbf {s} in Omega }operatorname {div} _{s}left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}Omega +{int }_{mathbf {s} in Omega }{frac {operatorname {div} _{mathbf {s} }[mathbf {M} ]}{|mathbf {r} -mathbf {s} |}} {d}^{3}Omega }{displaystyle {phi }_{M}(mathbf {r} )=-{int }_{mathbf {s} in Omega }operatorname {div} _{s}left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}Omega +{int }_{mathbf {s} in Omega }{frac {operatorname {div} _{mathbf {s} }[mathbf {M} ]}{|mathbf {r} -mathbf {s} |}} {d}^{3}Omega } （2-3-5）

右辺第一項に、ガウスの発散定理を適用すると、



∫s∈Ωdivs⁡[M(s)|r−s|] d3Ω=∫s∈∂ΩM(s)|r−s| d2(∂Ω){displaystyle {int }_{mathbf {s} in Omega }operatorname {div} _{s}left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}Omega ={int }_{mathbf {s} in partial Omega }{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}} {d}^{2}(partial Omega )}{displaystyle {int }_{mathbf {s} in Omega }operatorname {div} _{s}left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}Omega ={int }_{mathbf {s} in partial Omega }{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}} {d}^{2}(partial Omega )} （2-3-6）

さらに、ベクトル解析の公式を適用すると、



∫s∈∂ΩM(s)|r−s| d2(∂Ω)=∫s∈∂Ω<M(s)|n∂Ω>|r−s| |d2(∂Ω)|{displaystyle {int }_{mathbf {s} in partial Omega }{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}} {d}^{2}(partial Omega )={int }_{mathbf {s} in partial Omega }{frac {<mathbf {M} (mathbf {s} )|mathbf {n} _{partial Omega }>}{|mathbf {r} -mathbf {s} |}} |{d}^{2}(partial Omega )|}{displaystyle {int }_{mathbf {s} in partial Omega }{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}} {d}^{2}(partial Omega )={int }_{mathbf {s} in partial Omega }{frac {<mathbf {M} (mathbf {s} )|mathbf {n} _{partial Omega }>}{|mathbf {r} -mathbf {s} |}} |{d}^{2}(partial Omega )|} （2-3-7）

従って、



ϕM(r)=−∫s∈Ωdivs⁡[M(s)|r−s|] d3Ω +∫s∈∂Ω<M(s)|n∂Ω>|r−s| |d2(∂Ω)|{displaystyle {phi }_{M}(mathbf {r} )=-{int }_{mathbf {s} in Omega }operatorname {div} _{s}left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}Omega +{int }_{mathbf {s} in partial Omega }{frac {<mathbf {M} (mathbf {s} )|mathbf {n} _{partial Omega }>}{|mathbf {r} -mathbf {s} |}} |{d}^{2}(partial Omega )|}{displaystyle {phi }_{M}(mathbf {r} )=-{int }_{mathbf {s} in Omega }operatorname {div} _{s}left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}Omega +{int }_{mathbf {s} in partial Omega }{frac {<mathbf {M} (mathbf {s} )|mathbf {n} _{partial Omega }>}{|mathbf {r} -mathbf {s} |}} |{d}^{2}(partial Omega )|} （2-3-8）

ここで、∂Ω{displaystyle partial Omega }partial Omega は、領域Ω{displaystyle Omega }Omega の境界面を意味する。
また、n∂Ω{displaystyle mathbf {n} _{partial Omega }}{displaystyle mathbf {n} _{partial Omega }}は、∂Ω{displaystyle partial Omega }partial Omega の法線ベクトルを意味する。


今、体積磁荷密度と、表面磁荷密度を、




  • ρM:=−divr⁡[M]{displaystyle {rho }_{M}:=-operatorname {div} _{mathbf {r} }[mathbf {M} ]}{displaystyle {rho }_{M}:=-operatorname {div} _{mathbf {r} }[mathbf {M} ]} （体積磁荷密度） （2-3-10）


  • σM:=<M | n∂Ω>{displaystyle {sigma }_{M}:=<mathbf {M} | mathbf {n} _{partial Omega }>}{displaystyle {sigma }_{M}:=<mathbf {M} | mathbf {n} _{partial Omega }>} （表面磁荷密度） （2-3-11）


により定めると、



ϕM(r)=∫s∈ΩρM(s)|r−s| d3Ω +∫s∈∂ΩσM(s)|r−s| |d2(∂Ω)|{displaystyle {phi }_{M}(mathbf {r} )={int }_{mathbf {s} in Omega }{frac {{rho }_{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}} {d}^{3}Omega +{int }_{mathbf {s} in partial Omega }{frac {{sigma }_{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}} |{d}^{2}(partial Omega )|}{displaystyle {phi }_{M}(mathbf {r} )={int }_{mathbf {s} in Omega }{frac {{rho }_{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}} {d}^{3}Omega +{int }_{mathbf {s} in partial Omega }{frac {{sigma }_{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}} |{d}^{2}(partial Omega )|} （2-3-12）

を得る。


一方、ϕM(r){displaystyle {phi }_{M}(mathbf {r} )}{displaystyle {phi }_{M}(mathbf {r} )}の定義により、



HM(r)=−gradr⁡[ϕM(r)]{displaystyle {mathbf {H} }_{M}(mathbf {r} )=-operatorname {grad} _{mathbf {r} }[{phi }_{M}(mathbf {r} )]}{displaystyle {mathbf {H} }_{M}(mathbf {r} )=-operatorname {grad} _{mathbf {r} }[{phi }_{M}(mathbf {r} )]} （2-3-13）

であるため、



HM(r)=∫s∈Ωgradr⁡[−ρM(s)|r−s|] d3Ω +∫s∈∂Ωgradr⁡[−σM(s)|r−s|] |d2(∂Ω)|{displaystyle {mathbf {H} }_{M}(mathbf {r} )={int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {{-rho }_{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}Omega +{int }_{mathbf {s} in partial Omega }operatorname {grad} _{mathbf {r} }left[{frac {{-sigma }_{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] |{d}^{2}(partial Omega )|}{displaystyle {mathbf {H} }_{M}(mathbf {r} )={int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {{-rho }_{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}Omega +{int }_{mathbf {s} in partial Omega }operatorname {grad} _{mathbf {r} }left[{frac {{-sigma }_{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] |{d}^{2}(partial Omega )|} （2-3-14）

である。



磁化電流の導入


時間的に定常な磁化が作り出す磁気ベクトルポテンシャルを、別の側面から考察してみることにしよう。ここでは、「磁束密度の原因は、電流に帰される」という思想に従い、
だとすれば、「磁化と等価な効果を発揮する電流」がどのようなものかを検討することにする。


前記AM(r){displaystyle mathbf {A} _{M}(mathbf {r} )}{displaystyle mathbf {A} _{M}(mathbf {r} )}に、
ベクトル解析の公式を適用すると、


AM(r)=μ04π∫s∈Ω((M(s))×(r−s)|r−s|) d3s=μ04π∫s∈Ω(rots⁡[M(s)]|r−s|) d3s+μ04π∫s∈∂Ω(M(s)×n∂Ω|r−s|) |d2(∂Ω)|(2-4-1){displaystyle {begin{aligned}mathbf {A} _{M}(mathbf {r} )&={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }left({frac {(mathbf {M} (mathbf {s} ))times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right) {d}^{3}mathbf {s} \&={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }left({frac {operatorname {rot} _{mathbf {s} }[mathbf {M} (mathbf {s} )]}{|mathbf {r} -mathbf {s} |}}right) {d}^{3}mathbf {s} +{frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {M} (mathbf {s} )times mathbf {n} _{partial Omega }}{|mathbf {r} -mathbf {s} |}}right) |{d}^{2}(partial Omega )|qquad {text{(2-4-1)}}end{aligned}}}{displaystyle {begin{aligned}mathbf {A} _{M}(mathbf {r} )&={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }left({frac {(mathbf {M} (mathbf {s} ))times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right) {d}^{3}mathbf {s} \&={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }left({frac {operatorname {rot} _{mathbf {s} }[mathbf {M} (mathbf {s} )]}{|mathbf {r} -mathbf {s} |}}right) {d}^{3}mathbf {s} +{frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {M} (mathbf {s} )times mathbf {n} _{partial Omega }}{|mathbf {r} -mathbf {s} |}}right) |{d}^{2}(partial Omega )|qquad {text{(2-4-1)}}end{aligned}}}

となることが判る。ここで、rots{displaystyle operatorname {rot} _{mathbf {s} }}{displaystyle operatorname {rot} _{mathbf {s} }}
は、変数s{displaystyle mathbf {s} }{mathbf {s}}についての回転微分を意味する。∂Ω{displaystyle partial Omega }partial Omega は、領域Ω{displaystyle Omega }Omega の境界を意味する。
上式の右辺第二項の積分において、「∂Ω{displaystyle partial Omega }partial Omega の面素の絶対値」を意味し
所謂普通の面積分ではない（この「面積分」は、「ベクトル場からベクトルを作る」特殊なもの）ので注意が必要である。
この積分において、n∂Ω{displaystyle mathbf {n} _{partial Omega }}{displaystyle mathbf {n} _{partial Omega }}は、∂Ω{displaystyle partial Omega }partial Omega
の法線ベクトルを意味する。



実際、スカラー倍の回転微分の公式より、


rots⁡[M(s)|r−s|]=(grads⁡[M(s)|r−s|])×M(s)+rots⁡[M(s)]|r−s|=(r−s)×M(s)|r−s|3+rots⁡[M(s)]|r−s|(2-4-2){displaystyle {begin{aligned}operatorname {rot} _{mathbf {s} }left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right]&=left(operatorname {grad} _{mathbf {s} }left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right]right)times mathbf {M} (mathbf {s} )+{frac {operatorname {rot} _{mathbf {s} }left[mathbf {M} (mathbf {s} )right]}{|mathbf {r} -mathbf {s} |}}\&={frac {(mathbf {r} -mathbf {s} )times mathbf {M} (mathbf {s} )}{{|mathbf {r} -mathbf {s} |}^{3}}}+{frac {operatorname {rot} _{mathbf {s} }left[mathbf {M} (mathbf {s} )right]}{|mathbf {r} -mathbf {s} |}}qquad {text{(2-4-2)}}end{aligned}}}{displaystyle {begin{aligned}operatorname {rot} _{mathbf {s} }left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right]&=left(operatorname {grad} _{mathbf {s} }left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right]right)times mathbf {M} (mathbf {s} )+{frac {operatorname {rot} _{mathbf {s} }left[mathbf {M} (mathbf {s} )right]}{|mathbf {r} -mathbf {s} |}}\&={frac {(mathbf {r} -mathbf {s} )times mathbf {M} (mathbf {s} )}{{|mathbf {r} -mathbf {s} |}^{3}}}+{frac {operatorname {rot} _{mathbf {s} }left[mathbf {M} (mathbf {s} )right]}{|mathbf {r} -mathbf {s} |}}qquad {text{(2-4-2)}}end{aligned}}}

従って、



M(s)×(r−s)|r−s|3=rots⁡[M(s)]|r−s|−rots⁡[M(s)|r−s|]{displaystyle {frac {mathbf {M} (mathbf {s} )times (mathbf {r} -mathbf {s} )}{{|mathbf {r} -mathbf {s} |}^{3}}}={frac {operatorname {rot} _{mathbf {s} }left[mathbf {M} (mathbf {s} )right]}{|mathbf {r} -mathbf {s} |}}-operatorname {rot} _{mathbf {s} }left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right]}{displaystyle {frac {mathbf {M} (mathbf {s} )times (mathbf {r} -mathbf {s} )}{{|mathbf {r} -mathbf {s} |}^{3}}}={frac {operatorname {rot} _{mathbf {s} }left[mathbf {M} (mathbf {s} )right]}{|mathbf {r} -mathbf {s} |}}-operatorname {rot} _{mathbf {s} }left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right]} （2-4-3）

従って、



AM(r)=μ04π∫s∈Ωrots⁡[M(s)]|r−s| d3s −μ04π∫s∈Ωrots⁡[M(s)|r−s|] d3s {displaystyle mathbf {A} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }{frac {operatorname {rot} _{mathbf {s} }left[mathbf {M} (mathbf {s} )right]}{|mathbf {r} -mathbf {s} |}} {d}^{3}mathbf {s} -{frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }operatorname {rot} _{mathbf {s} }left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}mathbf {s} }{displaystyle mathbf {A} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }{frac {operatorname {rot} _{mathbf {s} }left[mathbf {M} (mathbf {s} )right]}{|mathbf {r} -mathbf {s} |}} {d}^{3}mathbf {s} -{frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }operatorname {rot} _{mathbf {s} }left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}mathbf {s} }　　（2-4-4）

を得る。さらに、上式の右辺第二項の積分にベクトル解析の公式を適用すると、



−μ04π∫s∈Ωrots⁡[M(s)|r−s|] d3s =μ04π∫s∈∂Ω(M(s)|r−s|)×d2s {displaystyle -{frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }operatorname {rot} _{mathbf {s} }left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}mathbf {s} ={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right)times {d}^{2}mathbf {s} }{displaystyle -{frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }operatorname {rot} _{mathbf {s} }left[{frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right] {d}^{3}mathbf {s} ={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right)times {d}^{2}mathbf {s} }　　（2-4-5）

を得る。さらに、右辺にベクトル解析の公式を適用すると、



μ04π∫s∈∂Ω(M(s)|r−s|)×d2s =μ04π∫s∈∂Ω(M(s)×n∂Ω|r−s|) |d2(∂Ω)| {displaystyle {frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right)times {d}^{2}mathbf {s} ={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {M} (mathbf {s} )times mathbf {n} _{partial Omega }}{|mathbf {r} -mathbf {s} |}}right) |{d}^{2}(partial Omega )| }{displaystyle {frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {M} (mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right)times {d}^{2}mathbf {s} ={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {M} (mathbf {s} )times mathbf {n} _{partial Omega }}{|mathbf {r} -mathbf {s} |}}right) |{d}^{2}(partial Omega )| } （2-4-6）

が得られる。以上から示すべき式が証明された。



今、iM(r), KM{displaystyle {boldsymbol {i}}_{M}(mathbf {r} ), {mathbf {K} }_{M}}{displaystyle {boldsymbol {i}}_{M}(mathbf {r} ), {mathbf {K} }_{M}}を、




  • iM(r):=rot⁡[M(r)]{displaystyle {boldsymbol {i}}_{M}(mathbf {r} ):=operatorname {rot} [mathbf {M} (mathbf {r} )]}{displaystyle {boldsymbol {i}}_{M}(mathbf {r} ):=operatorname {rot} [mathbf {M} (mathbf {r} )]} （体積磁化電流密度） （2-4-7）


  • KM(r):=M(r)×n∂Ω(r){displaystyle {mathbf {K} }_{M}(mathbf {r} ):=mathbf {M} (mathbf {r} )times mathbf {n} _{partial Omega }(mathbf {r} )}{displaystyle {mathbf {K} }_{M}(mathbf {r} ):=mathbf {M} (mathbf {r} )times mathbf {n} _{partial Omega }(mathbf {r} )} (表面磁化電流密度) （2-4-8）


と置くと、結局、



AM(r)=μ04π∫s∈Ω(iM(s)|r−s|) d3s+ μ04π∫s∈∂Ω(KM(s)|r−s|) |d2(∂Ω)| {displaystyle mathbf {A} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }left({frac {{boldsymbol {i}}_{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right) {d}^{3}mathbf {s} + {frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {K} _{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right) |{d}^{2}(partial Omega )| }{displaystyle mathbf {A} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in Omega }left({frac {{boldsymbol {i}}_{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right) {d}^{3}mathbf {s} + {frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {K} _{M}(mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right) |{d}^{2}(partial Omega )| } （2-4-9）

となる。


両辺の回転微分を取ると、



BM(r)=μ04π∫s∈R3(iM(s)×(r−s)|r−s|3)d3s + μ04π∫s∈∂Ω(KM(s)×(r−s)|r−s|3) |d2(∂Ω)|{displaystyle mathbf {B} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}int _{mathbf {s} in mathbb {R} ^{3}}left({frac {{boldsymbol {i}}_{M}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right)d^{3}mathbf {s} + {frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {K} _{M}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right) |{d}^{2}(partial Omega )|}{displaystyle mathbf {B} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}int _{mathbf {s} in mathbb {R} ^{3}}left({frac {{boldsymbol {i}}_{M}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right)d^{3}mathbf {s} + {frac {mu _{0}}{4pi }}{int }_{mathbf {s} in partial Omega }left({frac {mathbf {K} _{M}(mathbf {s} )times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |^{3}}}right) |{d}^{2}(partial Omega )|} （2-4-10）

が、判る。この見方は、特に、磁化が一様な場合（より一般にはrot[M]=0）といった特殊な場合に、特に威力を発揮する。





強制電流と磁化の両方が既知のとき(一般論)



全系の磁束密度


強制電流ifc (forced current)と磁化Mの両方が既知とし、これら以外に磁束密度の原因となるものがないとした場合、
全系の磁束密度Btot
は、強制電流ifcに起因する磁束密度の成分（式(1-2)から求められる）
Bfcと、磁化Mが作り出す磁束密度の成分（式(2-2-8)から求められる）
BMによって、



Btot=Bfc+BM (3-1-1)

と表される。



ここで、よく注意しておかないといけないことは、「強制電流ifcと磁化Mの両方が既知」という言葉の意味であるが、
仮に、強制電流ifc（や、他の磁化）がない状態での磁化Miniが既知したとして、
(磁化Miniが作り出す磁束密度BMini がどうなるかは(2-2-8)から計算できるが)



Btot=Bfc+BMini (間違った式)

は、よほど特殊な場合を除き成り立たない。要は間違いである。



式(3-1-1)より正確に書くならば、「外場の影響等により磁化が変化した後の磁化
Mcon
が作り出す磁束密度」
BMconを用いて、



Btot=Bfc+BMcon (正しい式)

である。結局以下の難所は残ったままである。


（難所）元々磁化Miniを帯びている物質に、外部から磁束密度Bext(r){displaystyle mathbf {B} _{text{ext}}(mathbf {r} )}{displaystyle mathbf {B} _{text{ext}}(mathbf {r} )}を印加したとき、最終的に、物質がどのようなMconを得るか？

この問題が、実のところは難しい。つまり、強制電流が作った磁束密度や、磁化自身が作り出す磁束密度により、物質の磁化が最初の磁化から変化してしまうという問題がある。
一般には。B-H曲線等の実測結果と、上記の拘束条件を考え合わせ、数値計算によって磁束密度や磁化が計算されるのである。但し、線形物質に関して言えば、上の難所は比較的簡単である。このような特殊な物質に関する問題については、次章で述べることにする。


再び式(3-1-1)について考えよう。ベクトル解析の公式から




div⁡[Bfc]=0{displaystyle operatorname {div} [mathbf {B} _{fc}]=0}{displaystyle operatorname {div} [mathbf {B} _{fc}]=0} (3-1-2a)


div⁡[BM]=0{displaystyle operatorname {div} [mathbf {B} _{M}]=0}{displaystyle operatorname {div} [mathbf {B} _{M}]=0} (3-1-2b)


div⁡[Btot]=0{displaystyle operatorname {div} [mathbf {B} _{text{tot}}]=0}{displaystyle operatorname {div} [mathbf {B} _{text{tot}}]=0} (3-1-2c)


であることが判る。即ち、磁束保存の式が満たされることが判る。



B=μ0(H+M)について


前節同様に、強制電流ifcと磁化Mの両方が既知とし、これら以外に磁束密度の原因となるものがないとした場合
について考える。式(3-1-1)のBtotに対し、新たな場Htotを、



Htot(r):=Btot(r)μ0−M(r){displaystyle mathbf {H} _{text{tot}}(mathbf {r} ):={frac {mathbf {B} _{text{tot}}(mathbf {r} )}{{mu }_{0}}}-mathbf {M} (mathbf {r} )}{displaystyle mathbf {H} _{text{tot}}(mathbf {r} ):={frac {mathbf {B} _{text{tot}}(mathbf {r} )}{{mu }_{0}}}-mathbf {M} (mathbf {r} )} （3-2-1）

と定める。式(3-2-1)は、所謂”B=μ0(H+M)”に他ならない。




さらに、強制電流ifc{displaystyle {boldsymbol {i}}_{fc}}{displaystyle {boldsymbol {i}}_{fc}}が作る磁束密度Bfc(r){displaystyle mathbf {B} _{fc}(mathbf {r} )}{displaystyle mathbf {B} _{fc}(mathbf {r} )}は、式(1-2)に
より、（周辺の物質の有無にかかわらず）定まるが、式(1-3)に倣い
「強制電流ifc{displaystyle {boldsymbol {i}}_{fc}}{displaystyle {boldsymbol {i}}_{fc}}が作る磁場Hfc(r){displaystyle mathbf {H} _{fc}(mathbf {r} )}{displaystyle mathbf {H} _{fc}(mathbf {r} )}」を、



Hfc(r):=1μ0Bfc(r){displaystyle mathbf {H} _{fc}(mathbf {r} ):={frac {1}{{mu }_{0}}}mathbf {B} _{fc}(mathbf {r} )}{displaystyle mathbf {H} _{fc}(mathbf {r} ):={frac {1}{{mu }_{0}}}mathbf {B} _{fc}(mathbf {r} )} （3-2-2）

によって定め、磁化による磁場HMを式（2-2-10）のように定めると、
式（3-2-1）と式（3-2-2）、式（2-2-10）より、


Htot(r)=Btot(r)μ0−M(r)=Bfc(r)+BM(r)μ0−M(r)=Hfc(r)+(HM(r)+M(r))−M(r)=Hfc(r)+HM(r)(3-2-3){displaystyle {begin{aligned}mathbf {H} _{text{tot}}(mathbf {r} )&={frac {mathbf {B} _{text{tot}}(mathbf {r} )}{{mu }_{0}}}-mathbf {M} (mathbf {r} )\&={frac {mathbf {B} _{fc}(mathbf {r} )+mathbf {B} _{M}(mathbf {r} )}{{mu }_{0}}}-mathbf {M} (mathbf {r} )\&=mathbf {H} _{fc}(mathbf {r} )+(mathbf {H} _{M}(mathbf {r} )+mathbf {M} (mathbf {r} ))-mathbf {M} (mathbf {r} )\&=mathbf {H} _{fc}(mathbf {r} )+mathbf {H} _{M}(mathbf {r} )qquad {text{(3-2-3)}}end{aligned}}}{displaystyle {begin{aligned}mathbf {H} _{text{tot}}(mathbf {r} )&={frac {mathbf {B} _{text{tot}}(mathbf {r} )}{{mu }_{0}}}-mathbf {M} (mathbf {r} )\&={frac {mathbf {B} _{fc}(mathbf {r} )+mathbf {B} _{M}(mathbf {r} )}{{mu }_{0}}}-mathbf {M} (mathbf {r} )\&=mathbf {H} _{fc}(mathbf {r} )+(mathbf {H} _{M}(mathbf {r} )+mathbf {M} (mathbf {r} ))-mathbf {M} (mathbf {r} )\&=mathbf {H} _{fc}(mathbf {r} )+mathbf {H} _{M}(mathbf {r} )qquad {text{(3-2-3)}}end{aligned}}}

が得られる。従って、



rot⁡[Htot(r)]=rot⁡[Hfc(r)]+rot⁡[HM(r)]=ifc{displaystyle operatorname {rot} [mathbf {H} _{text{tot}}(mathbf {r} )]=operatorname {rot} [mathbf {H} _{fc}(mathbf {r} )]+operatorname {rot} [mathbf {H} _{M}(mathbf {r} )]={boldsymbol {i}}_{fc}}{displaystyle operatorname {rot} [mathbf {H} _{text{tot}}(mathbf {r} )]=operatorname {rot} [mathbf {H} _{fc}(mathbf {r} )]+operatorname {rot} [mathbf {H} _{M}(mathbf {r} )]={boldsymbol {i}}_{fc}} （3-2-4）

が得られる。これは、即ち（静磁場の）アンペールの法則である。



式（3-2-4）を示そう。式（2-2-9） 及び以下の式(3-2-5)より、


rot grad=0 (3-2-5)

HM(r){displaystyle mathbf {H} _{M}(mathbf {r} )}{displaystyle mathbf {H} _{M}(mathbf {r} )}に対し、回転微分を作用させると、



rot⁡[HM(r)]=14π∫s∈Ωrotr⁡[gradr⁡[<M(s)|(r−s)>|r−s|3]] d3s =0{displaystyle operatorname {rot} [mathbf {H} _{M}(mathbf {r} )]={frac {1}{4pi }}{int }_{mathbf {s} in Omega }operatorname {rot} _{mathbf {r} }left[operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]right] {d}^{3}{s} =0}{displaystyle operatorname {rot} [mathbf {H} _{M}(mathbf {r} )]={frac {1}{4pi }}{int }_{mathbf {s} in Omega }operatorname {rot} _{mathbf {r} }left[operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]right] {d}^{3}{s} =0} (3-2-6)

が判る。従って、磁化に起因する磁場HM(r){displaystyle mathbf {H} _{M}(mathbf {r} )}{displaystyle mathbf {H} _{M}(mathbf {r} )}は、



rot⁡[HM(r)]=0{displaystyle operatorname {rot} [mathbf {H} _{M}(mathbf {r} )]=0}{displaystyle operatorname {rot} [mathbf {H} _{M}(mathbf {r} )]=0} (3-2-7)

を満たす。





式(3-2-7)と式(1-5)より、式（3-2-4）が示された。



透磁率の導入


透磁率なる概念を、新たに導入する。
即ち、3次正方行列に値を取る行列値関数μ(r){displaystyle {mu }(mathbf {r} )}{displaystyle {mu }(mathbf {r} )}（透磁率）を用いて、



Btot(r)=μ(r)Htot(r){displaystyle mathbf {B} _{text{tot}}(mathbf {r} )={mu }(mathbf {r} )mathbf {H} _{text{tot}}(mathbf {r} )}{displaystyle mathbf {B} _{text{tot}}(mathbf {r} )={mu }(mathbf {r} )mathbf {H} _{text{tot}}(mathbf {r} )} (3-3-1)

と書けるものとする。この透磁率は、磁化の概念の”すり替え”に過ぎない概念である。
尚、通常は、透磁率は、スカラー値関数と考えてよい場合が多いのだが、その場合は、
「対角成分が全部同じ値で、それ以外の成分が0の行列値関数と、スカラー値関数が同一視できる」ことを思い起こせばいい。


さらに、式(3-3-1) で見たような、透磁率（透磁率テンソル）を用いて、



Btot(r)=μ(r)Htot(r){displaystyle mathbf {B} _{text{tot}}(mathbf {r} )={mu }(mathbf {r} )mathbf {H} _{text{tot}}(mathbf {r} )}{displaystyle mathbf {B} _{text{tot}}(mathbf {r} )={mu }(mathbf {r} )mathbf {H} _{text{tot}}(mathbf {r} )} (3-3-2)

と書き表せ、かつ、μが全点で正則行列（逆行列を持つ）とする。(そう考えても、”あまり”一般性を失わない）逆に言えば、



Htot(r)=ν(r)Btot(r){displaystyle mathbf {H} _{text{tot}}(mathbf {r} )={nu }(mathbf {r} )mathbf {B} _{text{tot}}(mathbf {r} )}{displaystyle mathbf {H} _{text{tot}}(mathbf {r} )={nu }(mathbf {r} )mathbf {B} _{text{tot}}(mathbf {r} )} (3-3-3)

である。ここで、


ν:=μ-1 (3-3-4)

は、各点でμの逆行列を与えるような行列値関数である。即ち、


μ(s) ν(s)=μ(s) μ-1(s)=E3 (3-3-5)

を充たすような行列値関数である。（ここでE3は3次の単位行列)である。




これまでの議論では、電流素片や磁気双極子モーメントが作り出した磁束密度/磁場について論じてきたが、
式(3-1-2),式（3-2-4）より、このようにして作られた全系の磁場が、磁束保存の式と、アンペールの法則を
充たすことが判った。
ここからは、逆に、磁束保存の式(3-3-6)と、アンペールの法則(3-3-7)、即ち、




div⁡[Bfc]=0{displaystyle operatorname {div} [mathbf {B} _{fc}]=0}{displaystyle operatorname {div} [mathbf {B} _{fc}]=0} (3-3-6a)


div⁡[BM]=0{displaystyle operatorname {div} [mathbf {B} _{M}]=0}{displaystyle operatorname {div} [mathbf {B} _{M}]=0} (3-3-6b)


div⁡[Btot]=0{displaystyle operatorname {div} [mathbf {B} _{text{tot}}]=0}{displaystyle operatorname {div} [mathbf {B} _{text{tot}}]=0} (3-3-6c)


rot⁡[Htot]=ifc{displaystyle operatorname {rot} [mathbf {H} _{text{tot}}]={boldsymbol {i}}_{fc}}{displaystyle operatorname {rot} [mathbf {H} _{text{tot}}]={boldsymbol {i}}_{fc}} (3-3-7)


を出発点とし、全系の磁気ベクトルポテンシャルAtot{displaystyle mathbf {A} _{text{tot}}}{displaystyle mathbf {A} _{text{tot}}}が充たす微分方程式を導出する。


まず、磁束保存の式 (3-3-6c)は、強制電流に起因する成分 (3-3-6a)、磁化に起因する成分 (3-3-6b)
それぞれについて成り立つため、（ポアンカレの補助定理より）それぞれがベクトルポテンシャルを持ち、




rot⁡[Afc]=Bfc{displaystyle operatorname {rot} [mathbf {A} _{fc}]=mathbf {B} _{fc}}{displaystyle operatorname {rot} [mathbf {A} _{fc}]=mathbf {B} _{fc}} (3-3-8a)


rot⁡[AM]=BM{displaystyle operatorname {rot} [mathbf {A} _{M}]=mathbf {B} _{M}}{displaystyle operatorname {rot} [mathbf {A} _{M}]=mathbf {B} _{M}} (3-3-8b)


をみたすようなベクトル場Afc{displaystyle mathbf {A} _{fc}}{displaystyle mathbf {A} _{fc}}と、AM{displaystyle mathbf {A} _{M}}{displaystyle mathbf {A} _{M}}が存在する。
このようなベクトル場Afc{displaystyle mathbf {A} _{fc}}{displaystyle mathbf {A} _{fc}}と、AM{displaystyle mathbf {A} _{M}}{displaystyle mathbf {A} _{M}}は、ゲージ不定性を除き一意に定まるが、
本記事では、強制電流に起因する成分Afc{displaystyle mathbf {A} _{fc}}{displaystyle mathbf {A} _{fc}}は、式(1-1)のものを採用し、
と、磁化に起因する成分AM{displaystyle mathbf {A} _{M}}{displaystyle mathbf {A} _{M}}式(2-1-5）を採用することにする。


全系のベクトルポテンシャルを、



Atot:=AM+Afc{displaystyle mathbf {A} _{text{tot}}:=mathbf {A} _{M}+mathbf {A} _{fc}}{displaystyle mathbf {A} _{text{tot}}:=mathbf {A} _{M}+mathbf {A} _{fc}} (3-3-9)

と定めると、



rot⁡[Atot]=Btot{displaystyle operatorname {rot} [mathbf {A} _{text{tot}}]=mathbf {B} _{text{tot}}}{displaystyle operatorname {rot} [mathbf {A} _{text{tot}}]=mathbf {B} _{text{tot}}} (3-3-8c)

が成り立つ。これは、(3-3-6c)の、磁束保存の式をベクトルポテンシャルを用いて書いたものに他ならない。


さらに、式(3-3-8c)に、式 (3-3-2)と、アンペールの法則(3-3-7)を考え併せると、



rot⁡[νrot⁡[A]]−ifc=0{displaystyle operatorname {rot} [{nu }operatorname {rot} [mathbf {A} ]]-{boldsymbol {i}}_{fc}=0}{displaystyle operatorname {rot} [{nu }operatorname {rot} [mathbf {A} ]]-{boldsymbol {i}}_{fc}=0} (3-3-10)

が得られ、未知のベクトル場BやHが消え、未知のベクトル場はAのみとなる。これを
「ベクトルポテンシャルによる静磁場の方程式」（静磁場の支配方程式）という。



静磁場の境界条件


物質の境界において、式(3-3-6a)、式 (3-3-7)に、ガウスの発散定理や、ケルビンストークスの定理を適用すると、
磁場および磁束密度の境界条件が得られる。





静磁場のエネルギーと停留値関数


前節の前提条件において、
全系の磁気ベクトルポテンシャルAtot{displaystyle mathbf {A} _{text{tot}}}{displaystyle mathbf {A} _{text{tot}}}は、以下の汎関数Fの停留関数となることが判る。[4][8]


F[A]:=∫s∈R312< H(s) | B(s)>−< ifc(s)| Atot(s)> ds=∫s∈R312< ν(s)(rots⁡[A(s)]) | (rots⁡[A(s)])>−< ifc(s)| A(s)> ds(3-5-1){displaystyle {begin{aligned}F[mathbf {A} ]&:={int }_{mathbf {s} in mathbb {R} ^{3}}{frac {1}{2}}< mathbf {H} (mathbf {s} ) | mathbf {B} (mathbf {s} )>-< {boldsymbol {i}}_{fc}(mathbf {s} )| mathbf {A} _{text{tot}}(mathbf {s} )> dmathbf {s} \&={int }_{mathbf {s} in mathbb {R} ^{3}}{frac {1}{2}}< {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )])>-< {boldsymbol {i}}_{fc}(mathbf {s} )| mathbf {A} (mathbf {s} )> dmathbf {s} qquad {text{(3-5-1)}}end{aligned}}}{displaystyle {begin{aligned}F[mathbf {A} ]&:={int }_{mathbf {s} in mathbb {R} ^{3}}{frac {1}{2}}< mathbf {H} (mathbf {s} ) | mathbf {B} (mathbf {s} )>-< {boldsymbol {i}}_{fc}(mathbf {s} )| mathbf {A} _{text{tot}}(mathbf {s} )> dmathbf {s} \&={int }_{mathbf {s} in mathbb {R} ^{3}}{frac {1}{2}}< {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )])>-< {boldsymbol {i}}_{fc}(mathbf {s} )| mathbf {A} (mathbf {s} )> dmathbf {s} qquad {text{(3-5-1)}}end{aligned}}}

次元解析をすると、式(3-5-1)は、エネルギーの次元を持つことが判るが、実際に、式(3-5-1)は、全系のエネルギーとなっている。
ここで、< | >{displaystyle < | >}{displaystyle < | >}は、内積を表す。νの定義は、式(3-3-4)に記載のとおりである。




式(3-5-1)が、本当に、「磁気ベクトルポテンシャルの停留汎関数」ことを検証しよう。
Aに対し、微小な摂動δA[注釈 8]
を与えた際の第一変分δF[注釈 8]、即ち、


δF=F[A+δA]-F[A] (3-5-2)

を求める。


F[A+δA]は、以下の被積分関数を、全空間でｓについて積分したものである。



12< ν(s)(rots⁡[A(s)+δA(s)]) | (rots⁡[A(s)+δA(s)])> −< ifc(s)| A(s)+δA(s)>{displaystyle {frac {1}{2}}< {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )+delta mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )+delta mathbf {A} (mathbf {s} )])> -< {boldsymbol {i}}_{fc}(mathbf {s} )| mathbf {A} (mathbf {s} )+delta mathbf {A} (mathbf {s} )>}{displaystyle {frac {1}{2}}< {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )+delta mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )+delta mathbf {A} (mathbf {s} )])> -< {boldsymbol {i}}_{fc}(mathbf {s} )| mathbf {A} (mathbf {s} )+delta mathbf {A} (mathbf {s} )>} (3-5-3)

一方、F[A]は、以下の被積分関数を、全空間でｓについて積分したものである。



12< ν(s)(rots⁡[A(s)]) | (rots⁡[A(s)])> −< ifc(s)| A(s)>{displaystyle {frac {1}{2}}< {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )])> -< {boldsymbol {i}}_{fc}(mathbf {s} )| mathbf {A} (mathbf {s} )>}{displaystyle {frac {1}{2}}< {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )])> -< {boldsymbol {i}}_{fc}(mathbf {s} )| mathbf {A} (mathbf {s} )>} (3-5-4)

したがって、δFは、以下の(3-5-5)を、全空間でｓについて積分したものである。



< ν(s)(rots⁡[A(s)]) | (rots⁡[δA(s)])>+12< ν(s)(rots⁡[δA(s)]) | (rots⁡[δA(s)])> −< ifc(s)| δA(s)>{displaystyle < {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[delta mathbf {A} (mathbf {s} )])>+{frac {1}{2}}< {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[delta mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[delta mathbf {A} (mathbf {s} )])> -< {boldsymbol {i}}_{fc}(mathbf {s} )| delta mathbf {A} (mathbf {s} )>}{displaystyle < {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[delta mathbf {A} (mathbf {s} )])>+{frac {1}{2}}< {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[delta mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[delta mathbf {A} (mathbf {s} )])> -< {boldsymbol {i}}_{fc}(mathbf {s} )| delta mathbf {A} (mathbf {s} )>} (3-5-5)

式(3-5-5)から、二次の微小項を無視すると、Fの第”一”変分



δF[A]=∫s∈R3< ν(s)(rots⁡[A(s)]) | (rots⁡[δA(s)])> −< ifc(s)| δA(s)> d3s{displaystyle delta F[A]={int }_{mathbf {s} in mathbb {R} ^{3}}< {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[delta mathbf {A} (mathbf {s} )])> -< {boldsymbol {i}}_{fc}(mathbf {s} )| delta mathbf {A} (mathbf {s} )> {d}^{3}mathbf {s} }{displaystyle delta F[A]={int }_{mathbf {s} in mathbb {R} ^{3}}< {nu }(mathbf {s} )(operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )]) | (operatorname {rot} _{mathbf {s} }[delta mathbf {A} (mathbf {s} )])> -< {boldsymbol {i}}_{fc}(mathbf {s} )| delta mathbf {A} (mathbf {s} )> {d}^{3}mathbf {s} } (3-5-6)

を得る。
式(3-5-6) に、ベクトル解析の公式を適用する。



div⁡[X×Y]=< Y | rot⁡[X] >−< X | rot⁡[Y] >{displaystyle operatorname {div} [mathbf {X} times mathbf {Y} ]=< mathbf {Y} | operatorname {rot} [mathbf {X} ] >-< mathbf {X} | operatorname {rot} [mathbf {Y} ] >}{displaystyle operatorname {div} [mathbf {X} times mathbf {Y} ]=< mathbf {Y} | operatorname {rot} [mathbf {X} ] >-< mathbf {X} | operatorname {rot} [mathbf {Y} ] >} (3-5-7)

において


X:=ν(s)rot⁡[A(s)](3-5-8)Y:=δA(3-5-9){displaystyle {begin{aligned}mathbf {X} &:={nu }(mathbf {s} )operatorname {rot} [mathbf {A} (mathbf {s} )]&{text{(3-5-8)}}\mathbf {Y} &:=delta mathbf {A} &{text{(3-5-9)}}end{aligned}}}{displaystyle {begin{aligned}mathbf {X} &:={nu }(mathbf {s} )operatorname {rot} [mathbf {A} (mathbf {s} )]&{text{(3-5-8)}}\mathbf {Y} &:=delta mathbf {A} &{text{(3-5-9)}}end{aligned}}}

を代入すると、


δF[A]=∫s∈R3−divs⁡[(ν(s)rots⁡[A])×δA] d3s+∫s∈R3<(rots⁡[(ν(s))rots⁡[A(s)] ]−ifc(s) | δA(s)> d3s(3-5-10){displaystyle {begin{aligned}delta F[mathbf {A} ]&={int }_{mathbf {s} in mathbb {R} ^{3}}-operatorname {div} _{mathbf {s} }[({nu }(mathbf {s} )operatorname {rot} _{mathbf {s} }[mathbf {A} ])times delta mathbf {A} ] {d}^{3}mathbf {s} \&+{int }_{mathbf {s} in mathbb {R} ^{3}}<(operatorname {rot} _{mathbf {s} }[({nu }(mathbf {s} ))operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )] ]-{boldsymbol {i}}_{fc}(mathbf {s} ) | delta mathbf {A} (mathbf {s} )> {d}^{3}mathbf {s} qquad {text{(3-5-10)}}end{aligned}}}{displaystyle {begin{aligned}delta F[mathbf {A} ]&={int }_{mathbf {s} in mathbb {R} ^{3}}-operatorname {div} _{mathbf {s} }[({nu }(mathbf {s} )operatorname {rot} _{mathbf {s} }[mathbf {A} ])times delta mathbf {A} ] {d}^{3}mathbf {s} \&+{int }_{mathbf {s} in mathbb {R} ^{3}}<(operatorname {rot} _{mathbf {s} }[({nu }(mathbf {s} ))operatorname {rot} _{mathbf {s} }[mathbf {A} (mathbf {s} )] ]-{boldsymbol {i}}_{fc}(mathbf {s} ) | delta mathbf {A} (mathbf {s} )> {d}^{3}mathbf {s} qquad {text{(3-5-10)}}end{aligned}}}

式(3-5-10)の第一項については、式(3-4-6)と、「真空中ではdiv[H]=0であること」を考え併せると、
結局物質Ω内の効果しか寄与しないことが判る。さらにガウスの発散定理を考慮すると、式(3-5-10)の第一項は、



∫s∈Ω−divs⁡[(ν(s)rots⁡[A])×δA ] d3s=−∫s∈∂Ω(ν(s)rots⁡[A])×δA d2(∂Ω){displaystyle {int }_{mathbf {s} in Omega }-operatorname {div} _{mathbf {s} }[({nu }(mathbf {s} )operatorname {rot} _{mathbf {s} }[mathbf {A} ])times delta mathbf {A} ] {d}^{3}mathbf {s} =-{int }_{mathbf {s} in partial Omega }({nu }(mathbf {s} )operatorname {rot} _{mathbf {s} }[mathbf {A} ])times delta mathbf {A} {d}^{2}(partial Omega )}{displaystyle {int }_{mathbf {s} in Omega }-operatorname {div} _{mathbf {s} }[({nu }(mathbf {s} )operatorname {rot} _{mathbf {s} }[mathbf {A} ])times delta mathbf {A} ] {d}^{3}mathbf {s} =-{int }_{mathbf {s} in partial Omega }({nu }(mathbf {s} )operatorname {rot} _{mathbf {s} }[mathbf {A} ])times delta mathbf {A} {d}^{2}(partial Omega )} (3-5-11)

となる。従って、式(3-5-10)の第一項が、任意の摂動δAに対して0になるためには、ノイマン条件即ち、



(rot⁡[A])×n∂Ω = 0{displaystyle (operatorname {rot} [mathbf {A} ])times {n}_{partial Omega } = 0}{displaystyle (operatorname {rot} [mathbf {A} ])times {n}_{partial Omega } = 0} (3-5-12)

が満たされればよい。さらに、式(3-5-10)の第二項が、任意の摂動δAに対して0になるためには、



rot⁡[νrot⁡[A]]−ifc=0{displaystyle operatorname {rot} [{nu }operatorname {rot} [mathbf {A} ]]-{boldsymbol {i}}_{fc}=0}{displaystyle operatorname {rot} [{nu }operatorname {rot} [mathbf {A} ]]-{boldsymbol {i}}_{fc}=0} (3-5-12)

であればよい。これは、「ベクトルポテンシャルによる静磁場の方程式」（静磁場の支配方程式）に(式(3-3-10))他ならない。以上から、ノイマン条件下での、静磁場の支配方程式の解は、汎関数Fの停留関数となること判る。



ベクトル解析に関する補足


静磁場の解析には、静磁場の解析には、やや高度なベクトル解析の知識
(例えば 藤本[33]参照）
が要求される。同じ式変形が何度も繰り返されるのを避ける目的から、
一般的な知識は既知とする立場から、本記事の内容理解に必要と思われる事柄に限り、簡単に説明する。
が要求される。同じ式変形が何度も繰り返されるのを避ける目的から、
一般的な知識は既知とする立場から、本記事の内容理解に必要と思われる事柄に限り、簡単に説明する。



ベクトル場の代数演算と微分作用素


ベクトル場に代数演算を施したものに、微分作用素を作用させた場合に成り立つ公式について、本記事で用い、かつ、あまり本に載っていないものについて、簡潔にまとめる。



(1) <F|∇>について


F=(f1,f2,f3)をベクトル場とする。このとき、


⟨ F | ∇ ⟩:=f1∂∂x1+f2∂∂x2+f3∂∂x3{displaystyle leftlangle mathbf {F} | nabla rightrangle :={f}_{1}{partial over partial {x}_{1}}+{f}_{2}{partial over partial {x}_{2}}+{f}_{3}{partial over partial {x}_{3}}}{displaystyle leftlangle mathbf {F} | nabla rightrangle :={f}_{1}{partial over partial {x}_{1}}+{f}_{2}{partial over partial {x}_{2}}+{f}_{3}{partial over partial {x}_{3}}}

と定義する。ここで、∇は、


∇:=x^∂∂x+y^∂∂y+z^∂∂z{displaystyle nabla :=mathbf {hat {x}} {partial over partial x}+mathbf {hat {y}} {partial over partial y}+mathbf {hat {z}} {partial over partial z}}nabla :={mathbf {{hat {x}}}}{partial over partial x}+{mathbf {{hat {y}}}}{partial over partial y}+{mathbf {{hat {z}}}}{partial over partial z}

を意味する。<F|∇>のことを、”F・∇”と書くこともある。


Gを、ベクトル場としたとき、<F|∇>をGに作用させると、


⟨ F | ∇ ⟩G=f1∂g1∂x1+f2∂g2∂x2+f3∂g3∂x3=(J[G])⋅F{displaystyle leftlangle mathbf {F} | nabla rightrangle mathbf {G} ={f}_{1}{partial {g}_{1} over partial {x}_{1}}+{f}_{2}{partial {g}_{2} over partial {x}_{2}}+{f}_{3}{partial {g}_{3} over partial {x}_{3}}=(J[mathbf {G} ])cdot mathbf {F} }{displaystyle leftlangle mathbf {F} | nabla rightrangle mathbf {G} ={f}_{1}{partial {g}_{1} over partial {x}_{1}}+{f}_{2}{partial {g}_{2} over partial {x}_{2}}+{f}_{3}{partial {g}_{3} over partial {x}_{3}}=(J[mathbf {G} ])cdot mathbf {F} }

が成り立つ。ここで、J[G]{displaystyle J[mathbf {G} ]}{displaystyle J[mathbf {G} ]}は、Gのヤコビ行列（ヤコビアンではない）を意味する。



(2)rot とスカラー倍、ベクトル積


F,Gを、ベクトル場、fを、スカラー値関数とする。このとき、以下が成り立つ。


rot⁡[fF]=grad⁡[f]×F+f⋅rot⁡[F]{displaystyle operatorname {rot} [fmathbf {F} ]=operatorname {grad} [f]times mathbf {F} +fcdot operatorname {rot} [mathbf {F} ]}{displaystyle operatorname {rot} [fmathbf {F} ]=operatorname {grad} [f]times mathbf {F} +fcdot operatorname {rot} [mathbf {F} ]}

rot⁡[F×G]=⟨ G | ∇ ⟩F−⟨ F | ∇ ⟩G+Fdiv⁡[G]−Gdiv⁡[F]{displaystyle operatorname {rot} [mathbf {F} times mathbf {G} ]=leftlangle mathbf {G} | nabla rightrangle mathbf {F} -leftlangle mathbf {F} | nabla rightrangle mathbf {G} +mathbf {F} operatorname {div} [mathbf {G} ]-mathbf {G} operatorname {div} [mathbf {F} ]}{displaystyle operatorname {rot} [mathbf {F} times mathbf {G} ]=leftlangle mathbf {G} | nabla rightrangle mathbf {F} -leftlangle mathbf {F} | nabla rightrangle mathbf {G} +mathbf {F} operatorname {div} [mathbf {G} ]-mathbf {G} operatorname {div} [mathbf {F} ]}


(3)div とスカラー倍、ベクトル積


F,Gを、ベクトル場、fを、スカラー値関数とする。このとき、


div⁡[fF]=⟨F|grad⁡[f]⟩+fdiv⁡[F]{displaystyle operatorname {div} left[fmathbf {F} right]=leftlangle mathbf {F} |operatorname {grad} [f]rightrangle +foperatorname {div} [mathbf {F} ]}{displaystyle operatorname {div} left[fmathbf {F} right]=leftlangle mathbf {F} |operatorname {grad} [f]rightrangle +foperatorname {div} [mathbf {F} ]}


ディラックのデルタと、重積分


一般のn変数関数（関数はベクトル値関数であってよい)fの定義域をΩとしたとき、


f(r)=∫s∈Ωf(s)δn(r−s) dns =f∗δn{displaystyle f(mathbf {r} )={int }_{mathbf {s} in Omega }f(mathbf {s} ){delta }^{n}(mathbf {r} -mathbf {s} ) {d}^{n}s =f*{delta }^{n}}f({mathbf {r}})={int }_{{{mathbf {s}}in Omega }}f({mathbf {s}}){delta }^{{n}}({mathbf {r}}-{mathbf {s}}) {d}^{{n}}s =f*{delta }^{{n}}

であることが知られている。ここで、上式の"*"は、合成積(乗算ではない)である。また、δnは、n変数のδ関数である。



発散微分とディラックのデルタ


原点(r=0)を除いて、



divr⁡[r|r|]=3|r|3−3|r|2|r|5=0{displaystyle operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} }{|mathbf {r} |}}right]={frac {3}{|mathbf {r} |^{3}}}-{frac {3{|mathbf {r} |}^{2}}{{|mathbf {r} |}^{5}}}=0}{displaystyle operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} }{|mathbf {r} |}}right]={frac {3}{|mathbf {r} |^{3}}}-{frac {3{|mathbf {r} |}^{2}}{{|mathbf {r} |}^{5}}}=0}

であり、原点(r=0)を中心とする、球体BLに対し、ガウスの発散定理を用いると、
∫r∈BLdivr⁡[r|r|] d3r=4π{displaystyle {int }_{mathbf {r} in BL}operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} }{|mathbf {r} |}}right] d^{3}r=4pi }{displaystyle {int }_{mathbf {r} in BL}operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} }{|mathbf {r} |}}right] d^{3}r=4pi }


となるので、結局、
divr⁡[r|r|]=4πδ(r){displaystyle operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} }{|mathbf {r} |}}right]=4pi delta (mathbf {r} )}{displaystyle operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} }{|mathbf {r} |}}right]=4pi delta (mathbf {r} )}



であることが判る。(ここで、δは、3変数のδ関数である。

以上の議論を平行移動させると、


divr⁡[r−s|r−s|]=4πδ(r−s){displaystyle operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]=4pi delta (mathbf {r} -mathbf {s} )}{displaystyle operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]=4pi delta (mathbf {r} -mathbf {s} )}

であることが判る。



分数関数の偏微分


∂∂si[1|r−s|]=∂∂si[1((ri−si)2+const)1/2]=∂∂si[((ri−si)2+const)−(1/2)]{displaystyle {begin{aligned}{frac {partial }{partial mathbf {s} _{i}}}left[{frac {1}{|mathbf {r} -mathbf {s} |}}right]&={frac {partial }{partial {s}_{i}}}left[{frac {1}{{(({r}_{i}-{s}_{i})^{2}+{text{const}})}^{1/2}}}right]\&={frac {partial }{partial {s}_{i}}}left[{(({r}_{i}-{s}_{i})^{2}+{text{const}})}^{-(1/2)}right]end{aligned}}}{displaystyle {begin{aligned}{frac {partial }{partial mathbf {s} _{i}}}left[{frac {1}{|mathbf {r} -mathbf {s} |}}right]&={frac {partial }{partial {s}_{i}}}left[{frac {1}{{(({r}_{i}-{s}_{i})^{2}+{text{const}})}^{1/2}}}right]\&={frac {partial }{partial {s}_{i}}}left[{(({r}_{i}-{s}_{i})^{2}+{text{const}})}^{-(1/2)}right]end{aligned}}}


面積分に関する補足


ここでは、様々なベクトル場の面積分についてまとめる。
「絶対値による面積分」「外積面積分」という用語は一般的な用語ではないが、他に
適切な表現がないため、この場限りでそのような言い方をする。
本記事内での定義は、それぞれ以下の通り。


IをR2{displaystyle mathbb {R} ^{2}}{mathbb {R}}^{{2}}の閉長方形、
φ:I→R3{displaystyle varphi :Ito mathbb {R} ^{3}}{displaystyle varphi :Ito mathbb {R} ^{3}}は、
Iの近傍で、殆ど至る所区分的に滑らかかつ、非退化であり、かつ、Iの内部で単射なベクトル値関数、
S:=φ[I]{displaystyle S:=varphi [I]}{displaystyle S:=varphi [I]}を、R3{displaystyle mathbb {R} ^{3}}{mathbb {R}}^{{3}}内の曲面片とする。

また、Xを、Sの近傍で定義された、区分的に滑らかなベクトル場とする。このとき、



(1)普通の面積分


定義は、以下の通りである。
∫s∈SX(s) d2S:=∫(u1,u2)∈I⟨X(φ(u1,u2)) | (∂φ∂u1×∂φ∂u2)⟩ du1du2{displaystyle {int }_{mathbf {s} in S}mathbf {X} (mathbf {s} ) {d}^{2}S:={int }_{({u}_{1},{u}_{2})in I}langle mathbf {X} (varphi ({u}_{1},{u}_{2})) | left({frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right)rangle d{u}_{1}d{u}_{2}}{displaystyle {int }_{mathbf {s} in S}mathbf {X} (mathbf {s} ) {d}^{2}S:={int }_{({u}_{1},{u}_{2})in I}langle mathbf {X} (varphi ({u}_{1},{u}_{2})) | left({frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right)rangle d{u}_{1}d{u}_{2}}



上式の右辺は、(u1,u2){displaystyle ({u}_{1},{u}_{2})}{displaystyle ({u}_{1},{u}_{2})}についてのスカラー値関数
⟨X(φ(u1,u2)) | (∂φ∂u1×∂φ∂u2)⟩{displaystyle langle mathbf {X} (varphi ({u}_{1},{u}_{2})) | left({frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right)rangle }{displaystyle langle mathbf {X} (varphi ({u}_{1},{u}_{2})) | left({frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right)rangle }


を、区間I上で重積分したものを意味する。今、Sの単位法線ベクトルnS{displaystyle {mathbf {n} }_{S}}{displaystyle {mathbf {n} }_{S}}を、
nS:=1|∂φ∂u1×∂φ∂u2|(∂φ∂u1×∂φ∂u2){displaystyle {mathbf {n} }_{S}:={frac {1}{left|{frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right|}}left({frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right)}{displaystyle {mathbf {n} }_{S}:={frac {1}{left|{frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right|}}left({frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right)}


と定めると、



∫s∈SX(s) d2S=∫(u1,u2)∈I⟨X(φ(u1,u2)) | nS⟩|∂φ∂u1×∂φ∂u2| du1du2{displaystyle {int }_{mathbf {s} in S}mathbf {X} (mathbf {s} ) {d}^{2}S={int }_{({u}_{1},{u}_{2})in I}langle mathbf {X} (varphi ({u}_{1},{u}_{2})) | {mathbf {n} }_{S}rangle left|{frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right| d{u}_{1}d{u}_{2}}{displaystyle {int }_{mathbf {s} in S}mathbf {X} (mathbf {s} ) {d}^{2}S={int }_{({u}_{1},{u}_{2})in I}langle mathbf {X} (varphi ({u}_{1},{u}_{2})) | {mathbf {n} }_{S}rangle left|{frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right| d{u}_{1}d{u}_{2}}

である。



(2)絶対値による面積分


定義は、以下の通りである。
∫s∈SX(s) |d2S|:=∫(u1,u2)∈IX(φ(u1,u2)) |∂φ∂u1×∂φ∂u2| du1du2{displaystyle {int }_{mathbf {s} in S}mathbf {X} (mathbf {s} ) |{d}^{2}S|:={int }_{({u}_{1},{u}_{2})in I}mathbf {X} (varphi ({u}_{1},{u}_{2})) left|{frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right| d{u}_{1}d{u}_{2}}{displaystyle {int }_{mathbf {s} in S}mathbf {X} (mathbf {s} ) |{d}^{2}S|:={int }_{({u}_{1},{u}_{2})in I}mathbf {X} (varphi ({u}_{1},{u}_{2})) left|{frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right| d{u}_{1}d{u}_{2}}



本記事では、絶対値による面積分の場合は、 |d2S|{displaystyle |{d}^{2}S|}{displaystyle |{d}^{2}S|}のように、面素に絶対値記号をつけることにする。

右辺は、(u1,u2){displaystyle ({u}_{1},{u}_{2})}{displaystyle ({u}_{1},{u}_{2})}についてのベクトル値関数
X(φ(u1,u2)) |∂φ∂u1×∂φ∂u2|{displaystyle mathbf {X} (varphi ({u}_{1},{u}_{2})) left|{frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right|}{displaystyle mathbf {X} (varphi ({u}_{1},{u}_{2})) left|{frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right|}


を、成分ごとに区間I上で重積分したものを意味する。即ち、絶対値による面積分は「関数の面積分を各成分ごとにやる」というのと同じ意味である。



(3)外積面積分


定義は、以下の通りである。
∫s∈SX(s)×d2S:=∫(u1,u2)∈I(X(φ(u1,u2)))×(∂φ∂u1×∂φ∂u2) du1du2{displaystyle {int }_{mathbf {s} in S}mathbf {X} (mathbf {s} )times {d}^{2}S:={int }_{({u}_{1},{u}_{2})in I}left(mathbf {X} (varphi ({u}_{1},{u}_{2}))right)times left({frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right) d{u}_{1}d{u}_{2}}{displaystyle {int }_{mathbf {s} in S}mathbf {X} (mathbf {s} )times {d}^{2}S:={int }_{({u}_{1},{u}_{2})in I}left(mathbf {X} (varphi ({u}_{1},{u}_{2}))right)times left({frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{2}}}right) d{u}_{1}d{u}_{2}}



本記事では、外積面積分の場合は、× d2S{displaystyle times {d}^{2}S}{displaystyle times {d}^{2}S}のように、面素にの前に、×をつけることにする。

右辺は、(u1,u2){displaystyle ({u}_{1},{u}_{2})}{displaystyle ({u}_{1},{u}_{2})}についてのベクトル値関数
X(φ(u1,u2))×(∂φ∂u1×∂φ∂u1){displaystyle mathbf {X} (varphi ({u}_{1},{u}_{2}))times left({frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{1}}}right)}{displaystyle mathbf {X} (varphi ({u}_{1},{u}_{2}))times left({frac {partial varphi }{partial {u}_{1}}}times {frac {partial varphi }{partial {u}_{1}}}right)}



を、成分ごとに区間I上で重積分したものを意味する。即ち、(2)との関係でいえば、
∫s∈SX(s)×d2S=∫s∈SX(s)×nS|d2S|{displaystyle {int }_{mathbf {s} in S}mathbf {X} (mathbf {s} )times {d}^{2}S={int }_{mathbf {s} in S}mathbf {X} (mathbf {s} )times mathbf {n} _{S}|{d}^{2}S|}{displaystyle {int }_{mathbf {s} in S}mathbf {X} (mathbf {s} )times {d}^{2}S={int }_{mathbf {s} in S}mathbf {X} (mathbf {s} )times mathbf {n} _{S}|{d}^{2}S|}


となる。



脚注


[ヘルプ]


注釈




  1. ^
    尚、磁化がΩ{displaystyle Omega }Omega の外では0であるからといって、Ω{displaystyle Omega }Omega の外でAM{displaystyle mathbf {A} _{M}}{displaystyle mathbf {A} _{M}}が0となるとは限らない（大概の場合はΩ{displaystyle Omega }Omega の外でもAM{displaystyle mathbf {A} _{M}}{displaystyle mathbf {A} _{M}}は、0ではない）ことに注意されたい。


  2. ^
    なお、磁化は、Ω{displaystyle Omega }Omega の外では0なので、上式は、
    AM(r)=μ04π∫s∈R3((M(s))×(r−s)|r−s|) d3s{displaystyle mathbf {A} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in mathbb {R} ^{3}}left({frac {(mathbf {M} (mathbf {s} ))times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right) {d}^{3}mathbf {s} }{displaystyle mathbf {A} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in mathbb {R} ^{3}}left({frac {(mathbf {M} (mathbf {s} ))times (mathbf {r} -mathbf {s} )}{|mathbf {r} -mathbf {s} |}}right) {d}^{3}mathbf {s} }

    としてもよい（紛らわしいのであまり推奨しない）。


  3. ^
    尚、磁化がΩ{displaystyle Omega }Omega の外では0であるからといって、Ω{displaystyle Omega }Omega の外でBM{displaystyle mathbf {B} _{M}}{displaystyle mathbf {B} _{M}}が0となるとは限らない（大概の場合はΩ{displaystyle Omega }Omega の外でもBM{displaystyle mathbf {B} _{M}}{displaystyle mathbf {B} _{M}}は、0ではない）ことに注意されたい。


  4. ^
    なお、磁化は、Ω{displaystyle Omega }Omega の外では0なので、上式は、
    BM(r)=μ04πgradr⁡[<M(s)|(r−s)>|r−s|3]+divr⁡[r−s|r−s|]M(s) d3s{displaystyle mathbf {B} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]+operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {M} (mathbf {s} ) {d}^{3}mathbf {s} }{displaystyle mathbf {B} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right]+operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {M} (mathbf {s} ) {d}^{3}mathbf {s} }

    としてもよい（紛らわしいのであまり推奨しない）。


  5. ^
    なお、磁化は、Ω{displaystyle Omega }Omega の外では0なので、上式は、
    M(s) = ∫s∈R3 divr⁡[r−s|r−s|]M(s) d3s{displaystyle mathbf {M} (mathbf {s} ) = {int }_{mathbf {s} in mathbf {R} ^{3}} operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {M} (mathbf {s} ) {d}^{3}mathbf {s} }{displaystyle mathbf {M} (mathbf {s} ) = {int }_{mathbf {s} in mathbf {R} ^{3}} operatorname {div} _{mathbf {r} }left[{frac {mathbf {r} -mathbf {s} }{|mathbf {r} -mathbf {s} |}}right]mathbf {M} (mathbf {s} ) {d}^{3}mathbf {s} }

    でもある。


  6. ^
    なお、磁化は、Ω{displaystyle Omega }Omega の外では0なので、上式は、
    BM(r)=μ04π∫s∈R3gradr⁡[<M(s)|(r−s)>|r−s|3] d3s+μ0M(r){displaystyle mathbf {B} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in mathbf {R} ^{3}}operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right] {d}^{3}{s}+mu _{0}mathbf {M} (mathbf {r} )}{displaystyle mathbf {B} _{M}(mathbf {r} )={frac {mu _{0}}{4pi }}{int }_{mathbf {s} in mathbf {R} ^{3}}operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right] {d}^{3}{s}+mu _{0}mathbf {M} (mathbf {r} )}

    としてもよい（紛らわしいのであまり推奨しない）。


  7. ^
    なお、磁化は、Ω{displaystyle Omega }Omega の外では0なので、上式は、
    HM(r)=14π∫s∈Ωgradr⁡[<M(s)|(r−s)>|r−s|3] d3s{displaystyle mathbf {H} _{M}(mathbf {r} )={frac {1}{4pi }}{int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right] {d}^{3}{s}}{displaystyle mathbf {H} _{M}(mathbf {r} )={frac {1}{4pi }}{int }_{mathbf {s} in Omega }operatorname {grad} _{mathbf {r} }left[{frac {<mathbf {M} (mathbf {s} )|(mathbf {r} -mathbf {s} )>}{{|mathbf {r} -mathbf {s} |}^{3}}}right] {d}^{3}{s}}

    としてもよい（紛らわしいのであまり推奨しない）。

  8. ^ abここでいうδは、摂動あるいは変分を表し、ディラックのδと紛らわしいが全然別物である。



参考文献









  1. ^
    W.K.H. パノフスキー (著), M. フィリップス (著),林 忠四郎 (翻訳), 西田 稔 (翻訳);「新版 電磁気学〈上〉」吉岡書店; POD版 (2002/09)

  2. ^ abcd
    [★溝口 正 (著) 「電磁気学―SI UNITS 」裳華房 (2001/03)](特にP188付近)


  3. ^
    [★]竹山 説三 (著) 「電磁気学現象理論」丸善出版; 3版 (1949)

  4. ^ ab[★]P.P.シルベスタ(著),
    R.L.フェラーリ（著）,
    本間 利久（著）,
    田中 康博（著）;
    「有限要素法による電磁界解析 (Information & computing (26)) 」
    サイエンス社 (1988/09)

    [原書]Peter P. Silvester (Author), Ronald L. Ferrari (Author);
    「Finite Elements for Electrical Engineers」Cambridge University Press; 3版 (1996/9/5)
    [1]


  5. ^ 斎藤兆古,坂本禎智,藤原耕二；「磁気回路法と有限要素法の理論的関係」 (PDF) 電気学会マグネティックス研究会資料 MAG-03 号：54-64 ページ：5-10 発行年：2003年03月31日


  6. ^ 静磁場解析のための二次要素を用いる有限要素法の研究


  7. ^ 有限要素法による磁場解析

  8. ^ ab磁場が満たす偏微分方程式（秋田高専講義録）


  9. ^ 園田英徳;「大学院生のための基礎物理学」講談社（2011/09/29）ISBN 978-4-06-153277-9


  10. ^ 平川浩正;「電磁気学（新物理学シリーズ２）」培風館 (1986/04) ISBN 9784563024024


  11. ^ 守末 利弥「数値電磁気学のためのゲージ理論」森北出版 (1996/04) ISBN 978-4-627-71600-1


  12. ^ 依田 潔 (著) ;「Mathematicaによる電磁界シミュレーション入門 - ＰＯＤ版 (計算電気・電子工学シリーズ)」森北出版(2012/2/24) ISBN 978-4-627-71529-5


  13. ^
    電磁気学II(大阪大学 田中実教授の講義録） (PDF)


  14. ^ [2]立教大学講義ノート

    本講義の参考文献欄
    等、至る所に文字化けがあるが、たとえば、もじばけらった等によって解読できる。


  15. ^
    東京理科大学講義ノート (PDF)
    および、講義スライド (PDF)


  16. ^ 鹿児島 誠一 (著) ;「電磁気学 (パリティ物理学コース)」丸善(1997/01) ISBN 978-4-621-04277-9


  17. ^ 後藤 憲一(著), 山崎修一郎(著) ;「詳解 電磁気学演習」共立出版 (1970/12) ISBN 978-4-320-03022-0


  18. ^ 前田 和茂 (著) ,小林 俊雄(著);「ビジュアルアプローチ電磁気学」森北出版 (2009/12/5) ISBN 978-4-627-16221-1


  19. ^ Julius Adams Stratton;「Electromagnetic Theory」Wiley-IEEE Press; 1版 (2007/1/22) ISBN 978-0-470-13153-4


  20. ^ A. Pramanik (著) ;「Electro-Magnetism: Theory and Applications 」Prentice-Hall of India Pvt.Ltd (2004/8/15)
    [3]


  21. ^ Bo-nan Jiang(著) ;「The Least-Squares Finite Element Method: Theory and Applications in Computational Fluid Dynamics and Electromagnetics (Scientific Computation) 」Springer; 1998版 (1998/6/22) [4]


  22. ^ 電気磁気学特論（秋田高専講義録）


  23. ^ 高橋 康人 (早稲田大学):"高速多重極法を組み込んだ磁気モーメント法による 磁性体解析に関する基礎的検討" (PDF)


  24. ^ 東京電機大学 (編集)「入門 電磁気学」東京電機大学出版局 (2006/03) ISBN 9784501004200


  25. ^ 早川 義晴 (著) 「電気教科書 電験三種合格ガイド」翔泳社 (2011/2/25) ISBN 9784798126623


  26. ^
    磁気回路と電気回路 (PDF)


  27. ^ 坪井 一洋 (著) 「システムと微分方程式」三恵社 (2011/5/22) ISBN 9784883618248


  28. ^ 複数磁石による静磁場(簡易シミュレータ)


  29. ^ 近畿大学 講義ノート (PDF)


  30. ^
    等々力 二郎：「有限の太さの矩形断面ヘリカルコイルの磁場の計算」核融合研究,57(1987)318


  31. ^ 渡辺 二太:「空心円筒コイル群によって作られる磁場の計算」 核融合研究 Vol.35 (1976) No.3 P235-242


  32. ^ 渡辺 二太:「多様な形状のコイルに対する磁場計算法」 核融合研究 Vol.63 (1990) No.6 P482-507


  33. ^
    藤本 淳夫 (著) ベクトル解析 (現代数学レクチャーズ C- 1) 培風館 (1979/01)




関連項目



  • 磁場

  • 静電気学

  • ベクトルポテンシャル

  • 遅延ポテンシャル

  • ジェフィメンコ方程式

  • アンペールの法則

  • ビオ・サバールの法則

  • クーロンの法則


  • 電場の強さ(E)、磁束密度(B)、電束密度(D)

  • E-B対応とE-H対応

  • 地磁気

  • ホール素子

  • SQUID

  • ファラデー効果

  • マクスウェルの方程式

  • ローレンツ力

  • フレミングの法則



外部リンク




  • 磁化Mと磁場Hについて (PDF)


  • 磁場 (PDF)


  • 電気磁気学演習解答例 (PDF)




-

Popular posts from this blog

Coverage of Google Street View

Image
Coverage of Google Street View From Wikipedia, the free encyclopedia Jump to navigation Jump to search Countries and dependencies with:    mostly full coverage    partial coverage    full or partial coverage planned (official)    full or partial coverage planned (unofficial)    views of selected businesses and/or tourist attractions only    views of private businesses only    no current or planned coverage Google Street View was first introduced in the United States on May 25, 2007, and until November 26, 2008, featured camera icon markers, each representing at least one major city or area (such as a park), and usually the other nearby cities, towns, suburbs, and parks. Many areas that had coverage were represented by icons. Contents 1 Key additions 2 Timeline of introductions 2.1 2007 2.2 2008 2.3 2009 2.4 2010 2.5 2011 2.6 2012 2.7 2013 2.8 2014 2.9 20...
Read more

Full-time equivalent

Image
Full-time equivalent From Wikipedia, the free encyclopedia Jump to navigation Jump to search Full-time equivalent ( FTE ) or whole time equivalent ( WTE ) is a unit that indicates the workload of an employed person (or student) in a way that makes workloads or class loads comparable [1] across various contexts. FTE is often used to measure a worker's or student's involvement in a project, or to track cost reductions in an organization. An FTE of 1.0 is equivalent to a full-time worker or student, while an FTE of 0.5 signals half of a full work or school load. [2] Contents 1 U.S. Federal Government 2 In education 2.1 Example 3 Notes 4 References U.S. Federal Government [ edit ] In the U.S. Federal Government, FTE is defined by the Government Accountability Office (GAO) as the number of total hours worked divided by the maximum number of compensable hours in a full-time schedule as...
Read more

Surfing

Image
Surfing From Wikipedia, the free encyclopedia Jump to navigation Jump to search This article is about stand-up ocean surfing. For other uses, see Surfing (disambiguation). "Surfer" redirects here. For other uses, see Surfer (disambiguation). Surfing A surfer at the Cayucos Pier, Cayucos, California Highest governing body World Surf League (WSL) Characteristics Mixed gender Yes, separate competitions Presence Country or region Worldwide Olympic Will debut in 2020 Surfing on the Gold Coast, Queensland, Australia Surfing is a surface water sport in which the wave rider, referred to as a surfer, rides on the forward or deep face of a moving wave, which usually carries the surfer towards the shore. Waves suitable for surfing are primarily found in the ocean, but can also be found in lakes or rivers in the form of a standing wave or tidal bore. However, surfers can also utilize artifi...
Read more