次元解析とπ定理 — トラブルシューティングガイド
より充実した内容を dimensional-analysis.html でご覧いただけます。
単位系・スケーリングに起因するトラブル
次元解析がらみで実際に起きるトラブルって、どんなものがありますか？
単位系の不整合やスケーリングミスは、CFDで最も多い初歩的エラーの1つだ。具体的なケースを見ていこう。
1. 圧力値が桁違いにおかしい
症状: 大気圧条件なのに圧力が $10^9$ Pa のオーダーになる
原因: CAD形状がmm単位なのに、物性値をSI（m基準）で入力した。長さが1000倍になるので、$\Delta p \propto \rho U^2$ から圧力が $10^6$ 倍になる。
対策: Fluent の場合、Mesh > Scale で形状をm単位にスケーリングするか、整合する単位系（密度: tonne/mm$^3$、粘度: MPa$\cdot$s）を使う。
2. Re数が想定と合わない
症状: 層流のはずなのに乱流的な振動が出る、または逆に乱流のはずなのに層流解になる
原因: 代表長さや代表速度の定義が不適切。例えば管径 $D$ を使うべきところで半径 $r$ を使っている（Re数が2倍ずれる）。
対策:
- $\text{Re}$ を手計算で確認: $\text{Re} = \rho U D / \mu$
- 管内流の遷移: $\text{Re}_D \approx 2300$
- 外部流（平板）: $\text{Re}_x \approx 5 \times 10^5$
Re数を自分で計算して検算するのが大事なんですね。
3. 抗力係数 $C_D$ が文献値と合わない
症状: 球の $C_D$ が文献値から大幅にずれる
考えられる原因:
- Reference Valuesの面積が間違っている（投影面積 $\pi D^2/4$ vs. 表面積 $\pi D^2$）
- 速度の基準が自由流速でない
- 計算領域が狭く、ブロッケージ効果が出ている（推奨: 上流5D以上、下流15D以上、側方5D以上）
面積の定義が文献と合っているか確認するのがポイントですね。
$y^+$ 関連のトラブル
$y^+$ は壁面近傍メッシュの品質を示す最重要無次元数だ。
| 乱流モデル | 必要な $y^+$ | メッシュ戦略 |
|---|---|---|
| 標準壁関数 (Standard Wall Function) | 30 - 300 | 壁面から離して第一層配置 |
| Enhanced Wall Treatment | $\approx 1$ | 壁面に密着したプリズム層 |
| $k$-$\omega$ SST（壁面解像） | $\approx 1$ | 10層以上のプリズム層推奨 |
| SA (Spalart-Allmaras) | $\approx 1$ | 境界層内に十分な層数 |
$y^+$ の事前推定式:
平板の経験式 $C_f \approx 0.058\,\text{Re}_L^{-0.2}$ を使えば、解析前に第一層厚さを概算できる。
計算前に$y^+$を予測してメッシュ設計するんですね。
スケーリングミスの検出チェックリスト
計算結果を評価する際の「次元解析的」チェック項目をまとめよう。
- Re数を手計算で算出し、流れレジーム（層流/遷移/乱流）が想定通りか確認
- $C_D$, $C_L$ がオーダー的に妥当か（鈍頭物体の $C_D \sim O(1)$、流線形は $O(10^{-2})$）
- 圧力損失が $\Delta p \sim \frac{1}{2}\rho U^2$ のオーダーに収まるか
- Ma数が0.3未満であることを確認（非圧縮性仮定を使う場合）
- 温度変化がある場合、Gr/Re$^2$ の大きさで強制/自然対流の支配性を判断
次元解析の知識があると、結果の妥当性を素早く判断できるんですね。
その通り。計算結果を鵜呑みにせず、無次元数のオーダー推定で「おかしい」と気づけることが、CFDエンジニアの基本素養だよ。
Re数が合っても「同じ流れ」にならないことがある
「模型実験でRe数を実機と合わせれば相似則が成立する」——これは多くの場合正しいですが、落とし穴があります。例えば圧縮性効果（マッハ数）、表面張力（ウェーバー数）、重力（フルード数）が同時に重要な問題では、すべての無次元数を同時に一致させることは物理的に不可能です。風洞実験で高Re数を出すために高圧空気を使うと、今度はRe数は合うがマッハ数がズレる——という딜레마が起きます。次元解析のトラブルの本質は「どの無次元数を優先するか」という物理的な判断の問題です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——次元解析とπ定理の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告