intuitive_cfd

第12章:超声速流动设计

当飞行器速度超过声速时,流动的物理特性发生了根本性的改变。信息不再能够向上游传播,激波成为流场的主导特征,传统的亚声速设计理念完全失效。从V-2火箭到现代高超声速飞行器,超声速流动设计一直是航空航天工程的核心挑战。本章将介绍超声速流动中的基本设计原理,通过简单几何体的激波关系来理解复杂飞行器的设计思路。

12.1 楔形体与锥形体绕流

12.1.1 二维楔形体绕流

当超声速气流遇到楔形体时,会在楔形体前缘形成附体斜激波。这是超声速流动中最基本的流动模式。

        激波角 β
            ↓
    ━━━━━━╱━━━━━━
         ╱│
   Ma₁  ╱ │ Ma₂
       ╱  │
      ╱   │ 楔角 θ
     ╱    │
    ╱_____│

θ-β-M 关系(激波角关系式):

\[\tan \theta = 2 \cot \beta \frac{M_1^2 \sin^2 \beta - 1}{M_1^2 (\gamma + \cos 2\beta) + 2}\]

对于给定的来流马赫数 $M_1$ 和楔角 $θ$,存在两个可能的激波角:

自然界总是选择弱激波解,除非几何约束迫使形成强激波(脱体激波)。

最大楔角: 对于每个马赫数,存在一个最大楔角 $θ_{max}$,超过此角度将形成脱体激波:

Ma 1.5 2.0 3.0 5.0
θmax 12° 23° 34° 42°

12.1.2 轴对称锥形体绕流

锥形体绕流比楔形体复杂,因为流动具有三维特性。锥形激波是弯曲的,激波后的流动参数沿径向变化。

        锥形激波
           ╱╲
          ╱  ╲
    Ma₁  ╱    ╲  
        ╱ cone ╲
       ╱   θc   ╲
      ╱__________╲

Taylor-Maccoll 方程描述锥形流场,但工程上常用经验关系:

锥角-激波角关系的经验公式: 对于相同的偏转角,锥形体的激波角总是小于楔形体:

\[\beta_{cone} \approx \beta_{wedge} - \Delta\beta\]

其中修正量: \(\Delta\beta \approx \frac{3°~5°}{Ma_1 - 1}\)

这意味着:

12.1.3 工程应用实例

战斗机进气道设计: F-104 战斗机采用尖锐的进气锥,通过多级锥面逐步减速超声速气流:

超声速子弹设计

12.2 激波阻力的估算

12.2.1 波阻的物理本质

超声速飞行中,即使在无粘流动中也存在阻力——波阻。这是因为通过激波时熵增,动能不可逆地转化为热能。

波阻系数的基本形式

对于细长体,波阻系数可近似为:

\[C_{D,wave} = \frac{4}{\pi} \int_0^L \left(\frac{dS}{dx}\right)^2 dx\]

其中 $S(x)$ 是横截面积分布。

12.2.2 简单形状的波阻

二维楔形体: \(C_{D,wave} = \frac{4 \theta^2}{\sqrt{M_1^2 - 1}}\)

圆锥: \(C_{D,wave} = \frac{2 \theta_c^2}{\sqrt{M_1^2 - 1}}\)

钝头体: \(C_{D,wave} \approx 0.85 + \frac{0.35}{M_1^2}\)

12.2.3 面积律(Area Rule)

Whitcomb面积律(1950年代): 跨声速飞行器的波阻主要取决于横截面积的纵向分布,而不是具体形状。

理想的面积分布应该是Sears-Haack体

\[S(x) = S_{max} \left[1 - \left(\frac{2x}{L} - 1\right)^2\right]^{3/2}\]

工程应用

12.2.4 波阻的快速估算方法

经验公式集合

  1. 细长旋转体(长细比 > 5): \(C_{D,wave} \approx \frac{1.2 (d/L)^2}{\sqrt{M^2 - 1}}\)

  2. 机翼波阻(薄翼): \(C_{D,wave} \approx \frac{4 (t/c)^2}{\sqrt{M^2 - 1}}\)

  3. 组合体波阻(干扰因子): \(C_{D,total} = C_{D,body} + C_{D,wing} + K_{interference} \cdot C_{D,wing}\)

    其中干扰因子 $K_{interference} \approx 0.1 - 0.3$

12.3 喷管设计原理

12.3.1 拉瓦尔喷管基础

