静磁場

提供: testwiki
ナビゲーションに移動 検索に移動

静磁場(せいじば、テンプレート:Lang-en)とは、時間的に変動しない磁場のことである [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)の視点から、静磁場について述べる。

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

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


真空中に定常な(つまり時刻tに依存しない)電流密度 𝒊(𝐫)が与えられた とする。このとき、𝒊(𝐫) は、以下の磁気ベクトルポテンシャル 𝐀𝒊(𝐫)を空間内に作り出す。

𝐀𝒊(𝐫)=μ04π𝐬3(𝒊(𝐬)|𝐫𝐬|) d3𝐬 (1-1)

となる。

𝐁𝒊=rot[𝐀𝒊]を考え併せると、 𝒊(𝐫)が直接的に作り出す磁束密度𝐁𝒊は、

𝐁𝒊(𝐫)=μ04π𝐬3(𝒊(𝐬)×(𝐫𝐬)|𝐫𝐬|3)d3𝐬 (1-2)

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

上記の𝐁𝒊(𝐫)に対し、 新たな場𝐇𝒊を、

𝐇𝒊(𝐫):=1μ0𝐁𝒊(𝐫) (1-3)

と定義する。この場のことを、「電流密度iが作り出す磁場」と呼ぶ。ここでμ0は、真空の透磁率である。 尚、定義の上では、「電流密度iが作り出す磁場」𝐇𝒊は、 透磁率がμの場所でも、:𝐇𝒊(𝐫):=1μ0𝐁𝒊(𝐫)であることに特に注意されたい。

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

𝐇𝒊(𝐫)=14π𝐬3(𝒊(𝐬)×(𝐫𝐬)|𝐫𝐬|3)d3𝐬 (1-4)

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

rot𝐫[𝐇𝒊]=𝒊 (1-5)

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

rot𝐫[(𝒊(𝐬)×(𝐫𝐬)|𝐫𝐬|3)]=4π(div𝐫[(𝐫𝐬)|𝐫𝐬|3])𝒊(𝐬)=4πδ3(𝐫𝐬)𝒊(𝐬)(1-6)

従って、

rot𝐫[𝐇𝒊(𝐫)]=14π𝐬3rot𝐫[(𝒊(𝐬)×(𝐫𝐬)|𝐫𝐬|3)]d3𝐬=𝐬3δ3(𝐫𝐬)𝒊(𝐬)d3𝐬=𝒊(𝐫)(1-7)

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

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

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


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

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

物質内の各点𝐬Ωそれぞれに、それぞれ
𝐌(𝐬)d3𝐬 (2-1-1)
で与えられる磁気モーメントが配置されている

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

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

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

𝐚0(𝐫)=μ04π(𝐦×𝐫|𝐫|) (2-1-2)
である。

(2)従って、"(1)"を平行移動すれば、位置𝐬3に置かれた磁気モーメントmは、空間上の位置rに作り出す磁気ベクトルポテンシャルは、

𝐚s(𝐫)=μ04π(𝐦×(𝐫𝐬)|𝐫𝐬|) (2-1-3)
である。

従って、上記の磁気モーメント𝐌(𝐬)d𝐬それぞれは、磁気ベクトルポテンシャル

μ04π((𝐌(𝐬)d𝐬)×(𝐫𝐬)|𝐫𝐬|) (2-1-4)

を作り出す。

上記の磁気ベクトルポテンシャルそれぞれを、全ての𝐬Ωに渡って足し合わせると、

𝐀M(𝐫)=μ04π𝐬Ω((𝐌(𝐬))×(𝐫𝐬)|𝐫𝐬|) d3𝐬 (2-1-5)

を得る。即ち、物質の磁化𝐌(𝐫) は、上記の磁気ベクトルポテンシャル 𝐀Mを空間内に作り出す[2] [注釈 1] [注釈 2]

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

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

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

