← 返回知识卡片总目录

DS1号 · 自学教材系列 | 流体力学

第 4 章 理想流体动力学

这一章是流体力学的高潮:把牛顿第二定律用到流体上,先得到「欧拉方程」,再沿流线积分得到「伯努利方程」。伯努利方程不是背出来的公式,它是理想流体里「能量守恒」的化身。我们会完整推一遍,并讲清楚它为什么成立、什么时候失效。

🎯 本章学习目标

  1. 理解理想流体假设(忽略粘性)的物理含义与它的合理性。
  2. 从牛顿第二定律推导出欧拉方程,理解它「力 = 质量 × 加速度」的本质。
  3. 从欧拉方程沿流线积分,完整推导伯努利方程。
  4. 理解伯努利方程每一项的物理意义(压强能、动能、重力势能),及其「能量守恒」本质。
  5. 牢牢记住伯努利方程的适用条件,知道它什么时候不能乱用。

4.1理想流体:把粘性暂时「关掉」

真实流体都有粘性,粘性会让流体在流动中「内耗」、把机械能转成热能,这让数学变得极其复杂(那正是第 6、7 章要啃的硬骨头)。但很多流动中,粘性的影响其实很小,可以先忽略,抓住主要规律。

于是我们引入一个理想化模型:理想流体。

理想流体 假设粘性为零(μ = 0)的流体。它不产生切应力、不内耗、机械能守恒。真实流体没有一个是「理想」的,但这个模型能干净地揭示流动的核心规律,之后再叠加上粘性修正。

理想流体模型为什么合理?看一个例子:水管里高速流动的水,在离管壁较远的主流区,粘性切应力很小,流动几乎就是「无摩擦」的。只有在紧贴壁面的薄薄一层(边界层,第 9 章讲)里粘性才不可忽略。所以「先理想、后修正」是一个聪明的策略。

理想流体 + 不可压缩 + 定常流动,就是推导伯努利方程的三块基石。

先分清两个独立的假设:无粘 ≠ 不可压缩

初学者常把「理想流体」与「不可压缩流体」混为一谈,其实它们是两个互不相干的假设,各自回答不同的问题:

两条假设彼此独立,可以自由组合:水和低速空气两条都近似满足;接近声速的空气「无粘近似还行,但必须考虑压缩性」(第 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 章)、边界层与分离(第 9 章)、高粘油品的低速管输。这套「先理想、后修正」的两步法,是流体力学解决工程问题的标准姿势——就像结构分析先用刚体近似抓主骨架、再用变形体理论补细节一样:抽象模型是逼近真理的脚手架,而不是真理本身。

4.2欧拉方程:流体版的牛顿第二定律

牛顿第二定律说:F = ma(合力 = 质量 × 加速度)。把它用到流体微团上,就得到描述流体运动的方程——欧拉方程。

取一个微小的流体元

在流场中沿流线取一小段「流管微元」:截面积 dA、沿流线长度 ds,如图 4-1 所示。它的质量是 dm = ρ·dA·ds。现在分析它沿流线方向(s 方向)受到的力。

流管微元的受力分析 沿流线取微元,上游压力p,下游压力p+dp,重力分量为ρgdAdssinθ,加速度为vdv/ds。 流线 s 微元 dm p·dA (p+dp)·dA 重力分量 F = ma
图 4-1沿流线取微元:前后压力差与重力分量共同提供合力,等于微元质量乘以加速度。

微元沿流线方向受的力有三部分:

沿流线方向的合力是:

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) 在变,链式法则一展开,两股「变化率」就分了家:

a = dv/dt = ∂v/∂t + (∂v/∂s)·(ds/dt) = ∂v/∂t + v·(dv/ds)

一个直观类比:高速公路的车流是「定常」的——每个路段的限速常年不变,路况不随时间变;但当你从限速 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 章的拉瓦尔喷管,全都建立在这条「速度—压强」兑换关系上。

📖 人物:莱昂哈德·欧拉(Leonhard Euler) 瑞士数学家,史上最高产的数学家之一,一生写下八百多篇论文,晚年双目失明后仍靠心算与口述继续研究。1757 年,他把当时尚且年轻的微积分语言注入流体力学,写出以他命名的运动方程——连续介质力学从此有了第一副完整的数学骨架。有趣的是,欧拉方程比纳维-斯托克斯方程早了近七十年:粘性项是后人补上的,而「无粘的理想世界」这层抽象,正是欧拉的手笔。数学里处处是他的名字——欧拉公式、欧拉常数、欧拉示性数;而在流体力学教材里,他的方程与伯努利的名字永远连在一起。

4.3沿流线积分:推导伯努利方程

现在把欧拉方程沿流线从点 1 积分到点 2,就能得到伯努利方程。这是最关键的一步,跟着推一遍。

