バッキンガムπ定理シミュレーター — 次元解析と模型相似則

1. 実機のサイズ・流速を決める（D_p, U_p）：まず予測したい実機の直径と流速を設定します。既定値は D_p = 1m, U_p = 10m/s。これは「水中の直径1mの円柱に10m/sで水が当たる」状況に相当します。

2. 模型のサイズを決める（D_m）：試験に使う模型の直径を入力します。既定値は D_m = 0.10m（実機の1/10）。スケール比 λ = D_m/D_p が計算結果に表示されます。

3. 流体を選ぶ（ν_m/ν_p）：模型試験で使う流体の動粘度比です。同流体なら 1.0、空気→水なら約 0.07 など。

4. 結果を読む：「相似則による模型流速 U_m」が、Re を実機に一致させるために必要な模型試験の流速です。「力比 F_m/F_p」は、模型に働く力が実機の何倍になるかを示します。

ステップ1：効く変数を選ぶ
円柱周りの抗力に効くのは：抗力 F、直径 D、流速 U、流体の密度 ρ、流体の粘度 μ。表面粗さや流れの乱れなどは今は無視します。

ステップ2：基本次元で書く
F [MLT⁻²]、D [L]、U [LT⁻¹]、ρ [ML⁻³]、μ [ML⁻¹T⁻¹]。基本次元は M, L, T の3つ。

ステップ3：π群の数を求める
バッキンガムπ定理より、独立な無次元群の数は $n - r = 5 - 3 = 2$ 個。

ステップ4：π群を構成する
「繰り返し変数」として D, U, ρ を選び、残りの F と μ を無次元化します。
$\pi_1 = \dfrac{\mu}{\rho U D} = \dfrac{1}{Re}$ → 慣例的に逆数を取って $Re = \dfrac{\rho U D}{\mu}$
$\pi_2 = \dfrac{F}{\rho U^2 D^2} \propto C_D$（標準的な $C_D$ は分母に $\tfrac{1}{2}\rho U^2 A$ を使うが、形状係数を除けば本質は同じ）

ステップ5：関係式に落とす
π定理の結論として、抗力係数はレイノルズ数のみの関数：$C_D = f(Re)$。実験で1本の曲線を求めれば、あらゆる D, U, ρ, μ の組み合わせで抗力が予測できます。

ステップ6：模型相似則を導く
$Re_m = Re_p$ のとき $C_D$ も等しい。同流体（$\rho_m = \rho_p$）なら：

$\dfrac{U_m}{U_p} = \dfrac{D_p}{D_m}\cdot\dfrac{\nu_m}{\nu_p}, \qquad \dfrac{F_m}{F_p} = \left(\dfrac{U_m}{U_p}\right)^2 \left(\dfrac{D_m}{D_p}\right)^2$

既定値での手計算例：D_p=1.0m, U_p=10m/s, D_m=0.10m, ν_m/ν_p=1.0
スケール比 $\lambda = D_m/D_p = 0.100$
$U_m = U_p \cdot (D_p/D_m) \cdot (\nu_m/\nu_p) = 10 \times 10 \times 1 = 100$ m/s
水 20°C で $\nu_p = 1.0\times10^{-6}$ m²/s と仮定すると：$Re = U_p D_p / \nu_p = 10 \times 1 / 10^{-6} = 1.0\times10^7$
模型側でも $Re_m = U_m D_m / \nu_m = 100 \times 0.1 / 10^{-6} = 1.0\times10^7$ で一致 ✓
力比 $F_m/F_p = (100/10)^2 \times (0.1/1)^2 = 100 \times 0.01 = 1.000$

風洞試験による航空機・自動車設計：実機の航空機や自動車を直接試験するのは費用も時間もかかりすぎます。バッキンガムπ定理に基づく相似則を使えば、1/10〜1/50スケールの模型を風洞に入れ、Re や Ma（マッハ数）を合わせて試験することで、実機の抗力・揚力・モーメントを高精度に予測できます。F1チームが風洞模型で数百分の1秒を削るのも、この原理に支えられています。

船舶模型試験（曳航水槽）：船の造波抵抗にはフルード数 Fr が支配的に効きます。模型船を曳航水槽で動かし、Fr を実船と一致させることで、巨大なタンカーや空母の抵抗を実用的なサイズの模型で評価します。粘性抵抗（Re 支配）は別途補正する「フルード分割」が標準的手法です。

化学プラントのスケールアップ：実験室のフラスコ規模で開発した反応・撹拌プロセスを、商業生産規模のリアクター（数百〜数千倍の体積）に拡大する際、相似則が不可欠です。撹拌レイノルズ数、ヌッセルト数、ダムケラー数などの π 群を整理し、どれを優先するかで「幾何相似」「動的相似」「熱的相似」を選択します。

河川・港湾水理模型：洪水時の流れ、津波の遡上、港湾内の波浪——これらを数値計算だけでなく実物模型でも検証します。地形を縮尺し、フルード相似で流速と時間をスケーリングし、堤防の越流や防波堤の効果を実証的に評価します。

最も多い誤解は、「模型を小さくすれば模型試験は楽になる」と考えてしまうことです。実際にはちょうど逆で、Re を一致させるには模型流速を上げる必要があり、$U_m = U_p \cdot (D_p/D_m)$ なので、1/10 模型なら流速は 10 倍。1/100 模型なら 100 倍。シミュレーターで D_m スライダーを最小にしてみてください。U_m がとんでもない値になります。実機が水中 10 m/s なら、1/100 模型では水中 1000 m/s——もはや実験不可能です。だから流体を変える（ν を下げる）、圧力を上げる（高圧風洞）、別の支配パラメータを優先する、といった工夫が必要になります。

次に多いのが、「複数の無次元数を同時に一致させられる」と思い込むことです。船の自由表面流れでは Re と Fr の両方が効きますが、同じ流体で両者を同時に満たすには $\nu_m/\nu_p = (D_m/D_p)^{3/2}$ という現実離れした流体が必要で、実質的に不可能です。航空機の遷音速試験では Re と Ma が同時に効きます。実務では「支配的な π を1つ選んで合わせ、残りは補正する」という割り切りが常識です。バッキンガムπ定理が教えてくれるのは「変数を減らせる」ことであって「全部を魔法のように再現できる」ことではない、と覚えてください。

最後に、「次元解析だけで物理が決まる」と勘違いするミス。π定理は $C_D = f(Re)$ という関数形を保証しますが、$f$ の中身は教えてくれません。$C_D$ が Re とともにどう変化するか——層流域での Stokes 域、Re ≈ 10⁵ 付近の急減（抗力危機）、超臨界域でのほぼ一定値——これらは実験や CFD で曲線を実際に求める必要があります。シミュレーターは「相似条件」と「スケーリング」の骨格を示しますが、円柱抗力の数値そのものは出力していません。π定理は整理の道具であって、力学法則そのものを置き換えるものではない、と理解することが正しい使い方の出発点です。