𝐚0(𝐫)=μ04π(𝐦×𝐫|𝐫|) (2-2-1)
である。これの回転微分をとることで、原点に置かれた磁気モーメントmが、空間上の位置rに作り出す磁束密度は、
𝐛0(𝐫)=rot𝐫[μ04π(𝐦×𝐫|𝐫|)]=μ04πgrad𝐫[𝐦|𝐫|𝐫|3]+μ04πdiv𝐫[𝐫|𝐫|]𝐦 (2-2-2)
であることが判る。ここで、< | >は、内積を表す。

(2)従って、"(1)"を平行移動すれば、位置𝐬3に置かれた磁気モーメントmが、空間上の位置r に作り出す磁束密度が、

𝐛s(𝐫)=rot𝐫[μ04π(𝐦×(𝐫𝐬)|𝐫𝐬|)]=μ04π grad𝐫[𝐦|(𝐫𝐬)|𝐫𝐬|3]+μ04π div𝐫[𝐫𝐬|𝐫𝐬|]𝐦 (2-2-3)
であることが判る。

従って、上記の磁気モーメント𝐌(𝐬)d𝐬それぞれは、磁束密度

μ04πgrad𝐫[𝐌(𝐬)d3𝐬|(𝐫𝐬)|𝐫𝐬|3]+μ04πdiv𝐫[𝐫𝐬|𝐫𝐬|]𝐌(𝐬)d3𝐬 (2-2-4)

を作り出す。

上記の磁気ベクトルポテンシャルそれぞれを、全ての𝐬Ωに渡って足し合わせると、

𝐁M(𝐫)=μ04π𝐬Ωgrad𝐫[𝐌(𝐬)|(𝐫𝐬)|𝐫𝐬|3]+div𝐫[𝐫𝐬|𝐫𝐬|]𝐌(𝐬) d3𝐬 (2-2-5)

を得る。即ち、物質の磁化𝐌(𝐫) は、上記の磁気ベクトルポテンシャル 𝐁Mを空間内に作り出す [2] [注釈 3] [注釈 4]

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

𝐌(𝐬) = 𝐬Ω div𝐫[𝐫𝐬|𝐫𝐬|]𝐌(𝐬) d3s (2-2-6)

が判る。[注釈 5]

従って、上記の磁気モーメント𝐌(𝐬)d𝐬それぞれが作り出す磁束密度は、

μ04πgrad𝐫[𝐌(𝐬)d3𝐬|(𝐫𝐬)|𝐫𝐬|3]+μ0δ(𝐫𝐬)𝐌(𝐬)d3𝐬 (2-2-7)

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

𝐁M(𝐫)=μ04π𝐬Ωgrad𝐫[𝐌(𝐬)|(𝐫𝐬)|𝐫𝐬|3] d3s+μ0𝐌(𝐫) (2-2-8)

を得る。[注釈 6]

今、新たな場𝐇Mを、

𝐇M(𝐫):=14π𝐬Ωgrad𝐫[𝐌(𝐬)|(𝐫𝐬)|𝐫𝐬|3] d3s (2-2-9)

と定義すると [注釈 7]

𝐁M(𝐫)=μ0(𝐇M(𝐫)+𝐌(𝐫)) (2-2-10)

が得られる。

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

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

(難所)元々磁化Miniを帯びている物質に、外部から磁束密度𝐁ext(𝐫)を印加したとき、最終的に、物質がどのようなMconを得るか?

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

磁荷密度の導入

空間内の領域Ωに物質が置かれ、前記物質が、定常な(つまり時刻tに依存しない)磁化𝐌(𝐫)を帯びているとする。

先に定義した、(磁化による)磁場

𝐇M(𝐫)=14π𝐬Ωgrad𝐫[𝐌(𝐬)|(𝐫𝐬)|𝐫𝐬|3] d3s (2-3-1)

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

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

divs[𝐌(𝐬)|𝐫𝐬|]=grad𝐬[1|𝐫𝐬|]|𝐌(𝐬)+div𝐬[𝐌]|𝐫𝐬|=𝐫𝐬|𝐫𝐬|3|𝐌(𝐬)+div𝐬[𝐌]|𝐫𝐬|=𝐌(𝐬)|(𝐫𝐬)|𝐫𝐬|3+div𝐬[𝐌]|𝐫𝐬|(2-3-2)

