› Ergun 式 模拟器 返回
粉体·流动层模拟器

Ergun 式 模拟器 — 填充床压力损失

使用 Ergun 式 ΔP/L = 150μu(1−ε)²/(ε³d_p²) + 1.75ρu²(1−ε)/(ε³d_p) 实时计算球形颗粒填充床中流体通过时的单位长度压力损失。从空塔速度 u、颗粒径 d_p、空隙率 ε 和粘度 μ 计算 ΔP/L、颗粒雷诺数 Re_p、最小流动化速度 u_mf 和层流贡献比例，通过填充床侧面示意图和 ΔP/L − u 对数曲线可视化从固定床到流动化的行为。

参数设置
空塔速度 u
mm/s
颗粒径 d_p
mm
空隙率 ε
动粘性系数 μ
×10⁻³ Pa·s

默认参数：u=10 mm/s、d_p=2.0 mm、ε=0.40、μ=1.0×10⁻³ Pa·s（20 °C 水相当）、流体密度 ρ=1000 kg/m³ 固定。最小流动化速度 u_mf 以颗粒密度 ρ_p=2500 kg/m³（砂粒相当）计算。

计算结果
—
ΔP/L 压力梯度
—
颗粒 Re_p
—
最小流动化速度 u_mf
—
层流贡献（第 1 项）
填充床侧面示意图

圆柱形填充床的侧面示意。多个球形颗粒（直径 d_p）以空隙率 ε 填充，流体以空塔速度 u 从下向上通过。颗粒大小和密度与当前的 d_p 和 ε 值关联。

ΔP/L − u 曲线（对数坐标）

横轴为 u（mm/s、log10）、纵轴为 ΔP/L（Pa/m、log10）。低速处层流（斜率 1）、高速处湍流（斜率 2）的过渡。黄色标记为当前的 (u, ΔP/L)，虚线为各极限情况（层流/惯性）的渐近线。

理论和主要公式

球形颗粒填充床中流体通过时的单位长度压力损失由层流贡献和惯性贡献的总和给出：

$$\frac{\Delta P}{L} = \frac{150\,\mu\,u\,(1-\varepsilon)^2}{\varepsilon^3\,d_p^2} + \frac{1.75\,\rho\,u^2\,(1-\varepsilon)}{\varepsilon^3\,d_p}$$

其中 $u$ 为空塔速度（m/s）、$d_p$ 为球形颗粒径（m）、$\varepsilon$ 为空隙率、$\mu$ 为流体粘度（Pa·s）、$\rho$ 为流体密度（kg/m³）。第 1 项为 Carman-Kozeny 流（层流），第 2 项为 Burke-Plummer 流（惯性）贡献。颗粒雷诺数为：

$$\mathrm{Re}_p = \frac{\rho\,u\,d_p}{\mu\,(1-\varepsilon)}$$

层流近似下的最小流动化速度（颗粒密度 $\rho_p$）仅由第 1 项得出：

$$u_{mf} = \frac{\varepsilon^3\,d_p^2\,(\rho_p-\rho)\,g}{150\,\mu\,(1-\varepsilon)}$$

本工具以 $\rho = 1000$ kg/m³、$\rho_p = 2500$ kg/m³、$g = 9.81$ m/s² 进行固定计算。

Ergun 式模拟器说明

🙋
我在教科书中看到过"填充床压力损失"这个概念，但还是不太理解。这是什么场景下的问题？
🎓
一个好的例子是咖啡滴漏过滤器或过滤装置。当你在粉或颗粒层上方倒水时，水通过粒子层需要克服"压力"。在化工厂中，常见应用包括触媒填充塔、吸附塔、干燥塔和过滤层。在土木工程中，砾石过滤器和玻璃珠试验柱也是代表例。Ergun 式用来计算流体通过填充层需要的"单位长度压力 ΔP/L"，它取决于空塔速度 u、颗粒径 d_p、空隙率 ε、流体粘度 μ 和密度 ρ。在我们工具的默认参数下（水在 20 °C、d_p=2 mm、ε=0.40、u=10 mm/s），ΔP/L 约为 2929 Pa/m，也就是说要让 1 m 的填充层有流体通过，需要约 0.03 气压的压力。
🙋
公式中有两项，它们分别代表什么？
🎓
第一项是层流贡献（Carman-Kozeny），在流速缓慢时起主导作用。它将颗粒间的缝隙看作细小管道，然后用 Hagen-Poiseuille 流的叠加来推导。第二项是惯性贡献（Burke-Plummer），在流速较快时起主导作用。此时流体在孔隙内流动变得湍流化，颗粒对流体的碰撞损失变成主要因素。两项通过颗粒雷诺数 Re_p = ρu d_p/(μ(1−ε)) 来竞争。在我们工具的默认参数下，Re_p ≈ 33.3，处在"层流贡献约 70%、惯性贡献约 30%"的过渡区。当你把 u 改到 1 mm/s 时，层流贡献上升到 95% 以上；改到 100 mm/s 时，惯性贡献就变得支配性了。
🙋
右下方的对数图上，低速和高速区的斜率不同，这是为什么？
🎓
在对数坐标上，第一项与 u 成正比，所以斜率是 1；第二项与 u² 成正比，所以斜率是 2。低速时绿色虚线（第一项）主导，高速时红色虚线（第二项）主导。我们工具的右下图清楚地显示了这种过渡，黄色标记点指示当前的 u 所对应的位置。
🙋
"最小流动化速度 u_mf"是什么意思？和流化床有关吗？
🎓
完全正确。u_mf 是指当上升流体的阻力等于颗粒层见清自重的那个空塔速度，超过这个速度后颗粒层就会"浮起来"进入流动化状态。在层流近似下，我们只用 Ergun 式的第一项来估算 u_mf = ε³ d_p² (ρ_p − ρ)g / [150μ(1−ε)]（颗粒密度取 ρ_p = 2500 kg/m³，模拟砂粒）。在我们工具的默认参数下，u_mf ≈ 41.9 mm/s，所以当前的操作流速 10 mm/s 还在固定床状态。在高 Re_p 情况下需要用 Wen-Yu 或 Grace 相关式来补正。
🙋
改变空隙率 ε 时，ΔP/L 变化很大，这是为什么？
🎓
因为 ε³ 出现在分母中。当 ε 从 0.40 降到 0.30 时，ε³ 从 0.064 降到 0.027，约降低了 2.4 倍，所以 ΔP/L 会上升 2 到 3 倍。反过来，当 ε 升到 0.50 时，ε³ 从 0.064 升到 0.125，ΔP/L 约降半。在实际操作中，"颗粒球形度低、微粉混入、偏析发生"都会导致局部 ε 下降，引起意想不到的压力上升。设计时通常会对 ε 的不确定性（±0.03）留出安全系数。你可以用我们工具的滑块将 ε 在 0.30 到 0.60 范围内变化，感受一下 ΔP/L 对 ε 的敏感程度。

常见问题

Ergun 式是一个半经验式，用于计算球形颗粒填充床中流体通过时的单位长度压力损失 ΔP/L。公式为 ΔP/L = 150μu(1−ε)²/(ε³d_p²) + 1.75ρu²(1−ε)/(ε³d_p)，其中第 1 项为 Carman-Kozeny 层流贡献，第 2 项为 Burke-Plummer 惯性贡献。在本工具的默认参数下（u=10 mm/s、d_p=2 mm、ε=0.40、μ=1×10⁻³ Pa·s、ρ=1000 kg/m³），ΔP/L ≈ 2929 Pa/m、Re_p ≈ 33.3、层流贡献率约 72.0%。
填充床的颗粒雷诺数定义为 Re_p = ρ u d_p / (μ (1−ε))。当 Re_p < 10 时为层流（第 1 项支配），当 Re_p > 1000 时为完全湍流（第 2 项支配），介于两者之间为过渡域，两项同程度贡献。本工具默认参数下 Re_p ≈ 33.3，位于过渡域，层流贡献约 72%、惯性贡献约 28%。当 u 降至 1 mm/s 时，Re_p ≈ 3.3；当 u 升至 100 mm/s 时，Re_p ≈ 333。
最小流动化速度 u_mf 是指上升流体的阻力与颗粒层表观自重达到平衡时的空塔速度。超过此值时，颗粒层进入流动化状态。在层流近似下，仅用 Ergun 式第 1 项计算：u_mf = ε³ d_p² (ρ_p − ρ) g / [150 μ (1−ε)]（颗粒密度 ρ_p = 2500 kg/m³）。本工具默认参数下 u_mf ≈ 41.9 mm/s，当前操作流速 10 mm/s 处于固定床范围。高 Re_p 需用 Wen-Yu 或 Grace 相关式补正。
在低速层流域，Ergun 式第 1 项支配，ΔP/L 与 u 成正比（对数图斜率 1），源自孔隙内粘性流（Hagen-Poiseuille），是 Carman-Kozeny 式的本质。在高速湍流域，第 2 项支配，ΔP/L 与 u² 成正比（斜率 2），源自惯性损失和流体与颗粒碰撞。本工具右下图清楚显示这种过渡，黄色标记表示当前 u 的位置。

实际应用

固定床催化反应器压力设计：在石油精炼的氢化脱硫（HDS）、合成气改质、汽车催化转换器等固定床反应器中，压力损失直接影响泵的动力消耗。例如，若使用 d_p=3 mm 的颗粒催化剂，ε=0.40、u=50 mm/s 运行，本工具计算出 ΔP/L ≈ 8 kPa/m。对于塔高 5 m 的反应器，总压损约 40 kPa（约 0.4 气压），是压缩机的主要能耗。虽然将颗粒径增大到 5 mm 可以使 ΔP/L 降至原来的约 1/2.8，但会恶化内部物质传递，降低反应速率，所以设计时必须权衡"压损 vs 反应性"。

吸附·干燥塔和穿透曲线：空气干燥（分子筛、硅胶）和 VOC 除去用活性炭塔通常在 u=100～300 mm/s 下运行。将 d_p=3 mm、ε=0.38、u=200 mm/s 输入本工具，得 ΔP/L ≈ 30 kPa/m，此时 Re_p ≈ 600 处于完全湍流域。实际操作中需要同时评估穿透时间、物质转移区域宽度和压损，并据此选定塔径和颗粒径。由于微粉混入会导致局部 ε 下降，ΔP/L 急剧上升，所以再生时的颗粒破碎也是寿命设计的关键。

流动床反应器启动和 u_mf：流化接触分解（FCC）、气化炉、流动床锅炉等在超过 u_mf 的速度下运行。将 d_p=0.5 mm、ε=0.45、μ=2×10⁻⁵ Pa·s（高温空气）代入本工具，可得 u_mf 为数厘米/秒量级，实际运行速度通常为其 3～10 倍。层流近似的 u_mf 是出发点，实际操作需根据 Geldart 颗粒分类（A、B、C、D）应用相应的相关式（Wen-Yu、Grace、Saxena-Vogel）进行补正。本工具对于理解固定床到流动化过渡边界在颗粒径和空隙率变化下的敏感性特别有用。

砾石过滤和地下水净化：上水处理厂的快速过滤池（砂层 d_p=0.5～1 mm）、地下水净化用砾石过滤、雨水渗透设施中的碎石层等，Ergun 式（或 Forchheimer 式）都是基础方程。将 d_p=0.8 mm、ε=0.40、u=2 mm/s（典型过滤速度）代入，得 ΔP/L ≈ 1～2 kPa/m，此时 Re_p ≈ 4 处于层流域。反洗时将 u 提高到 50 mm/s 以上超过 u_mf，使砂层膨胀消除堵塞。使用本工具可以确认设计流速与 u_mf 的比值，从而估算反洗泵的规格。

常见误解和注意事项

