曲線座標におけるテンソル

曲線座標はテンソル計算で定式化することができ物理学工学、特に流体力学連続体力学における物理量の輸送や物質の変形を記述するために重要な用途があります

3次元曲線座標におけるベクトル代数とテンソル代数

曲線座標系における初等的なベクトル代数とテンソル代数は、力学物理学の古い科学文献の一部で用いられており、1900年代初頭から中期にかけての研究、例えばGreenとZernaの著書[1]を理解する上で不可欠なものとなり 得る。本節では、曲線座標系におけるベクトル代数と2階テンソル代数における有用な関係式をいくつか示す。表記法と内容は主にOgden、[2] Naghdi、[3] Simmonds、[4] GreenとZerna、[1] BasarとWeichert、[5] Ciarlet [6]に拠っている。

座標変換

座標変数 と を持つ2つの座標系を考えてみましょう。これらはそれぞれ と略記し、添え字は常に1から3までとします。これらの座標系は3次元ユークリッド空間に埋め込まれていると仮定します。座標 と は互いに説明するために使用できます。なぜなら、一方の座標系で座標線に沿って移動すると、もう一方の座標系を使用して位置を記述できるからです。このように、座標 と は互いの関数です 。 ( Z 1 , Z 2 , Z 3 ) {\displaystyle (Z^{1},Z^{2},Z^{3})} ( Z 1 ´ , Z 2 ´ , Z 3 ´ ) {\displaystyle (Z^{\acute {1}},Z^{\acute {2}},Z^{\acute {3}})} Z i {\displaystyle Z^{i}} Z i ´ {\displaystyle Z^{\acute {i}}} i {\displaystyle i} Z i {\displaystyle Z^{i}} Z i ´ {\displaystyle Z^{\acute {i}}} Z i {\displaystyle Z^{i}} Z i ´ {\displaystyle Z^{\acute {i}}}

Z i = f i ( Z 1 ´ , Z 2 ´ , Z 3 ´ ) {\displaystyle Z^{i}=f^{i}(Z^{\acute {1}},Z^{\acute {2}},Z^{\acute {3}})} のために i = 1 , 2 , 3 {\displaystyle i=1,2,3}

これは次のように書ける。

Z i = Z i ( Z 1 ´ , Z 2 ´ , Z 3 ´ ) = Z i ( Z i ´ ) {\displaystyle Z^{i}=Z^{i}(Z^{\acute {1}},Z^{\acute {2}},Z^{\acute {3}})=Z^{i}(Z^{\acute {i}})} のために i ´ , i = 1 , 2 , 3 {\displaystyle {\acute {i}},i=1,2,3}

これら3つの方程式は、 からへ の座標変換とも呼ばれます。この変換を と表記します。したがって、座標変数を含む座標系から座標を含む座標系への変換は、次のように表されます Z i ´ {\displaystyle Z^{\acute {i}}} Z i {\displaystyle Z^{i}} T {\displaystyle T} Z i ´ {\displaystyle Z^{\acute {i}}} Z i {\displaystyle Z^{i}}

Z = T ( z ´ ) {\displaystyle Z=T({\acute {z}})}

同様にを の関数として次のように表すことができます Z i ´ {\displaystyle Z^{\acute {i}}} Z i {\displaystyle Z^{i}}

Z i ´ = g i ´ ( Z 1 , Z 2 , Z 3 ) {\displaystyle Z^{\acute {i}}=g^{\acute {i}}(Z^{1},Z^{2},Z^{3})} のために i ´ = 1 , 2 , 3 {\displaystyle {\acute {i}}=1,2,3}

そして自由方程式は次のように簡潔に書くことができる。

Z i ´ = Z i ´ ( Z 1 , Z 2 , Z 3 ) = Z i ´ ( Z i ) {\displaystyle Z^{\acute {i}}=Z^{\acute {i}}(Z^{1},Z^{2},Z^{3})=Z^{\acute {i}}(Z^{i})} のために i ´ , i = 1 , 2 , 3 {\displaystyle {\acute {i}},i=1,2,3}

これら3つの方程式は、 からへの座標変換とも呼ばれます。この変換を と表記します。座標変数を含む座標系から座標を含む座標系への変換は、次のように表されます。 Z i {\displaystyle Z^{i}} Z i ´ {\displaystyle Z^{\acute {i}}} S {\displaystyle S} Z i {\displaystyle Z^{i}} Z i ´ {\displaystyle Z^{\acute {i}}}

z ´ = S ( z ) {\displaystyle {\acute {z}}=S(z)}

変換が全単射である場合、変換の像、すなわち を の許容座標の集合と呼びますが線形である場合、座標系はアフィン座標系と呼ばれ、そうでない場合は曲線座標系と呼ばれます T {\displaystyle T} Z i {\displaystyle Z^{i}} Z i ´ {\displaystyle Z^{\acute {i}}} T {\displaystyle T} Z i {\displaystyle Z^{i}} Z i {\displaystyle Z^{i}}

ヤコビアン

座標 と は互いの関数であることがわかったので、座標変数 を座標変数 に関して微分することができます Z i {\displaystyle Z^{i}} Z i ´ {\displaystyle Z^{\acute {i}}} Z i {\displaystyle Z^{i}} Z i ´ {\displaystyle Z^{\acute {i}}}

考慮する

Z i Z i ´ = d e f J i ´ i {\displaystyle {\frac {\partial {Z^{i}}}{\partial {Z^{\acute {i}}}}}\;{\overset {\underset {\mathrm {def} }{}}{=}}\;J_{\acute {i}}^{i}} の場合、これらの導関数は行列 に配置することができ、例えば は- 行- 列 の要素である。 i ´ , i = 1 , 2 , 3 {\displaystyle {\acute {i}},i=1,2,3} J {\displaystyle J} J i ´ i {\displaystyle J_{\acute {i}}^{i}} i {\displaystyle i} i ´ {\displaystyle {\acute {i}}}

J = ( J 1 ´ 1 J 2 ´ 1 J 3 ´ 1 J 1 ´ 2 J 2 ´ 2 J 3 ´ 2 J 1 ´ 3 J 2 ´ 3 J 3 ´ 3 ) = ( Z 1 Z 1 ´ Z 1 Z 2 ´ Z 1 Z 3 ´ Z 2 Z 1 ´ Z 2 Z 2 ´ Z 2 Z 3 ´ Z 3 Z 1 ´ Z 3 Z 2 ´ Z 3 Z 3 ´ ) {\displaystyle J={\begin{pmatrix}J_{\acute {1}}^{1}&J_{\acute {2}}^{1}&J_{\acute {3}}^{1}\\J_{\acute {1}}^{2}&J_{\acute {2}}^{2}&J_{\acute {3}}^{2}\\J_{\acute {1}}^{3}&J_{\acute {2}}^{3}&J_{\acute {3}}^{3}\end{pmatrix}}={\begin{pmatrix}{\partial {Z^{1}} \over \partial {Z^{\acute {1}}}}&{\partial {Z^{1}} \over \partial {Z^{\acute {2}}}}&{\partial {Z^{1}} \over \partial {Z^{\acute {3}}}}\\{\partial {Z^{2}} \over \partial {Z^{\acute {1}}}}&{\partial {Z^{2}} \over \partial {Z^{\acute {2}}}}&{\partial {Z^{2}} \over \partial {Z^{\acute {3}}}}\\{\partial {Z^{3}} \over \partial {Z^{\acute {1}}}}&{\partial {Z^{3}} \over \partial {Z^{\acute {2}}}}&{\partial {Z^{3}} \over \partial {Z^{\acute {3}}}}\end{pmatrix}}}

結果として得られる行列はヤコビ行列と呼ばれます。

曲線座標のベクトル

を3次元ユークリッド空間の任意の基底とする。一般に、基底ベクトルは単位ベクトルでも互いに直交でもない。しかし、それらは線形独立であることが要求される。すると、ベクトルは[4]のように表される。27  成分はベクトルの反変成分である ( b 1 , b 2 , b 3 ) {\displaystyle (\mathbf {b} _{1},\mathbf {b} _{2},\mathbf {b} _{3})} v {\displaystyle \mathbf {v} } v = v k b k {\displaystyle \mathbf {v} =v^{k}\,\mathbf {b} _{k}} v k {\displaystyle v^{k}} v {\displaystyle \mathbf {v} }

基底 は関係式[4]によって定義される:28–29  ここでクロネッカーのデルタである ( b 1 , b 2 , b 3 ) {\displaystyle (\mathbf {b} ^{1},\mathbf {b} ^{2},\mathbf {b} ^{3})} b i b j = δ j i {\displaystyle \mathbf {b} ^{i}\cdot \mathbf {b} _{j}=\delta _{j}^{i}} δ j i {\displaystyle \delta _{j}^{i}}

ベクトルは逆基底で表すこともできます。 成分はベクトルの共変成分です v {\displaystyle \mathbf {v} } v = v k   b k {\displaystyle \mathbf {v} =v_{k}~\mathbf {b} ^{k}} v k {\displaystyle v_{k}} v {\displaystyle \mathbf {v} }

曲線座標における2階テンソル

2 次テンソルは次のように表現できます。 成分は、 2 次テンソルの 反変成分、混合右共変成分、混合左共変成分、成分と呼ばれます。 S = S i j   b i b j = S   j i   b i b j = S i   j   b i b j = S i j   b i b j {\displaystyle {\boldsymbol {S}}=S^{ij}~\mathbf {b} _{i}\otimes \mathbf {b} _{j}=S_{~j}^{i}~\mathbf {b} _{i}\otimes \mathbf {b} ^{j}=S_{i}^{~j}~\mathbf {b} ^{i}\otimes \mathbf {b} _{j}=S_{ij}~\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}} S i j {\displaystyle S^{ij}} S   j i {\displaystyle S_{~j}^{i}} S i   j {\displaystyle S_{i}^{~j}} S i j {\displaystyle S_{ij}}

計量テンソルと成分間の関係

量はのように定義される[4] :39  g i j {\displaystyle g_{ij}} g i j {\displaystyle g^{ij}}

g i j = b i b j = g j i   ;     g i j = b i b j = g j i {\displaystyle g_{ij}=\mathbf {b} _{i}\cdot \mathbf {b} _{j}=g_{ji}~;~~g^{ij}=\mathbf {b} ^{i}\cdot \mathbf {b} ^{j}=g^{ji}} 上記の式から、 v i = g i k   v k   ;     v i = g i k   v k   ;     b i = g i j   b j   ;     b i = g i j   b j {\displaystyle v^{i}=g^{ik}~v_{k}~;~~v_{i}=g_{ik}~v^{k}~;~~\mathbf {b} ^{i}=g^{ij}~\mathbf {b} _{j}~;~~\mathbf {b} _{i}=g_{ij}~\mathbf {b} ^{j}}

ベクトルの成分は[4]によって関係付けられる:30–32  また、 v b i = v k   b k b i = v k   δ k i = v i {\displaystyle \mathbf {v} \cdot \mathbf {b} ^{i}=v^{k}~\mathbf {b} _{k}\cdot \mathbf {b} ^{i}=v^{k}~\delta _{k}^{i}=v^{i}} v b i = v k   b k b i = v k   δ i k = v i {\displaystyle \mathbf {v} \cdot \mathbf {b} _{i}=v_{k}~\mathbf {b} ^{k}\cdot \mathbf {b} _{i}=v_{k}~\delta _{i}^{k}=v_{i}} v b i = v k   b k b i = g k i   v k {\displaystyle \mathbf {v} \cdot \mathbf {b} _{i}=v^{k}~\mathbf {b} _{k}\cdot \mathbf {b} _{i}=g_{ki}~v^{k}} v b i = v k   b k b i = g k i   v k {\displaystyle \mathbf {v} \cdot \mathbf {b} ^{i}=v_{k}~\mathbf {b} ^{k}\cdot \mathbf {b} ^{i}=g^{ki}~v_{k}}