超声速流动的产生需要收缩-扩张型喷管(de Laval nozzle)。这个看似简单的装置蕴含着深刻的物理原理。

     收缩段    喉道   扩张段
        ╲      │      ╱
         ╲     │     ╱
    ━━━━━━╲    │    ╱━━━━━━
           ╲   │   ╱
            ╲  │  ╱   Ma > 1
     Ma < 1  ╲ │ ╱
              ╲│╱
               ○  Ma = 1

面积-马赫数关系

\[\frac{A}{A^*} = \frac{1}{M} \left[\frac{2}{\gamma + 1} \left(1 + \frac{\gamma - 1}{2} M^2\right)\right]^{\frac{\gamma + 1}{2(\gamma - 1)}}\]

其中 $A^*$ 是喉道面积(声速点)。

12.3.2 喷管轮廓设计

理想喷管轮廓应满足:

  1. 无激波(等熵膨胀)
  2. 出口流动均匀平行
  3. 长度最短

特征线法(Method of Characteristics): 喷管扩张段可分为两个区域:

  1. 初始膨胀区:从喉道下游开始的快速膨胀
  2. 矫直区:将径向速度分量转为轴向

实用设计公式(Rao喷管):

喷管半径分布: \(r(x) = r_t + (r_e - r_t) \left(\frac{x}{L}\right)^n\)

其中指数 $n$ 的选择:

12.3.3 非设计工况

实际喷管经常工作在非设计压比下:

欠膨胀(背压 < 设计压力):

过膨胀(背压 > 设计压力):

流动分离准则(Summerfield准则): \(\frac{p_{sep}}{p_0} \approx 0.4 \left(\frac{p_a}{p_0}\right)^{0.85}\)

12.3.4 特殊喷管设计

塞式喷管(Aerospike)

矢量喷管

12.4 超声速进气道

12.4.1 进气道的功能与挑战

超声速进气道必须将高速气流减速至亚声速(Ma ≈ 0.3-0.5)供发动机使用,同时:

12.4.2 进气道类型

皮托进气道(正激波进气道):

    ━━━━┃━━━━
   Ma>1 ┃ Ma<1
    ━━━━┃━━━━
        正激波

外压进气道(斜激波系):

       ╱━━━━━
      ╱╱━━━━━
     ╱╱╱━━━━━
    ╱╱╱╱

混压进气道

12.4.3 起动问题

超声速进气道存在”起动”和”不起动”两种状态:

不起动状态

起动条件(Kantrowitz准则): \(\frac{A_{throat}}{A_{capture}} > \left(\frac{A}{A^*}\right)_{M_1}\)

改善起动的方法

  1. 可变几何喉道
  2. 边界层吸除
  3. 辅助进气门

12.4.4 进气道-发动机匹配

流量匹配: \(\dot{m}_{inlet} = \dot{m}_{engine}\)

不匹配会导致:

动态响应: 进气道激波系统对扰动敏感,需要:

12.5 历史人物:冯·布劳恩与V-2火箭

12.5.1 早期梦想与挫折

韦纳·冯·布劳恩(Wernher von Braun, 1912-1977)从小就梦想太空旅行。1930年,18岁的他加入柏林工业大学,师从火箭先驱赫尔曼·奥伯特。他的博士论文《液体推进剂火箭发动机的设计、理论和实验研究》奠定了现代火箭推进的理论基础。

12.5.2 V-2火箭的技术突破

V-2火箭(1944)是第一种实用的弹道导弹,其超声速设计包含多项创新:

空气动力学设计

关键参数

12.5.3 技术遗产

V-2的设计理念影响了后续所有弹道导弹和运载火箭:

  1. 级间分离技术
  2. 涡轮泵供给系统
  3. 陀螺稳定平台
  4. 再入空气动力学

冯·布劳恩后来主持了美国土星五号火箭的研制,将人类送上月球。他的名言:”基础研究是我不知道自己在做什么的时候所做的事情”,体现了探索精神。

12.6 高级话题:乘波体设计与高超声速飞行器

12.6.1 乘波体概念

乘波体(Waverider)是一种利用自身产生的激波面作为”空气动力学下表面”的飞行器构型。

     顶视图              侧视图
    ╱─────╲            ━━━━━━━━
   ╱       ╲              ╲___
  ╱         ╲          激波附体

设计原理

  1. 从已知流场(如楔形流或锥形流)出发
  2. 沿流线追踪设计表面
  3. 前缘与激波面重合
  4. 激波完全附体,无泄漏

12.6.2 性能优势

高升阻比

设计马赫数的选择

12.6.3 实际挑战

非设计点性能

热防护

粘性效应

12.6.4 发展趋势

概念验证飞行器

