穆迪图模拟器 返回
管道流动模拟器

穆迪图模拟器 — 管道摩擦系数

基于雷诺数与相对粗糙度,在穆迪图上可视化圆管内流动的摩擦系数 f。Swamee-Jain 显式公式即时给出湍流的 f,再由达西-韦斯巴赫方程算出压降与水头损失。

参数设置
log₁₀ Re
log
Re = 1.00e+5
log₁₀ ε/D
log
ε/D = 1.00e-3
长径比 L/D
平均流速 V
m/s

假定流体为水(ρ = 1000 kg/m³,g = 9.81 m/s²)。Re 与 ε/D 滑块采用对数刻度。

计算结果
摩擦系数 f
压降 ΔP
水头损失 h_L
流动区域
穆迪图(f vs Re)

横轴=Re(对数)/纵轴=f(对数)/曲线=相对粗糙度 ε/D 对应的摩擦系数/黄点=当前工况点

流速 V 与压降 ΔP

在当前 Re、ε/D、L/D 下变化 V 时的压降曲线(黄点=当前值)

理论与主要公式

达西-韦斯巴赫方程统一表示了圆管内由摩擦引起的压降。摩擦系数 $f$ 是雷诺数 $Re$ 与相对粗糙度 $\varepsilon/D$ 的函数,层流区使用解析式,湍流区使用 Swamee-Jain 显式公式。

达西-韦斯巴赫方程(压降与水头损失):

$$\Delta P = f\,\frac{L}{D}\,\frac{\rho V^2}{2},\qquad h_L = \frac{\Delta P}{\rho g}$$

层流($Re \lt 2300$)的摩擦系数:

$$f = \frac{64}{Re}$$

湍流($Re \geq 4000$)的 Swamee-Jain 公式:

$$f = \frac{0.25}{\left[\log_{10}\!\left(\dfrac{\varepsilon/D}{3.7} + \dfrac{5.74}{Re^{0.9}}\right)\right]^2}$$

$\rho$ 为密度 [kg/m³],$V$ 为平均流速 [m/s],$L$ 为管长,$D$ 为管径,$\varepsilon$ 为管壁粗糙度 [m],$g$ 为重力加速度 [m/s²]。

什么是穆迪图模拟器

🙋
我想算一下管道的压降,可教科书只说"从穆迪图上读 f",从图上读太麻烦了,有没有公式可以直接算?
🎓
原始的方程是 Colebrook 式,$\frac{1}{\sqrt{f}} = -2\log_{10}\!\left(\frac{\varepsilon/D}{3.7}+\frac{2.51}{Re\sqrt{f}}\right)$,但 f 同时出现在两边,是隐函数,需要迭代求解。实务上常用 Swamee-Jain 显式式,一次代入就能直接算出 f,误差大概 1% 左右。这个工具用的就是它。
🙋
原来如此。那层流的时候呢?
🎓
层流(Re < 2300)就更简单,$f = 64/Re$,直接是哈根-泊肃叶流动的解析解。它完全由黏性应力决定,所以粗糙度 ε/D 一点都不起作用。在模拟器里把 log₁₀(Re) 调到 3 以下试试,你会看到工况点正好落在橙色直线上(双对数下斜率为 −1)。
🙋
过渡区(Re=2300〜4000)是怎么算的呢?
🎓
说实话,过渡区"f 不能唯一确定"才是真相,流动会在层流和湍流之间反复切换。设计上为了偏安全,通常用湍流公式(Colebrook 或 Swamee-Jain)来评估。本工具在过渡区也是直接套用湍流公式,并在图上用红色带状区域标出"过渡区"。
🙋
"相对粗糙度 ε/D"具体应该取多少呢?
🎓
典型值:拉拔加工的铜管、PVC 大约 ε ≈ 1.5 μm,商用钢管约 45 μm,镀锌钢管约 150 μm,老旧铸铁管约 260 μm。除以管内径 D 就得到 ε/D。比如 D=100 mm 的钢管,ε/D ≈ 4.5e-4,把滑块设到 log₁₀(ε/D)=−3.35 附近,就会落到对应的曲线上。

常见问题

