曲率テンソル (2回目)

2023/8/4(金)
曲率テンソル (2回目)
 
(Curvature tensor)
 
リーマン・クリストッフェル曲率テンソル
(Riemann-Christoffel curvature tensor)
 
球面の
曲率テンソルRnm,αβ 
リッチテンソルRmβ 
スカラー曲率R
を求める
 
■ 定義
▼ 曲率テンソル
基底ベクトルの2階共変微分の
1回目と2回目の微分の順番を入れ替えたときの
差の係数をリーマン曲率テンソルとする
Rnm,αβen = ∇α(∇βem) - ∇β(∇αem) = -Rnm,βα 
 
Rnm,αβ 
= (∂/∂xαnβm - (∂/∂xβnαm 
+ ΓnαkΓkβm  - ΓnβkΓkαm 
 
▼ 極座標系のクリストッフェル記号
 
https://ulprojectmail.blogspot.com/2023/07/christoffel-2.html
クリストッフェル記号 (2回目)
より
 
クリストッフェル記号
Γγαβeγ = (∂/∂xβ)eα = ∇βeα =
Γγβαeγ 
 
Γθθ = -r
Γφφ = -rsin2θ
Γθrθ = Γθθr = Γφrφ = Γφφr = 1/r
Γθφφ = -sinθcosθ
Γφθφ = Γφφθ = cosθ/sinθ
 
▼ 球面のクリストッフェル記号
半径rの球面の変数はθ,φなので
rを固定するとrに関する微分は0になる
よって0以外の成分は次の通り
 
Γθφφ = -sinθcosθ
Γφθφ = cosθ/sinθ
Γφφθ = cosθ/sinθ
 
 
■ (∂/∂xαnβm 
▼ 導出
(∂/∂θ)Γθφφ = (∂/∂θ)(-sinθcosθ)
= -cosθcosθ + sinθsinθ = -cos2θ+sin2θ
= -(cos2θ-sin2θ)
 
(∂/∂θ)Γφθφ = (∂/∂θ)Γφφθ =
(∂/∂θ)(cosθ/sinθ)
= -sinθ/sinθ + cosθ(-cosθ/sin2θ)
= -(sin2θ + cos2θ)/sin2θ
= -1/sin2θ
 
▼ まとめ
(∂/∂θ)Γθφφ =  sin2θ-cos2θ
(∂/∂θ)Γφθφ = -1/sin2θ
(∂/∂θ)Γφφθ = -1/sin2θ
 
■ 微分項Dの計算
D = (∂/∂xαnβm - (∂/∂xβnαm 
α⇔β なら -D、α=β なら D = 0 になる
 
▼ n=θ
(∂/∂θ)Γθφφ =  sin2θ-cos2θ
D = (∂/∂xαθβm - (∂/∂xβθαm 
α=θ,β=φ,m=φ,D =  sin2θ-cos2θ
 
▼ n=φ
(∂/∂θ)Γφθφ = -1/sin2θ
(∂/∂θ)Γφφθ = -1/sin2θ
D = (∂/∂xαφβm - (∂/∂xβφαm 
α=θ,β=φ,m=θ,D = -1/sin2θ
α=θ,β=θ,m=φ,D = 0
 
▼ D
n=θ,α=θ,β=φ,m=φ,D =  sin2θ-cos2θ
n=φ,α=θ,β=φ,m=θ,D = -1/sin2θ
 
 
■ 通常項Eの計算
E = ΓnαkΓkβm  - ΓnβkΓkαm 
α⇔β なら -D
α=β なら D = 0 になる
 
▼ n=θ,k=φ
E = ΓnαkΓkβm  - ΓnβkΓkαm 
Γθφφ = -sinθcosθ
Γφθφ = cosθ/sinθ
Γφφθ = cosθ/sinθ
α=φ
E = ΓθφφΓφβm  - ΓθβφΓφφm 
β=θ,m=φ,E = -sinθcosθcosθ/sinθ = -cos2θ
 
n=θ,k=φ,α=φ,β=θ,m=φ,E = -cos2θ
 
▼ n=φ,k=θ
E = ΓnαkΓkβm  - ΓnβkΓkαm 
Γθφφ = -sinθcosθ
Γφθφ = cosθ/sinθ
Γφφθ = cosθ/sinθ
α=φ
E = ΓφφθΓθβm  - ΓφβθΓθφm 
β=φ,m=φ,E = 0
 
n=φ,k=θ,α=φ,β=φ,m=φ,E = 0
 
▼ n=φ,k=φ
E = ΓnαkΓkβm  - ΓnβkΓkαm 
Γθφφ = -sinθcosθ
Γφθφ = cosθ/sinθ
Γφφθ = cosθ/sinθ
α=θ
E = ΓφθφΓφβm  - ΓφβφΓφθm 
β=θ,m=φ,E = 0
β=φ,m=θ,E = (cosθ/sinθ)(cosθ/sinθ) = cos2θ/sin2θ
 
n=φ,k=φ,α=θ,β=φ,m=θ,E = cos2θ/sin2θ
 
▼ Eのkの縮約(Σ[k=r,θ,φ])
E = ΓnαkΓkβm  - ΓnβkΓkαm 
n=θ,k=φ,α=φ,β=θ,m=φ,E = -cos2θ
n=φ,k=φ,α=θ,β=φ,m=θ,E = cos2θ/sin2θ
 
n=θ,α=φ,β=θ,m=φ,E = -cos2θ
n=φ,α=θ,β=φ,m=θ,E = cos2θ/sin2θ
 
 
■ 球面の曲率テンソルRnm,αβ 
▼ D,E
D = (∂/∂xαn - (∂/∂xβn 
n=θ,α=θ,β=φ,m=φ,D =  sin2θ-cos2θ
n=φ,α=θ,β=φ,m=θ,D = -1/sin2θ
 
E = ΓnαkΓkβm  - ΓnβkΓkαm 
n=θ,α=φ,β=θ,m=φ,E = -cos2θ
n=φ,α=θ,β=φ,m=θ,E = cos2θ/sin2θ
 
▼ K=D+E
Rnm,αβ 
 
n=θ,α=θ,β=φ,m=φ,D =  sin2θ-cos2θ
n=θ,α=φ,β=θ,m=φ,E = -cos2θ
n=θ,α=θ,β=φ,m=φ,K =  sin2θ-cos2θ+cos2θ
=  sin2θ
 
Rθφ,θφ =  sin2θ
Rθφ,φθ = -sin2θ
 
n=φ,α=θ,β=φ,m=θ,D = -1/sin2θ
n=φ,α=θ,β=φ,m=θ,E = cos2θ/sin2θ
n=φ,α=θ,β=φ,m=θ,K = -1/sin2θ+cos2θ/sin2θ
=  (cos2θ-1)/sin2θ = -sin2θ/sin2θ = -1
 
Rφθ,θφ = -1
Rφθ,φθ =  1
 
 
■ リッチテンソルRmβ の導出
曲率テンソルを縮約したもの
R = Rαm,αβ (Σ[α=θ,φ])
 
球面の曲率テンソルRnm,αβ 
Rnm,αβ 
= (∂/∂xαnβm - (∂/∂xβnαm 
+ ΓnαkΓkβm  - ΓnβkΓkαm 
 
Rθφ,θφ =  sin2θ
Rθφ,φθ = -sin2θ
Rφθ,θφ = -1
Rφθ,φθ =  1
 
Rθθ = Rθθ,θθ + Rφθ,φθ = 0 + 1
Rθφ = Rθθ,θφ + Rφθ,φφ = 0 + 0
Rφθ = Rθφ,θθ + Rφφ,φθ = 0 + 0
Rφφ = Rθφ,θφ + Rφφ,φφ = sin2θ + 0
より
 
Rθθ = 1
Rφφ = sin2θ
 
 
■ スカラー曲率(リッチスカラー)
gαβ = diag(1, 1/r2, 1/(r2sin2θ))
  |1 0    0           |
= |0 1/r2 0           |
  |1 0    1/(r2sin2θ)|
Rθθ = 1
Rφφ = sin2θ
 
R = gR 
= gθθRθθ + gφφRφφ 
= (1/r2)1 + {1/(r2sin2θ)}sin2θ
= 1/r2 + 1/r2 
= 2/r2 
 
 
■ 結果
▼ 球面の曲率テンソルRnm,αβ 
Rnm,αβ 
= (∂/∂xαnβm - (∂/∂xβnαm 
+ ΓnαkΓkβm  - ΓnβkΓkαm 
 
Rθφ,θφ =  sin2θ
Rθφ,φθ = -sin2θ
Rφθ,θφ = -1
Rφθ,φθ =  1
 
▼ リッチテンソル
曲率テンソルを縮約したもの
R = Rαm,αβ (Σ[α=θ,φ])
 
Rθθ = 1
Rφφ = sin2θ
 
▼ スカラー曲率(リッチスカラー)
gαβ = diag(1, 1/r2, 1/(r2sin2θ))
  |1 0    0           |
= |0 1/r2 0           |
  |1 0    1/(r2sin2θ)|
Rθθ = 1
Rφφ = sin2θ
 
R = gR = gθθRθθ + gφφRφφ 
= 2/r2 


このブログの人気の投稿

NEWS

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