第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)

在亚声速流动中，扰动能向各个方向传播，整个流场相互影响。特征包括：

  1. 椭圆型方程特性：下游条件能影响上游流动
  2. 压力恢复良好：流动减速时压力能有效恢复
  3. 无激波：压力变化连续
  4. 可用不可压假设：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
                    \___________/

特征现象：

  1. 局部超声速区：即使来流亚声速，翼型上表面可能超声速
  2. 激波形成：超声速区通过激波回到亚声速
  3. 激波-边界层干扰：可能导致分离
  4. 阻力急剧增加：激波阻力出现

10.3.3 超声速流动 (Ma > 1.2)

超声速流动具有独特的特征：

  1. 双曲型方程特性：上游不受下游影响
  2. 激波和膨胀波主导：压力通过间断面变化
  3. 波系干扰复杂：激波相交、反射形成复杂流场
  4. 气动加热显著：动能转化为热能

典型应用：

  • 战斗机（最大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 的点

阻力增加的原因：

  1. 激波阻力：激波造成总压损失
  2. 激波诱导分离：激波后逆压梯度导致边界层分离
  3. 形状阻力增加：分离区扩大

10.4.3 延迟阻力发散的方法

工程上采用多种方法推迟阻力发散：

  1. 超临界翼型： - 上表面平坦，延缓加速 - 下表面后部反弯，提供升力 - 可将Ma_DD提高0.1-0.15

  2. 后掠翼： - 有效马赫数 $Ma_{eff} = Ma_{\infty} \cos(\Lambda)$ - 30°后掠可将Ma_cr提高约15%

  3. 翼型厚度控制： - 减小厚弦比 - 最大厚度位置后移

  4. 层流控制： - 保持层流边界层 - 减小速度峰值

10.5 面积法则

10.5.1 惠特科姆面积法则

理查德·惠特科姆在1952年发现的面积法则是跨声速空气动力学的重大突破：

核心思想：跨声速飞行器的波阻主要取决于横截面积沿轴向的分布，而不是具体形状。

理想的面积分布应该平滑变化，避免突变。这就像声音通过变截面管道：截面积突变会产生反射波。

10.5.2 面积法则的应用