2階テンソルの成分は次のように関係している。 S i j = g i k   S k   j = g j k   S   k i = g i k   g j l   S k l {\displaystyle S^{ij}=g^{ik}~S_{k}^{~j}=g^{jk}~S_{~k}^{i}=g^{ik}~g^{jl}~S_{kl}}

交代テンソル

直交右手基底では、三階交代テンソルは次のように定義されます。 一般の曲線基底では、同じテンソルは次のように表すことができます。 したがって 同様 に、 E = ε i j k   e i e j e k {\displaystyle {\boldsymbol {\mathcal {E}}}=\varepsilon _{ijk}~\mathbf {e} ^{i}\otimes \mathbf {e} ^{j}\otimes \mathbf {e} ^{k}} E = E i j k   b i b j b k = E i j k   b i b j b k {\displaystyle {\boldsymbol {\mathcal {E}}}={\mathcal {E}}_{ijk}~\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}\otimes \mathbf {b} ^{k}={\mathcal {E}}^{ijk}~\mathbf {b} _{i}\otimes \mathbf {b} _{j}\otimes \mathbf {b} _{k}} E i j k = [ b i , b j , b k ] = ( b i × b j ) b k   ;     E i j k = [ b i , b j , b k ] {\displaystyle {\mathcal {E}}_{ijk}=\left[\mathbf {b} _{i},\mathbf {b} _{j},\mathbf {b} _{k}\right]=(\mathbf {b} _{i}\times \mathbf {b} _{j})\cdot \mathbf {b} _{k}~;~~{\mathcal {E}}^{ijk}=\left[\mathbf {b} ^{i},\mathbf {b} ^{j},\mathbf {b} ^{k}\right]} b i × b j = J   ε i j p   b p = g   ε i j p   b p {\displaystyle \mathbf {b} _{i}\times \mathbf {b} _{j}=J~\varepsilon _{ijp}~\mathbf {b} ^{p}={\sqrt {g}}~\varepsilon _{ijp}~\mathbf {b} ^{p}} E i j k = J   ε i j k = g   ε i j k {\displaystyle {\mathcal {E}}_{ijk}=J~\varepsilon _{ijk}={\sqrt {g}}~\varepsilon _{ijk}} E i j k = 1 J   ε i j k = 1 g   ε i j k {\displaystyle {\mathcal {E}}^{ijk}={\cfrac {1}{J}}~\varepsilon ^{ijk}={\cfrac {1}{\sqrt {g}}}~\varepsilon ^{ijk}}

ベクトル演算

アイデンティティマップ

によって定義される恒等写像は次のように示される: [4] : 39  I {\displaystyle \mathbf {I} } I v = v {\displaystyle \mathbf {I} \cdot \mathbf {v} =\mathbf {v} }

I = g i j b i b j = g i j b i b j = b i b i = b i b i {\displaystyle \mathbf {I} =g^{ij}\mathbf {b} _{i}\otimes \mathbf {b} _{j}=g_{ij}\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}=\mathbf {b} _{i}\otimes \mathbf {b} ^{i}=\mathbf {b} ^{i}\otimes \mathbf {b} _{i}}

スカラー(ドット)積

曲線座標における2つのベクトルのスカラー積は[4] :32 である。

u v = u i v i = u i v i = g i j u i v j = g i j u i v j {\displaystyle \mathbf {u} \cdot \mathbf {v} =u^{i}v_{i}=u_{i}v^{i}=g_{ij}u^{i}v^{j}=g^{ij}u_{i}v_{j}}

ベクトル(クロス)積

2つのベクトルの外積次のように与えられる: [4] : 32–34 

u × v = ε i j k u j v k e i {\displaystyle \mathbf {u} \times \mathbf {v} =\varepsilon _{ijk}u_{j}v_{k}\mathbf {e} _{i}}

ここで置換記号は直交座標基底ベクトルです。曲線座標では、等価な式は次のようになります。 ε i j k {\displaystyle \varepsilon _{ijk}} e i {\displaystyle \mathbf {e} _{i}}

u × v = [ ( b m × b n ) b s ] u m v n b s = E s m n u m v n b s {\displaystyle \mathbf {u} \times \mathbf {v} =[(\mathbf {b} _{m}\times \mathbf {b} _{n})\cdot \mathbf {b} _{s}]u^{m}v^{n}\mathbf {b} ^{s}={\mathcal {E}}_{smn}u^{m}v^{n}\mathbf {b} ^{s}}

ここでは3階交代テンソルです。2つのベクトルの 外積は次のように与えられます。 E i j k {\displaystyle {\mathcal {E}}_{ijk}}

u × v = ε i j k u ^ j v ^ k e i {\displaystyle \mathbf {u} \times \mathbf {v} =\varepsilon _{ijk}{\hat {u}}_{j}{\hat {v}}_{k}\mathbf {e} _{i}}

ここでは置換記号でありは直交座標基底ベクトルである。したがって、 ε i j k {\displaystyle \varepsilon _{ijk}} e i {\displaystyle \mathbf {e} _{i}}

e p × e q = ε i p q e i {\displaystyle \mathbf {e} _{p}\times \mathbf {e} _{q}=\varepsilon _{ipq}\mathbf {e} _{i}}

そして

b m × b n = x q m × x q n = ( x p e p ) q m × ( x q e q ) q n = x p q m x q q n e p × e q = ε i p q x p q m x q q n e i . {\displaystyle \mathbf {b} _{m}\times \mathbf {b} _{n}={\frac {\partial \mathbf {x} }{\partial q^{m}}}\times {\frac {\partial \mathbf {x} }{\partial q^{n}}}={\frac {\partial (x_{p}\mathbf {e} _{p})}{\partial q^{m}}}\times {\frac {\partial (x_{q}\mathbf {e} _{q})}{\partial q^{n}}}={\frac {\partial x_{p}}{\partial q^{m}}}{\frac {\partial x_{q}}{\partial q^{n}}}\mathbf {e} _{p}\times \mathbf {e} _{q}=\varepsilon _{ipq}{\frac {\partial x_{p}}{\partial q^{m}}}{\frac {\partial x_{q}}{\partial q^{n}}}\mathbf {e} _{i}.}

したがって、

( b m × b n ) b s = ε i p q x p q m x q q n x i q s {\displaystyle (\mathbf {b} _{m}\times \mathbf {b} _{n})\cdot \mathbf {b} _{s}=\varepsilon _{ipq}{\frac {\partial x_{p}}{\partial q^{m}}}{\frac {\partial x_{q}}{\partial q^{n}}}{\frac {\partial x_{i}}{\partial q^{s}}}}

ベクトル積に戻り、次の関係を使用します。

u ^ j = x j q m u m , v ^ k = x k q n v n , e i = x i q s b s , {\displaystyle {\hat {u}}_{j}={\frac {\partial x_{j}}{\partial q^{m}}}u^{m},\quad {\hat {v}}_{k}={\frac {\partial x_{k}}{\partial q^{n}}}v^{n},\quad \mathbf {e} _{i}={\frac {\partial x_{i}}{\partial q^{s}}}\mathbf {b} ^{s},}

結果は次のようになります:

u × v = ε i j k u ^ j v ^ k e i = ε i j k x j q m x k q n x i q s u m v n b s = [ ( b m × b n ) b s ] u m v n b s = E s m n u m v n b s {\displaystyle \mathbf {u} \times \mathbf {v} =\varepsilon _{ijk}{\hat {u}}_{j}{\hat {v}}_{k}\mathbf {e} _{i}=\varepsilon _{ijk}{\frac {\partial x_{j}}{\partial q^{m}}}{\frac {\partial x_{k}}{\partial q^{n}}}{\frac {\partial x_{i}}{\partial q^{s}}}u^{m}v^{n}\mathbf {b} ^{s}=[(\mathbf {b} _{m}\times \mathbf {b} _{n})\cdot \mathbf {b} _{s}]u^{m}v^{n}\mathbf {b} ^{s}={\mathcal {E}}_{smn}u^{m}v^{n}\mathbf {b} ^{s}}

テンソル演算

アイデンティティマップ

によって定義される恒等写像は 次のように示される[4] :39  I {\displaystyle {\mathsf {I}}} I v = v {\displaystyle {\mathsf {I}}\cdot \mathbf {v} =\mathbf {v} }

I = g i j b i b j = g i j b i b j = b i b i = b i b i {\displaystyle {\mathsf {I}}=g^{ij}\mathbf {b} _{i}\otimes \mathbf {b} _{j}=g_{ij}\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}=\mathbf {b} _{i}\otimes \mathbf {b} ^{i}=\mathbf {b} ^{i}\otimes \mathbf {b} _{i}}

2階テンソルのベクトルへの作用

この作用は曲線座標で次のように表される。 v = S u {\displaystyle \mathbf {v} ={\boldsymbol {S}}\mathbf {u} }

v i b i = S i j u j b i = S j i u j b i ; v i b i = S i j u j b i = S i j u j b i {\displaystyle v^{i}\mathbf {b} _{i}=S^{ij}u_{j}\mathbf {b} _{i}=S_{j}^{i}u^{j}\mathbf {b} _{i};\qquad v_{i}\mathbf {b} ^{i}=S_{ij}u^{j}\mathbf {b} ^{i}=S_{i}^{j}u_{j}\mathbf {b} ^{i}}

2つの2階テンソルの内積

2つの2階テンソルの内積曲線座標で次のように表すことができます。 U = S T {\displaystyle {\boldsymbol {U}}={\boldsymbol {S}}\cdot {\boldsymbol {T}}}

U i j b i b j = S i k T . j k b i b j = S i . k T k j b i b j {\displaystyle U_{ij}\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}=S_{ik}T_{.j}^{k}\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}=S_{i}^{.k}T_{kj}\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}}

あるいは、

U = S i j T . n m g j m b i b n = S . m i T . n m b i b n = S i j T j n b i b n {\displaystyle {\boldsymbol {U}}=S^{ij}T_{.n}^{m}g_{jm}\mathbf {b} _{i}\otimes \mathbf {b} ^{n}=S_{.m}^{i}T_{.n}^{m}\mathbf {b} _{i}\otimes \mathbf {b} ^{n}=S^{ij}T_{jn}\mathbf {b} _{i}\otimes \mathbf {b} ^{n}}

2階テンソルの行列式

が2階テンソルの場合、行列式は次の関係で定義される。 S {\displaystyle {\boldsymbol {S}}}

[ S u , S v , S w ] = det S [ u , v , w ] {\displaystyle \left[{\boldsymbol {S}}\mathbf {u} ,{\boldsymbol {S}}\mathbf {v} ,{\boldsymbol {S}}\mathbf {w} \right]=\det {\boldsymbol {S}}\left[\mathbf {u} ,\mathbf {v} ,\mathbf {w} \right]}

ここで任意のベクトルであり、 u , v , w {\displaystyle \mathbf {u} ,\mathbf {v} ,\mathbf {w} }

[ u , v , w ] := u ( v × w ) . {\displaystyle \left[\mathbf {u} ,\mathbf {v} ,\mathbf {w} \right]:=\mathbf {u} \cdot (\mathbf {v} \times \mathbf {w} ).}

曲線基底ベクトルと直交基底ベクトルの関係

