第2章：无量纲数与相似准则

在流体力学中，无量纲数就像是流动现象的"指纹"——它们能够跨越尺度、跨越流体种类，揭示流动的本质特征。想象一下，为什么鲸鱼游泳的姿态与小鱼截然不同？为什么昆虫能在空中悬停而鸟类不能？为什么风洞中的缩比模型能预测真实飞机的性能？答案都隐藏在几个关键的无量纲数中。本章将介绍流体力学中最重要的无量纲数，让你能够快速判断流动的类型，估算流动特征，并理解相似准则在工程实践中的应用。

2.1 雷诺数：层流与湍流的分界

雷诺数的定义与物理意义

雷诺数（Reynolds number, $Re$）是流体力学中最基础也是最重要的无量纲数，它表征了惯性力与粘性力的相对大小：

$$Re = \frac{\rho V L}{\mu} = \frac{VL}{\nu}$$ 其中：

• $\rho$ - 流体密度
• $V$ - 特征速度
• $L$ - 特征长度
• $\mu$ - 动力粘度
• $\nu = \mu/\rho$ - 运动粘度

物理直觉：

• $Re \ll 1$：粘性主导，流动像蜂蜜一样缓慢粘稠
• $Re \sim 1$：惯性与粘性相当
• $Re \gg 1$：惯性主导，流动容易失稳产生湍流

日常生活中的雷诺数

让我们通过具体例子来建立对雷诺数量级的直觉：

微观世界（$Re < 1$）：

• 精子游动：$Re \sim 10^{-2}$
• 细菌运动：$Re \sim 10^{-5}$
• 在这个世界里，没有惯性，停止划动就立即停止

低雷诺数流动（$Re \sim 1-100$）：

• 蜂蜜从勺子滴落：$Re \sim 0.1$
• 雨滴下落：$Re \sim 100-1000$（取决于大小）

中等雷诺数（$Re \sim 10^3-10^5$）：

• 人游泳：$Re \sim 10^6$（基于身长）
• 鸟类飞行：$Re \sim 10^4-10^5$
• 汽车行驶：$Re \sim 10^6-10^7$

高雷诺数（$Re > 10^6$）：

• 商用飞机：$Re \sim 10^7-10^8$
• 潜艇航行：$Re \sim 10^8-10^9$

临界雷诺数与转捩

不同几何形状的流动有不同的临界雷诺数：

圆管内流动：

• $Re < 2300$：层流
• $2300 < Re < 4000$：过渡区
• $Re > 4000$：湍流

平板边界层：

• $Re_x < 5 \times 10^5$：层流（基于距前缘距离）
• $Re_x > 5 \times 10^5$：开始转捩

圆柱绕流：

• $Re < 5$：无分离的爬流
• $5 < Re < 40$：对称的分离涡
• $40 < Re < 200$：层流涡街
• $200 < Re < 3 \times 10^5$：湍流涡街
• $Re > 3 \times 10^5$：边界层转捩，阻力危机

        层流                过渡                 湍流
    ============        ~~~~~~~~~~~~        ≈≈≈≈≈≈≈≈≈≈≈≈
    平滑有序            不稳定振荡           混沌随机
    可预测              间歇性              统计描述
    低混合              中等混合             强混合

雷诺数的工程应用

管道设计：

• 层流：摩擦系数 $f = 64/Re$
• 湍流：需要用Moody图或Colebrook公式

换热器设计：

• 湍流换热系数比层流高10-100倍
• 但压降也相应增加

减阻设计：

• 游泳衣的疏水表面维持层流
• 高尔夫球凹坑促进转捩降低压差阻力

2.2 马赫数：可压缩性的度量

声速与马赫数定义

马赫数（Mach number, $Ma$）定义为流速与当地声速之比： $$Ma = \frac{V}{a}$$ 其中声速 $a = \sqrt{\gamma RT}$（理想气体）

• $\gamma$ - 比热比（空气约1.4）
• $R$ - 气体常数
• $T$ - 绝对温度

标准大气条件下的声速：

• 海平面（15°C）：340 m/s
• 10 km高空（-50°C）：295 m/s

