次元解析

出典: フリー百科事典『ウィキペディア(Wikipedia)』
移動: 案内検索

次元解析(じげんかいせき、: dimensional analysis)とは、物理量における、長さ質量時間電荷などの次元から、複数の物理量の間の関係を予測することである。

物理的な関係を表す数式においては、両辺の次元が一致しなくてはならない。この規則を逆に利用すると、既知の量を組み合わせ、求めたい未知の物理量の次元に一致するように式を立てれば、それは正しい関係式になっている可能性が高い。

次元解析を用いると、一般解を得ることが困難な(ときには不可能な)現象に対して、物理量間の関係を推測することができる。

次元一致の原理[編集]

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

数学的表現[編集]

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

Q=F(x_1,\dots,x_n)

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

F(x_1,\dots,x_n) = \prod_i [X_i]^{a_i} \times F^*(x_1^*,\dots,x_n^*)

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

バッキンガムのΠ定理[編集]

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

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

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

定式化[編集]

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

f(q_1, q_2, \dots, q_n) = 0

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

F(\pi_1, \pi_2, \dots, \pi_p) = 0

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

\pi_i = q_1^{a_1} \, q_2^{a_2} \cdots q_n^{a_n}, \quad i=1,\dots,p

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

証明[編集]

概要[編集]

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

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

正式な証明[編集]

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

M = \begin{pmatrix} a_1 \\ \vdots \\ a_n \end{pmatrix}

は物理変数

q_1^{a_1} \, q_2^{a_2} \cdots q_n^{a_n}

の次元行列である。

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

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

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

[編集]

調和振動[編集]

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

式に含まれるであろう定数は、物体の質量m [\mathsf{M}]、ばね定数k [\mathsf{M T}^{-2}]、初期変位x [\mathsf{L}]の3つで、求めるべき周期T次元\mathsf{T}である(長さの次元を\mathsf{L}、質量の次元を\mathsf{M}、時間の次元を\mathsf{T}とした)。この中で長さの次元を含んでいるのは初期変位のみなので、式に含めることができない。なぜなら式の左辺と右辺では次元が一致しなくてはならず、初期変位を含めるならば両辺に同じだけかける必要があり、それならば無くても同じだからである。

次元が\mathsf{T}になるようにmk を組み合わせる方法は一つしかない。結果次の式が求まる。

T = A \sqrt{\frac{m}{k}}

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

T=2\pi\sqrt{\frac{m}{k}}

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

バッキンガムのΠ定理にしたがって考えると、物理量がmkx およびT の4つで、次元が\mathsf{M}, \mathsf{T}, \mathsf{L}の3種類なので、次元行列は

M = \begin{pmatrix} 1 & 1 & 0 & 0  \\  0 & -2 & 0 & 1  \\  0 & 0 & 1 & 0 \end{pmatrix}

となる(列がmkxT に、行が\mathsf{M}, \mathsf{T}, \mathsf{L}に対応する)。null M = 1から、1個の無次元量があることが分かる。関係式はすなわちこの無次元量が定数ということである。

\frac{T}{\sqrt{m/k}} = A (=2\pi)

減衰振動[編集]

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

 m \ddot{x} = -\gamma \dot{x} - k x

式に現れる定数は、物体の質量 m [\mathsf{M}]粘性抵抗の比例係数γ [\mathsf{M T}^{-1}]ばね定数 k [\mathsf{M T}^{-2}]の3つである。

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

  • 減衰時間: \tau = \frac{m}{\gamma}
  • 固有周期: \frac{1}{\omega} = \sqrt{\frac{m}{k}}

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

Π定理からも、物理量が m 、γ、k の3つで、次元が\mathsf{M}, \mathsf{T}の2種類である(調和振動のときと同じ理由によって、初期変位は入れなくても良い)から、次元行列が

M = \begin{pmatrix} 1 & 1 & 0  \\  0 & -1 & 1 \end{pmatrix}

となる。したがって1つの無次元量τωの大きさで振動を定性的に分類できることがわかる。

実際、運動方程式を解析的に解くと、τω > 1/2のとき減衰振動、τω = 1/2のとき臨界減衰、τω < 1/2のとき過減衰となり、運動が定性的にも変化する。

流体機械[編集]

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

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

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

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

P = A \rho N^2 D^2 \left(\frac{Q}{ND^3}\right)^\alpha
L = B \rho N^3 D^5 \left(\frac{Q}{ND^3}\right)^\beta

ここで、AB、α、βは次元解析から求めることはできないが、条件で考慮していない流体の粘度や機械の各部寸法バランスなどに依存する無次元量である。この場合の次元行列は

M = \begin{pmatrix} 1&0&0&0&1&1  \\  0&0&-1&-1&-2&-3  \\  -3&1&0&3&-1&2 \end{pmatrix}

で無次元数は null M = 3つ存在する:

\phi=\frac{Q}{ND^3},\quad\psi=\frac{P}{\rho N^2 D^2},\quad\tau=\frac{L}{\rho N^3 D^5}

無次元の関係式 fg で表すと

\psi = f(\phi),\quad\tau = g(\phi)

となる。

原子構造[編集]

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

水素原子電子クーロン力で惑星のように陽子に束縛されている。その軌道の半径 [\mathsf{L}] は、

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

a = \frac{\epsilon_0 h^2}{m e^2}

が導かれる。これはボーア半径のπ倍である。

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

方向性次元解析[編集]

長さの次元 \mathsf{L} に対して、3方向 (x , y , z ) を区別して次元解析することを方向性次元解析vectorial dimensional analysis)という[7]重力境界層など、特別な方向をもつ物理現象に対しては方向性次元解析が有効になる場合がある。

参考文献[編集]

  1. ^ 化学工学会編 『化学工学』 (3版) 槇書店、2006年、6頁。ISBN 4-8375-0690-9 
  2. ^ 大野克嗣 『非線形な世界』 東京大学出版会、2009年ISBN 978-4-13-063352-9 
  3. ^ 五十嵐保; 杉山均 『流体工学と伝熱工学のための次元解析活用法』 共立出版、2013年、6頁。ISBN 978-4-320-07189-6 
  4. ^ 白樫正高 (1994). 次元解析再考. http://lib.nagaokaut.ac.jp/kiyou/data/study/k16/K16_14.pdf. 
  5. ^ 山本鎮男・曽根彰・芦野隆一・守本晃 『ダイナミカルシステムの数理 基礎』 共立出版。
  6. ^ 大野克嗣 『非線形な世界』 東京大学出版会、2009年、165頁。ISBN 978-4-13-063352-9
  7. ^ 五十嵐保; 杉山均 『流体工学と伝熱工学のための次元解析活用法』 共立出版、2013年、104頁。ISBN 978-4-320-07189-6 

関連項目[編集]

外部リンク[編集]