← 返回知识卡片总目录

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

第 6 章 粘性流体与流动阻力

前几章我们「关掉」了粘性,得到了漂亮的伯努利方程。这一章把粘性重新打开:真实流体流动时,为什么会分「层流」和「湍流」两种截然不同的状态?雷诺数这个神奇的无量纲数为什么能一眼判断流动状态?它的物理意义是什么?

🎯 本章学习目标

  1. 理解粘性流体流动中存在的两种流态:层流与湍流。
  2. 掌握雷诺数的定义,并深刻理解它的物理意义——惯性力与粘性力之比。
  3. 知道临界雷诺数,会判断圆管流动是层流还是湍流。
  4. 理解层流与湍流的内部结构差异,以及湍流为何阻力大。
  5. 把粘性切应力(牛顿内摩擦)重新放回流动分析中,理解「阻力」从哪来。

6.1雷诺实验:两种流态的发现

1883 年,英国物理学家奥斯本·雷诺做了一个后来载入史册的实验:让水流过一根透明玻璃管,并在管入口处滴入染色液体,观察染色线的形态。他发现了两种截然不同的流动状态。

雷诺实验:层流与湍流 上:层流,染色线保持直线清晰;中:过渡;下:湍流,染色线散开搅浑。 层流(流速小) 染色线保持直线、清晰 → 分层流动 过渡(流速中等) 染色线开始摆动、波浪状 湍流(流速大)
图 6-1雷诺实验:流速小时染色线平直(层流),流速增大后染色线摆动、最终散开搅浑(湍流)。

📖 人物:奥斯本·雷诺(Osborne Reynolds) 雷诺是 19 世纪英国的工程师兼物理学家,是「把流体力学带进工程」的关键人物之一。他不仅发现了两种流态,还提出了那个以他命名的判据——雷诺数,一举解决了「什么条件下是层流、什么条件下是湍流」这个困扰工程师已久的问题。

实验设计的三个门道:为什么雷诺能「看见」流态

这个实验看似简单,实则处处是功夫,门道至少有三条:

先纠一个直觉偏差:层流不等于「慢」

日常语言里「水流很急」几乎等于湍流,这不算全错,但工程上必须精确:流态不取决于速度的绝对大小,只取决于雷诺数(下一节)。同样的 1.2 m/s:在 50 mm 管里输水,Re 高达 6 万,必然湍流(例题 2);换输送粘稠的润滑油,Re 只有约 260,是从容的层流。反过来,大江大河虽然「缓慢」,但尺度巨大,Re 达百万量级,照样湍流。「细管快流」多半湍流,「粗管慢流」且流体粘稠时可以是层流——流态是 ρ、v、L、μ 四个量的「团体赛」,不是速度一个人的独角戏。

📐 例题 1 家用软管接水:你的水管是什么流态?

家用软管内径 d = 16 mm,水温 20 ℃(ν = 1.0×10⁻⁶ m²/s)。(a) 水龙头全开,管内平均流速约 0.8 m/s,判断流态;(b) 把流速调小到 0.10 m/s 呢?(c) 若想让水在管内保持层流,流速最多多大?此时接满 1 L 水要多久?

思路 Re = vd/ν,与 2300 比较;(c) 由 Re = 2300 反解流速,再用 Q = vA 与 t = V/Q 换时间。

解 (a) Re = 0.8×0.016/1.0×10⁻⁶ = 12800 ≫ 4000——湍流。(b) Re = 0.10×0.016/10⁻⁶ = 1600 < 2300——层流。(c) vmax = 2300×10⁻⁶/0.016 ≈ 0.14 m/s;流量 Q = vA = 0.14×(π/4×0.016²) ≈ 2.9×10⁻⁵ m³/s = 0.029 L/s,接满 1 L 约 35 秒——比全开龙头(湍流状态下不到 10 秒)慢得多。

答:(a) Re≈12800,湍流;(b) Re≈1600,层流;(c) 最大约 0.14 m/s,接 1 L 水约 35 s。 ⚠ 易错点:①家用软管的标称尺寸常是外径,壁厚吃掉 3–4 mm,别拿标称值当内径硬算;②Re 落在过渡区(2300–4000)时,工程上按湍流处理才保险。

转捩不是「开灯」,是「黎明」

染色线从「笔直」到「搅浑」之间,雷诺其实观察到了一段连续的中间过程:先是染色线出现规则的正弦式摆动——微弱的行波在管中成长;随后摆动加剧、局部崩出零星的「湍流斑」——一小片一小片的湍流像麻疹一样在层流里冒头、随流向下游扩张;最后湍流斑连成一片,流动完全湍流。这个从层流到湍流的演变叫转捩。它提醒我们:临界雷诺数 2300 不是一堵墙,更像一条「黄昏线」——真实管流里,扰动的强弱(入口毛刺、弯头、泵的脉动)会把这条线左右挪动几百甚至上千。工程上按 2300 保守判断,正是因为转捩位置本身是「看扰动脸色」的。

