范诺流（摩擦管道流）模拟器 — 摩擦壅塞与出口马赫数计算

参数设置

入口马赫数 M₁

仅限亚音速流动（0.05～0.90）

管径 D

管长 L

达西摩擦系数 f

取自穆迪图。商用钢管约 0.015～0.030

比热比 γ

空气＝1.40，天然气≈1.31，CO₂≈1.30

计算结果

—

入口 fL*/D

—

可用 fL/D

—

出口 fL*/D

—

出口马赫数 M₂

—

压力比 P₂/P₁

—

壅塞判定

—

摩擦管道流 — 动画

亚音速气体进入管道后因壁面摩擦而加速，马赫数沿程上升。当出现壅塞时，管道中部会出现 M=1 的壁。

范诺函数 fL*/D 与马赫数

出口马赫数与管长

理论与主要公式

$$\frac{fL^{*}}{D}=\frac{1-M^{2}}{\gamma M^{2}}+\frac{\gamma+1}{2\gamma}\ln\!\frac{(\gamma+1)M^{2}}{2+(\gamma-1)M^{2}}$$

范诺函数 fL*/D 是从马赫数 M 通过摩擦到达 M=1（壅塞）所需的无量纲摩擦长度。亚音速时摩擦使 M 上升至 1；超音速时摩擦使 M 下降至 1。

$$\frac{fL_{1\to 2}}{D}=\left(\frac{fL^{*}}{D}\right)_{\!M_{1}}-\left(\frac{fL^{*}}{D}\right)_{\!M_{2}}$$

入口到出口的实际摩擦长度 fL/D 等于入口与出口范诺函数值之差。当出口值降至 0 时，发生摩擦壅塞。

$$\frac{P}{P^{*}}=\frac{1}{M}\sqrt{\frac{\gamma+1}{2+(\gamma-1)M^{2}}}$$

范诺线上的静压比 P/P*。出口/入口压力比为 P₂/P₁＝(P/P*)|_{M₂} / (P/P*)|_{M₁}。亚音速流动中静压沿程降低。

什么是范诺流（摩擦管道流）

🙋

"范诺流"是计算长管道中气体的吗？为什么不能直接用伯努利方程？

🎓

好问题。伯努利方程的前提是"不可压缩、无摩擦"，对于短水管这没问题。但当空气或天然气通过几十米以上的细管时，密度沿程变化的"可压缩性"和壁面摩擦损失都会变得很重要。范诺流就是把"绝热（无热量进出）＋等截面＋有壁面摩擦"这三个条件同时认真求解后得到的模型。

🙋

有意思的是，把左边马赫数设为 0.3 或 0.5，再拖动管长滑块，出口马赫数会一直上升。摩擦反而让流动"变快"，这与直觉相反。

🎓

这正是范诺流最精彩的地方。不可压缩直觉认为"摩擦＝阻力＝变慢"，但气体不一样。摩擦消耗能量使静温度 T 下降，由状态方程密度 ρ 随之降低，等截面下 ρuA＝常数迫使速度 u 上升，连锁反应。再加上 M＝u/√(γRT) 中分子 u 上升、分母 √T 下降，马赫数双重提升。这就是亚音速流向 M＝1 加速的原因。

🙋

如果管长再加长，屏幕上"壅塞"会变红。这时候发生了什么？

🎓

这叫"摩擦壅塞"，是范诺流最著名的结论。入口范诺函数 fL*/D 表示"从此处通过摩擦到达 M＝1 所需的无量纲摩擦长度"，当实际管道的摩擦长度 fL/D 达到这个值时，管内某处就出现 M＝1。比此更长的管道物理上不可能存在，上游流量会被自动节流。所以设计长距离管道时，最重要的第一道检查就是"f·L/D 是否超过入口的 fL*/D"。

🙋

右边的范诺函数图中，M=1 附近曲线趋近 0，这表示"M 越接近 1 越临近壅塞"吗？

🎓

完全正确。曲线在 M＝0 处发散、M＝1 处归零，是单调递减的。所以入口马赫数越大，fL*/D 越小，距壅塞的余量越少。比如 M₁＝0.3 时 fL*/D≈5.3，M₁＝0.5 时约 1.07，M₁＝0.7 时骤降到约 0.21。所以"要在不壅塞的前提下推大流量"，技巧是加粗管径，或用增压机提高入口压力，使入口马赫数保持低位。

🙋

实际工程中出口压力比 P₂/P₁ 也很重要吧？这里显示约 0.84。

🎓

没错——管道设计的目的就是"出口剩多少压力"，P₂/P₁ 是最终判断标准。亚音速范诺流中，静压 P 沿程总是下降，本例约降 16%。而瑞利流（加热）中根据加热量大小，压力可升可降。长距离天然气管道就是典型的范诺问题，每隔一段设置压缩机站重新升压。

常见问题

在不可压缩直觉中，摩擦总是让流动变慢，但在绝热等截面的可压缩管道中情况正好相反。摩擦消耗能量，使静温度、静压与密度都下降，由等截面下的质量守恒（ρuA＝常数），速度 u 必须上升。马赫数 M＝u/√(γRT) 中分子的 u 上升、分母的 √T 下降，因此亚音速流动总是向 M＝1 加速。超音速范诺流则反向减速至 M＝1。无论亚音速还是超音速，终点都是 M＝1，这正是摩擦壅塞的本质。

