次元解析

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

次元解析(じげんかいせき、テンプレート:Lang-en-short)とは、物理量における、長さ質量時間電荷などの次元から、複数の物理量の間の関係を予測することである。

物理的な関係を表す数式においては、両辺や各項の次元が一致しなくてはならない。この規則を逆に利用すると、既知の量を組み合わせ、求めたい未知の物理量の次元に一致するように式を立てれば、それは正しい関係式になっている可能性が高い。これをレイリーの方法テンプレート:En)という[1]

次元解析を用いると、一般解を得ることが困難な(ときには不可能な)現象に対して、物理量間の関係を推測することができる。一方で、未知の物理量を決めるのに既知の物理量では不十分な場合にそれとわかることもある[2]。また、次元の不一致といったミスの防止にも役立つ。

次元一致の原理

数式の左右両辺の各項の次元が等しい式は次元的に健全[3]または次元的に斉一テンプレート:En[4]であると呼ばれる。物理法則に基いて理論的に導かれる理論式は次元的に健全であり、次元的に健全な式のみ物理では意味があると考える。すなわち物理現象を支配する物理方程式の各項の次元は次元的に健全でなければならない。この原理を次元一致の原理テンプレート:En)という[5]

数学的表現

物理量Qn 個の物理量xi によって決定されるとき、それらの関係を表す式

Q=F(x1,,xn)

が次元的に健全であるということは、次のように変形できることを意味する[6]

F(x1,,xn)=i[Xi]ai×F*(x1*,,xn*)

ここで[Xi]は物理量xiの単位または次元、*付きの変数は無次元量を意味する。

バッキンガムのπ定理

バッキンガムのπ定理Buckingham Π theorem)とは、数理物理学の分野において、次元解析の基礎となる理論である。大雑把に言うと、物理的な関係式が物理変数n 個含み、それらの変数がk 種類の独立な基本単位を持つならば、その式は元の物理変数で構成されるp = n - k 個の無次元パラメータを含む式と等価であるという定理である。この定理により、与えられた物理変数から、たとえ関係式の形が不明であっても無次元パラメータを求めることができる。物理量を無次元量で書き直せば、式の次元の一致・不一致をチェックする必要がなくなり、解析が簡単になる。ただし、無次元パラメータの選び方は一意ではない。バッキンガムのΠ定理は無次元パラメータを求める方法を与えるだけであり、物理的に意味のあるものを選ぶわけではない。

2つの物理的な系の無次元パラメータが一致するとき、それらの系は相似であるという(大きさのみが異なる三角形を相似と呼ぶのと同様である)。これらの系は数学的には等価であるため、解析をするために便利な(実験などがしやすい)系を選ぶことができる。

より正確に表現すると、無次元パラメータの個数p は次元行列M退化次数 null M に等しく、k はその階数 rank M に等しい。物理的に異なる系に対して、無次元パラメータが等しくなるなら、それらの系は数学的に等価である。

定式化

次のような物理的な関係式があるとする:

f(q1,q2,,qn)=0

ここでq1, ..., qnn 個の物理変数であり、k 種類の独立な基本単位で表されている。このとき、上式は次の数学的に等価な式に書き換えることができる:

F(π1,π2,,πp)=0

ここでπ1, ..., πpq1, ..., qn で構成されるp = n - k 個の 無次元パラメータである:

πi=q1a1q2a2qnan,i=1,,p

ここで指数ai有理数である(適当にべき乗すれば常に整数としてよい)。

証明

テンプレート:節スタブ

概要

前提として、与えられた基本単位は有理数体上のベクトル空間(物理次元ベクトル空間と呼ぶ)の基底であり、物理単位の積はベクトルの和で表され、べき乗はスカラー倍を表すとする。有次元の物理変数を必要な基本単位の指数の組で表す(現れない基本単位に対しては指数はゼロとする)。例えば、重力加速度g𝖫𝖳2(長さ÷時間2)の次元を持つ。したがってこれは基底(長さ,時間)に関してベクトル(1, -2)で表される。

物理的単位を物理的関係式の両辺で一致させることは、物理次元ベクトル空間で線形従属性を課すこととみなすことができる。

正式な証明