を、対象となるユークリッド空間の通常の直交座標基底ベクトルとし、を に写像する2階変換テンソルとします すると この関係 式から、 あることが示せます 。を のヤコビアンとします。すると、行列式の定義から、 となるので、 上記 関係式を用いて、多くの興味深い結果を導き出すことができます。 ( e 1 , e 2 , e 3 ) {\displaystyle (\mathbf {e} _{1},\mathbf {e} _{2},\mathbf {e} _{3})} b i = F e i {\displaystyle \mathbf {b} _{i}={\boldsymbol {F}}\mathbf {e} _{i}} F {\displaystyle {\boldsymbol {F}}} e i {\displaystyle \mathbf {e} _{i}} b i {\displaystyle \mathbf {b} _{i}} b i e i = ( F e i ) e i = F ( e i e i ) = F   . {\displaystyle \mathbf {b} _{i}\otimes \mathbf {e} _{i}=({\boldsymbol {F}}\mathbf {e} _{i})\otimes \mathbf {e} _{i}={\boldsymbol {F}}(\mathbf {e} _{i}\otimes \mathbf {e} _{i})={\boldsymbol {F}}~.} b i = F T e i   ;     g i j = [ F 1 F T ] i j   ;     g i j = [ g i j ] 1 = [ F T F ] i j {\displaystyle \mathbf {b} ^{i}={\boldsymbol {F}}^{-{\rm {T}}}\mathbf {e} ^{i}~;~~g^{ij}=[{\boldsymbol {F}}^{-{\rm {1}}}{\boldsymbol {F}}^{-{\rm {T}}}]_{ij}~;~~g_{ij}=[g^{ij}]^{-1}=[{\boldsymbol {F}}^{\rm {T}}{\boldsymbol {F}}]_{ij}} J := det F {\displaystyle J:=\det {\boldsymbol {F}}} [ b 1 , b 2 , b 3 ] = det F [ e 1 , e 2 , e 3 ]   . {\displaystyle \left[\mathbf {b} _{1},\mathbf {b} _{2},\mathbf {b} _{3}\right]=\det {\boldsymbol {F}}\left[\mathbf {e} _{1},\mathbf {e} _{2},\mathbf {e} _{3}\right]~.} [ e 1 , e 2 , e 3 ] = 1 {\displaystyle \left[\mathbf {e} _{1},\mathbf {e} _{2},\mathbf {e} _{3}\right]=1} J = det F = [ b 1 , b 2 , b 3 ] = b 1 ( b 2 × b 3 ) {\displaystyle J=\det {\boldsymbol {F}}=\left[\mathbf {b} _{1},\mathbf {b} _{2},\mathbf {b} _{3}\right]=\mathbf {b} _{1}\cdot (\mathbf {b} _{2}\times \mathbf {b} _{3})}

まず、 を考えます。 次に 、 同様に、 が示せます。 したがって、 という事実を用いて g := det [ g i j ] {\displaystyle g:=\det[g_{ij}]} g = det [ F T ] det [ F ] = J J = J 2 {\displaystyle g=\det[{\boldsymbol {F}}^{\rm {T}}]\cdot \det[{\boldsymbol {F}}]=J\cdot J=J^{2}} det [ g i j ] = 1 J 2 {\displaystyle \det[g^{ij}]={\cfrac {1}{J^{2}}}} [ g i j ] = [ g i j ] 1 {\displaystyle [g^{ij}]=[g_{ij}]^{-1}} g g i j = 2   J   J g i j = g   g i j {\displaystyle {\cfrac {\partial g}{\partial g_{ij}}}=2~J~{\cfrac {\partial J}{\partial g_{ij}}}=g~g^{ij}}

もう一つの興味深い関係式が以下で導出されます。 定数ですが、まだ定まっていません。 この観察から、関係式が導かれます。 指数表記では、 通常の置換記号です。 b i b j = δ j i b 1 b 1 = 1 ,   b 1 b 2 = b 1 b 3 = 0 b 1 = A   ( b 2 × b 3 ) {\displaystyle \mathbf {b} ^{i}\cdot \mathbf {b} _{j}=\delta _{j}^{i}\quad \Rightarrow \quad \mathbf {b} ^{1}\cdot \mathbf {b} _{1}=1,~\mathbf {b} ^{1}\cdot \mathbf {b} _{2}=\mathbf {b} ^{1}\cdot \mathbf {b} _{3}=0\quad \Rightarrow \quad \mathbf {b} ^{1}=A~(\mathbf {b} _{2}\times \mathbf {b} _{3})} A {\displaystyle A} b 1 b 1 = A   b 1 ( b 2 × b 3 ) = A J = 1 A = 1 J {\displaystyle \mathbf {b} ^{1}\cdot \mathbf {b} _{1}=A~\mathbf {b} _{1}\cdot (\mathbf {b} _{2}\times \mathbf {b} _{3})=AJ=1\quad \Rightarrow \quad A={\cfrac {1}{J}}} b 1 = 1 J ( b 2 × b 3 )   ;     b 2 = 1 J ( b 3 × b 1 )   ;     b 3 = 1 J ( b 1 × b 2 ) {\displaystyle \mathbf {b} ^{1}={\cfrac {1}{J}}(\mathbf {b} _{2}\times \mathbf {b} _{3})~;~~\mathbf {b} ^{2}={\cfrac {1}{J}}(\mathbf {b} _{3}\times \mathbf {b} _{1})~;~~\mathbf {b} ^{3}={\cfrac {1}{J}}(\mathbf {b} _{1}\times \mathbf {b} _{2})} ε i j k   b k = 1 J ( b i × b j ) = 1 g ( b i × b j ) {\displaystyle \varepsilon _{ijk}~\mathbf {b} ^{k}={\cfrac {1}{J}}(\mathbf {b} _{i}\times \mathbf {b} _{j})={\cfrac {1}{\sqrt {g}}}(\mathbf {b} _{i}\times \mathbf {b} _{j})} ε i j k {\displaystyle \varepsilon _{ijk}}

変換テンソルの明示的な表現は示していない。これは、曲線基底と直交基底間の写像の別の形式の方がより有用であるためである。写像が十分な滑らかさを持つと仮定すると(そして多少の表記法の乱用)、次の式が得られる。 同様に、 これらの結果から、次の式が得られる F {\displaystyle {\boldsymbol {F}}} b i = x q i = x x j   x j q i = e j   x j q i {\displaystyle \mathbf {b} _{i}={\cfrac {\partial \mathbf {x} }{\partial q^{i}}}={\cfrac {\partial \mathbf {x} }{\partial x_{j}}}~{\cfrac {\partial x_{j}}{\partial q^{i}}}=\mathbf {e} _{j}~{\cfrac {\partial x_{j}}{\partial q^{i}}}} e i = b j   q j x i {\displaystyle \mathbf {e} _{i}=\mathbf {b} _{j}~{\cfrac {\partial q^{j}}{\partial x_{i}}}} e k b i = x k q i x k q i   b i = e k ( b i b i ) = e k {\displaystyle \mathbf {e} ^{k}\cdot \mathbf {b} _{i}={\frac {\partial x_{k}}{\partial q^{i}}}\quad \Rightarrow \quad {\frac {\partial x_{k}}{\partial q^{i}}}~\mathbf {b} ^{i}=\mathbf {e} ^{k}\cdot (\mathbf {b} _{i}\otimes \mathbf {b} ^{i})=\mathbf {e} ^{k}} b k = q k x i   e i {\displaystyle \mathbf {b} ^{k}={\frac {\partial q^{k}}{\partial x_{i}}}~\mathbf {e} ^{i}}

3次元曲線座標におけるベクトルとテンソルの計算

シモンズ[4]はテンソル解析に関する著書の中でアルバート・アインシュタインの言葉を引用している[7]。

この理論の魔法は、それを真に理解した人なら誰でも感じ取らずにはいられないでしょう。これは、ガウス、リーマン、リッチ、レヴィ=チヴィタによって確立された絶対微分積分法の真の勝利を表しています。

一般曲線座標系におけるベクトルおよびテンソル計算は、一般相対性理論における4次元曲線多様体のテンソル解析、[8] 、曲面シェルの力学[6] 、メタマテリアル[9] [10]で関心を集めているマクスウェル方程式不変性の調査その他多くの分野で使用されています。

この節では、曲線座標系におけるベクトルと2階テンソルの微積分における有用な関係式をいくつか示す。表記法と内容は主にOgden, [2] Simmonds, [4] Green and Zerna, [1] Basar and Weichert, [5] Ciarlet [6]に拠る。

基本的な定義

空間内の点の位置を3つの座標変数で特徴付けるとします ( q 1 , q 2 , q 3 ) {\displaystyle (q^{1},q^{2},q^{3})}

座標曲線は 、とが定数となる曲線を表します点の位置ベクトルをある原点からの相対的な位置ベクトルとします。そして、そのような写像とその逆写像が存在し、それらが連続であると仮定すると、 [2] のように書くことができます。55  これらの体は、曲線座標系曲線座標関数と呼ばれます q 1 {\displaystyle q^{1}} q 2 {\displaystyle q^{2}} q 3 {\displaystyle q^{3}} x {\displaystyle \mathbf {x} } x = φ ( q 1 , q 2 , q 3 )   ;     q i = ψ i ( x ) = [ φ 1 ( x ) ] i {\displaystyle \mathbf {x} ={\boldsymbol {\varphi }}(q^{1},q^{2},q^{3})~;~~q^{i}=\psi ^{i}(\mathbf {x} )=[{\boldsymbol {\varphi }}^{-1}(\mathbf {x} )]^{i}} ψ i ( x ) {\displaystyle \psi ^{i}(\mathbf {x} )} ψ ( x ) = φ 1 ( x ) {\displaystyle {\boldsymbol {\psi }}(\mathbf {x} )={\boldsymbol {\varphi }}^{-1}(\mathbf {x} )}

座標曲線は、 が固定されによって与えられる 1 パラメータ関数族によって定義されます q i {\displaystyle q^{i}} x i ( α ) = φ ( α , q j , q k )   ,     i j k {\displaystyle \mathbf {x} _{i}(\alpha )={\boldsymbol {\varphi }}(\alpha ,q^{j},q^{k})~,~~i\neq j\neq k} q j {\displaystyle q^{j}} q k {\displaystyle q^{k}}

座標曲線の接線ベクトル

点における曲線の接線ベクトル(または点における座標曲線の接線ベクトル)は、 x i {\displaystyle \mathbf {x} _{i}} x i ( α ) {\displaystyle \mathbf {x} _{i}(\alpha )} q i {\displaystyle q_{i}} x {\displaystyle \mathbf {x} } d x i d α x q i {\displaystyle {\cfrac {\rm {{d}\mathbf {x} _{i}}}{\rm {{d}\alpha }}}\equiv {\cfrac {\partial \mathbf {x} }{\partial q^{i}}}}

勾配

スカラー場

を空間上のスカラー場とする。すると、 の勾配はで定義される。 ここで は 任意の定数ベクトルである。成分を と定義すると f ( x ) {\displaystyle f(\mathbf {x} )} f ( x ) = f [ φ ( q 1 , q 2 , q 3 ) ] = f φ ( q 1 , q 2 , q 3 ) {\displaystyle f(\mathbf {x} )=f[{\boldsymbol {\varphi }}(q^{1},q^{2},q^{3})]=f_{\varphi }(q^{1},q^{2},q^{3})} f {\displaystyle f} [ f ( x ) ] c = d d α f ( x + α c ) | α = 0 {\displaystyle [{\boldsymbol {\nabla }}f(\mathbf {x} )]\cdot \mathbf {c} ={\cfrac {\rm {d}}{\rm {{d}\alpha }}}f(\mathbf {x} +\alpha \mathbf {c} ){\biggr |}_{\alpha =0}} c {\displaystyle \mathbf {c} } c i {\displaystyle c^{i}} c {\displaystyle \mathbf {c} } q i + α   c i = ψ i ( x + α   c ) {\displaystyle q^{i}+\alpha ~c^{i}=\psi ^{i}(\mathbf {x} +\alpha ~\mathbf {c} )} [ f ( x ) ] c = d d α f φ ( q 1 + α   c 1 , q 2 + α   c 2 , q 3 + α   c 3 ) | α = 0 = f φ q i   c i = f q i   c i {\displaystyle [{\boldsymbol {\nabla }}f(\mathbf {x} )]\cdot \mathbf {c} ={\cfrac {\rm {d}}{\rm {{d}\alpha }}}f_{\varphi }(q^{1}+\alpha ~c^{1},q^{2}+\alpha ~c^{2},q^{3}+\alpha ~c^{3}){\biggr |}_{\alpha =0}={\cfrac {\partial f_{\varphi }}{\partial q^{i}}}~c^{i}={\cfrac {\partial f}{\partial q^{i}}}~c^{i}}