摩擦壅塞指在没有截面变化、没有热量进出的条件下，仅靠壁面摩擦就在管道中途到达 M＝1 的现象。判定指标是无量纲摩擦长度 fL/D，与入口马赫数 M1 对应的范诺函数 fL*/D（从 M1 到 M＝1 所需的摩擦长度）比较。当 fL/D 达到 fL*/D 时，M＝1 出现，再长的管道在物理上无法存在。工程实践中建议按 fL/D ≤ 0.7·(fL*/D) 设计，预留余量以应对温度与流量波动。

两者都是等截面一维可压缩流，但驱动机制不同。范诺流是"绝热＋摩擦"，熵增的同时 M 趋于 1。瑞利流是"忽略摩擦＋热量进出"，亚音速下加热使 M 上升、冷却使 M 下降。两者在亚音速下都向 M＝1 加速，但范诺流是温度下降同时 M 上升，瑞利流（加热时）则是温度上升同时 M 上升。长距离管道以摩擦为主时用范诺流；燃烧室以热加为主时用瑞利流。

可通过穆迪图或科尔布鲁克方程确定。商用钢管的相对粗糙度 ε/D ≈ 5×10⁻⁵～5×10⁻⁴，雷诺数为 10⁵～10⁷ 的完全湍流区，f ≈ 0.015～0.030。常见初始值：长距离天然气主管线（D≈500mm）取 f≈0.012；HVAC 钢板风管（D≈300mm）取 f≈0.020；橡胶/塑料软管取 f≈0.030。范宁摩擦系数 f_F 与达西系数关系为 f＝4·f_F，请注意文献定义。本工具采用达西-韦斯巴赫摩擦系数。

实际应用

长距离天然气管道：横贯西伯利亚、北美大陆数千公里的管道工程，其核心正是范诺流分析——既用于压力损失核算，也用于壅塞余量校核。每隔数百公里设置压缩机站重新升压，操作上正是"在出口范诺函数值降至 0 之前送回上游"。运营商使用直径 600～1200mm 的管道，摩擦系数 f≈0.012，并将入口马赫数压低到 0.05～0.10，以最大化压缩机站间距。

建筑暖通空调（HVAC）风管设计：大型办公楼或工厂的送风管长度超过 100m，含有许多弯头与分支，摩擦损失占主导。运行流速远低于音速（M≈0.1～0.2），但对于某些直径与长度的组合，仍可能达到摩擦壅塞。利用范诺流分析可确定风机所需的静压及最终送风量。大型数据中心冷却风管也用同样的方法保证机柜间气流均衡。

燃气轮机叶片内部冷却通道：现代航空发动机涡轮入口温度超过 1500℃，叶片内部贯穿直径 1～3mm 的细冷却通道，由压气机抽气进行对流冷却。细通道中摩擦很大，流动接近范诺流。设计者需要匹配通道直径与长度，确保出口不会壅塞，从而保证冷却气量。一旦壅塞，叶片冷却不足将导致蠕变断裂，这是关键分析。

气力（粉体）输送：水泥厂、食品厂用气流输送粉体的系统是典型的长距离高速气流。粉体加入会提高等效摩擦系数，使系统处于比纯气体更紧张的范诺范围。若入口压力不足，粉体会在管道中途堵塞——"阻塞"故障。设计者通过范诺流反推"在多长处达到 M＝1"，据此布置增压风机的位置。

常见误解与注意事项

最常见的误解是把"摩擦总是让流动变慢"的不可压缩直觉照搬到可压缩流。亚音速范诺流中，摩擦让流动变快、马赫数向 1 攀升。正确的直觉是"气体在温度下降时收缩，等截面下同样的质量流必须流得更快"。超音速范诺流则相反，摩擦使流动减速至 M＝1，这同样源自"M＝1 是摩擦的终点"这一普适原理。"亚音速与超音速均向 M＝1 收敛"是要点，只回答"加快"或"减慢"都只对了一半。

其次是摩擦系数 f 的定义混淆。范诺流文献中"达西摩擦系数 f"与"范宁摩擦系数 f_F"并存，关系为 f＝4·f_F。式 fL*/D 中的 f 是达西系数（本工具也采用达西），但基于范宁系数的教科书会写成 4fL*/D。两者混用会导致致命错误：会在管长仅为实际 1/4 处误判出壅塞。穆迪图采用达西，许多国外化工手册采用范宁，请务必核实文献。

最后是忘记"范诺流严格假设绝热"。实际管道与环境存在热交换，尤其是长距离管道埋于地下时会向土壤温度冷却（或加热）。本工具结果基于"绝热、等截面、一维、理想气体"的理想化条件。若工程精度要求高于 10%，则需要将范诺（摩擦）与瑞利（热）组合的广义可压缩管流数值解（分段积分或 CFD）。即便如此，范诺分析作为"摩擦是否占主导""壅塞余量多少"的筛选工具仍极为有效，至今仍是正式 CFD 分析前不可或缺的方法。