› 埃尔根方程模拟器 返回
颗粒与流化床模拟器

埃尔根方程模拟器 — 填充床的压降

基于埃尔根方程 ΔP/L = 150μu(1−ε)²/(ε³d_p²) + 1.75ρu²(1−ε)/(ε³d_p) 实时计算流体通过球形颗粒填充床时的单位长度压降。从空塔流速 u、颗粒直径 d_p、空隙率 ε、粘度 μ 求 ΔP/L、颗粒雷诺数 Re_p、最小流化速度 u_mf 及层流贡献占比，并以填充床侧视示意图与 log-log Δ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 曲线（log-log）

横轴 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² 计算。

什么是埃尔根方程模拟器

🙋
"填充床的压降"这个词在教科书上看过，但还是没什么概念。具体是什么场景？
🎓
最直观的例子是咖啡滤杯或砂滤器：水从颗粒层上方流下时，要克服压力才能穿过缝隙。化工里有催化剂填充塔、吸附塔、干燥塔、过滤层；土木里有砂砾过滤器、玻璃珠试验柱。埃尔根方程就是用空塔流速 u、颗粒直径 d_p、空隙率 ε、流体粘度 μ 与密度 ρ 估算"通过床所需的单位长度压力 ΔP/L"的半经验公式。本工具默认值（20 °C 水、d_p=2 mm、ε=0.40、u=10 mm/s）下 ΔP/L ≈ 2.93 kPa/m，相当于每米床约 0.03 个大气压。
🙋
公式里有两项，它们各自代表什么？
🎓
第 1 项是层流贡献（Carman-Kozeny），慢流时占主导，可以把颗粒间隙看成细管，由 Hagen-Poiseuille 流叠加得出。第 2 项是惯性贡献（Burke-Plummer），快流时占主导：流体穿过缝隙后撞到下一颗颗粒，损失动能。两者的强弱由颗粒雷诺数 Re_p = ρu d_p/(μ(1−ε)) 决定。本工具默认值下 Re_p ≈ 33.3，处于过渡区——层流约 7 成、惯性约 3 成。把 u 调到 1 mm/s，层流占 95% 以上；调到 100 mm/s，惯性占主导。
🙋
右下方 log-log 图里，低速和高速的斜率不同，为什么？
🎓
log-log 看，第 1 项 ΔP/L ∝ u，斜率 1；第 2 项 ΔP/L ∝ u²，斜率 2。低速时埃尔根曲线贴近绿色虚线（第 1 项），高速时贴近红色虚线（第 2 项）。设计实务先确认目标 Re_p 落在哪个区——层流要做粘度修正、湍流要做颗粒撞击损失修正，再决定操作条件。比如反应器催化剂层在高速下运行时 ΔP/L 会按二次方暴增，泵的能耗也指数式飙升，需特别注意。
🙋
"最小流化速度 u_mf"是什么？是不是和流化床有关？
🎓
没错，它就是流化床的入口指标。当上行流体对颗粒的曳力刚好等于颗粒视重（重力减浮力）时，颗粒床被"轻轻抬起"开始流化，此时的空塔流速即 u_mf。在仅取埃尔根方程第 1 项的层流近似下，u_mf = ε³ d_p² (ρ_p−ρ)g / [150μ(1−ε)]。本工具固定砂粒（ρ_p=2500 kg/m³）计算，默认值下 u_mf ≈ 41.9 mm/s——所以默认 u=10 mm/s 是固定床；u=50 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.125，ΔP/L 几乎减半。实务中"颗粒球度差、混入细粉、偏析"会使局部 ε 下降，导致差压意外升高。设计时通常预留 ε 的不确定度 ±0.03 作安全裕度。本工具中把 ε 在 0.30〜0.60 间扫一下，即可直观感受 ΔP/L 的灵敏度。

常见问题

埃尔根方程（Ergun equation）是用空塔流速 u、颗粒直径 d_p、空隙率 ε、流体粘度 μ 与密度 ρ 估算流体通过球形颗粒填充床时单位长度压降 ΔP/L 的半经验公式。ΔP/L = 150μu(1−ε)²/(ε³d_p²) + 1.75ρu²(1−ε)/(ε³d_p)。第 1 项为 Carman-Kozeny 层流贡献，第 2 项为 Burke-Plummer 惯性贡献。本工具默认值下 ΔP/L ≈ 2.93 kPa/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；调到 100 mm/s，Re_p ≈ 333。
u_mf 是上行流体曳力恰等于颗粒床表观重量时的空塔流速，超过该值颗粒床即开始流化。层流近似下 u_mf = ε³ d_p² (ρ_p − ρ) g / [150 μ (1−ε)]，其中颗粒密度 ρ_p = 2500 kg/m³。本工具默认值下 u_mf ≈ 41.9 mm/s，故当前 u=10 mm/s 仍处于固定床区域。Re_p 较高时常用 Wen-Yu 或 Grace 关联式。
低速层流区下，第 1 项支配，ΔP/L ∝ u（log-log 斜率 1），由颗粒间隙的粘性流（Hagen-Poiseuille 流）所致，这是 Carman-Kozeny 的本质。高速湍流区下，第 2 项支配，ΔP/L ∝ u²（斜率 2），由惯性损失与射流撞击主导，由 Burke-Plummer 公式描述。本工具右下方 log-log 图清晰呈现这一过渡，黄色标记为当前 u。