欧拉方程各项乘以 ds:

(1/ρ)·dp + v·dv + g·dz = 0

对不可压缩流体,密度 ρ 是常数,可以拿到积分号外。沿流线从点 1 到点 2 积分:

(1/ρ)∫₁² dp + ∫₁² v dv + g∫₁² dz = 0

逐项积分:

代入整理:

(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 年的《流体动力学》中提出了这个方程。有意思的是,他的父亲约翰·伯努利也研究过类似问题,父子之间还有过「谁先发现」的竞争。这个方程后来成为流体力学乃至航空、水利工程最常用的公式之一。

把每一步的「依据」摆上桌面

上面六行推导行云流水,但每一步都有明确依据,值得逐条点名:

换个视角看,整个推导其实是把一个全微分积出来:对不可压缩流体,(1/ρ)dp + v dv + g dz = d(p/ρ + v²/2 + gz) = 0。一个函数沿路径的全微分为零,说明它沿这条路径取常值——伯努利「常数」不是假设,是推导的必然产物。这也解释了方程为什么「长成这样」:三项全都能塞进同一个全微分的括号里。

「常数」对谁不变?——适用范围与两个推广

推导里的积分是沿一条流线做的,所以 p/ρ + v²/2 + gz 只对这条流线守恒;换一条流线,常数一般不同——这正面回答了思考题 4:为什么不能跨流线乱连两点。

但有两个重要的推广,让方程的适用面大幅拓宽:

📐 例题 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² + ρgzPa(单位体积能量)气动计算:动压、总压都是这套单位,皮托管读数直接对接
除以质量 ρp/ρ + v²/2 + gzm²/s²(单位质量能量)热力学分析、可压缩流动(第 11 章的滞止焓就是它的推广)
除以重量 ρgp/ρg + v²/2g + zm(水头,单位重量能量)水力机械与管网:泵扬程、管路损失全用「米」记账

三种形式完全等价,高手会按问题场景挑最顺手的一种:算飞机空速用压强形式,算水管网用水头形式。读工程文献时看到不同的写法,先别慌——看看除的是 ρ 还是 ρg,就知道对方在用哪套账本。

4.4伯努利方程的能量意义与几何意义

伯努利方程不只是个公式,它背后有清晰的物理图像。我们分两层理解。

能量意义:三项 = 三种能量

把方程 p + ρv²/2 + ρgz = 常数 理解为「单位体积流体的能量」:

伯努利方程三项的物理意义(单位体积能量)
项名称物理意义
p压强能(压能)单位体积流体因受压力而「储存」的能量,相当于被压缩的气体/液体的势能
(1/2)ρv²动能(动压)单位体积流体因运动而具有的动能
ρgz重力势能(位能)单位体积流体因处于高度 z 而具有的重力势能

所以伯努利方程本质就是机械能守恒:理想流体(无摩擦、不内耗)沿流线流动时,压强能、动能、重力势能三者可以相互转化,但总和不变。这也反过来解释了为什么「理想流体」才满足它——真实流体有粘性,会不断把机械能耗散成热能,总和就守不住了(第 6、7 章会引入「损失项」修正)。

几何意义:三种「水头」

把方程各项都除以 ρg(即除以单位体积重量),得到长度量纲的形式:

p/(ρg) + v²/(2g) + z = 常数

三项现在都有「高度(米)」的量纲,各叫一个「水头」:

三者之和叫总水头,沿流线守恒。画成图,就是一条水平的「总水头线」:位置水头 + 压强水头 + 速度水头三者此消彼长,总高度不变。这幅图是工程师分析管道的核心工具,第 5 章会反复用到。

总水头线:三种水头此消彼长 沿流线总水头线水平不变,位置水头、压强水头、速度水头三者相互转化。 总水头线(守恒) 流线 v²/2g p/ρg z 速度水头变小 压强水头变大 位置水头 z
图 4-2总水头线守恒:三种水头(速度、压强、位置)此消彼长,总和不变。

🧠 一句话记忆 伯努利方程 = 理想流体的能量守恒:压强能 + 动能 + 重力势能 = 常数。速度大→压强小(在高度不变时),这是它最常用的推论。

两条线要分清:总水头线与测压管水头线

工程师画「水力坡度图」时常用两条线,初学者容易混为一谈:

两条线的间距 = 速度水头:管径越小、流速越大,两线分得越开。图 4-3 用一组竖直测压管把这个间距「画」在了管道上:收缩段流速增大,测压管液面明显下跌——压强水头换成了速度水头,而虚线画出的总水头始终水平。这幅图也是文丘里流量计(第 5 章)的原理草图。

收缩管道上的测压管组 水平变径管道上方装三根竖直测压管,粗段液面高,喉部液面明显下降,扩张段部分回升;总水头线保持水平虚线。 总水头线(水平,守恒) 液面高 喉部液面大跌 部分回升 v 小 v 大 粗段:流速小、压强水头高 喉部:压强水头换速度水头,总水头不变
图 4-3收缩管道上的测压管组:喉部流速增大,压强水头下降、换成速度水头;总水头线(虚线)保持水平——三种水头此消彼长。

📐 例题 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 章讲绕流时会触及)。所以记住一句话:伯努利方程是理想化的、局部的、沿流线的近似,不是放之四海皆准的普适规律。

