Processing math: 100%

四方山話

目次を表示

平面と球面での計算(2025年5月11日)

具体例で少し計算しておきましょう:R を求めるためには Rij,kl が必要で

Rij,kl=kΓijllΓijkΓimlΓmjk+ΓimkΓmjl

であり

¯rjk=Γijkri

Γijk の定義式です。ここで ¯rij は、ベクトル

rij=jri=jir

を面に接するベクトルと垂直なベクトルの和に分解するときの前者を表しました。

 まず Γijk を計算しやすい形に書き直すことを考えてみます。

 rjk の面ベクトル成分を求めることが絶対必要なように思いますが、面に垂直な成分と rl との内積が 0 であることに気付くなら、次のような少しずるい方法でも何とかなるのでは、と考えてしまいます。(1) 式

¯rjk=Γijkri

の両辺と rl との内積をとるのです。

rjkrl=Γijkgil=Γl,ij

が得られます。これは

(krj)rl=Γl,ij

を意味します。左辺は、3次元ユークリッド空間の原点 O と曲面上の点 (x,y) を結ぶベクトル r(x,y) がわかれば、

ri=ir

は求まるはずですから、計算できそうです。従って

Γkij=gklΓl,ij

によって Γkij を計算できます。

 ただ、数学者としてはもう少し整理した形で書いてみたくなります。単なる美意識の問題と言えば言えるのですが、おそらくそうしておく方が後で役に立つと思われるからです。

 ここまでの議論では gij が基本的な役割を果たしてきました。ここでもそれを表に出す工夫をしてみます。

k(rjrl)=(krj)rl+rj(krl)

です。従って

kgjl=Γl,jk+Γj,lk

です。右辺の Γl,jk だけを残したいので、等式を増やすために添え字を順ぐりに回すと

jglk=Γk,lj+Γl,kj

lgkj=Γj,kl+Γk,jl

が得られます。ここで

Γi,jk=Γi,kj

に注意するならば

(2)+(3)(4)

によって右辺に Γl,jk だけが残るようにできます。つまり

12(lgjk+jgkl+kglj)=Γl,jk

です。ただ、gij=gji を使って添え字を少し見やすく配置し直しました。よって

Γijk=12gil(lgjk+jgkl+kglj)

が得られます。

この形なら gij さえ求めれば良いから、確かに見通しがいいですね。具体例について Γijk を計算してみましょう。平面、球面の2例を取り上げます:

 平面の場合がもちろん最も簡単でガウス座標としては、直交座標系を1つ決めて (x,y) とします。このとき

r(x,y)=(x,y,0)

としてよいので

r1(x,y)=xr(x,y)=(1,0,0),r2(x,y)=yr(x,y)=(0,1,0)

であり、よって

g11=g22=1,g12=0

が得られます。よって、(6) によれば、すべての i,j,k について

Γijk=12gil(lgjk+jgkl+kglj)=0

です。

 よって、Rij,kl の定義式

Rij,kl=kΓijllΓijkΓimlΓmjk+ΓimkΓmjl

から、これらすべて Rij,kl=0 です。よって

Rij=Rki,kj=0,R=gijRij=0

でもあります。

 半径 r の球面上の点(但し、北極と南極は除外しますが)の座標としては、Z 軸からの角度 θX 軸からの角度 ϕ の組 (θ,ϕ) を用います。このとき

r(θ,ϕ)=(rsinθcosϕ,rsinθsinϕ,rcosθ)

であって、

rθ=(rcosθcosϕ,rcosθsinϕ,rsinθ)rϕ=(rsinθsinϕ,rsinθcosϕ,0)

となります。よって

gθθ=r2,gθϕ=gϕθ=0,gϕϕ=r2sin2θ

です。

 これらを (5) 式

Γl,jk=12(lgjk+jgkl+kglj)

へ代入すれば

Γθ,ϕϕ=r2sinθcosθ,Γϕ,θϕ=Γϕ,ϕθ=r2sinθcosθ,0

が得られます。

 また、行列

(gθθgθϕgϕθgϕϕ)=(r200r2sin2θ)

の逆行列を求めることにより

gθθ=r2,gθϕ=gϕθ=0,gϕϕ=r2sin2θ

となり、よって (5) 式

Γl,jk=12(lgjk+jgkl+kglj)

を丹念に計算することにより、最終的に

Γθϕϕ=sinθcosθ,Γϕθϕ=Γϕϕθ=cosθsinθ,0

が得られます。

 これらから

R12,12=r2sin2θ

が導かれ

g=r4sin2θ

であり、

R=2r2

です。従って

K=12R=1r2

となります。これは最初の出発点で選んだ K の値に一致しています。

目次を表示

このページの先頭へ