📚 课程目录第3章管式反应器设计方程推导

3. 管式反应器设计方程推导

管式反应器（PFR）的设计方程不是凭空给出的公式，它来自对反应器内一个无限小体积微元所做的物料衡算。理解这条推导线索，比记住最终公式更重要——因为空时、空速、转化率之间的关联，以及方程的适用条件，都藏在这条推导线索里。

一、微元物料衡算的建立

先约定参照系：取物料流动方向为轴向正方向，反应器入口为 z=0，出口为 z=L；反应物 A 的摩尔流量记为 F_A，转化率 X_A 定义为已反应的 A 量占进入反应器 A 量的分率，因此入口处 X_A=0，沿轴向单调增大。这一约定在后续公式与图示中保持一致。

在轴向位置 z 处截取长度 dz、截面积为 A_c 的微元。对组分 A 做稳态物料衡算，进入量减去离开量再减去反应消耗量应等于零（稳态下微元内无累积）：

F_A|z  −  F_A|(z+dz)  −  (−r_A)·A_c·dz  =  0

其中 (−r_A) 为 A 的消耗速率（mol·m⁻³·s⁻¹），取正值。把 F_A|(z+dz) 在 z 处展开并保留一阶项，得到 F_A|(z+dz) ≈ F_A|z + (dF_A/dz)dz，代入上式消去同类项：

− dF_A/dz  =  (−r_A)·A_c

这就是 PFR 的微分形式物料衡算。它说明：沿轴向摩尔流量的下降速率，正比于当地反应速率与截面积。

二、引入转化率，得到设计方程

把 F_A 用入口流量 F_A0 与转化率表示：F_A = F_A0(1 − X_A)，于是 dF_A = −F_A0dX_A。代入微分衡算：

F_A0 · dX_A/dz  =  (−r_A)·A_c

沿整个反应器长度从 z=0（X_A=0）积分到 z=L（X_A=X_A,out），并注意 A_c·L 即反应器体积 V：

V / F_A0  =  ∫₀^(X_A,out)  dX_A / (−r_A)

这就是 PFR 的积分形式设计方程。它的物理含义非常清晰：单位入口摩尔流量所需要的反应器体积，等于转化率从 0 增到出口值过程中 1/(−r_A) 对 X_A 的积分面积。反应速率越低（1/(−r_A) 越大），积分面积越大，所需体积越大——这正是 PFR 设计方程最实用的判读方式。

三、空时与空速：把体积换算成操作强度

设计方程给出的是 V/F_A0，工程上更常用空时 τ 描述反应器负荷。空时定义为反应器体积与入口体积流量之比：

τ = V / v_0        （单位：s 或 min）
空速 SV = 1/τ = v_0 / V   （单位：s⁻¹ 或 h⁻¹）

注意 τ 与物料在反应器内的真实停留时间 t 并不总是相等：当反应伴随摩尔数变化（如 A→2B）或温度、压力沿轴向变化时，体积流量 v 会改变，真实停留时间需按 t = ∫dV/v 计算。只有在恒容、恒温、恒压且无摩尔数变化的条件下，τ 才等于真实停留时间。这是初学者最容易混淆的一点。

把 F_A0 = C_A0·v₀ 代入设计方程，可得用空时表达的等价形式：

τ = C_A0 · ∫₀^(X_A,out) dX_A / (−r_A)

可见空时 τ 与入口浓度 C_A0、目标转化率 X_A,out 以及速率表达式共同决定。空速越大，意味着同样体积下处理量越大，但达到同一转化率所需的 τ 不变——空速与转化率之间不是简单的线性关系，而是通过速率积分耦合在一起。

四、示例：一级反应下的定量计算

设液相一级反应 A→P，速率 (−r_A)=kC_A，恒容下 C_A=C_A0(1−X_A)，k=0.05 s⁻¹，要求 X_A,out=0.90。以下为教学示例参数，用于演示计算流程，非实测数据。

∫₀^0.9 dX_A / [k(1−X_A)] = (1/k)·ln[1/(1−0.9)]
= (1/0.05)·ln(10) ≈ 20 × 2.303 ≈ 46.05 s

即 τ = 46.05 s，若 v_0 = 2 L/s，则 V = τ·v_0 ≈ 92.1 L

若把目标转化率提高到 0.99，则 τ = (1/0.05)·ln(100) ≈ 92.1 s，体积翻倍。这一结果直观说明：高转化率区段 1/(−r_A) 急剧上升，积分面积迅速增大，追求最后几个百分点的转化率代价很高。这正是设计方程几何判读的实际价值。

五、方程的适用条件与常见误用

PFR 设计方程的推导隐含了若干前提，使用时必须核对：

稳态操作：微元内无累积项，开停车、负荷突变阶段不适用。
理想平推流：同一截面上浓度、温度均匀，无轴向返混与径向梯度；实际管式反应器需用停留时间分布或轴向扩散模型修正。
沿轴向连续积分：反应速率随 X_A（以及随温度）变化，必须用当地速率积分，不能拿入口或出口速率取平均代替。
体积流量处理正确：气相变摩尔反应需按 v=v₀(1+εX_A)(T/T₀)(P₀/P) 修正浓度，否则 τ 与真实停留时间混淆，转化率会被高估或低估。
单一反应或指定关键组分：多反应并存时，需对每个组分分别列衡算，不能只用一个转化率方程概括。

常见错误：直接把 τ=V/v₀ 当作真实停留时间代入一级反应积分式。对气相变摩尔反应，正确做法是把 C_A 写成 X_A 的函数后再积分，得到的 τ 表达式会多出与膨胀因子 ε 相关的项。

校核思路：先用设计方程算出理论 τ 与 V，再用停留时间分布实验（脉冲示踪）测得的平均停留时间与方差，判断实际反应器偏离平推流的程度；若方差偏大，说明存在返混，需用轴向扩散模型或串联釜模型重新核算转化率，而不是直接套用理想 PFR 公式。

把上述推导串起来看：微元衡算给出微分方程，转化率定义把它变成可积分形式，积分得到 V/F_A0 与 1/(−r_A) 的面积关系，再用空时把体积换算成操作强度。每一步都在回答同一个问题——要达到指定转化率，反应器需要多大、负荷能开多高。掌握这条主线，后续的变温操作、非理想流动修正和误差校核才有稳固的出发点。