intuitive_cfd

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

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

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

雷诺数的定义与物理意义

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

\[Re = \frac{\rho V L}{\mu} = \frac{VL}{\nu}\]

其中：

物理直觉：

日常生活中的雷诺数

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

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

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

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

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

临界雷诺数与转捩

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

圆管内流动：

平板边界层：

圆柱绕流：

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

雷诺数的工程应用

管道设计：

换热器设计：

减阻设计：

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

声速与马赫数定义

马赫数（Mach number, $Ma$）定义为流速与当地声速之比：

\[Ma = \frac{V}{a}\]

其中声速 $a = \sqrt{\gamma RT}$（理想气体）

标准大气条件下的声速：

流动分类

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

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

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

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

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

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

生活中的马赫数效应

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

高铁过隧道：

子弹飞行：

压缩性效应的估算

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

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

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

涡街现象与斯特劳哈尔数

斯特劳哈尔数（Strouhal number, $St$）描述了非定常流动中的特征频率：

\[St = \frac{fL}{V}\]

其中：

对于圆柱绕流，在很宽的雷诺数范围内（$200 < Re < 10^5$），斯特劳哈尔数近似为常数： \(St \approx 0.2\)

这意味着涡脱落频率与流速成正比： \(f = 0.2 \frac{V}{D}\)

日常生活中的涡街

电线在风中的歌唱：

烟囱涡激振动：

潜望镜尾流：

斯特劳哈尔数的变化规律

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

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

工程应用

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

桥梁设计：

高层建筑：

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

三种相似性

几何相似： 模型与原型的所有对应线性尺寸成比例： \(\frac{L_m}{L_p} = \lambda\)

运动相似： 对应点的速度方向相同，大小成比例： \(\frac{V_m}{V_p} = \frac{L_m/t_m}{L_p/t_p}\)

动力相似： 对应点的力的方向相同，大小成比例。这要求所有相关无量纲数相等：

风洞实验的缩比原则

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

要保持雷诺数相等： \(\frac{V_m L_m}{\nu_m} = \frac{V_p L_p}{\nu_p}\)

如果 $L_m = L_p/10$，在相同流体中：

实际解决方案：

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

相似准则的限制

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

船模实验的困境：

解决方法：

尺度效应：

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

相似分析的实例

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

关注点：压差阻力为主

测量修正： \(C_{D,full} = C_{D,model} \times (1 - 0.025 \log_{10}\frac{Re_{model}}{Re_{full}})\)

建筑风洞（1:300模型）：

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

时间尺度： \(\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'\)

工程遗产：

实验的现代意义

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

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

工程应用：

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

Buckingham π定理

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

示例：圆柱绕流阻力

涉及的物理量：

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

因此：$C_D = f(Re)$

不完全相似与自模化

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

圆柱阻力系数：

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

量纲分析的巧妙应用

泰勒-库埃特不稳定性：

两同心圆柱间的旋转流动，内圆柱旋转角速度$\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$
    • 用于涡激振动分析

相似准则的应用

模型实验三要素：

实际限制：

关键公式速查

无量纲数 定义 物理意义 典型应用
雷诺数 $\frac{\rho VL}{\mu}$ 惯性/粘性 流态判断
马赫数 $\frac{V}{a}$ 流速/声速 压缩性判断
斯特劳哈尔数 $\frac{fL}{V}$ 频率特征 涡激振动
弗劳德数 $\frac{V}{\sqrt{gL}}$ 惯性/重力 船舶阻力
韦伯数 $\frac{\rho V^2 L}{\sigma}$ 惯性/表面张力 液滴破碎

经验法则

  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。

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

答案 $Re = \frac{VL}{\nu} = \frac{2 \times 1.7}{10^{-6}} = 3.4 \times 10^6$ 这是高雷诺数湍流。运动员周围会形成湍流尾流，这就是为什么跳水时会产生大量气泡和水花。

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

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

答案 $V = 350 \text{ km/h} = 97.2 \text{ m/s}$ $Ma = 97.2/340 = 0.286$ 虽然马赫数较低，但在受限空间（隧道）中，列车推动的空气无法侧向逸出，形成压力波。压力波以声速传播，在隧道出口产生微压波，即"隧道爆破"。 解决方案：隧道入口设置缓冲结构，列车头部优化为细长形。

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

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

答案 $f = St \times \frac{V}{D} = 0.2 \times \frac{10}{0.005} = 400$ Hz 这个频率在人耳敏感范围内（中音区），所以风铃声音悦耳。不同直径的管产生不同频率，形成和声。

挑战题

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

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

答案 实机：$Re_p = 10^8$ 模型要求相同雷诺数：$\frac{V_m L_m}{\nu_m} = \frac{V_p L_p}{\nu_p}$ 由于$L_m = 0.1L_p$，需要： - 方案1：提高速度到2500 m/s（不现实） - 方案2：降低粘度到1/25（使用加压或低温） - 方案3：接受雷诺数不匹配，在自模化区测试 **实际方案**： 1. **加压风洞**：5倍大气压，$\nu$降至1/5，可达$Re = 2 \times 10^7$ 2. **低温风洞**：-180°C液氮冷却，$\nu$降至1/10 3. **变雷诺数外推**：测试多个雷诺数，找出趋势 优选方案3，成本最低，结合理论修正。

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

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

答案 主动脉： $Re_{aorta} = \frac{0.3 \times 0.025}{4 \times 10^{-6}} = 1875$ 接近层流-湍流转换，脉动时可能局部湍流，有助于混合但不会损伤血管壁。 毛细血管： $Re_{cap} = \frac{0.0005 \times 8 \times 10^{-6}}{4 \times 10^{-6}} = 0.001$ 极低雷诺数，完全层流，扩散主导。这保证了氧气和营养物质的充分交换时间。 **生理意义**： - 大血管的适度湍流防止血栓 - 毛细血管的层流保证物质交换 - 病理狭窄处雷诺数增加，易形成涡流和血栓

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

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

答案 雷诺数： $Re = \frac{1 \times 0.001}{1.5 \times 10^{-5}} = 67$ 斯特劳哈尔数： $St = \frac{200 \times 0.001}{1} = 0.2$ **分析**： 1. 低雷诺数（$Re < 100$）：粘性效应显著，常规翼型理论（基于高Re）失效 2. 高斯特劳哈尔数：强非定常效应，每个拍动周期都在起动和停止 3. 前缘涡（LEV）机制：低Re下前缘涡不脱落，提供额外升力 4. 拍合机制（clap-and-fling）：两翅拍合后快速分开，产生强环量 这就是为什么昆虫能产生远超稳态理论预测的升力系数（可达4-6）。

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

参考方案 **实验设计：弹性带振动** 材料： - 弹性带（不同宽度） - 风扇（可调速） - 手机慢动作摄像 - 尺子 步骤： 1. 固定弹性带一端，另一端自由 2. 用风扇吹，调节风速直到稳定振动 3. 慢动作拍摄，数振动频率 4. 测量风速（用轻纸片飘移距离估算） 5. 改变带宽，重复实验 预期结果： - $St = fD/V ≈ 0.15-0.2$（平板） - 频率与风速成正比 - 频率与带宽成反比 **延伸**： - 比较不同截面形状（圆棒、方棒） - 研究张力的影响 - 探索共振条件

常见陷阱与错误

1. 特征长度的选择

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

示例：

2. 雷诺数的误用

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

3. 相似准则的局限

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

4. 马赫数效应的忽视

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

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

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

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

最佳实践检查清单

问题定义阶段

实验设计阶段

数据分析阶段

工程应用阶段