と設定すると、 なので、 となり、これは ベクトル の反変成分を抽出する手段を提供します f ( x ) = ψ i ( x ) {\displaystyle f(\mathbf {x} )=\psi ^{i}(\mathbf {x} )} q i = ψ i ( x ) {\displaystyle q^{i}=\psi ^{i}(\mathbf {x} )} [ ψ i ( x ) ] c = ψ i q j   c j = c i {\displaystyle [{\boldsymbol {\nabla }}\psi ^{i}(\mathbf {x} )]\cdot \mathbf {c} ={\cfrac {\partial \psi ^{i}}{\partial q^{j}}}~c^{j}=c^{i}} c {\displaystyle \mathbf {c} }

が点における共変(または自然)基底であり、 がその点における反変(または逆数)基底である場合 のセクションでこの基底選択の簡単な根拠を示します。 b i {\displaystyle \mathbf {b} _{i}} b i {\displaystyle \mathbf {b} ^{i}} [ f ( x ) ] c = f q i   c i = ( f q i   b i ) ( c i   b i ) f ( x ) = f q i   b i {\displaystyle [{\boldsymbol {\nabla }}f(\mathbf {x} )]\cdot \mathbf {c} ={\cfrac {\partial f}{\partial q^{i}}}~c^{i}=\left({\cfrac {\partial f}{\partial q^{i}}}~\mathbf {b} ^{i}\right)\left(c^{i}~\mathbf {b} _{i}\right)\quad \Rightarrow \quad {\boldsymbol {\nabla }}f(\mathbf {x} )={\cfrac {\partial f}{\partial q^{i}}}~\mathbf {b} ^{i}}

ベクトル場

同様のプロセスを使用して、ベクトル場 の勾配を求めることができます。勾配は、次のように与えられます。 位置ベクトル場 の勾配を考えると、次が示されます。 ベクトル場 は座標曲線に接し、曲線上の各点で自然な基底を形成します。この記事の冒頭で説明したように、この基底は共変曲線基底とも呼ばれます。逆基底、または反変曲線基底を定義することもできます。テンソル代数のセクションで説明したように、基底ベクトル間のすべての代数関係は、各点 における自然な基底とその逆数に適用されます f ( x ) {\displaystyle \mathbf {f} (\mathbf {x} )} [ f ( x ) ] c = f q i   c i {\displaystyle [{\boldsymbol {\nabla }}\mathbf {f} (\mathbf {x} )]\cdot \mathbf {c} ={\cfrac {\partial \mathbf {f} }{\partial q^{i}}}~c^{i}} r ( x ) = x {\displaystyle \mathbf {r} (\mathbf {x} )=\mathbf {x} } c = x q i   c i = b i ( x )   c i   ;     b i ( x ) := x q i {\displaystyle \mathbf {c} ={\cfrac {\partial \mathbf {x} }{\partial q^{i}}}~c^{i}=\mathbf {b} _{i}(\mathbf {x} )~c^{i}~;~~\mathbf {b} _{i}(\mathbf {x} ):={\cfrac {\partial \mathbf {x} }{\partial q^{i}}}} b i {\displaystyle \mathbf {b} _{i}} q i {\displaystyle q^{i}} b i {\displaystyle \mathbf {b} ^{i}} x {\displaystyle \mathbf {x} }

は任意なので、次のように書くことができる。 c {\displaystyle \mathbf {c} } f ( x ) = f q i b i {\displaystyle {\boldsymbol {\nabla }}\mathbf {f} (\mathbf {x} )={\cfrac {\partial \mathbf {f} }{\partial q^{i}}}\otimes \mathbf {b} ^{i}}

反変基底ベクトルは定数の面に垂直であり、次のように与えられる ことに注意する。 b i {\displaystyle \mathbf {b} ^{i}} ψ i {\displaystyle \psi ^{i}} b i = ψ i {\displaystyle \mathbf {b} ^{i}={\boldsymbol {\nabla }}\psi ^{i}}

第一種クリストッフェル記号

第一種クリストッフェル記号は次のように定義される。 を 用い表す 、 となる。 となる ので、 となる。これらを用いて上記の関係式を整理すると、 b i , j = b i q j := Γ i j k   b k b i , j b l = Γ i j l {\displaystyle \mathbf {b} _{i,j}={\frac {\partial \mathbf {b} _{i}}{\partial q^{j}}}:=\Gamma _{ijk}~\mathbf {b} ^{k}\quad \Rightarrow \quad \mathbf {b} _{i,j}\cdot \mathbf {b} _{l}=\Gamma _{ijl}} Γ i j k {\displaystyle \Gamma _{ijk}} g i j {\displaystyle g_{ij}} g i j , k = ( b i b j ) , k = b i , k b j + b i b j , k = Γ i k j + Γ j k i g i k , j = ( b i b k ) , j = b i , j b k + b i b k , j = Γ i j k + Γ k j i g j k , i = ( b j b k ) , i = b j , i b k + b j b k , i = Γ j i k + Γ k i j {\displaystyle {\begin{aligned}g_{ij,k}&=(\mathbf {b} _{i}\cdot \mathbf {b} _{j})_{,k}=\mathbf {b} _{i,k}\cdot \mathbf {b} _{j}+\mathbf {b} _{i}\cdot \mathbf {b} _{j,k}=\Gamma _{ikj}+\Gamma _{jki}\\g_{ik,j}&=(\mathbf {b} _{i}\cdot \mathbf {b} _{k})_{,j}=\mathbf {b} _{i,j}\cdot \mathbf {b} _{k}+\mathbf {b} _{i}\cdot \mathbf {b} _{k,j}=\Gamma _{ijk}+\Gamma _{kji}\\g_{jk,i}&=(\mathbf {b} _{j}\cdot \mathbf {b} _{k})_{,i}=\mathbf {b} _{j,i}\cdot \mathbf {b} _{k}+\mathbf {b} _{j}\cdot \mathbf {b} _{k,i}=\Gamma _{jik}+\Gamma _{kij}\end{aligned}}} b i , j = b j , i {\displaystyle \mathbf {b} _{i,j}=\mathbf {b} _{j,i}} Γ i j k = Γ j i k {\displaystyle \Gamma _{ijk}=\Gamma _{jik}} Γ i j k = 1 2 ( g i k , j + g j k , i g i j , k ) = 1 2 [ ( b i b k ) , j + ( b j b k ) , i ( b i b j ) , k ] {\displaystyle \Gamma _{ijk}={\frac {1}{2}}(g_{ik,j}+g_{jk,i}-g_{ij,k})={\frac {1}{2}}[(\mathbf {b} _{i}\cdot \mathbf {b} _{k})_{,j}+(\mathbf {b} _{j}\cdot \mathbf {b} _{k})_{,i}-(\mathbf {b} _{i}\cdot \mathbf {b} _{j})_{,k}]}

第二種クリストッフェル記号

第二種クリストッフェル記号は次のように定義 れる Γ i j k = Γ j i k {\displaystyle \Gamma _{ij}^{k}=\Gamma _{ji}^{k}}

b i q j = Γ i j k   b k {\displaystyle {\cfrac {\partial \mathbf {b} _{i}}{\partial q^{j}}}=\Gamma _{ij}^{k}~\mathbf {b} _{k}}

これは、 以下の他の関係が Γ i j k = b i q j b k = b i b k q j {\displaystyle \Gamma _{ij}^{k}={\cfrac {\partial \mathbf {b} _{i}}{\partial q^{j}}}\cdot \mathbf {b} ^{k}=-\mathbf {b} _{i}\cdot {\cfrac {\partial \mathbf {b} ^{k}}{\partial q^{j}}}} b i q j = Γ j k i   b k   ;     b i = Γ i j k   b k b j   ;     b i = Γ j k i   b k b j {\displaystyle {\cfrac {\partial \mathbf {b} ^{i}}{\partial q^{j}}}=-\Gamma _{jk}^{i}~\mathbf {b} ^{k}~;~~{\boldsymbol {\nabla }}\mathbf {b} _{i}=\Gamma _{ij}^{k}~\mathbf {b} _{k}\otimes \mathbf {b} ^{j}~;~~{\boldsymbol {\nabla }}\mathbf {b} ^{i}=-\Gamma _{jk}^{i}~\mathbf {b} ^{k}\otimes \mathbf {b} ^{j}}

クリストッフェル記号が計量テンソルとその導関数のみに依存することを示すもう一つの特に有用な関係は、 Γ i j k = g k m 2 ( g m i q j + g m j q i g i j q m ) {\displaystyle \Gamma _{ij}^{k}={\frac {g^{km}}{2}}\left({\frac {\partial g_{mi}}{\partial q^{j}}}+{\frac {\partial g_{mj}}{\partial q^{i}}}-{\frac {\partial g_{ij}}{\partial q^{m}}}\right)}

ベクトル場の勾配の明示的な表現

曲線座標におけるベクトル場の勾配を表す次の式は非常に便利です。 v = [ v i q k + Γ l k i   v l ]   b i b k = [ v i q k Γ k i l   v l ]   b i b k {\displaystyle {\begin{aligned}{\boldsymbol {\nabla }}\mathbf {v} &=\left[{\cfrac {\partial v^{i}}{\partial q^{k}}}+\Gamma _{lk}^{i}~v^{l}\right]~\mathbf {b} _{i}\otimes \mathbf {b} ^{k}\\[8pt]&=\left[{\cfrac {\partial v_{i}}{\partial q^{k}}}-\Gamma _{ki}^{l}~v_{l}\right]~\mathbf {b} ^{i}\otimes \mathbf {b} ^{k}\end{aligned}}}

物理的なベクトル場を表現する

ベクトル場は、場の共変成分、物理的な成分、(合計なし) は正規化された反変基底ベクトル として表すことができます 。 v {\displaystyle \mathbf {v} } v = v i   b i = v ^ i   b ^ i {\displaystyle \mathbf {v} =v_{i}~\mathbf {b} ^{i}={\hat {v}}_{i}~{\hat {\mathbf {b} }}^{i}} v i {\displaystyle v_{i}} v ^ i {\displaystyle {\hat {v}}_{i}} b ^ i = b i g i i {\displaystyle {\hat {\mathbf {b} }}^{i}={\cfrac {\mathbf {b} ^{i}}{\sqrt {g^{ii}}}}}

2階テンソル場

2階テンソル場の勾配も同様に次のように表される。 S = S q i b i {\displaystyle {\boldsymbol {\nabla }}{\boldsymbol {S}}={\frac {\partial {\boldsymbol {S}}}{\partial q^{i}}}\otimes \mathbf {b} ^{i}}

勾配の明示的な表現