従って、

𝐌(𝐬)|(𝐫𝐬)|𝐫𝐬|3=divs[𝐌(𝐬)|𝐫𝐬|]div𝐬[𝐌]|𝐫𝐬| (2-3-3)

今、スカラー値関数ϕMを、

ϕM(𝐫):=𝐬Ω𝐌(𝐬)|(𝐫𝐬)|𝐫𝐬|3 d3Ω (2-3-4)

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

ϕM(𝐫)=𝐬Ωdivs[𝐌(𝐬)|𝐫𝐬|] d3Ω+𝐬Ωdiv𝐬[𝐌]|𝐫𝐬| d3Ω (2-3-5)

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

𝐬Ωdivs[𝐌(𝐬)|𝐫𝐬|] d3Ω=𝐬Ω𝐌(𝐬)|𝐫𝐬| d2(Ω) (2-3-6)

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

𝐬Ω𝐌(𝐬)|𝐫𝐬| d2(Ω)=𝐬Ω𝐌(𝐬)|𝐧Ω|𝐫𝐬| |d2(Ω)| (2-3-7)

従って、

ϕM(𝐫)=𝐬Ωdivs[𝐌(𝐬)|𝐫𝐬|] d3Ω +𝐬Ω𝐌(𝐬)|𝐧Ω|𝐫𝐬| |d2(Ω)| (2-3-8)

ここで、Ωは、領域Ωの境界面を意味する。 また、𝐧Ωは、Ωの法線ベクトルを意味する。

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

  • ρM:=div𝐫[𝐌] (体積磁荷密度) (2-3-10)
  • σM:=<𝐌 | 𝐧Ω> (表面磁荷密度) (2-3-11)

により定めると、

ϕM(𝐫)=𝐬ΩρM(𝐬)|𝐫𝐬| d3Ω +𝐬ΩσM(𝐬)|𝐫𝐬| |d2(Ω)| (2-3-12)

を得る。

一方、ϕM(𝐫)の定義により、

𝐇M(𝐫)=grad𝐫[ϕM(𝐫)] (2-3-13)

であるため、

𝐇M(𝐫)=𝐬Ωgrad𝐫[ρM(𝐬)|𝐫𝐬|] d3Ω +𝐬Ωgrad𝐫[σM(𝐬)|𝐫𝐬|] |d2(Ω)| (2-3-14)

である。

磁化電流の導入

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

前記𝐀M(𝐫)に、 ベクトル解析の公式を適用すると、

𝐀M(𝐫)=μ04π𝐬Ω((𝐌(𝐬))×(𝐫𝐬)|𝐫𝐬|) d3𝐬=μ04π𝐬Ω(rot𝐬[𝐌(𝐬)]|𝐫𝐬|) d3𝐬+μ04π𝐬Ω(𝐌(𝐬)×𝐧Ω|𝐫𝐬|) |d2(Ω)|(2-4-1)

となることが判る。ここで、rot𝐬 は、変数𝐬についての回転微分を意味する。Ωは、領域Ωの境界を意味する。 上式の右辺第二項の積分において、「Ωの面素の絶対値」を意味し 所謂普通の面積分ではない(この「面積分」は、「ベクトル場からベクトルを作る」特殊なもの)ので注意が必要である。 この積分において、𝐧Ωは、Ω の法線ベクトルを意味する。

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

rot𝐬[𝐌(𝐬)|𝐫𝐬|]=(grad𝐬[𝐌(𝐬)|𝐫𝐬|])×𝐌(𝐬)+rot𝐬[𝐌(𝐬)]|𝐫𝐬|=(𝐫𝐬)×𝐌(𝐬)|𝐫𝐬|3+rot𝐬[𝐌(𝐬)]|𝐫𝐬|(2-4-2)

従って、

𝐌(𝐬)×(𝐫𝐬)|𝐫𝐬|3=rot𝐬[𝐌(𝐬)]|𝐫𝐬|rot𝐬[𝐌(𝐬)|𝐫𝐬|] (2-4-3)

