参数设置
假设完全气体、绝热、一维定常流动。f 为达西摩擦系数（与穆迪图相同的定义）。当实际 4fL/D 超过入口处的 4fL*/D 时流动壅塞，出口马赫数被锁定在 M=1。
管道示意图
左＝状态 1（入口）／右＝状态 2（出口）／箭头长度＝局部速度／壁面波线＝摩擦／颜色渐变指示相对温度
法诺曲线 (T-s 图)
纵轴＝T/T*／横轴＝(s − s*)/c_p／上分支＝亚音速，下分支＝超音速／蓝点＝状态 1，红点＝状态 2／橙线＝1 → 2 路径
理论与主要公式
法诺流动是绝热、带壁面摩擦、等截面管道中可压缩流动的理想模型。摩擦总是把流动推向 M=1：亚音速分支被加速，超音速分支被减速。
由当前马赫数 M 到 M=1 的无量纲壅塞距离：
$$\frac{4 f L^*}{D_h} = \frac{1 - M^2}{\gamma M^2} + \frac{\gamma+1}{2\gamma}\ln\!\left(\frac{(\gamma+1)M^2}{2 + (\gamma-1)M^2}\right)$$
出口马赫数 M_2 由管道长度 L 通过下式数值求解：
$$\frac{4 f L^*_2}{D_h} = \frac{4 f L^*_1}{D_h} - \frac{4 f L}{D_h}$$
温度比与压力比以声速参考态表示：
$$\frac{T}{T^*} = \frac{\gamma+1}{2 + (\gamma-1)M^2}, \qquad \frac{P}{P^*} = \frac{1}{M}\sqrt{\frac{\gamma+1}{2 + (\gamma-1)M^2}}$$
T_2/T_1 与 P_2/P_1 即上述参考形式之比。当实际 4fL/D 超过 4fL*_1/D 时流动壅塞，出口被锁定在 M=1。
法诺流动模拟器是什么
🙋
老听人说"天然气管道越长，质量流量就越上不去"，这背后到底是什么物理？
🎓
这就是法诺流动。在绝热、等截面、带壁面摩擦的管道里，摩擦总是把流动推向 M=1。亚音速入口越往下游越快，足够长时出口就变成 M=1，再延长也没用，这叫'壅塞 choke'。用模拟器默认值（M_1 = 0.30，γ = 1.4，f = 0.020，L/D = 50）跑一下，出口大约 M_2 ≈ 0.474。
🙋
从 M_1 = 0.30 只加速到 M_2 = 0.47 啊。继续加大 L/D 真的能让出口逼近 1 吗？
🎓
对。L/D 一大，4fL/D 就向入口的 4fL*/D（这里大约 5.30）靠近。当 L/D ≈ 66 时出口正好 M=1，管道就壅塞了。再想送更长，就只能降低入口压力或入口马赫数重新设计。把 L/D 滑块拉到很大，可以看到 M_2 卡在 1.000，并出现红色 CHOKED 徽标。
🙋
右边 T-s 图里 M=1 在最右边，上下还分成两支，这是什么意思？
🎓
那条曲线就是法诺曲线。纵轴是 T/T*，横轴是无量纲熵；上分支是亚音速、下分支是超音速，两支在 M=1 处汇合，正好是熵的最大点。绝热条件下熵只能增加，所以状态点只能沿曲线向 M=1 方向走。仔细看，蓝色入口点经橙色虚线移到红色出口点，就是沿着亚音速分支朝右下走的。
🙋
如果入口换成超音速，比如 M_1 = 2.0，行为会反过来吗？
🎓
试试就知道。这次摩擦让 M_1 = 2.0 减速向 M=1 靠近，T-s 图上状态点沿着超音速下分支，从右下出发往左上走。看上去和亚音速的方向相反，但'熵升高、状态趋向 M=1'这个本质完全相同，这正是法诺流动最迷人的对称性。
常见问题
本工具的 f 是'达西摩擦系数（Darcy friction factor）'，即穆迪图直接读出的值。部分教材使用'Fanning 摩擦系数 f_F'，二者关系是 4·f_F = f_D。本工具的公式 4fL*/D 默认采用达西定义，将 4·f_D·L/D 直接代入即可。湍流光滑管的典型值 f_D ≈ 0.02，工业粗糙管约 0.04 ~ 0.06。
出口达到 M=1 之后，该截面的质量流量被锁定在'上游条件允许的最大值'，再降低背压也无法让出口越过 1，所以流量保持不变。这与收缩喷嘴的临界流是同一物理。实际工程会通过加大管径 D 或降低入口马赫数来事先避免壅塞，并以理论值为上限取一定安全裕度。
由于法诺流动是绝热的，总焓 h_0 = h + V²/2 守恒。亚音速向下游加速时速度 V 增大，静温 T 必然下降。默认参数下下降很小（T_2/T_1 ≈ 0.974），但若入口马赫数很低、出口逼近 M=1，T_2/T_1 可降到约 0.83。在超音速分支上则相反：减速使 T 上升。
法诺流动是'绝热 + 带摩擦'的理想模型；瑞利流动是'无摩擦 + 带热量交换'的理想模型。两者都适用于等截面可压缩管道流动，且都把流动推向 M=1，但 T-s 图上的曲线形状与熵的演化不同。实际燃烧室与换热器同时含两种效应，工程上按主导效应近似选用其中之一。
实际工程应用
长输管道设计：横跨数十至数百公里的天然气、压缩空气长输管道，必须用法诺流动来确定管径。从所需的出口流量与入口压力反推不会壅塞的最小管径，并据此判断是否需要中途增压站。即便严格的绝热假设需要补充与外界的换热修正，法诺流动仍是计算的起点。
火箭与喷气发动机管道：液体火箭的推进剂供给管路、燃气轮机的引气冷却管、喷气发动机的中冷器管道等，都用法诺流动评估给定长度下出口马赫数能上升到何处。冷却引气管一旦壅塞，下游冷却空气将不足，可能引发严重事故，因此设计上需在所有工况下保持 4fL*/D 显著大于实际 4fL/D。
超音速风洞扩压器：超音速风洞试验段后接扩压器，把超音速气流减速回收。摩擦的作用类似法诺流动，扩压器长度与出口马赫数受壅塞极限约束。过长会壅塞、破坏试验条件；过短又增加下游真空泵容量。这种取舍是典型的法诺流动设计问题。
HVAC 与高速通风管道：建筑高速通风管道（马赫 0.1 ~ 0.3）严格做压降计算时也包含法诺流动修正。低马赫管道用不可压缩的 Darcy-Weisbach 公式即可，但工厂的高速排风、厨房抽油烟管等可能需要可压缩法诺分析，否则风机选型容易偏小导致风量不足。
常见误区与注意事项
最常见的误区是认为"摩擦总是让流动减速"。但在亚音速法诺流动中，摩擦实际上让气体'加速'。这看起来违反直觉，却直接来自'绝热假设 + 质量守恒'的组合。模拟器默认情况（亚音速 M_1 = 0.30 加速至 M_2 = 0.474）就直接证明了这一点。只有超音速分支才是减速的。
第二个误区是认为'管道想加多长就加多长'。法诺流动有 4fL*_1/D 这个硬性上限，越过就壅塞，出口被锁在 M=1。把模拟器的 L/D 滑块拉大，会看到 CHOKED 徽标出现。工程上要先算这个极限再取 1.3 ~ 2 的安全系数。
第三点要小心达西（Darcy）和 Fanning 摩擦系数的区别。本工具采用穆迪图的达西定义；使用 Fanning 系数 f_F 的教材满足 4·f_F = f_D，混淆会让压降计算偏小 4 倍。引用文献或教材数据时务必先确认采用的是哪一种定义。