5.3 拉瓦尔喷管与一维流动

5.3 拉瓦尔喷管与一维流动

本节摘要：拉瓦尔喷管是先缩后张的管道，亚音速气流在收缩段加速、在喉道达到声速、在扩张段膨胀为超音速。它能否在设计状态工作，取决于喉道壅塞与出口背压的配合。本节推导面积比公式，用 Python 求解双解马赫数，并把火箭发动机与蒸汽涡轮的现场数据带进来复算。

火箭尾焰那道明亮的锥形喷流，穿过一个精心计算的几何形状：拉瓦尔喷管。1888 年瑞典工程师拉瓦尔把它用在蒸汽涡轮上，半个世纪后成为所有火箭发动机的标准件。它的原理浓缩了上一节全部结论，是可压缩流动最漂亮的一件工程作品。

一、面积比公式：几何决定马赫数

等熵流动中，任一截面面积 A 与喉道面积 A* 之比只取决于当地马赫数（γ 取 1.4）：

A/A* = (1/Ma) 乘以 [(1 + 0.2 Ma²)/1.2] 的三次方。

注意这个方程一个面积比对应两个马赫数（一亚一超）——同一个几何形状，喷管可以是亚音速扩压器也可以是超音速喷管，几何不决定流态，边界条件（背压）才决定。这是拉瓦尔喷管一切"怪脾气"的根源。

import numpy as np
from scipy.optimize import brentq

def area_ratio(Ma, g=1.4):
    return (1/Ma) * ((1 + (g-1)/2*Ma**2) / ((g+1)/2))**((g+1)/(2*(g-1)))

def mach_from_area(AR):
    sub = brentq(lambda m: area_ratio(m) - AR, 1e-3, 1.0)
    sup = brentq(lambda m: area_ratio(m) - AR, 1.0, 30.0)
    return sub, sup

# 火箭喷管 出口与喉道面积比 8
sub, sup = mach_from_area(8.0)
print(f"面积比 8 : 亚音速解 Ma={sub:.3f}  超音速解 Ma={sup:.2f}")

# 喉道直径 80 mm 求出口直径
print(f"出口直径约 {80*np.sqrt(8):.0f} mm")
​

面积比 8 对应超音速解约 Ma 2.94（这正是 5.1 节算过的"总压 36.7 倍静压"工况）。设计一个喷管的流程就是反过来：想要多高的出口马赫数，查面积比定几何。

二、壅塞与工况图谱：背压说了算

喉道一旦达到声速，质量流量就到达最大值——壅塞。此后无论背压再怎么降，喉道上游毫无感觉（扰动过不了声速喉），流量不再增加。下游则随背压 p_b 与出口设计压 p_e 的关系演出一整幕剧：

背压从设计值升高时，一道正激波从出口向喉道推进：背压越高，激波位置越靠上游（当地马赫数低、激波弱），直到把它推回喉道、流动退化为亚音速。5.2 节的激波计算器在这里就是位置解算器。

现场算例：小型液体火箭发动机推力室

燃烧室总压 2.0 MPa、总温 3000 K、喉道直径 60 mm、面积比 4（γ 取 1.2 近似燃气）。估算海平面推力：

import numpy as np

g, R, T0, p0 = 1.2, 380.0, 3000.0, 2.0e6
At = np.pi * 0.06**2 / 4

# 壅塞流量公式 特征速度法
Gamma = np.sqrt(g) * (2/(g+1))**((g+1)/(2*(g-1)))
c_star = np.sqrt(R*T0) / Gamma
mdot = p0 * At / c_star
print(f"特征速度 {c_star:.0f} m/s  质量流量 {mdot:.2f} kg/s")

# 面积比4的推力系数 Cf 约 1.4
Cf = 1.40
F = Cf * p0 * At
print(f"海平面推力约 {F/1e3:.1f} kN")
​

质量流量约 2.7 kg/s、推力约 4 kN——一台小型上面级发动机的量级。喷管设计点应选在主要工作高度的环境压力：海平面发射的发动机面积比不敢做大（出口压过低会被大气压瘪、诱发流动分离与侧向载荷），真空级则可以做到面积比上百。

三、其他现场的拉瓦尔原理

• 蒸汽涡轮：喷嘴把高压蒸汽膨胀为高速射流打向动叶片，拉瓦尔的原始应用；
• 超音速风洞：拉瓦尔喷管建立试验段超音速流，下游再用第二喉道与扩压器回收压力；
• 超声速进气道：反着用——超音速来流经激波系减速增压到亚音速进入压气机，几何同样遵循面积-马赫数关系，只是方向倒过来；
• 日常彩蛋：拧松气球口让空气冲出、高压锅限压阀排气，出口射流的收缩-膨胀结构里都有拉瓦尔喷管的影子，只是雷诺数与马赫数都小得多。

💡 关键直觉：喉道是"信息单向阀"。声速喉道隔绝上下游，这就是壅塞的物理本质——上游对下游的一切变化免疫。设计变工况部件（可调喷管、可调进气道）本质上都是在挪动这道阀门。

四、过膨胀与欠膨胀：出口波系的读法

设计点之外的两种偏离都写在天上的羽流形状里。过膨胀（出口压低于背压，常见于海平面试车的真空设计喷管）：出口外先被压缩、出现斜激波与马赫盘，羽流边界先内收再外扩，照片上是一串亮珠似的马赫盘。欠膨胀（出口压高于背压，真空级发动机在海平面）：气流出口外继续膨胀，羽流鼓成胖梭形。读羽流形状反推出口压比，是火箭试验台的日常判读技能，也是判断喷管是否分离、侧向载荷风险的依据。

工程处理：海平面试车的高面积比喷管常采用双钟形或延伸喷管，起飞段用小面积比、高空段展开大面积比——用可动几何同时伺候两个背压世界，呼应了本节"边界条件决定流态"的主线。

本节要点回顾

• 面积比公式连接几何与马赫数，一个面积比对应亚/超两个解，流态由背压决定；
• 壅塞：喉道达声速后流量封顶，下游降压不再影响上游；
• 背压图谱：从全亚音速到管内激波、设计点、过膨胀，一张表覆盖全部工况；
• 推力估算用流量公式加推力系数，喷管面积比按主要工作高度选取；
• 拉瓦尔几何在涡轮、风洞、进气道中反复变装出现，识别"先缩后张加声速喉道"即可认出它。

高速世界到此收束。下一章回到低速，贴着壁面看粘性最精彩的表演：边界层与湍流。