💡 深挖:染色线为什么层流不散、湍流秒散?——两种「扩散」差六个量级 染料分子自己也会散开,靠的是分子扩散(布朗运动),水中的扩散系数 D ≈ 10⁻⁹ m²/s——一秒钟只推开 √(2Dt) ≈ 0.05 mm,一分钟也才 0.4 mm 上下,这就是层流里染色线能保持「笔直」一整个管长的原因:不是不散,是散得太慢看不见。湍流里换成涡扩散:涡团横冲直撞地搬运染料,等效「扩散系数」高达 10⁻³–10⁻² m²/s,比分子扩散快六个数量级——一秒钟搅匀几厘米,整管染料瞬间花掉。同一根染色线的两种命运,背后是两种扩散机制的 10⁶ 倍差距;风洞里烟线之所以能清晰勾画流谱,也是同一个道理:只要流动是层流,烟线就是「冻结」在流线上的。

工程上谁在乎流态?

在乎流态的工程师比想象中多。输油管线:原油粘度大、常为层流,沿程损失远小于同条件的湍流,故「加热降粘」「加减阻剂」都是围绕流态与阻力做文章;液压系统:设计在层流区工作——噪声低、能耗稳、控制精确;换热器:恰恰相反,刻意制造湍流,用涡团掺混把冷热流体「搅」在一起,换热效率成倍提高;生物医学工程:血流大部分区段是层流(主动脉 Re≈2000,贴着临界值),瓣膜后方的湍流区恰是血栓与心脏杂音的温床。同一个雷诺数,四种行业、四种悲喜。

💡 拓展:现代流场可视化——雷诺思路的技术延续 雷诺用一根染色线「看见」了流态,一百多年后的工程师仍在用同样的思路,只是工具换了代:烟线法——风洞里从细金属丝上释放白烟,把绕机翼的流动拍成可见的流谱;氢气泡法——水槽中用通电细丝电解出氢气泡,给水流打上「时间刻度线」;油膜法——在模型表面涂薄油层,气流剪切让油膜流成条纹,直接显示壁面流动方向(转捩线、分离线的位置一目了然);PIV 粒子图像测速——撒入示踪小粒子,用激光片光照亮一个切面、连续拍两帧,由粒子位移反算出整个切面的速度场——等于把雷诺的「一条染色线」升级成「千万条染色线」+「自动读数」。风洞试验报告里那一张张漂亮的流谱图,祖师爷都是 1883 年那根染色线。

6.2雷诺数:惯性力与粘性力的较量

雷诺实验之后,下一个问题自然浮现:到底由什么决定流动是层流还是湍流?是流速?是管径?是流体的粘度?雷诺发现,答案不是一个量,而是一个组合起来的无量纲数。

雷诺数的定义

雷诺数 Re 定义为:

Re = ρvL / μ = vL / ν

其中 ρ 是密度、v 是特征流速、L 是特征长度(圆管流里取管直径 d)、μ 是动力粘度、ν 是运动粘度。注意 Re 是个无量纲数——分子分母单位约掉了,所以它的大小与用什么单位制无关。

物理意义:惯性力 / 粘性力

这是理解雷诺数最关键的一点,务必弄懂它的「出身」。

两者之比:

惯性力 / 粘性力 ≈ (ρv²L²) / (μvL) = ρvL/μ = Re

量纲推导:这两个「力」是怎么估出来的

上面的量级不是拍脑袋,各推一步就明白。取边长为 L 的流体微团,速度 v,特征时间 t ≈ L/v(流过自身尺度的时间):

有了这层出身,就能理解为什么 Re 小时流动「稳」:粘性相对强大,任何扰动刚冒头就被抹平,像果冻里戳一下立刻复原;Re 大时惯性独大,扰动被流动放大成涡团、滚雪球成湍流——像推倒一排多米诺。

💡 细节:特征长度 L 怎么选? Re 里的 L 取「最能代表该流动的几何尺度」:圆管取内径 d(不是半径——临界值 2300 就是对 d 定义的,换成半径数字就全变了);沿平板流动取离前缘的距离 x(边界层沿程增厚,Rex 沿程增大);非圆截面用当量直径 de = 4A/χ(A 为过流面积,χ 为湿周)。原则一句话:取「速度梯度发生方向上的代表性尺度」。用错特征长度,Re 的数值就失去可比性,临界值也跟着作废。

核心理解 雷诺数 = 惯性力 / 粘性力。Re 小,粘性力占上风,扰动被「按住」,流动保持层流;Re 大,惯性力占上风,扰动被「放大」,流动失稳变成湍流。所以雷诺数刻画的是两种力量的此消彼长。

这就是雷诺数的魅力:它用一个数字,把「流速、尺度、粘度」对流动状态的影响打包成一个判据。不管是一根细管里慢流的水,还是一条大河里奔腾的水,只要 Re 相同,流动的「性格」就相似(这正是第 8 章相似原理的思想源头)。

临界雷诺数与判断

对圆管流动,实验给出:

工程上常取 Rec ≈ 2300 作为判断圆管层流/湍流的临界值。注意这个值是下临界雷诺数(层流一定稳定的上限);实际从层流转到湍流的「上临界值」可以高得多,取决于扰动大小,但工程上用 2300 保守判断。

