ベルヌーイの定理
ベルヌーイの定理の理論基礎
概要
先生、ベルヌーイの定理ってよく聞くんですけど、何がそんなに大事なんですか？
定常・非粘性・非圧縮の流れにおけるエネルギー保存則だ。流体力学で最も基本的かつ強力なツールの一つだよ。流速と圧力の関係を直接結びつけるから、配管設計やピトー管、ベンチュリ管など至るところで使われる。
支配方程式の導出
どうやって導出するんですか？
オイラー方程式（非粘性流体の運動方程式）の流線に沿った成分を考える。定常流で流線に沿って積分すると、ベルヌーイの式が得られる。
オイラー方程式は次の通りだ。
定常かつ流線に沿って積分すると、ベルヌーイの式になる。
各項はそれぞれ何を表してるんですか？
第1項 $p$ が静圧（static pressure）、第2項 $\frac{1}{2}\rho v^2$ が動圧（dynamic pressure）、第3項 $\rho gz$ が位置圧（hydrostatic pressure）だ。三者の和を全圧（total pressure）$p_0$ と呼ぶ。
水頭の形で書くこともありますよね？
その通り。両辺を $\rho g$ で割ると水頭形式になる。
ここで $\frac{p}{\rho g}$ が圧力水頭、$\frac{v^2}{2g}$ が速度水頭、$z$ が位置水頭、$H$ が全水頭だ。
適用条件と制約
どんな流れでも使えるんですか？
いや、厳密な適用条件がある。
| 条件 | 意味 | 違反した場合 |
|---|---|---|
| 定常流 | $\partial/\partial t = 0$ | 非定常ベルヌーイ式が必要 |
| 非粘性 | $\mu = 0$（摩擦なし） | 損失水頭項を追加 |
| 非圧縮 | $\rho = \text{const}$ | 圧縮性ベルヌーイ式が必要 |
| 同一流線上 | エネルギー保存は流線ごと | 渦なし流なら全空間で成立 |
粘性損失がある場合は拡張形を使う。
ここで $h_L$ は損失水頭（摩擦やバルブ、曲がりなどの圧力損失）だ。
なるほど、実務では損失水頭をいかに正確に見積もるかが鍵なんですね。
その通り。管路設計ではDarcy-Weisbachの式 $h_L = f \frac{L}{D}\frac{v^2}{2g}$ と組み合わせるのが定石だ。
「翼が揚力を生む理由」の教科書はたぶん間違っている
「上面が長いので空気が速く流れ圧力が下がる」——中学理科や入門書でよく見る翼の揚力説明ですが、これには大きな穴があります。「上側の空気と下側の空気は同時に後縁に到達しなければならない」という前提が実は根拠なしです。実際には上下の空気の「等時刻到達」は起きておらず、上面の流速はその前提より大幅に速い。本当の揚力は循環（Circulation）の理論で説明されます。ベルヌーイの定理は「速いほど圧力が低い」という道具として正しいですが、なぜ速くなるかの説明が教科書でサボられているのです。
ベルヌーイの定理の数値計算手法
ベルヌーイ式の数値的活用
CFDでベルヌーイの式ってどう使うんですか？ 直接解くわけじゃないですよね？
いい質問だ。CFDではNavier-Stokes方程式を解くのが基本だから、ベルヌーイ式を直接離散化して解くわけではない。ただし、境界条件の設定、結果の検証、1Dネットワーク解析ではベルヌーイが核心的な役割を果たす。
境界条件への応用
境界条件としてはどう使うんですか？
例えば、圧力入口（pressure inlet）と圧力出口（pressure outlet）を設定するとき、全圧と静圧の関係が必要になる。
Ansys Fluentでpressure-inletを設定する場合、Total Pressure（全圧）を指定するが、内部でこの式を使って速度場と整合性を取っている。全圧と静圧を混同すると、流量が想定と大きくずれるので注意が必要だ。
1Dネットワーク解析
3D CFDでなく、1Dで済む場合もあるんですか？
配管ネットワークの予備検討ではベルヌーイ式＋損失項の1D解析が非常に有効だ。節点で圧力バランス、枝管で流量バランスを連立する。
ここで $K$ は局所損失係数（バルブ、エルボ等）、$f$ はDarcy摩擦係数だ。FlownexやAFT Fathomなどの配管解析ツールがこれを実装している。
検証への活用
CFD結果の検証にも使えますか？
もちろん。ベンチュリ管や縮流部のCFD結果をベルヌーイの理論値と比較するのは、基本的な検証手法だ。
| 検証項目 | ベルヌーイ理論値 | CFD結果との比較ポイント |
|---|---|---|
| 縮流部の速度増加 | $v_2 = v_1 \sqrt{A_1/A_2}$（連続の式併用） | 中心流速で比較 |
| 縮流部の圧力降下 | $\Delta p = \frac{1}{2}\rho(v_2^2 - v_1^2)$ | 断面平均圧力で比較 |
| ピトー管の全圧 | $p_0 = p + \frac{1}{2}\rho v^2$ | 淀み点圧力と比較 |
Re数が十分高く（$Re > 10^4$程度）、流れの剥離がなければ、理論値との差は数%以内に収まるはずだ。大きくずれていれば、メッシュか境界条件に問題がある。
つまりベルヌーイはCFDの「サニティチェック」に最適ってことですね。
ベルヌーイで設計できる範囲、できない範囲
ベルヌーイの定理が使える条件は「非粘性・定常・非圧縮・同一流線上」の4つ。実務でよくハマるのが「同一流線上」という条件を忘れることです。例えばポンプ前後の圧力差をベルヌーイで計算する場合、入口と出口の流線が同一かどうかを確認しないと間違いになります。配管設計では「損失ヘッド（損失水頭）」をベルヌーイに加えた拡張形を使うのが標準。摩擦損失・曲管損失・急拡大損失などを損失係数として足し込むことで、粘性効果も近似的に扱えます。本来は「理想流体」の式ですが、うまく拡張して使うのがエンジニアの知恵です。
ベルヌーイの定理の実務適用
実践ガイド
実務でベルヌーイの定理を使う場面って、具体的にはどんなときですか？
代表的な応用例を挙げよう。
ベンチュリ流量計
管路の断面積を絞り、上流と喉部の圧力差から流量を求める。
ここで $C_d$ は流量係数（通常0.95〜0.99）、$A_1$ が上流断面積、$A_2$ が喉部断面積。工業用ベンチュリでは $C_d \approx 0.985$ が典型値だ。
ピトー管
ピトー管はどうですか？
淀み点の全圧と静圧の差から流速を求める。
航空機の対気速度計はまさにこの原理。$p_0 - p = 3000$ Pa、空気密度 $\rho = 1.225$ kg/m$^3$ なら $v \approx 70$ m/s だ。
タンク流出問題（トリチェリの定理）
大きなタンクの底の小さな穴からの流出速度は次の通り。
これはベルヌーイ式で、タンク水面の速度を0（$A_{\text{tank}} \gg A_{\text{hole}}$）とした特別な場合だ。
CFD実務での注意点
CFDで流量計の検証をするとき、気をつけるべきことは？
重要なポイントをまとめよう。
| 項目 | 推奨事項 |
|---|---|
| 上流助走区間 | 管径Dの10倍以上確保 |
| メッシュ | 喉部に10層以上の境界層メッシュ |
| 乱流モデル | SST k-omegaが一般的に良好 |
| 圧力タップ位置 | 壁面圧力を使用、ISO 5167準拠の位置で比較 |
| 収束判定 | 残差 $10^{-5}$ 以下、流量モニタの安定 |
ベルヌーイの式だけで設計するのは危険ですか？
予備設計の段階では十分有用だが、詳細設計では必ずCFDまたは実験で検証すべきだ。特に以下のケースではベルヌーイの仮定が崩れる。
- 急拡大・急縮小: 流れの剥離が発生し損失が大きい
- 高速流: Ma > 0.3 では圧縮性の影響が無視できない
- 旋回流: 遠心力の効果でベルヌーイ式の適用が困難
まずベルヌーイで概算、次にCFDで詳細確認という2段階が実務の王道ですね。
キャビテーションはベルヌーイが引き起こす
ポンプや船のスクリューで問題になる「キャビテーション（空洞現象）」は、ベルヌーイの定理の直接的な帰結です。流速が上がると局所圧力が下がり、水の蒸気圧以下になると気泡が発生する。その気泡が高圧部で潰れるときに発生する衝撃波は、金属材料をえぐり取るほどの破壊力があります。実務でキャビテーション解析をする際は、ベルヌーイで「最低圧力点がどこか」を事前に見積もり、そこが蒸気圧を下回るかどうかを確認するのが設計の第一歩です。船のスクリューデザインがあれほど複雑なのも、このキャビテーション対策が理由のひとつです。
ベルヌーイの定理のソフトウェア比較
ベルヌーイ関連機能の商用ツール対応
ベルヌーイの定理に関連する解析って、各CFDソフトではどう扱われてるんですか？
ベルヌーイそのものは理論式だが、全圧・静圧・動圧の扱い、圧力境界条件の設定、ポスト処理での水頭計算など、各ツールの実装に差がある。
ソルバー別の圧力境界条件
| ソルバー | 圧力入口の指定方法 | 全圧/静圧の扱い | ||
|---|---|---|---|---|
| Ansys Fluent | Pressure Inlet: 全圧を指定 | 内部でベルヌーイ的に静圧を算出 | ||
| Ansys CFX | Total Pressure指定が標準 | Coupled solverで速度と同時に解く | ||
| STAR-CCM+ | Stagnation Inlet: 全圧指定 | 等エントロピー関係も選択可能 | ||
| OpenFOAM | totalPressure BC | $p_0 = p + 0.5\rho | U | ^2$ を内部で適用 |
全圧と静圧を間違えると大変なことになりそうですね。
その通り。典型的なミスは、全圧を指定すべき場所に静圧を入れること。これで流量が全然違ってくる。
1Dシステム解析ツール
ベルヌーイを直接使う1Dツールもあるんですか？
配管ネットワーク解析では以下のツールがベルヌーイ式ベースの解析を行う。
| ツール | 開発元 | 特徴 |
|---|---|---|
| Flownex | Flownex SE | 1Dシステム解析、CFD連成対応 |
| AFT Fathom | Applied Flow Technology | 配管ネットワーク専用 |
| Flowmaster (Simcenter) | Siemens | 自動車・航空の配管系統 |
| RELAP5 | Idaho National Lab | 原子力プラント向け熱水力 |
これらは内部でベルヌーイ式+損失相関を節点ごとに連立して解いている。3D CFDの前段として配管系の圧力分布を把握するのに非常に有効だ。
可視化とポスト処理
CFD結果から全圧分布を見るにはどうすればいいですか？
多くのソルバーで全圧（Total Pressure）は標準出力変数として用意されている。
- Fluent: Report > Surface Integrals で全圧の面平均値を取得
- STAR-CCM+: Field Function で
$TotalPressureを選択 - ParaView (OpenFOAM): Calculator フィルタで
p + 0.5mag(U)^2rhoを計算
全圧のコンター図で損失の大きい箇所が一目瞭然なんですね。
そう。全圧の低下が大きい領域が圧力損失の主因だ。設計改善のターゲットを特定するのに全圧分布は非常に有効なツールだよ。
1D解析ツールとCFDの役割分担
配管システムの設計では、ベルヌーイ+損失係数ベースの1D解析ツール（例：AFT Fathom、FluidFlow）が今も広く使われています。プラント全体の数百ノードの配管ネットワークをCFDで全部解くのは現実的でないため、まず1Dで系全体を評価し、問題箇所だけを3D CFDで詳細解析する——という役割分担が実務標準です。この「1D→3D連成」の設計プロセスを知っていると、ツール選定の議論がずっと具体的になります。ベルヌーイの式は古くて単純に見えますが、システムレベルの最初のふるい分けとして今も第一線で動いています。
ベルヌーイの定理の先端研究
非定常ベルヌーイ方程式
定常流以外にもベルヌーイの式は拡張できるんですか？
もちろん。非定常のオイラー方程式を流線に沿って積分すると、非定常ベルヌーイ式が得られる。
ここで $\phi$ は速度ポテンシャル（$\mathbf{u} = \nabla\phi$）。非定常項 $\partial\phi/\partial t$ が追加される。水撃作用（ウォーターハンマー）の解析や、心臓弁の非定常流にも応用される。
圧縮性ベルヌーイ式
圧縮性流れの場合はどうなりますか？
等エントロピー流れを仮定すると、圧縮性ベルヌーイ式が得られる。
または、全圧と静圧の関係として次の等エントロピー関係式が使われる。
空気（$\gamma = 1.4$）でMa = 0.3の場合、$p_0/p \approx 1.064$ で、非圧縮（$p_0/p \approx 1.065$）との差は0.1%程度。Ma > 0.3で差が顕著になり始める。
回転系への拡張
ターボ機械のような回転する流れではどうですか？
回転座標系でのベルヌーイ式はロタルピー（rothalpy）の保存として表される。
ここで $h$ はエンタルピー、$W$ は相対速度、$U$ は周速。これはポンプやタービン翼の設計で基礎となる関係式だ。
乱流場への適用限界
乱流ではベルヌーイ式は使えないんですか？
時間平均したRANS場では、平均流に対して「修正ベルヌーイ式」を使うことがある。
ここで $k = \frac{1}{2}\overline{u_i'u_i'}$ が乱流エネルギー。変動速度のエネルギーを考慮しないと全圧が過小評価される。これがFluent等で全圧出力に乱流エネルギー項を含めるオプションがある理由だ。
古典的な式でも、現代の数値流体力学と深く結びついてるんですね。
楽器の音もベルヌーイ方程式で動いている
フルートやオルガンパイプの発音原理は「エアリード（空気リード）」と呼ばれ、ベルヌーイ効果が中心的な役割を果たします。吹き口に当たった空気が管の内外でジェット振動を起こし、その周波数が管の共振と合致することで音が出ます。このジェットの不安定性と圧力変動の連成は、非常に複雑な非線形現象です。近年では楽器設計にCFDが使われるようになり、フルート奏者の口形状の違いによる音色変化をシミュレーションする研究も登場しています。ベルヌーイの射程は意外に広いです。
ベルヌーイの定理のトラブル対応
トラブルシューティング
ベルヌーイの定理に関連した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結果と手計算のベルヌーイ値がズレると「計算がおかしい」と思いがちですが、たいていの場合「ベルヌーイの仮定が成立していない」のが原因です。実務でよくある落とし穴は①粘性損失を無視できない管内流れ（損失ヘッドの見落とし）、②剥離や再循環が起きている非定常領域にベルヌーイを適用してしまうこと、③密度変化が大きい高速流れへの誤適用です。ベルヌーイは非常に便利な式ですが「理想流体の特定の流線上」でのみ成立します。まずその前提を一つひとつ確認するのが、結果の整合性チェックの基本です。
関連トピック
なった
詳しく
報告