テンソルの表現を反変基底で考えると、次 のようにも書ける。 S = q k [ S i j   b i b j ] b k = [ S i j q k Γ k i l   S l j Γ k j l   S i l ]   b i b j b k {\displaystyle {\boldsymbol {\nabla }}{\boldsymbol {S}}={\frac {\partial }{\partial q^{k}}}[S_{ij}~\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}]\otimes \mathbf {b} ^{k}=\left[{\frac {\partial S_{ij}}{\partial q^{k}}}-\Gamma _{ki}^{l}~S_{lj}-\Gamma _{kj}^{l}~S_{il}\right]~\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}\otimes \mathbf {b} ^{k}} S = [ S i j q k + Γ k l i   S l j + Γ k l j   S i l ]   b i b j b k = [ S   j i q k + Γ k l i   S   j l Γ k j l   S   l i ]   b i b j b k = [ S i   j q k Γ i k l   S l   j + Γ k l j   S i   l ]   b i b j b k {\displaystyle {\begin{aligned}{\boldsymbol {\nabla }}{\boldsymbol {S}}&=\left[{\cfrac {\partial S^{ij}}{\partial q^{k}}}+\Gamma _{kl}^{i}~S^{lj}+\Gamma _{kl}^{j}~S^{il}\right]~\mathbf {b} _{i}\otimes \mathbf {b} _{j}\otimes \mathbf {b} ^{k}\\[8pt]&=\left[{\cfrac {\partial S_{~j}^{i}}{\partial q^{k}}}+\Gamma _{kl}^{i}~S_{~j}^{l}-\Gamma _{kj}^{l}~S_{~l}^{i}\right]~\mathbf {b} _{i}\otimes \mathbf {b} ^{j}\otimes \mathbf {b} ^{k}\\[8pt]&=\left[{\cfrac {\partial S_{i}^{~j}}{\partial q^{k}}}-\Gamma _{ik}^{l}~S_{l}^{~j}+\Gamma _{kl}^{j}~S_{i}^{~l}\right]~\mathbf {b} ^{i}\otimes \mathbf {b} _{j}\otimes \mathbf {b} ^{k}\end{aligned}}}

物理的な2階テンソル場を表現する

2階テンソル場の物理的成分は、正規化された反変基底、すなわち ハット基底ベクトルが正規化された基底を用いることで得られる。これは(これもまた和は不要である) S = S i j   b i b j = S ^ i j   b ^ i b ^ j {\displaystyle {\boldsymbol {S}}=S_{ij}~\mathbf {b} ^{i}\otimes \mathbf {b} ^{j}={\hat {S}}_{ij}~{\hat {\mathbf {b} }}^{i}\otimes {\hat {\mathbf {b} }}^{j}}

S ^ i j = S i j   g i i   g j j {\displaystyle {\hat {S}}_{ij}=S_{ij}~{\sqrt {g^{ii}~g^{jj}}}}

発散

ベクトル場

ベクトル場の発散 、曲線基底に関する成分に関して 次のように定義される。 v {\displaystyle \mathbf {v} } div   v = v = tr ( v ) {\displaystyle \operatorname {div} ~\mathbf {v} ={\boldsymbol {\nabla }}\cdot \mathbf {v} ={\text{tr}}({\boldsymbol {\nabla }}\mathbf {v} )} v = v i q i + Γ i i   v = [ v i q j Γ j i   v ]   g i j {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\cfrac {\partial v^{i}}{\partial q^{i}}}+\Gamma _{\ell i}^{i}~v^{\ell }=\left[{\cfrac {\partial v_{i}}{\partial q^{j}}}-\Gamma _{ji}^{\ell }~v_{\ell }\right]~g^{ij}}

ベクトル場の発散を表す別の方程式が頻繁に用いられます。この関係式を導くには、次の関係を思い出してください。 ここで、 の対称性により、次の関係が成り立ち ます が を成分とする行列である 場合、その逆行列は となります。逆行列は で与えられます。 ここで、はを成分とする余因子行列です。行列代数から、 次の関係が成り立ちます。 したがって、 この関係式を発散の式に代入すると、次の式が得られます。 少し操作を加えると、よりコンパクトな形式になります。 v = v i q i + Γ i i   v {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\frac {\partial v^{i}}{\partial q^{i}}}+\Gamma _{\ell i}^{i}~v^{\ell }} Γ i i = Γ i i = g m i 2 [ g i m q + g m q i g i l q m ] {\displaystyle \Gamma _{\ell i}^{i}=\Gamma _{i\ell }^{i}={\cfrac {g^{mi}}{2}}\left[{\frac {\partial g_{im}}{\partial q^{\ell }}}+{\frac {\partial g_{\ell m}}{\partial q^{i}}}-{\frac {\partial g_{il}}{\partial q^{m}}}\right]} g {\displaystyle {\boldsymbol {g}}} g m i   g m q i = g m i   g i q m {\displaystyle g^{mi}~{\frac {\partial g_{\ell m}}{\partial q^{i}}}=g^{mi}~{\frac {\partial g_{i\ell }}{\partial q^{m}}}} v = v i q i + g m i 2   g i m q   v {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\frac {\partial v^{i}}{\partial q^{i}}}+{\cfrac {g^{mi}}{2}}~{\frac {\partial g_{im}}{\partial q^{\ell }}}~v^{\ell }} [ g i j ] {\displaystyle [g_{ij}]} g i j {\displaystyle g_{ij}} [ g i j ] 1 = [ g i j ] {\displaystyle [g_{ij}]^{-1}=[g^{ij}]} [ g i j ] = [ g i j ] 1 = A i j g   ;     g := det ( [ g i j ] ) = det g {\displaystyle [g^{ij}]=[g_{ij}]^{-1}={\cfrac {A^{ij}}{g}}~;~~g:=\det([g_{ij}])=\det {\boldsymbol {g}}} A i j {\displaystyle A^{ij}} g i j {\displaystyle g_{ij}} g = det ( [ g i j ] ) = i g i j   A i j g g i j = A i j {\displaystyle g=\det([g_{ij}])=\sum _{i}g_{ij}~A^{ij}\quad \Rightarrow \quad {\frac {\partial g}{\partial g_{ij}}}=A^{ij}} [ g i j ] = 1 g   g g i j {\displaystyle [g^{ij}]={\cfrac {1}{g}}~{\frac {\partial g}{\partial g_{ij}}}} v = v i q i + 1 2 g   g g m i   g i m q   v = v i q i + 1 2 g   g q   v {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\frac {\partial v^{i}}{\partial q^{i}}}+{\cfrac {1}{2g}}~{\frac {\partial g}{\partial g_{mi}}}~{\frac {\partial g_{im}}{\partial q^{\ell }}}~v^{\ell }={\frac {\partial v^{i}}{\partial q^{i}}}+{\cfrac {1}{2g}}~{\frac {\partial g}{\partial q^{\ell }}}~v^{\ell }} v = 1 g   q i ( v i   g ) {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\cfrac {1}{\sqrt {g}}}~{\frac {\partial }{\partial q^{i}}}(v^{i}~{\sqrt {g}})}

2階テンソル場

2階テンソル場の発散は、任意の定数ベクトルを用いて定義される [ 11 ]曲線 座標では、 ( S ) a = ( S a ) {\displaystyle ({\boldsymbol {\nabla }}\cdot {\boldsymbol {S}})\cdot \mathbf {a} ={\boldsymbol {\nabla }}\cdot ({\boldsymbol {S}}\mathbf {a} )} a {\displaystyle \mathbf {a} } S = [ S i j q k Γ k i l   S l j Γ k j l   S i l ]   g i k   b j = [ S i j q i + Γ i l i   S l j + Γ i l j   S i l ]   b j = [ S   j i q i + Γ i l i   S   j l Γ i j l   S   l i ]   b j = [ S i   j q k Γ i k l   S l   j + Γ k l j   S i   l ]   g i k   b j {\displaystyle {\begin{aligned}{\boldsymbol {\nabla }}\cdot {\boldsymbol {S}}&=\left[{\cfrac {\partial S_{ij}}{\partial q^{k}}}-\Gamma _{ki}^{l}~S_{lj}-\Gamma _{kj}^{l}~S_{il}\right]~g^{ik}~\mathbf {b} ^{j}\\[8pt]&=\left[{\cfrac {\partial S^{ij}}{\partial q^{i}}}+\Gamma _{il}^{i}~S^{lj}+\Gamma _{il}^{j}~S^{il}\right]~\mathbf {b} _{j}\\[8pt]&=\left[{\cfrac {\partial S_{~j}^{i}}{\partial q^{i}}}+\Gamma _{il}^{i}~S_{~j}^{l}-\Gamma _{ij}^{l}~S_{~l}^{i}\right]~\mathbf {b} ^{j}\\[8pt]&=\left[{\cfrac {\partial S_{i}^{~j}}{\partial q^{k}}}-\Gamma _{ik}^{l}~S_{l}^{~j}+\Gamma _{kl}^{j}~S_{i}^{~l}\right]~g^{ik}~\mathbf {b} _{j}\end{aligned}}}

ラプラシアン

スカラー場

スカラー場のラプラシアンは 次の ように定義されます。 ベクトル場の発散の別の表現を使うと次のようになります 。したがって、 φ ( x ) {\displaystyle \varphi (\mathbf {x} )} 2 φ := ( φ ) {\displaystyle \nabla ^{2}\varphi :={\boldsymbol {\nabla }}\cdot ({\boldsymbol {\nabla }}\varphi )} 2 φ = 1 g   q i ( [ φ ] i   g ) {\displaystyle \nabla ^{2}\varphi ={\cfrac {1}{\sqrt {g}}}~{\frac {\partial }{\partial q^{i}}}([{\boldsymbol {\nabla }}\varphi ]^{i}~{\sqrt {g}})} φ = φ q l   b l = g l i   φ q l   b i [ φ ] i = g l i   φ q l {\displaystyle {\boldsymbol {\nabla }}\varphi ={\frac {\partial \varphi }{\partial q^{l}}}~\mathbf {b} ^{l}=g^{li}~{\frac {\partial \varphi }{\partial q^{l}}}~\mathbf {b} _{i}\quad \Rightarrow \quad [{\boldsymbol {\nabla }}\varphi ]^{i}=g^{li}~{\frac {\partial \varphi }{\partial q^{l}}}} 2 φ = 1 g   q i ( g l i   φ q l   g ) {\displaystyle \nabla ^{2}\varphi ={\cfrac {1}{\sqrt {g}}}~{\frac {\partial }{\partial q^{i}}}\left(g^{li}~{\frac {\partial \varphi }{\partial q^{l}}}~{\sqrt {g}}\right)}

ベクトル場の回転

共変曲線座標におけるベクトル場の回転はのように表される。 v {\displaystyle \mathbf {v} } × v = E r s t v s | r   b t {\displaystyle {\boldsymbol {\nabla }}\times \mathbf {v} ={\mathcal {E}}^{rst}v_{s|r}~\mathbf {b} _{t}} v s | r = v s , r Γ s r i   v i {\displaystyle v_{s|r}=v_{s,r}-\Gamma _{sr}^{i}~v_{i}}

直交曲線座標