🧠 使用前必查清单 理想?定常?不可压?同一流线?只有重力?五条都满足才能放心用伯努利方程。

逐条把「为什么」补齐

表格列了五条条件,每一条背后都有机理,值得逐条说透——知道了「为什么」,才能在工程现场灵活判断「能不能用」:

一个能算出数字的反例:突扩管

「有粘性则机械能损耗」这句话,可以算出具体数字来体会。图 4-4 画的是突扩管——小管突然接进大管:流体冲进大管后「刹不住车」,在突扩的角区分离出一片旋涡区,主流先收缩再慢慢扩张。旋涡区里剧烈的掺混把机械能不可逆地搅成了热,于是伯努利方程在这里「多算」了压强回升。例题 5 把「理想账」和「真实账」并排摊开。

突扩管中的分离与耗散 小管突然接入大管,角区形成旋涡分离区,主流先收缩再扩张;理想伯努利预言的压强回升高于实际值,差额被旋涡耗散。 v₁ v₂ 分离区:旋涡把机械能搅成热能(伯努利在此失效) 截面1 截面2 理想预言:p₂−p₁ = ½ρ(v₁²−v₂²) | 实际回升更小,差额 ½ρ(v₁−v₂)² 被旋涡耗散
图 4-4突扩管:角区旋涡耗散机械能,伯努利方程高估了压强回升——「理想流体」条件的定量账单。

📐 例题 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 章的主角),再减一项。写成通用形式:

z₁ + p₁/ρg + v₁²/2g + H泵 = z₂ + p₂/ρg + v₂²/2g + hf + H轮机

原版方程是这台「记账机」的零账本特例:没有泵、没有轮机、没有损失。全书后面几章的管路计算,本质上都是在学习怎么把 hf 算准——沿程损失、局部损失(第 7 章),以及为什么损失是那个形式(第 6 章)。带上这一版方程,就等于把后续两章的路线图提前装进了口袋。

衔接:从「守恒」到「记账」

本章建立的是一个无损失的守恒世界:总水头沿流线一成不变。第 5 章马上把这套守恒律变成测量仪器(皮托管、文丘里管)与计算工具(托里拆利、虹吸);第 6、7 章再把粘性的账算细:先弄清能量为什么漏(层流/湍流的机制),再给出漏多少的工程公式(达西公式、局部损失)。一句话记住全书主线:理想流体给骨架,粘性给修正,能量这本账贯穿始终。

假设的『记账单』:违反哪条,欠什么账 把伯努利方程的适用条件记成一本账,审计思路立刻清晰。欠「无粘」的账:真实流体沿途耗散机械能,管路越长欠得越多——第 6 章的水头损失正是还账的账本;欠「不可压」的账:流速逼近音速三分之一时密度开始变化,压强-速度换算失真——第 11 章的可压缩流动接手(红线已在 4.5 节划过);欠「定常」的账:方程里的常数随时间漂移,启动、阀门突关这类瞬态工况必须把时间项请回来;欠「同一流线」的账:两条流线各有各的常数,跨线直接套用等于拿别家的存折取钱。记账的目的不是吓退使用,而是遇到可疑结果时,第一时间知道去哪本账上找误差来源。

本章小结

  1. 理想流体是忽略粘性(μ=0)的模型,用于先抓住流动的核心规律。
  2. 欧拉方程 (1/ρ)dp/ds + v·dv/ds + g·dz/ds = 0 是把牛顿第二定律用到流体微团的结果,是理想流体的运动方程。
  3. 沿流线积分欧拉方程(不可压缩)得到伯努利方程:p + (1/2)ρv² + ρgz = 常数。
  4. 伯努利方程的三项分别是压强能、动能、重力势能,本质是机械能守恒;几何上对应压强水头、速度水头、位置水头,总水头守恒。
  5. 伯努利方程只在理想、定常、不可压缩、沿同一流线、只受重力时成立;真实粘性流体会耗散机械能,需要加损失项修正。

读完本章,你应该能……

思考题

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):此时可以证明全流场共用同一个常数,直管均匀段近似满足;有旋时(如强弯曲流道)常数逐线不同。工程计算沿主流线取点,执行的就是这条纪律。