従って、

𝐀M(𝐫)=μ04π𝐬Ωrot𝐬[𝐌(𝐬)]|𝐫𝐬| d3𝐬 μ04π𝐬Ωrot𝐬[𝐌(𝐬)|𝐫𝐬|] d3𝐬   (2-4-4)

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

μ04π𝐬Ωrot𝐬[𝐌(𝐬)|𝐫𝐬|] d3𝐬 =μ04π𝐬Ω(𝐌(𝐬)|𝐫𝐬|)×d2𝐬   (2-4-5)

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

μ04π𝐬Ω(𝐌(𝐬)|𝐫𝐬|)×d2𝐬 =μ04π𝐬Ω(𝐌(𝐬)×𝐧Ω|𝐫𝐬|) |d2(Ω)|  (2-4-6)

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

今、𝒊M(𝐫), 𝐊Mを、

  • 𝒊M(𝐫):=rot[𝐌(𝐫)] (体積磁化電流密度) (2-4-7)
  • 𝐊M(𝐫):=𝐌(𝐫)×𝐧Ω(𝐫) (表面磁化電流密度) (2-4-8)

と置くと、結局、

𝐀M(𝐫)=μ04π𝐬Ω(𝒊M(𝐬)|𝐫𝐬|) d3𝐬+ μ04π𝐬Ω(𝐊M(𝐬)|𝐫𝐬|) |d2(Ω)|  (2-4-9)

となる。

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

