🙋
SpaceX Merlin 的"针栓喷射器"和普通火箭发动机喷射器有什么区别？从照片上看，中间似乎立着一根粗针。
🎓
没错，那根"中间的针"就是 pintle（针栓），LOX 从这里沿轴向喷出。外侧套筒的燃料 RP-1 沿径向喷出，正好以90°碰撞。这就形成对撞喷雾，展开成锥形。普通同轴喷射器要并排数百个喷嘴，但针栓只需一根就能承载全部推力。这种设计在1969年的阿波罗月着陆船发动机 LMDE 中得到验证，SpaceX 在 Merlin 和 Dragon 的 SuperDraco 中复兴了它。
🙋
就一根针来承担全部功能，不会很危险吗？常理来看多个分散的设计更安全。
🎓
虽然直觉上这样想，但针栓设计与燃烧室的声学模式结合弱，高频燃烧不稳定很难出现。这与 F-1（土星五号）为了解决声学问题而采用挡板设计的痛苦不同。而且针栓能轻易实现 10:1 的节流比。LMDE 能达到 12.5:1，正因为如此才能在着陆时精确控制推力。Falcon 9 一级能垂直着陆，也是因为 Merlin 允许深度节流。简而言之，针栓设计一次性实现了"强、快、便宜"。
🙋
明白了！设计中最关键的参数是什么？TMR 经常被提到…
🎓
TMR (Total Momentum Ratio) = ṁ_f V_f / ṁ_o V_o 最关键。它表示"燃料和氧化剂的动量平衡"，直接通过 Cheng 公式 θ = 2·90°/(1+1/√TMR) 决定喷雾锥角。TMR 小时（Merlin 的 0.3～0.4）喷雾窄，TMR=1 时约90°，TMR=4 时约120°。Hutt & Cramer (1996) 的稳定性基准认为 0.7～2.5 是安全域，但 Merlin 故意在较低的 TMR 运行，配合 O/F=2.34 的偏氧化剂设置，目的是降低燃烧室热负荷。
🎓
SMD（Sauter平均粒径）是液滴的代表尺度，蒸发时间 ∝ SMD² 与燃烧室长度直接关联。越小蒸发越快，可以缩短 L*，但太小会激发声学不稳定。Merlin 级的目标是10～50 μm。本工具使用 SMD ∝ d_p·We⁻⁰·⁴·TMR⁻⁰·² 的 Mehegan 系列相关式近似。左边滑块调节针栓直径，可以看到 SMD 如何变化，这有助于获得设计直觉。
🙋
最后一个问题，Raptor（Starship）也是这种针栓喷射器吗？
🎓
Raptor 采用全流分级燃烧循环 + LOX/CH4，喷射元件据说是多元同轴旋流型，不是针栓。但 Merlin 在节流能力和燃烧稳定性上积累的经验也用在了 Raptor 的燃烧室设计上。本工具可以切换到 Raptor 的典型值（推力2300 kN、Pc=300 bar、O/F=3.6），能看出 TMR 和锥角如何变化，体会推进剂选择对喷射器设计的影响。
针栓喷射器通过中心轴的针栓和外周套筒以90°角进行推进剂对撞喷射，具有三个优点：(1) 易于实现 10:1 以上的节流比，可在着陆时进行大幅推力调节；(2) 仅需一个喷射元件，结构极简、成本低廉；(3) 与燃烧室声学模式结合弱，不易出现高频不稳定。阿波罗 LMDE（12.5:1 节流）曾验证该方案，SpaceX 的 Merlin（845 kN，LOX/RP-1）和 Raptor（LOX/CH4）也采用了类似布局。
TMR = (ṁ_f V_f) / (ṁ_o V_o) 是针栓设计的最关键参数，表示燃料侧和氧化剂侧的动量比。Hutt-Cramer 稳定性准则认为 TMR 0.7～2.5 为稳定域，0.3 以下或 4 以上属于不稳定风险区。实际中 Merlin 运行在 0.3～0.4 范围内，与 O/F=2.34 的偏氧化剂设计相符。TMR 过低会导致燃料被推向中心，混合恶化；过高则导致过多壁面喷溅，热负荷增加。
锥体全角 θ 由 Cheng 等人(2017)的经验式 θ = 2·90°/(1 + 1/√TMR) 近似，TMR 越高喷雾越发散。TMR=1 时 θ≈90°，TMR=4 时 θ≈120°，TMR=0.25 时 θ≈60°。锥角决定了对燃烧室壁面的热负荷和循环流模式，应与燃烧室 L*（特征长度）一起设计。过宽导致壁面烧损，过窄则中心留下未燃区。
SMD 是液滴的体积/表面积比，代表蒸发和燃烧速率。本工具使用 SMD ≈ d_pintle·50·We⁻⁰·⁴·TMR⁻⁰·² 的简化相关式。SMD 越小蒸发越快，可缩短燃烧室，但过小会激发声学不稳定。Merlin 级别以 10～50 μm 为目标，超过此范围会出现未燃损失（c* 效率下降）或高频不稳定风险。
SpaceX Merlin（Falcon 9 一级）：海平面推力 845 kN、Pc=95 bar、LOX/RP-1、O/F=2.34 的燃气发生器循环。单级 9 台并联运行，单个针栓喷嘴的高结构裕度和 40～100% 的深度节流范围使得 Falcon 9 一级能够垂直着陆。本工具的默认值再现了 Merlin 1D 的特性。
Apollo LMDE（月着陆船下降发动机）：1969 年成功实现了人类首次月球着陆的 45 kN 级针栓发动机（N2O4/Aerozine 50）。在当时首次实现了 12.5:1 的节流，使宇航员能在着陆最后阶段手动调节推力。本工具的 N2O4/MMH 预设 + 低推力可以看到相当于 LMDE 的运行区域的 TMR、锥角。
SpaceX Dragon SuperDraco（紧急脱离发动机）：8 台 71 kN N2O4/MMH 针栓发动机组，为 Crew Dragon 提供发射时紧急脱离（LES）和未来动力着陆能力。100% 到 20% 的瞬时节流得以实现，正是因为针栓的特性。
CAE / CFD 前期分析：进行详细反应性 LES（大涡模拟）之前，用本工具这样的相关式快速估算 TMR、锥角、SMD，结合燃烧室 L*（典型 0.8～1.5 m），缩小设计空间。反之，若 CFD 结果与此估算相差很大，应检查边界条件或雾化模型设置。
最大的陷阱是认为"提高 TMR 必然提高燃烧效率"。确实 TMR 升高会拓宽锥角、改善混合，但过宽的锥会让未燃液滴直接打在燃烧室壁上，导致热负荷和壁面烧损。Merlin 之所以在 TMR 0.3～0.4 这样较低的值运行，是因为 O/F 2.34 的偏氧化剂设计使燃烧温度控制在 3670 K，这个温度范围内再生冷却就能管理。"TMR 稳定域 0.7～2.5"的教科书规则不能盲目套用，要理解实机为何要突破这个范围。
其次，别以为"比推力 Isp 只由推进剂组合决定"。本工具简化假设 Isp ≈ 300 s 量级（c_star × η_c* = 1.5×3000 m/s），但实际设计中 c* 效率在 92～98% 之间波动，直接影响 Isp。TMR 降至 0.5 以下时 c* 效率显著下降，推力、Isp 会有个位数百分比的衰减。工具单独输出"混合效率"正是这个原因——有时牺牲 c* 效率来降低壁面热负荷，这个权衡由再生冷却能力决定。
最后，小心"SMD 越小越好"的过度简化。若把 SMD 压至 5 μm 以下，蒸发虽快，但液滴群的声学响应会与燃烧室纵向或切向模态耦合，引发"高频不稳定（screech）"。Apollo F-1 正是在这个区间吃过亏。实机通常把 SMD 控制在 10～50 μm，必要时配合挡板或声学阻尼器额外衰减。"细化雾化"的决策必须与燃烧室声学模态分析（Helmholtz/纵向）联动。