📚 学習ノート中学・高校・大学の学習内容に沿った自学習サイト

次元解析と相似則——無次元量から何がわかるか

単位を変えても物理法則が変わらないことを出発点に、指数の決定、バッキンガムのΠ定理、模型と実物の相似条件を導く。比例係数や未知の関数、変数の選び漏れなど、次元だけでは決まらない部分も見分ける。

この単元の目標
  • 単位と次元を区別し、式の各項の次元を点検できる。
  • 次元の指数を連立させて尺度を求め、無次元係数が残る理由を説明できる。
  • Π定理の変数の数と次元行列の階数から、独立な無次元量の数を数えられる。
  • レイノルズ数とフルード数の意味から、模型でそろえるべき条件と相似の限界を判断できる。
📑 目次
  1. 1．物差しを取り替えても、現象は変わらない
  2. 2．指数を決める——それでも定数は残る
  3. 3．無次元化は「何倍か」という物差しを作る
  4. 4．Π定理——五つの量を二つの組へ
  5. 5．レイノルズ数は競争する効果の比
  6. 6．模型は形だけでなく、支配する比もそろえる
  7. 7．一つの無次元数が合っても十分でない場合
  8. 8．手元の問題で使う順序

1．物差しを取り替えても、現象は変わらない

机の長さを1 mとも100 cmとも書けますが、机自体が100倍になるわけではありません。物理法則にも、単位の選び方に左右されない骨組みがあるはずです。その骨組みを長さ・時間・質量の種類で点検するのが次元解析です。数値を当てはめる前に「その足し算や比例式は、量の種類として可能か」を調べます。

次元は量の種類、単位はその量を測る基準です。長さの次元をL、時間をT、質量をM、量qの次元を[q]と書きます。ここでTは時間の次元記号であり、振り子の周期は混同を避けてPと表します。mとcmは異なる単位ですが同じ次元Lです。力学の範囲では次のようにM、L、Tの積で表せます。

量次元SI単位
速度UL T⁻¹m/s
加速度gL T⁻²m/s²
力FM L T⁻²N = kg·m/s²
密度ρM L⁻³kg/m³
粘性係数μ（動粘度ではない）M L⁻¹ T⁻¹Pa·s = kg/(m·s)
動粘度ν = μ/ρL² T⁻¹m²/s
[a+b]=[a]=[b],[ab]=[a][b],[a/b]=[a]/[b][a+b]=[a]=[b],\qquad [ab]=[a][b],\qquad [a/b]=[a]/[b]

最初の式は、足し算が物理量として意味をもつための必要条件です。移動距離を「初速度×時間＋加速度×時間」と書くと、第1項は長さ、第2項は速度なので不適切。加速度の項には時間の2乗が必要です。ただし、これだけでは係数が1/2だとはわかりません。数値計算では、同じ次元であってもmとcmを統一してから足します。

次元が一致することは、式が正しいことの必要条件であって十分条件ではありません。エネルギーとトルクはともにM L² T⁻²ですが同じ物理量ではありません。また、指数関数や対数、三角関数の引数は無次元にします。減衰をe^(−t/τ)と書くのは、tと時定数τの比が無次元だからです。角度も弧長/半径なので無次元ですが、角度という意味を失うわけではありません。

2．指数を決める——それでも定数は残る

長いブランコがゆっくり揺れることを、運動方程式を解く前に予想できるでしょうか。固定支点、質点のおもり、一様な重力、摩擦なしを仮定し、おもりを支点から固定距離ℓに保つ質量のない剛棒で拘束した理想振り子を考えます。鉛直下向きからの角度θ₀で静止から解放します。糸ではθ₀ > π/2で解放した直後に弛み、固定距離の拘束が保たれません。剛棒は押す力も伝えられるため、この問題を避けられます。長さℓ、重力加速度g、質量mが周期Pを決めると仮定し、まずPをこれらのべきの積で作ります。振幅の影響はこの後で点検します。

🔍 周期の次元Tを作る

  1. べきの積を仮定する
    P=Cℓagbmc,[C]=1P=C\ell^a g^b m^c,\qquad [C]=1

    Cは無次元係数。ℓ、g、mはいずれも正とする。

  2. 次元の指数を比較する
    T=La+bT−2bMc⟹a+b=0, −2b=1, c=0T=L^{a+b}T^{-2b}M^c\quad\Longrightarrow\quad a+b=0,\ -2b=1,\ c=0

    長さと質量の指数は0、時間の指数は1でなくてはならない。

  3. 許される尺度を得る
    a=12,b=−12,c=0,P=Cℓga=\frac12,\quad b=-\frac12,\quad c=0,\qquad P=C\sqrt{\frac{\ell}{g}}

    この変数集合では質量依存は残らず、長さを4倍にすると周期は2倍になる。ただしCが同じ場合の比較である。