「下临界」与「上临界」值得多说两句。下临界雷诺数 ≈ 2300:Re 低于它,无论怎么扰动,层流都稳如泰山——这是「层流一定安全的下界」。上临界雷诺数:把装置做到极致安静(消除振动、光滑入口、过滤杂质)时,层流甚至能维持到 Re 上万;但轻轻一碰立刻崩塌成湍流——这是「层流苟延残喘的上界」,数值因装置而异,没有工程价值。所以 2300 是「保底判据」:低于它放心按层流算;高于它按湍流处理,最坏只是保守。工程判断永远站在保守一侧,这个思想后面会反复出现。

雷诺数的量级地图(各领域典型值)
流动场景典型参数Re 量级流态与解读
毛细血管d≈8 μm,v≈0.5 mm/s(血液 ν≈3.5×10⁻⁶)≈10⁻³粘性绝对碾压,蠕动般层流
主动脉d≈25 mm,v≈0.3 m/s≈2000贴着临界值,剧烈运动时可转湍(对应心脏杂音)
家用水管d≈25 mm,v≈1 m/s10⁴–10⁵日常用水几乎总是湍流
大江大河水深≈5 m,v≈2 m/s10⁶–10⁷天然水域必然湍流,掺混带来强自净能力
客机机翼绕流弦长≈5 m,v≈250 m/s(空气)10⁷–10⁸极高 Re:边界层薄,转捩与分离主导气动设计

读表的姿势:Re 跨越 11 个数量级,流体世界的「性格光谱」全在这一个数上——毛细血管里粘性是暴君,机翼上方惯性是国王。同一门力学,尺度决定脾气。

💡 深挖:无量纲数 = 流动「性格」的通用编码 为什么工程师如此钟爱无量纲数?因为它与单位制无关——用 SI 还是英制算出同一个 Re,性格判断完全一致;更因为它是「力与力之比」,天然适合做模型实验:想用 1:10 的缩比模型预测真实飞行器的流动,只要把模型试验调到相同 Re,两种尺度下的惯性/粘性较量比例相同,流态就相似(第 8 章相似原理的核心)。但这里有个昂贵的技术难题:模型缩小 10 倍,要保持 Re 不变,就得把流速×10 或把运动粘度 ν 降到 1/10——把风「吹快十倍」会撞上空气压缩性(马赫数限制),把 ν 降下来则催生了增压风洞、低温风洞这类造价惊人的设备:往风洞里灌低温高压的氮气,把 ν 压到常温空气的几分之一,只为让缩比模型「活得像真机」。一个无量纲数,牵着一代代风洞的预算。

💡 拓展:Re ≪ 1 的世界——粘性是绝对的国王 量级地图的左端(Re ≪ 1)是一个我们直觉陌生的世界,叫爬流(蠕动流):惯性小到可以整体忽略,粘性成为唯一的力学语言。这里的流动极端「温顺」:扰动了立刻平息,没有分离涡,没有湍流,永远层流。球在爬流中受的阻力有精确解——斯托克斯律 F = 6πμrv(r 为球半径、v 为球速),阻力与速度严格线性。这是微米世界的物理:微珠在液体中的沉降、粉尘颗粒在空气中的飘落、微流控芯片里的流动、机器轴与轴承间油膜的挤压流动,全部活在 Re ≪ 1 的疆域。用人类的直觉去想象它:像在糖浆里游泳,惯性毫无用处,你停下来流动就瞬间停住——连划水都「存不下」任何动量。

📐 例题 2 同样的管、同样的速度:水与油两种命运

同一条 d = 50 mm 的管道,平均流速 v = 1.2 m/s。(a) 输送 20 ℃ 的水(ν = 1.0×10⁻⁶ m²/s),判断流态;(b) 改输送润滑油(μ = 0.2 Pa·s,ρ = 870 kg/m³)呢?(c) 若仍输水,流速小到多少才能保持层流?

思路 Re = vd/ν;(b) 先换算 ν = μ/ρ 再算;(c) 令 Re = 2300 反解 v。

解 (a) Re = 1.2×0.05/1.0×10⁻⁶ = 60000 ≫ 4000——湍流。(b) ν = 0.2/870 ≈ 2.3×10⁻⁴ m²/s;Re = 1.2×0.05/2.3×10⁻⁴ ≈ 261 ≪ 2300——层流。(c) vmax = 2300×10⁻⁶/0.05 ≈ 0.046 m/s——水要保持层流,只能以每秒不到 5 厘米的速度「爬」,几乎不实用;而油以 1.2 m/s「快跑」仍是从容的层流。

答:水 Re=60000 湍流;油 Re≈261 层流;水保持层流需 v ≤ 0.046 m/s。 ⚠ 易错点:①μ 与 ν 差一个密度,代公式前先核对单位;②「流速快→湍流」只在同一种流体里成立,跨流体比较必须看 Re。

