参数设置
假设完全气体、绝热、一维、稳态流。f 采用穆迪图相同的达西摩擦系数定义。当实际管道的 4fL/D 超过入口侧的 4fL*/D 时,发生堵塞现象,出口马赫固定为 M=1。
管道示意图
左=状态 1(入口)/右=状态 2(出口)/箭头长度=速度/下方波线=壁面摩擦/温度、压力用颜色表示相对值
范诺线图 (T-s)
纵轴=T/T*/横轴=无量纲熵 (s−s*)/cp/上分支=亚音速,下分支=超音速/蓝圆=状态 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 通过以下关系的数值求解获得:
$$\frac{4 f L^*_2}{D_h} = \frac{4 f L^*_1}{D_h} - \frac{4 f L}{D_h}$$
温度比和压力比用以 M=1 为基准的参考量表示:
$$\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 分别由 T/T* 和 P/P* 的比值获得。当实际管道的 4fL/D 超过 4fL*_1/D 时发生堵塞,出口固定在 M=1。
范诺流模拟器简介
🙋
我听说天然气管道"越长流量越受限",这是什么物理原理呢?
🎓
这正是范诺流的体现。在绝热、有摩擦、恒定截面的管道中,摩擦作用总是使流动趋向 M=1。亚音速流体进入时会沿程加速,管道足够长时出口达到 M=1,此后流动发生『堵塞(choking)』,无法再加速。默认参数(M_1=0.30, γ=1.4, f=0.020, L/D=50)下,出口应显示 M_2 ≈ 0.474。
🙋
从 M_1=0.3 只加速到 M_2=0.47,并没加速多少呢。把 L/D 增加的话能继续加速吗?
🎓
可以的。增加 L/D,当 4fL/D 等于 4fL*_1/D(现在约为 5.30)时,出口会恰好堵塞在 M=1。在目前参数下,当 L/D ≈ 66 时发生堵塞。再想增加长度的话,只能降低上游压力或减小入口马赫重新设计。你可以用模拟器调节 L/D 滑块,看 M_2 逼近 1 时『堵塞』标签的出现。
🙋
右边 T-s 线图上,M=1 的点在右端,上下各有一条分支,这是什么?
🎓
那就是范诺曲线。纵轴是 T/T*,横轴是无量纲熵,上分支是亚音速、下分支是超音速。M=1 时达到最大熵,这是范诺曲线的特征。在绝热摩擦条件下,熵必然增加,所以流动沿曲线只能向 M=1 的方向(向右)移动。现在的入口(蓝点)到出口(红点)的箭头,应该沿着亚音速分支向右下方移动。
🙋
那如果我把入口改成超音速,比如 M_1=2.0,会怎样?
🎓
用滑块调节 M_1=2.0 试试看。这次摩擦会让流动『减速』,M_2 会从 2.0 向 1 下降。在 T-s 线图上,入口会在下分支(超音速)的右下方,然后沿下分支向左上方(M=1)移动。看起来与亚音速的行为完全相反,但本质上都是『趋向 M=1』,这就是范诺流最有趣的地方。
常见问题
本工具中的 f 是『达西摩擦系数(Darcy friction factor)』,直接输入穆迪图上显示的数值即可。某些教科书使用『范宁摩擦系数(Fanning friction factor)』,两者关系为 4f_F = f_D,请注意区分。本工具的 4fL*/D 公式以达西 f 为基础,直接使用 4·f_D·L/D。乱流光滑管典型值 f_D ≈ 0.02,粗糙管约 0.04~0.06。
出口达到 M=1 时,该出口截面的质量流量就达到了『在上游条件下允许的最大值』。无论如何降低下游压力,出口马赫都突破不了 1,因此流量保持不变。这与收敛喷嘴的『临界流』物理机制完全相同。实际工程中提前计算堵塞极限,通过加大管径或降低入口马赫来避免堵塞,设计时通常留 1.3~2.0 倍的安全系数。
范诺流是绝热的,所以总焓 h_0 = h + V²/2 得以保存。亚音速加速时,速度 V 增加,静温 T 必然下降以维持总焓守恒。默认参数下 T_2/T_1 ≈ 0.974,只下降少许。但若入口马赫很小而出口接近 M=1,温度比可能降至 0.83 左右。相反,超音速减速时温度会上升。
范诺流是『绝热 + 有摩擦』模型,瑞利流是『无摩擦 + 有热交换』模型。两者都是恒定截面管道可压缩流的理想模型,都会使流动趋向 M=1,但熵变和 T-s 线图形状不同。实际燃烧室和热交换器中两种效应并存,因此需要灵活选用作为近似模型。
实际应用
长距离管道设计:天然气或压缩空气运输数十至数百公里的长距离管道设计时,根据下游所需流量和上游压力反推,选择不会发生堵塞的管径。范诺分析确定『该压力下可以运输多长的距离』以及『需要多长时再加装压缩机』,是设计初期的关键计算。实际工程中还需加入与环境热交换的修正。
火箭和喷气发动机的管道:液体火箭推进剂供应管线、燃气轮机的引气冷却管、喷气发动机的涡轮冷却气管路等,都需用范诺流评估给定管道长度下出口马赫能达到多少。特别是引气管若发生堵塞,冷却气供应不足会导致严重事故,设计必须在 4fL*/D 基础上留足余量。
超音速风洞扩压段:超音速风洞在测试段后接扩压器以减速并回收气体。摩擦效应会导致范诺流式的减速,扩压段长度和出口马赫的关系需提前预测以避免堵塞。太长易堵塞影响试验条件,太短则下游真空泵需要更大容量,存在设计权衡。
建筑换气管道:高速建筑排气管道(马赫 0.1~0.3 量级)也会产生范诺型压力损失。低马赫下用非压缩达西-韦斯巴赫公式足够,但高速厨房排风或工业排气管可能需要可压缩性修正。压损估算偏低会导致送风机规格不足。
常见误区与注意事项
最常见的误解是认为『有摩擦的流动必然减速』。而事实是亚音速范诺流中,摩擦反而『加速』流动。这违反直觉,但结合绝热条件和质量守恒就能得出这个结论。模拟器默认参数(亚音速 M_1=0.30 → M_2=0.474 加速)就能清楚演示。减速只在超音速分支发生。
另一常见错误是『无限延长管道就能无限加速』。范诺流有上限 4fL*_1/D,超过它就发生堵塞,出口固定在 M=1。用模拟器拉动 L/D 滑块上升,看『堵塞』标签出现就能印证。实际设计需提前计算该上限,并预留 1.3~2 倍安全系数。
还要注意达西(Darcy)和范宁(Fanning)摩擦系数的区别。本工具用的是达西系数,穆迪图上的数值直接使用。某些教科书用范宁系数,两者满足 4f_F = f_D,搞混会导致压损计算差 4 倍。查阅论文或教材时务必确认定义。