無次元の変数は、次元の指数比較には見えません。振り子には初期角度θ₀もあります。剛棒拘束で静止から解放する往復振動の範囲0 < θ₀ < π（rad）なら、正しくはCがθ₀の関数になりえます。次元解析だけでは「振幅に無関係」も「C = 2π」も証明できません。棒の質量や空気抵抗、支点の摩擦が重要なら、変数集合そのものを見直します。
P=ℓg f(θ0),θ0≪1 ⟹ f(θ0)≈2πP=\sqrt{\frac{\ell}{g}}\,f(\theta_0),\qquad \theta_0\ll1\ \Longrightarrow\ f(\theta_0)\approx2\pi

小角近似では運動方程式を線形化すると角振動数が√(g/ℓ)となり、1往復は位相2πなので係数2πが得られます。これは次元解析に加えた力学的情報です。角度を変えず長さだけ変える比較なら、未知のfが比の中で消えます。この「係数を知らなくても倍率がわかる」が相似則の強みです。

例題1：絶対値より先に倍率を予測する

同じ重力加速度g = 9.8 m/s²の場所で、質量のない剛棒で拘束した理想振り子Aの長さは0.25 m、Bは1.00 m。両者を同じ小さい初期角度で静止から解放する。周期比を求め、小角近似で各周期も計算せよ。どの結果が次元解析だけで得られるか。

解答・解説

比はP_B/P_A = √(1.00/0.25) = 2。係数が共通なので、周期比2は次元解析と同じ初期角度という条件から得られます。小角近似の2πを追加するとP_A ≈ 2π√(0.25/9.8) ≈ 1.00 s、P_B ≈ 2π√(1.00/9.8) ≈ 2.01 s。丸める前の値は約1.00354 sと2.00709 sです。絶対値には運動方程式から得た係数が必要です。

3．無次元化は「何倍か」という物差しを作る

模型の1秒と実物の1秒をそのまま比べても、同じ運動段階とは限りません。例えば振り子なら自然な時間尺度τ = √(ℓ/g)で時刻を割り、t* = t/τとすれば「その系の時間尺度の何倍が経ったか」を比較できます。無次元量は単位を消すだけの細工ではなく、基準に対する相対的な大きさを表します。基準の選び方と物理的意味をセットで述べることが重要です。

t∗=tℓ/g,d2θdt∗2+sin⁡θ=0t^*=\frac{t}{\sqrt{\ell/g}},\qquad \frac{d^2\theta}{d{t^*}^2}+\sin\theta=0

剛棒拘束の理想振り子の運動方程式をこの時刻で書けば、ℓとgは式から消えます。初期角度と無次元の初期角速度まで同じなら、θをt*に対して描く曲線は同じです。一方、角度が違えば無次元化後の初期条件が違うので別の曲線です。「式から単位が消えた」だけでは相似は完成しません。

4．Π定理——五つの量を二つの組へ

物体が流体中で受ける抗力の大きさFを考えます。形状を固定し、代表長さℓ、流速U、密度ρ、粘性係数μのほかは影響しないと仮定します。非圧縮のニュートン流体、一定の物性値を想定し、自由表面・圧縮性・壁の影響などは除いたモデルです。乱れを伴う場合のFは十分に平均した抗力とします。五つの量を五本の軸で整理する代わりに、単位を消した少数の組でまとめます。

バッキンガムのΠ定理は、単位の取り方に依存しない関係にn個の変数があり、その次元行列の階数がrなら、n−r個の独立な無次元の積Πで関係を書けるというものです。次元行列とは、各変数のM、L、Tの指数を縦に並べた行列です。rは独立な次元の方向の数であり、常に機械的に3とするのではありません。無次元の入力変数もnに数えます。

🔍 抗力に使う無次元量を構成する

  1. 繰り返し変数を選ぶ
    [ρ]=ML−3,[U]=LT−1,[ℓ]=L[\rho]=ML^{-3},\quad [U]=LT^{-1},\quad [\ell]=L

    ρ、U、ℓの次元は独立でM、L、Tを張る。すべて正とする。五変数の階数は3なので独立なΠは二つ。

  2. 力を無次元化する
    [FρaUbℓc]=M1+aL1−3a+b+cT−2−b=1⟹ΠF=FρU2ℓ2[F\rho^a U^b\ell^c]=M^{1+a}L^{1-3a+b+c}T^{-2-b}=1\quad\Longrightarrow\quad \Pi_F=\frac{F}{\rho U^2\ell^2}

    指数をすべて0にするとa = −1、b = −2、c = −2。分母は力の尺度で、単位はN。

  3. 粘性係数を無次元化する
    Πμ=μρUℓ=1Re,Re=ρUℓμ=Uℓν\Pi_\mu=\frac{\mu}{\rho U\ell}=\frac1{\mathrm{Re}},\qquad \mathrm{Re}=\frac{\rho U\ell}{\mu}=\frac{U\ell}{\nu}

    μの指数を同様に消すとこの比になる。逆数も無次元なので、慣例的にレイノルズ数Reを用いる。

  4. 未知の関数にまとめる
    ΠF=Φ(Re),F=ρU2ℓ2Φ(Re)\Pi_F=\Phi(\mathrm{Re}),\qquad F=\rho U^2\ell^2\Phi(\mathrm{Re})

    Fが他の変数で一意に定まるモデルを仮定した表示。Π定理は関数Φの形や数値を決めない。