最常见的误解是"Ergun 式对所有填充床都适用"。Ergun 式是在"球形颗粒、窄颗粒径分布、ε 在 0.35～0.55、壁面影响小"等条件下的半经验式，对于针状、薄片状催化剂或颗粒径分布宽的砂体，以及颗粒直径与塔径比（D_t/d_p）< 10 的情况，精度会明显下降。对于球形度 φ_s < 1 的颗粒，应使用 Sauter 平均径 d_{32}、形状补正项（如 Carman 的球形度修正），或切换到 Forchheimer 式或广义 Carman-Kozeny 式。本工具的计算基于"球形、单分散、连续体假设"的理想条件，实际应用时应留出 20～30% 的不确定性余地。

第二个常见误解是"空隙率 ε 用一阶补正就够了"。因为 ε³ 在分母中，Ergun 式对 ε 的敏感性很强：ε=0.40 ± 0.02 的不确定性就会导致 ΔP/L 有 ±15～20% 的误差。实际操作中，充填方式（自然倒入、振动、蒸汽抽出）可以使 ε 变化超过 0.05，因此要分别测量初始空隙率和运行时空隙率。使用本工具的滑块将 ε 在 0.38、0.40、0.42 之间变化，体感 ΔP/L 的敏感性后，再决定安全系数的取值。

第三个误解是"层流近似的 u_mf 总是偏于安全一侧"。层流近似公式 u_mf = ε³ d_p² (ρ_p − ρ)g / [150μ(1−ε)] 仅在 Re_p,mf < 10 时有效。对于 d_p > 1 mm 的大颗粒（Geldart B、D 分类），Re_p,mf 会超过 10，本工具计算的 u_mf 会高估（实际 u_mf 应该从层流项+惯性项的交点求取）。推荐使用 Wen-Yu 相关式 Re_p,mf = √(33.7² + 0.0408·Ar) − 33.7（Ar 为阿基米德数），本工具的 u_mf 应视为"层流近似出发点"。对 Geldart D 颗粒（d_p 数毫米）尤其要注意补正。

相关工具

HVAC 矩形风管尺寸计算器 — 压力损失计算 HVAC 矩形风管的尺寸和压力损失计算模拟器。从风量、宽度、高度和风管长度计算水力直径 D_h、平均风速 V、摩擦系数 f、单位长度压力损失... 管道内压力损失计算机（Darcy-Weisbach） 使用 Darcy-Weisbach 式的管道压力损失计算工具。自动判断层流/湍流，计算摩擦系数和弯头、阀门等局部损失。用于管道设计和流体系统... 管路压力损失计算机 免费在线管路压力损失计算工具。使用 Darcy-Weisbach 式和 Colebrook 式精确计算，考虑管件和阀门损失。实时可视化运行状态，直观展示... 过滤压力·滤料阻力 计算机 NovaSolver 免费在线计算工具，计算过滤压力损失和滤料阻力。基于 Darcy 定律和 Ruth 过滤方程，支持填充层、纤维、膜、颗粒等各种媒体... 管内流动 模拟器 NovaSolver 免费在线计算工具模拟管内流动。计算内径 50 mm、长度 100 m 的钢管中水以 3 m/s 流速流动时的压力损失。考虑管件和阀门...

使用指南

  1. 使用滑块设置空塔速度 u（mm/s）。触媒填充床通常在 0.05～0.5 m/s 范围内。
  2. 输入颗粒径 d_p（mm）。硅胶吸附塔通常为 3～5 mm，FCC 触媒再生器通常为 50～150 μm。
  3. 设置空隙率 ε（0～1）和流体粘度 μ×10⁻³ Pa·s，即时显示压力梯度 ΔP/L、颗粒雷诺数 Re_p 和最小流动化速度 u_mf。

具体计算示例

玻璃珠填充床：d_p=2 mm、ε=0.40、u=0.15 m/s、μ=1.2×10⁻³ Pa·s（水溶液），计算结果为 ΔP/L=580 Pa/m、Re_p=250、层流贡献率 65%。而 FCC 再生器中焦炭燃烧情况：d_p=80 μm、ε=0.45、u=0.8 m/s、μ=2.0×10⁻⁵ Pa·s（110 °C 空气），ΔP/L=18500 Pa/m、Re_p=0.32，此时湍流项占支配地位，最小流动化速度 u_mf=0.21 m/s。

实务中的注意事项