实际应用

固定床催化反应器的差压设计：石油精炼的加氢脱硫（HDS）、合成气重整、汽车尾气催化转化器等固定床反应器，其差压直接决定运行泵的能耗。例如在 d_p=3 mm 颗粒催化剂、ε=0.40、u=50 mm/s 条件下，本工具给出 ΔP/L ≈ 8 kPa/m，对 5 m 床而言全压损约 40 kPa（约 0.4 atm），是压缩机能耗的主要来源。把 d_p 从 3 mm 增大到 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（完全湍流区）。实际装置同时评估穿透时间、传质区长度与压降，决定塔径与颗粒粒径选择。混入细粉（fines）局部 ε 下降会使 ΔP/L 急升，再生工况下颗粒破碎也是寿命设计关键。

流化床反应器的启动与 u_mf：流化催化裂化（FCC）、气化炉、流化床锅炉等装置在 u_mf 之上运行。代入 d_p=0.5 mm、ε=0.45、μ=2×10⁻⁵ Pa·s（高温空气），u_mf 在数 cm/s 量级，实际操作速度约为其 3〜10 倍。埃尔根方程的层流近似 u_mf 是设计起点，实际中按 Geldart 颗粒分类（A、B、C、D）选用相应关联式（Wen-Yu、Grace、Saxena-Vogel）作修正。本工具最适合通过粒径、空隙率扫描直观感受"固定床→流化"边界的灵敏度。

砂砾过滤器与地下水净化：城市供水的快速过滤池（砂层 d_p=0.5〜1 mm）、地下水净化的砂砾过滤器、雨水渗透设施的碎石层等土木场景，埃尔根方程（或 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 之比，可估算反洗泵规格。

常见误解与注意点

最常见的误解是"埃尔根方程对所有填充床都适用"。该方程是在"近球形颗粒、窄粒度分布、ε=0.35〜0.55、壁效应小"这一范围内导出的半经验公式，对针状／片状催化剂、粒度分布宽的砂、粒径柱径比（D_t/d_p）<10 等情况精度下降。球度 φ_s < 1 的颗粒应使用 Sauter 平均粒径 d_32、引入形状修正项（d_p → φ_s·d_p），或改用 Forchheimer 方程或非牛顿广义形式。本工具基于"球形、单分散、连续介质假设"的理想条件给值，实际数据通常需预留 20〜30% 的不确定度。

其次是"空隙率 ε 仅作零阶修正即可"的误解。由于 ε³ 在分母里，埃尔根方程对 ε 极敏感：ε=0.40 ± 0.02 的不确定可换算成 ΔP/L 的 ±15〜20% 误差。实际装填操作（降装、振动、排汽）会让 ε 变动 0.05 以上，初始空隙率与运行空隙率必须分开测量。本工具中将 ε 在 0.38、0.40、0.42 间扫一下感受灵敏度后再决定安全系数，是稳妥做法。

最后是"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 为阿基米德数）。Geldart D 颗粒（d_p > 数 mm）下尤其需要修正。本工具的 u_mf 值仅作"层流近似的起点"使用。

相关工具

HVAC 风管尺寸模拟器 — 矩形风管的压损计算 HVAC 矩形风管尺寸与压损模拟器：由风量、宽度、高度、风管长度即时计算水力直径 D_h、平均风速 V、摩擦系数 f、单位长压降 ΔP/L 与总压降 ΔP，并以 Swamee-Jai... 管道压降计算器（达西-韦斯巴赫） 使用达西-韦斯巴赫公式在线计算管道压降。输入管径、长度、流量与流体属性，实时获取雷诺数、摩擦系数、沿程与局部损失（弯头/阀门）。工具在莫迪图上动态标注工作点，帮助工程师快速分析流体系... 阀门压降计算工具 阀门压降计算工具：基于Kv流量系数，在线计算闸阀、截止阀、球阀、蝶阀、止回阀的压降。支持多阀门串联，自动合成Kv与总压降，提供流速预警。适用于管道系统与流体工程设计的免费工具。 管道压降计算器 使用Darcy-Weisbach方程与Colebrook-White公式精确计算管道压降。本工具集成K系数法叠加局部损失，支持实时莫迪图标注工作点，适用于湍流区摩擦系数的迭代求解，为... 过滤压力与滤材阻力计算器 使用达西定律、科兹尼-卡曼方程与Ruth过滤方程，实时计算过滤压力损失、透过率及滤饼比阻。本计算器提供直观的滤饼积累动画，帮助工程师深入理解并优化固液分离过程，适用于过滤工艺设计、性... 管道流动压力损失计算器（Moody图） 管道流动压力损失计算器：基于Colebrook-White方程与Moody图，在线计算雷诺数、摩擦系数和压力损失。输入管径、流量等参数，实时获取结果并可视化ΔP-Q曲线。免费使用，助...