電磁気学 (4回目)

2023/9/6(水)
電磁気学 (4回目)
 
(Electro magnetics)
 
(特殊)相対論的マクスウェルの方程式の導出
(Maxwell's equation)
 
■ 前提
▼ 定義
ナブラ∇ = (∂/∂x,∂/∂y,∂/∂z)
grad f = ∇f , div E = ∇・E , rot E = ∇×E 
ラプラシアン∇・∇ = ∇2 = Δ
ダランべルシアン
□ = ∂μμ = -(1/c2)(∂2/∂t2) + Δ
 
E  :電場(N/C)
B  :磁束密度(T) = (Wb/m2) = (Vs/m2) = (N/(A・m))
ρ :電荷密度[総量Q:電荷(C)]
ε0:真空中の誘電率(F/m)[ε:誘電率]
μ0:真空中の透磁率(N/A2)[μ:透磁率]
j  :電流密度[総量I:電流(A)]
φ :スカラーポテンシャル
A  :ベクトルポテンシャル
c  :真空中の光速度(m/s)[c = 1/√(ε0μ0)導出略]
 
▼ ポテンシャルの方程式
E = -gradφ - ∂A/∂t
B = rotA 
 
{Δ - ε0μ0(∂2/∂t2)}A 
- grad{divA + ε0μ0(∂φ/∂t)} = -μ0j  … ①
Δφ + div(∂A/∂t) = -ρ/ε0  … ②
 
 
■ 導出
▼ 式の変形
c = 1/√(ε0μ0)  … 光速度
より
ε0μ0 = 1/c2 
を使って
{Δ - ε0μ0(∂2/∂t2)}A 
- grad{divA + ε0μ0(∂φ/∂t)} = -μ0j  … ①

{Δ - (1/c2)(∂2/∂t2)}A 
- grad{divA + (1/c2)(∂φ/∂t)} = -μ0j  … ①'
 
Δφ + div(∂A/∂t) = -ρ/ε0  … ②

(1/c2)(∂2/∂t2)φを引いて足すと
{Δ - (1/c2)(∂2/∂t2)}φ
+ div(∂A/∂t) + (1/c2)(∂2φ/∂t2) = -ρ/ε0  
{Δ - (1/c2)(∂2/∂t2)}φ
+ (∂/∂t){divA + (1/c2)(∂φ/∂t)} = -ρ/ε0  
{Δ - (1/c2)(∂2/∂t2)}(φ/c)
+ (1/c)(∂/∂t){divA + (1/c)(∂(φ/c)/∂t)}
= -ρ/(ε0c)  … ②'
に変形
 
{Δ - (1/c2)(∂2/∂t2)}A 
- grad{divA + (1/c2)(∂φ/∂t)} = -μ0j  … ①'
{Δ - (1/c2)(∂2/∂t2)}(φ/c)
+ (1/c)(∂/∂t){divA + (1/c)(∂(φ/c)/∂t)}
= -ρ/(ε0c)  … ②'
 
▼ 電磁ポテンシャルと電流密度の定義
四元ポテンシャルAμの定義(ゲージ場)
Aμ = (A0, A1, A2, A3) = (φ/c, A)
四元電流密度jμの定義
jμ = (j0, j1, j2, j3) = (cρ, j)
 
②'式
{Δ - (1/c2)(∂2/∂t2)}(φ/c)
+ (∂/c∂t){divA + (∂(φ/c)/c∂t)}
= -ρ/(ε0c)  … ②'

{Δ - (1/c2)(∂2/∂t2)}A0 
+ (1/c)(∂/∂t){divA + (1/c)(∂/∂t)A0}
= -cρ/(ε0c2) = -j0/{ε0/(ε0μ0)} = -μ0j0 
 
①'式
{Δ - (1/c2)(∂2/∂t2)}A 
- grad{divA + (1/c2)(∂φ/∂t)} = -μ0j  … ①'

k = 1,2,3として
{Δ - (1/c2)(∂2/∂t2)}Ak 
- (∂/∂xk){divA + (1/c)(∂/∂t)A0} = -μ0jk 
 
▼ 微分の定義
ラプラシアン
Δ = ∂2/∂x2 + ∂2/∂y2 + ∂2/∂z2 
 
ミンコフスキー計量ημν = diag(-1,1,1,1)
 
ナブラ∇
μ = {(1/c)(∂/∂t), ∇}
= {(1/c)(∂/∂t), ∂/∂x, ∂/∂y, ∂/∂z}
μ = ημνμ = {-(1/c)(∂/∂t), ∇}
= {-(1/c)(∂/∂t), ∂/∂x, ∂/∂y, ∂/∂z}
 
ダランべルシアン(d'Alembertian)□
□ = ∂μμ 
= -(1/c2)∂2/∂t2 + ∂2/∂x2 + ∂2/∂y2 + ∂2/∂z2 
= -(1/c2)(∂2/∂t2) + Δ
 
クアブラ(quabla)□
2 = ∂μμ と書く事もある
 
▼ 式をまとめる
{Δ - (1/c2)(∂2/∂t2)}A0 
+ (1/c)(∂/∂t){divA + (1/c)(∂/∂t)A0} = -μ0j0 

k = 1,2,3として
{Δ - (1/c2)(∂2/∂t2)}Ak 
- (∂/∂xk){divA + (1/c)(∂/∂t)A0} = -μ0jk 

divA + (1/c)(∂/∂t)A0 
= (1/c)(∂/∂t)A0 + (∂/∂xk)Ak 
= ∂νAν 
Δ - (1/c2)(∂2/∂t2) = □
を代入して
□A0 + (1/c)(∂/∂t)(∂νAν) = -μ0j0 
□Ak - (∂/∂xk)(∂νAν) = -μ0jk 
これをまとめて
□Aμ - ∂μ(∂νAν) = -μ0jμ 
 
 
▼ ローレンツ条件(次回解説)
divA + ε0μ0(∂φ/∂t) = 0  … ローレンツ条件

νAν = 0  … ローレンツ条件
と書ける
 
 
▼ 電磁場テンソル
□ = ∂μμ 
□Aμ - ∂μ(∂νAν) = -μ0jμ 
より
ννAμ - ∂μ(∂νAν) = -μ0jμ 
ν(∂νAμ - ∂μAν) = -μ0jμ 
ν(∂μAν - ∂νAμ) = μ0jμ 
 
ここで、電磁場テンソルを
Fμν = ∂μAν - ∂νAμ 
と定義する
 
νFμν = μ0jμ 
 
因みに
Fμν = ∂μAν - ∂νAμ = -(∂νAμ - ∂μAν) = -Fνμ 
Fμμ = ∂μAμ - ∂μAμ = 0
 
以下(k=1,2,3)とする
B = rotA = ∇×A 
= (Bx, By, Bz) = (B1, B2, B3) = Bk = ∂k×Bk 
= (∂2A3 - ∂3A2, ∂3A1 - ∂1A3, ∂1A2 - ∂2A1)
E = -gradφ - ∂A/∂t = -∇φ - ∂A/∂t
= -∂A/∂t - ∇φ
Ek = -∂Ak/∂t - ∂φ/∂xk = (Ex, Ey, Ez)
を使って
 
F0k = ∂0Ak-∂kA0 = -(1/c)∂Ak/∂t - ∂φ/∂xk = Ek/c
F12 = ∂1A2-∂2A1 =  B3 
F13 = ∂1A3-∂3A1 = -B2 
F23 = ∂2A3-∂3A2 =  B1 
よって
 
Fμν =
|0      Ex/c  Ey/c  Ez/c|
|-Ex/c   0     Bz   -By |
|-Ey/c  -Bz    0     Bx |
|-Ez/c   By   -Bx    0  |
 
 
■ 結果
▼ 定義
E  :電場(N/C)
B  :磁束密度(T) = (Wb/m2) = (Vs/m2) = (N/(A・m))
ρ :電荷密度[総量Q:電荷(C)]
ε0:真空中の誘電率(F/m)[ε:誘電率]
μ0:真空中の透磁率(N/A2)[μ:透磁率]
j  :電流密度[総量I:電流(A)]
φ :スカラーポテンシャル(V)
A  :ベクトルポテンシャル
c  :真空中の光速度(m/s)[c = 1/√(ε0μ0)導出略]
 
▼ ポテンシャル
E = -gradφ - ∂A/∂t
B = rotA 
 
▼ 四元ベクトル
電磁ポテンシャルAμの定義
Aμ = (A0, A1, A2, A3) = (φ/c, A)
四元電流密度jμの定義
jμ = (j0, j1, j2, j3) = (cρ, j)
 
ナブラ∇ = (∂/∂x, ∂/∂y, ∂/∂z)
ミンコフスキー計量ημν = diag(-1,1,1,1)
μ = {(1/c)(∂/∂t), ∇}
μ = ημνμ = {-(1/c)(∂/∂t), ∇}
ラプラシアンΔ = ∂2/∂x2 + ∂2/∂y2 + ∂2/∂z2 
ダランべルシアン
□ = ∂μμ = -(1/c2)(∂2/∂t2) + Δ
 
クアブラ(quabla)□
2 = ∂μμ と書く事もある
 
▼ 相対論的マクスウェルの方程式
□Aμ - ∂μ(∂νAν) = -μ0jμ 
または
νFμν = μ0jμ 
 
▼ ローレンツ条件(次回解説)
νAν = 0
 
▼ 電磁場テンソルFμν 
Fμν = ∂μAν - ∂νAμ 
 
(Bx, By, Bz) = B = rotA 
(Ex, Ey, Ez) = E = -gradφ - ∂A/∂t
 
Fμν =
|0      Ex/c  Ey/c  Ez/c|
|-Ex/c   0     Bz   -By |
|-Ey/c  -Bz    0     Bx |
|-Ez/c   By   -Bx    0  |
 

このブログの人気の投稿

NEWS

N88-BASICでゲーム (1回目)