この節では、曲線座標系が直交、すなわち 、 あるいは同値で であると仮定する。前述と同様に、は共変基底ベクトルでありは反変基底ベクトルである。また、 は背景の固定直交基底とする。直交曲線座標系の一覧は以下の通りである。 b i b j = { g i i if  i = j 0 if  i j , {\displaystyle \mathbf {b} _{i}\cdot \mathbf {b} _{j}={\begin{cases}g_{ii}&{\text{if }}i=j\\0&{\text{if }}i\neq j,\end{cases}}} b i b j = { g i i if  i = j 0 if  i j , {\displaystyle \mathbf {b} ^{i}\cdot \mathbf {b} ^{j}={\begin{cases}g^{ii}&{\text{if }}i=j\\0&{\text{if }}i\neq j,\end{cases}}} g i i = g i i 1 {\displaystyle g^{ii}=g_{ii}^{-1}} b i , b j {\displaystyle \mathbf {b} _{i},\mathbf {b} _{j}} b i {\displaystyle \mathbf {b} ^{i}} b j {\displaystyle \mathbf {b} ^{j}} ( e 1 , e 2 , e 3 ) {\displaystyle (\mathbf {e} ^{1},\mathbf {e} ^{2},\mathbf {e} ^{3})}

直交曲線座標における計量テンソル

座標系の原点に対する点の位置ベクトルをと する。表記は=とすることで簡略化できる。各点において、小さな線分要素 を作図できる。線分要素の長さの2乗はスカラー積であり、空間計量と呼ばれる。曲線座標系について話す場合、対象となる空間はユークリッド空間であると仮定することを思い出してほしい。位置ベクトルを、背景となる固定された直交座標基底、すなわち を用いて表すとしよう。 r ( x ) {\displaystyle \mathbf {r} (\mathbf {x} )} x {\displaystyle \mathbf {x} } x {\displaystyle \mathbf {x} } r ( x ) {\displaystyle \mathbf {r} (\mathbf {x} )} d x {\displaystyle \mathrm {d} \mathbf {x} } d x d x {\displaystyle \mathrm {d} \mathbf {x} \cdot \mathrm {d} \mathbf {x} } x = i = 1 3 x i   e i {\displaystyle \mathbf {x} =\sum _{i=1}^{3}x_{i}~\mathbf {e} _{i}}

連鎖律を用いると、 3次元直交曲線座標系で次のように表すことができます。したがって 、計量は次のように与えられます。 d x {\displaystyle \mathrm {d} \mathbf {x} } ( q 1 , q 2 , q 3 ) {\displaystyle (q^{1},q^{2},q^{3})} d x = i = 1 3 j = 1 3 ( x i q j   e i ) d q j {\displaystyle \mathrm {d} \mathbf {x} =\sum _{i=1}^{3}\sum _{j=1}^{3}\left({\cfrac {\partial x_{i}}{\partial q^{j}}}~\mathbf {e} _{i}\right)\mathrm {d} q^{j}} d x d x = i = 1 3 j = 1 3 k = 1 3 x i q j   x i q k   d q j   d q k {\displaystyle \mathrm {d} \mathbf {x} \cdot \mathrm {d} \mathbf {x} =\sum _{i=1}^{3}\sum _{j=1}^{3}\sum _{k=1}^{3}{\cfrac {\partial x_{i}}{\partial q^{j}}}~{\cfrac {\partial x_{i}}{\partial q^{k}}}~\mathrm {d} q^{j}~\mathrm {d} q^{k}}

対称量は 、曲線座標における ユークリッド空間基本テンソル(または計量テンソル) と呼ばれます。 g i j ( q i , q j ) = k = 1 3 x k q i   x k q j = b i b j {\displaystyle g_{ij}(q^{i},q^{j})=\sum _{k=1}^{3}{\cfrac {\partial x_{k}}{\partial q^{i}}}~{\cfrac {\partial x_{k}}{\partial q^{j}}}=\mathbf {b} _{i}\cdot \mathbf {b} _{j}}

また、 Lamé 係数が どこにあるかにも注意してください。 g i j = x q i x q j = ( k h k i   e k ) ( m h m j   e m ) = k h k i   h k j {\displaystyle g_{ij}={\cfrac {\partial \mathbf {x} }{\partial q^{i}}}\cdot {\cfrac {\partial \mathbf {x} }{\partial q^{j}}}=\left(\sum _{k}h_{ki}~\mathbf {e} _{k}\right)\cdot \left(\sum _{m}h_{mj}~\mathbf {e} _{m}\right)=\sum _{k}h_{ki}~h_{kj}} h i j {\displaystyle h_{ij}}

スケール係数 を定義すると、基本 テンソルとラメ係数の関係が得られます。 h i {\displaystyle h_{i}} b i b i = g i i = k h k i 2 =: h i 2 | x q i | = | b i | = g i i = h i {\displaystyle \mathbf {b} _{i}\cdot \mathbf {b} _{i}=g_{ii}=\sum _{k}h_{ki}^{2}=:h_{i}^{2}\quad \Rightarrow \quad \left|{\cfrac {\partial \mathbf {x} }{\partial q^{i}}}\right|=\left|\mathbf {b} _{i}\right|={\sqrt {g_{ii}}}=h_{i}}

例: 極座標

の極座標を考えるときは曲線座標であり、変換のヤコビ行列式は であることに注意してください R 2 {\displaystyle \mathbb {R} ^{2}} ( x , y ) = ( r cos θ , r sin θ ) {\displaystyle (x,y)=(r\cos \theta ,r\sin \theta )} ( r , θ ) {\displaystyle (r,\theta )} ( r , θ ) ( r cos θ , r sin θ ) {\displaystyle (r,\theta )\to (r\cos \theta ,r\sin \theta )} r {\displaystyle r}

直交基底ベクトルは、です。正規化され基底ベクトルはであり、スケール係数は 、 です基礎テンソルは、です b r = ( cos θ , sin θ ) {\displaystyle \mathbf {b} _{r}=(\cos \theta ,\sin \theta )} b θ = ( r sin θ , r cos θ ) {\displaystyle \mathbf {b} _{\theta }=(-r\sin \theta ,r\cos \theta )} e r = ( cos θ , sin θ ) {\displaystyle \mathbf {e} _{r}=(\cos \theta ,\sin \theta )} e θ = ( sin θ , cos θ ) {\displaystyle \mathbf {e} _{\theta }=(-\sin \theta ,\cos \theta )} h r = 1 {\displaystyle h_{r}=1} h θ = r {\displaystyle h_{\theta }=r} g 11 = 1 {\displaystyle g_{11}=1} g 22 = r 2 {\displaystyle g_{22}=r^{2}} g 12 = g 21 = 0 {\displaystyle g_{12}=g_{21}=0}

線積分と面積分

ベクトル解析の計算に曲線座標系を用いる場合、線積分、面積分、体積積分の計算において調整が必要です。簡潔にするため、ここでも議論は3次元と直交曲線座標系に限定します。ただし、座標系が直交でない場合は式にいくつかの追加項が加わりますが、同じ議論は3次元問題にも当てはまります。 n {\displaystyle n}

線積分

通常、線積分の計算では、直交座標系における が媒介変数化される 場所を計算することに関心がある 。曲線座標系では、 C f d s = a b f ( x ( t ) ) | x t | d t {\displaystyle \int _{C}f\,ds=\int _{a}^{b}f(\mathbf {x} (t))\left|{\partial \mathbf {x} \over \partial t}\right|\;dt} x ( t ) {\displaystyle \mathbf {x} (t)} C {\displaystyle C}

| x t | = | i = 1 3 x q i q i t | {\displaystyle \left|{\partial \mathbf {x} \over \partial t}\right|=\left|\sum _{i=1}^{3}{\partial \mathbf {x} \over \partial q^{i}}{\partial q^{i} \over \partial t}\right|}

連鎖律により。そしてラメ係数の定義から、

x q i = k h k i   e k {\displaystyle {\partial \mathbf {x} \over \partial q^{i}}=\sum _{k}h_{ki}~\mathbf {e} _{k}}

そしてこうして

| x t | = | k ( i h k i   q i t ) e k | = i j k h k i   h k j q i t q j t = i j g i j   q i t q j t {\displaystyle {\begin{aligned}\left|{\partial \mathbf {x} \over \partial t}\right|&=\left|\sum _{k}\left(\sum _{i}h_{ki}~{\cfrac {\partial q^{i}}{\partial t}}\right)\mathbf {e} _{k}\right|\\[8pt]&={\sqrt {\sum _{i}\sum _{j}\sum _{k}h_{ki}~h_{kj}{\cfrac {\partial q^{i}}{\partial t}}{\cfrac {\partial q^{j}}{\partial t}}}}={\sqrt {\sum _{i}\sum _{j}g_{ij}~{\cfrac {\partial q^{i}}{\partial t}}{\cfrac {\partial q^{j}}{\partial t}}}}\end{aligned}}}

さて、いつから、私たちは 、そして正常に進めることができるようになりました。 g i j = 0 {\displaystyle g_{ij}=0} i j {\displaystyle i\neq j} | x t | = i g i i   ( q i t ) 2 = i h i 2   ( q i t ) 2 {\displaystyle \left|{\partial \mathbf {x} \over \partial t}\right|={\sqrt {\sum _{i}g_{ii}~\left({\cfrac {\partial q^{i}}{\partial t}}\right)^{2}}}={\sqrt {\sum _{i}h_{i}^{2}~\left({\cfrac {\partial q^{i}}{\partial t}}\right)^{2}}}}

面積分

同様に、面積分に興味がある場合、直交座標系における表面のパラメータ化を用いた関連する計算は次のようになる。 また、曲線座標系では、 となり 、曲線座標系の定義を再び利用して、 を得る。 S f d S = T f ( x ( s , t ) ) | x s × x t | d s d t {\displaystyle \int _{S}f\,dS=\iint _{T}f(\mathbf {x} (s,t))\left|{\partial \mathbf {x} \over \partial s}\times {\partial \mathbf {x} \over \partial t}\right|\,ds\,dt} | x s × x t | = | ( i x q i q i s ) × ( j x q j q j t ) | {\displaystyle \left|{\partial \mathbf {x} \over \partial s}\times {\partial \mathbf {x} \over \partial t}\right|=\left|\left(\sum _{i}{\partial \mathbf {x} \over \partial q^{i}}{\partial q^{i} \over \partial s}\right)\times \left(\sum _{j}{\partial \mathbf {x} \over \partial q^{j}}{\partial q^{j} \over \partial t}\right)\right|} x q i q i s = k ( i = 1 3 h k i   q i s ) e k   ;     x q j q j t = m ( j = 1 3 h m j   q j t ) e m {\displaystyle {\partial \mathbf {x} \over \partial q^{i}}{\partial q^{i} \over \partial s}=\sum _{k}\left(\sum _{i=1}^{3}h_{ki}~{\partial q^{i} \over \partial s}\right)\mathbf {e} _{k}~;~~{\partial \mathbf {x} \over \partial q^{j}}{\partial q^{j} \over \partial t}=\sum _{m}\left(\sum _{j=1}^{3}h_{mj}~{\partial q^{j} \over \partial t}\right)\mathbf {e} _{m}}

したがって、 順列記号です | x s × x t | = | k m ( i = 1 3 h k i   q i s ) ( j = 1 3 h m j   q j t ) e k × e m | = | p k m E k m p ( i = 1 3 h k i   q i s ) ( j = 1 3 h m j   q j t ) e p | {\displaystyle {\begin{aligned}\left|{\partial \mathbf {x} \over \partial s}\times {\partial \mathbf {x} \over \partial t}\right|&=\left|\sum _{k}\sum _{m}\left(\sum _{i=1}^{3}h_{ki}~{\partial q^{i} \over \partial s}\right)\left(\sum _{j=1}^{3}h_{mj}~{\partial q^{j} \over \partial t}\right)\mathbf {e} _{k}\times \mathbf {e} _{m}\right|\\[8pt]&=\left|\sum _{p}\sum _{k}\sum _{m}{\mathcal {E}}_{kmp}\left(\sum _{i=1}^{3}h_{ki}~{\partial q^{i} \over \partial s}\right)\left(\sum _{j=1}^{3}h_{mj}~{\partial q^{j} \over \partial t}\right)\mathbf {e} _{p}\right|\end{aligned}}} E {\displaystyle {\mathcal {E}}}

行列式の形式では、曲線座標に関する外積は次のようになります。 | e 1 e 2 e 3 i h 1 i q i s i h 2 i q i s i h 3 i q i s j h 1 j q j t j h 2 j q j t j h 3 j q j t | {\displaystyle {\begin{vmatrix}\mathbf {e} _{1}&\mathbf {e} _{2}&\mathbf {e} _{3}\\&&\\\sum _{i}h_{1i}{\partial q^{i} \over \partial s}&\sum _{i}h_{2i}{\partial q^{i} \over \partial s}&\sum _{i}h_{3i}{\partial q^{i} \over \partial s}\\&&\\\sum _{j}h_{1j}{\partial q^{j} \over \partial t}&\sum _{j}h_{2j}{\partial q^{j} \over \partial t}&\sum _{j}h_{3j}{\partial q^{j} \over \partial t}\end{vmatrix}}}

勾配、回転、除算、ラプラシアン

3 次元の 直交曲線座標では、スカラー場またはベクトル場勾配を次のように 表すことができます。 直交基底の場合、 ベクトル場の発散は次のように記述できます。 また したがって 、 次の式に注意することで、同様の方法で ラプラシアン の式を取得できます。 次に、次の式が得られます。 勾配、発散、およびラプラシアンの式は、直接-次元に拡張できます。 b i = k g i k   b k   ;     g i i = 1 g i i = 1 h i 2 {\displaystyle \mathbf {b} ^{i}=\sum _{k}g^{ik}~\mathbf {b} _{k}~;~~g^{ii}={\cfrac {1}{g_{ii}}}={\cfrac {1}{h_{i}^{2}}}} φ = i φ q i   b i = i j φ q i   g i j   b j = i 1 h i 2   f q i   b i   ;     v = i 1 h i 2   v q i b i {\displaystyle \nabla \varphi =\sum _{i}{\partial \varphi \over \partial q^{i}}~\mathbf {b} ^{i}=\sum _{i}\sum _{j}{\partial \varphi \over \partial q^{i}}~g^{ij}~\mathbf {b} _{j}=\sum _{i}{\cfrac {1}{h_{i}^{2}}}~{\partial f \over \partial q^{i}}~\mathbf {b} _{i}~;~~\nabla \mathbf {v} =\sum _{i}{\cfrac {1}{h_{i}^{2}}}~{\partial \mathbf {v} \over \partial q^{i}}\otimes \mathbf {b} _{i}} g = g 11   g 22   g 33 = h 1 2   h 2 2   h 3 2 g = h 1 h 2 h 3 {\displaystyle g=g_{11}~g_{22}~g_{33}=h_{1}^{2}~h_{2}^{2}~h_{3}^{2}\quad \Rightarrow \quad {\sqrt {g}}=h_{1}h_{2}h_{3}} v = 1 h 1 h 2 h 3   q i ( h 1 h 2 h 3   v i ) {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\cfrac {1}{h_{1}h_{2}h_{3}}}~{\frac {\partial }{\partial q^{i}}}(h_{1}h_{2}h_{3}~v^{i})} v i = g i k   v k v 1 = g 11   v 1 = v 1 h 1 2   ;     v 2 = g 22   v 2 = v 2 h 2 2   ;     v 3 = g 33   v 3 = v 3 h 3 2 {\displaystyle v^{i}=g^{ik}~v_{k}\quad \Rightarrow v^{1}=g^{11}~v_{1}={\cfrac {v_{1}}{h_{1}^{2}}}~;~~v^{2}=g^{22}~v_{2}={\cfrac {v_{2}}{h_{2}^{2}}}~;~~v^{3}=g^{33}~v_{3}={\cfrac {v_{3}}{h_{3}^{2}}}} v = 1 h 1 h 2 h 3   i q i ( h 1 h 2 h 3 h i 2   v i ) {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\cfrac {1}{h_{1}h_{2}h_{3}}}~\sum _{i}{\frac {\partial }{\partial q^{i}}}\left({\cfrac {h_{1}h_{2}h_{3}}{h_{i}^{2}}}~v_{i}\right)} g l i   φ q l = { g 11   φ q 1 , g 22   φ q 2 , g 33   φ q 3 } = { 1 h 1 2   φ q 1 , 1 h 2 2   φ q 2 , 1 h 3 2   φ q 3 } {\displaystyle g^{li}~{\frac {\partial \varphi }{\partial q^{l}}}=\left\{g^{11}~{\frac {\partial \varphi }{\partial q^{1}}},g^{22}~{\frac {\partial \varphi }{\partial q^{2}}},g^{33}~{\frac {\partial \varphi }{\partial q^{3}}}\right\}=\left\{{\cfrac {1}{h_{1}^{2}}}~{\frac {\partial \varphi }{\partial q^{1}}},{\cfrac {1}{h_{2}^{2}}}~{\frac {\partial \varphi }{\partial q^{2}}},{\cfrac {1}{h_{3}^{2}}}~{\frac {\partial \varphi }{\partial q^{3}}}\right\}} 2 φ = 1 h 1 h 2 h 3   i q i ( h 1 h 2 h 3 h i 2   φ q i ) {\displaystyle \nabla ^{2}\varphi ={\cfrac {1}{h_{1}h_{2}h_{3}}}~\sum _{i}{\frac {\partial }{\partial q^{i}}}\left({\cfrac {h_{1}h_{2}h_{3}}{h_{i}^{2}}}~{\frac {\partial \varphi }{\partial q^{i}}}\right)} n {\displaystyle n}

