第11章:激波与膨胀波

当超声速战斗机从头顶呼啸而过,那震耳欲聋的音爆让人印象深刻——这就是激波现象最直观的体现。激波不仅存在于军用航空,在火箭发射、高速列车隧道通过、甚至开香槟时都能观察到。本章将揭示这些剧烈压缩现象背后的物理机制,以及与之对应的膨胀过程。通过学习激波与膨胀波,你将理解超声速流动中最基本也最重要的物理现象,掌握估算激波强度、预测流动偏转角度的方法,并能识别实际工程中的激波问题。

11.1 正激波关系式

激波的形成机理

想象一列高速列车驶入隧道的瞬间。列车前方的空气来不及让开,被剧烈压缩,形成一道压力突跃——这就是激波的本质。在超声速流动中,扰动信息的传播速度(声速)小于流体速度,导致压缩波不断叠加,最终形成厚度仅几个分子自由程的间断面。

亚声速(Ma < 1):         超声速(Ma > 1):

   扰动波可以向上游传播        扰动波被流动带向下游
        ←  ←  ←                    → → →
    ————————————→           ————————————→
        流动方向                 流动方向
                              压缩波叠加形成激波

从微观角度看,激波是分子剧烈碰撞的区域。在激波厚度(约3-4个分子平均自由程)内,分子速度分布严重偏离麦克斯韦分布,动能转化为内能,导致温度急剧上升。

朗肯-雨果尼奥条件

激波前后必须满足三个守恒定律:

  1. 质量守恒:$\rho_1 V_1 = \rho_2 V_2$
  2. 动量守恒:$p_1 + \rho_1 V_1^2 = p_2 + \rho_2 V_2^2$
  3. 能量守恒:$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 激波-边界层干扰

激波诱导分离的物理机制

当激波与边界层相遇时,会产生复杂的相互作用。激波造成的逆压梯度可能导致边界层分离:

激波-边界层干扰示意图:

入射激波 ╱
        ╱
       ╱  分离激波
      ╱  ╱
     ╱  ╱ 再附激波
    ╱  ╱ ╱
━━━━━━━━━━━━━━━━━━━━━━
    分离泡

分离过程:

  1. 激波产生强逆压梯度
  2. 边界层内低动量流体无法克服逆压
  3. 流动分离形成回流区
  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$

控制方法

被动控制

  1. 边界层抽吸
    激波位置
       ↓
━━━━━━┬━━━━━━━
      └─ 抽吸缝
  1. 涡流发生器
  ∧ ∧ ∧  涡流发生器
━━━━━━━━━━━━━━
  增强边界层动量
  1. 多孔壁面: - 减弱激波强度 - 降低逆压梯度 - 延迟分离

主动控制

  1. 合成射流: - 周期性吹吸 - 增强掺混 - 抑制分离

  2. 等离子体控制: - 局部加热 - 改变激波结构 - 减小分离区

工程实例

超声速进气道

  • 唇口激波与边界层干扰
  • 可能导致进气道不起动
  • 设计泄流槽避免分离

火箭喷管

  • 过膨胀时产生分离
  • 可能导致侧向力
  • 采用高度补偿喷管设计

高超声速飞行器

  • 控制面铰链处热流集中
  • 采用间隙密封设计
  • 热防护系统重点区域

设计准则

分离判据: 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$)
  • 化学反应耦合
  • 爆轰波传播

这些高级话题代表着激波研究的前沿,对高超声速飞行和空间探索至关重要。

本章小结

激波与膨胀波是超声速流动中最基本的物理现象。通过本章学习,你应该掌握:

核心概念

  1. 正激波:垂直于流动方向的激波,造成最强的压缩和总压损失
  2. 斜激波:与流动成一定角度,可保持下游超声速,存在弱解和强解
  3. 膨胀波:等熵膨胀过程,形成扇形马赫波区域
  4. 激波-边界层干扰:产生分离、热流集中等复杂现象

