直交曲線座標

数学において、直交曲線座標、直交座標とは、座標超曲面同士が互いに直交するようなd個の座標q = の組として定義される。ある座標//q////k//に対する座標超曲面とは、//q////k//が定数となる超曲面のことである。たとえば、3次元のデカルト座標系  では「x = 定数」、「y = 定数」、「z = 定数」は座標超曲面であるが、これらが互いに直角に交るので、直交座標系である。直交曲線座標は曲線座標の特殊な例である。

1. 動機
-----

ベクトル同士の演算や物理法則の導出は、通常、デカルト座標系  で行うのが最も簡単であるが、量子力学における場の理論、流体力学、等角性を保持する地図投影、電気力学、プラズマ物理学、化学種や熱の拡散等において生じるような境界値問題においては、デカルト座標ではない直交座標がよく用いられる。非デカルト直交座標の利点は、問題の対称性に合わせて座標を選ぶことができる点である。例えば、地面から遠く離れた場所での爆発による圧力波は、デカルト座標では3次元空間に依存するが、球座標では問題はほぼ1次元となる。デカルト座標では偏微分方程式を含む2次元の境界値問題を解かなければならないが、円筒座標では偏微分方程式を用いずとも、常微分方程式で表現可能1次元の問題に帰着される。

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

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

即ち、:

ここで、//d// は次次元を表す。 また、スケーリング関数

は、計量テンソルの対角成分の平方根に等しい。これらのスケーリング関数 //h////i// は新しい座標における微分演算子、例え 勾配、ラプラシアン、発散や回転を計算する上でも使用される。

2次元の直交曲線座標の一例を生成する簡単な方法として、標準的なデカルト座標}}が定める2次元格子の共形写像による方法がある。非ゼロの複素微分を持つ正則関数 //w// = //f// は共形写像を生成する。得られた複素数を1=//w// = //u// + //iv//}}と書くと、元の定数 //x// と //y// の直線と同じように定数 //u// と //v// の曲線は直交する。

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

2. 基底ベクトル
---------

2.1. 共変基底
---------

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

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

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

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

単位ベクトルをハット付きで表記し、これは、上記の を、その長さで割ることで得られる。

ベクトル場は「基底ベクトル」 に対する成分で指定される場合と、「正規化された基底ベクトル」 に対する成分で指定される場合があり、どちらの場合を指しているのかを確認する必要がある。正規化基底の成分は、数量を明確にする目的では最も一般的に使われる。。しかし、微分する場合にはより複雑になるため正規化基底はあまり一般的に使わない。

2.2. 反変基底
---------

前節に示した基底ベクトルは共変基底ベクトルといわれるが、それはベクトルと「共変」するからである。直交曲線座標の場合、反変基底ベクトル  は、共変ベクトルと同じ方向になるので、簡単に求められる即ち、

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

また、

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

添字の位置は成分の計算方法を表している。なお、すべての基底ベクトルに対する和を示す記号Σと和の範囲は、しばしば省略されることがある。それぞれの基底における成分同士の関係は、以下のようになる。

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

3. ベクトル代数
---------

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

3.1. 内積
-------

デカルト座標系における内積においては、単純に成分の積の和になる。同様に、直交曲線座標でも、2つのベクトル x と y の内積は、ベクトルの成分を正規化基底で表示すると、このような馴染みのある形になる。

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

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

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

3.2. 外積
-------

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

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

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

さらに展開すれば、 が右手系であるという仮定の下で

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

4. ベクトル解析
---------

4.1. 微分
-------

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

定義によれば、関数の勾配は以下を満たさなければならない。

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

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

4.2. 基底ベクトルの式
------------------------------------

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

{| class="wikitable"|- !scope=col width="10px"| Differential element !scope=col width="200px"| Vectors !scope=col width="200px"| Scalars |- |線要素 ||Tangent vector to coordinate curve //qi//: || Infinitesimal length |- |面積要素 || Normal to coordinate surface //qk// = constant: || Infinitesimal surface |- |体積要素 ||N/A ||Infinitesimal volume |- |}

ここで、

はヤコビ行列式で、これは「デカルト座標における無限小の立方体 d//x//d//y//d//z//」から、「無限小の曲った立方体」への体積の変形という幾何学的解釈を持つものである。ただしここで ヤコビ行列式は正と仮定してあることに注意する 。以下ではヤコビ行列式が正の場合のみ考える。

4.3. 積分
-------

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

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

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

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

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

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

5. Differential operators in three dimensions
---------------------------------------------

これらの演算は応用上共通なので、本節ではすべてのベクトル成分を正規化基底を用いて以下のように示す。.

{| class="wikitable"|- ! Operator ! Expression |- | Gradient of a scalar field | |- |Divergence of a vector field | |- |Curl of a vector field | |- | Laplacian of a scalar field | |}

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

{| class="wikitable"|- ! Operator ! Expression |- | Gradient of a scalar field | |- |Divergence of a vector field | |- |Curl of a vector field  | |- | Laplacian of a scalar field | |}

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

upon a change of basis:

where the rotation and scaling matrices are:

6. 直交曲線座標の表
-----------

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

{| class="wikitable"|- !scope=col width="20px"| Curvillinear coordinates  !scope=col width="200px"| Transformation from cartesian  !scope=col width="200px"| Scale factors |- | Spherical polar coordinates | | |- | Cylindrical polar coordinates

| | |- | Parabolic cylindrical coordinates | | |- | Parabolic coordinates | | |- | Paraboloidal coordinates | where | |- | Ellipsoidal coordinates | where | |- | Elliptic cylindrical coordinates | | |- | Prolate spheroidal coordinates | | |- | Oblate spheroidal coordinates | | |- | Bipolar cylindrical coordinates | | |- | Toroidal coordinates | | |- | Bispherical coordinates | | |- | Conical coordinates | | |- |}

7. 関連項目
-------

Geodetic coordinates
テンソル
ベクトル場
Skew coordinates

8. 脚注
-----

9. 参考文献
-------

Korn GA and Korn TM.  //Mathematical Handbook for Scientists and Engineers//, McGraw-Hill, pp.   164–182.
Margenau H. and Murphy GM.  //The Mathematics of Physics and Chemistry//, 2nd. ed., Van Nostrand, pp.   172–192.
Leonid P. Lebedev and Michael J. Cloud  //Tensor Analysis//, pp.   81 – 88.