Swamee-Jain 公式(1976)是 Colebrook 方程的显式近似,在适用范围 5e3 ≤ Re ≤ 1e8 且 1e-6 ≤ ε/D ≤ 1e-2 内,与 Colebrook 解的相对误差大致在 1% 以内。对工程压降计算来说精度完全够用,又免去了对隐式 Colebrook 方程的迭代,可直接在电子表格中使用,因此在水泵选型与管道设计实务中被广泛采用。
本工具使用的是达西摩擦系数 f_D(也称穆迪摩擦系数),与化工领域常用的范宁摩擦系数 f_F 之间满足 f_D = 4·f_F。例如层流时 f_D = 64/Re,但范宁记法写成 f_F = 16/Re。阅读文献时务必确认所用公式是按达西基准(与达西-韦斯巴赫式配对)还是按范宁基准书写。
严格意义上不能直接用,但作为近似可以用"水力直径"D_h = 4A/P(A = 流通截面积,P = 湿周)替代圆管直径 D,从而把穆迪图扩展到矩形风管、环形流道等场合。湍流下这种代换的精度足以满足工程使用,但层流下需要按截面形状的修正系数(如矩形约 0.89)。这是通风管道与换热器管侧计算中常用的近似方法。
本工具仅计算直管段的摩擦损失(沿程损失),不包含弯头、三通、阀门、突扩等带来的局部损失。实务上常用每个管件的损失系数 K,按 ΔP_minor = K·(ρV²/2) 累加,或换算成等效长度 L_eq 加到直管长度 L 中。在长管系中沿程摩擦占主导,但在短管或管件较多的系统中,局部损失所占比例不可忽略。

实际应用

水泵与输液系统的功率计算:在估算水泵所需扬程与功率时,先从穆迪图读取 f,再用达西-韦斯巴赫方程计算管系摩擦损失。从过程工业的长距离管线、楼宇给水与空调水系统,到核电厂一回路冷却水回路,所有"输送流体"的设计计算几乎都从这里出发。f 翻一番,水泵功率几乎也翻一番,因此管径选择与粗糙度管理直接关系到经济性。

HVAC 与通风风管设计:空调、通风风管在用上述水力直径换算成圆管之后,使用同一张穆迪图。它常用于风机静压估算,ASHRAE 手册中给出了各种风管材料的 ε(如镀锌铁皮约 0.15 mm)。把长风管做细会带来风量大幅下降与能耗骤增,因此设计时要在图上比较"等效"摩擦损失再决定路径与尺寸。

石油与天然气管线:长距离输送中摩擦损失占成本的大部分,因此常采用低粗糙度内涂层(如环氧)把 ε/D 下降一个数量级,从而减少所需的泵站/压气站数量。在输送效率评估中,常把 Colebrook 式(或 Swamee-Jain 式)与黏度的温度修正结合,并利用高 Re 区"完全粗糙流"中 f 的稳定区域做鲁棒设计。

核电、火电厂的热工水力:反应堆一回路、蒸汽发生器、冷凝器等含有大量管束的系统,各支路的流量分配由摩擦损失决定,因此把穆迪图给出的 f 嵌入到一维热工水力程序(RELAP、TRACE 等)中。事故工况下的冷却剂流动分析中,依赖 Re 与 ε/D 的 f 也是安全评估的关键参数。

常见误区与注意点

最常见的误解是"管道越粗糙摩擦总是越大"。看穆迪图就能发现,在低 Re 区(层流〜过渡区)几乎看不出 ε/D 的影响,所有粗糙度的曲线都汇合到 64/Re 的同一条直线上。粗糙度起作用是在湍流区,而且 Re 越高,效果越显著。可以在模拟器中把 log₁₀(Re) 设为较低值(例如 3.0),改变 ε/D 时观察 f 几乎不变,这就是"水力光滑"状态——黏性底层把粗糙度凸起完全覆盖住了。

第二个常见误解是认为过渡区(Re=2300〜4000)的 f 可以唯一确定。实际上该区域内流动在层流与湍流之间来回切换,处于不稳定状态,f 随时间波动且重复性差,因此教科书的穆迪图常以虚线或留白方式标注此带。设计上的常规做法是用湍流公式做偏安全评估,或者尽可能选择管径与流速,使工况避开过渡区。

最后一个要警惕的是把达西-韦斯巴赫式中的 f 与"范宁摩擦系数"搞混。化工类文献常用范宁记法 f_F = f_D/4,层流公式写成 f_F = 16/Re。如果只看 f 的数值就套用公式,会出现 4 倍的差异。必须把手头的公式和 f 的定义成对确认——这是横向阅读流体工学教材时必须遵守的铁律。本工具自始至终都使用达西基准。