第11章:激波与膨胀波
当超声速战斗机从头顶呼啸而过,那震耳欲聋的音爆让人印象深刻——这就是激波现象最直观的体现。激波不仅存在于军用航空,在火箭发射、高速列车隧道通过、甚至开香槟时都能观察到。本章将揭示这些剧烈压缩现象背后的物理机制,以及与之对应的膨胀过程。通过学习激波与膨胀波,你将理解超声速流动中最基本也最重要的物理现象,掌握估算激波强度、预测流动偏转角度的方法,并能识别实际工程中的激波问题。
11.1 正激波关系式
激波的形成机理
想象一列高速列车驶入隧道的瞬间。列车前方的空气来不及让开,被剧烈压缩,形成一道压力突跃——这就是激波的本质。在超声速流动中,扰动信息的传播速度(声速)小于流体速度,导致压缩波不断叠加,最终形成厚度仅几个分子自由程的间断面。
亚声速(Ma < 1): 超声速(Ma > 1):
扰动波可以向上游传播 扰动波被流动带向下游
← ← ← → → →
————————————→ ————————————→
流动方向 流动方向
压缩波叠加形成激波
从微观角度看,激波是分子剧烈碰撞的区域。在激波厚度(约3-4个分子平均自由程)内,分子速度分布严重偏离麦克斯韦分布,动能转化为内能,导致温度急剧上升。
朗肯-雨果尼奥条件
激波前后必须满足三个守恒定律:
- 质量守恒:$\rho_1 V_1 = \rho_2 V_2$
- 动量守恒:$p_1 + \rho_1 V_1^2 = p_2 + \rho_2 V_2^2$
- 能量守恒:$h_1 + \frac{V_1^2}{2} = h_2 + \frac{V_2^2}{2}$
其中下标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$ |
| $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 |
经验法则:
- 当$Ma_1 = 2$时,压力增加约4.5倍,这是超声速进气道设计的重要参考点
- 密度比有上限:$\frac{\rho_2}{\rho_1} \leq \frac{\gamma+1}{\gamma-1} = 6$(空气)
- 强激波后流动总是亚声速($Ma_2 < 1$)
激波的不可逆性
激波过程伴随着熵增,是典型的不可逆过程: $$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)$$ 工程影响:
- 超声速进气道的总压恢复系数通常为0.7-0.9
- 激波越强,总压损失越大,这是超声速飞行效率的主要限制因素
11.2 斜激波与激波角
楔形体产生的斜激波
当超声速气流遇到楔形体时,气流必须偏转以绕过物体。如果偏转角不太大,会形成附体的斜激波:
斜激波
/
/β(激波角)
/
━━━━━━━━━━━━━→ Ma₁
/
/θ(偏转角)
/━━━━━━━━━━━━
楔形体
斜激波的特点:
- 激波面与来流方向成一定角度(激波角β)
- 流动通过激波后偏转θ角
- 激波强度小于相同马赫数的正激波
- 可以保持下游流动为超声速
θ-β-M关系
偏转角θ、激波角β和上游马赫数$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 |
| $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}$时,斜激波无法附体,会形成脱体激波。
脱体激波与附体激波
附体激波:
- 激波起始于楔形体顶点
- 激波后流动方向与楔面平行
- 激波强度相对较弱
- 常见于尖锐前缘的超声速飞行器
脱体激波:
- 激波离开物体表面一定距离
- 形成弓形激波
- 驻点处为正激波,向外逐渐变为斜激波
- 常见于钝头体或大偏转角情况
附体激波: 脱体激波:
/ ╭─╮
/ │ │
/ │ │
/━━━━ ╰━━━╯
弓形激波
激波反射与相交
规则反射: 当斜激波遇到固壁时,会发生反射以满足壁面边界条件(流动与壁面平行):
入射激波 反射激波
\ /
\ /
\/
━━━━━━━━━━━━━━━━ 壁面
马赫反射: 当入射角过大时,规则反射无法实现,会形成马赫反射:
入射激波
\
\ 马赫杆
\│
\│
━━━━━━━━━━━━━━━━ 壁面
滑移线
马赫反射的特征:
- 出现垂直于壁面的马赫杆(正激波)
- 形成滑移线(速度间断面)
- 三激波交汇于三叉点
激波-激波相互作用
两道激波相交会产生复杂的波系结构:
同族激波相交(如两道压缩激波):
- 产生更强的激波
- 可能产生滑移线
- 应用:激波聚焦产生高温高压
异族激波相交(压缩波与膨胀波):
- 相互削弱
- 可能完全抵消
- 应用:激波消除技术
工程应用实例
超声速进气道: 利用多道斜激波逐级减速,比单道正激波总压损失小:
╱──╱──╱━━━━━━
╱ ╱ ╱ 燃烧室
╱ ╱ ╱
╱ ╱ ╱
━━━━━━━━━
多级斜激波压缩
楔形减阻: 超声速飞行器前缘采用楔形设计,产生斜激波而非脱体激波,显著降低波阻。
激波角估算法则:
- 弱激波:$\beta \approx \theta + \frac{20°}{Ma_1}$
- 音速线位置:激波后马赫数为1的点约在$\beta = 90°/\sqrt{Ma_1}$
11.3 普朗特-迈耶膨胀
膨胀波的形成
与激波相反,当超声速流动遇到凸角时,需要膨胀以适应边界条件的变化。由于膨胀是等熵过程,不会形成间断面,而是产生一系列膨胀波(马赫波):
超声速流动绕凸角的膨胀:
━━━━━━━━━━━━━→ 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 | ∞ |
| 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膨胀计算喷管型线:
收缩段 喉部 扩张段
╲ │ ╱
╲ │ ╱
╲ │ ╱
╲ │ ╱
╲│╱
设计要点:
- 喉部马赫数 = 1
- 出口马赫数由面积比确定
- 型线由P-M函数确定
超声速风洞设计:
- 利用膨胀角设计拐角
- 产生均匀超声速流场
- 典型设计:$Ma = 2-5$的风洞喷管扩张角15-25°
激波消除技术: 战斗机进气道利用膨胀波削弱激波:
╱━━━╲ 膨胀槽
╱ ╲╲╲
╱ ╲╲╲ 膨胀波
╱ 减弱的激波
━━━━━━━━━━━━━━━━
简化估算方法
小偏转角近似($\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}$$ 这些近似在初步设计中非常有用。
11.4 激波-边界层干扰
激波诱导分离的物理机制
当激波与边界层相遇时,会产生复杂的相互作用。激波造成的逆压梯度可能导致边界层分离:
激波-边界层干扰示意图:
入射激波 ╱
╱
╱ 分离激波
╱ ╱
╱ ╱ 再附激波
╱ ╱ ╱
━━━━━━━━━━━━━━━━━━━━━━
分离泡
分离过程:
- 激波产生强逆压梯度
- 边界层内低动量流体无法克服逆压
- 流动分离形成回流区
- 产生分离激波和再附激波
λ激波结构
强激波与边界层干扰时形成特征性的λ形激波系:
λ激波结构:
主激波
╱
╱
╱─── 分离激波
╱
━━━━━━━━━━━━━━━━
分离区
λ激波特点:
- 分离激波角度较小(弱激波)
- 主激波与分离激波在边界层外缘相交
- 交点高度约为边界层厚度的3-5倍
- 壁面压力呈台阶状分布
压力分布特征
激波-边界层干扰区的压力分布:
壁面压力分布:
p/p₁
4 ┤ ┌─────── 激波后压力
│ ╱│
3 ┤ ╱ │
│ ╱ │
2 ┤ ──╱ │
│ │
1 ┤────────┘
└──────────────→ x
上游 干扰区 下游
特征长度:
- 上游影响距离:$L_u \approx 3\delta$($\delta$为边界层厚度)
- 分离区长度:$L_s \approx (2-5)\delta$
- 压力上升距离:$L_p \approx (5-10)\delta$
热流集中现象
激波-边界层干扰导致局部热流急剧增加:
热流分布:
q/q₀
6 ┤ ∧ 热流峰值
│ ╱ ╲
4 ┤ ╱ ╲
│ ╱ ╲
2 ┤──╱ ╲──
│
1 ┤
└──────────────→ x
热流增强机理:
- 再附点附近薄边界层
- 湍流度增加
- 压力和温度同时升高
- 峰值热流可达无干扰时的3-8倍
干扰类型分类
根据激波来源和边界层状态:
II型干扰(入射激波):
- 外部激波入射到壁面
- 常见于进气道唇口
- 分离区相对较大
III型干扰(激波反射):
- 壁面产生的激波反射
- 常见于喷管和扩压器
- 可能产生多次反射
IV型干扰(激波相交):
- 多道激波相交于边界层
- 最复杂的干扰形式
- 可能产生严重的热流集中
影响因素
马赫数影响:
- $Ma > 2$:开始出现明显分离
- $Ma = 3-4$:分离区最大
- $Ma > 5$:真实气体效应显著
雷诺数影响:
- 高雷诺数:分离区较小
- 低雷诺数:容易产生大规模分离
- 转捩位置影响分离特性
激波强度: 压力比阈值(产生分离):
- 层流:$p_2/p_1 > 1.5-2.0$
- 湍流:$p_2/p_1 > 2.5-3.0$
控制方法
被动控制:
- 边界层抽吸:
激波位置
↓
━━━━━━┬━━━━━━━
└─ 抽吸缝
- 涡流发生器:
∧ ∧ ∧ 涡流发生器
━━━━━━━━━━━━━━
增强边界层动量
- 多孔壁面: - 减弱激波强度 - 降低逆压梯度 - 延迟分离
主动控制:
-
合成射流: - 周期性吹吸 - 增强掺混 - 抑制分离
-
等离子体控制: - 局部加热 - 改变激波结构 - 减小分离区
工程实例
超声速进气道:
- 唇口激波与边界层干扰
- 可能导致进气道不起动
- 设计泄流槽避免分离
火箭喷管:
- 过膨胀时产生分离
- 可能导致侧向力
- 采用高度补偿喷管设计
高超声速飞行器:
- 控制面铰链处热流集中
- 采用间隙密封设计
- 热防护系统重点区域
设计准则
分离判据: 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)$$ 这些经验公式在初步设计中广泛应用。
历史人物:朗肯与雨果尼奥的激波理论
威廉·约翰·麦夸恩·朗肯(1820-1872)
苏格兰工程师和物理学家朗肯是热力学的奠基人之一。1870年,他在研究声波传播时,首次推导出了激波前后的跃变关系。朗肯的贡献不仅在于数学推导,更在于他认识到激波是一种物理上可能存在的间断面。
朗肯的主要成就:
- 建立了激波的守恒方程组
- 预言了激波的不可逆性
- 首次计算了激波前后的参数关系
有趣的是,朗肯的激波理论最初受到质疑,因为当时还没有产生激波的实验装置。直到20年后,恩斯特·马赫的超声速弹丸实验才证实了朗肯的理论。
皮埃尔·亨利·雨果尼奥(1851-1887)
法国工程师雨果尼奥独立于朗肯,在1887年研究爆炸波时得出了类似的结论。他的独特贡献在于引入了"雨果尼奥曲线"概念,这成为研究激波和爆轰波的基本工具。
雨果尼奥的创新:
- 发展了图解法分析激波
- 研究了爆轰波的传播
- 建立了激波稳定性理论
不幸的是,雨果尼奥36岁就英年早逝,但他的工作深刻影响了后来的激波研究。法国至今仍将激波关系式称为"雨果尼奥关系式"以纪念他。
激波理论的实验验证
1886年,恩斯特·马赫使用他发明的纹影技术,首次拍摄到了超声速弹丸产生的激波。这些照片不仅验证了朗肯-雨果尼奥理论,还揭示了斜激波、脱体激波等复杂现象。
马赫的实验装置巧妙而简单:
- 使用火花放电产生瞬时照明
- 利用密度变化造成的光线偏折
- 首次实现激波的可视化
这些先驱者的工作奠定了现代激波理论的基础,他们的名字永远镌刻在流体力学的历史上。
高级话题:非平衡激波与真实气体效应
非平衡现象
在强激波($Ma > 5$)和高空稀薄大气中,激波内部出现显著的非平衡效应:
振动激发滞后:
- 平动温度瞬间升高
- 振动温度缓慢达到平衡
- 弛豫区长度可达激波厚度的1000倍
温度分布:
T
│ 平动温度
│ ╱────────
│ ╱
│ ╱ 振动温度
│ ╱ ╱─────────
│╱ ╱
└──────────────→ x
激波 弛豫区
离解与电离:
- $T > 2000K$:氧分子开始离解
- $T > 4000K$:氮分子离解
- $T > 9000K$:开始电离
这些化学反应吸收能量,影响激波后的温度和压力。
真实气体效应
高温高压下,理想气体假设失效:
比热变化: $$\gamma = \gamma(T, p)$$ 高温下$\gamma$从1.4降至1.2甚至更低,影响激波强度。
高密度效应: 分子间作用力不可忽略,状态方程偏离理想气体: $$p = \rho RT(1 + B\rho + C\rho^2 + ...)$$
辐射效应
强激波后的高温气体产生强烈辐射:
辐射机制:
- 连续谱辐射(轫致辐射)
- 线谱辐射(原子跃迁)
- 分子带辐射
辐射预热: 辐射能量传递到激波前方,改变来流性质:
- 降低激波强度
- 增加激波厚度
- 改变激波结构
应用实例
再入飞行器:
- 激波层温度超过10000K
- 强烈的辐射加热
- 需要烧蚀热防护
高超声速风洞:
- 真实气体效应模拟
- 需要高焓设备
- 激波管、电弧风洞
爆炸与爆轰:
- 极强激波($Ma > 10$)
- 化学反应耦合
- 爆轰波传播
这些高级话题代表着激波研究的前沿,对高超声速飞行和空间探索至关重要。
本章小结
激波与膨胀波是超声速流动中最基本的物理现象。通过本章学习,你应该掌握:
核心概念
- 正激波:垂直于流动方向的激波,造成最强的压缩和总压损失
- 斜激波:与流动成一定角度,可保持下游超声速,存在弱解和强解
- 膨胀波:等熵膨胀过程,形成扇形马赫波区域
- 激波-边界层干扰:产生分离、热流集中等复杂现象
关键公式
- 正激波压力比:$\frac{p_2}{p_1} = \frac{2\gamma Ma_1^2 - (\gamma-1)}{\gamma+1}$
- θ-β-M关系:连接偏转角、激波角和马赫数
- P-M函数:$\nu(Ma)$描述膨胀过程
- 分离判据:$\frac{p_2}{p_1} < 1 + 0.3Ma_1\sqrt{\frac{C_f}{2}}$
工程应用要点
- 超声速进气道设计:多级斜激波压缩
- 喷管设计:利用膨胀波加速到设计马赫数
- 热防护:关注激波-边界层干扰区的热流集中
- 激波控制:涡流发生器、边界层抽吸等方法
经验法则
- 密度比上限:空气中约为6
- 弱激波角估算:$\beta \approx \theta + 20°/Ma_1$
- 热流放大:再附点可达3-8倍
- 分离压力比阈值:层流1.5-2.0,湍流2.5-3.0
练习题
基础题
11.1 空气流($Ma_1 = 2.5$)通过正激波,计算激波后的马赫数、压力比、温度比和密度比。
提示
使用正激波关系式或查正激波表。注意空气的比热比$\gamma = 1.4$。
答案
$Ma_2 = 0.513$,$p_2/p_1 = 7.125$,$T_2/T_1 = 2.138$,$\rho_2/\rho_1 = 3.333$
11.2 超声速气流($Ma_1 = 3.0$)绕过15°楔角,求弱激波解的激波角β和下游马赫数$Ma_2$。
提示
使用θ-β-M关系或查斜激波图表。弱激波解对应较小的β值。
答案
$\beta \approx 32.2°$,$Ma_2 \approx 2.28$
11.3 马赫数2.0的均匀流经过10°凸角膨胀,计算膨胀后的马赫数和压力比。
提示
使用P-M函数:先计算$\nu(Ma_1)$,加上偏转角得$\nu(Ma_2)$,反求$Ma_2$。
答案
$Ma_2 \approx 2.38$,$p_2/p_1 \approx 0.58$
11.4 判断以下情况是否会产生激波诱导的边界层分离: 湍流边界层,$Ma_1 = 2.5$,入射激波造成的压力比$p_2/p_1 = 3.5$。
提示
湍流边界层分离的临界压力比约为2.5-3.0。
答案
会产生分离。压力比3.5超过了湍流边界层的典型分离阈值(约3.0)。
挑战题
11.5 设计一个二级斜激波进气道,将$Ma = 3.0$的来流减速到$Ma = 1.5$。第一道激波偏转角10°,求第二道激波的偏转角。比较与单道正激波的总压恢复系数。
提示
分步计算:第一道斜激波后的状态→第二道斜激波的偏转角→总压损失。记住总压恢复系数是各级的乘积。
答案
第一道激波后:$Ma \approx 2.3$,第二道激波偏转角约11°。 二级斜激波总压恢复:约0.78 单道正激波总压恢复:约0.33 二级设计显著改善了总压恢复。
11.6 火箭喷管出口$Ma = 3.5$,压力50 kPa,环境压力10 kPa。描述喷流结构,估算第一个马赫盘的位置。
提示
这是欠膨胀喷流。喷流先膨胀,然后通过激波调整到环境压力。马赫盘位置约为喷管直径的0.67倍乘以压力比的平方根。
答案
欠膨胀喷流形成钻石形激波结构。压力比5:1导致强烈膨胀。 第一个马赫盘距离:$x/D \approx 0.67\sqrt{5} \approx 1.5$ 喷流呈现周期性的膨胀-压缩结构。
11.7 高超声速飞行器($Ma = 6$)的控制面偏转20°。考虑真实气体效应,定性分析激波结构和热流分布的变化。
提示
$Ma = 6$时开始出现显著的真实气体效应。考虑振动激发、离解对激波强度的影响,以及激波-边界层干扰造成的热流集中。
答案
20°偏转在$Ma = 6$时仍能形成附体斜激波,但接近脱体极限。 真实气体效应:
- 振动激发和离解吸收能量,激波后温度低于理想气体预测
- 有效$\gamma$降低,激波角增大
- 激波-边界层干扰产生λ激波
- 再附点热流可达层流值的5-10倍
- 需要特殊的热防护设计
11.8 超声速风洞试验段要求$Ma = 2.5$的均匀流。喷管喉部到试验段的膨胀角是多少?如果喉部宽度10 cm,估算试验段宽度。
提示
使用P-M函数计算从$Ma = 1$到$Ma = 2.5$的总膨胀角。面积比可用$A/A^* = \frac{1}{Ma}[(1+\frac{\gamma-1}{2}Ma^2)/(\frac{\gamma+1}{2})]^{\frac{\gamma+1}{2(\gamma-1)}}$
答案
P-M膨胀角:$\nu(2.5) - \nu(1) = 39.1° - 0° = 39.1°$ 单侧膨胀角:19.55° 面积比:$A/A^* \approx 2.64$ 假设二维喷管,试验段宽度约26.4 cm 实际设计会采用渐变型线以获得均匀流场。
常见陷阱与错误
概念误区
-
混淆激波与声波 - ❌ 错误:激波是大振幅声波 - ✅ 正确:激波是非线性间断,声波是线性小扰动
-
忽视激波解的多值性 - ❌ 错误:斜激波只有唯一解 - ✅ 正确:存在弱解和强解,需根据下游条件确定
-
膨胀波的理解偏差 - ❌ 错误:膨胀也会产生间断 - ✅ 正确:膨胀是等熵连续过程,通过马赫波扇实现
计算错误
-
单位不一致 - ❌ 错误:混用度和弧度 - ✅ 正确:P-M函数通常用度,三角函数用弧度,需要转换
-
忽略总压损失累积 - ❌ 错误:多级激波的总压恢复是各级之和 - ✅ 正确:总压恢复系数是各级的乘积
-
激波表使用错误 - ❌ 错误:用激波前马赫数直接查下游参数 - ✅ 正确:斜激波需要先计算法向马赫数分量
工程判断失误
-
低估热流集中 - ❌ 错误:用平均热流设计热防护 - ✅ 正确:考虑激波干扰区局部热流可达5-10倍
-
忽视非设计工况 - ❌ 错误:只考虑设计马赫数 - ✅ 正确:检查起动特性和非设计点性能
-
边界层影响 - ❌ 错误:用无粘流理论设计进气道 - ✅ 正确:考虑激波-边界层干扰,预留分离裕度
调试技巧
-
激波位置诊断 - 压力分布突跃 - 纹影/阴影照相 - 压力敏感涂料
-
分离检测 - 壁面压力平台 - 表面油流显示 - 热流分布异常
-
性能问题排查 - 总压恢复低→检查激波系配置 - 流动不均匀→检查激波反射/相交 - 振荡/喘振→检查激波-边界层干扰
最佳实践检查清单
设计阶段
激波系设计
- [ ] 选择合适的激波配置(正激波/斜激波/混合)
- [ ] 验证所有偏转角小于最大偏转角
- [ ] 计算总压损失并满足要求
- [ ] 检查激波相交和反射模式
膨胀区设计
- [ ] 使用P-M方法设计膨胀轮廓
- [ ] 避免过度膨胀(考虑真实气体极限)
- [ ] 检查膨胀波反射影响
- [ ] 验证出口流场均匀性
激波-边界层干扰
- [ ] 识别所有潜在干扰区域
- [ ] 估算分离区大小
- [ ] 计算局部热流峰值
- [ ] 设计必要的控制措施
分析阶段
性能评估
- [ ] 计算设计点性能
- [ ] 分析非设计点特性
- [ ] 评估起动/不起动边界
- [ ] 检查瞬态响应
载荷评估
- [ ] 计算压力载荷分布
- [ ] 识别热流集中区域
- [ ] 评估非定常载荷
- [ ] 考虑疲劳和热应力
验证阶段
数值验证
- [ ] 网格无关性检查(激波捕捉)
- [ ] 湍流模型敏感性
- [ ] 与理论/经验公式对比
- [ ] 极限工况检查
试验验证
- [ ] 风洞试验方案设计
- [ ] 测量技术选择(压力、纹影、PSP等)
- [ ] 雷诺数/马赫数相似性
- [ ] 不确定度分析
优化阶段
性能优化
- [ ] 总压恢复最大化
- [ ] 流场畸变最小化
- [ ] 起动裕度优化
- [ ] 宽速域性能平衡
结构优化
- [ ] 热防护优化
- [ ] 结构重量优化
- [ ] 控制系统简化
- [ ] 制造可行性
通过系统地应用这些检查项,可以避免常见错误,确保超声速流动系统的设计质量。处