💡 细节:「充分发展」与入口段——公式都写在成熟区 本章的抛物线剖面、泊肃叶公式,前提都是「充分发展」:剖面形状沿程不再变化。刚进管子的流体可不是这样——入口处剖面近似均匀,边壁的粘性作用像「冻层」一样逐渐向管心蔓延,走一段距离后剖面才长成稳定的抛物线。这段「成长期」叫入口段,层流时长约 0.05·Re·d(管径的几十到上百倍),湍流时短得多,约 10–60 倍管径。工程含义有二:其一,测压孔、流量计要装在充分发展段之后,否则读到的是「未成年」的剖面,标定作废;其二,例题 1 的软管(湍流,入口段约 40 倍管径 ≈ 0.6 m)恰好整根都接近充分发展——家用尺度下这几乎是自动满足的,大装置里则必须认真留出这个距离。

🧠 一句话记忆 Re = ρvL/μ = 惯性力/粘性力。Re 小→粘性镇压→层流;Re 大→惯性造反→湍流。圆管 Re≈2300 是分界线。

辨析:判流态要用『最不利工况』的粘度 雷诺数里的 ν 不是常数,它随温度大幅变动——液体粘度对温度极其敏感,20 ℃ 粘稠的机械油到 50 ℃ 能稀一大截。同一条 d = 50 mm 管路、同一流速 1.2 m/s:20 ℃ 时 ν ≈ 10⁻⁴ m²/s,Re ≈ 600,从容层流;升温到 50 ℃ 时 ν ≈ 2×10⁻⁵ m²/s,Re ≈ 3000,一脚跨进转捩区。冬天层流的管路,夏天可能悄悄换成湍流的脾气——阻力特性整个换挡。工程纪律由此而来:判流态按「最热工况、最高流速、最小管径」的组合来算——阻力算小了是事故,算大了只是保守;粘度值必须标注温度,脱离温度报粘度,等于报股价不报日期。

6.3层流与湍流的本质区别

层流和湍流不只是「整齐 vs 杂乱」的外观差别,它们有深刻的结构和物理机制差异。

速度剖面

圆管里,由于粘性(贴壁流体被「粘住」,速度为零),截面上的流速分布不均匀:

层流与湍流的速度剖面 层流速度剖面为抛物线,中心速度约为平均2倍;湍流剖面较饱满、中心平坦,壁面附近梯度大。 层流:抛物线剖面 湍流:饱满剖面 中心平坦
图 6-2层流是抛物线速度剖面,湍流是「饱满、中心平坦」的剖面,但壁面附近梯度极大。

层流的抛物线剖面:一次力平衡推导

「层流剖面是抛物线」不是实验拟合,它可以由力平衡严格导出。条件:等直圆管、充分发展(速度沿程不再变化,u 只是半径 r 的函数)、定常。取一个与管同轴、半径 r、长度 L 的圆柱形流体元:

梯度与半径成正比、方向指向减速,积分并代入壁面条件 u(R) = 0(无滑移):

u(r) = Δp/(4μL) · (R² − r²)

这正是抛物线。三个立刻能读出的推论:轴线速度最大 umax = ΔpR²/(4μL);平均速度是它的一半 v̄ = ΔpR²/(8μL) = umax/2——所以层流里测了管心速度就知道流量;壁面切应力 τw = μ·|du/dr|w = ΔpR/(2μL) = 8μv̄/d——阻力的第一现场,6.4 节直接拿去用。

📐 例题 3 油管层流:从轴线速度到压降

输油管 d = 20 mm(R = 10 mm),油 μ = 0.2 Pa·s、ρ = 870 kg/m³(ν = μ/ρ ≈ 2.3×10⁻⁴ m²/s),管内为充分发展层流,测得轴线最大速度 umax = 1.0 m/s。(a) 验证流态确为层流;(b) 求断面平均流速;(c) 求壁面速度梯度与壁面切应力;(d) 求 20 m 管长上的压降。

思路 (a) Re = v̄d/ν;(b) 层流 umax = 2v̄;(c) 壁面梯度 |du/dr|w = 4v̄/R,再乘 μ 得 τw;(d) 由 τw = ΔpR/(2L) 反解,或直接用 Δp/L = 8μv̄/R²。

解 (a) v̄ = umax/2 = 0.5 m/s;Re = 0.5×0.02/2.3×10⁻⁴ ≈ 43.5 ≪ 2300,层流无疑。(b) v̄ = 0.5 m/s。(c) |du/dr|w = 4v̄/R = 4×0.5/0.01 = 200 s⁻¹;τw = 0.2×200 = 40 Pa。(d) Δp/L = 8μv̄/R² = 8×0.2×0.5/0.01² = 8000 Pa/m;20 m 管长 Δp = 160 kPa。复核:τw = ΔpR/(2L) = 160000×0.01/40 = 40 Pa ✓ 与 (c) 一致。

答:Re≈43.5 层流;v̄ = 0.5 m/s;壁面梯度 200 s⁻¹、τw = 40 Pa;20 m 压降 160 kPa。 ⚠ 易错点:①umax = 2v̄ 只对圆管充分发展层流成立,湍流约 1.2 倍,别张冠李戴;②R 与 d 差一倍,公式里该用谁就用谁——本例 |du/dr|w = 4v̄/R 用半径,Δp/L = 8μv̄/R² 也是半径平方。