关键公式

  • 正激波压力比:$\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 实际设计会采用渐变型线以获得均匀流场。

常见陷阱与错误

概念误区

  1. 混淆激波与声波 - ❌ 错误:激波是大振幅声波 - ✅ 正确:激波是非线性间断,声波是线性小扰动

  2. 忽视激波解的多值性 - ❌ 错误:斜激波只有唯一解 - ✅ 正确:存在弱解和强解,需根据下游条件确定

  3. 膨胀波的理解偏差 - ❌ 错误:膨胀也会产生间断 - ✅ 正确:膨胀是等熵连续过程,通过马赫波扇实现

计算错误

  1. 单位不一致 - ❌ 错误:混用度和弧度 - ✅ 正确:P-M函数通常用度,三角函数用弧度,需要转换

  2. 忽略总压损失累积 - ❌ 错误:多级激波的总压恢复是各级之和 - ✅ 正确:总压恢复系数是各级的乘积

  3. 激波表使用错误 - ❌ 错误:用激波前马赫数直接查下游参数 - ✅ 正确:斜激波需要先计算法向马赫数分量

工程判断失误

  1. 低估热流集中 - ❌ 错误:用平均热流设计热防护 - ✅ 正确:考虑激波干扰区局部热流可达5-10倍

  2. 忽视非设计工况 - ❌ 错误:只考虑设计马赫数 - ✅ 正确:检查起动特性和非设计点性能

  3. 边界层影响 - ❌ 错误:用无粘流理论设计进气道 - ✅ 正确:考虑激波-边界层干扰,预留分离裕度

调试技巧

  1. 激波位置诊断 - 压力分布突跃 - 纹影/阴影照相 - 压力敏感涂料

  2. 分离检测 - 壁面压力平台 - 表面油流显示 - 热流分布异常

  3. 性能问题排查 - 总压恢复低→检查激波系配置 - 流动不均匀→检查激波反射/相交 - 振荡/喘振→检查激波-边界层干扰

最佳实践检查清单

设计阶段

激波系设计

  • [ ] 选择合适的激波配置(正激波/斜激波/混合)
  • [ ] 验证所有偏转角小于最大偏转角
  • [ ] 计算总压损失并满足要求
  • [ ] 检查激波相交和反射模式

膨胀区设计

  • [ ] 使用P-M方法设计膨胀轮廓
  • [ ] 避免过度膨胀(考虑真实气体极限)
  • [ ] 检查膨胀波反射影响
  • [ ] 验证出口流场均匀性

激波-边界层干扰

  • [ ] 识别所有潜在干扰区域
  • [ ] 估算分离区大小
  • [ ] 计算局部热流峰值
  • [ ] 设计必要的控制措施

分析阶段

性能评估

  • [ ] 计算设计点性能
  • [ ] 分析非设计点特性
  • [ ] 评估起动/不起动边界
  • [ ] 检查瞬态响应

载荷评估

  • [ ] 计算压力载荷分布
  • [ ] 识别热流集中区域
  • [ ] 评估非定常载荷
  • [ ] 考虑疲劳和热应力

验证阶段

数值验证

  • [ ] 网格无关性检查(激波捕捉)
  • [ ] 湍流模型敏感性
  • [ ] 与理论/经验公式对比
  • [ ] 极限工况检查

试验验证

  • [ ] 风洞试验方案设计
  • [ ] 测量技术选择(压力、纹影、PSP等)
  • [ ] 雷诺数/马赫数相似性
  • [ ] 不确定度分析

优化阶段

性能优化

  • [ ] 总压恢复最大化
  • [ ] 流场畸变最小化
  • [ ] 起动裕度优化
  • [ ] 宽速域性能平衡

结构优化

  • [ ] 热防护优化
  • [ ] 结构重量优化
  • [ ] 控制系统简化
  • [ ] 制造可行性

通过系统地应用这些检查项,可以避免常见错误,确保超声速流动系统的设计质量。处