𝐁M(𝐫)=μ04π𝐬3(𝒊M(𝐬)×(𝐫𝐬)|𝐫𝐬|3)d3𝐬 + μ04π𝐬Ω(𝐊M(𝐬)×(𝐫𝐬)|𝐫𝐬|3) |d2(Ω)| (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を帯びている物質に、外部から磁束密度𝐁ext(𝐫)を印加したとき、最終的に、物質がどのようなMconを得るか?

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

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

div[𝐁fc]=0 (3-1-2a)
div[𝐁M]=0 (3-1-2b)
div[𝐁tot]=0 (3-1-2c)

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

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

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

𝐇tot(𝐫):=𝐁tot(𝐫)μ0𝐌(𝐫) (3-2-1)

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


さらに、強制電流𝒊fcが作る磁束密度𝐁fc(𝐫)は、式(1-2)に より、(周辺の物質の有無にかかわらず)定まるが、式(1-3)に倣い 「強制電流𝒊fcが作る磁場𝐇fc(𝐫)」を、

𝐇fc(𝐫):=1μ0𝐁fc(𝐫) (3-2-2)

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

𝐇tot(𝐫)=𝐁tot(𝐫)μ0𝐌(𝐫)=𝐁fc(𝐫)+𝐁M(𝐫)μ0𝐌(𝐫)=𝐇fc(𝐫)+(𝐇M(𝐫)+𝐌(𝐫))𝐌(𝐫)=𝐇fc(𝐫)+𝐇M(𝐫)(3-2-3)

が得られる。従って、

rot[𝐇tot(𝐫)]=rot[𝐇fc(𝐫)]+rot[𝐇M(𝐫)]=𝒊fc (3-2-4)

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

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

rot grad=0 (3-2-5)

𝐇M(𝐫)に対し、回転微分を作用させると、

rot[𝐇M(𝐫)]=14π𝐬Ωrot𝐫[grad𝐫[𝐌(𝐬)|(𝐫𝐬)|𝐫𝐬|3]] d3s =0 (3-2-6)

が判る。従って、磁化に起因する磁場𝐇M(𝐫)は、

rot[𝐇M(𝐫)]=0 (3-2-7)

を満たす。


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

透磁率の導入

透磁率なる概念を、新たに導入する。 即ち、3次正方行列に値を取る行列値関数μ(𝐫)(透磁率)を用いて、

𝐁tot(𝐫)=μ(𝐫)𝐇tot(𝐫) (3-3-1)

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

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

𝐁tot(𝐫)=μ(𝐫)𝐇tot(𝐫) (3-3-2)

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

𝐇tot(𝐫)=ν(𝐫)𝐁tot(𝐫) (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[𝐁fc]=0 (3-3-6a)
div[𝐁M]=0 (3-3-6b)
div[𝐁tot]=0 (3-3-6c)
rot[𝐇tot]=𝒊fc (3-3-7)

を出発点とし、全系の磁気ベクトルポテンシャル𝐀totが充たす微分方程式を導出する。

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

rot[𝐀fc]=𝐁fc (3-3-8a)
rot[𝐀M]=𝐁M (3-3-8b)

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

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

𝐀tot:=𝐀M+𝐀fc (3-3-9)

と定めると、

rot[𝐀tot]=𝐁tot (3-3-8c)

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

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

rot[νrot[𝐀]]𝒊fc=0 (3-3-10)

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

静磁場の境界条件

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


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

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

F[𝐀]:=𝐬312< 𝐇(𝐬) | 𝐁(𝐬)>< 𝒊fc(𝐬)| 𝐀tot(𝐬)> d𝐬=𝐬312< ν(𝐬)(rot𝐬[𝐀(𝐬)]) | (rot𝐬[𝐀(𝐬)])>< 𝒊fc(𝐬)| 𝐀(𝐬)> d𝐬(3-5-1)

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


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

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

を求める。

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

12< ν(𝐬)(rot𝐬[𝐀(𝐬)+δ𝐀(𝐬)]) | (rot𝐬[𝐀(𝐬)+δ𝐀(𝐬)])> < 𝒊fc(𝐬)| 𝐀(𝐬)+δ𝐀(𝐬)> (3-5-3)

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

12< ν(𝐬)(rot𝐬[𝐀(𝐬)]) | (rot𝐬[𝐀(𝐬)])> < 𝒊fc(𝐬)| 𝐀(𝐬)> (3-5-4)

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

< ν(𝐬)(rot𝐬[𝐀(𝐬)]) | (rot𝐬[δ𝐀(𝐬)])>+12< ν(𝐬)(rot𝐬[δ𝐀(𝐬)]) | (rot𝐬[δ𝐀(𝐬)])> < 𝒊fc(𝐬)| δ𝐀(𝐬)> (3-5-5)

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

δF[A]=𝐬3< ν(𝐬)(rot𝐬[𝐀(𝐬)]) | (rot𝐬[δ𝐀(𝐬)])> < 𝒊fc(𝐬)| δ𝐀(𝐬)> d3𝐬 (3-5-6)

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

div[𝐗×𝐘]=< 𝐘 | rot[𝐗] >< 𝐗 | rot[𝐘] > (3-5-7)

において

𝐗:=ν(𝐬)rot[𝐀(𝐬)](3-5-8)𝐘:=δ𝐀(3-5-9)

を代入すると、

δF[𝐀]=𝐬3div𝐬[(ν(𝐬)rot𝐬[𝐀])×δ𝐀] d3𝐬+𝐬3<(rot𝐬[(ν(𝐬))rot𝐬[𝐀(𝐬)] ]𝒊fc(𝐬) | δ𝐀(𝐬)> d3𝐬(3-5-10)

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

𝐬Ωdiv𝐬[(ν(𝐬)rot𝐬[𝐀])×δ𝐀 ] d3𝐬=𝐬Ω(ν(𝐬)rot𝐬[𝐀])×δ𝐀 d2(Ω) (3-5-11)

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

(rot[𝐀])×nΩ = 0 (3-5-12)

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

rot[νrot[𝐀]]𝒊fc=0 (3-5-12)

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

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

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

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

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

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

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

 𝐅 |  :=f1x1+f2x2+f3x3

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

:=𝐱^x+𝐲^y+𝐳^z

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

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

 𝐅 |  𝐆=f1g1x1+f2g2x2+f3g3x3=(J[𝐆])𝐅

が成り立つ。ここで、J[𝐆]は、Gのヤコビ行列(ヤコビアンではない)を意味する。

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

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

rot[f𝐅]=grad[f]×𝐅+frot[𝐅]
rot[𝐅×𝐆]= 𝐆 |  𝐅 𝐅 |  𝐆+𝐅div[𝐆]𝐆div[𝐅]

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

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

div[f𝐅]=𝐅|grad[f]+fdiv[𝐅]

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

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

f(𝐫)=𝐬Ωf(𝐬)δn(𝐫𝐬) dns =f*δn

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

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

原点(r=0)を除いて

div𝐫[𝐫|𝐫|]=3|𝐫|33|𝐫|2|𝐫|5=0
であり、原点(r=0)を中心とする、球体BLに対し、ガウスの発散定理を用いると、
𝐫BLdiv𝐫[𝐫|𝐫|] d3r=4π
となるので、結局、
div𝐫[𝐫|𝐫|]=4πδ(𝐫)

であることが判る。(ここで、δは、3変数のδ関数である。
以上の議論を平行移動させると、

div𝐫[𝐫𝐬|𝐫𝐬|]=4πδ(𝐫𝐬)

であることが判る。

分数関数の偏微分

𝐬i[1|𝐫𝐬|]=si[1((risi)2+const.)1/2]=si[((risi)2+const.)(1/2)]

面積分に関する補足

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

Iを2の閉長方形、
φ:I3は、 Iの近傍で、殆ど至る所区分的に滑らかかつ、非退化であり、かつ、Iの内部で単射なベクトル値関数、
S:=φ[I]を、3内の曲面片とする。
また、Xを、Sの近傍で定義された、区分的に滑らかなベクトル場とする。このとき、

(1)普通の面積分

定義は、以下の通りである。
𝐬S𝐗(𝐬) d2S:=(u1,u2)I𝐗(φ(u1,u2)) | (φu1×φu2) du1du2
上式の右辺は、(u1,u2)についてのスカラー値関数
𝐗(φ(u1,u2)) | (φu1×φu2)
を、区間I上で重積分したものを意味する。今、Sの単位法線ベクトル𝐧Sを、
𝐧S:=1|φu1×φu2|(φu1×φu2)
と定めると、
𝐬S𝐗(𝐬) d2S=(u1,u2)I𝐗(φ(u1,u2)) | 𝐧S|φu1×φu2| du1du2
である。

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

定義は、以下の通りである。
𝐬S𝐗(𝐬) |d2S|:=(u1,u2)I𝐗(φ(u1,u2)) |φu1×φu2| du1du2
本記事では、絶対値による面積分の場合は、 |d2S|のように、面素に絶対値記号をつけることにする。
右辺は、(u1,u2)についてのベクトル値関数
𝐗(φ(u1,u2)) |φu1×φu2|
を、成分ごとに区間I上で重積分したものを意味する。即ち、絶対値による面積分は「関数の面積分を各成分ごとにやる」というのと同じ意味である。

(3)外積面積分

定義は、以下の通りである。
𝐬S𝐗(𝐬)×d2S:=(u1,u2)I(𝐗(φ(u1,u2)))×(φu1×φu2) du1du2
本記事では、外積面積分の場合は、× d2Sのように、面素にの前に、×をつけることにする。
右辺は、(u1,u2)についてのベクトル値関数
𝐗(φ(u1,u2))×(φu1×φu1)
を、成分ごとに区間I上で重積分したものを意味する。即ち、(2)との関係でいえば、
𝐬S𝐗(𝐬)×d2S=𝐬S𝐗(𝐬)×𝐧S|d2S|

となる。

脚注

テンプレート:脚注ヘルプ

注釈

  1. 尚、磁化がΩの外では0であるからといって、Ωの外で𝐀Mが0となるとは限らない(大概の場合はΩの外でも𝐀Mは、0ではない)ことに注意されたい。
  2. なお、磁化は、Ωの外では0なので、上式は、
    𝐀M(𝐫)=μ04π𝐬3((𝐌(𝐬))×(𝐫𝐬)|𝐫𝐬|) d3𝐬
    としてもよい(紛らわしいのであまり推奨しない)。
  3. 尚、磁化がΩの外では0であるからといって、Ωの外で𝐁Mが0となるとは限らない(大概の場合はΩの外でも𝐁Mは、0ではない)ことに注意されたい。
  4. なお、磁化は、Ωの外では0なので、上式は、
    𝐁M(𝐫)=μ04πgrad𝐫[𝐌(𝐬)|(𝐫𝐬)|𝐫𝐬|3]+div𝐫[𝐫𝐬|𝐫𝐬|]𝐌(𝐬) d3𝐬
    としてもよい(紛らわしいのであまり推奨しない)。
  5. なお、磁化は、Ωの外では0なので、上式は、
    𝐌(𝐬) = 𝐬𝐑3 div𝐫[𝐫𝐬|𝐫𝐬|]𝐌(𝐬) d3𝐬
    でもある。
  6. なお、磁化は、Ωの外では0なので、上式は、
    𝐁M(𝐫)=μ04π𝐬𝐑3grad𝐫[𝐌(𝐬)|(𝐫𝐬)|𝐫𝐬|3] d3s+μ0𝐌(𝐫)
    としてもよい(紛らわしいのであまり推奨しない)。
  7. なお、磁化は、Ωの外では0なので、上式は、
    𝐇M(𝐫)=14π𝐬Ωgrad𝐫[𝐌(𝐬)|(𝐫𝐬)|𝐫𝐬|3] d3s
    としてもよい(紛らわしいのであまり推奨しない)。
  8. 8.0 8.1 ここでいうδは、摂動あるいは変分を表し、ディラックのδと紛らわしいが全然別物である。

参考文献

テンプレート:参照方法 テンプレート:Reflist

関連項目

外部リンク

  1. W.K.H. パノフスキー (著), M. フィリップス (著),林 忠四郎 (翻訳), 西田 稔 (翻訳);「新版 電磁気学〈上〉」吉岡書店; POD版 (2002/09)
  2. 2.0 2.1 2.2 2.3 [★溝口 正 (著) 「電磁気学―SI UNITS 」裳華房 (2001/03)](特にP188付近)
  3. [★]竹山 説三 (著) 「電磁気学現象理論」丸善出版; 3版 (1949)
  4. 4.0 4.1 [★]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. テンプレート:PDFlink 電気学会マグネティックス研究会資料 MAG-03 号:54-64 ページ:5-10 発行年:2003年03月31日
  6. 静磁場解析のための二次要素を用いる有限要素法の研究
  7. 有限要素法による磁場解析
  8. 8.0 8.1 磁場が満たす偏微分方程式(秋田高専講義録)
  9. 園田英徳;「大学院生のための基礎物理学」講談社(2011/09/29)ISBN 978-4-06-153277-9
  10. 平川浩正;「電磁気学(新物理学シリーズ2)」培風館 (1986/04) ISBN 9784563024024
  11. 守末 利弥「数値電磁気学のためのゲージ理論」森北出版 (1996/04) ISBN 978-4-627-71600-1
  12. 依田 潔 (著) ;「Mathematicaによる電磁界シミュレーション入門 - POD版 (計算電気・電子工学シリーズ)」森北出版(2012/2/24) ISBN 978-4-627-71529-5
  13. テンプレート:PDFlink
  14. [2]立教大学講義ノート 本講義の参考文献欄 等、至る所に文字化けがあるが、たとえば、もじばけらった等によって解読できる。
  15. テンプレート:PDFlink テンプレート:PDFlink
  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. 高橋 康人 (早稲田大学):テンプレート:PDFlink
  24. 東京電機大学 (編集)「入門 電磁気学」東京電機大学出版局 (2006/03) ISBN 9784501004200
  25. 早川 義晴 (著) 「電気教科書 電験三種合格ガイド」翔泳社 (2011/2/25) ISBN 9784798126623
  26. テンプレート:PDFlink
  27. 坪井 一洋 (著) 「システムと微分方程式」三恵社 (2011/5/22) ISBN 9784883618248
  28. 複数磁石による静磁場(簡易シミュレータ)
  29. テンプレート:PDFlink
  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)