流动分类

根据马赫数，流动可分为：

不可压缩流（$Ma < 0.3$）：

• 密度变化 < 5%
• 汽车、火车、低速飞机
• 日常生活中的大部分流动

亚声速流（$0.3 < Ma < 0.8$）：

• 需考虑压缩性效应
• 商用客机巡航（$Ma \sim 0.8$）

跨声速流（$0.8 < Ma < 1.2$）：

• 局部出现超声速区
• 激波开始形成
• 波阻急剧增加

超声速流（$1.2 < Ma < 5$）：

• 激波成为主导特征
• 战斗机、超音速客机

高超声速流（$Ma > 5$）：

• 真实气体效应
• 弹道导弹、航天器再入

生活中的马赫数效应

鞭子抽打声： 鞭梢速度可超过声速（$Ma > 1$），产生的"啪"声实际上是微型音爆。

高铁过隧道：

• 列车速度：300 km/h = 83 m/s
• $Ma = 83/340 = 0.24$
• 虽然是低马赫数，但在隧道中压力波传播会产生"隧道爆破"现象

子弹飞行：

• 手枪子弹：$Ma \sim 1.0-1.5$
• 步枪子弹：$Ma \sim 2.0-3.0$
• 产生特征的激波锥

压缩性效应的估算

压力系数修正（Prandtl-Glauert）： $$C_p = \frac{C_{p,0}}{\sqrt{1-Ma^2}}$$ （适用于 $Ma < 0.7$）

阻力发散： 临界马赫数附近，波阻急剧增加：

• 典型翼型：$Ma_{crit} \sim 0.7-0.8$
• 超临界翼型：$Ma_{crit} \sim 0.8-0.85$

2.3 斯特劳哈尔数：涡脱落频率

涡街现象与斯特劳哈尔数

斯特劳哈尔数（Strouhal number, $St$）描述了非定常流动中的特征频率： $$St = \frac{fL}{V}$$ 其中：

• $f$ - 特征频率（如涡脱落频率）
• $L$ - 特征长度
• $V$ - 流速

对于圆柱绕流，在很宽的雷诺数范围内（$200 < Re < 10^5$），斯特劳哈尔数近似为常数： $$St \approx 0.2$$ 这意味着涡脱落频率与流速成正比： $$f = 0.2 \frac{V}{D}$$

日常生活中的涡街

电线在风中的歌唱：

• 10 mm直径电线，风速 10 m/s
• 频率：$f = 0.2 \times 10/0.01 = 200$ Hz
• 正好在人耳可听范围内

烟囱涡激振动：

• 高烟囱必须考虑涡激振动
• 当涡脱落频率接近结构固有频率时会发生共振
• 工程对策：螺旋条纹破坏涡的相关性

潜望镜尾流：

• 潜艇潜望镜伸出水面会产生特征涡街
• 可被反潜机探测
• 对策：非圆形截面设计

斯特劳哈尔数的变化规律

St
^
|
0.3 |    _______________
    |   /               \___
0.2 |  /                    \________ ≈ 0.2
    | /                              
0.1 |/
    |__________________________|_____> log(Re)
    10²  10³   10⁴   10⁵   10⁶

不同形状的斯特劳哈尔数：

• 圆柱：$St \approx 0.2$
• 方柱：$St \approx 0.13$
• 平板（垂直来流）：$St \approx 0.15$

工程应用

流量计设计： 涡街流量计利用稳定的斯特劳哈尔数关系：

• 测量涡脱落频率
• 反算流速：$V = f \cdot D / St$
• 精度可达 ±1%

桥梁设计：

• 塔科马大桥事故（1940）：涡激振动导致
• 现代斜拉桥拉索采用螺旋线、凹坑等措施抑制涡街

高层建筑：

• 摩天大楼顶部的涡激摆动
• 上海中心大厦的螺旋造型减少涡激力24%

2.4 相似准则在模型实验中的应用

三种相似性