💡 实用:层流的「测速红利」——只测管心就能定流量 umax = 2v̄ 这个精确关系给现场测量留了个便宜:把探针(皮托管针头或热线探头)只放到管轴中心读一个 umax,流量立刻得手——Q = A·v̄ = (πd²/4)·(umax/2) = πd²·umax/8。一次单点测量顶一次全剖面扫描。但这个红利仅层流有效:湍流剖面 umax/v̄ ≈ 1.17–1.26 且比值随 Re 缓慢漂移,单点读数必须配手册系数修正。这也是工程里「先判流态、再选测量方法」的顺序逻辑——不知道 Re,连该测哪个点都不知道。

为什么湍流阻力大?

湍流阻力远大于层流,原因在于动量传递机制不同:

所以同样流速下,湍流的流动阻力(能量损失)比层流大很多。这也是为什么工程上有时要刻意「保持层流」来省能。

本质区别 层流靠分子粘性传递动量,湍流靠涡团掺混传递动量。后者效率高得多,所以湍流的阻力和能量耗散都远大于层流。

湍流的时均化:把「乱流」拆成平均 + 脉动

湍流虽然乱,但不是没有章法。处理它的标准姿势是时均化:把瞬时速度写成

u(t) = ū + u′(时均值 + 随机脉动)

ū 是足够长时间内的平均,u′ 是围绕它的快速涨落(图 6-3)。脉动的强弱用湍流度衡量(√u′²/ū),风洞品质、大气颠簸、发动机进气畸变,说的都是它。这一拆分让「不可解」的湍流变得可工程化:工程师关心的是平均流速剖面、平均阻力——脉动部分则统计地藏进「湍流切应力」里(6.4 节)。顺带解释了湍流剖面为何「中心平坦、壁面陡峭」:掺混像一把抹刀,把中心区的动量摊匀;而无滑移条件钉死了壁面速度为零,梯度只能挤在极薄的近壁层里完成。

湍流的时均化:瞬时速度 = 时均值 + 脉动 一条围绕水平均值线上下无规则起伏的速度时间序列,虚线为时均速度,起伏部分标注为脉动u′。 u(t) 时间 t 时均速度 ū 瞬时速度 u(t) = ū + u′ 脉动 u′(正负交替、随机) 脉动越强 → 湍流度越高;时均化之后,「乱流」才能进入工程计算。
图 6-3湍流的时均化:瞬时速度拆成「时均值 ū + 脉动 u′」,平均交给计算,脉动交给统计。
层流 vs 湍流:一张表看全
对比项层流湍流
速度剖面抛物线(光滑)饱满、中心平坦,近壁梯度极大
umax/v̄2(精确)≈1.17–1.26(随 Re 缓变)
动量传递者分子(粘性)分子 + 涡团(掺混主导)
切应力规律τ = μ(du/dy) 严格成立分子粘性 + 湍流(雷诺)应力
沿程压降随速度Δp ∝ v̄(线性)Δp ∝ v̄1.75–2(近似平方)
能量损失小大得多(掺混是昂贵的)
工程偏好液压、输油、层流翼型想留住它换热器、燃烧室、搅拌槽想要它

💡 辨析:牛顿 vs 非牛顿——雷诺数判据的适用边界 本章的 τ = μ(du/dy)(μ 为常数)只对牛顿流体成立:水、空气、轻质油都是。但流体世界还有三类「脾气古怪」的成员:剪切变稀(粘度随剪切率增大而下降——油漆、血液:刷得越快越顺滑);剪切变稠(越搅越稠——浓淀粉糊);宾汉体(先要越过「屈服应力」才开始流动——牙膏挤才动、润滑脂能「挂」在轴承上不流走)。对它们,「粘度」变成随剪切率变化的表观粘度,雷诺数判据只能取某一工况的表观粘度近似使用。航空语境下的燃油、液压油、滑油在常温下都近似牛顿流体,可放心使用本章结论;润滑脂例外,它是典型的宾汉体。

6.4粘性切应力与流动阻力

现在把粘性正式纳入受力分析。第 1 章已经给出牛顿内摩擦定律 τ = μ(du/dy),它正是「流动阻力」的微观来源。

从切应力到总阻力

流体沿管道流动时,紧贴壁面的一层流体速度为 0,稍靠里的一层有速度,于是壁面附近有速度梯度 du/dy,产生切应力 τw(壁面切应力)。流体要维持流动,就必须不断克服这个切应力,也就是「摩擦阻力」。

在层流里,τ = μ(du/dy) 是严格的;在湍流里,壁面附近的时均切应力是「分子粘性 + 湍流涡团」两部分之和,比层流大得多。总阻力就是对壁面切应力沿整个壁面积分。

层流的定量账:从 τw 到哈根-泊肃叶公式

