直交曲線座標

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

テンプレート:Redirect テンプレート:翻訳中途 数学において、直交曲線座標(ちょっこうきょくせんざひょう)、直交座標(ちょっこうざひょう、テンプレート:Lang-en-short)とは、座標超曲面同士が互いに直交するようなd個の座標q = (q1, q2, ..., qd)の組として定義される(注:上付き添え字は指数ではなく添え字 (Einstein notation) を意味する)。ある座標qkに対する座標超曲面とは、qkが定数となる超曲面(場合によっては曲線、曲面)のことである。たとえば、3次元のデカルト座標系 (x, y, z) では「x = 定数」、「y = 定数」、「z = 定数」は座標超曲面であるが、これらが互いに直角に交るので、直交座標系である。直交曲線座標は曲線座標の特殊な例である。

動機

矩形グリッドに施される等角写像。曲がったグリッドの直交性は保持されている。

ベクトル同士の演算や物理法則の導出は、通常、デカルト座標系 (テンプレート:Lang-en-short) で行うのが最も簡単であるが、量子力学における場の理論、流体力学、等角性を保持する地図投影正角図法)、電気力学プラズマ物理学化学種拡散等において生じるような境界値問題においては、デカルト座標ではない直交座標(非デカルト直交座標、即ち直交座標)がよく用いられる。

非デカルト直交座標の利点は、問題の対称性に合わせて座標を選ぶことができる点である。例えば、地面(あるいは他の障害物)から遠く離れた場所での爆発による圧力波は、デカルト座標では3次元空間に依存するが、球座標では問題はほぼ1次元となる(圧力波は時間と中心からの距離のみに依存するため)。デカルト座標では偏微分方程式を含む2次元の境界値問題を解かなければならないが、円筒座標では偏微分方程式を用いずとも、常微分方程式で表現可能1次元の問題に帰着される。

一般的なテンプレート:仮リンク (テンプレート:Lang-en-short) ではなく、直交曲線座標を好まれる理由は、これを用いたほうが単純であるからである。直交しない座標では多くの複雑な問題が発生する。例えば、直交曲線座標では、多くの問題が変数分離によって解決されることがある。変数分離とは、複雑なd次元の問題を、「既知の関数で解くことができるd個の1次元の問題」に変換する数学的手法である。多くの方程式は、ラプラス方程式ヘルムホルツ方程式に還元することができる。ラプラス方程式は、下表13番に示す座標系(トロイダル座標を除く下表14番に示す座標)で変数分離可能であり、ヘルムホルツ方程式は、下表11番の座標系で変数分離可能である[1][2]

直交曲線座標は、計量テンソルの非対角項を決して持たない。つまり、無限小の2乗距離、即ち ds2 は、常に、「無限小の座標変位の2乗の総和」として書くことができる。

即ち、:ds2=k=1d(hkdqk)2

ここで、d は次次元を表す。 また、スケーリング関数(スケールファクター、h-パラメータ)

hk(𝐪) =def gkk(𝐪)=|𝐞k|

は、計量テンソルの対角成分の平方根(後述する局所基底ベクトル (テンプレート:Lang-en-short) 𝐞kの長さ)に等しい。これらのスケーリング関数 hi は新しい座標における微分演算子、例え 勾配ラプラシアン発散回転を計算する上でも使用される。

2次元の直交曲線座標の一例を生成する簡単な方法として、標準的なデカルト座標テンプレート:Nowrapが定める2次元格子の共形写像による方法がある。非ゼロの複素微分を持つ正則関数 w = f(z) は共形写像を生成する。得られた複素数をテンプレート:Nowrapと書くと、元の定数 xy の直線と同じように定数 uv の曲線は直交する。

3次元以上の直交曲線座標の一例を生成する方法の一つとして、直交する2次元座標系から、新しい次元に投影する(「円筒座標」)か、2次元座標系をその対称軸の1つを中心に回転させる方法がある。 しかし、2次元座標系を射影したり回転させたりしても得られない3次元の直交曲線座標系もあり、例えばテンプレート:仮リンクはそのような例である。より一般的な直交曲線座標は、いくつかの必要な座標面から出発し、その直交軌道を考えることで得られることがある。