几何相似： 模型与原型的所有对应线性尺寸成比例： $$\frac{L_m}{L_p} = \lambda$$ 运动相似： 对应点的速度方向相同，大小成比例： $$\frac{V_m}{V_p} = \frac{L_m/t_m}{L_p/t_p}$$ 动力相似： 对应点的力的方向相同，大小成比例。这要求所有相关无量纲数相等：

• 雷诺数相等：$(Re)_m = (Re)_p$
• 马赫数相等：$(Ma)_m = (Ma)_p$（可压缩流）
• 弗劳德数相等：$(Fr)_m = (Fr)_p$（自由液面）

风洞实验的缩比原则

飞机模型测试（1:10缩比）：

要保持雷诺数相等： $$\frac{V_m L_m}{\nu_m} = \frac{V_p L_p}{\nu_p}$$ 如果 $L_m = L_p/10$，在相同流体中：

• 需要 $V_m = 10 V_p$
• 实际飞机 250 m/s → 模型需要 2500 m/s（不现实！）

实际解决方案：

1. 增压风洞：提高压力，降低运动粘度
2. 低温风洞：降低温度，降低粘度
3. 变雷诺数测试：测试多个雷诺数，外推结果
4. 部分模拟：只模拟关键区域的雷诺数

相似准则的限制

不可能同时满足所有相似准则：

船模实验的困境：

• 保持弗劳德数：$Fr = V/\sqrt{gL}$ → $V_m = V_p \sqrt{\lambda}$
• 保持雷诺数：$Re = VL/\nu$ → $V_m = V_p/\lambda$
• 两者矛盾！

解决方法：

• 阻力分解：摩擦阻力（依赖$Re$）+ 兴波阻力（依赖$Fr$）
• 分别修正

尺度效应：

某些现象在小尺度下消失：

• 转捩位置改变
• 分离特性不同
• 表面粗糙度的相对影响

相似分析的实例

汽车风洞测试（1:5模型）：

关注点：压差阻力为主

• 保证 $Re > Re_{crit}$（进入自模化区）
• 地面边界层模拟（移动地板）
• 轮胎旋转效应

测量修正： $$C_{D,full} = C_{D,model} \times (1 - 0.025 \log_{10}\frac{Re_{model}}{Re_{full}})$$ 建筑风洞（1:300模型）：

关注点：风压分布、涡激振动

• 大气边界层模拟（粗糙元、尖塔）
• 斯特劳哈尔数相似
• 风压系数直接适用（$C_p$无量纲）

时间尺度： $$\frac{f_m}{f_p} = \frac{St_m V_m/L_m}{St_p V_p/L_p} = \frac{V_m L_p}{V_p L_m}$$

历史人物：奥斯本·雷诺与1883年的染料实验

划时代的实验

1883年，英国曼彻斯特大学的奥斯本·雷诺（Osborne Reynolds, 1842-1912）进行了流体力学史上最著名的实验之一。他在一根水平玻璃管中注入染料，观察不同流速下的流动形态。

实验装置：

    染料注入
       ↓
[水箱]===╤═══════════════════════[阀门]
         │←——————玻璃管——————→│

低速：━━━━━━━━━━━━━━━━━━  层流
高速：～～～～～～～～～～～  湍流

关键发现：

1. 低速时，染料形成一条直线（层流）
2. 提高流速，染料线开始波动
3. 超过临界速度，染料迅速扩散（湍流）
4. 临界速度与管径、流体粘度有关

无量纲数的诞生

雷诺的天才在于他认识到，转捩不是由单一参数决定，而是由一个无量纲组合决定： $$\text{雷诺数} = \frac{\text{惯性力}}{\text{粘性力}} = \frac{\rho V D}{\mu}$$ 这个发现的革命性在于：

• 统一了不同尺度、不同流体的流动规律
• 使模型实验成为可能
• 奠定了相似理论的基础

雷诺的其他贡献

雷诺应力： 湍流中的附加应力项，源于速度脉动： $$\tau_{turb} = -\rho \overline{u'v'}$$ 雷诺分解： 将瞬时量分解为平均值和脉动值： $$u = \bar{u} + u'$$ 工程遗产：

• 管道设计的摩擦系数图（Moody图）
• 传热学中的雷诺类比
• 湍流模型的基础

实验的现代意义