有次元の物理変数n 個で表される系を考える。基本単位はk 種類とする。次元行列 M ∈ Rk×n を(i , j )成分がj 番目の物理変数のi 番目の基本単位の指数である行列とする。例えば

M=(a1an)

は物理変数 テンプレート:Math の次元行列である。

無次元量は単位のべきが全てゼロとなる(すなわち次元がない)組み合わせであり、次元行列の零空間に相当する。無次元変数は有次元変数間の単位の線型結合である。

階数・退化次数公式により、k 個の(必要な)次元を持つn 個のベクトルから成る系は関係のp (= n - k )-次元空間を満足する。任意の基底の選択はp 個の無次元数の要素を持つ。

無次元変数は(分母を払うことで)いつも有次元変数の整数の組み合わせになるように取られる。不自然な有次元数の選択が数学的にはある。いくつかの無次元変数の選択は物理的により意味があり、理想的に使われるものがある。

調和振動

例としてばねにつないだ物体の振動運動について考える。水平面上に質量 テンプレート:Math の物体をおき、垂直に立った壁と物体との間をばね定数 テンプレート:Mathばねで結ぶ。ばねの自然長の状態から物体を テンプレート:Math だけずらし、静かに手を離すと物体は振動運動を始める。このときの振動の周期(1振動にかかる時間)テンプレート:Math を与える式を推測する。水平面との摩擦や空気抵抗は考えない。

式に含まれるであろう定数は、物体の質量 テンプレート:Math、ばね定数 テンプレート:Math、初期変位 テンプレート:Math の3つである。長さの次元を𝖫、質量の次元𝖬、時間の次元を𝖳とすれば、それぞれの定数および周期 テンプレート:Math の次元は[m]=𝖬,[k]=𝖬𝖳2,[x]=𝖫,[T]=𝖳である。この中で長さの次元𝖫を含んでいるのは初期変位 テンプレート:Math のみなので、式に含めることができない。なぜなら式の左辺と右辺では次元が一致しなくてはならず、初期変位を含めるならば両辺に同じだけかける必要があり、それならば無くても同じだからである。

次元が𝖳になるように テンプレート:Mathテンプレート:Math を組み合わせる方法は一つしかない。結果次の式が求まる。

T=Amk

比例係数 テンプレート:Math無次元量の定数で次元解析から求めることはできない。この運動の運動方程式を直接解くと周期は

T=2πmk

となり、テンプレート:Math のもとで両者は見事に一致している(固有振動も参照)。このように簡単な問題ならば次元を考えるだけで見通しが立つ。式の次元が合うことは必須の要請であるので、式の間違いをチェックする場合にも使える。

バッキンガムのΠ定理にしたがって考えると、物理量が テンプレート:Math および テンプレート:Math の4つで、次元が𝖬,𝖳,𝖫の3種類なので、次元行列は

M=(mkxT𝖬1100𝖳0201𝖫0010)

となる(便宜的に列が テンプレート:Math 、行が𝖬,𝖳,𝖫に対応していることを明記しているが、本来の次元行列には含まれない)。テンプレート:Math から、1個の無次元量があることが分かる。関係式はすなわちこの無次元量が定数ということである。

Tm/k=A(=2π)

減衰振動

ばねにつながれた物体が、速度に比例した大きさの抵抗(粘性抵抗力)を受けながら一次元運動することを考える。運動方程式は以下である[7]テンプレート:Main

mx¨=cx˙kx

式に現れる定数は、物体の質量 テンプレート:Math粘性抵抗の比例係数 テンプレート:Mathばね定数 テンプレート:Math の3つで、それぞれの次元は[m]=𝖬,[c]=𝖬𝖳1,[k]=[𝖬𝖳2]である。

この運動では、特徴的な時間尺度 (characteristic time scale) が2つ存在する。即ち、

  • 減衰時間:τ=mc
  • 固有周期:1ω=mk

の2つの時間が現象を特徴づけており、時間尺度の競合が起こる。つまり テンプレート:Mathテンプレート:Math の大きさのバランスによって運動の様子が変わることが予想される。

Π定理からは、物理量が テンプレート:Math の3つで次元が𝖬,𝖳の2種類である(調和振動のときと同じ理由によって初期変位は入れなくても良い)から、次元行列が

M=(mck𝖬111𝖳012)