基底ベクトル

共変基底(Covariant basis)

デカルト座標では、基底ベクトルは固定(一定)である。より一般的なテンプレート:仮リンクでは、座標によって空間の点が指定され、そのような点ごとに基底ベクトルの集合が束ねられるが、それは一般に一定ではない。直交曲線座標の特徴は、基底ベクトルが変化しても、互いに対して常に直交していることである。言い換えれば

𝐞i𝐞j=0ifij

これらの基底ベクトル「ある座標を変化させ、他の座標を固定して得られる曲線の接ベクトル」として定義される。即ち、

𝐞i=𝐫qi
2次元直交座標の可視化。1つの座標以外を一定にして得られる曲線が、基底ベクトルとともに示されている。基底ベクトルは長さが等しくないことに注意すること:等しい必要はなく、直交していればよい。

ここで 'r は何らかの点を表し、q i は基底ベクトルを抽出した座標である。つまり、1つの座標以外を固定して曲線を得、固定しない座標をパラメトリック曲線のように変化させ、パラメータ(変化する座標)に対する曲線の微分を、その座標の基底ベクトルとする。

なお、ベクトルは必ずしも等しい長さとは限らない。座標のスケールファクターとして知られる便利な関数は、単に基底ベクトル 𝐞i の長さ hi (下表参照) である。スケールファクターはLamé coefficientsと呼ばれることもあるが、弾性論におけるラメ定数と混同しないように注意のこと。

単位ベクトル(正規化基底ベクトル)をハット付きで表記し、これは、上記の 𝐞i を、その長さで割ることで得られる。

𝐞^i=𝐞ihi=𝐞i|𝐞i|

ベクトル場は「基底ベクトル」(𝐞i) に対する成分で指定される場合と、「正規化された基底ベクトル」(𝐞^i) に対する成分で指定される場合があり、どちらの場合を指しているのかを確認する必要がある。正規化基底の成分は、数量を明確にする目的では最も一般的に使われる。(例えば、接線速度にスケールファクターをかけたものではなく、接線速度を扱いたい場合がある)。しかし、微分する場合にはより複雑になるため正規化基底はあまり一般的に使わない。

反変基底(Contravariant basis)

前節に示した基底ベクトル(𝐞^i)は共変基底ベクトルといわれるが、それはベクトルと「共変」するからである。直交曲線座標の場合、反変基底ベクトル (𝐞i) は、共変ベクトルと同じ方向(但し逆長)になるので、簡単に求められる(このため、2組の基底ベクトルは互いに対して逆であると言われる)即ち、

𝐞i=𝐞^ihi=𝐞ihi2

クロネッカーのデルタを使うと、𝐞i=δij となることに注意のこと。

また、

𝐞^i=𝐞ihi=hi𝐞i=𝐞^i

我々は、直交曲線座標上の「ベクトル」を記述するためによく使われる3つの異なる基底セット、 即ち、共変基底ei、反変基底ei、正規化基底êiの3つの基底を持つ。「ベクトル」はobjective quantity,であり、その同一性はどの座標系にも依存しないが、「ベクトル」の成分はそのベクトルがどの基底で表現されるかに依存する。

𝐱=ixi𝐞i=ixi𝐞i

添字の位置は成分の計算方法を表している(上付の添字を指数と混同してはいけない)。なお、すべての基底ベクトル(i=1, 2, ..., d)に対する和を示す記号Σ(大文字のSigma)と和の範囲は、しばしば省略(アインシュタイン表記)されることがある。それぞれの基底における成分同士の関係は、以下のようになる。

hi2xi=xi

正規化基底に関するベクトルの成分を指定するために広く使われている表記法はない。本稿では、ベクトル成分には添え字を用い、成分が正規化基底で計算されていることに着目する。

ベクトル代数

