ベルヌーイの定理
理論と物理
概要
先生、ベルヌーイの定理ってよく聞くんですけど、何がそんなに大事なんですか？
定常・非粘性・非圧縮の流れにおけるエネルギー保存則だ。流体力学で最も基本的かつ強力なツールの一つだよ。流速と圧力の関係を直接結びつけるから、配管設計やピトー管、ベンチュリ管など至るところで使われる。
支配方程式の導出
どうやって導出するんですか？
オイラー方程式（非粘性流体の運動方程式）の流線に沿った成分を考える。定常流で流線に沿って積分すると、ベルヌーイの式が得られる。
オイラー方程式は次の通りだ。
定常かつ流線に沿って積分すると、ベルヌーイの式になる。
各項はそれぞれ何を表してるんですか？
第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）の理論で説明されます。ベルヌーイの定理は「速いほど圧力が低い」という道具として正しいですが、なぜ速くなるかの説明が教科書でサボられているのです。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
ベルヌーイ式の数値的活用
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つ。実務でよくハマるのが「同一流線上」という条件を忘れることです。例えばポンプ前後の圧力差をベルヌーイで計算する場合、入口と出口の流線が同一かどうかを確認しないと間違いになります。配管設計では「損失ヘッド（損失水頭）」をベルヌーイに加えた拡張形を使うのが標準。摩擦損失・曲管損失・急拡大損失などを損失係数として足し込むことで、粘性効果も近似的に扱えます。本来は「理想流体」の式ですが、うまく拡張して使うのがエンジニアの知恵です。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
実務でベルヌーイの定理を使う場面って、具体的にはどんなときですか？
代表的な応用例を挙げよう。
ベンチュリ流量計
管路の断面積を絞り、上流と喉部の圧力差から流量を求める。
ここで $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の解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
ベルヌーイ関連機能の商用ツール対応
ベルヌーイの定理に関連する解析って、各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連成」の設計プロセスを知っていると、ツール選定の議論がずっと具体的になります。ベルヌーイの式は古くて単純に見えますが、システムレベルの最初のふるい分けとして今も第一線で動いています。
選定で最も重要な3つの問い
- 「何を解くか」：ベルヌーイの定理に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
非定常ベルヌーイ方程式
定常流以外にもベルヌーイの式は拡張できるんですか？
もちろん。非定常のオイラー方程式を流線に沿って積分すると、非定常ベルヌーイ式が得られる。
ここで $\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結果と手計算のベルヌーイ値がズレると「計算がおかしい」と思いがちですが、たいていの場合「ベルヌーイの仮定が成立していない」のが原因です。実務でよくある落とし穴は①粘性損失を無視できない管内流れ（損失ヘッドの見落とし）、②剥離や再循環が起きている非定常領域にベルヌーイを適用してしまうこと、③密度変化が大きい高速流れへの誤適用です。ベルヌーイは非常に便利な式ですが「理想流体の特定の流線上」でのみ成立します。まずその前提を一つひとつ確認するのが、結果の整合性チェックの基本です。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ベルヌーイの定理の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告