ベルヌーイの定理 — トラブルシューティングガイド
トラブルシューティング
ベルヌーイの定理に関連した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結果と突き合わせること。これがデバッグの最短ルートだ。
F1と空力の戦い
F1マシンは時速300kmで走ると、車重と同じくらいのダウンフォース（下向きの空力的な力）を発生します。つまり理論上、天井に貼り付けて走れる！ チームは数千CPU時間のCFDシミュレーションを毎週実行し、フロントウィングの角度を0.1°単位で最適化しています。F1はCAEの技術力がそのまま順位に直結する世界です。
トラブル解決の考え方
デバッグのイメージ
CFDのデバッグは「水道管の詰まり修理」に似ている。まず「どこで詰まっているか」（どの残差が下がらないか）を特定し、次に「何が詰まっているか」（メッシュ品質？境界条件？乱流モデル？）を調べ、最後に「どう直すか」（メッシュ修正？緩和係数？）を判断する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ベルヌーイの定理の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
ベルヌーイの定理の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →