力学 · 从零学起 第 2 / 5 节

第 2 章 · 粘性流体与雷诺数:一个参数统治一切

加入粘性项后方程变成纳维-斯托克斯方程,而雷诺数一个无量纲数立刻接管了流态的分类权——层流还是湍流、相似不相似,都由它说了算。

纳维-斯托克斯方程

$$\partial_t\mathbf{v} + (\mathbf{v}\cdot\nabla)\mathbf{v} = -\frac{1}{\rho}\nabla p + \frac{\eta}{\rho}\nabla^2\mathbf{v}$$

粘性系数 $\eta$(动力)与 $\nu = \eta/\rho$(运动)把动量扩散进来。量纲分析产出雷诺数:

$$Re = \frac{UL}{\nu} = \frac{\text{惯性}}{\text{粘性}}$$

$Re$ 小:粘性统治,流动温顺(斯托克斯流);$Re$ 大:惯性统治,粘性被挤到边界层与湍流里。同 $Re$ 的流动彼此相似——风洞缩比实验的合法性证书。

泊肃叶与斯托克斯:两个可解标杆

两者共同的教训:可解析解的 N-S 问题屈指可数,每一个都被反反复复使用。

  • 圆管层流(泊肃叶):$Q = \frac{\pi R^4}{8\eta}\frac{\Delta p}{L}$——流量正比半径四次方,血管、管道工程的定量基础;
  • 低 $Re$ 绕流(斯托克斯阻力):$F = 6\pi\eta a v$——小颗粒沉降、胶体、微流控的公式。

涡旋的粘性演化

涡量方程 $\partial_t\boldsymbol\omega = \nu\nabla^2\boldsymbol\omega$(无膨胀项时)是扩散方程——粘性让涡量向外弥散。突然启动的平板带来"扩散时间 $t \sim L^2/\nu$"的基本估计:粘性影响的长、慢、广,全在这个平方标度里。

原因要点:第 2 章给出应力张量的完整推导($\sigma_{ik} = -p\delta_{ik} + \eta(\partial_i v_k + \partial_k v_i) + \zeta\delta_{ik}\nabla\cdot\mathbf{v}$),两个粘性系数(剪切 $\eta$、体积 $\zeta$)来自各向同性的最一般线性本构;相似律一节把无量纲化作为方法论正式立档。

老师的提醒

  • 任何流动问题开工先算 $Re$:它决定你该用哪个近似(斯托克斯/边界层/湍流模型)、该查哪一章。
  • 泊肃叶的 $R^4$ 定律的工程直觉:管道半径减半,泵功涨 16 倍——设计里最贵的平方立方惩罚之一。

自检

  • 雷诺数的物理含义与相似律?
  • 推导泊肃叶流的速度剖面与流量公式。
  • 涡量扩散方程说明粘性的什么作用?

更深一层:雷诺数是"相似性的外交护照"

雷诺数 $Re = UL/\nu$ 的深层身份是无量纲化的外交成果:把方程无量纲化后,纳维-斯托克斯方程只含一个参数 $Re$——$Re$ 相同的两个流动,解(无量纲形式)完全相同。这张外交护照的工程价值:风洞里 1:10 模型的结果可以外推到真机(只要 $Re$ 匹配——大型风洞为凑 $Re$ 而造的原因)、管道阻力系数表对一切流体通用、微流控的 $Re \sim 1$ 世界与航空的 $Re \sim 10^8$ 世界遵循不同物理。无量纲数是流体力学的外交语言,每个新无量纲数(马赫数——压缩性、弗劳德数——重力、普朗特数——热扩散比动量扩散)都是一份新的相似性条约。

泊肃叶与斯托克斯两个"可解孤岛"的工程地位也值得标注:血管流动(泊肃叶 + Womersley 脉动修正)、微流控芯片(层流扩散混合设计)、气溶胶沉降(斯托克斯阻力 + Cunningham 滑移修正)——低 $Re$ 世界是"解析解统治的世界"。

推导细读:涡量扩散方程——"粘性是涡的稀释剂"

涡量输运方程(不可压二维)$\partial_t\boldsymbol\omega + (\mathbf{v}\cdot\nabla)\boldsymbol\omega = \nu\nabla^2\boldsymbol\omega$ 的读法:涡量像热量一样扩散。突然启动的平板(Stokes 第一问题)在流体中制造的涡量层厚度 $\delta \sim \sqrt{\nu t}$——一秒后空气中约 4 mm、水中约 1 mm。$\sqrt{\nu t}$ 的扩散标度同时控制边界层生长、涡环衰减、刹车盘的热渗透——扩散问题全家共用一条标度律。实验对应:突然拔出板上的奶油,几秒后奶油痕迹的厚度就是 $\sqrt{\nu t}$ 的肉眼测量。