未来应用

  1. 高超声速巡航导弹:1小时全球打击
  2. 可重复使用空天飞机:降低发射成本
  3. 高速民航客机:2小时跨太平洋

关键技术

本章小结

超声速流动设计的核心在于理解和控制激波。从最简单的楔形体和锥形体绕流出发,我们建立了以下关键概念:

基本原理

  1. 激波关系:θ-β-M关系决定了激波角度,自然界偏好弱激波解
  2. 波阻本质:熵增导致的不可逆能量损失,与物体形状和马赫数密切相关
  3. 面积定律:喉道处达到声速是产生超声速流的必要条件
  4. 流动匹配:进气道、发动机、喷管必须协调工作

设计准则

关键公式汇总

参数 公式 适用条件
波阻系数(楔) $C_D = 4\theta^2/\sqrt{M^2-1}$ 小偏转角
波阻系数(锥) $C_D = 2\theta_c^2/\sqrt{M^2-1}$ 细长锥
面积比 $A/A^* = f(M)$ 等熵流动
总压恢复 $p_{02}/p_{01} = f(M_1)$ 正激波

工程经验值

练习题

基础题

12.1 激波角计算 一个半楔角为15°的楔形体在Ma=2.5的气流中飞行。计算: a) 附体激波的激波角 b) 激波后的马赫数 c) 总压损失

提示 使用θ-β-M关系式或查表。记住自然界选择弱激波解。总压损失可用正激波关系式的斜激波分量计算。
答案 a) 激波角 β ≈ 45° b) 激波后马赫数 Ma₂ ≈ 1.87 c) 总压损失约 5% 激波后的马赫数仍为超声速,这是弱激波的特征。

12.2 喷管设计 设计一个将气流从Ma=0.3加速到Ma=3.0的拉瓦尔喷管。如果喉道直径为10cm,计算出口直径。

提示 使用面积-马赫数关系。注意入口处也需要用此关系式。面积比是出口面积与喉道面积之比。
答案 - 喉道面积:A* = π(5cm)² = 78.5 cm² - Ma=3.0时,A/A* = 2.64 - 出口面积:A_exit = 2.64 × 78.5 = 207 cm² - 出口直径:d_exit = 16.2 cm 入口处Ma=0.3,A/A* = 2.04,所以入口直径约14.3cm。

12.3 波阻估算 比较以下三种头部形状在Ma=2时的波阻系数: a) 半锥角20°的圆锥 b) 半楔角20°的楔形体 c) 半球形钝头体

提示 使用本章给出的波阻公式。注意角度要转换为弧度。钝头体使用经验公式。
答案 a) 圆锥:$C_D = 2(20°π/180°)²/\sqrt{4-1} = 0.14$ b) 楔形:$C_D = 4(20°π/180°)²/\sqrt{4-1} = 0.28$ c) 钝头:$C_D = 0.85 + 0.35/4 = 0.94$ 圆锥波阻最小,钝头体波阻最大,相差近7倍!

12.4 进气道起动 一个超声速进气道在Ma=2.0时捕获面积为1.0m²。如果要保证起动,喉道面积至少应该多大?

提示 使用Kantrowitz起动准则。查Ma=2时的A/A*值。实际设计需要留有裕度。
答案 Ma=2.0时,(A/A*) = 1.69 根据起动准则: $A_{throat}/A_{capture} > 1/1.69 = 0.59$ 因此 $A_{throat} > 0.59 m²$ 实际设计通常取0.65-0.70 m²以确保可靠起动。

挑战题

12.5 多级压缩设计 设计一个三级楔形压缩系统,将Ma=3的气流通过三道斜激波减速。要求:

提示 等转角设计不一定是最优的。可以先假设每级转角10°,然后迭代调整。考虑每道激波的压力损失累积。
答案 优化设计(非等角): - 第一级:θ₁ = 12°, Ma₁ = 3.0 → Ma₂ = 2.35 - 第二级:θ₂ = 10°, Ma₂ = 2.35 → Ma₃ = 1.85 - 第三级:θ₃ = 8°, Ma₃ = 1.85 → Ma₄ = 1.35 总转角:30° 总压恢复:约85% 相比单个30°楔形(会产生脱体激波),性能大幅提升。

12.6 喷管非设计工况分析 一个火箭喷管设计用于海平面(背压101kPa),出口压力50kPa,燃烧室压力5MPa。当火箭上升到10km高度(背压26kPa)时: a) 判断流动状态 b) 估算推力变化 c) 会出现什么现象?

