伯努利定理
伯努利定理的理论基础
概述
老师，伯努利定理是我经常听到的，为什么这么重要呢？
定常、非粘性、非压缩流中的能量守恒定律。流体力学中最基础、最强大的工具之一。直接关联流速与压力，在管道设计、皮托管、文丘里管等处广泛应用。
支配方程导出
怎样推导呢？
欧拉方程（非粘性流体运动方程）沿流线方向的分量。定常流沿流线积分，得到伯努利式。
定常且沿流线积分，得到伯努利式。
各项分别表示什么？
第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基本上解纳维尔-斯托克斯方程，不直接离散求解伯努利式。但边界条件设置、结果检证、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的"健全性检查"最适合的工具啊。
伯努利能设计的范围和不能的范围
伯努利定理的适用条件是"非粘性、定常、非压缩、同一流线上"四个。实务中常见的失误是忽略"同一流线上"。比如计算泵前后的压力差用伯努利时，入口和出口的流线是否一致没有确认就算了，结果就会错。配管设计的标准做法是对伯努利式加入"损失水头（损失高度）"扩展形。摩擦损失、弯管损失、急扩大损失等作为损失系数加进去，这样就能近似处理粘性效应。本来是"理想流体"的式子，但通过巧妙的扩展而成为实用工具，这就是工程师的智慧。
伯努利定理的实务应用
实践指南
实务中什么时候会用到伯努利定理？具体有什么例子？
列举代表性应用。
文丘里流量计
管道的截面积缩小，用上游与喉部的压力差求流量。
其中 $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精确检证，这是实务的标准做法吧。
空蚀是伯努利定理引起的
泵和船舶螺旋桨出现的"空蚀（cavitation）"现象正是伯努利定理的直接后果。流速上升导致局部压力下降，低于水的蒸气压时产生气泡。这些气泡在高压部塌陷时产生的冲击波能摧毁金属材料。设计时用伯努利先找"最低压点在哪儿"，检查它是否降至蒸气压以下——这是第一步。船舶螺旋桨设计那么复杂，空蚀防护是重要原因之一。
伯努利定理的软件对比
伯努利相关功能的商用工具对应
伯努利定理相关的解析，各CFD软件怎么处理的？
伯努利本身是理论式，但全压、静压、动压的处理，压力边界条件设置，后处理中水头计算，各工具有所不同。
求解器的压力边界条件
| 求解器 | 压力入口的指定方法 | 全压/静压的处理 |
|---|---|---|
| Ansys Fluent | Pressure Inlet：指定全压 | 内部用伯努利计算静压 |
| Ansys CFX | Total Pressure指定为标准 | 耦合求解器同时解速度和压力 |
| 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.5*mag(U)^2*rho
全压等值线图能一眼看出损失大的地方啊。
对。全压降大的区域就是压力损失的主因。设计改进的目标定位用全压分布图非常高效。
1D解析工具与CFD的角色分工
配管系统设计中，伯努利+损失系数基础的1D解析工具（如AFT Fathom、FluidFlow）至今广泛使用。整个工厂数百个节点的配管网络全部CFD求解在现实中不可行，标准做法是先用1D评估整体系统，再对问题部分进行3D CFD详细分析——这就是"1D→3D连成"的设计流程。这个"1D→3D"的分层思路掌握好了，工具选择讨论会变得具体有力。伯努利式看似古老单纯，但作为系统级评估的第一道门槛，今天仍在第一线发挥作用。
伯努利定理的前沿研究
非定常伯努利方程
定常流以外能扩展伯努利式吗？
可以。对非定常欧拉方程沿流线积分，得到非定常伯努利式。
其中 $\phi$ 是速度势（$\mathbf{u} = \nabla\phi$）。加入了非定常项 $\partial\phi/\partial t$。水锤效应（管道中的压力波）解析和心脏瓣膜的非定常流也能应用。
压缩性伯努利式
压缩流怎么处理？
等熵流假设下，得到压缩性伯努利式。
或者，全压与静压的等熵关系为：
空气（$\gamma = 1.4$）Ma = 0.3时，$p_0/p \approx 1.064$，与非压缩时差约0.1%。Ma > 0.3开始明显差异。
旋转系的扩展
旋转的流像泵那样怎么办？
旋转坐标系下的伯努利式表现为转焓（rothalpy）守恒。
其中 $h$ 是焓，$W$ 是相对速度，$U$ 是周速。这是泵和涡轮设计的基础关系式。
湍流场的适用限界
湍流中伯努利式不适用吗？
时均RANS场中对平均流可用"修正伯努利式"。
其中 $k = \frac{1}{2}\overline{u_i'u_i'}$ 是湍流能。不考虑速度脉动的能量，全压会被低估。所以Fluent等的全压输出有选项加湍流能项。
古老的式子也与现代数值流体力学深度联系啊。
乐器的音也由伯努利方程驱动
长笛和风琴管这样的乐器发音原理是"气流簧"（air reed），伯努利效应起中心作用。吹口的气流在管内外产生喷流振荡，当频率与管道共振相匹配时就发出声音。这种喷流的不稳定性与压力变化耦合是非常复杂的非线性现象。最近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+不在壁函数适用范围（应在30～300）
- 湍流模型不当（误用层流计算）
对策：设置壁面粗糙度参数（钢管约 $k_s \approx 0.045$ mm）。调整y+到30～300。根据Re数选用合适的湍流模型。
结局就是，有伯努利理论值的话，必须和CFD对比。这是找问题最近的路。
完全同意。理论可以预测的值一定要与CFD结果对标。这是调试的最短路径。
"伯努利不符"往往是假设破裂
CFD结果与手算的伯努利值差异大时，通常会怀疑"计算有问题"，但多数情况其实是"伯努利的假设失效"。实务中常见的陷阱：①粘性损失不能忽视的管内流（损失水头遗漏）、②剥离或再循环的非定常区域硬用伯努利、③密度变化大的高速流误用。伯努利是强大的工具，但"理想流体的特定流线上"才成立。一个一个确认前提，是结果整合的基础。
关联主题
价值
更详细
错误