层流里这一切都能算到底。6.3 节已导出壁面切应力 τw = 8μv̄/d。对长 L 的管段,把整根液柱当作隔离体再写一次力平衡——压差推力 Δp·(πd²/4) 等于湿壁上的总摩擦力 τw·πdL:

Δp·(πd²/4) = τw·πdL ⟹ Δp = 32μLv̄/d² ⟹ Q = πΔpd⁴/(128μL)

最后一步就是著名的哈根-泊肃叶公式:流量与压差成正比、与管径四次方成正比。两个立即要记住的推论:其一,层流压降与速度线性(Δp ∝ v̄)——速度翻倍、压降翻倍;其二,管径的杠杆极大——管径减半,同样压差下流量只剩 1/16,这就是为什么血管变细一点点,供血就吃紧。注意此式只对层流成立:湍流里 Δp 近似随速度平方增长(下一章的 λ 定量处理),沿用泊肃叶会把损失严重低估。

💡 深挖:湍流的切应力——雷诺应力与涡粘度 把 6.3 的时均化用在动量上,湍流的时均切应力分成两部分:τ = μ(dū/dy) − ρ·u′v′(u′、v′ 为脉动分量的时均值乘积)。第二项就是雷诺应力——涡团横穿流层「搬运」动量产生的等效摩擦,通常远大于分子项。工程上模仿分子粘性的形式,把它写成 ρ·εm(dū/dy),其中 εm 叫涡粘度;水流里涡粘度可达分子运动粘度的成百上千倍,这正是「湍流阻力大」的定量注脚。但 εm 不是物性,它随流动本身变化——湍流理论百年来苦苦攻坚的,就是这一个小小的「不是常数」。

💡 深挖:漏掉的机械能,最后变成了什么? 答案是热——全部、不可逆地变成热。湍流的能量账本叫「能量级串」:大涡从主流汲取能量,逐级传递给小涡,直到小到某个尺度时分子粘性终于接管、把剩余动能直接耗散为内能。给例题 4 的油管算一笔绝热小账:泵功率 145 W 全部变成热,质量流量 2.05 kg/s、油比热约 2000 J/(kg·K),油温每循环只升约 0.04 K——短期无感,但长输管线数百万瓦的摩擦热积累起来就是可观温升,所以长距离输油的水力计算(算压降)与热力计算(算温度、算粘度)必须耦合迭代:温度影响粘度、粘度影响压降、压降影响发热,环环相扣。

💡 预告:近壁区的三层结构 「壁面速度梯度极大」的湍流近壁区,按主导力量还能细分:紧贴壁面的粘性底层——涡团被壁面压得抬不起头,分子粘性说了算,τ = μ(du/dy) 仍然成立;往外是缓冲层——两种机制犬牙交错;再往外对数律区——涡团掺混完全主导,速度呈对数分布。下一章的沿程阻力系数 λ(穆迪图)为什么是那个形状、光滑管与粗糙管的分野在哪,物理答案全部埋在这三层里。本章先把「分子粘性只在极薄一层当家」的图景记牢即可。

📐 例题 4 输油管层流全链条:切应力 → 阻力 → 压降 → 泵功率

输油管 d = 50 mm、长 L = 20 m,油 μ = 0.2 Pa·s、ρ = 870 kg/m³,平均流速 v̄ = 1.2 m/s(例题 2 已判明 Re ≈ 261,层流)。(a) 求壁面切应力;(b) 分别用「切应力 × 湿壁面积」与「压差 × 断面」两条路求管壁总阻力;(c) 求 20 m 管长的压降;(d) 泵需要多大功率;(e) 折算成多少米油柱的水头损失。

思路 (a) τw = 8μv̄/d;(b) 两条路径应一致(力平衡自检);(c) Δp = 32μLv̄/d²;(d) P = Δp·Q;(e) hf = Δp/(ρg)。

解 (a) τw = 8×0.2×1.2/0.05 = 38.4 Pa。(b) 路径一:F = τw·πdL = 38.4×π×0.05×20 ≈ 120.6 N。路径二:Δp = 32μLv̄/d² = 32×0.2×20×1.2/0.05² = 61440 Pa;断面 A = π/4×0.05² = 1.96×10⁻³ m²;F = Δp·A = 61440×1.96×10⁻³ ≈ 120.6 N ✓ 两条路吻合。(c) Δp = 61.4 kPa。(d) Q = v̄A = 1.2×1.96×10⁻³ = 2.36×10⁻³ m³/s(≈141 L/min);P = Δp·Q = 61440×2.36×10⁻³ ≈ 145 W。(e) hf = 61440/(870×9.8) ≈ 7.2 m 油柱。

答:τw = 38.4 Pa;总阻力 ≈ 120.6 N;Δp ≈ 61.4 kPa;泵功率 ≈ 145 W;水头损失 ≈ 7.2 m 油柱。 ⚠ 易错点:①此题全程层流,泊肃叶公式可用;若把 Re≈261 误判成湍流、套平方规律,压降会算错好几倍;②功率是 Δp×Q 而不是 τ×Q——切应力作用在壁面,泵功率结算的是断面上的压差做功。

