📖 人物:莱昂哈德·欧拉(Leonhard Euler) 瑞士数学家,史上最高产的数学家之一,一生写下八百多篇论文,晚年双目失明后仍靠心算与口述继续研究。1757 年,他把当时尚且年轻的微积分语言注入流体力学,写出以他命名的运动方程——连续介质力学从此有了第一副完整的数学骨架。有趣的是,欧拉方程比纳维-斯托克斯方程早了近七十年:粘性项是后人补上的,而「无粘的理想世界」这层抽象,正是欧拉的手笔。数学里处处是他的名字——欧拉公式、欧拉常数、欧拉示性数;而在流体力学教材里,他的方程与伯努利的名字永远连在一起。
🎯 本章学习目标
- 理解理想流体假设(忽略粘性)的物理含义与它的合理性。
- 从牛顿第二定律推导出欧拉方程,理解它「力 = 质量 × 加速度」的本质。
- 从欧拉方程沿流线积分,完整推导伯努利方程。
- 理解伯努利方程每一项的物理意义(压强能、动能、重力势能),及其「能量守恒」本质。
- 牢牢记住伯努利方程的适用条件,知道它什么时候不能乱用。
4.1理想流体:把粘性暂时「关掉」
真实流体都有粘性,粘性会让流体在流动中「内耗」、把机械能转成热能,这让数学变得极其复杂(那正是第 6、7 章要啃的硬骨头)。但很多流动中,粘性的影响其实很小,可以先忽略,抓住主要规律。
于是我们引入一个理想化模型:理想流体。
理想流体 假设粘性为零(μ = 0)的流体。它不产生切应力、不内耗、机械能守恒。真实流体没有一个是「理想」的,但这个模型能干净地揭示流动的核心规律,之后再叠加上粘性修正。
理想流体模型为什么合理?看一个例子:水管里高速流动的水,在离管壁较远的主流区,粘性切应力很小,流动几乎就是「无摩擦」的。只有在紧贴壁面的薄薄一层(边界层,第 9 章讲)里粘性才不可忽略。所以「先理想、后修正」是一个聪明的策略。
理想流体 + 不可压缩 + 定常流动,就是推导伯努利方程的三块基石。
先分清两个独立的假设:无粘 ≠ 不可压缩
初学者常把「理想流体」与「不可压缩流体」混为一谈,其实它们是两个互不相干的假设,各自回答不同的问题:
- 无粘(μ = 0)回答的是力的问题——流体内部有没有内摩擦?切应力是否为零?它决定机械能会不会在流体内部被耗散成热。
- 不可压缩(ρ = 常数)回答的是物性的问题——流体受力后密度变不变?体积缩不缩?它决定压强功能不能写成 (1/ρ)∫dp 这样简洁的形式。
两条假设彼此独立,可以自由组合:水和低速空气两条都近似满足;接近声速的空气「无粘近似还行,但必须考虑压缩性」(第 11 章);蜂蜜几乎不可压缩、粘性却极强;糖浆被快速搅拌时又粘又可压。推导伯努利方程时两条都用——所以 4.5 节的适用条件表里,「理想流体」与「不可压缩」是并列的两行,而不是同一件事的两种说法。
还有一个容易被忽略的「副作用」值得点破:真实流体在固体壁面必须满足无滑移条件——贴壁那层流体与壁面同速,速度从壁面的零一层层被「拖」起来。而理想流体完全没有切应力,允许流体沿壁面滑移。换句话说,理想流体模型天生描述不了紧贴壁面的那一层:它把全场最难啃的边界层(第 9 章)整个「抹掉」了。这不是模型的缺陷,而是分工——理想流体负责主流,边界层理论负责壁面,两者拼起来才是完整图景。判断某处该用哪个模型,钥匙是第 6 章的雷诺数。
📐 例题 1 主流区切应力的量级估算:为什么敢忽略粘性?
水管内水流平均速度 v = 1 m/s,水的动力粘度 μ = 1.0×10−3 Pa·s,密度 ρ = 1000 kg/m³。(a) 主流区速度剖面很平缓,取速度梯度 du/dy ≈ 1 s⁻¹,估算粘性切应力;(b) 壁面附近速度在约 0.5 mm 的薄层内从 0 涨到主流速度,估算该处的速度梯度与切应力;(c) 把两处切应力与动压 ½ρv² 比较,说明粘性在哪些区域可以忽略。
思路 切应力用牛顿内摩擦定律 τ = μ(du/dy) 逐点估算;壁面附近用 du/dy ≈ v/δ 估算(δ 为速度从 0 到主流速度的薄层厚度);最后与动压比较、给结论。
解 (a) 主流区:τ = 1.0×10⁻³ × 1 = 1.0×10⁻³ Pa——只有千分之一帕斯卡,微不足道。(b) 壁面附近:du/dy ≈ 1 / 0.0005 = 2000 s⁻¹,τ ≈ 1.0×10⁻³ × 2000 = 2 Pa,是主流区的一千倍。(c) 动压 ½ρv² = 0.5×1000×1² = 500 Pa。主流区 τ/动压 ≈ 2×10⁻⁶——粘性完全可忽略;壁面区 τ/动压 ≈ 4×10⁻³——绝对值仍然不大,但比主流区大了三个数量级,而且越贴近壁面梯度越陡,粘性正是在这里「东山再起」。
答:主流区 τ ≈ 10⁻³ Pa 可忽略;壁面区 τ ≈ 2 Pa,粘性影响集中在壁面薄层。 ⚠ 易错点:「忽略粘性」不是全场开关,而是分区域的近似——主流里忽略、壁面附近不忽略,这正是边界层思想的雏形。
「理想流体」的适用地图
哪些场景能放心用理想流体模型?工程经验给出一张粗略的地图:
- 绕流的主流区:机翼表面以外的绝大部分流场、建筑物周围的风场——粘性只通过薄薄的边界层起作用,主流几乎是「无摩擦」的,用理想流体算出的表面压强分布与实测相差无几。
- 远离壁面的管道流动:管网初步水力计算先按理想流体估压差,再叠加沿程损失修正(第 7 章的「水头损失」就是给理想模型打的补丁)。
- 水波与明渠流动的远场:自由表面流动的能量估算,粘性几乎插不上手。
反过来,粘性必须「在场」的问题也一目了然:管壁附近的速度剖面与流动阻力(第 7 章)、边界层与分离(第 9 章)、高粘油品的低速管输。这套「先理想、后修正」的两步法,是流体力学解决工程问题的标准姿势——就像结构分析先用刚体近似抓主骨架、再用变形体理论补细节一样:抽象模型是逼近真理的脚手架,而不是真理本身。
4.2欧拉方程:流体版的牛顿第二定律
牛顿第二定律说:F = ma(合力 = 质量 × 加速度)。把它用到流体微团上,就得到描述流体运动的方程——欧拉方程。
取一个微小的流体元
在流场中沿流线取一小段「流管微元」:截面积 dA、沿流线长度 ds,如图 4-1 所示。它的质量是 dm = ρ·dA·ds。现在分析它沿流线方向(s 方向)受到的力。
微元沿流线方向受的力有三部分:
- 上游压力:p·dA,沿 s 正方向。
- 下游压力:(p + dp)·dA,沿 s 负方向。净压力 = p·dA − (p+dp)·dA = −dp·dA。
- 重力分量:重力 dm·g 竖直向下,它沿流线方向的分量是 −dm·g·sinθ = −ρ·dA·ds·g·(dz/ds),其中 dz/ds = sinθ 是流线沿高度 z 的变化率。
沿流线方向的合力是:
Fs = −dp·dA − ρ·dA·ds·g·(dz/ds)
加速度是多少?
微团沿流线的加速度。定常流动中,微团速度随位置变化(对流加速度):微团从 s 运动到 s+ds,速度从 v 变到 v+dv,加速度 a = dv/dt = (dv/ds)·(ds/dt) = v·(dv/ds)。
由牛顿第二定律 Fs = dm·a:
−dp·dA − ρ·dA·ds·g·(dz/ds) = ρ·dA·ds·v·(dv/ds)
两边同除以 ρ·dA·ds,移项整理得:
欧拉方程(沿流线的一维形式) (1/ρ)·(dp/ds) + v·(dv/ds) + g·(dz/ds) = 0
这就是流体版的 F=ma。它把「压力变化、速度变化、高度变化」三者沿流线联系了起来。注意:欧拉方程因为忽略了粘性,所以是「理想流体」的运动方程;如果把粘性项加回来,就得到更复杂的纳维-斯托克斯(N-S)方程(第 6 章再提)。
定常流动里哪来的加速度?——质点导数
正文那句「微团速度随位置变化(对流加速度)」是初学者最容易卡住的地方:流场明明定常——每一点的速度读数永远不变——微团凭什么还有加速度?答案藏在质点导数的视角里。加速度的定义是「跟着同一个微团」看速度的时间变化率 dv/dt;定常流场中速度是位置的函数 v = v(s),而微团在动、位置 s = s(t) 在变,链式法则一展开,两股「变化率」就分了家:
- 第一项 ∂v/∂t叫局部加速度(时变项):流场本身随时间变化贡献的加速度,定常流动中它恒为零。
- 第二项 v·(dv/ds)叫对流加速度(迁移项):微团「换地方」贡献的加速度——它从流速 2 m/s 的地段挪进流速 6 m/s 的地段,自然越走越快,哪怕任何一点的速度读数从头到尾纹丝不动。
一个直观类比:高速公路的车流是「定常」的——每个路段的限速常年不变,路况不随时间变;但当你从限速 60 的匝道汇入限速 120 的主路时,你踩了油门,你就有加速度。变的不是路况,而是你所在的位置。流体微团穿过收缩段被「提速」,正是同一回事:不是哪一点在变快,而是它在「换更快的点」。
非定常时两项都在。比如管道末端的阀门突然关闭,全管流速整体随时间衰减,∂v/∂t 席卷全场,管内压强剧烈波动——那是水锤分析的舞台(第 7 章)。所以欧拉方程的一般形式是 (1/ρ)(∂p/∂s) + ∂v/∂t + v(dv/ds) + g(dz/ds) = 0,正文推导的定常形式只是它的特例;4.5 节「定常」这条适用条件的根源,就在这里。
📐 例题 2 收缩喷管的对流加速度与压强降落
水平放置的收缩喷管(类似灭火水枪、燃油喷嘴的前端),水流沿流线从点 1 流到点 2,两点流速分别为 2 m/s 与 6 m/s,沿流线相距 0.5 m。水的密度 ρ = 1000 kg/m³。(a) 估算这段路程上微团的加速度;(b) 用欧拉方程算该段上的压强降落;(c) 用伯努利方程复核,并说明两法为何必然一致。
思路 (a) 用平均速度梯度 dv/ds ≈ Δv/Δs,配合中点速度代入 a = v(dv/ds);(b) 欧拉方程给出 dp/ds = −ρ(v·dv/ds + g·dz/ds),水平流动 dz/ds = 0;(c) 伯努利直接算 Δp = ½ρ(v₁² − v₂²)。
解 (a) 平均速度梯度 dv/ds ≈ (6−2)/0.5 = 8 s⁻¹;取中点流速 v ≈ (2+6)/2 = 4 m/s,得 a ≈ 4 × 8 = 32 m/s²——约 3.3 倍重力加速度,微团在半米之内被猛烈提速。(b) dp/ds = −ρ·a = −1000 × 32 = −3.2×10⁴ Pa/m;该段压降 Δp = −3.2×10⁴ × 0.5 = −1.6×10⁴ Pa = −16 kPa。(c) 伯努利(水平、同一流线):Δp = ½ρ(v₁²−v₂²) = 500×(4−36) = −16 kPa ✓ 与 (b) 完全一致——伯努利方程本来就是欧拉方程沿流线的积分,两者是同一条定律的「微分版」与「积分版」。
答:a ≈ 32 m/s²;喷管内压强沿程降落 16 kPa。 ⚠ 易错点:①加速度不是 (6−2)/0.5 = 8「m/s²」——8 s⁻¹ 是速度梯度,还要再乘流速本身才是加速度;②压强沿收缩方向降低,与「流速大压强低」方向一致,但因果是「压差推动加速」,不是速度把压强「吸」小了。
压力与重力之外:完整的受力盘点
正文的受力清单只写了压力与重力,严格来说还欠两句交代。其一,微元侧表面上也有压强作用,但它们沿侧向互相抵消,对沿流线方向的合力没有贡献——所以沿 s 方向的动量方程里确实只剩前后压差与重力分量,这不是偷懒,是矢量投影的严谨。其二,重力是体积力、压强是表面力,两者在方程里以「单位质量的力」口径统一记账,最后才化简得只剩三项。理解了「谁在方程里、谁不在、为什么」,公式就不再是黑箱。
再看三维的欧拉方程:ρ[∂v/∂t + (v·∇)v] = −∇p + ρg。左边是「密度×加速度」((v·∇)v 就是三维化的对流加速度),右边是压力梯度与体积力。正文的流线一维形式,是矢量方程沿流线方向 s 投影、再除以 ρ 的结果。看懂这层「投影」关系,才能明白 4.3 节末尾的重要推广——无旋流动里伯努利常数全场同值——为什么成立:那时积分路径可以任选,不再被流线「锁死」。
工程里的欧拉方程:压强是「布置」出来的
欧拉方程在工程上最直接的用途,是压强分布的推演与设计。涡喷发动机的进气道在低速飞行时是收缩通道:流速沿程增大,按欧拉方程压强必须沿程降低——空气被「抽」进来,为发动机提供进气动量;反过来,进气道的扩压段让流速降低、压强回升,把动能「压」回压力能,为燃烧室增压。同一份数学,一收一扩、方向相反:加速靠压强换,增压靠流速换。第 5 章的皮托管、文丘里管,第 11 章的拉瓦尔喷管,全都建立在这条「速度—压强」兑换关系上。
4.3沿流线积分:推导伯努利方程
现在把欧拉方程沿流线从点 1 积分到点 2,就能得到伯努利方程。这是最关键的一步,跟着推一遍。
欧拉方程各项乘以 ds:
(1/ρ)·dp + v·dv + g·dz = 0
对不可压缩流体,密度 ρ 是常数,可以拿到积分号外。沿流线从点 1 到点 2 积分:
(1/ρ)∫₁² dp + ∫₁² v dv + g∫₁² dz = 0
逐项积分:
- ∫ dp = p₂ − p₁
- ∫ v dv = v₂²/2 − v₁²/2
- ∫ dz = z₂ − z₁
代入整理:
(p₂ − p₁)/ρ + (v₂² − v₁²)/2 + g(z₂ − z₁) = 0
把点 1 的项放到一边、点 2 的项放到另一边,得到:
伯努利方程(不可压缩、理想、定常,沿流线) p + (1/2)ρv² + ρgz = 常数,或写成两点形式:p₁ + (1/2)ρv₁² + ρgz₁ = p₂ + (1/2)ρv₂² + ρgz₂。
这就是著名的伯努利方程。它说:沿同一条流线,这个组合量「p + ρv²/2 + ρgz」保持不变。
📖 人物:丹尼尔·伯努利(Daniel Bernoulli) 瑞士伯努利家族是科学史上罕见的「科学世家」,出了好几代杰出数学家物理学家。丹尼尔·伯努利在 1738 年的《流体动力学》中提出了这个方程。有意思的是,他的父亲约翰·伯努利也研究过类似问题,父子之间还有过「谁先发现」的竞争。这个方程后来成为流体力学乃至航空、水利工程最常用的公式之一。
把每一步的「依据」摆上桌面
上面六行推导行云流水,但每一步都有明确依据,值得逐条点名:
- 「乘以 ds」:欧拉方程沿流线处处成立(每一点都满足),所以能沿曲线积分;ds 是沿流线的微元弧长,与方程同方向,不会「积歪」。
- 「ρ 提出积分号」:不可压缩假设正式登场的地方——ρ 为常数才提得出来。若可压缩,(1/ρ)∫dp 必须先知道 p 与 ρ 的关系(等温?等熵?)才能积分,那是第 11 章可压缩流动的入口。
- 「∫v dv = v²/2」:凑微分 v dv = d(v²/2)。
- 「∫g dz = gz」:重力是保守力,g dz = d(gz),所以它有资格被积成一个「势」。
换个视角看,整个推导其实是把一个全微分积出来:对不可压缩流体,(1/ρ)dp + v dv + g dz = d(p/ρ + v²/2 + gz) = 0。一个函数沿路径的全微分为零,说明它沿这条路径取常值——伯努利「常数」不是假设,是推导的必然产物。这也解释了方程为什么「长成这样」:三项全都能塞进同一个全微分的括号里。
「常数」对谁不变?——适用范围与两个推广
推导里的积分是沿一条流线做的,所以 p/ρ + v²/2 + gz 只对这条流线守恒;换一条流线,常数一般不同——这正面回答了思考题 4:为什么不能跨流线乱连两点。
但有两个重要的推广,让方程的适用面大幅拓宽:
- 无旋推广:若流动无旋(全场旋度为零——均匀来流绕过物体之前的远场、水箱刚放水时的初始阶段都近似如此),可以证明每条流线的常数都相同。伯努利方程从「沿一条流线」升级为「全流场通用」。工程上大量绕流问题近似满足无旋条件,这正是这个方程被用到「泛滥」的原因。
- 缓变流截面近似:管流工程里,在缓变流截面(流线近似平行、弯曲很小的横截面)上,各点的 p/ρg + z 近似相等,于是可以把整个截面「当成一个点」,与下游另一个截面连用伯努利方程——第 7 章管路计算的日常操作。注意两条截面的连线穿过的是「同一个缓变流场的两个截面」,本质仍是沿流动方向的能量核算。
📐 例题 3 变径又爬升的水管:两点伯努利全流程
输水管从点 1(管径 100 mm,流速 2 m/s,表压 200 kPa,标高 2 m)流到点 2(管径 50 mm,标高 5 m),流动定常,忽略粘性损失,水的密度 ρ = 1000 kg/m³,g = 9.8 m/s²。求点 2 的流速与表压。
思路 两步走:①连续性方程定 v₂(同一管流、质量守恒);②沿流线写两点伯努利解出 p₂。表压是相对大气的读数,两点同用表压,大气压项两边同时约掉。
解 ① v₂ = v₁A₁/A₂ = v₁(d₁/d₂)² = 2×(100/50)² = 8 m/s。② 伯努利:p₁ + ½ρv₁² + ρgz₁ = p₂ + ½ρv₂² + ρgz₂,解出 p₂ = p₁ + ½ρ(v₁²−v₂²) + ρg(z₁−z₂) = 200000 + 500×(4−64) + 9800×(2−5) = 200000 − 30000 − 29400 = 140600 Pa ≈ 140.6 kPa(表压)。
答:v₂ = 8 m/s,p₂ ≈ 140.6 kPa。动压增大吃掉 30 kPa,爬升 3 m 又吃掉 29.4 kPa——压强能是唯一的「提款机」。 ⚠ 易错点:①面积比是直径比的平方,直接代 100/50 会差 4 倍;②爬升项是 z₁−z₂ = −3 m,符号代错答案会偏出 58.8 kPa;③表压与绝对压两套体系不能混用——两点必须同用表压(或同用绝对压)才能把大气压约掉。
辨析:「流速大所以压强小」错在哪
口语常把伯努利方程总结成「流速大、压强小」,这句话只在同一流线、高度不变、理想定常的前提下成立,而且极易把因果说反。严格的逻辑链是:截面积缩小 → 连续性要求流速增大 → 动能增加需要能量 → 能量只能从压强能里来 → 压强必然降低。是压差在推动加速,而不是「速度大 somehow 把压强吸小了」。理顺因果,公式的方向感就有了:已知压差求速度增益,或已知速度变化反推压差,都不会代错。
顺着这条逻辑还能预先看出方程的「死穴」:如果流动里有旋涡、分离、阀门、弯头这类把机械能搅成热的地方,能量账就不平了——第 4.5 节的突扩例题会给出具体的「亏空数字」。此外,若两点间有泵或涡轮(外界加入/抽走能量),也得在方程两端加减「设备功」——水泵扬程、水轮机水头,说的都是这件事。
同一个方程的三种「记账单位」
细心的读者可能发现:4.3 节写成 p + ½ρv² + ρgz,4.4 节又写成 p/ρg + v²/2g + z——它们是同一个方程的三种「记账单位」,除的参考量不同而已,用途各有侧重:
| 除以什么 | 形式 | 单位 | 典型使用者 |
|---|---|---|---|
| 除以体积(不除) | p + ½ρv² + ρgz | Pa(单位体积能量) | 气动计算:动压、总压都是这套单位,皮托管读数直接对接 |
| 除以质量 ρ | p/ρ + v²/2 + gz | m²/s²(单位质量能量) | 热力学分析、可压缩流动(第 11 章的滞止焓就是它的推广) |
| 除以重量 ρg | p/ρg + v²/2g + z | m(水头,单位重量能量) | 水力机械与管网:泵扬程、管路损失全用「米」记账 |
三种形式完全等价,高手会按问题场景挑最顺手的一种:算飞机空速用压强形式,算水管网用水头形式。读工程文献时看到不同的写法,先别慌——看看除的是 ρ 还是 ρg,就知道对方在用哪套账本。
4.4伯努利方程的能量意义与几何意义
伯努利方程不只是个公式,它背后有清晰的物理图像。我们分两层理解。
能量意义:三项 = 三种能量
把方程 p + ρv²/2 + ρgz = 常数 理解为「单位体积流体的能量」:
| 项 | 名称 | 物理意义 |
|---|---|---|
| p | 压强能(压能) | 单位体积流体因受压力而「储存」的能量,相当于被压缩的气体/液体的势能 |
| (1/2)ρv² | 动能(动压) | 单位体积流体因运动而具有的动能 |
| ρgz | 重力势能(位能) | 单位体积流体因处于高度 z 而具有的重力势能 |
所以伯努利方程本质就是机械能守恒:理想流体(无摩擦、不内耗)沿流线流动时,压强能、动能、重力势能三者可以相互转化,但总和不变。这也反过来解释了为什么「理想流体」才满足它——真实流体有粘性,会不断把机械能耗散成热能,总和就守不住了(第 6、7 章会引入「损失项」修正)。
几何意义:三种「水头」
把方程各项都除以 ρg(即除以单位体积重量),得到长度量纲的形式:
p/(ρg) + v²/(2g) + z = 常数
三项现在都有「高度(米)」的量纲,各叫一个「水头」:
- 压强水头 p/(ρg):把流体压到某个压强所需的水柱高度。
- 速度水头 v²/(2g):以速度 v 竖直上抛能到达的高度(v² = 2gh 反推)。
- 位置水头 z:流体所处的高度。
三者之和叫总水头,沿流线守恒。画成图,就是一条水平的「总水头线」:位置水头 + 压强水头 + 速度水头三者此消彼长,总高度不变。这幅图是工程师分析管道的核心工具,第 5 章会反复用到。
🧠 一句话记忆 伯努利方程 = 理想流体的能量守恒:压强能 + 动能 + 重力势能 = 常数。速度大→压强小(在高度不变时),这是它最常用的推论。
两条线要分清:总水头线与测压管水头线
工程师画「水力坡度图」时常用两条线,初学者容易混为一谈:
- 总水头线:H = p/ρg + v²/(2g) + z 沿程的连线。理想流体中它是水平线;真实流体中它沿程下降(第 7 章的「水头损失」就体现在这条线的坡度上)。
- 测压管水头线(水力坡度线):Hp = p/ρg + z 的连线,处处比总水头线低一个速度水头 v²/(2g)。明渠流里它就是水面线;有压管道里,它就是「在管壁接一根竖直测压管、液面能升到的高度」。
两条线的间距 = 速度水头:管径越小、流速越大,两线分得越开。图 4-3 用一组竖直测压管把这个间距「画」在了管道上:收缩段流速增大,测压管液面明显下跌——压强水头换成了速度水头,而虚线画出的总水头始终水平。这幅图也是文丘里流量计(第 5 章)的原理草图。
📐 例题 4 水头核算:总水头真的不变
接例题 3:点 1(表压 200 kPa,v₁ = 2 m/s,z₁ = 2 m)与点 2(表压 140.6 kPa,v₂ = 8 m/s,z₂ = 5 m),水 ρ = 1000 kg/m³,g = 9.8 m/s²。(a) 分别计算两点的压强水头、速度水头、位置水头与总水头,验证守恒;(b) 计算两点测压管水头的落差,说明它去哪了。
思路 水头是「单位重量流体的能量」:逐点求 p/ρg、v²/2g、z 三项之和;测压管水头 = p/ρg + z。
解 (a) 点 1:p₁/ρg = 200000/9800 = 20.41 m,v₁²/2g = 4/19.6 = 0.20 m,z₁ = 2 m,总水头 H₁ = 22.61 m。点 2:p₂/ρg = 140600/9800 = 14.35 m,v₂²/2g = 64/19.6 = 3.27 m,z₂ = 5 m,总水头 H₂ = 22.61 m。两点完全相等 ✓。(b) 测压管水头:Hp1 = 20.41+2 = 22.41 m;Hp2 = 14.35+5 = 19.35 m;落差 = 3.06 m——恰好等于速度水头的增量(3.27−0.20 = 3.06 m)✓。
答:两点总水头同为 22.61 m(守恒验证通过);测压管水头线从 22.41 m 跌到 19.35 m,跌掉的 3.06 m 全部换成了速度水头。 ⚠ 易错点:总水头线与测压管水头线是两条不同的线——把测压管水头误当「总能量」,会在变径管处得出「能量不守恒」的错误结论。
工程视角:为什么工程师爱说「水头」不爱说「压强」
翻开水泵样本或管网设计手册,工程师满口都是「扬程 32 m」「水头损失 4 m」,原因有二。其一,水头是长度量纲,直观可比——「损失 4 m」就是每吨水白白丢掉 4 m 高度的能量,人人有体感;其二,不同设备、不同密度介质(清水、液压油、航空燃油)的能量账用「米」记就能横向比较:泵的扬程减去管路的总损失,剩下的就是还能爬多高、喷多远,加减法直接做。这套「水头语言」的语法基础正是本节的几何解释——能量折算成高度,账目一目了然。第 5 章的水泵—管路计算、第 7 章的损失计算,都将用这门语言书写。
收尾预演:驻点——速度归零,压强封顶 沿流线走到物体正前方,流体被逼停的那一点 v = 0,伯努利方程只剩压强一项——压强升到这条流线的天花板,这就是驻点,该处压强叫驻点压强(总压)。拿数字找手感:风以 25 m/s 正对告示牌吹来,空气密度取 1.225 kg/m³,驻点表压 = ½×1.225×25² ≈ 383 Pa,相当于牌面每平方米多承受约 383 牛——台风天广告牌最先从中心开裂,物理根源就在这个「压强峰」上。第 5 章的皮托管干脆把驻点当探头:量出总压与静压之差,流速就能反算出来——驻点是本章理论通向下一章仪器的桥。
4.5适用条件:什么时候能用、什么时候不行
伯努利方程好用,但滥用它是初学者最常犯的错。它有一串严格的成立条件,必须逐条核对:
| 条件 | 含义 | 不满足时的后果 |
|---|---|---|
| 理想流体 | 无粘性,μ = 0 | 有粘性则机械能损耗,需加损失项 |
| 定常流动 | 流场不随时间变化 | 非定常时需加非定常项 |
| 不可压缩 | ρ 为常数 | 可压缩时需用气体动力学处理 |
| 沿同一条流线 | 两点必须在同一流线上 | 跨流线一般不再守恒 |
| 只受重力 | 无其他体积力 | 有其他力(如电磁力)需另加项 |
为什么「理想流体」才满足?
这是最值得追问的一条。回顾 4.3 的推导:我们沿流线积分时,只写了压力差和重力两个力,完全没写粘性切应力——因为理想流体假设 μ=0,切应力 τ = μ(du/dy) = 0。如果流体有粘性,沿流线还会有一个「摩擦力」项,它做的功会把机械能变成热,于是 p + ρv²/2 + ρgz 就不再守恒,而会沿流向不断减小。所以:
核心理解 伯努利方程「机械能守恒」的成立,根源于没有粘性耗散。理想流体是它成立的物理前提;真实流体的粘性越强、流得越慢、越靠近壁面,伯努利方程的误差就越大。
一个常见的反例
两块平行板之间的流动(如第 1 章图 1-3 那样的 Couette 流)里,粘性占主导,伯努利方程根本不能用。再比如,机翼升力虽然常用伯努利方程解释,但那也是「近似」——真实的升力还要靠粘性导致的环量(第 10 章讲绕流时会触及)。所以记住一句话:伯努利方程是理想化的、局部的、沿流线的近似,不是放之四海皆准的普适规律。
🧠 使用前必查清单 理想?定常?不可压?同一流线?只有重力?五条都满足才能放心用伯努利方程。
逐条把「为什么」补齐
表格列了五条条件,每一条背后都有机理,值得逐条说透——知道了「为什么」,才能在工程现场灵活判断「能不能用」:
- 定常条件:推导中用了 a = v(dv/ds),悄悄丢掉了 ∂v/∂t 项。阀门突开突关、活塞泵往复推动的管流、海洋潮汐通道都是非定常的,必须补上非定常项(沿流线对 ∂v/∂t 积分)——水锤(第 7 章)正是「含时间伯努利」的经典应用。
- 不可压缩条件:ρ 提出积分号的代价。气体流速超过约 100 m/s(马赫数超过 0.3)后密度变化超过百分之几,(1/ρ)∫dp 积不出来,必须换用等熵关系(第 11 章)。
- 沿流线条件:常数「一线一值」。两个补救:无旋流动全场同值(4.3 节);管流取缓变流截面近似连线(第 7 章的日常做法)。
- 只受重力条件:推导里体积力只写了 g。若还有电磁力(磁流体、等离子体)或旋转机械里的离心力场,方程要加项。
- 理想流体条件:核心中的核心,下面用突扩管的例子给这条「记账」。
一个能算出数字的反例:突扩管
「有粘性则机械能损耗」这句话,可以算出具体数字来体会。图 4-4 画的是突扩管——小管突然接进大管:流体冲进大管后「刹不住车」,在突扩的角区分离出一片旋涡区,主流先收缩再慢慢扩张。旋涡区里剧烈的掺混把机械能不可逆地搅成了热,于是伯努利方程在这里「多算」了压强回升。例题 5 把「理想账」和「真实账」并排摊开。
📐 例题 5 突扩管的「理想账」与「真实账」
水管从 d₁ = 50 mm 突然扩大到 d₂ = 100 mm,小管流速 v₁ = 4 m/s(水,ρ = 1000 kg/m³,g = 9.8 m/s²)。(a) 求大管流速 v₂;(b) 若突扩前后无任何损失(理想伯努利),压强能回升多少?(c) 实验与动量方程(第 12 章)给出真实回升 p₂ − p₁ = ρv₂(v₁ − v₂),求真实值与「漏掉」的能量;(d) 把漏掉的能量折算成水头损失。
思路 (a) 连续性方程;(b) 水平流动、同一流线,用伯努利两点式;(c) 按题给的真实关系式计算,差额即损失;(d) 损失压强除以 ρg。
解 (a) v₂ = v₁(d₁/d₂)² = 4×(50/100)² = 1 m/s。(b) 理想回升 = ½ρ(v₁²−v₂²) = 500×(16−1) = 7500 Pa = 7.5 kPa。(c) 真实回升 = ρv₂(v₁−v₂) = 1000×1×3 = 3000 Pa = 3 kPa,只有理想值的 40%;「漏掉」7500−3000 = 4500 Pa,恰等于 ½ρ(v₁−v₂)² = 500×9 = 4500 Pa——这就是著名的博尔达–卡诺(Borda–Carnot)损失。(d) 水头损失 = 4500/9800 ≈ 0.46 m:每吨水流过这个突扩,白白丢掉近半米高程的机械能。
答:v₂ = 1 m/s;理想回升 7.5 kPa,真实仅 3 kPa,损失 4.5 kPa(合 0.46 m 水头)。 ⚠ 易错点:①突扩前后两点不能直接连伯努利——中间隔着分离旋涡区;②真实回升式来自动量方程,不是经验猜测;③同样的扩径若改成 5°–7° 锥角的渐扩管,分离消失、损失降到一成以下——文丘里管的扩散段就是这么设计的(第 5 章)。
伯努利方程的「亲戚」:加项与减项
工程实践里很少用到「原汁原味」的伯努利方程,更常用的是它的「带账本版」:两点之间若有水泵(外界对流体做功,单位重量流体获得扬程 H泵),方程左边加一项;有水轮机/涡轮(流体对外做功),减一项;有粘性损失 hf(第 7 章的主角),再减一项。写成通用形式:
原版方程是这台「记账机」的零账本特例:没有泵、没有轮机、没有损失。全书后面几章的管路计算,本质上都是在学习怎么把 hf 算准——沿程损失、局部损失(第 7 章),以及为什么损失是那个形式(第 6 章)。带上这一版方程,就等于把后续两章的路线图提前装进了口袋。
衔接:从「守恒」到「记账」
本章建立的是一个无损失的守恒世界:总水头沿流线一成不变。第 5 章马上把这套守恒律变成测量仪器(皮托管、文丘里管)与计算工具(托里拆利、虹吸);第 6、7 章再把粘性的账算细:先弄清能量为什么漏(层流/湍流的机制),再给出漏多少的工程公式(达西公式、局部损失)。一句话记住全书主线:理想流体给骨架,粘性给修正,能量这本账贯穿始终。
假设的『记账单』:违反哪条,欠什么账 把伯努利方程的适用条件记成一本账,审计思路立刻清晰。欠「无粘」的账:真实流体沿途耗散机械能,管路越长欠得越多——第 6 章的水头损失正是还账的账本;欠「不可压」的账:流速逼近音速三分之一时密度开始变化,压强-速度换算失真——第 11 章的可压缩流动接手(红线已在 4.5 节划过);欠「定常」的账:方程里的常数随时间漂移,启动、阀门突关这类瞬态工况必须把时间项请回来;欠「同一流线」的账:两条流线各有各的常数,跨线直接套用等于拿别家的存折取钱。记账的目的不是吓退使用,而是遇到可疑结果时,第一时间知道去哪本账上找误差来源。
本章小结
- 理想流体是忽略粘性(μ=0)的模型,用于先抓住流动的核心规律。
- 欧拉方程 (1/ρ)dp/ds + v·dv/ds + g·dz/ds = 0 是把牛顿第二定律用到流体微团的结果,是理想流体的运动方程。
- 沿流线积分欧拉方程(不可压缩)得到伯努利方程:p + (1/2)ρv² + ρgz = 常数。
- 伯努利方程的三项分别是压强能、动能、重力势能,本质是机械能守恒;几何上对应压强水头、速度水头、位置水头,总水头守恒。
- 伯努利方程只在理想、定常、不可压缩、沿同一流线、只受重力时成立;真实粘性流体会耗散机械能,需要加损失项修正。
读完本章,你应该能……
- 说出理想流体的定义,以及为什么它能「先抓住主要规律」。
- 独立推导欧拉方程,并说明它如何体现 F=ma。
- 独立从欧拉方程积分出伯努利方程,一步不跳。
- 解释伯努利方程三项各自的能量意义和水头意义。
- 完整列出五条适用条件,并解释为什么「理想流体」是机械能守恒的前提。
- 区分总水头线与测压管水头线,并能画图说明两者间距的含义。
- 用突扩管反例定量说明「滥用伯努利」会高估多少压强回升。
- 在压强、单位质量、水头三种归一形式之间自由换算,并知道各自适合的工程场景。
思考题
1. 为什么伯努利方程只对「理想流体」严格成立?粘性会破坏什么?
提示:推导中沿流线只写了压力差和重力,没写粘性切应力。有粘性时切应力做功把机械能变热,p+ρv²/2+ρgz 沿流向不断减小,不再守恒。
2. 用能量观点解释:为什么水管变细(高度不变)时,流速变大、压强变小?
提示:连续性方程要求细处流速大;伯努利方程(z 不变)要求 p + ρv²/2 守恒,v 增大则 p 必然减小。动能增、压能减,能量相互转化。
3. 伯努利方程里,p/(ρg)、v²/(2g)、z 三项都有长度量纲,分别叫什么「水头」?物理上各代表什么?
提示:压强水头(把流体压到该压强所需水柱高)、速度水头(以此速度上抛能到的最大高度)、位置水头(所在高度)。
4. 两个点分别位于两条不同的流线上,能直接用伯努利方程连它们吗?为什么?
提示:不能。伯努利方程是「沿同一条流线」积分的结果,常数只对同一条流线成立;不同流线上的常数一般不同,不能跨流线乱用。
5. 高速公路每个路段的限速常年不变(「路况定常」),你从限速 60 的匝道汇入限速 120 的主路时仍有加速度。类比流体,解释「定常流动中微团为何仍有加速度」。
提示:定常说的是「场不随时间变」,加速度说的是「微团的速度随时间变」。微团在换位置、穿过不均匀的速度场,于是有了迁移加速度 v(dv/ds);只有场既定常又均匀时,它才为零。
6. 突扩管的小管截面 1 与大管截面 2,能直接用两点伯努利方程吗?应该怎么处理?
提示:不能。突扩产生分离旋涡,机械能耗散为热,p₁+½ρv₁² 不再等于 p₂+½ρv₂²。处理办法:在下游侧减去损失项 ½ρ(v₁−v₂)²,或改用动量方程求真实压差(第 12 章)。
7. 沿同一条流线的伯努利「常数」,数值由什么决定?两条流线的常数一定不同吗?
提示:常数由该流线的积分起点(边界条件,比如上游水箱液面)决定;一般流动中两条流线的常数可以不同,但无旋流动中全流场同值——这是伯努利方程能全场通用的例外情形,也是它被广泛使用的深层原因。
8. 伯努利方程有三种写法:p + ½ρv² + ρgz、p/ρ + v²/2 + gz、p/ρg + v²/2g + z。它们是什么关系?各自适合什么场景?
提示:同一个方程除以不同参考量——体积、质量、重量——得到三种归一形式,数值不同、信息等价。压强形式(Pa)对接气动与皮托管,水头形式(m)对接泵扬程与管路损失,单位质量形式在可压缩流动(第 11 章)里最自然。
9. 风以 25 m/s 吹向一块大告示牌正中,空气密度 1.225 kg/m³。求驻点相对来流的表压,并换算成每平方米多少牛。
提示:驻点处 v = 0,动压全部转成压强:p_stag − p∞ = ½ρv² = 0.5×1.225×25² ≈ 383 Pa,即每平方米约 383 N。牌面其余位置的压强低于驻点,这个「中心压强峰」正是迎风面受力的主力来源。
10. 同一横截面上,靠管壁与管心两条不同流线上的点,能直接共用同一个伯努利常数吗?什么情况下可以全流场共用?
提示:一般不能——常数沿流线成立,两条流线各有各的值。特例是无旋流动(第 3 章的涡量 = 0):此时可以证明全流场共用同一个常数,直管均匀段近似满足;有旋时(如强弯曲流道)常数逐线不同。工程计算沿主流线取点,执行的就是这条纪律。