📖 人物:奥斯本·雷诺(Osborne Reynolds) 雷诺是 19 世纪英国的工程师兼物理学家,是「把流体力学带进工程」的关键人物之一。他不仅发现了两种流态,还提出了那个以他命名的判据——雷诺数,一举解决了「什么条件下是层流、什么条件下是湍流」这个困扰工程师已久的问题。
🎯 本章学习目标
- 理解粘性流体流动中存在的两种流态:层流与湍流。
- 掌握雷诺数的定义,并深刻理解它的物理意义——惯性力与粘性力之比。
- 知道临界雷诺数,会判断圆管流动是层流还是湍流。
- 理解层流与湍流的内部结构差异,以及湍流为何阻力大。
- 把粘性切应力(牛顿内摩擦)重新放回流动分析中,理解「阻力」从哪来。
6.1雷诺实验:两种流态的发现
1883 年,英国物理学家奥斯本·雷诺做了一个后来载入史册的实验:让水流过一根透明玻璃管,并在管入口处滴入染色液体,观察染色线的形态。他发现了两种截然不同的流动状态。
- 流速较小时:染色线是一条稳定、笔直的细线,说明流体分层流动、各层互不掺混——这是层流。
- 流速增大:染色线开始波浪状摆动,进入过渡状态。
- 流速再增大:染色线迅速散开、与周围水搅浑,说明流体杂乱无章、剧烈掺混——这是湍流(紊流)。
实验设计的三个门道:为什么雷诺能「看见」流态
这个实验看似简单,实则处处是功夫,门道至少有三条:
- 染料要与水「等密度」:染料密度若比水大,会自己下沉、把实验搅浑;雷诺特意调配染料密度与水接近,保证看到的运动全来自「水带着染料走」,而不是「染料自己走」。
- 入口要顺、上游要稳:水箱保持恒定水位(定常),管入口做成流线喇叭口,避免装置自己制造扰动——这样观察到的失稳才来自流动本身,而不是设备的杂音。
- 染色线的身份是「示踪迹线」:定常流动中迹线与流线重合,一条染色线就是一束可见的流线。层流里它笔直清晰——分层不掺混;湍流里它被涡团反复折叠、拉伸、撕碎——掺混的直接影像证据。用「示踪法」而非「测压法」研究流态是雷诺的高明之处:压差反映结果,示踪反映结构。
先纠一个直觉偏差:层流不等于「慢」
日常语言里「水流很急」几乎等于湍流,这不算全错,但工程上必须精确:流态不取决于速度的绝对大小,只取决于雷诺数(下一节)。同样的 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(τ × 面积)。它抑制扰动、把流动「拉平」。
两者之比:
惯性力 / 粘性力 ≈ (ρv²L²) / (μvL) = ρvL/μ = Re
量纲推导:这两个「力」是怎么估出来的
上面的量级不是拍脑袋,各推一步就明白。取边长为 L 的流体微团,速度 v,特征时间 t ≈ L/v(流过自身尺度的时间):
- 惯性项:质量 m ≈ ρL³;加速度量级 a ≈ v/t = v²/L。惯性力 Fin ≈ ρL³·(v²/L) = ρv²L²。
- 粘性项:牛顿内摩擦 τ = μ(du/dy),速度在 L 的尺度内变化 v,故 du/dy ≈ v/L;作用面积 ≈ L²。粘性力 Fvis ≈ μ(v/L)·L² = μvL。
- 相除:Fin/Fvis = ρv²L²/(μvL) = ρvL/μ = Re。雷诺数不是被硬性定义出来的,它是从「两个力谁大谁小」这个问题里长出来的。
有了这层出身,就能理解为什么 Re 小时流动「稳」:粘性相对强大,任何扰动刚冒头就被抹平,像果冻里戳一下立刻复原;Re 大时惯性独大,扰动被流动放大成涡团、滚雪球成湍流——像推倒一排多米诺。
💡 细节:特征长度 L 怎么选? Re 里的 L 取「最能代表该流动的几何尺度」:圆管取内径 d(不是半径——临界值 2300 就是对 d 定义的,换成半径数字就全变了);沿平板流动取离前缘的距离 x(边界层沿程增厚,Rex 沿程增大);非圆截面用当量直径 de = 4A/χ(A 为过流面积,χ 为湿周)。原则一句话:取「速度梯度发生方向上的代表性尺度」。用错特征长度,Re 的数值就失去可比性,临界值也跟着作废。
核心理解 雷诺数 = 惯性力 / 粘性力。Re 小,粘性力占上风,扰动被「按住」,流动保持层流;Re 大,惯性力占上风,扰动被「放大」,流动失稳变成湍流。所以雷诺数刻画的是两种力量的此消彼长。
这就是雷诺数的魅力:它用一个数字,把「流速、尺度、粘度」对流动状态的影响打包成一个判据。不管是一根细管里慢流的水,还是一条大河里奔腾的水,只要 Re 相同,流动的「性格」就相似(这正是第 8 章相似原理的思想源头)。
临界雷诺数与判断
对圆管流动,实验给出:
- Re < 2300:层流。
- 2300 < Re < 4000:过渡区(不稳定)。
- Re > 4000:湍流。
工程上常取 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/s | 10⁴–10⁵ | 日常用水几乎总是湍流 |
| 大江大河 | 水深≈5 m,v≈2 m/s | 10⁶–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 倍),从中心到壁面平滑递减。
- 湍流:因为强烈掺混,中心区的速度被「抹平」了,速度剖面比较饱满、中心区较平坦,但紧贴壁面处速度梯度极大(速度在极薄一层内从很大骤降到 0)。
层流的抛物线剖面:一次力平衡推导
「层流剖面是抛物线」不是实验拟合,它可以由力平衡严格导出。条件:等直圆管、充分发展(速度沿程不再变化,u 只是半径 r 的函数)、定常。取一个与管同轴、半径 r、长度 L 的圆柱形流体元:
- 推动它的力:两端压差 × 端面积 = Δp·πr²。
- 阻碍它的力:圆柱侧表面上的粘性力 = 2πrL·μ·|du/dr|。
- 定常 = 不加速:两力平衡 → Δp·πr² = 2πrLμ·|du/dr|,化简得 du/dr = −Δp·r/(2μL)。
梯度与半径成正比、方向指向减速,积分并代入壁面条件 u(R) = 0(无滑移):
这正是抛物线。三个立刻能读出的推论:轴线速度最大 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′ 是围绕它的快速涨落(图 6-3)。脉动的强弱用湍流度衡量(√u′²/ū),风洞品质、大气颠簸、发动机进气畸变,说的都是它。这一拆分让「不可解」的湍流变得可工程化:工程师关心的是平均流速剖面、平均阻力——脉动部分则统计地藏进「湍流切应力」里(6.4 节)。顺带解释了湍流剖面为何「中心平坦、壁面陡峭」:掺混像一把抹刀,把中心区的动量摊匀;而无滑移条件钉死了壁面速度为零,梯度只能挤在极薄的近壁层里完成。
| 对比项 | 层流 | 湍流 |
|---|---|---|
| 速度剖面 | 抛物线(光滑) | 饱满、中心平坦,近壁梯度极大 |
| 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 ∝ 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 章留下的伏笔:真实流体有粘性,粘性切应力做「负功」,把机械能不可逆地转成热能。所以沿流动方向,总机械能(总水头)会不断下降。这个「下降」就是能量损失(水头损失),是下一章管内流动计算的核心对象。
至此,第 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² 猛涨。低速系统的「省电红利」正源于这条一次方律,长输管线刻意保温降粘、把流动按在层流段,就是在给这条定律发工资。
本章小结
- 雷诺实验揭示了两种流态:层流(分层、互不掺混)与湍流(杂乱、剧烈掺混)。
- 雷诺数 Re = ρvL/μ = vL/ν,其物理意义是惯性力与粘性力之比。
- Re 小(粘性主导)→ 层流;Re 大(惯性主导)→ 湍流;圆管临界 Re ≈ 2300。
- 层流速度剖面是抛物线,湍流剖面饱满;湍流因涡团掺混传递动量,阻力和能耗远大于层流。
- 粘性切应力 τ = μ(du/dy) 是流动阻力的微观来源,粘性导致机械能沿流不可逆地耗散为热(水头损失)。
读完本章,你应该能……
- 描述雷诺实验,说出层流、过渡、湍流三种状态的特征。
- 默写雷诺数定义,并用自己的话解释「惯性力/粘性力」的物理意义。
- 用 Re ≈ 2300 判断圆管流动的流态。
- 画出并对比层流与湍流的速度剖面。
- 解释为什么湍流的流动阻力远大于层流。
- 推导圆管层流抛物线剖面,记住 umax = 2v̄ 与 τw = 8μv̄/d。
- 用哈根-泊肃叶公式算层流管流的流量与压降,并先核对 Re。
- 说出层流与湍流在剖面、动量传递、损失规律上的三大差异。
- 会估算水管、油管、血管、机翼等典型流动的雷诺数量级,先判量级、再动笔细算。
思考题
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 是账单基数,折扣变浅、基数翻倍,总账自然翻倍。