4.1 欧拉方程与伯努利方程
4.1 欧拉方程与伯努利方程
本节摘要：忽略粘性后 N-S 方程退化为欧拉方程；沿流线积分（定常、不可压缩、无体积力做功路径）得到伯努利方程 p + ½ρU² + ρgz = 常数，即压能、动能、位置势能之和守恒。本节用皮托管、文丘里流量计、孔口出流三个仪表级应用演练这条方程，并明确列出它的失效场景。
对着两张纸的中间吹气，纸会向中间靠拢而不是被吹开。这是伯努利效应最便宜的家庭实验：流速高的地方静压低。1700 年代丹尼尔·伯努利把这件事写成了能量账本，如今每架飞机的空速管、每座水厂的流量计、每支化学实验室的抽滤泵，都还在用这本账。
一、从欧拉方程到能量守恒
把 N-S 方程的粘性项 μ∇²u 划掉，剩下欧拉方程。在"定常、不可压缩、沿一条流线"三个条件下积分，得到：
p + ½ρU² + ρgz = 常数（沿流线）。
三项分别是静压（压能）、动压（动能）与位置水头（势能）。物理图像就一句话：无粘、无损失的流动里，这三种能量可以互相兑换，总额不变。水管由粗变细，流速上升，静压必须下降——不是"被抽走了"，是能量换了形态。可压缩或非定常时各需修正（前者换等熵关系，后者补当地加速度项），但核心思想不变。
动压是"驻点账户"里的钱
把手指挡在水管口，水在你指腹前停驻，动能全额转成压力：驻点压力 = 静压 + 动压。这正是皮托管的原理——测总压与静压之差，反推速度：
import numpy as np # 皮托管测风速：压差 500 Pa 对应多少风速 dp = 500.0 # 总压减静压 Pa rho = 1.225 # 空气密度 U = np.sqrt(2 * dp / rho) print(f"风速 {U:.1f} m/s") # 反过来 飞机 100 m/s 飞行时动压 U_true = 100.0 print(f"100 m/s 飞行的动压 {0.5*rho*U_true**2/1e3:.1f} kPa")
500 Pa 压差对应约 28.6 m/s 的风速；飞机 100 m/s 飞行时动压约 6.1 kPa——足够让空速膜盒产生可观的机械变形，机械式空速表不需要任何电子元件。压差传感器加一根小管，就是从伯努利方程到航空仪表的完整产业链。
二、文丘里流量计与孔口出流
文丘里管是收缩-扩张的缩放管：喉道处流速升高、静压下降，测收缩段与喉道的压差反推流量。联合连续性方程消去速度比，得到流量公式 Q = C·A₂·√(2Δp/(ρ(1−β⁴)))，β 为喉道与入口直径比，C 为考虑粘性损失的流量系数（通常 0.95–0.99）。
孔口出流：水箱侧壁开孔，水面到孔心高差 h，出流速度 v = √(2gh)——托里拆利公式，伯努利方程的免费赠品。注意真实孔流要乘流量系数（锐孔约 0.61），因为射流在孔外先收缩（缩脉）。
# 文丘里管：水路 D=100mm d=50mm 压差 8 kPa import numpy as np D, d = 0.10, 0.05 beta = d / D dp, rho = 8e3, 998.0 Cd = 0.98 A2 = np.pi * d**2 / 4 Q = Cd * A2 * np.sqrt(2*dp / (rho * (1 - beta**4))) print(f"流量 {Q*1e3:.2f} L/s 喉道流速 {Q/A2:.2f} m/s") # 水箱孔口 水面高差 3 m h, C_or = 3.0, 0.61 print(f"孔口理论出流速度 {np.sqrt(2*9.81*h):.1f} m/s 实际流量系数 {C_or}")
出流速度约 7.7 m/s——与"从 3 米高处自由落体"完全一致，因为动能正是从同样高的势能兑换来的。伯努利方程在这里把落体公式平移到了流体上。
三、伯努利方程的四个"不许用"
⚠️ 常见坑清单（每条都有真实事故背景）：
- 跨流线乱用：方程常数只在同一流线上（或无旋的全场）。把上游慢流线的常数套到下游快流线上，结果是能量凭空多出来；
- 有泵/风机段：设备对流体做功，必须加功率项，否则算出的扬程永远不够；
- 分离区：流线脱离壁面、能量损失严重的地方（突扩、阀门后）不满足"无损失"，要按损失系数或突扩公式另算；
- 可压缩高速流：密度变化大时改用等熵关系（第 5 章），ρ 当常数的伯努利会低估动压。
经典工程事故类型：按理想伯努利选泵，忽略管路沿程损失与局部损失，结果实际流量只有设计值的一半。伯努利是"理想情况的下限账本"，工程使用必须叠上损失项（第 3 章达西公式）。
四、虹吸的完整账单
虹吸管是伯努利方程最好的综合练习：从高处水池翻过坝顶流向低处。取水池自由面与管顶最高点两截面列伯努利（同一条流线、两端通大气），最高点的压强为
p_top = p_atm − ρg·h（h 为最高点高出水面的距离）。
rho, g = 998.0, 9.81 for h in [2.0, 5.0, 9.0]: p_top = 101325 - rho*g*h print(f"顶点高出水面 {h:.0f} m 顶点绝对压 {p_top/1e3:6.1f} kPa")
顶点高出 9 m 时绝对压已降到约 13 kPa，溶解气析出会让管内积气断流——虹吸管越高越容易"断"，工程虹吸灌溉与虹吸雨水排放系统的顶高都控制在约 8 m 以内，并在顶点设排气阀。这个例子把"压力能与势能互换"演到了极限工况：伯努利方程不止算流量，还算得出系统里压力最低、最先出事的点——找最低压力点是泵、水轮机、喷油器设计里反复出现的例行步骤（汽蚀判定）。
变式思考
伯努利方程还有一条常被忽略的推论：沿流线"总压"守恒意味着能量一旦被损失吃掉就回不来。风道里一片滤网、一只阀门的损失，都要靠风机额外补回——这也是暖通设计中"低阻力滤网省的是风机电费"的定量出处。把每个部件的损失系数折算成电费，节能改造的优先级排序就一目了然。
本节要点回顾
- 欧拉方程是 N-S 方程去掉粘性项的版本，沿流线积分得伯努利方程：压能、动能、势能之和守恒；
- 驻点压力 = 静压 + 动压，皮托管测总静压差反推速度，是空速测量的根基；
- 文丘里计、孔口出流分别把能量互换用于计量与泄流，各有必须乘的流量系数；
- 四大禁用场景：跨流线、含功设备、分离区、高速可压缩——越界即错；
- 理想无粘结果做一阶近似，粘性损失当修正项叠加，是工程流体力学的标准工作流。
一条流线上的账算清了，下一节把整张流场当拼图来拼。