ここでのΠ_FはF/(ρU²ℓ²)という定義です。よく使う抗力係数C_D = F/(ρU²A/2)とは規格化が異なり、Aは基準面積です。係数を比較する前に分母を確認しましょう。また、Fの前にU²が見えても「抗力は常に速度の2乗に比例する」とは言えません。Uを変えるとReも変わり、Φも変わりえます。

5．レイノルズ数は競争する効果の比

流体の慣性は運動を保とうとし、粘性は隣り合う部分の速度差をならします。代表長さℓにわたり速度がU程度変わると考えると、単位体積当たりの慣性項の尺度はρU²/ℓ、粘性項の尺度はμU/ℓ²。その比がReです。粘性係数μの次元も「粘性によるせん断応力 = μ×速度勾配」から確かめられます。

ρU2/ℓμU/ℓ2=Re\frac{\rho U^2/\ell}{\mu U/\ell^2}=\mathrm{Re}

Reが小さければこの尺度比較では粘性が強く、Reが大きければ慣性が強いと読めます。ただし大きなReでも壁の近くなどでは局所的に粘性が重要です。どのReで乱流になるかも形状や流入状態に依存します。すべての流れに共通の一つの境界値があるわけではありません。

6．模型は形だけでなく、支配する比もそろえる

写真を縮小すれば見た目は似ますが、水中の模型を同じ速度で動かしても流れまで似るとは限りません。幾何学的相似は、全長・幅・曲率半径など対応する長さの比が一定であること。動的相似には、それに加えて支配する無次元量と無次元の境界・初期条件をそろえます。以下は前節の抗力モデル内での理想的な比較です。実際に測る場合は表面の相対粗さや周囲の壁までの相対距離も点検します。

例題2：1/4模型でレイノルズ数を合わせる

幾何学的に相似な実物pと模型mを同じ流体（ρ = 1000 kg/m³、ν = 1.0×10⁻⁶ m²/s）中で比較する。実物はℓ_p = 4.0 m、U_p = 0.50 m/s、平均抗力F_p = 40 Nという架空の条件。模型はℓ_m = 1.0 m。前節のモデルでReを一致させる模型の速度と、模型の抗力を求めよ。

解答・解説

同じνなのでU_mℓ_m = U_pℓ_p。したがってU_m = 2.0 m/sで、両者ともRe = 2.0×10⁶。Reと形状などが同じならΠ_Fが同じなので、F_m/F_p = (U_m/U_p)²(ℓ_m/ℓ_p)² = 4²×(1/4)² = 1、F_m = 40 N。実物の力の尺度ρU_p²ℓ_p²は4000 N、模型も4000 Nで、Π_Fはともに0.010。小さい模型だから抗力も小さい、とは限りません。

例題2で模型を実物と同じ0.50 m/sにすれば、Reは実物の1/4です。Π_Fが変わりうるため、長さ比の2乗だけでF_m = 40/16 Nと断定することはできません。あるReの範囲でΦがほぼ一定という追加情報がある場合に限り、そのような近似が可能になります。

7．一つの無次元数が合っても十分でない場合

水面の波が重要な模型では、重力gを変数から除けません。慣性項の尺度ρU²/ℓと重力項ρgの比はU²/(gℓ)です。その平方根であるフルード数Frは、「流れの速さ」と「重力と長さから作る速度尺度」の比になります。次元の階数は3のまま変数が一つ増え、独立な無次元量も一つ増えます。

Fr=Ugℓ,ΠF=Ψ(Re,Fr)\mathrm{Fr}=\frac{U}{\sqrt{g\ell}},\qquad \Pi_F=\Psi(\mathrm{Re},\mathrm{Fr})