在真实世界里:雷诺数的"选型地图"

按 $Re$ 分区的工程选型地图:$Re \ll 1$ → 斯托克斯流(微机电、气溶胶、血细胞);$Re \sim 1\text{–}1000$ → 中间区(数值为主,实验辅助);$Re \gg 1$ → 边界层 + 位势流的分区世界(翼型、管道湍流、船体)。同一个几何形状在不同 $Re$ 区的阻力系数可以差 100 倍(圆柱:低 $Re$ 的 $\propto 1/Re$ 到高 $Re$ 的 $\approx 1$)——这就是"先算 $Re$ 再查表"的纪律来源。

专题深入:雷诺数是"相似性的外交护照"完整版

雷诺数 $Re = UL/\nu$ 的深层身份是无量纲化的外交成果:把方程无量纲化后,纳维-斯托克斯方程只含一个参数 $Re$——$Re$ 相同的两个流动,解(无量纲形式)完全相同。这张外交护照的工程价值:风洞里 1:10 模型的结果可以外推到真机(只要 $Re$ 匹配——大型风洞为凑 $Re$ 而造的原因)、管道阻力系数表对一切流体通用、微流控的 $Re \sim 1$ 世界与航空的 $Re \sim 10^8$ 世界遵循不同物理。无量纲数是流体力学的外交语言,每个新无量纲数(马赫数——压缩性、弗劳德数——重力、普朗特数——热扩散比动量扩散)都是一份新的相似性条约。

泊肃叶与斯托克斯两个"可解孤岛"的工程地位也值得标注:血管流动(泊肃叶 + Womersley 脉动修正)、微流控芯片(层流扩散混合设计)、气溶胶沉降(斯托克斯阻力 + Cunningham 滑移修正)——低 $Re$ 世界是"解析解统治的世界"。

数学细节展开:涡量扩散方程——"粘性是涡的稀释剂"完整版

涡量输运方程(不可压二维)$\partial_t\boldsymbol\omega + (\mathbf{v}\cdot\nabla)\boldsymbol\omega = \nu\nabla^2\boldsymbol\omega$ 的读法:涡量像热量一样扩散。突然启动的平板(Stokes 第一问题)在流体中制造的涡量层厚度 $\delta \sim \sqrt{\nu t}$——一秒后空气中约 4 mm、水中约 1 mm。$\sqrt{\nu t}$ 的扩散标度同时控制边界层生长、涡环衰减、刹车盘的热渗透——扩散问题全家共用一条标度律。实验对应:突然拔出板上的奶油,几秒后奶油痕迹的厚度就是 $\sqrt{\nu t}$ 的肉眼测量。

工程案例细节:雷诺数的"选型地图"工程版

按 $Re$ 分区的工程选型地图:$Re \ll 1$ → 斯托克斯流(微机电、气溶胶、血细胞);$Re \sim 1\text{–}1000$ → 中间区(数值为主,实验辅助);$Re \gg 1$ → 边界层 + 位势流的分区世界(翼型、管道湍流、船体)。同一个几何形状在不同 $Re$ 区的阻力系数可以差 100 倍(圆柱:低 $Re$ 的 $\propto 1/Re$ 到高 $Re$ 的 $\approx 1$)——这就是"先算 $Re$ 再查表"的纪律来源。

例题

水($\nu = 10^{-6}\ \mathrm{m^2/s}$)以 $v = 1\ \mathrm{m/s}$ 流过直径 $D = 5\ \mathrm{cm}$ 的管道。求雷诺数并判断流态;若换成同尺寸的甘油($\nu = 1000\times$ 水)呢?

解答

$$Re = \frac{vD}{\nu} = \frac{1\times0.05}{10^{-6}} = 5\times10^4$$

管流临界 $Re \approx 2300$ ⟹ 充分发展湍流。

甘油:$Re = 50 < 2300$ ⟹ 层流——泊肃叶公式直接可用(同管道压差损失比湍流可预测得多)。

  • 一句话记住:粘性越大/速度越低/尺寸越小 → 越层流。所有流动问题的第一笔账就是这一个除法。

本页是《流体力学(理论物理教程第六卷)》第 2 节的学习笔记:为便于初学,讲解、例题与要点概括均为本站原创整理,公式与章节定位供对照原书使用。

📑 本书目录 经典研读