当超声速战斗机从头顶呼啸而过，那震耳欲聋的音爆让人印象深刻——这就是激波现象最直观的体现。激波不仅存在于军用航空，在火箭发射、高速列车隧道通过、甚至开香槟时都能观察到。本章将揭示这些剧烈压缩现象背后的物理机制，以及与之对应的膨胀过程。通过学习激波与膨胀波，你将理解超声速流动中最基本也最重要的物理现象，掌握估算激波强度、预测流动偏转角度的方法，并能识别实际工程中的激波问题。
想象一列高速列车驶入隧道的瞬间。列车前方的空气来不及让开，被剧烈压缩，形成一道压力突跃——这就是激波的本质。在超声速流动中，扰动信息的传播速度（声速）小于流体速度，导致压缩波不断叠加，最终形成厚度仅几个分子自由程的间断面。
亚声速（Ma < 1）： 超声速（Ma > 1）：
扰动波可以向上游传播 扰动波被流动带向下游
← ← ← → → →
————————————→ ————————————→
流动方向 流动方向
压缩波叠加形成激波
从微观角度看，激波是分子剧烈碰撞的区域。在激波厚度（约3-4个分子平均自由程）内，分子速度分布严重偏离麦克斯韦分布，动能转化为内能，导致温度急剧上升。
激波前后必须满足三个守恒定律：
其中下标1表示激波前（上游），2表示激波后（下游）。这三个方程构成了朗肯-雨果尼奥条件，是分析所有激波问题的基础。
对于理想气体的正激波，可以导出以下关系式：
压力比： \(\frac{p_2}{p_1} = \frac{2\gamma Ma_1^2 - (\gamma-1)}{\gamma+1}\)
密度比： \(\frac{\rho_2}{\rho_1} = \frac{(\gamma+1)Ma_1^2}{(\gamma-1)Ma_1^2 + 2}\)
温度比： \(\frac{T_2}{T_1} = \left[1 + \frac{2(\gamma-1)}{(\gamma+1)^2}(Ma_1^2-1)\right]\left[\frac{2\gamma Ma_1^2}{\gamma+1} - \frac{\gamma-1}{\gamma+1}\right]\)
下游马赫数： \(Ma_2^2 = \frac{Ma_1^2 + \frac{2}{\gamma-1}}{2\frac{\gamma}{\gamma-1}Ma_1^2 - 1}\)
工程实践中，通常使用正激波表快速查询激波前后的参数关系。以下是空气（$\gamma=1.4$）的部分正激波表：
| $Ma_1$ | $p_2/p_1$ | $\rho_2/\rho_1$ | $T_2/T_1$ | $Ma_2$ |
|---|---|---|---|---|
| 1.5 | 2.46 | 1.86 | 1.32 | 0.70 |
| 2.0 | 4.50 | 2.67 | 1.69 | 0.58 |
| 2.5 | 7.13 | 3.33 | 2.14 | 0.51 |
| 3.0 | 10.3 | 3.86 | 2.68 | 0.48 |
经验法则：
激波过程伴随着熵增，是典型的不可逆过程：
\[s_2 - s_1 = c_v \ln\left[\frac{p_2/p_1}{(\rho_2/\rho_1)^\gamma}\right] > 0\]熵增导致总压损失： \(\frac{p_{02}}{p_{01}} = \left(\frac{p_2}{p_1}\right)^{\frac{1}{\gamma}}\left(\frac{\rho_1}{\rho_2}\right)\)
工程影响：
当超声速气流遇到楔形体时，气流必须偏转以绕过物体。如果偏转角不太大，会形成附体的斜激波：
斜激波
/
/β（激波角）
/
━━━━━━━━━━━━━→ Ma₁
/
/θ（偏转角）
/━━━━━━━━━━━━
楔形体
斜激波的特点：
偏转角θ、激波角β和上游马赫数$Ma_1$之间存在以下关系：
\[\tan\theta = 2\cot\beta \frac{Ma_1^2\sin^2\beta - 1}{Ma_1^2(\gamma + \cos2\beta) + 2}\]这个方程有两个解：
激波角β随偏转角θ的变化（Ma₁=2.0，空气）：
β(度)
90 ┤ 强激波解
│ ╱
70 ┤ ╱
│ ╱
50 ┤ ╱━━━━━━━━━ 弱激波解
│ ╱
30 ┤━━━━
└─────────────────────→ θ(度)
0 10 20 30 θmax
对于给定的上游马赫数，存在最大偏转角$\theta_{max}$：
| $Ma_1$ | 1.5 | 2.0 | 2.5 | 3.0 | 4.0 | 5.0 |
|---|---|---|---|---|---|---|
| $\theta_{max}$(度) | 11.8 | 23.0 | 32.2 | 39.3 | 48.0 | 53.5 |
当偏转角超过$\theta_{max}$时，斜激波无法附体，会形成脱体激波。
附体激波：
脱体激波：
附体激波： 脱体激波：
/ ╭─╮
/ │ │
/ │ │
/━━━━ ╰━━━╯
弓形激波
规则反射： 当斜激波遇到固壁时，会发生反射以满足壁面边界条件（流动与壁面平行）：
入射激波 反射激波
\ /
\ /
\/
━━━━━━━━━━━━━━━━ 壁面
马赫反射： 当入射角过大时，规则反射无法实现，会形成马赫反射：
入射激波
\
\ 马赫杆
\│
\│
━━━━━━━━━━━━━━━━ 壁面
滑移线
马赫反射的特征：
两道激波相交会产生复杂的波系结构：
同族激波相交（如两道压缩激波）：
异族激波相交（压缩波与膨胀波）：
超声速进气道： 利用多道斜激波逐级减速，比单道正激波总压损失小：
╱──╱──╱━━━━━━
╱ ╱ ╱ 燃烧室
╱ ╱ ╱
╱ ╱ ╱
━━━━━━━━━
多级斜激波压缩
楔形减阻： 超声速飞行器前缘采用楔形设计，产生斜激波而非脱体激波，显著降低波阻。
激波角估算法则：
与激波相反，当超声速流动遇到凸角时，需要膨胀以适应边界条件的变化。由于膨胀是等熵过程，不会形成间断面，而是产生一系列膨胀波（马赫波）：
超声速流动绕凸角的膨胀：
━━━━━━━━━━━━━→ Ma₁
╲╲╲╲╲
╲╲╲╲╲ 膨胀扇
╲╲╲╲╲
╲╲╲╲╲
━━━━━━━━━→ Ma₂ > Ma₁
θ
膨胀波的特征：
P-M函数描述了马赫数与流动偏转角的关系：
\[\nu(Ma) = \sqrt{\frac{\gamma+1}{\gamma-1}} \arctan\sqrt{\frac{\gamma-1}{\gamma+1}(Ma^2-1)} - \arctan\sqrt{Ma^2-1}\]对于空气（$\gamma=1.4$）：
| Ma | 1.0 | 1.5 | 2.0 | 2.5 | 3.0 | 4.0 | 5.0 | ∞ |
|---|---|---|---|---|---|---|---|---|
| ν(°) | 0 | 11.9 | 26.4 | 39.1 | 49.8 | 65.8 | 76.9 | 130.5 |
膨胀角计算： \(\theta = \nu(Ma_2) - \nu(Ma_1)\)
理论上，流动可以从任意马赫数膨胀到$Ma = \infty$，对应的最大偏转角为：
\(\theta_{max} = \nu(\infty) - \nu(Ma_1) = 130.5° - \nu(Ma_1)\)（空气）
实际限制：
马赫角： 第一道和最后一道马赫波的角度： \(\mu_1 = \arcsin(1/Ma_1)\) \(\mu_2 = \arcsin(1/Ma_2)\)
膨胀扇的几何结构：
μ₁╲
╲
━━━━━━╲━━━━━━━
╲
╲μ₂
╲━━━━━━━
θ
中心膨胀： 所有马赫波起源于凸角顶点，形成扇形区域。
火箭发动机喷管出口常见的膨胀模式：
完全膨胀（设计工况）：
欠膨胀（高空工况）：
欠膨胀喷流：
━━━━┓
┃╲╲╲╲╲╲
━━━━┛ ╲╲╲╲╲╲ 膨胀扇
╲╲╲╲╲╲
过膨胀（低空工况）：
膨胀波穿过激波：
膨胀波与激波相互作用：
激波 ╱ 减弱的激波
╱╲╲╲╲╲╱
╱ ╲╲╲╲╲╱
╱ 膨胀波
膨胀波反射： 膨胀波遇到固壁反射为压缩波：
入射膨胀波 ╲╲╲
╲╲╲
━━━━━━━━━━━━━━━━━━ 壁面
╱╱╱
反射压缩波 ╱╱╱
超声速喷管设计： 利用P-M膨胀计算喷管型线：
收缩段 喉部 扩张段
╲ │ ╱
╲ │ ╱
╲ │ ╱
╲ │ ╱
╲│╱
设计要点：
超声速风洞设计：
激波消除技术： 战斗机进气道利用膨胀波削弱激波：
╱━━━╲ 膨胀槽
╱ ╲╲╲
╱ ╲╲╲ 膨胀波
╱ 减弱的激波
━━━━━━━━━━━━━━━━
小偏转角近似（$\theta < 15°$）： \(\frac{\Delta p}{p} \approx -\gamma Ma^2 \frac{\theta}{57.3}\)
马赫数变化： \(\frac{\Delta Ma}{Ma} \approx \frac{\theta}{57.3} \frac{Ma^2}{Ma^2-1}\)
这些近似在初步设计中非常有用。
当激波与边界层相遇时，会产生复杂的相互作用。激波造成的逆压梯度可能导致边界层分离：
激波-边界层干扰示意图：
入射激波 ╱
╱
╱ 分离激波
╱ ╱
╱ ╱ 再附激波
╱ ╱ ╱
━━━━━━━━━━━━━━━━━━━━━━
分离泡
分离过程：
强激波与边界层干扰时形成特征性的λ形激波系：
λ激波结构：
主激波
╱
╱
╱─── 分离激波
╱
━━━━━━━━━━━━━━━━
分离区
λ激波特点：
激波-边界层干扰区的压力分布：
壁面压力分布：
p/p₁
4 ┤ ┌─────── 激波后压力
│ ╱│
3 ┤ ╱ │
│ ╱ │
2 ┤ ──╱ │
│ │
1 ┤────────┘
└──────────────→ x
上游 干扰区 下游
特征长度：
激波-边界层干扰导致局部热流急剧增加：
热流分布：
q/q₀
6 ┤ ∧ 热流峰值
│ ╱ ╲
4 ┤ ╱ ╲
│ ╱ ╲
2 ┤──╱ ╲──
│
1 ┤
└──────────────→ x
热流增强机理：
根据激波来源和边界层状态：
II型干扰（入射激波）：
III型干扰（激波反射）：
IV型干扰（激波相交）：
马赫数影响：
雷诺数影响：
激波强度： 压力比阈值（产生分离）：
被动控制：
激波位置
↓
━━━━━━┬━━━━━━━
└─ 抽吸缝
∧ ∧ ∧ 涡流发生器
━━━━━━━━━━━━━━
增强边界层动量
主动控制：
超声速进气道：
火箭喷管：
高超声速飞行器：
分离判据： Korkegi准则： \(\frac{p_2}{p_1} < 1 + 0.3Ma_1\sqrt{\frac{C_f}{2}}\)
其中$C_f$为摩擦系数。
热流估算： 再附点热流放大系数： \(\frac{q_{peak}}{q_0} \approx 3 + 0.5(Ma_1 - 2)\)
这些经验公式在初步设计中广泛应用。
苏格兰工程师和物理学家朗肯是热力学的奠基人之一。1870年，他在研究声波传播时，首次推导出了激波前后的跃变关系。朗肯的贡献不仅在于数学推导，更在于他认识到激波是一种物理上可能存在的间断面。
朗肯的主要成就：
有趣的是，朗肯的激波理论最初受到质疑，因为当时还没有产生激波的实验装置。直到20年后，恩斯特·马赫的超声速弹丸实验才证实了朗肯的理论。
法国工程师雨果尼奥独立于朗肯，在1887年研究爆炸波时得出了类似的结论。他的独特贡献在于引入了”雨果尼奥曲线”概念，这成为研究激波和爆轰波的基本工具。
雨果尼奥的创新：
不幸的是，雨果尼奥36岁就英年早逝，但他的工作深刻影响了后来的激波研究。法国至今仍将激波关系式称为”雨果尼奥关系式”以纪念他。
1886年，恩斯特·马赫使用他发明的纹影技术，首次拍摄到了超声速弹丸产生的激波。这些照片不仅验证了朗肯-雨果尼奥理论，还揭示了斜激波、脱体激波等复杂现象。
马赫的实验装置巧妙而简单：
这些先驱者的工作奠定了现代激波理论的基础，他们的名字永远镌刻在流体力学的历史上。
在强激波（$Ma > 5$）和高空稀薄大气中，激波内部出现显著的非平衡效应：
振动激发滞后：
温度分布：
T
│ 平动温度
│ ╱────────
│ ╱
│ ╱ 振动温度
│ ╱ ╱─────────
│╱ ╱
└──────────────→ x
激波 弛豫区
离解与电离：
这些化学反应吸收能量，影响激波后的温度和压力。
高温高压下，理想气体假设失效：
比热变化： \(\gamma = \gamma(T, p)\)
高温下$\gamma$从1.4降至1.2甚至更低，影响激波强度。
高密度效应： 分子间作用力不可忽略，状态方程偏离理想气体： \(p = \rho RT(1 + B\rho + C\rho^2 + ...)\)
强激波后的高温气体产生强烈辐射：
辐射机制：
辐射预热： 辐射能量传递到激波前方，改变来流性质：
再入飞行器：
高超声速风洞：
爆炸与爆轰：
这些高级话题代表着激波研究的前沿，对高超声速飞行和空间探索至关重要。
激波与膨胀波是超声速流动中最基本的物理现象。通过本章学习，你应该掌握：
11.1 空气流（$Ma_1 = 2.5$）通过正激波，计算激波后的马赫数、压力比、温度比和密度比。
11.2 超声速气流（$Ma_1 = 3.0$）绕过15°楔角，求弱激波解的激波角β和下游马赫数$Ma_2$。
11.3 马赫数2.0的均匀流经过10°凸角膨胀，计算膨胀后的马赫数和压力比。
11.4 判断以下情况是否会产生激波诱导的边界层分离： 湍流边界层，$Ma_1 = 2.5$，入射激波造成的压力比$p_2/p_1 = 3.5$。
11.5 设计一个二级斜激波进气道，将$Ma = 3.0$的来流减速到$Ma = 1.5$。第一道激波偏转角10°，求第二道激波的偏转角。比较与单道正激波的总压恢复系数。
11.6 火箭喷管出口$Ma = 3.5$，压力50 kPa，环境压力10 kPa。描述喷流结构，估算第一个马赫盘的位置。
11.7 高超声速飞行器（$Ma = 6$）的控制面偏转20°。考虑真实气体效应，定性分析激波结构和热流分布的变化。
11.8 超声速风洞试验段要求$Ma = 2.5$的均匀流。喷管喉部到试验段的膨胀角是多少？如果喉部宽度10 cm，估算试验段宽度。
1. 混淆激波与声波
2. 忽视激波解的多值性
3. 膨胀波的理解偏差
4. 单位不一致
5. 忽略总压损失累积
6. 激波表使用错误
7. 低估热流集中
8. 忽视非设计工况
9. 边界层影响
10. 激波位置诊断
11. 分离检测
12. 性能问题排查
通过系统地应用这些检查项，可以避免常见错误，确保超声速流动系统的设计质量。处