同じ流体・同じ重力で長さ比λ = ℓ_m/ℓ_pとすると、Re一致には速度比U_m/U_p = 1/λ、Fr一致にはU_m/U_p = √λが必要です。両方を同時に満たす正のλは1だけ。したがって縮小模型で両者を厳密に合わせることはできません。どちらを重視し、無視した効果がどれほど残るかを別に評価します。表面張力や圧縮性まで重要なら、さらに無次元量が必要です。

「無次元数が同じなら何でも同じ」は誤りです。結論は選んだ支配方程式、変数集合、形状、無次元の初期・境界条件が同じ範囲に限られます。次元解析は、未知の係数・関数の形、解の安定性や一意性、モデルが現実に十分かどうかまでは決めません。相似が崩れたら、計算ミスだけでなく省いた変数も疑います。

8．手元の問題で使う順序

  1. 求めたい量と影響する入力を列挙し、除く効果とその理由を書く。角度・形状比などの無次元入力も忘れない。
  2. 各量の次元をそろえて書き、指数の連立またはΠ定理で尺度と無次元量を作る。
  3. 残った係数・関数のうち、法則の解や観測によって決める部分を明示する。
  4. 比較する二つの系で無次元条件が一致するか調べてから、元の単位へ戻す。極限・単位換算・数値比でも検算する。
練習問題：尺度を作り、条件を点検する

  1. 拡散係数Dの次元はL² T⁻¹である。長さℓとDだけから時間尺度τを、無次元係数を1として作れ。ℓ = 0.020 m、D = 1.0×10⁻⁶ m²/sで計算せよ。この時間が必ず「拡散が完了する時刻」なのか。

    ヒント

    ℓの2乗をDで割ると何の単位が残るか。

    解答・解説

    τ = ℓ²/D = 0.000400/0.000001 = 400 s。これは代表的な時間尺度であって、完了時刻ではありません。具体的な到達割合を満たす時刻は形状・境界条件・初期分布・到達の定義で変わり、次元解析だけでは決まりません。

  2. 前節の抗力モデルで、ρ = 1000 kg/m³、U = 0.20 m/s、ℓ = 0.10 m、μ = 0.010 Pa·s、F = 0.12 Nとする。ReとΠ_Fを求めよ。

    解答・解説

    Re = 1000×0.20×0.10/0.010 = 2000。力の尺度はρU²ℓ² = 1000×0.040×0.010 = 0.400 NなのでΠ_F = 0.12/0.400 = 0.30。どちらも無次元であり、後者はこのReにおけるΦの一つの値を与えるだけです。

  3. 例題2の1/4模型で水面の波を重視し、同じgの下でFrを合わせることにした。必要な模型速度と、そのときの模型Reを求めよ。

    解答・解説

    速度比は√(1/4) = 1/2なのでU_m = 0.25 m/s。Re_m = 0.25×1.0/(1.0×10⁻⁶) = 2.5×10⁵で、実物の2.0×10⁶の1/8。Frは一致しますがReは一致しないため、例題2の抗力40 Nをそのまま使うことはできません。

  4. 質量のない剛棒で拘束した理想振り子を、0 < θ₀ < πの角度で静止から解放する。周期P、長さℓ、g、質量m、初期角度θ₀の五変数を扱い、次元行列の階数は3とする。独立な無次元量を二つ挙げ、次元解析だけで周期が初期角度によらないと言えるか答えよ。

    解答・解説

    五変数−階数3 = 二つで、P√(g/ℓ)とθ₀を選べます。関係はP√(g/ℓ) = f(θ₀)となるので、角度に依存しないとは言えません。質量mは他の変数で質量次元を打ち消せないため、この変数集合では周期に影響する無次元の組に入りません。

確認クイズ：わかることと決まらないこと

問1. 式の左右で次元が一致していた。この事実だけから言えるのはどれか。

問2. n = 6個の変数の次元行列の階数がr = 2だった。Π定理で作る独立な無次元量はいくつか。

問3. F = ρU²ℓ²Φ(Re)で、ρ、ℓ、μを固定しUを2倍にした。次元解析だけで抗力比はどうなるか。

問4. 同じ流体・同じ重力で、全長を1/4にした相似模型のReとFrを厳密に合わせたい。正しい説明はどれか。

根拠と範囲：次元の指数比較とΠ定理、支点の周りの慣性モーメントがmℓ²である質量のない剛棒拘束の理想振り子の運動方程式、非圧縮ニュートン流体の慣性・粘性・重力項の尺度比較に基づく教材です。数値は独立計算用の架空条件であり、特定の船・装置の実測性能や設計値ではありません。抗力の具体的な関数ΦやΨを実測なしに与えたものではなく、模型の結果を実物へ外挿するには採用したモデルの妥当性の確認が別途必要です。
📝 メモこの端末のみに保存
メモ欄はページ読み込み後に使えるようになります。