📖 人物:哈根与泊肃叶——素不相识的同题共振 1839 年,德国水利工程师哈根用毛细管做实验,发现流量「与压差成正比、与管径四次方反常地敏感」;次年,法国医生泊肃叶(Jean Léonard Marie Poiseuille)为了研究血压与血流,独立地用更精密的毛细管装置得到了同样规律——他原本想当数学家,因病转攻医学,却以一条流体力学公式名留青史。两人素不相识、动机各异(一个修水管、一个量血压),却在同一年摸到同一条四次方定律——科学史上反复上演的「同题共振」。今天这条公式前并排挂着两个人的名字:哈根-泊肃叶。顺带一提:血液在毛细血管里正是 Re ≪ 1 的爬流,泊肃叶医生当年选的课题,本身就是流体力学与医学的一次握手。

工程应用:管线的「粘度经济学」

例题 4 的链条延伸出去就是真实的输油工程:层流压降与粘度成正比,所以「降粘」直接省钱——原油加热站的存在意义就在于此(把 μ 从 0.4 降到 0.2,泵功率立省一半)。但若管线里是湍流,加热的收益就小得多(湍流压降对 μ 不敏感)——投不投加热炉,先算 Re。同样的账在液压系统里倒着算:液压油粘度不能一味选小,因为太稀容易泄漏、太稠则损失大,粘度与流态的选择是液压设计的核心权衡之一。粘度、流速、管径、温度四个旋钮,转的都是雷诺数这把锁。

把视野从管道抬到飞行中的机翼上:贴着翼面有一层几毫米到几十毫米厚的边界层——它就是管流近壁区在曲面上的翻版。前缘附近雷诺数低(离前缘距离 x 还小),边界层是层流,摩擦阻力小;往后 Rex 增大,流动转捩成湍流,摩擦阻力陡增。设计者的两难在于:层流阻力小但「体质弱」——逆压梯度一来容易分离;湍流阻力大但「身板硬」——高动量的掺混能顶住逆压、推迟分离。想让机翼「滑」(推迟转捩、保层流)还是想要「抗」(早转捩、防分离),取决于飞行阶段更怕阻力还是更怕失速——高亚声速巡航的运输类飞行器普遍把推迟转捩当作减阻的重要抓手,因为巡航阻力中摩擦阻力占大头,省下百分之一都是真金白银的燃油。这个权衡,正是本章「流态决定阻力」在航空工程里最贵的一次应用。

有粘性后,能量会「漏」

这正是第 4 章留下的伏笔:真实流体有粘性,粘性切应力做「负功」,把机械能不可逆地转成热能。所以沿流动方向,总机械能(总水头)会不断下降。这个「下降」就是能量损失(水头损失),是下一章管内流动计算的核心对象。

总水头线沿程下降:能量「漏」的几何表达 理想流体总水头线是水平线;粘性流体的总水头线沿程倾斜下降,与测压管水头线平行,两者间竖直距离为速度水头。 管道(等径) 理想流体:总水头线水平(能量不漏) 真实流体总水头线 H:沿程下降 测压管水头线(与 H 平行,间距 = 速度水头) 水头损失 h_f 断面(1) 断面(2) 粘性切应力持续做负功 → 机械能沿流不可逆地漏成热。
图 6-4水头线语言下的能量损失:理想流体的总水头线水平;真实流体的总水头线一路下斜,两端的高差就是 hf。

至此,第 4 章的「理想账本」补上了最后一栏:伯努利方程算的是能量怎么转化,本章讲清了能量为什么漏、漏多快——壁面切应力(层流 8μv̄/d、湍流还要乘涡团的账)是漏水的孔径,雷诺数决定孔径的大小量级。下一章管内流动计算,就是把 hf 精确到每个管段、每个弯头,工程上每一次「泵选大了一号」或「管线压不够」,背后都是这笔账的得与失。

🧠 一句话记忆 粘性切应力 τ = μ(du/dy) 是阻力的微观来源;有粘性就有水头损失,机械能沿流不断「漏」成热。

λ = 64/Re:层流的『出厂设定』 圆管层流的摩擦系数不用查表,理论上直接锁死:λ = 64/Re。Re = 500 时 λ ≈ 0.128,Re = 2000 时 λ ≈ 0.032——雷诺数越大,粘性相对越弱,摩擦系数越小。更妙的是把它代回第 7 章预告的达西公式 h_f = λ(L/d)(v²/2g):Re 中的 v 把 v² 抵消掉一个,层流的水头损失与流速成正比。实算一条百米长、内径 50 mm 的油管(ν = 10⁻⁴ m²/s):流速 1 m/s 损失约 13.0 m 油柱,2 m/s 时恰为 26.1 m——翻倍即翻倍;而湍流段 λ 近似常数,损失按 v² 猛涨。低速系统的「省电红利」正源于这条一次方律,长输管线刻意保温降粘、把流动按在层流段,就是在给这条定律发工资。