ベクトルの加算とマイナスは、デカルト座標と同様に成分毎に行うことが出来、複雑な操作は不要である。他のベクトル演算については、特別な配慮が必要な場合がある。 ただし、これらの演算はすべてベクトル場の2つのベクトルが同じ点に束縛されている(言い換えれば、ベクトルの尾が一致している)ことを前提としていることに注意のこと。基底ベクトルは一般に直交曲線座標で変化するため、空間上の異なる点で計算された成分を持つ2つのベクトルを足し合わせる場合、基底ベクトルの違いを考慮する必要がある。

内積(Dot product)

デカルト座標系における内積においては、(ユークリッド空間に直交基底を考えると)単純に成分の積の和になる。同様に、直交曲線座標でも、2つのベクトル xy の内積は、ベクトルの成分を正規化基底で表示すると、このような馴染みのある形になる。

𝐱𝐲=ixi𝐞^ijyj𝐞^j=ixiyi

これは、ある点での正規化基底がデカルト座標系を形成できるという事実の直接的な帰結である。この基底は正規直交基底である。

𝐱𝐲=ihi2xiyi=ixiyihi2=ixiyi=ixiyi

これは、ベクトルを成分形式で書き出し、基底ベクトルを正規化し、内積を取ることで容易に導き出すことができる。例えば、2Dの場合、

𝐱𝐲=(x1𝐞1+x2𝐞2)(y1𝐞1+y2𝐞2)=(x1h1𝐞^1+x2h2𝐞^2)(y1𝐞^1h1+y2𝐞^2h2)=x1y1+x2y2

ここでは、正規化された共変基底と反変基底が等しいことが利用されている。

外積(Cross product)

3次元デカルト座標における外積は、以下の通りである。

𝐱×𝐲=(x2y3x3y2)𝐞^1+(x3y1x1y3)𝐞^2+(x1y2x2y1)𝐞^3

そして、直交曲線座標系でも、成分を正規化した基準で計算すれば、上記の式は有効である。

直交曲線座標において、共変基底あるいは反変基底を考えた場合の外積を構成するには、やはり基底ベクトルを正規化する必要がある。例えば、

𝐱×𝐲=ixi𝐞i×jyj𝐞j=ixihi𝐞^i×jyjhj𝐞^j

さらに展開すれば、(𝐞1,𝐞2,𝐞3) が右手系であるという仮定の下で

𝐱×𝐲=(x2y3x3y2)h2h3h1𝐞1+(x3y1x1y3)h1h3h2𝐞2+(x1y2x2y1)h1h2h3𝐞3

直交しない座標や高次元への一般化を単純化するために、外積の簡潔な表記がレビ・チビタテンソルで可能であるが、スケールファクターがすべて1に等しくない場合、0と1以外の成分を持つことになる。

ベクトル解析

微分

ある点からの無限小の変位を見てみると、明らかに以下が成り立つ。

d𝐫=i𝐫qidqi=i𝐞idqi

定義によれば、関数の勾配は以下を満たさなければならない(この定義はƒが任意のテンソルであっても真である)。

df=fd𝐫df=fi𝐞idqi

従って、ナブラ演算子は必ず、以下を満たさねばならないことになる。

=i𝐞iqi

これは、これは直交曲線座標に限らない一般的な曲線座標の場合にも当てはまる。勾配ラプラシアンのような演算子は、この演算子を適切に適用することで得られるものである。

基底ベクトルの式(Basis vector formulae)

d'r と正規化基底ベクトル êiから、次のように構成できる[3][4]

Differential element Vectors Scalars
線要素 Tangent vector to coordinate curve qi:

d=hidqi𝐞^i=𝐫qidqi

Infinitesimal length

d=d𝐫d𝐫=(h1dq1)2+(h2dq2)2+(h3dq3)2

面積要素 Normal to coordinate surface qk = constant:

d𝐒=(hidqi𝐞^i)×(hjdqj𝐞^j)=dqidqj(𝐫qi×𝐫qj)=hihjdqidqj𝐞^k

Infinitesimal surface

dSk=hihjdqidqj

体積要素 N/A Infinitesimal volume

dV=|(h1dq1𝐞^1)(h2dq2𝐞^2)×(h3dq3𝐞^3)|=|𝐞^1𝐞^2×𝐞^3|h1h2h3dq1dq2dq3=h1h2h3dq1dq2dq3=Jdq1dq2dq3

