第10章：声速与马赫数效应
当一架战斗机以接近声速飞行时，机翼上会突然出现一团白雾；当鞭子甩动时末端会发出清脆的爆响；当高铁进入隧道时会感受到明显的压力波动——这些现象都与流体的可压缩性密切相关。本章将探讨声速的物理本质，以及当流动速度接近或超过声速时出现的各种有趣现象。我们将学习如何通过马赫数来判断流动特性，理解为什么音障曾经被认为是不可逾越的，以及现代飞行器如何巧妙地处理跨声速流动问题。
10.1 声速的物理意义
10.1.1 声速的本质：压力扰动的传播
声速并不是声音特有的速度，而是小扰动在介质中的传播速度。想象在静止的池塘中投入一颗石子，涟漪以固定速度向外扩散——这就是压力波的传播。在空气中，这个速度就是声速。
从微观角度看，声速反映了分子热运动的平均速度。温度越高，分子运动越剧烈，信息传递越快，声速也就越大。这就是为什么：
\[a = \sqrt{\gamma RT}\]
其中 $a$ 是声速，$\gamma$ 是比热比（空气约为1.4），$R$ 是气体常数，$T$ 是绝对温度。
在标准海平面条件下（15°C），空气中的声速约为 340 m/s。有趣的是，声速只与温度有关，与压力无关——这是因为压力增加时，密度也成比例增加，两者的效应相互抵消。
10.1.2 声速的测量与估算
工程上常用的声速估算公式：
\[a \approx 20\sqrt{T} \text{ (m/s，T单位为K)}\]
这个简化公式误差在1%以内，非常实用。例如：
- 0°C (273K): $a \approx 20\sqrt{273} = 331$ m/s
- 20°C (293K): $a \approx 20\sqrt{293} = 343$ m/s
- 高空-50°C (223K): $a \approx 20\sqrt{223} = 299$ m/s
10.1.3 不同介质中的声速
声速在不同介质中差异很大：
- 水中：约 1480 m/s（比空气快4倍多）
- 钢铁中：约 5000 m/s
- 氦气中：约 970 m/s（所以吸入氦气后声音变尖）
这种差异源于介质的可压缩性和密度。液体和固体的可压缩性远小于气体，所以声速更快。
10.2 马赫数与流动regime
10.2.1 马赫数的定义与物理意义
马赫数定义为流动速度与当地声速之比：
\[Ma = \frac{V}{a}\]
马赫数告诉我们流动相对于信息传播的快慢。这就像在河流中的船只：
- Ma < 1：船速小于水波速度，涟漪能向前传播
- Ma = 1：船速等于水波速度，波浪在船头堆积
- Ma > 1：船速超过水波速度，形成V形尾迹
10.2.2 马赫锥与信息传播
当物体以超声速运动时，它产生的扰动被限制在一个锥形区域内，称为马赫锥：
物体运动方向 →
·
/|\
/ | \
/ | \ 马赫锥
/ | \
/ α | \
/ | \
/______|______\
sin(α) = 1/Ma
马赫角 α = arcsin(1/Ma)，马赫数越大，锥角越小。这解释了为什么：
- 超声速飞机飞过头顶时，先看到飞机，后听到声音
- 子弹的弹道冲击波呈锥形
- 超声速飞行器的激波附体在尖锐前缘
10.3 亚声速、跨声速、超声速流动特征
10.3.1 亚声速流动 (Ma < 0.8)
在亚声速流动中，扰动能向各个方向传播，整个流场相互影响。特征包括：
- 椭圆型方程特性：下游条件能影响上游流动
- 压力恢复良好：流动减速时压力能有效恢复
- 无激波：压力变化连续
- 可用不可压假设：Ma < 0.3时，密度变化小于5%
日常例子：
- 汽车行驶（通常Ma < 0.15）
- 民航客机巡航（Ma ≈ 0.85）
- 风力发电机叶片（叶尖Ma < 0.5）
10.3.2 跨声速流动 (0.8 < Ma < 1.2)
跨声速区是最复杂的流动区域，局部同时存在亚声速和超声速区域：
翼型上的跨声速流动：
来流 Ma=0.85 → ___________
/ \___
| 超声速区 激波
| Ma>1.0 ↓
============================================
| 亚声速区
| Ma<1.0
\___________/
特征现象：
- 局部超声速区：即使来流亚声速，翼型上表面可能超声速
- 激波形成：超声速区通过激波回到亚声速
- 激波-边界层干扰：可能导致分离
- 阻力急剧增加：激波阻力出现
10.3.3 超声速流动 (Ma > 1.2)
超声速流动具有独特的特征：
- 双曲型方程特性：上游不受下游影响
- 激波和膨胀波主导：压力通过间断面变化
- 波系干扰复杂：激波相交、反射形成复杂流场
- 气动加热显著：动能转化为热能
典型应用：
- 战斗机（最大Ma ≈ 2.5）
- 超声速客机（协和式Ma = 2.0）
- 弹道导弹（再入Ma > 20）
10.4 临界马赫数与阻力发散
10.4.1 临界马赫数的概念
临界马赫数 $Ma_{cr}$ 是翼型上第一次出现声速点（Ma = 1.0）时的来流马赫数。对于典型的亚声速翼型：
- 对称翼型：$Ma_{cr} \approx 0.7-0.75$
- 传统翼型（有弯度）：$Ma_{cr} \approx 0.65-0.70$
- 超临界翼型：$Ma_{cr} \approx 0.75-0.80$
临界马赫数可以通过普朗特-格劳特法则估算：
\[C_p = \frac{C_{p,incomp}}{\sqrt{1-Ma^2}}\]
当最大负压系数处首次达到声速时，即为临界条件。
10.4.2 阻力发散现象
当飞行马赫数超过临界马赫数后，阻力系数急剧增加：
CD (阻力系数)
↑
| /
| / 阻力发散
| /
| /
|__/_______________
Ma_cr Ma_DD Ma →
Ma_DD: 阻力发散马赫数
定义为 dCD/dMa = 0.1 的点
阻力增加的原因：
- 激波阻力：激波造成总压损失
- 激波诱导分离：激波后逆压梯度导致边界层分离
- 形状阻力增加：分离区扩大
10.4.3 延迟阻力发散的方法
工程上采用多种方法推迟阻力发散：
- 超临界翼型：
- 上表面平坦，延缓加速
- 下表面后部反弯，提供升力
- 可将Ma_DD提高0.1-0.15
- 后掠翼：
- 有效马赫数 $Ma_{eff} = Ma_{\infty} \cos(\Lambda)$
- 30°后掠可将Ma_cr提高约15%
- 翼型厚度控制：
- 层流控制：
10.5 面积法则
10.5.1 惠特科姆面积法则
理查德·惠特科姆在1952年发现的面积法则是跨声速空气动力学的重大突破：
核心思想：跨声速飞行器的波阻主要取决于横截面积沿轴向的分布，而不是具体形状。
理想的面积分布应该平滑变化，避免突变。这就像声音通过变截面管道：截面积突变会产生反射波。
10.5.2 面积法则的应用
实际应用中的”蜂腰”设计：
普通机身：
========[===]======== 机翼处截面积突增
面积法则修形后：
======)(===)(====== 机身收缩补偿机翼面积
蜂腰
经典案例：
- F-102战斗机：应用面积法则后，最大速度从Ma 0.98提升到Ma 1.22
- B-58轰炸机：发动机吊舱位置精心设计以满足面积法则
- 波音747：机身-机翼连接处的整流设计
10.5.3 面积法则的定量应用
跨声速波阻可以近似为：
\[C_{D,wave} \propto \int_0^L \left(\frac{d^2A}{dx^2}\right)^2 dx\]
其中A(x)是横截面积分布。这表明：
- 面积二阶导数越小越好
- 光滑过渡优于突变
- 总体积相同时，细长体阻力更小
实践中的设计准则：
- 面积分布应接近Sears-Haack体
- 最大截面积位置在50-60%机身长度处
- 面积增长率控制在特定范围内
10.6 历史人物：恩斯特·马赫与超声速摄影
10.6.1 从哲学家到物理学家
恩斯特·马赫（Ernst Mach, 1838-1916）是一位独特的科学家——他既是物理学家、哲学家，也是心理学家。他最初研究多普勒效应，试图通过实验验证这一理论，这把他引向了声学和激波的研究。
马赫的哲学思想深刻影响了他的科学研究。他坚持”思维经济原理”——科学理论应该以最简洁的方式描述现象。这种思想later影响了爱因斯坦的相对论发展。
10.6.2 开创性的激波摄影
1887年，马赫和他的学生彼得·萨尔切尔（Peter Salcher）首次拍摄到子弹飞行产生的激波。他们使用的纹影技术（Schlieren photography）至今仍是流动可视化的重要方法：
光源 → 凹面镜 → 测试区 → 凹面镜 → 刀口 → 相机
↑
密度变化
折射光线
马赫的照片清晰显示了：
这些照片不仅证实了理论预测，更重要的是提供了直观的物理图像。
10.6.3 马赫数的命名
有趣的是，马赫本人从未使用过”马赫数”这个概念。1929年，瑞士工程师雅各布·阿克莱特（Jakob Ackeret）首次提出用马赫的名字命名这个无量纲数。这个选择非常恰当，因为：
- 马赫首次直观展示了超声速现象
- 他的工作奠定了高速空气动力学基础
- 他强调了声速作为特征速度的重要性
10.6.4 马赫的科学遗产
马赫的贡献远超空气动力学：
- 马赫带：视觉感知中的边缘增强现象
- 马赫原理：惯性源于宇宙中所有物质的相互作用
- 感觉分析：开创了实验心理学
- 科学哲学：影响了维也纳学派和逻辑实证主义
他的名言”描述事实，仅此而已”体现了他的科学态度——这正是我们研究流体现象时应该秉持的精神。
10.7 高级话题：跨声速相似律与临界点理论
10.7.1 卡门-钱学森公式
跨声速相似律是理解高亚声速流动的关键。1939年，冯·卡门提出，1946年钱学森完善了这一理论：
对于细长体，压力系数修正为：
\[C_p = \frac{C_{p0}}{\sqrt{1-Ma_{\infty}^2} + \frac{Ma_{\infty}^2}{1+\frac{\gamma-1}{2}Ma_{\infty}^2}\frac{C_{p0}}{2}}\]
其中 $C_{p0}$ 是不可压流中的压力系数。
这个公式的物理意义：
- 第一项：普朗特-格劳特修正（线性理论）
- 第二项：非线性效应（局部超声速修正）
10.7.2 跨声速小扰动方程
跨声速流动的控制方程具有混合型特征：
\[\left(1-Ma_{\infty}^2 - \frac{\gamma+1}{U_{\infty}}Ma_{\infty}^2\phi_x\right)\phi_{xx} + \phi_{yy} + \phi_{zz} = 0\]
这个方程在不同区域表现不同：
- 亚声速区：椭圆型（所有二阶导数同号）
- 超声速区：双曲型（存在特征线）
- 声速线上：抛物型（过渡区）
10.7.3 临界点理论与跨声速解的非唯一性
跨声速流动最令人困惑的特征之一是解的非唯一性。给定相同的边界条件，可能存在多个稳定解：
升力系数 CL
↑
| ___B
| / \
| / \___C
| /
| /A
|/_____________
Ma →
A-B: 连续解
B-C: 不连续跳跃
C点后: 另一解支
这种现象的物理根源：
- 混合型方程：信息传播方式在声速点改变
- 激波位置不定：小扰动可能导致激波大幅移动
- 分离敏感性：激波诱导分离的临界特性
10.7.4 现代跨声速设计工具
现代CFD方法处理跨声速流动的策略：
- 自适应网格：
- 人工粘性：
- 多重网格方法：
- 伴随方法优化：
10.7.5 跨声速风洞的特殊考虑
跨声速风洞设计必须处理独特挑战：
- 壁面干扰：
- 雷诺数效应：
- 动态测试：
10.8 本章小结
本章探讨了可压缩性效应如何从根本上改变流动特性。关键概念包括：
- 声速的物理本质：
- 小扰动传播速度：$a = \sqrt{\gamma RT}$
- 只与温度相关，与压力无关
- 反映分子热运动速度
- 马赫数的核心地位：
- 定义流动regime的关键参数
- 决定信息传播模式
- 控制方程类型的转变
- 跨声速流动的复杂性：
- 阻力发散与缓解策略：
- 面积法则的工程价值：
- 横截面积分布决定波阻
- 蜂腰设计的物理基础
- Sears-Haack理想体
关键公式汇总：
- 声速：$a = \sqrt{\gamma RT} \approx 20\sqrt{T}$ m/s
- 马赫数：$Ma = V/a$
- 马赫角：$\sin\alpha = 1/Ma$
- 普朗特-格劳特修正：$C_p = C_{p,incomp}/\sqrt{1-Ma^2}$
- 有效马赫数（后掠翼）：$Ma_{eff} = Ma_{\infty}\cos\Lambda$
10.9 练习题
题目1（基础）
一架民航客机在11000米高空巡航，该高度温度约为-56°C。如果巡航速度为850 km/h，计算飞行马赫数。该飞机是否会遇到跨声速效应？
提示：先计算该温度下的声速，注意单位转换。
答案
1. 温度：T = -56 + 273 = 217 K
2. 声速：a = 20√217 = 20 × 14.73 = 295 m/s
3. 飞行速度：V = 850 km/h = 850/3.6 = 236 m/s
4. 马赫数：Ma = 236/295 = 0.80
结论：Ma = 0.80，刚好在跨声速区边缘，翼型上可能出现局部超声速区。
题目2（基础）
战斗机以Ma = 1.5飞行，计算马赫锥的半角。如果飞机在你头顶10公里处飞过，声音传到地面时飞机已经飞了多远？
提示：使用马赫角公式，考虑几何关系。
答案
1. 马赫角：sin α = 1/1.5 = 0.667，α = 41.8°
2. 马赫锥半角即为41.8°
3. 几何关系：飞机飞行距离 = 10 km × tan(90° - 41.8°) = 10 × tan(48.2°) = 11.2 km
当声音传到地面时，飞机已经向前飞行了11.2公里。
题目3（基础）
某超临界翼型的临界马赫数为0.75，普通翼型为0.65。如果要在Ma = 0.85巡航，采用30°后掠角是否足够避免激波产生？
提示：计算有效马赫数，与临界马赫数比较。
答案
1. 有效马赫数：Ma_eff = 0.85 × cos(30°) = 0.85 × 0.866 = 0.736
2. 对于超临界翼型：Ma_eff = 0.736 < Ma_cr = 0.75，理论上可避免强激波
3. 对于普通翼型：Ma_eff = 0.736 > Ma_cr = 0.65，仍会产生激波
结论：30°后掠角配合超临界翼型可以缓解激波问题，但对普通翼型不够。
题目4（挑战）
F-104战斗机被称为”有人驾驶的导弹”，其机翼极薄（厚弦比仅3.4%）且几乎无后掠。请分析这种设计的空气动力学考虑，以及在不同马赫数下的优缺点。
提示：考虑激波强度与翼型厚度的关系，以及不同速度regime的需求。
答案
设计分析：
1. **超声速优化**（Ma > 1.5）：
- 薄翼减小激波角度和强度
- 直机翼在超声速下阻力反而较小
- 激波阻力 ∝ (t/c)²
2. **跨声速问题**（0.8 < Ma < 1.2）：
- 无后掠导致临界马赫数很低
- 阻力发散严重
- 加速通过跨声速区需要强大推力
3. **亚声速性能差**：
- 薄翼升力线斜率小
- 失速特性恶劣
- 着陆速度高（>300 km/h）
4. **结构考虑**：
- 薄翼结构效率低
- 需要特殊材料和制造工艺
这是典型的单点设计——为高超声速性能牺牲其他性能。
题目5（挑战）
高速列车进入隧道时会产生压缩波。如果列车速度300 km/h，隧道横截面积100 m²，列车横截面积12 m²，估算初始压缩波的压力升高。该现象与航空中的什么问题类似？
提示：考虑活塞效应和面积堵塞比。
答案
分析：
1. **基本参数**：
- V = 300 km/h = 83.3 m/s
- Ma = 83.3/340 = 0.245
- 堵塞比 = 12/100 = 0.12
2. **压力升高估算**：
使用活塞理论：Δp/p ≈ γ·Ma·(A_train/A_tunnel)
Δp/p ≈ 1.4 × 0.245 × 0.12 = 0.041
压力升高约4.1%
3. **实际考虑**：
- 列车头部形状影响压力梯度
- 隧道入口设计可缓解压力波
- 多重反射可能放大效应
4. **航空类比**：
类似于风洞堵塞效应和跨声速风洞壁面干扰问题。都涉及有限空间内的压缩波传播。
题目6（开放思考）
“音障”曾被认为不可突破，但现在超声速飞行已经常规化。请讨论：
a) 为什么早期会有”音障”的错误认识？
b) 现代是否存在类似的”障碍”？（如”热障”）
c) 这对工程创新有什么启示？
答案
讨论要点：
a) **"音障"误解的原因**：
- 线性理论在Ma→1时发散（1/√(1-Ma²)）
- 早期飞机在跨声速区确实遇到严重问题
- 缺乏对激波物理的理解
- 风洞实验的局限性（壁面干扰）
b) **现代的"障碍"**：
- **热障**（Ma > 5）：材料耐温极限
- **控制障碍**：高超声速下的等离子体鞘
- **经济障碍**：超声速客运的成本问题
- **环境障碍**：声爆对地面的影响
c) **工程启示**：
- 理论局限性不等于物理不可能
- 突破需要新概念（如面积法则）
- 多学科综合（材料、控制、推进）
- 渐进式创新vs.颠覆式创新
关键认识：每个"障碍"都促进了新技术发展。
题目7（综合应用）
设计一个Ma = 0.9巡航的公务机机身-机翼结合部。已知机翼根部弦长4米，厚度0.4米，机身直径3米。如何应用面积法则优化设计？绘制横截面积分布草图。
提示：考虑面积补偿和过渡区设计。
答案
设计方案：
1. **问题分析**：
- 机翼增加横截面积：~1.6 m²（0.4×4）
- 需要在机翼区域减少机身面积
- Ma = 0.9处于跨声速区，面积法则关键
2. **具体措施**：
- 机翼前：机身直径从3.0米逐渐收缩到2.6米
- 机翼处：最小直径约2.4米（"蜂腰"）
- 机翼后：逐渐恢复到3.0米
3. **面积分布**（示意）：
```
面积(m²)
8 | ___
7 | / \___
6 | / \
5 |__/ \__
|__________________|
0 前 翼 后 x
```
4. **细节优化**：
- 过渡长度约2倍机翼弦长
- 避免二阶导数突变
- 整流罩采用双曲线型面
- 考虑结构和系统布置约束
5. **预期效果**：
- 波阻降低20-30%
- Ma_DD提高0.02-0.03
- 燃油经济性改善3-5%
10.10 常见陷阱与错误
10.10.1 概念理解错误
- “声速是常数”
- 错误：认为声速在所有条件下都是340 m/s
- 正确：声速随温度变化，高空声速明显降低
- “马赫数只与速度有关”
- 错误：忽略了当地声速的变化
- 正确：同样速度在不同高度对应不同马赫数
- “超声速就一定有激波”
- 错误：认为Ma > 1就必然产生激波
- 正确：等熵加速到超声速可以无激波（拉瓦尔喷管）
10.10.2 工程应用陷阱
- 忽略局部马赫数
- 陷阱：只考虑来流马赫数
- 正确：翼型上局部Ma可能比来流高50%
- 线性理论的误用
- 陷阱：在跨声速区使用普朗特-格劳特公式
- 正确：Ma > 0.7后需要考虑非线性效应
- 面积法则的片面理解
- 陷阱：只在机身-机翼结合处应用
- 正确：考虑所有部件（发动机、尾翼等）的综合效应
10.10.3 设计决策失误
- 过度追求层流
- 问题：层流翼型往往临界马赫数较低
- 平衡：需要在层流benefits和跨声速性能间权衡
- 忽视跨声速配平
- 问题：激波移动导致压心大幅变化
- 解决：需要考虑配平阻力和操纵性
- 雷诺数效应
- 陷阱：风洞试验的跨声速特性可能与实际不符
- 原因：转捩位置影响激波-边界层干扰
10.11 最佳实践检查清单
设计审查要点
初步设计阶段
详细设计阶段
性能评估
试验验证
运行考虑