本章小结

  1. 雷诺实验揭示了两种流态:层流(分层、互不掺混)与湍流(杂乱、剧烈掺混)。
  2. 雷诺数 Re = ρvL/μ = vL/ν,其物理意义是惯性力与粘性力之比。
  3. Re 小(粘性主导)→ 层流;Re 大(惯性主导)→ 湍流;圆管临界 Re ≈ 2300。
  4. 层流速度剖面是抛物线,湍流剖面饱满;湍流因涡团掺混传递动量,阻力和能耗远大于层流。
  5. 粘性切应力 τ = μ(du/dy) 是流动阻力的微观来源,粘性导致机械能沿流不可逆地耗散为热(水头损失)。

读完本章,你应该能……

思考题

1. 雷诺数为什么能「一眼判断」流动是层流还是湍流?它的分子、分母各代表什么?

提示:Re = ρvL/μ,分子量级是惯性力(促使扰动放大、失稳),分母量级是粘性力(抑制扰动、稳定流动)。两者之比决定谁占上风。

2. 同样的水流,在细管里可能是层流,在粗管里却是湍流。用雷诺数解释。

提示:Re ∝ L(管径 d)。管径增大,Re 增大,超过临界值就变成湍流。流动状态由 Re 决定,而不是单看流速。

3. 为什么湍流的能量损失比层流大得多?

提示:层流动量靠分子粘性传递(低效);湍流靠涡团大尺度掺混传递动量(高效),产生远大于分子粘性的等效切应力,耗散更多能量。

4. 雷诺数是一个无量纲数,这个性质有什么实际意义?

提示:无量纲意味着它与单位制无关、普适;更重要的是,几何相似的流动只要 Re 相同,流态就相似——这是第 8 章相似原理和模型实验的基础。

5. 用同一根吸管分别吸水和吸蜂蜜,为什么吸蜂蜜费力得多?这笔「力」是惯性还是粘性主导?

提示:蜂蜜 μ ≈ 10 Pa·s(ρ≈1400 → ν≈7×10⁻³ m²/s),管径几毫米、流速几厘米每秒 → Re ≈ 0.04 ≪ 1,粘性力碾压惯性,费力的大头是粘性剪切;水在同样条件下 Re 可达数千,惯性占比高得多。「浓稠 = 难吸」的本质是 Re 极小、粘性做功。

6. 为什么湍流剖面的中心区平坦,而近壁区的速度梯度反而比层流更陡?

提示:中心区——涡团掺混把动量摊匀;近壁区——无滑移条件钉死壁面速度为零,而湍流平均速度更高(剖面饱满),「从 0 拉到接近 umax」只能在极薄的近壁层里完成,梯度必然更陡。

7. 油漆的「表观粘度」是什么?对它使用雷诺数判据要注意什么?

提示:油漆是剪切变稀流体,表观粘度 μapp = τ/(du/dy) 随剪切率变化——刷得快就变稀。用雷诺数时只能取工作点附近的 μapp 近似计算,结论只在该工况附近有效;工况一变(如喷涂 vs 静置),粘度就换了值。

8. 长输原油管线为什么要设加热站?为什么流量大增后加热的「省钱效果」反而变差?

提示:低速层流时 Δp ∝ μv̄(泊肃叶),加热降粘收益与降粘比例直接挂钩;流量增大进入湍流后,压降近似 ∝ v̄² 且对 μ 不敏感,同样加热换不回多少压降——所以加热站的经济性随流态改变,设计时要按运行工况的 Re 逐段算。

9. 血管里的血流是层流还是湍流?剧烈运动时听到的心脏杂音说明什么?

提示:主动脉 Re ≈ 2000,平时是层流;剧烈运动时心输出量可翻几倍,Re 越过临界,瓣膜后方出现湍流区——杂音正是湍流脉动引起的振动。用流态「听诊」,是本章内容在临床中最直观的应用之一。同样的道理也适用于动脉狭窄段:管径变细、局部流速升高,Re 与脉动同步增大——狭窄后方的杂音强度正是医生判断狭窄程度的听诊线索之一。

10. 内径 100 mm 的管道,流速 2 m/s。(a) 20 ℃ 水(ν = 10⁻⁶ m²/s)的 Re 是多少?(b) 换成 ν = 10⁻⁴ m²/s 的冷油呢?各按什么流态处理?

提示:Re = vd/ν。(a) 2×0.10/10⁻⁶ = 200000,铁定湍流;(b) 2×0.10/10⁻⁴ = 2000,正压在临界线上——按「宁可保守」的原则当湍流处理。同一条管、同一流速,粘度差百倍,流态是两个世界。

11. 层流中 Re 从 800 升到 1600,λ = 64/Re 怎么变?层流的水头损失又随流速怎么涨?两个结论矛盾吗?

提示:λ 从 0.080 降到 0.040,减半;但 Re 翻倍源于流速翻倍,而层流损失 h_f = λ(L/d)(v²/2g) 代入 λ = 64/Re 后与 v 成正比——流速翻倍,损失照样翻倍。「系数变小」与「账单变大」并存不矛盾:λ 是折扣率,v²/2g 是账单基数,折扣变浅、基数翻倍,总账自然翻倍。