雷诺实验看似简单，却包含深刻的物理：

非线性动力学： 层流-湍流转捩是典型的分岔现象，预示着混沌理论的诞生。

工程应用：

• 输油管道：保持层流减少泵功
• 血管流动：动脉粥样硬化与局部湍流相关
• 微流控芯片：利用低雷诺数实现精确控制

高级话题：π定理与量纲分析的深层应用

Buckingham π定理

设有$n$个物理量，涉及$k$个基本量纲，则可组成$(n-k)$个独立的无量纲数。

示例：圆柱绕流阻力

涉及的物理量：

• 阻力 $F$ [$MLT^{-2}$]
• 直径 $D$ [$L$]
• 流速 $V$ [$LT^{-1}$]
• 密度 $\rho$ [$ML^{-3}$]
• 粘度 $\mu$ [$ML^{-1}T^{-1}$]

5个物理量，3个基本量纲（M, L, T），可组成2个无量纲数：

• $\Pi_1 = \frac{F}{\rho V^2 D^2}$ （阻力系数）
• $\Pi_2 = \frac{\rho V D}{\mu}$ （雷诺数）

因此：$C_D = f(Re)$

不完全相似与自模化

自模化区域： 当雷诺数足够高时，某些流动特征不再随雷诺数变化：

圆柱阻力系数：

• $Re < 5$：$C_D \propto 1/Re$（Stokes流）
• $10^3 < Re < 10^5$：$C_D \approx 1.2$（亚临界）
• $Re > 10^6$：$C_D \approx 0.3$（超临界）

部分相似： 飞机起落架的缩比测试：

• 主体部分：确保 $Re > 10^5$（自模化）
• 细节部分：局部雷诺数匹配
• 缝隙流动：保持几何相似即可

量纲分析的巧妙应用

泰勒-库埃特不稳定性：

两同心圆柱间的旋转流动，内圆柱旋转角速度$\Omega$：

关键无量纲数（泰勒数）： $$Ta = \frac{\Omega^2 r_i d^3}{\nu^2}$$ 其中$d = r_o - r_i$是间隙宽度。

临界泰勒数：$Ta_{crit} \approx 1700$

应用：

• 轴承润滑
• 离心分离器
• 大气环流模拟

多重尺度问题

边界层的双重尺度：

外部尺度：$L$（物体特征长度） 内部尺度：$\delta \sim L/\sqrt{Re}$（边界层厚度）

匹配原则：

• 外部解在壁面奇异
• 内部解在远场发散
• 匹配区域：两解都有效

湍流的多重尺度：

能量注入尺度：$L$ 惯性子区：$\eta < \ell < L$ 耗散尺度（Kolmogorov尺度）： $$\eta = \left(\frac{\nu^3}{\varepsilon}\right)^{1/4}$$

尺度分离：$L/\eta \sim Re^{3/4}$

本章小结

核心概念

三大无量纲数：

1. 雷诺数 $Re = \rho VL/\mu$：惯性力与粘性力之比 - 决定层流/湍流状态 - 管流临界值 ≈ 2300 - 平板边界层临界值 ≈ 5×10⁵

2. 马赫数 $Ma = V/a$：流速与声速之比 - $Ma < 0.3$：不可压缩 - $Ma ≈ 1$：跨声速，激波形成 - $Ma > 1$：超声速流动

3. 斯特劳哈尔数 $St = fL/V$：非定常特征 - 圆柱涡脱落：$St ≈ 0.2$ - 用于涡激振动分析

相似准则的应用

模型实验三要素：

• 几何相似：形状缩放
• 运动相似：流线图案相同
• 动力相似：无量纲数相等

实际限制：

• 不能同时满足所有相似准则
• 利用自模化区域简化
• 分项修正补偿

关键公式速查

| 无量纲数 | 定义 | 物理意义 | 典型应用 |

经验法则

1. 雷诺数量级估算： - 日常物体在空气中：$Re \sim 10^4-10^7$ - 在水中：$Re$约为空气中的15倍

2. 临界雷诺数记忆： - 圆管：2300 - 平板：500,000 - 圆球：200,000