实际应用中的"蜂腰"设计：

    普通机身：
    ========[===]========  机翼处截面积突增

    面积法则修形后：
    ======)(===)(======  机身收缩补偿机翼面积
          蜂腰

经典案例：

  1. F-102战斗机：应用面积法则后，最大速度从Ma 0.98提升到Ma 1.22
  2. B-58轰炸机：发动机吊舱位置精心设计以满足面积法则
  3. 波音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）首次提出用马赫的名字命名这个无量纲数。这个选择非常恰当，因为：

  1. 马赫首次直观展示了超声速现象
  2. 他的工作奠定了高速空气动力学基础
  3. 他强调了声速作为特征速度的重要性

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点后: 另一解支

这种现象的物理根源：

  1. 混合型方程：信息传播方式在声速点改变
  2. 激波位置不定：小扰动可能导致激波大幅移动
  3. 分离敏感性：激波诱导分离的临界特性

10.7.4 现代跨声速设计工具

现代CFD方法处理跨声速流动的策略：

  1. 自适应网格： - 激波附近加密 - 基于压力梯度的自适应准则

  2. 人工粘性： - 捕捉激波 - 避免数值振荡

  3. 多重网格方法： - 加速收敛 - 处理不同尺度

  4. 伴随方法优化： - 敏感性分析 - 形状优化

10.7.5 跨声速风洞的特殊考虑

跨声速风洞设计必须处理独特挑战：

  1. 壁面干扰： - 开孔壁/槽壁减少堵塞 - 主动壁面控制

  2. 雷诺数效应： - 变压风洞提高雷诺数 - 低温风洞（如ETW）

  3. 动态测试： - 跨声速颤振 - 非定常激波运动

10.8 本章小结

本章探讨了可压缩性效应如何从根本上改变流动特性。关键概念包括：

  1. 声速的物理本质： - 小扰动传播速度：$a = \sqrt{\gamma RT}$ - 只与温度相关，与压力无关 - 反映分子热运动速度

  2. 马赫数的核心地位： - 定义流动regime的关键参数 - 决定信息传播模式 - 控制方程类型的转变

  3. 跨声速流动的复杂性： - 局部超声速区与激波 - 混合型方程特征 - 解的非唯一性

  4. 阻力发散与缓解策略： - 临界马赫数概念 - 超临界翼型设计 - 后掠翼的作用

  5. 面积法则的工程价值： - 横截面积分布决定波阻 - 蜂腰设计的物理基础 - 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

  1. 细节优化： - 过渡长度约2倍机翼弦长 - 避免二阶导数突变 - 整流罩采用双曲线型面 - 考虑结构和系统布置约束

  2. 预期效果： - 波阻降低20-30% - Ma_DD提高0.02-0.03 - 燃油经济性改善3-5%

10.10 常见陷阱与错误

10.10.1 概念理解错误

  1. "声速是常数" - 错误：认为声速在所有条件下都是340 m/s - 正确：声速随温度变化，高空声速明显降低

  2. "马赫数只与速度有关" - 错误：忽略了当地声速的变化 - 正确：同样速度在不同高度对应不同马赫数

  3. "超声速就一定有激波" - 错误：认为Ma > 1就必然产生激波 - 正确：等熵加速到超声速可以无激波（拉瓦尔喷管）

10.10.2 工程应用陷阱

  1. 忽略局部马赫数 - 陷阱：只考虑来流马赫数 - 正确：翼型上局部Ma可能比来流高50%

  2. 线性理论的误用 - 陷阱：在跨声速区使用普朗特-格劳特公式 - 正确：Ma > 0.7后需要考虑非线性效应

  3. 面积法则的片面理解 - 陷阱：只在机身-机翼结合处应用 - 正确：考虑所有部件（发动机、尾翼等）的综合效应

10.10.3 设计决策失误

  1. 过度追求层流 - 问题：层流翼型往往临界马赫数较低 - 平衡：需要在层流benefits和跨声速性能间权衡

  2. 忽视跨声速配平 - 问题：激波移动导致压心大幅变化 - 解决：需要考虑配平阻力和操纵性

  3. 雷诺数效应 - 陷阱：风洞试验的跨声速特性可能与实际不符 - 原因：转捩位置影响激波-边界层干扰

10.11 最佳实践检查清单

设计审查要点

初步设计阶段

  • [ ] 确定巡航马赫数和高度包线
  • [ ] 计算各工况下的雷诺数范围
  • [ ] 识别临界部件（最先遇到跨声速效应）
  • [ ] 评估是否需要超临界翼型
  • [ ] 确定后掠角需求

详细设计阶段

  • [ ] 验证临界马赫数满足要求
  • [ ] 检查横截面积分布符合面积法则
  • [ ] 优化机翼-机身、发动机-机翼结合部
  • [ ] 考虑激波位置的移动范围
  • [ ] 评估激波-边界层干扰风险

性能评估

  • [ ] 计算阻力发散特性
  • [ ] 确认加速通过跨声速区的推力裕度
  • [ ] 检查跨声速配平特性
  • [ ] 评估抖振边界
  • [ ] 验证操纵面效率

试验验证

  • [ ] 规划跨声速风洞试验
  • [ ] 考虑壁面干扰修正
  • [ ] 安排压力分布测量
  • [ ] 准备纹影或PSP流场显示
  • [ ] 制定飞行试验逐步推进计划

运行考虑

  • [ ] 确定经济巡航马赫数
  • [ ] 制定超速保护逻辑
  • [ ] 考虑高空/低温对马赫数限制的影响
  • [ ] 评估阵风响应特性
  • [ ] 准备机组培训材料（跨声速特性）