2. Langmuir-Hinshelwood动力学在固定床反应器中的实现
反应体系与速率模型设定
本例以气相甲醇脱水制二甲醚为教学对象，在固定床反应器上实现 Langmuir-Hinshelwood（L-H）型速率表达式。反应式取 2 CH3OH → CH3OCH3 + H2O，催化剂为 γ-Al₂O₃。L-H 机理假设甲醇在催化剂表面单点吸附，表面反应为速率控制步骤，二甲醚与水竞争吸附较弱，可忽略其吸附项，于是速率式写成：
r = k·K_M·p_M / (1 + K_M·p_M)^2 [kmol/(kg_cat·s)]
k = A·exp(-Ea/(R·T)) Arrhenius 表面反应速率常数
K_M = K0·exp(-ΔH_ads/(R·T)) 甲醇吸附平衡常数（放热吸附，ΔH_ads < 0）关键点在于：k 与 K_M 都随温度变化，但方向相反——升温使 k 增大，却使 K_M 减小。因此 L-H 速率对温度常出现非单调行为，这正是本例要考察吸附常数影响的原因。Aspen Plus 中 L-H 型速率不能直接用幂律形式录入，需在反应动力学中选用 Langmuir-Hinshelwood 速率类型，或通过 General 表达式自定义吸附项。
在 Aspen Plus 中构建 L-H 速率
操作路径（以 Aspen Plus V12 的常规界面逻辑描述，不同版本菜单文字可能略有差异）：
- 物性方法选
IDEAL（本例为气相中压体系，教学简化；物性方法匹配问题在后续案例单独讨论）。 - 定义组分：
CH3OH、CH3OCH3、H2O。 - 在 Reactions → Kinetic 中新建反应，反应类型选 LHHW（Langmuir-Hinshelwood-Hougen-Watson）。
- 在 Kinetic 页签填写 Arrhenius 参数：指前因子 A、活化能 Ea；在 Adsorption 页签为甲醇添加吸附项，填写 K0 与 ΔH_ads，吸附形式选单点吸附、指数 1。
- 驱动力表达式按化学计量写为
p_M^2 - p_DME·p_H2O/K_eq的简化形式；教学示例中令平衡项可忽略，驱动力取p_M^2，吸附项取(1+K_M·p_M)^2。
教学示例参数与 RPlug 设置
下表为示例假设参数（非实测值，仅用于演示趋势与操作流程）：
| 参数 | 符号 | 示例值 | 单位 |
|---|---|---|---|
| 指前因子 | A | 3.0×10⁶ | kmol/(kg_cat·s·bar) |
| 活化能 | Ea | 8.5×10⁴ | J/mol |
| 吸附指前因子 | K0 | 2.0×10⁻³ | 1/bar |
| 吸附焓 | ΔH_ads | -4.5×10⁴ | J/mol |
| 反应器长度 | L | 2.0 | m |
| 床层直径 | D | 0.05 | m |
| 床层空隙率 | ε | 0.4 | — |
| 进口温度 | T_in | 523 | K |
| 进口压力 | P | 5.0 | bar |
| 甲醇进料 | F_M | 1.0×10⁻³ | kmol/s |
RPlug 设置：反应器类型选 Kinetic，反应集挂接上面建立的 LHHW 反应；床层按上表填写长度、直径、空隙率；热模式先选 Isothermal（等温）以隔离温度效应，便于单独观察吸附常数的影响；压力沿程取恒定 5 bar。
吸附常数对出口转化率的影响
保持 A、Ea、T、P、进料不变，仅改变 K0（等价于改变吸附强度），在 523 K 下得到如下趋势（示例计算结果，用于说明规律，非实测）：
| K0 (1/bar) | K_M @523K (1/bar) | 出口甲醇转化率 | 现象 |
|---|---|---|---|
| 2.0×10⁻⁴ | ≈0.0065 | ≈0.21 | 吸附弱，表面覆盖度低，速率受吸附控制 |
| 2.0×10⁻³ | ≈0.065 | ≈0.58 | 覆盖度适中，速率较高 |
| 2.0×10⁻² | ≈0.65 | ≈0.47 | 吸附过强，活性位被占满，速率下降 |
规律：转化率随 K_M 先升后降，存在最优吸附强度。原因是速率式分母 (1+K_M·p_M)^2 随 K_M 增大而迅速增大，当 K_M·p_M ≫ 1 时速率近似与 1/K_M 成正比，即「强吸附中毒」。这一非单调性正是 L-H 模型区别于幂律模型的核心特征，也是参数录入后必须做敏感性检查的原因。
结果判读与常见错误
- 非单调是正常现象：若转化率随 K0 单调上升，多半是吸附项符号写反（ΔH_ads 取了正值）或吸附项未进入分母。
- 量纲核对：LHHW 中吸附常数单位须与分压单位一致（本例用 bar），若压力单位设为 Pa 而 K0 按 1/bar 录入，K_M·p_M 会放大约 10⁵ 倍，直接落入强吸附抑制区。
- 等温先验证：先用 Isothermal 模式确认速率式与参数无误，再切换 Adiabatic 观察床层温升；绝热下升温会同时抬高 k、压低 K_M，出口转化率对 K0 的敏感性会被温度效应叠加。
- 与幂律对照：同一组进出口条件下，幂律模型无法复现「吸附过强导致速率下降」的现象，这正是选用 L-H 模型的判据。