3. 马赫数效应： - $Ma < 0.3$：忽略压缩性 - $Ma > 0.7$：必须考虑压缩性 - $Ma ≈ 0.85$：典型客机巡航

练习题

基础题

习题2.1：游泳池跳水 一位跳水运动员（身高1.7m）以2 m/s的速度入水。估算雷诺数并判断流动类型。水的运动粘度 $\nu = 10^{-6}$ m²/s。

提示：使用身高作为特征长度。

习题2.2：高铁隧道压力波 高铁以350 km/h通过隧道，估算马赫数。讨论为什么会产生"隧道爆破"现象。

提示：将速度转换为m/s，标准声速340 m/s。

习题2.3：风铃频率 直径5mm的风铃管在10 m/s的风中发出声音。估算声音频率。

提示：使用斯特劳哈尔数 $St = 0.2$。

挑战题

习题2.4：缩比模型的困境 要在风洞中测试1:10的飞机模型，实机巡航速度250 m/s，雷诺数$10^8$。风洞最大风速100 m/s。提出三种可能的解决方案并讨论优缺点。

提示：考虑改变流体性质或测试条件。

习题2.5：血管中的流动 人体主动脉直径25mm，血流速度0.3 m/s，血液运动粘度$4 \times 10^{-6}$ m²/s。毛细血管直径8μm，流速0.5 mm/s。分别计算雷诺数并讨论生理意义。

提示：考虑不同雷诺数下的混合效率。

习题2.6：昆虫飞行之谜 果蝇翅膀弦长1mm，拍动频率200 Hz，翅尖速度1 m/s。计算雷诺数和斯特劳哈尔数，解释为什么昆虫飞行不能用常规翼型理论解释。

提示：空气运动粘度$1.5 \times 10^{-5}$ m²/s。

习题2.7：开放性思考题 设计一个实验验证斯特劳哈尔数的普适性。要求：材料易得，现象明显，可定量测量。

常见陷阱与错误

1. 特征长度的选择

错误：随意选择特征长度 正确：根据物理过程选择

示例：

• 管流：直径D
• 平板：距前缘距离x
• 翼型：弦长c
• 钝体：迎风面宽度

2. 雷诺数的误用

错误：认为高雷诺数一定是湍流 正确：需要考虑扰动水平和入口条件

• 精心设计的风洞可在$Re = 10^7$保持层流
• 粗糙表面可在$Re = 1000$触发湍流

3. 相似准则的局限

错误：认为无量纲数相等就完全相似 正确：还需考虑：

• 边界条件相似
• 初始条件相似
• 非定常过程的时间尺度
• 表面粗糙度的缩放

4. 马赫数效应的忽视

错误：在$Ma = 0.2$时完全忽略压缩性 正确：

密度变化：$\Delta\rho/\rho ≈ Ma^2/2$

• $Ma = 0.2$：2%变化
• $Ma = 0.3$：4.5%变化
• 精确测量时需要修正

5. 斯特劳哈尔数的适用范围

错误：所有涡脱落都是$St = 0.2$ 正确：

• 仅对圆柱在特定$Re$范围
• 形状改变，$St$改变
• 受限空间中$St$会变化

最佳实践检查清单

问题定义阶段

• [ ] 明确主导物理过程（惯性、粘性、压缩性、重力...）
• [ ] 识别所有相关无量纲数
• [ ] 估算各无量纲数的量级
• [ ] 判断哪些效应可以忽略

实验设计阶段

• [ ] 选择正确的特征长度和速度
• [ ] 计算模型的无量纲数
• [ ] 检查是否在自模化区
• [ ] 确定需要的修正方法
• [ ] 验证测量精度要求

数据分析阶段

• [ ] 用无量纲形式整理数据
• [ ] 检查无量纲数的依赖关系
• [ ] 识别临界值和转换区
• [ ] 与已知关联式对比
• [ ] 评估尺度效应的影响

工程应用阶段

• [ ] 从模型到原型的正确缩放
• [ ] 考虑部分相似的补偿
• [ ] 留出安全裕度
• [ ] 记录假设和限制
• [ ] 建议验证方法