ベクトル場回転は で与えられ、 レヴィ・チヴィタ記号です × v = 1 h 1 h 2 h 3 i = 1 n e i j k ε i j k h i ( h k v k ) q j {\displaystyle \nabla \times \mathbf {v} ={\frac {1}{h_{1}h_{2}h_{3}}}\sum _{i=1}^{n}\mathbf {e} _{i}\sum _{jk}\varepsilon _{ijk}h_{i}{\frac {\partial (h_{k}v_{k})}{\partial q^{j}}}} ε i j k {\displaystyle \varepsilon _{ijk}}

例: 円筒極座標

円筒 座標 ( x 1 , x 2 , x 3 ) = x = φ ( q 1 , q 2 , q 3 ) = φ ( r , θ , z ) = { r cos θ , r sin θ , z } {\displaystyle (x_{1},x_{2},x_{3})=\mathbf {x} ={\boldsymbol {\varphi }}(q^{1},q^{2},q^{3})={\boldsymbol {\varphi }}(r,\theta ,z)=\{r\cos \theta ,r\sin \theta ,z\}} { ψ 1 ( x ) , ψ 2 ( x ) , ψ 3 ( x ) } = ( q 1 , q 2 , q 3 ) ( r , θ , z ) = { x 1 2 + x 2 2 , tan 1 ( x 2 / x 1 ) , x 3 } {\displaystyle \{\psi ^{1}(\mathbf {x} ),\psi ^{2}(\mathbf {x} ),\psi ^{3}(\mathbf {x} )\}=(q^{1},q^{2},q^{3})\equiv (r,\theta ,z)=\{{\sqrt {x_{1}^{2}+x_{2}^{2}}},\tan ^{-1}(x_{2}/x_{1}),x_{3}\}} 0 < r <   ,     0 < θ < 2 π   ,     < z < {\displaystyle 0<r<\infty ~,~~0<\theta <2\pi ~,~~-\infty <z<\infty }

すると、共変基底ベクトルと反変基底ベクトルは、 方向の単位ベクトルになり ます b 1 = e r = b 1 b 2 = r   e θ = r 2   b 2 b 3 = e z = b 3 {\displaystyle {\begin{aligned}\mathbf {b} _{1}&=\mathbf {e} _{r}=\mathbf {b} ^{1}\\\mathbf {b} _{2}&=r~\mathbf {e} _{\theta }=r^{2}~\mathbf {b} ^{2}\\\mathbf {b} _{3}&=\mathbf {e} _{z}=\mathbf {b} ^{3}\end{aligned}}} e r , e θ , e z {\displaystyle \mathbf {e} _{r},\mathbf {e} _{\theta },\mathbf {e} _{z}} r , θ , z {\displaystyle r,\theta ,z}

計量テンソルの成分は、 基底が直交していることを示すものであることに注意してください。 g i j = g i j = 0 ( i j )   ;     g 11 = 1 ,   g 22 = 1 r ,   g 33 = 1 {\displaystyle g^{ij}=g_{ij}=0(i\neq j)~;~~{\sqrt {g^{11}}}=1,~{\sqrt {g^{22}}}={\cfrac {1}{r}},~{\sqrt {g^{33}}}=1}

第二種クリストッフェル記号の非零成分は Γ 12 2 = Γ 21 2 = 1 r   ;     Γ 22 1 = r {\displaystyle \Gamma _{12}^{2}=\Gamma _{21}^{2}={\cfrac {1}{r}}~;~~\Gamma _{22}^{1}=-r}

物理的なベクトル場を表現する

円筒極座標における正規化された反変基底ベクトルは 、ベクトルの物理的成分は、 b ^ 1 = e r   ;     b ^ 2 = e θ   ;     b ^ 3 = e z {\displaystyle {\hat {\mathbf {b} }}^{1}=\mathbf {e} _{r}~;~~{\hat {\mathbf {b} }}^{2}=\mathbf {e} _{\theta }~;~~{\hat {\mathbf {b} }}^{3}=\mathbf {e} _{z}} v {\displaystyle \mathbf {v} } ( v ^ 1 , v ^ 2 , v ^ 3 ) = ( v 1 , v 2 / r , v 3 ) =: ( v r , v θ , v z ) {\displaystyle ({\hat {v}}_{1},{\hat {v}}_{2},{\hat {v}}_{3})=(v_{1},v_{2}/r,v_{3})=:(v_{r},v_{\theta },v_{z})}

スカラー場の勾配

円筒座標におけるスカラー場の勾配は、曲線座標における一般的な式から計算することができ、次の形をとる。 f ( x ) {\displaystyle f(\mathbf {x} )} f = f r   e r + 1 r   f θ   e θ + f z   e z {\displaystyle {\boldsymbol {\nabla }}f={\cfrac {\partial f}{\partial r}}~\mathbf {e} _{r}+{\cfrac {1}{r}}~{\cfrac {\partial f}{\partial \theta }}~\mathbf {e} _{\theta }+{\cfrac {\partial f}{\partial z}}~\mathbf {e} _{z}}

ベクトル場の勾配

同様に、円筒座標におけるベクトル場の勾配は次のように表される。 v ( x ) {\displaystyle \mathbf {v} (\mathbf {x} )} v = v r r   e r e r + 1 r ( v r θ v θ )   e r e θ + v r z   e r e z + v θ r   e θ e r + 1 r ( v θ θ + v r )   e θ e θ + v θ z   e θ e z + v z r   e z e r + 1 r v z θ   e z e θ + v z z   e z e z {\displaystyle {\begin{aligned}{\boldsymbol {\nabla }}\mathbf {v} &={\cfrac {\partial v_{r}}{\partial r}}~\mathbf {e} _{r}\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left({\cfrac {\partial v_{r}}{\partial \theta }}-v_{\theta }\right)~\mathbf {e} _{r}\otimes \mathbf {e} _{\theta }+{\cfrac {\partial v_{r}}{\partial z}}~\mathbf {e} _{r}\otimes \mathbf {e} _{z}\\[8pt]&+{\cfrac {\partial v_{\theta }}{\partial r}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left({\cfrac {\partial v_{\theta }}{\partial \theta }}+v_{r}\right)~\mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }+{\cfrac {\partial v_{\theta }}{\partial z}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{z}\\[8pt]&+{\cfrac {\partial v_{z}}{\partial r}}~\mathbf {e} _{z}\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}{\cfrac {\partial v_{z}}{\partial \theta }}~\mathbf {e} _{z}\otimes \mathbf {e} _{\theta }+{\cfrac {\partial v_{z}}{\partial z}}~\mathbf {e} _{z}\otimes \mathbf {e} _{z}\end{aligned}}}

ベクトル場の発散

曲線座標系におけるベクトル場の発散の式を用いると、円筒座標系における発散は次のように示される。 v = v r r + 1 r ( v θ θ + v r ) + v z z {\displaystyle {\begin{aligned}{\boldsymbol {\nabla }}\cdot \mathbf {v} &={\cfrac {\partial v_{r}}{\partial r}}+{\cfrac {1}{r}}\left({\cfrac {\partial v_{\theta }}{\partial \theta }}+v_{r}\right)+{\cfrac {\partial v_{z}}{\partial z}}\end{aligned}}}

スカラー場のラプラシアン

