参数设置
假设空气 (γ = 1.4, R = 287 J/(kg·K))·等熵·一维·定常流动。
收敛喷嘴的断面和流动
左=上游（P_0, T_0）／中央缩颈=喉部 A*／右=背压区域 P_b。颜色代表速度（蓝=低速·红=高速），红线表示窒息时的音速面。
质量流量 ṁ 与背压比 P_b/P_0
当 P_b/P_0 ≤ 0.528 时窒息（ṁ = ṁ_max）／以上时背压升高导致 ṁ 降低。黄色竖线=当前 P_b/P_0。
理论与主要公式
假设完全气体、等熵、准一维定常流动，收敛喷嘴的临界条件仅取决于 γ。喉部马赫数 M = 1 时的临界压力比：
$$\frac{P^*}{P_0} = \left(\frac{2}{\gamma+1}\right)^{\gamma/(\gamma-1)}$$
当 P_b/P_0 ≤ P*/P_0 时发生窒息，质量流量达到最大值并保持不变：
$$\dot{m}_{\max} = A^*\,P_0\,\sqrt{\frac{\gamma}{R\,T_0}}\,\left(\frac{2}{\gamma+1}\right)^{(\gamma+1)/(2(\gamma-1))}$$
窒息时的喉部温度和音速：
$$T^* = T_0\,\frac{2}{\gamma+1}, \qquad a^* = \sqrt{\gamma\,R\,T^*}$$
对于 γ=1.4 的空气，P*/P_0 ≈ 0.528，T*/T_0 ≈ 0.833。降低背压时，喉部音速面形成"信息屏障"，因此 ṁ_max 不再增加。
窒息流模拟器说明
🙋
压缩空气配管中听说"发生了窒息，所以流量不能再增加"，这是什么意思？
🎓
简单地说，就是喷嘴或节流部最狭窄处的气体速度达到了音速。看一下模拟器的初始值：P_0 = 500 kPa，P_b = 200 kPa。背压比是 0.4，低于临界压力比 0.528，所以质量流量 ṁ = 0.117 kg/s 达到了极限，无法继续增加。即使进一步降低 P_b，流量也不会增加一克。这种现象在实际工程中的安全阀、破裂片和小孔节流孔中很常见。
🙋
为什么流速达到音速流量就停止增加呢？听说与信息传播有关？
🎓
完全正确。背压的信息通过声波向上游传播，但当喉部流速等于音速时，下游的声波再也无法传回上游。也就是说，喉部以下的压力变化对喷嘴入口的流动完全没有影响。从准一维流方程 dA/A = (M²−1) dV/V 来看，当 M < 1 时，缩小断面会加速流动；但当 M = 1 时，断面变化必须为零。因此收敛喷嘴的喉部是音速的极限。如果要加速到超音速，必须使用拉瓦尔（收敛-扩张）喷嘴。
🎓
分界点是 P_b/P_0 = 0.528。对于默认的 P_0 = 500 kPa，当 P_b > 264 kPa 时窒息状态就消失，进入"亚音速流"。试试把 P_b 调到 270 kPa，你会看到流量从 0.117 开始逐渐下降。继续增加到 P_b = 500 kPa（等于 P_0），流量就变成 0，流动停止。点击播放按钮扫描 P_b，你能清楚地看到这条 S 形曲线。
🙋
显示喉部速度是 317 m/s，这就是音速吗？
🎓
没错。在 T_0 = 300 K 的情况下，喉部音速 a* = √(γ R T*) 中，T* = 2T_0/(γ+1) = 250 K。所以 317 m/s 低于地面常温音速 343 m/s，原因是断热膨胀使喉部温度下降。试试用滑块增加 T_0，你会看到 a* 随之变大。比如 T_0 = 600 K 时，a* ≈ 448 m/s。
🎓
看公式 ṁ_max = A* P_0 √(γ/(R T_0)) (2/(γ+1))^((γ+1)/(2(γ-1)))，在窒息状态下只有三种办法：(1) 增加上游全压 P_0，(2) 放大喉部面积 A*，(3) 降低上游全温度 T_0（相同温度下密度更高）。背压 P_b 已经失效了。实际工程中常见做法是充分利用设计压力 P_0，然后根据所需流量 ṁ 反推出 A*。
常见问题解答
在完全气体、等熵流动的假设下，当喉部马赫数 M = 1 时，用滞止点和喉部的关系式 P_0/P* = (1 + (γ−1)/2)^(γ/(γ-1))，可推导出 P*/P_0 = (2/(γ+1))^(γ/(γ-1))。对空气 γ = 1.4，得 (2/2.4)^3.5 = 0.5283。对于 He（γ = 1.67）是 0.487，CO_2（γ = 1.30）是 0.546，不同气体略有差异。本模拟器使用的是空气的值。
本模拟器基于完全气体、等熵、准一维定常流动的假设。实际的孔板或控制阀需要乘以流出系数 C_d（通常 0.6～0.95）：ṁ_实际 = C_d · ṁ_理想。对于低温高压液化气（如 LNG），蒸发和焦耳-汤姆逊效应不可忽视，需要调用 NIST REFPROP 等状态方程库进行数值积分。API 520、ISO 4126 等安全阀设计标准也规定了包含 C_d 的实用公式。
拉瓦尔喷嘴的喉部临界条件和 ṁ_max 公式保持不变。区别在于：喉部下游有扩张部分，在设计压力比下，流体可以加速到超音速。若背压高于设计值，扩张部会产生垂直激波，下游恢复为亚音速。但只要发生窒息，喉部的 ṁ 就相同，所以本模拟器的"ṁ 与 P_b/P_0"曲线在窒息区（左半部分）对收敛喷嘴和收敛-扩张喷嘴的行为是完全相同的。
因为放出流量不依赖下游条件，仅由容器压力 P_0 和喉部面积 A* 决定，设计简单且重现性强。下游无论连接什么——配管、消音器、洗气瓶等——放出能力都有保证。API 520 和 ASME Section VIII 的安全阀设计公式正是以窒息状态为前提计算 ṁ_max 的。防止容器在火灾或反应暴走时爆裂的核心做法，就是根据最大需排流量反推设计喷嘴面积。
实际应用
压力容器的安全阀与破裂片：化工厂、液化天然气罐、空气瓶等在内压异常上升时需要立即泄压。核心是"始终保持窒息状态设计"——背压变化（下游配管变化）不影响泄压能力，容器压力能被迅速降低。API 520 和 ISO 4126 标准的 ṁ_max 公式就是本模拟器的公式乘以流出系数 C_d，从所需流量反向计算需要的面积 A*。比如用模拟器设置 P_0 = 1000 kPa，A* = 100 mm²，一个安全阀能排出约 0.23 kg/s 的空气。
孔板流量计（临界流喷嘴）：广泛用作气体流量标准器的 CFV（临界流文丘里）喷嘴正是为了窒息而设计的。ISO 9300 规范了这种方法，用 ṁ = C_d · ṁ_max 的简单关系可达到 0.5% 以下的精度。下游连什么都不影响流量，所以可以同时校准多个流量计，是计量领域的利器。
火箭和喷气发动机喷嘴：固体火箭、液体火箭、涡轮喷气发动机的燃烧气必然经过窒息状态。燃烧室压力（相当于 P_0）确定，喉部 ṁ_max 就唯一确定，推力 F = ṁ · V_e + (P_e − P_atm)·A_e 的设计变得直观。土星 V 火箭的 F-1 发动机在燃烧室压力 7 MPa、喉部直径 0.95 m 条件下实现 ṁ ≈ 2.6 t/s 的巨大窒息流量。用模拟器将单位改成 MPa 和 m²，能直观体感这个数量级。
真空系统泄漏评估：真空腔室或半导体工艺设备中"真空度升不上来"的原因排查时，常用窒息流模型来评估漏孔。室温大气压空气从 0.1 mm 直径孔泄漏，A* ≈ 0.008 mm² 对应 ṁ ≈ 9 μg/s。从设计压力和面积推算最坏泄漏量，反推所需排气速度的计算中，本模拟器的公式被直接应用。这在半导体、航天、粒子物理的真空配管设计中是常见课题。
常见误解与注意事项
最普遍的误解是"背压越低流量越大"的直觉。但在窒息区，背压的信息被音速壁垒阻挡，无法传回上游。模拟器中从 P_b = 200 kPa 降到 10 kPa，ṁ 仍然固定在 0.117 kg/s，纹丝不动。"真空吸就能增大流量"这个想法仅在 P_b/P_0 > 0.528 的亚音速区有效，把这个直觉带进设计会导致容量不足。
第二个常见错误是"ṁ_max 只由喉部状态决定"的理解偏差。公式 ṁ_max = A* P_0 √(γ/(R T_0)) (2/(γ+1))^((γ+1)/(2(γ-1))) 中，出现的是上游的 P_0、T_0 和喉部面积 A* 与气体常数，喉部压力 P* 和温度 T* 在中途被消去。实现代码时别想着"用喉部压力除"或"单独测喉部温度"，那样容易出错。模拟器直观演示了：仅改变 P_0 和 T_0，ṁ_max 就直接变化，而不需要计算中间状态。
还有一个误解是"收敛喷嘴也能加速到超音速"。从 dA/A = (M²−1) dV/V 看，M < 1 时缩小断面升速，但 M = 1 处必然 dA = 0，要使 M > 1 必须反向加大断面。收敛喷嘴的喉部音速是物理上的极限。无论背压多低，收敛喷嘴出口速度最高就是音速。要产生超音速气流，必须使用拉瓦尔（收敛-扩张）形状。