2. 流体力学基础回顾
在进入除尘脱硫脱硝系统的阻力核算之前，有必要把流体力学中几条最基础、也最常用的规律重新梳理一遍。系统阻力本质上就是烟气在管道、设备、塔器内流动时不断损失机械能的过程，而描述这一过程的工具无非三件：能量守恒（伯努利方程）、阻力分类（沿程与局部）、以及决定阻力大小的流动状态参数（雷诺数与摩擦系数）。本节把这些工具逐一展开，并给出可直接用于后续章节的管道阻力计算方法。
一、伯努利方程：能量守恒在流动中的表达
对于不可压缩流体的定常流动，沿同一流线取两个截面1和2，伯努利方程可写为：
p₁/ρg + z₁ + v₁²/2g = p₂/ρg + z₂ + v₂²/2g + h_f式中 p 为静压（Pa），z 为位置高度（m），v 为流速（m/s），ρ 为流体密度（kg/m³），g 取 9.81 m/s²，h_f 为两截面间的总水头损失（m）。三项分别代表压力能、位能和动能，h_f 则是流动过程中不可逆地转化为热能的机械能损失。这里约定：沿流动方向取截面1在上游、截面2在下游，h_f 恒为正，表示能量减少。
对除尘脱硫脱硝系统而言，烟气密度随温度变化显著，严格来说属于可压缩流动，但在风机全压不超过约 15 kPa、马赫数远小于 0.3 的常规工况下，按不可压缩处理误差很小，工程上普遍采用这一简化。真正需要注意的是密度取值：应以系统实际运行温度下的烟气密度代入，而不是常温空气密度。例如 150 ℃ 烟气密度约为 0.83 kg/m³，仅为常温空气的 0.68 倍，同样的流速下动压明显偏低。
伯努利方程在阻力核算中的直接用途是“由压差反推阻力”或“由阻力预估所需压头”。当两截面标高相同、流速相近时，方程简化为 Δp = ρg·h_f，即阻力损失几乎全部体现为静压下降。这正是现场用U形管或微压计测量设备进出口静压差来估算阻力的理论依据。
二、沿程阻力与局部阻力
系统总阻力由两部分构成。沿程阻力是流体在直管段内因内摩擦产生的损失，与管长成正比；局部阻力是流体经过弯头、变径、阀门、三通、进出口等局部构件时，因流线急剧改变、产生涡流而造成的损失。二者叠加即为总阻力：
Δp = Δp_沿程 + Δp_局部 = (λ·L/d + Σζ)·(ρv²/2)其中 λ 为沿程阻力系数，L 为管道长度（m），d 为管道内径（m），ζ 为局部阻力系数，ρv²/2 为动压（Pa）。这个式子说明一个重要事实：阻力与流速的平方成正比。流速提高一倍，阻力约变为四倍。因此在系统阻力超标时，降低局部流速往往比更换风机更经济。
局部阻力系数 ζ 是查表或由厂家提供的经验值，常见取值示例如下（示例数据，供教学参考，实际应以设计手册或设备样本为准）：
| 局部构件 | 典型 ζ 值 | 说明 |
|---|---|---|
| 90° 标准弯头 | 0.75～1.0 | 曲率半径越小取值越大 |
| 突然扩大（A₁/A₂=0.5） | 约 0.56 | 按 (1−A₁/A₂)² 估算 |
| 突然缩小（A₁/A₂=0.5） | 约 0.35 | 入口圆滑可显著降低 |
| 全开闸阀 | 0.1～0.2 | 半开时急剧增大 |
| 管道入口（锐边） | 0.5 | 喇叭口可降至 0.05 |
| 管道出口 | 1.0 | 动能全部损失 |
工程上常把局部阻力折算成“当量长度”，即用一段直管的沿程阻力等效某个局部构件，便于统一累加。但更推荐的做法是分别列出沿程项与局部项，因为局部构件的数量在系统布置阶段是可以优化的，分开列账才能看出哪一部分值得削减。
三、雷诺数与摩擦系数
流动状态由雷诺数判定：Re = ρvd/μ = vd/ν，其中 μ 为动力黏度，ν 为运动黏度。圆管流动以 Re ≈ 2300 和 4000 为界，分为层流、过渡区和湍流。除尘脱硫脱硝系统的烟气管道内 Re 通常在 10⁵～10⁶ 量级，几乎全部处于湍流粗糙区或过渡粗糙区，层流公式基本用不上，但理解 Re 的物理意义仍然重要——它表征惯性力与黏性力之比，Re 越大，惯性作用越强，涡流和掺混越剧烈，阻力系数趋于稳定。
湍流区的沿程阻力系数 λ 取决于 Re 和相对粗糙度 ε/d。常用计算式有：
- 光滑管区（Blasius 式）：λ = 0.3164·Re⁻⁰·²⁵，适用于 Re < 10⁵；
- 湍流粗糙区（Colebrook-White 式，需迭代或查 Moody 图）：1/√λ = −2lg[ε/(3.7d) + 2.51/(Re√λ)]；
- 完全粗糙区（Shifrinson 式）：λ = 0.11·(ε/d)⁰·²⁵，此时 λ 仅与相对粗糙度有关，与 Re 无关。
对碳钢烟道，绝对粗糙度 ε 一般取 0.1～0.2 mm；内衬玻璃鳞片或衬胶的管道可取 0.05 mm 左右。可以看出，在完全粗糙区，管道材质和锈蚀状况对阻力的影响比流速之外的其他因素更直接，这也是为什么老旧烟道阻力往往高于设计值。
四、管道阻力计算示例
下面用一个教学示例走通计算流程。设某段烟气管道：内径 d = 1.2 m，长度 L = 30 m，烟气温度 150 ℃，密度 ρ = 0.83 kg/m³，运动黏度 ν ≈ 2.4×10⁻⁵ m²/s，流速 v = 12 m/s，碳钢管道 ε = 0.15 mm，管段含 3 个 90° 弯头（ζ 各取 0.9）、1 个全开闸阀（ζ = 0.15）、1 个锐边入口（ζ = 0.5）。
动压：ρv²/2 = 0.83 × 12² / 2 ≈ 59.8 Pa
雷诺数：Re = vd/ν = 12 × 1.2 / 2.4e-5 = 6.0×10⁵
相对粗糙度：ε/d = 0.15/1200 ≈ 1.25×10⁻⁴
按 Colebrook 式迭代得 λ ≈ 0.0175
沿程阻力：Δp_沿 = λ·(L/d)·动压 = 0.0175 × (30/1.2) × 59.8 ≈ 26.2 Pa
局部系数和：Σζ = 3×0.9 + 0.15 + 0.5 = 3.35
局部阻力：Δp_局 = 3.35 × 59.8 ≈ 200.3 Pa
该管段总阻力：Δp ≈ 226.5 Pa这个结果传递出一个很实际的信号：在短管段中，局部阻力往往远大于沿程阻力（本例约 7.6 倍）。因此系统阻力核算时，弯头、变径、进出口的数量与形式必须逐项统计，不能只按管长估算。若把 3 个弯头改为曲率半径更大的弯管（ζ 降至 0.4），局部阻力可降到约 130 Pa，整段阻力下降三成以上。
五、为系统阻力分析做的三点准备
把上述工具用到除尘脱硫脱硝系统上，需要提前明确三件事。第一，统一基准与符号约定：全课程以流动方向为正，上游截面为1、下游截面为2，阻力损失一律取正值，压差测量以“进口静压减出口静压”为正。第二，分段累加：把系统拆成烟道段、除尘器、脱硫塔、脱硝反应器、风机进出口等若干段，每段分别计算沿程与局部阻力，再串联求和，这样任何一段超标都能定位。第三，密度与流速取实际值：温度变化会引起体积流量和密度变化，同一质量流量的烟气在高温段流速高、密度低，动压计算必须用当地实际参数，不能全系统套用同一个密度。
需要提醒的是，本节给出的 λ、ζ 数值和算例结果均为教学示例，用于说明计算方法和量级关系。实际工程中应以设计手册、设备厂家样本和现场烟气参数为准，并在系统投运后通过实测压差对计算值进行校核。掌握了伯努利方程、阻力分类、雷诺数与摩擦系数这几件工具，后续分析除尘器、脱硫塔、脱硝反应器及风机全压时，就能把“阻力从哪来、有多大、怎么降”这条线索贯穿起来。