提示 比较出口压力与环境压力。欠膨胀会形成膨胀扇,过膨胀可能导致激波和分离。使用推力公式考虑压力项。
答案 a) 海平面:过膨胀(50kPa < 101kPa),喷管内有激波 10km高度:欠膨胀(50kPa > 26kPa),出口有膨胀扇 b) 推力增加约15-20%(压力项贡献) c) 现象: - 海平面:可能出现流动分离,侧向载荷 - 高空:羽流扩张,但性能接近设计值 - 这解释了为什么第一级火箭常用海平面喷管

12.7 乘波体概念设计 为Ma=6巡航设计一个简单的乘波体。基准流场采用15°楔形流,估算: a) 理论最大升阻比 b) 如果在Ma=4飞行,性能如何变化? c) 主要的工程挑战是什么?

提示 使用L/D = 4(M+3)/(M+1)估算理论极限。非设计点性能下降显著。考虑热、结构、控制等实际问题。
答案 a) Ma=6时理论最大L/D = 4(6+3)/(6+1) = 5.14 考虑粘性,实际L/D ≈ 3.5-4.0 b) Ma=4时: - 激波不再完全附体于前缘 - 升力损失约30% - 阻力增加约20% - 实际L/D降至2.5-3.0 c) 工程挑战: - 前缘热流:~10 MW/m²,需要主动冷却 - 尖锐前缘制造困难且易损坏 - 横向稳定性差,需要主动控制 - 有效容积小,装载困难

12.8 综合设计问题 设计一个Ma=2.5巡航导弹的气动外形。要求射程最大化,考虑: a) 头部形状选择 b) 进气道类型 c) 弹体截面积分布 d) 估算升阻比和航程

提示 平衡波阻、容积、结构重量。考虑发动机集成。使用Breguet航程公式。
答案 优化设计方案: a) 头部:轴对称冯·卡门曲线(介于尖锐与钝头之间) - 波阻适中 CD ≈ 0.05 - 可容纳导引头 b) 进气道:下颌式混压进气道 - 两道斜激波 + 正激波 - 总压恢复 ~90% c) 面积分布:修正Sears-Haack体 - 考虑进气道凸起的补偿 - 长细比 L/D ≈ 12 d) 性能估算: - L/D ≈ 4.5(巡航状态) - 航程 R = (L/D) × (V/SFC) × ln(W₀/Wf) - 典型值:500-800 km 关键权衡:更尖锐→波阻小但容积小;更钝→容积大但阻力大。

常见陷阱与错误

1. 激波角度混淆

错误:将偏转角直接当作激波角

错误:15°楔形 → 15°激波角 ❌
正确:15°楔形 → 约45°激波角(Ma=2.5时)✓

记忆方法:激波角总是大于偏转角,且随马赫数增加而减小。

2. 强弱激波解选择

错误:随意选择强激波或弱激波解

正确理解

3. 喷管设计误区

错误:认为喷管越长效率越高

实际情况

错误:忽视非设计工况

正确做法

4. 进气道起动问题

错误:设计时只考虑巡航马赫数

问题

解决方案

5. 波阻计算错误

常见错误

a) 单位混淆

错误CD = 4θ²/(M²-1)  # θ用度数 ❌
正确CD = 4θ²/(M²-1)  # θ用弧度 ✓

b) 忽视干扰

c) 基准面积选择

6. 乘波体设计陷阱

理想vs现实

理论假设 实际情况 性能损失
无粘流动 边界层厚 L/D降低30%
尖锐前缘 钝化处理 升力损失10%
设计马赫数 变马赫数飞行 非设计点L/D降低50%
二维流动 三维效应 横向失稳

7. 热问题低估

错误:只考虑空气动力,忽视气动加热

驻点温度估算: \(T_0 = T_∞(1 + 0.2M²)\)

Ma=3时,驻点温度可达750K!

后果

8. 尺度效应忽视

风洞试验陷阱

修正方法

调试技巧

1. 快速判断激波类型

2. 性能快速评估

3. 实验观察要点

最佳实践检查清单

概念设计阶段

任务定义

初步布局

激波系统设计

外部压缩

进气道设计

推进系统集成

喷管设计

推力矢量(如适用)

性能评估

阻力分解

关键参数计算

非设计点考虑

变马赫数性能

变攻角效应

热管理

气动加热评估

结构温度限制

试验验证

风洞试验规划

测量技术

制造可行性

几何约束

成本考虑

安全性评估

失效模式

应急预案

优化迭代

性能提升潜力

验证方法

文档记录

设计依据

经验教训


记住:超声速设计是激波管理的艺术。每个设计决策都要考虑激波的影响,成功的设计在于找到性能、可靠性和可制造性之间的最佳平衡点。