ベルヌーイの定理 — トラブルシューティングガイド
より充実した内容を bernoulli-equation.html でご覧いただけます。
トラブルシューティング
ベルヌーイの定理に関連したCFDの問題って、どういうものがありますか？
全圧・静圧に関わるトラブルは実務で非常に多い。代表的なものを整理しよう。
よくある問題と対処法
1. 全圧と静圧の取り違え
症状: 圧力境界条件で流量が想定の2倍以上になる、または全く流れない。
原因: Pressure Inletに静圧を入れてしまう（または逆）。全圧 $p_0$ = 静圧 $p$ + 動圧 $\frac{1}{2}\rho v^2$ であり、例えば流速50 m/sの空気流では動圧は約1530 Paにもなる。
対策: 入口条件が全圧か静圧か、ソルバーのマニュアルで必ず確認する。Fluentの Pressure Inlet は全圧、Pressure Outletは静圧がデフォルトだ。
2. 全圧が増加する異常
CFD結果で下流の全圧が上流より高くなることがあるんですが、これっておかしいですよね？
症状: 流線に沿って全圧が増加している。
考えられる原因:
- メッシュが粗すぎて数値拡散が逆方向の「数値的エネルギー注入」を引き起こしている
- 境界条件の設定ミス（入口と出口の圧力が逆）
- 回転体や外力源がある場合は正当（ポンプ等）
対策: メッシュを2倍に細分化して全圧分布が変わるか確認。粘性流れでは全圧は必ず減少するはずだ。
3. ベンチュリ管のCFD結果が理論値と合わない
| 乖離の程度 | 主な原因 | 対策 |
|---|---|---|
| 5%以内 | 粘性損失（理論では無視） | 正常。$C_d$ で補正 |
| 5〜20% | メッシュ不足、乱流モデル不適切 | 喉部のメッシュ細分化 |
| 20%以上 | 境界条件の誤り、流れの剥離 | 入口の速度分布プロファイル見直し |
4. 圧力損失の過小評価
CFDの圧損が実験値より明らかに小さいことがあるんですけど…
症状: 圧損がDarcy-Weisbachの理論値や実験値の50%以下。
考えられる原因:
- 壁面粗さ（wall roughness）が未設定（デフォルトは滑面）
- y+が壁関数の適用範囲外
- 乱流モデルが不適切（層流で計算してしまっている）
対策: 壁面粗さパラメータを設定（鋼管なら $k_s \approx 0.045$ mm）。y+を30〜300の範囲に調整。Re数に応じた乱流モデルを選択する。
結局、ベルヌーイの理論値との比較が問題発見の第一歩なんですね。
その通り。理論的に予測できる値があるなら、必ずCFD結果と突き合わせること。これがデバッグの最短ルートだ。
「ベルヌーイが合わない」のはほぼ仮定の破れ
CFD結果と手計算のベルヌーイ値がズレると「計算がおかしい」と思いがちですが、たいていの場合「ベルヌーイの仮定が成立していない」のが原因です。実務でよくある落とし穴は①粘性損失を無視できない管内流れ（損失ヘッドの見落とし）、②剥離や再循環が起きている非定常領域にベルヌーイを適用してしまうこと、③密度変化が大きい高速流れへの誤適用です。ベルヌーイは非常に便利な式ですが「理想流体の特定の流線上」でのみ成立します。まずその前提を一つひとつ確認するのが、結果の整合性チェックの基本です。
トラブル解決の考え方
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ベルヌーイの定理の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告