2. 间歇反应器转化率与反应时间计算

问题设定与物料衡算方程

间歇反应器（BR）在投料后密闭运行，无进料也无出料，因此对反应物 A 做组分衡算时，累积项等于反应消耗项。设反应器内液相体积 V 恒定，A 的初始浓度为 C_A0，转化率定义为 X = (C_A0 − C_A)/C_A0，则衡算方程可写为：

−dC_A/dt = (−r_A) = k·C_Aⁿ

用转化率表示，C_A = C_A0(1 − X)，dC_A = −C_A0dX，代入得间歇反应器的设计方程：

t = C_A0 ∫₀^X dX / (−r_A)

本例取一级反应 n = 1，−r_A = k·C_A，积分后得到解析解：

t = −(1/k)·ln(1 − X)

教学示例参数：k = 0.05 min⁻¹（30 ℃ 下的假设值），C_A0 = 2.0 mol/L。注意一级反应的 t 与 C_A0 无关，这是本例需要验证的一个重要结论。

积分求解与时间计算

将 X = 0.90 代入解析式：

t = −(1/0.05)·ln(1 − 0.90)
  = −20·ln(0.10)
  = −20 × (−2.302585)
  = 46.05 min

若目标转化率提高到 X = 0.99，则 t = −20·ln(0.01) = 92.10 min。可见一级反应中，转化率从 90% 提到 99%，反应时间恰好翻倍——因为每提升一个数量级的剩余反应物浓度，都需要相同的 ln10/k 时间增量。

下面用数值积分（梯形法）独立校核，避免只依赖解析式：

import math

k = 0.05          # min^-1
CA0 = 2.0         # mol/L
X_target = 0.90

# 解析解
t_analytic = -math.log(1 - X_target) / k

# 数值积分：t = CA0 * ∫ dX / (k*CA0*(1-X)) = ∫ dX / (k*(1-X))
N = 100000
dX = X_target / N
integral = 0.0
for i in range(N):
    X1 = i * dX
    X2 = (i + 1) * dX
    f1 = 1.0 / (k * (1 - X1))
    f2 = 1.0 / (k * (1 - X2))
    integral += 0.5 * (f1 + f2) * dX

print(f"解析解 t = {t_analytic:.3f} min")
print(f"数值解 t = {integral:.3f} min")
print(f"相对偏差 = {abs(integral - t_analytic)/t_analytic*100:.4f}%")

预期输出：解析解 46.052 min，数值解约 46.052 min，相对偏差在 10⁻⁴% 量级（梯形法在 10⁵ 分段下的截断误差）。若偏差明显偏大，应检查被积函数是否写成了 1/(k·C_A0·(1−X)) 而漏乘 C_A0 的抵消关系。

转化率随时间的变化趋势

由 X = 1 − exp(−kt) 可知，转化率随时间呈指数趋近 1 的曲线：初期反应物浓度高、速率快，X 上升陡；后期 C_A 降低，速率变慢，X 上升趋缓。下表给出若干时刻的示例计算值（k = 0.05 min⁻¹）：

工程判断：若工艺要求 X ≥ 0.90，则单批反应时间至少 46 min；再考虑装料、卸料、清洗等辅助时间（假设共 30 min），单批周期约 76 min，据此可估算日批次数与产能。追求 99% 转化率会使反应时间翻倍而产能近乎减半，需在转化率收益与时间成本之间权衡。

结果校核与工程判读

• 量纲校核：t = −ln(1−X)/k，k 的单位 min⁻¹，故 t 单位为 min，正确。
• 极限校核：X → 0 时 t → 0；X → 1 时 t → ∞，符合一级反应理论上无法达到 100% 转化率的规律。
• 与 C_A0 无关性：把 C_A0 从 2.0 改为 5.0 mol/L，解析式结果仍为 46.05 min，验证一级反应时间与初始浓度无关；若换成二级反应（−r_A = kC_A²），则 t = X/[kC_A0(1−X)]，此时 C_A0 会显著影响时间。
• 数值与解析一致性：两者偏差小于 0.01%，说明积分区间划分足够细，可放心用于工艺时间估算。