ここで、

J=𝐫q1(𝐫q2×𝐫q3)=(x,y,z)(q1,q2,q3)=h1h2h3>0

ヤコビ行列式で、これは「デカルト座標における無限小の立方体 dxdydz」から、「無限小の曲った立方体」への体積の変形という幾何学的解釈を持つものである。ただしここで ヤコビ行列式は正と仮定してあることに注意する (一般にはヤコビ行列式は0でなければよいので負の値も許される)。以下ではヤコビ行列式が正の場合のみ考える。

積分

上に示した線素を用いると、ベクトル F の経路 𝒫 に沿った線積分は次のようになる。

𝒫𝐅d𝐫=𝒫iFi𝐞ij𝐞jdqj=i𝒫Fidqi

1つの座標qkを一定にして記述した面の面積の無限小要素は、以下のように変換され、

dAk=ikdsi=ikhidqi

同様に、体積要素も以下のように変換される。

dV=idsi=ihidqi

ここで、大きな記号Π(πの大文字)は、総乗を示す。即ち、すべてのスケールファクターの積はヤコビ行列式に等しいことを意味している。

例として、3次元のq1 = 定数で定まる面𝒮上のベクトル値関数F面積分は次のようになる。

𝒮𝐅d𝐀=𝒮𝐅𝐧^ dA=𝒮𝐅𝐞^1 dA=𝒮F1h2h3h1dq2dq3

ただし、F1/h1は、Fの、この表面に垂直な成分である。

Differential operators in three dimensions

テンプレート:Main これらの演算は応用上共通なので、本節ではすべてのベクトル成分を正規化基底を用いて以下のように示す。 Fi=𝐅𝐞^i.

Operator Expression
Gradient of a scalar field ϕ=𝐞^1h1ϕq1+𝐞^2h2ϕq2+𝐞^3h3ϕq3
Divergence of a vector field 𝐅=1h1h2h3[q1(F1h2h3)+q2(F2h3h1)+q3(F3h1h2)]
Curl of a vector field ×𝐅=𝐞^1h2h3[q2(h3F3)q3(h2F2)]+𝐞^2h3h1[q3(h1F1)q1(h3F3)]+𝐞^3h1h2[q1(h2F2)q2(h1F1)]=1h1h2h3|h1𝐞^1h2𝐞^2h3𝐞^3q1q2q3h1F1h2F2h3F3|
Laplacian of a scalar field 2ϕ=1h1h2h3[q1(h2h3h1ϕq1)+q2(h3h1h2ϕq2)+q3(h1h2h3ϕq3)]

上記の式は、レヴィ=チヴィタ記号を用いてより簡潔に書くことができる。ϵijkとヤコビ行列式J=h1h2h3>0で、繰り返し添字に対する和を考える。

Operator Expression
Gradient of a scalar field ϕ=𝐞^khkϕqk
Divergence of a vector field 𝐅=1Jqk(JhkFk)
Curl of a vector field (3D only) ×𝐅=hk𝐞^kJϵijkqi(hjFj)
Laplacian of a scalar field 2ϕ=1Jqk(Jhk2ϕqk)

また、スカラー場の勾配は正準偏導関数を含むヤコビ行列式 J で表現できることに注意。

𝐉=[ϕq1,ϕq2,ϕq3]

upon a change of basis:

ϕ=𝐒𝐑T𝐉T

where the rotation and scaling matrices are:

𝐑=[𝐞1,𝐞2,𝐞3]
𝐒=diag([h11,h21,h31]).

直交曲線座標の表

通常の直交曲線座標の他に、いくつかのやや珍しい直交曲線座標を以下に表に示す[4]Interval notation is used for compactness in the coordinates column.

Curvillinear coordinates (q1, q2, q3) Transformation from cartesian (x, y, z) Scale factors
Spherical polar coordinates

(r,θ,ϕ)[0,)×[0,π]×[0,2π)

x=rsinθcosϕy=rsinθsinϕz=rcosθ h1=1h2=rh3=rsinθ
Cylindrical polar coordinates

