次元解析とπ定理 — トラブルシューティングガイド

単位系の不整合やスケーリングミスは、CFDで最も多い初歩的エラーの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）を使う。

症状: 層流のはずなのに乱流的な振動が出る、または逆に乱流のはずなのに層流解になる

原因: 代表長さや代表速度の定義が不適切。例えば管径 $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$

症状: 球の $C_D$ が文献値から大幅にずれる

考えられる原因:

• Reference Valuesの面積が間違っている（投影面積 $\pi D^2/4$ vs. 表面積 $\pi D^2$）
• 速度の基準が自由流速でない
• 計算領域が狭く、ブロッケージ効果が出ている（推奨: 上流5D以上、下流15D以上、側方5D以上）

$y^+$ は壁面近傍メッシュの品質を示す最重要無次元数だ。

$y^+$ の事前推定式:

平板の経験式 $C_f \approx 0.058\,\text{Re}_L^{-0.2}$ を使えば、解析前に第一層厚さを概算できる。

計算結果を評価する際の「次元解析的」チェック項目をまとめよう。

• 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エンジニアの基本素養だよ。