4. 气相流动压降的通用关联式
干填料层的压降计算是填料塔流体力学分析的基础。所谓“干塔”，是指塔内只有填料、没有液体喷淋的状态，此时气相通道完全开放，压降全部来自气体绕流填料几何结构所产生的形阻与摩擦。工程上最通用的做法，是把填料层视为一束当量直径固定的不规则通道，用 Ergun 型方程统一描述层流黏性损失与湍流惯性损失两部分贡献。掌握这一关联式，就能在缺少厂商实测曲线时，用手算快速估算干塔压降，并判断操作点是否落在合理区间。
一、Ergun 型方程的基本形式
Ergun 方程最初针对颗粒填充床提出，其核心思想是：单位床层高度的压降由黏性项和惯性项叠加而成。写成通用形式为：
ΔP/Z = A · μ·u / d_p² · (1-ε)²/ε³ + B · ρ·u² / d_p · (1-ε)/ε³式中 ΔP 为填料层压降（Pa），Z 为填料层高度（m），μ 为气相动力黏度（Pa·s），ρ 为气相密度（kg/m³），u 为以空塔截面计的表观气速（m/s），d_p 为填料当量直径（m），ε 为床层空隙率。A、B 为经验常数，对散装颗粒床通常取 A≈150、B≈1.75。对于规整填料或开孔率较高的散装填料，直接套用颗粒床常数会明显高估压降，因此工程上更常用以填料因子 F_p 或比表面积 a 表达的修正形式。
把 d_p 用比表面积 a（m²/m³）替换，利用 d_p ≈ 6(1-ε)/a 的关系，可得到更贴近填料实际的表达式：
ΔP/Z = A' · μ·u·a²·(1-ε)²/ε³ + B' · ρ·u²·a·(1-ε)/ε³这一形式的好处是：a 和 ε 都是填料样本的标准几何参数，厂商样本中通常直接给出，无需再反算当量直径。实际手算时，常把两项合并写成摩擦因子形式 ΔP/Z = f · (ρu²/2) · (a/ε³)，其中 f 为修正摩擦因子，随 Re 变化。
二、Reynolds 数与流动分区
填料层内的流动状态用颗粒 Reynolds 数判别：
Re_p = ρ·u·d_p / μ （或 Re = ρ·u / (a·μ)）不同 Re 区间，黏性项与惯性项的权重不同，摩擦因子的变化规律也不同。可以划分为三个区：
- 层流区（Re_p < 约 10）：黏性项主导，ΔP/Z 近似与 u 成正比，摩擦因子随 Re 增大而下降，斜率约为 -1。
- 过渡区（Re_p 约 10～1000）：两项均不可忽略，ΔP/Z 与 u 呈 1～2 次方之间的幂律关系，摩擦因子下降趋缓。散装填料塔的常规操作大多落在此区。
- 湍流区（Re_p > 约 1000）：惯性项主导，ΔP/Z 近似与 u² 成正比，摩擦因子趋于常数。此时提高气速压降上升很快，是液泛预警的重要信号。
三、摩擦因子与填料几何参数的关联
把 Ergun 方程无量纲化，可提炼出摩擦因子 f 与 Re 的关系。以 f = (ΔP/Z)·d_p/(ρu²) 定义时，方程变为 f = A/Re_p + B，即一条双曲线叠加常数的形式。这说明：小 Re 时 f 由 A/Re_p 控制，大 Re 时 f 趋近常数 B。对填料而言，A、B 并非普适常数，而与填料的形状因子、开孔率、表面粗糙度相关。规整波纹填料的通道规则、转折少，其 B 值明显低于同比表面积的散装拉西环；而比表面积 a 越大、空隙率 ε 越小，同样气速下压降越高。因此选型时不能只看 a，还要看 ε 与填料因子 F_p——F_p 本质上就是把 a 与 ε 综合成一个压降指示参数。
四、干塔压降手算示例
下面用一个教学示例演示完整手算过程。设某散装填料塔，气相为常温常压空气，ρ = 1.2 kg/m³，μ = 1.8×10⁻⁵ Pa·s；填料比表面积 a = 200 m²/m³，空隙率 ε = 0.72，床层高度 Z = 3 m；表观气速 u = 0.8 m/s。取修正常数 A′ = 150、B′ = 1.75（示例取值，实际应按填料类型选取）。
第一步，算当量直径：d_p ≈ 6(1-ε)/a = 6×0.28/200 = 0.0084 m。第二步，算 Reynolds 数：Re_p = 1.2×0.8×0.0084 / 1.8×10⁻⁵ ≈ 448，落在过渡区，两项都要保留。第三步，分别计算两项：
黏性项 = 150 × 1.8e-5 × 0.8 × 200² × 0.28² / 0.72³
≈ 150 × 1.8e-5 × 0.8 × 40000 × 0.0784 / 0.373
≈ 227 Pa/m
惯性项 = 1.75 × 1.2 × 0.8² × 200 × 0.28 / 0.72³
≈ 1.75 × 1.2 × 0.64 × 200 × 0.28 / 0.373
≈ 202 Pa/m
ΔP/Z ≈ 227 + 202 = 429 Pa/m
ΔP ≈ 429 × 3 ≈ 1287 Pa ≈ 1.29 kPa结果约 1.3 kPa，两项贡献接近，说明确实处于过渡区，不能只留一项。若把气速提高到 1.6 m/s，Re_p 翻倍至约 896，仍在过渡区但更靠近湍流，惯性项按 u² 增长到约 808 Pa/m，黏性项只翻倍到约 454 Pa/m，总压降约 3.8 kPa，增幅远超气速增幅。这正是操作中要警惕的“压降陡升”现象。
五、手算练习建议
建议按以下顺序练习，形成稳定的计算习惯：先由 a、ε 反算 d_p，再算 Re_p 定区；然后按分区决定保留项，逐项代入并检查量纲；最后用 ΔP/Z 乘以 Z 得总压降，并与同类型填料的经验范围对照。练习时可改变单一变量做敏感性分析——例如固定其他参数只改 ε，观察压降如何随空隙率下降而快速上升，这比死记公式更能建立对填料几何参数与压降关系的直觉。掌握这一关联式后，后续引入喷淋液后的湿塔压降与载点、泛点判断，就有了可靠的干塔基准。