(r,ϕ,z)[0,)×[0,2π)×(,)

x=rcosϕy=rsinϕz=z h1=h3=1h2=r
Parabolic cylindrical coordinates

(u,v,z)(,)×[0,)×(,)

x=12(u2v2)y=uvz=z h1=h2=u2+v2h3=1
Parabolic coordinates

(u,v,ϕ)[0,)×[0,)×[0,2π)

x=uvcosϕy=uvsinϕz=12(u2v2) h1=h2=u2+v2h3=uv
Paraboloidal coordinates

(λ,μ,ν)[0,b2)×(b2,a2)×(a2,)b2<a2

x2qia2+y2qib2=2z+qi

where (q1,q2,q3)=(λ,μ,ν)

hi=12(qjqi)(qkqi)(a2qi)(b2qi)
Ellipsoidal coordinates

(λ,μ,ν)[0,c2)×(c2,b2)×(b2,a2)λ<c2<b2<a2,c2<μ<b2<a2,c2<b2<ν<a2,

x2a2qi+y2b2qi+z2c2qi=1

where (q1,q2,q3)=(λ,μ,ν)

hi=12(qjqi)(qkqi)(a2qi)(b2qi)(c2qi)
Elliptic cylindrical coordinates

(u,v,z)[0,)×[0,2π)×(,)

x=acoshucosvy=asinhusinvz=z h1=h2=asinh2u+sin2vh3=1
Prolate spheroidal coordinates

(ξ,η,ϕ)[0,)×[0,π]×[0,2π)

x=asinhξsinηcosϕy=asinhξsinηsinϕz=acoshξcosη h1=h2=asinh2ξ+sin2ηh3=asinhξsinη
Oblate spheroidal coordinates

(ξ,η,ϕ)[0,)×[π2,π2]×[0,2π)

x=acoshξcosηcosϕy=acoshξcosηsinϕz=asinhξsinη h1=h2=asinh2ξ+sin2ηh3=acoshξcosη
Bipolar cylindrical coordinates

(u,v,z)[0,2π)×(,)×(,)

x=asinhvcoshvcosuy=asinucoshvcosuz=z h1=h2=acoshvcosuh3=1
Toroidal coordinates

(u,v,ϕ)(π,π]×[0,)×[0,2π)

x=asinhvcosϕcoshvcosuy=asinhvsinϕcoshvcosuz=asinucoshvcosu h1=h2=acoshvcosuh3=asinhvcoshvcosu
Bispherical coordinates

(u,v,ϕ)(π,π]×[0,)×[0,2π)

x=asinucosϕcoshvcosuy=asinusinϕcoshvcosuz=asinhvcoshvcosu h1=h2=acoshvcosuh3=asinucoshvcosu
Conical coordinates

(λ,μ,ν)ν2<b2<μ2<a2λ[0,)

x=λμνaby=λa(μ2a2)(ν2a2)a2b2z=λb(μ2b2)(ν2b2)b2a2 h1=1h22=λ2(μ2ν2)(μ2a2)(b2μ2)h32=λ2(μ2ν2)(ν2a2)(ν2b2)

関連項目

脚注

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

参考文献

  • Korn GA and Korn TM. (1961) Mathematical Handbook for Scientists and Engineers, McGraw-Hill, pp. 164–182.
  • テンプレート:Cite journal
  • Margenau H. and Murphy GM. (1956) The Mathematics of Physics and Chemistry, 2nd. ed., Van Nostrand, pp. 172–192.
  • Leonid P. Lebedev and Michael J. Cloud (2003) Tensor Analysis, pp. 81 – 88.
  1. テンプレート:Mathworld
  2. テンプレート:Harvnb
  3. Mathematical Handbook of Formulas and Tables (3rd edition), S. Lipschutz, M.R. Spiegel, J. Liu, Schuam's Outline Series, 2009, テンプレート:ISBN2.
  4. 4.0 4.1 Vector Analysis (2nd Edition), M.R. Spiegel, S. Lipschutz, D. Spellman, Schaum’s Outlines, McGraw Hill (USA), 2009, テンプレート:ISBN2