となる。したがって1つの無次元量でこの現象を特徴づけられることがわかる。この無次元量には通常、減衰比と呼ばれる

ζ=1/2τω=c/2mk

が用いられ、実際に運動方程式を解析的に解くと、テンプレート:Math のとき減衰振動、テンプレート:Math のとき臨界減衰、テンプレート:Math のとき過減衰となり、運動が定性的にも変化する。

流体機械

ポンプ送風機発電用水車などのターボ機械は内部流れが複雑であるため、その挙動を表すナビエ-ストークス方程式を直接解くことができない。しかしその運転状態は以下の条件を与えるとおおよそ決まることが分かっている:

このとき、次の未知量を推測する:

この場合は物理量は6つ、次元が3種類である。

次元が一致するように各変数のべきを調整すると、(変数が多いので一意ではないが)以下のように関係式を推測できる:

P=AρN2D2(QND3)α
L=BρN3D5(QND3)β

ここで、テンプレート:Math は次元解析から求めることはできないが、条件で考慮していない流体の粘度や機械の各部寸法バランスなどに依存する無次元量である。

この場合の次元行列は

M=(ρDNQPL𝖬100011𝖳001123𝖫310312)

であるため無次元数は テンプレート:Math3つ存在する。よく用いられるのはそれぞれ流量係数、圧力係数、出力係数と呼ばれる以下の3つである:

ϕ=QND3,ψ=PρN2D2,τ=LρN3D5

無次元の関係式 テンプレート:Math で表すと

ψ=f(ϕ),τ=g(ϕ)

となる。

原子構造

原子構造を古典物理学が説明できないということも次元解析から理解できる[8]

水素原子電子クーロン力で惑星のように陽子に束縛されている。その軌道の半径 テンプレート:Math[a]=𝖫)は、

で表されると考えられる。ここで、𝖬 は質量、𝖫 は長さ、𝖳 は時間、𝖨電流の次元を表す。ところが、これらの量をどう組み合わせても、長さの次元 𝖫 を持った量を構成することができない。すなわち、水素原子は一定の大きさをとることができない。そこでニールス・ボーアは、このようなミクロの世界では次元が 𝖬𝖫2𝖳1プランク定数 h が関係していると考えた。以上の4つの物理量を組み合わせて長さの次元を持つ量を作ると、

a=ϵ0h2me2

が導かれる。これはボーア半径テンプレート:Math 倍である。

以上の次元解析的議論により、ボーアは テンプレート:Math が必須であることを確信した。

拡張

次元解析を行う際に用いる次元は国際単位系の基本単位に対応する7つの次元に限る必要はなく、扱う問題に応じて独立した次元を選ぶことができる[9]。たとえば加速度のない流れでは質量、長さ、時間に加えて力を独立次元とみなすことでより厳密な情報が得られるというブリッジマン(1921)に由来する方法がある。

また長さの次元 𝖫 に対して、3方向 テンプレート:Math を区別して次元解析してもよい。この方法はHuntley(1955)に由来し[9]方向性次元解析テンプレート:En[10]またはCosta(1971)によって指向解析 (directional analysis)[9]と呼ばれる。重力境界層など、特別な方向をもつ物理現象に対しては方向性次元解析が有効になる場合がある。

例として、流れの中に、流れに平行に置かれた平板が受ける抗力の問題を考える[9]。抗力 テンプレート:Math、平板の面積 テンプレート:Math、流速 テンプレート:Math、流体の密度 テンプレート:Math、粘性 テンプレート:Math、平板前縁から流れに沿って測った距離を テンプレート:Math とする。独立次元として𝖬𝖫𝖳を用いる通常の次元解析では2つの無次元数:抗力係数 テンプレート:Math とレイノルズ数 テンプレート:Math

f=FSρu2,Re=xuρμ

が得られるが、これらの間に成り立つ関係式の具体形は分からない。しかし平板に平行な2方向 テンプレート:Math の長さの次元𝖫と、平板に直交する テンプレート:Math 方向の長さの次元𝖫zを独立と考えることによって、層流の場合には

fRe1/2=const.

という、より詳細な関係式を得ることができる。

また、Moran(1967)によって群論的方法との関連も論じられている[9]

脚注

テンプレート:Reflist

関連項目

外部リンク

テンプレート:Physics-stub テンプレート:Authority control