ラプラシアンは、 とすればより簡単に計算できる。円筒極座標では 、したがって、 2 f = f {\displaystyle {\boldsymbol {\nabla }}^{2}f={\boldsymbol {\nabla }}\cdot {\boldsymbol {\nabla }}f} v = f = [ v r     v θ     v z ] = [ f r     1 r f θ     f z ] {\displaystyle \mathbf {v} ={\boldsymbol {\nabla }}f=\left[v_{r}~~v_{\theta }~~v_{z}\right]=\left[{\cfrac {\partial f}{\partial r}}~~{\cfrac {1}{r}}{\cfrac {\partial f}{\partial \theta }}~~{\cfrac {\partial f}{\partial z}}\right]} v = 2 f = 2 f r 2 + 1 r ( 1 r 2 f θ 2 + f r ) + 2 f z 2 = 1 r [ r ( r f r ) ] + 1 r 2 2 f θ 2 + 2 f z 2 {\displaystyle {\boldsymbol {\nabla }}\cdot \mathbf {v} ={\boldsymbol {\nabla }}^{2}f={\cfrac {\partial ^{2}f}{\partial r^{2}}}+{\cfrac {1}{r}}\left({\cfrac {1}{r}}{\cfrac {\partial ^{2}f}{\partial \theta ^{2}}}+{\cfrac {\partial f}{\partial r}}\right)+{\cfrac {\partial ^{2}f}{\partial z^{2}}}={\cfrac {1}{r}}\left[{\cfrac {\partial }{\partial r}}\left(r{\cfrac {\partial f}{\partial r}}\right)\right]+{\cfrac {1}{r^{2}}}{\cfrac {\partial ^{2}f}{\partial \theta ^{2}}}+{\cfrac {\partial ^{2}f}{\partial z^{2}}}}

物理的な2階テンソル場を表現する

2階テンソル場の物理的成分は、テンソルを正規化された反変基底で表したときに得られる成分である。円筒極座標では、これらの成分は以下の通りである。

S ^ 11 = S 11 =: S r r , S ^ 12 = S 12 r =: S r θ , S ^ 13 = S 13 =: S r z S ^ 21 = S 21 r =: S θ r , S ^ 22 = S 22 r 2 =: S θ θ , S ^ 23 = S 23 r =: S θ z S ^ 31 = S 31 =: S z r , S ^ 32 = S 32 r =: S z θ , S ^ 33 = S 33 =: S z z {\displaystyle {\begin{aligned}{\hat {S}}_{11}&=S_{11}=:S_{rr},&{\hat {S}}_{12}&={\frac {S_{12}}{r}}=:S_{r\theta },&{\hat {S}}_{13}&=S_{13}=:S_{rz}\\[6pt]{\hat {S}}_{21}&={\frac {S_{21}}{r}}=:S_{\theta r},&{\hat {S}}_{22}&={\frac {S_{22}}{r^{2}}}=:S_{\theta \theta },&{\hat {S}}_{23}&={\frac {S_{23}}{r}}=:S_{\theta z}\\[6pt]{\hat {S}}_{31}&=S_{31}=:S_{zr},&{\hat {S}}_{32}&={\frac {S_{32}}{r}}=:S_{z\theta },&{\hat {S}}_{33}&=S_{33}=:S_{zz}\end{aligned}}}

2階テンソル場の勾配

上記の定義を用いると、円筒極座標における2次テンソル場の勾配は次のように表せることがわかる。 S = S r r r   e r e r e r + 1 r [ S r r θ ( S θ r + S r θ ) ]   e r e r e θ + S r r z   e r e r e z + S r θ r   e r e θ e r + 1 r [ S r θ θ + ( S r r S θ θ ) ]   e r e θ e θ + S r θ z   e r e θ e z + S r z r   e r e z e r + 1 r [ S r z θ S θ z ]   e r e z e θ + S r z z   e r e z e z + S θ r r   e θ e r e r + 1 r [ S θ r θ + ( S r r S θ θ ) ]   e θ e r e θ + S θ r z   e θ e r e z + S θ θ r   e θ e θ e r + 1 r [ S θ θ θ + ( S r θ + S θ r ) ]   e θ e θ e θ + S θ θ z   e θ e θ e z + S θ z r   e θ e z e r + 1 r [ S θ z θ + S r z ]   e θ e z e θ + S θ z z   e θ e z e z + S z r r   e z e r e r + 1 r [ S z r θ S z θ ]   e z e r e θ + S z r z   e z e r e z + S z θ r   e z e θ e r + 1 r [ S z θ θ + S z r ]   e z e θ e θ + S z θ z   e z e θ e z + S z z r   e z e z e r + 1 r   S z z θ   e z e z e θ + S z z z   e z e z e z {\displaystyle {\begin{aligned}{\boldsymbol {\nabla }}{\boldsymbol {S}}&={\frac {\partial S_{rr}}{\partial r}}~\mathbf {e} _{r}\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{rr}}{\partial \theta }}-(S_{\theta r}+S_{r\theta })\right]~\mathbf {e} _{r}\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{rr}}{\partial z}}~\mathbf {e} _{r}\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{z}\\[8pt]&+{\frac {\partial S_{r\theta }}{\partial r}}~\mathbf {e} _{r}\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{r\theta }}{\partial \theta }}+(S_{rr}-S_{\theta \theta })\right]~\mathbf {e} _{r}\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{r\theta }}{\partial z}}~\mathbf {e} _{r}\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{z}\\[8pt]&+{\frac {\partial S_{rz}}{\partial r}}~\mathbf {e} _{r}\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{rz}}{\partial \theta }}-S_{\theta z}\right]~\mathbf {e} _{r}\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{rz}}{\partial z}}~\mathbf {e} _{r}\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{z}\\[8pt]&+{\frac {\partial S_{\theta r}}{\partial r}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{\theta r}}{\partial \theta }}+(S_{rr}-S_{\theta \theta })\right]~\mathbf {e} _{\theta }\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{\theta r}}{\partial z}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{z}\\[8pt]&+{\frac {\partial S_{\theta \theta }}{\partial r}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{\theta \theta }}{\partial \theta }}+(S_{r\theta }+S_{\theta r})\right]~\mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{\theta \theta }}{\partial z}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{z}\\[8pt]&+{\frac {\partial S_{\theta z}}{\partial r}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{\theta z}}{\partial \theta }}+S_{rz}\right]~\mathbf {e} _{\theta }\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{\theta z}}{\partial z}}~\mathbf {e} _{\theta }\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{z}\\[8pt]&+{\frac {\partial S_{zr}}{\partial r}}~\mathbf {e} _{z}\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{zr}}{\partial \theta }}-S_{z\theta }\right]~\mathbf {e} _{z}\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{zr}}{\partial z}}~\mathbf {e} _{z}\otimes \mathbf {e} _{r}\otimes \mathbf {e} _{z}\\[8pt]&+{\frac {\partial S_{z\theta }}{\partial r}}~\mathbf {e} _{z}\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{z\theta }}{\partial \theta }}+S_{zr}\right]~\mathbf {e} _{z}\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{z\theta }}{\partial z}}~\mathbf {e} _{z}\otimes \mathbf {e} _{\theta }\otimes \mathbf {e} _{z}\\[8pt]&+{\frac {\partial S_{zz}}{\partial r}}~\mathbf {e} _{z}\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{r}+{\cfrac {1}{r}}~{\frac {\partial S_{zz}}{\partial \theta }}~\mathbf {e} _{z}\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{\theta }+{\frac {\partial S_{zz}}{\partial z}}~\mathbf {e} _{z}\otimes \mathbf {e} _{z}\otimes \mathbf {e} _{z}\end{aligned}}}

2階テンソル場の発散

円筒極座標における2階テンソル場の発散は、勾配の式から、2項積における2つの外側のベクトルのスカラー積が非ゼロとなる項を集めることによって得られる。したがって、 S = S r r r   e r + S r θ r   e θ + S r z r   e z + 1 r [ S r θ θ + ( S r r S θ θ ) ]   e r + 1 r [ S θ θ θ + ( S r θ + S θ r ) ]   e θ + 1 r [ S θ z θ + S r z ]   e z + S z r z   e r + S z θ z   e θ + S z z z   e z {\displaystyle {\begin{aligned}{\boldsymbol {\nabla }}\cdot {\boldsymbol {S}}&={\frac {\partial S_{rr}}{\partial r}}~\mathbf {e} _{r}+{\frac {\partial S_{r\theta }}{\partial r}}~\mathbf {e} _{\theta }+{\frac {\partial S_{rz}}{\partial r}}~\mathbf {e} _{z}\\[8pt]&+{\cfrac {1}{r}}\left[{\frac {\partial S_{r\theta }}{\partial \theta }}+(S_{rr}-S_{\theta \theta })\right]~\mathbf {e} _{r}+{\cfrac {1}{r}}\left[{\frac {\partial S_{\theta \theta }}{\partial \theta }}+(S_{r\theta }+S_{\theta r})\right]~\mathbf {e} _{\theta }+{\cfrac {1}{r}}\left[{\frac {\partial S_{\theta z}}{\partial \theta }}+S_{rz}\right]~\mathbf {e} _{z}\\[8pt]&+{\frac {\partial S_{zr}}{\partial z}}~\mathbf {e} _{r}+{\frac {\partial S_{z\theta }}{\partial z}}~\mathbf {e} _{\theta }+{\frac {\partial S_{zz}}{\partial z}}~\mathbf {e} _{z}\end{aligned}}}

参照

参考文献

注記
  1. ^ abc Green, AE; Zerna, W. (1968).理論的弾性. オックスフォード大学出版局. ISBN 0-19-853486-8
  2. ^ abc Ogden, RW (2000).非線形弾性変形. ドーバー.
  3. ^ Naghdi, PM (1972). 「殻と板の理論」. S. Flügge (編).物理学ハンドブック. 第6a巻第2号. pp.  425– 640.
  4. ^ abcdefghijk Simmonds, JG (1994).テンソル解析の概要. Springer. ISBN 0-387-90639-8
  5. ^ ab Basar, Y.; Weichert, D. (2000).固体の数値連続体力学:基本概念と展望. Springer.
  6. ^ abc Ciarlet, PG (2000).シェルの理論. 第1巻. エルゼビア・サイエンス.
  7. ^ アインシュタイン、A. (1915). 「一般相対性理論への貢献」. ラチョス、C. (編). 『アインシュタインの10年』 p. 213.
  8. ^ ミスナー, CW; ソーン, KS; ウィーラー, JA (1973). 『重力』 WHフリーマン・アンド・カンパニーISBN 0-7167-0344-0
  9. ^ Greenleaf, A.; Lassas, M.; Uhlmann, G. (2003). 「EITでは検出できない異方性導電率」.生理学的測定. 24 (2): 413– 419. doi :10.1088/0967-3334/24/2/353. PMID  12812426. S2CID  250813768.
  10. ^ Leonhardt, U.; Philbin, TG (2006). 「電気工学における一般相対性理論」. New Journal of Physics . 8 (10): 247. arXiv : cond-mat/0607418 . Bibcode :2006NJPh....8..247L. doi :10.1088/1367-2630/8/10/247. S2CID  12100599.
  11. ^ 「テンソル場の発散」弾性/テンソル入門.ウィキバーシティ. 2010年11月26日閲覧
さらに読む
  • シュピーゲル、MR (1959)。ベクトル解析。ニューヨーク: シャウムのアウトライン シリーズ。ISBN 0-07-084378-3 {{cite book}}: ISBN / Date incompatibility (help)
  • アーフケン、ジョージ(1995年)『物理学者のための数学的手法』アカデミック・プレス、ISBN 0-12-059877-9
  • 曲線座標における単位ベクトルの導出
  • MathWorldの曲線座標に関するページ
  • R. ブランノン教授の曲線座標に関する電子書籍
Retrieved from "https://en.wikipedia.org/w/index.php?title=Tensors_in_curvilinear_coordinates&oldid=1322697727"