暖气管里结锈瘤比长锈层更费泵吗:同样一堆锈摊平阻力涨1.96倍,结成6.3毫米瘤涨2.97倍,碎成3.2毫米小瘤反而只涨1.7倍,临界点是4.07毫米

钢管用上十几年,内壁一定会长东西。切开看,有的是一层均匀发黄鼓起的锈壳,有的却是一颗一颗鼓出来的锈瘤——像墙上起的疹子,最大的能有绿豆那么大。
直觉上,这两种"生锈"只是外观不同,反正都是把管子堵细了。但把阻力算出来会发现,它们走的完全是两套物理:均匀锈层增加的是摩阻,锈瘤增加的是形阻。而形阻这一项,此前的算法里根本没被单独拎出来过。
结论先摆在这里:在 DN20、270 L/h、60 ℃ 的常规支管工况下,同样是"平均锈层厚度 0.5 毫米"那么一堆锈——摊平成锈壳,阻力涨 1.96 倍;结成半径 6.3 毫米的锈瘤,涨 2.97 倍;但如果碎成半径 3.2 毫米的密集小瘤,只涨 1.70 倍,比摊平还省。分界线是瘤半径 4.07 毫米。
换句话说:决定阻力的不只是"锈有多少",还有"锈长成什么样"。同样一堆锈,聚成大瘤比摊匀了费泵,碎成小瘤反而比摊匀了省泵。
一、先分清两件事:摩阻和形阻不是一回事
流体在管里损失压力,只有两条路:
第一条:摩阻(skin friction)。水贴着管壁走,壁面粗糙,一层一层的水互相拖拽。这一项由沿程摩阻系数 λ描述,λ 由雷诺数 Re 和相对粗糙度 ε/d 决定,用 Colebrook 隐式式求解。均匀锈层走的就是这条路——它把整根管子的 ε 抬高,同时把通径从 d 缩到 d−2δ。
第二条:形阻(form drag)。水流撞上一个凸起,在凸起前面减速增压、后面分离出一个尾流涡,动能白白耗散掉。这一项由局部阻力系数 ζ描述,只跟凸起的形状与迎风面积有关,跟 Re 基本无关。锈瘤走的是这条路——它几乎不改变通径,但每颗瘤都是一个小钝体。
这个区别是全部结论的支点。摩阻是"整根管子均匀变糙",形阻是"管子里摆了一串障碍物"。前一篇《暖气管内壁变粗糙会更费泵吗》已经把第一条路用 Colebrook 完整走了一遍,也明确写了自己没有算第二条:"只算了均匀锈层——实际锈瘤是局部的、离散分布的,局部凸起产生的额外形阻与均匀粗糙度不等价,本文未建模。"本文就是补上这一条。
二、基准工况:先把数字跟前作对齐
为了能和前作互相校验,本文沿用完全相同的基准:
| 参数 | 取值 | 说明 |
|---|---|---|
| 管子 | DN20 钢管,内径 d = 20.0 mm | 户内支管常见规格 |
| 水温 | 60 ℃ | ρ = 983 kg/m³,ν = 4.78×10⁻⁷ m²/s |
| 流速 | v₀ = 0.239 m/s | 对应流量 Q₀ = 270.3 L/h |
| 雷诺数 | Re₀ = 10 000 | 与前作完全一致 |
| 新钢管粗糙度 | ε₀ = 0.1 mm | 《实用供热空调设计手册》取值 |
| 沿程阻力 | Δp₀/L = 52.82 Pa/m | λ₀ = 0.03763 |
三个交叉校验点全部对上:
- Colebrook 在 ε→0 的极限下退化为 Prandtl-Kármán 光滑管式,本文 Re=10 000 算得 λ = 0.03089,前作同为 0.03089;
- 同 Re 下 Blasius 式给 0.03164,两篇一致;
- 新钢管 0.1 mm 粗糙度使 λ 从 0.03089 抬到 0.03763,吃掉 21.8% 摩阻,与前文"已吃掉两成"完全吻合。
基准站住了,下面加锈。
三、一颗锈瘤值多少阻力:Borda-Carnot 突扩
把一颗锈瘤理想化为贴在管壁上的半球,底半径 r、高度 k = r。它在水流方向上的迎风投影是一个半圆盘:
Af = πr²/2
相对管子截面积 A₀ = πd²/4,单颗瘤的阻塞比为:
a = Af/A₀ = 2r²/d²
水流过这颗瘤时,自由截面从 A₀ 缩到 A₀−Af,流速从 v 提到 v/(1−a),绕过瘤之后又扩回 A₀。这个"缩了又扩"的过程,不可逆损失几乎全在扩张那一步 —— 按 Borda-Carnot 突扩公式,以管内流速 v 为基准:
ζ = (1/(1−a) − 1)² = (a/(1−a))²
这条式子有个重要性质:ζ 与流速无关。不像 λ 会随 Re 变化,ζ 是一个纯几何量。这一点在第十节会成为判据的突破口,也是它的死穴。
还有一个物理上限必须守住:半球的高度 k = r 不能超过管子半径 d/2,否则瘤就长到对面壁上去了。代入 a = 2r²/d² 得:
r ≤ d/2 = 10 mm ⇒ a ≤ 0.50
单颗锈瘤最多只能堵掉 50% 的截面。这个上限后面会用到。
四、把"同样一堆锈"这句话变成方程
要让"摊平"和"结瘤"可比,必须锁死一个量:锈的总体积。
定义 δ 为摊平厚度——把管壁上所有锈刮下来、均匀糊回整个内表面 πdL 上,得到的厚度。这样 δ 就是现场能测、能说的量(内窥镜量、割管量、甚至用除下来的锈重反算):
δ = V总 / (πdL)
若这些锈结成半球瘤,每颗体积 Vb = 2πr³/3,每米有 n 颗,则每米锈体积为 n·Vb = πdδ,于是:
n = πdδ / Vb = 3dδ / (2r³)
再配合 r = d·√(a/2),δ 一旦定死,a 和 n 就被一根绳子拴住了:瘤越大(a 大、r 大),个数必然越少(n 小);瘤越小,个数必然越多。这条约束是后面所有反直觉结论的来源。
于是单位长度的总阻力可以拆成两块:
Δp/L = λ₀/d · ρv²/2 (沿程摩阻,管径未变,仍用新钢 λ₀)
+ n · ζ · ρv²/2 (形阻,逐颗累加)
两边一除,得到一个非常干净的比值:
Δp形/Δp摩 = n·d·ζ / λ₀
把 n、r、ζ 全部用 a 和 η = δ/d 代掉,整理后有闭式:
Δp形/Δp摩 = (3η/2) · 23/2 · a1/2 / [ (1−a)² · λ ]
五、阻力跟着"瘤有多大"单调上升
把 DN20、δ = 0.5 mm 代进去,只改变瘤的大小(同时按体积守恒自动改变个数),结果如下:
| 阻塞比 a | 瘤半径 r (mm) | 每米个数 n | ζ | 形阻/摩阻 | 阻力倍数 | 定压差流量剩 |
|---|---|---|---|---|---|---|
| 0.01 | 1.41 | 5303 | 0.0001 | 0.29 | 1.29 | 88% |
| 0.02 | 2.00 | 1875 | 0.0004 | 0.42 | 1.42 | 84% |
| 0.05 | 3.16 | 474 | 0.003 | 0.70 | 1.70 | 77% |
| 0.10 | 4.47 | 168 | 0.012 | 1.10 | 2.10 | 69% |
| 0.15 | 5.48 | 91 | 0.031 | 1.51 | 2.51 | 63% |
| 0.20 | 6.32 | 59 | 0.062 | 1.97 | 2.97 | 58% |
| 0.30 | 7.75 | 32 | 0.184 | 3.15 | 4.15 | 49% |
| 0.40 | 8.94 | 21 | 0.444 | 4.95 | 5.95 | 41% |
| 0.50 | 10.00 | 15 | 1.000 | 7.97 | 8.97 | 33% |
看最后一列:同样是那么多锈,阻力从 1.29 倍一路涨到 8.97 倍,跨度接近 7 倍。唯一的变量就是锈"结"成什么尺寸。
最反直觉的一段在最上面:把锈碎成半径 3.2 毫米以下的小瘤,阻力倍数低于摊平锈壳的 1.96 倍。也就是说,在某些情况下,"结瘤"反而比"摊匀"更省泵。
六、为什么:ζ 涨得比 n 掉得快
这个单调性不是巧合,可以从三个指数上直接看出来。
由体积守恒,n ∝ 1/r³ ∝ a−3/2。而 ζ = (a/(1−a))²,在 a 不大时 ζ ≈ a²。两者相乘:
n·ζ ∝ a−3/2 · a² = a1/2
总形阻正比于 a 的二分之一次方,单调递增。物理上讲就是:把锈从 100 颗合并成 1 颗,颗数少了 100 倍,但单颗的阻力涨了不止 100 倍——因为 ζ 对阻塞比是超线性的(近似平方),阻塞比翻倍、阻力翻四倍。
再看 (1−a)² 那一项:当 a 接近 0.5 的上限时它开始发威,把曲线往上狠狠一抬。a 从 0.4 到 0.5,个数只减了 29%,阻力却从 5.95 倍跳到 8.97 倍。所以真正致命的不是"有很多锈瘤",而是"有几颗特别大的锈瘤"。
七、临界瘤半径:4.07 毫米
把"摊平"这条基准线也算出来,就能给出一条可以直接用的判据。
摊平锈壳按前作的模型处理:通径 d′ = d − 2δ,粗糙度取 ε = δ,定流量下重新解 Colebrook:
| δ (mm) | d′ (mm) | v′ (m/s) | Re′ | λ′ | 阻力倍数 | 定压差流量剩 |
|---|---|---|---|---|---|---|
| 0.1 | 19.8 | 0.244 | 10 101 | 0.0376 | 1.05 | 97.5% |
| 0.2 | 19.6 | 0.249 | 10 204 | 0.0432 | 1.27 | 88.7% |
| 0.5 | 19.0 | 0.265 | 10 526 | 0.0572 | 1.96 | 71.4% |
| 1.0 | 18.0 | 0.295 | 11 111 | 0.0771 | 3.47 | 53.7% |
| 2.0 | 16.0 | 0.373 | 12 500 | 0.1167 | 9.46 | 32.5% |
(δ = 0.5 mm 处,定压差流量掉到 71.4%,即衰减 28.6%。前作在同一锈厚下含泵曲线算出衰减 34.1%,同量级;差额来自泵特性曲线与本文未计的局部项。)
现在求"离散 = 均匀"的临界瘤半径 r*。由于离散阻力随 a 单调上升,这个交叉点唯一,用二分即可:
| δ (mm) | 摊平阻力倍数 | 临界阻塞比 a* | 临界瘤半径 r* | r*/d |
|---|---|---|---|---|
| 0.1 | 1.05 | 0.0081 | 1.27 mm | 0.064 |
| 0.2 | 1.27 | 0.0476 | 3.09 mm | 0.154 |
| 0.5 | 1.96 | 0.0828 | 4.07 mm | 0.203 |
| 1.0 | 3.47 | 0.1168 | 4.83 mm | 0.242 |
| 2.0 | 9.46 | 0.2145 | 6.55 mm | 0.327 |
读法很直接:锈的平均厚度是 0.5 毫米时,如果锈瘤半径小于 4.07 毫米,结瘤比摊平省泵;大于 4.07 毫米,结瘤开始比摊平费泵。锈越厚,这条门槛越宽松(因为摊平侧的惩罚涨得更快)——但这不代表结瘤变安全了,只是"摊平"这条参照线本身恶化得更快。
八、同样的锈,摊平和结瘤差多少:一张对照表
| δ (mm) | 摊平 | a=0.02 | a=0.05 | a=0.10 | a=0.20 | a=0.30 | a=0.40 | a=0.50 |
|---|---|---|---|---|---|---|---|---|
| 0.1 | 1.05 | 1.14 | 1.14 | 1.22 | 1.39 | 1.63 | 1.99 | 2.59 |
| 0.2 | 1.27 | 1.17 | 1.28 | 1.44 | 1.79 | 2.26 | 2.98 | 4.19 |
| 0.5 | 1.96 | 1.42 | 1.70 | 2.10 | 2.97 | 4.15 | 5.95 | 8.97 |
| 1.0 | 3.47 | 1.83 | 2.40 | 3.20 | 4.94 | 7.30 | 10.90 | 16.94 |
| 2.0 | 9.46 | 2.66 | 3.79 | 5.40 | 8.88 | 13.60 | 20.81 | 32.89 |
这张表最值得看的是左半边:在锈还不算厚(δ ≤ 1 mm)时,只要锈是碎成小瘤的形态,阻力倍数一直低于摊平值。这意味着"管壁上有鼓包"并不自动等于"阻力超标"——它也可能是锈长得比较"碎"。
九、管径的惩罚:形阻是 d⁻⁶,不是 d⁻⁴·⁷⁵
前作给出了一个很好用的标度:定流量下沿程阻力 Δp ∝ d−4.75(λ ∝ d0.25 与 Δp ∝ λd−5 相乘)。这个标度在结瘤时会失效。
推导形阻项的管径标度。定流量下 v ∝ d−2,故 v² ∝ d−4;而 n·ζ ∝ δ·a1/2/d²,于是:
Δp形 ∝ (n·ζ)·v² ∝ d−2 · d−4 = d−6
用四个规格数值验证(δ = 0.5 mm,a = 0.20,定流量 270 L/h):
| 规格 | d (mm) | v (m/s) | 摊平倍数 | 结瘤倍数 | 形阻项 (Pa/m) |
|---|---|---|---|---|---|
| DN15 | 15.7 | 0.388 | 2.28 | 3.48 | 444.6 |
| DN20 | 20.0 | 0.239 | 1.96 | 2.97 | 104.0 |
| DN25 | 26.0 | 0.141 | 1.70 | 2.50 | 21.6 |
| DN32 | 34.0 | 0.083 | 1.50 | 2.12 | 4.3 |
对最后一列做对数拟合,指数为 −6.00,与理论 d−6 完全吻合;同组沿程摩阻项拟合得 −4.96,落在前作 d−4.75 附近(偏差来自本文 Re 已部分偏离纯 Blasius 区)。
工程含义:结瘤对细管的惩罚比锈层更狠。DN15 上,形阻项高达 444.6 Pa/m,是 DN20 的 4.3 倍;而同一组锈摊平时,DN15 只比 DN20 高 16%。所以"细管先堵"这句老经验,在结瘤情形下要被放大好几倍。
十、现场怎么判断:X 判据
现在给一条能落地的判据。定义:
X = 实测阻力 ÷ 按实测平均锈厚摊平算出的阻力
X 大于 1 的部分,就是形阻——也就是"锈是结成瘤的"这一事实带来的额外代价。
| δ (mm) | a=0.05 | a=0.10 | a=0.20 | a=0.35 | a=0.50 |
|---|---|---|---|---|---|
| 0.2 | 1.01 | 1.13 | 1.41 | 2.03 | 3.29 |
| 0.5 | 0.86 | 1.07 | 1.51 | 2.52 | 4.57 |
| 1.0 | 0.69 | 0.92 | 1.42 | 2.56 | 4.88 |
用法:内窥镜或割管量出平均锈厚 δ,按第七节的摊平表读出"本该涨到几倍",再拿实测压差去除。X 在 1.0 附近 = 锈基本是摊平的(均匀腐蚀);X 在 1.5~2.5 = 明显结瘤;X 超过 3 = 已经长成少数几颗大瘤,这时候靠加泵档是补不回来的,该考虑换管或化学清洗。
注意 X < 1 也是可能的(表左上角),它的含义是"锈长得比摊平更有利",不是测错了。
十一、为什么"变流量测压差"这条老办法在这里失灵
常见做法是测两个流量下的压差,看 Δp ∝ vⁿ 的指数 n,用 n 判断流态。混合阻力下:
n = (1.75 + 2R)/(1 + R), R = 形阻/摩阻
| R | 0 | 0.5 | 1 | 2 | 3 | 5 | 10 |
|---|---|---|---|---|---|---|---|
| n | 1.750 | 1.833 | 1.875 | 1.917 | 1.938 | 1.958 | 1.977 |
| Δp(2v)/Δp(v) | 3.364 | 3.564 | 3.668 | 3.775 | 3.830 | 3.886 | 3.937 |
看起来 n 从 1.75 走到 1.98,跨度够大。问题是:重锈的均匀锈层也会把 n 推到同一个位置。δ = 0.5 mm 的摊平工况,Re 从 10 526 降到 5 263 时 λ 只从 0.0572 变到 0.0599,n ≈ 1.93;而 a = 0.20 的结瘤工况算出来 n ≈ 1.94。两者完全重叠。
所以用 n 只能判断"有没有脱离光滑管摩阻",判断不了"是锈层还是锈瘤"。要分开这两者,只能回到第十节的 X 判据——它需要的不是流量扫描,而是一个独立的锈厚读数。这条边界值得写清楚,免得被拿去误用。
十二、敏感性:把不确定度摊开
形阻系数 ζ = (a/(1−a))² 取的是纯 Borda-Carnot 突扩,对应"瘤后面完全分离"的上界;真实锈瘤若是平缓的鼓包,尾流更弱,ζ 会小一些。引入修正系数 κ,取 0.5 / 1.0 / 1.5 / 2.0 扫描(δ = 0.5 mm):
| a | κ=0.5 | κ=1.0 | κ=1.5 | κ=2.0 |
|---|---|---|---|---|
| 0.05 | 1.35 | 1.70 | 2.05 | 2.40 |
| 0.10 | 1.55 | 2.10 | 2.65 | 3.20 |
| 0.20 | 1.98 | 2.97 | 3.95 | 4.94 |
| 0.30 | 2.58 | 4.15 | 5.73 | 7.30 |
| 0.50 | 4.99 | 8.97 | 12.96 | 16.94 |
即使把形阻砍半(κ = 0.5),"阻力随瘤尺寸单调上升"和"存在临界瘤半径 r*"这两条定性结论一条都不变,只是 a = 0.20 处的阻力倍数从 2.97 降到 1.98,恰好回到摊平线附近。所以 κ 影响的是"分界线在哪里",不影响"分界线存在"。
十三、本文不算的部分
- 只用半球:真实锈瘤是球冠、是连成片的脊、是带尖的结核,形状千差万别。本文的 ζ 只由阻塞比 a 决定,没有区分陡峭度。扁而缓的鼓包真实 ζ 会低于本文值,用 κ 打折即可。
- 瘤是均匀分布的:实际锈瘤成簇、成带,局部密度远高于平均。成簇会抬高局部流速、让尾流互相干扰,真实阻力可能更高。
- 没有算锈瘤对摩阻的反作用:本文假设非瘤区管壁仍是新钢(ε = 0.1 mm)。实际上瘤周围往往伴随一片粗糙区,沿程项会被抬高,本文给的是下界。
- 没算锈瘤之间的相互遮蔽:当 a = 0.5、每米 15 颗时,间距已接近瘤径,尾流会打到下一颗瘤上,ζ 不再可简单相加。
- 单相、等温:未考虑含气两相流、生物粘泥的流变、以及锈瘤表面松散层被冲刷剥离的非稳态过程。
- 没算水锤:锈瘤改变了局部波阻抗,但那是瞬变问题,与本文的稳态阻力是两回事。
- 不涉及"该不该换管"的经济账——那是另一个问题。
十四、一句话收尾
管子里长锈这件事,直觉只问"长了多少"。但算完会发现,还得再问一句"长成了什么样"。
同样一堆锈,摊匀了阻力涨 1.96 倍,结成 6.3 毫米的瘤涨 2.97 倍,碎成 3.2 毫米的小瘤只涨 1.70 倍——三者相差五倍多,而锈的量一模一样。分界点在瘤半径 4.07 毫米:小于它,结瘤是好事;大于它,结瘤才是麻烦。
而最该记住的其实是最上面那一行:让人害怕的从来不是"管壁上有鼓包",而是"有那么几颗特别大的鼓包"。一颗堵到 50% 的瘤,单颗的 ζ 就等于 1.0,一颗顶得上十米直管。除锈、换管、加泵的优先顺序,应该按最大的那颗瘤来排,而不是按平均锈厚来排。
相关产品推荐
结合本文场景,瑞华特可提供以下采暖产品和配套选型建议:
- PPR采暖主管 DN25/DN32:用于壁挂炉、分水器和采暖主管连接,适合地暖与暖气片混合系统配套。
- 铜铝复合暖气片 TAL-800:铜铝复合暖气片TAL-800适合集中供暖、自采暖和水质复杂场景,兼顾耐腐蚀、快散热和长期稳定。核心参数为材质铜管+铝合金,高度800mm,工作压力1.2MPa。
- PERT-II阻氧地暖管 DN20:适用于家庭地暖盘管和工程地暖系统,柔韧性好,适合低温热水地面辐射采暖。
需要确认规格、数量或工程采购清单,可联系瑞华特暖气片厂家获取免费选型指导。



