数学上的区别在于反应源项的位置和单位。 对于均相反应，反应速率直接出现在微分物种守恒方程中，作为体积源项 (R_{Vi})，单位为 (\mathrm{mol,m^{-3},s^{-1}})。对于发生在表面或相界面的非均相反应，其速率作为边界通量 (R_{Si}) 施加，单位为 (\mathrm{mol,m^{-2},s^{-1}})，或者在采用反应器尺度平均模型时转换为等效体积源项。

均相反应改变整个相体积内的物种浓度；非均相反应改变表面或界面处的物种通量。 正确的数学形式取决于反应是在表面局部解析，还是在中试规模反应器模型中表示为平均源项。

物种衡算如何表示反应

通用守恒方程

对于物种 (i)，摩尔物种衡算可以写为

[ \frac{\partial C_i}{\partial t}

-\nabla \cdot \mathbf{N}i + R{Vi} ]

其中 (C_i) 是物种浓度，(\mathbf{N}i) 是其摩尔通量，而 (R{Vi}) 是每单位体积的净生成速率。

散度项表示进出微分控制体的输运。反应项表示该体积内化学反应的生成或消耗。

均相反应速率

均相反应发生在单一流体相内，例如液相溶液或气相混合物。其速率直接进入控制微分方程：

[ \frac{\partial C_i}{\partial t}

-\nabla \cdot \mathbf{N}i + R{Vi} ]

对于包含 (j=1,\ldots,N_r) 个反应的反应机理，

[ R_{Vi}

\sum_{j=1}^{N_r} \nu_{ij} r_{Vj} ]

其中 (\nu_{ij}) 是物种 (i) 在反应 (j) 中的化学计量系数，而 (r_{Vj}) 是每单位流体体积的反应速率。

典型单位为

[ r_{Vj}, R_{Vi} ;:; \mathrm{mol,m^{-3},s^{-1}} ]

负值表示净消耗，正值表示净生成。

非均相反应如何进入衡算

作为边界条件的表面反应

当反应仅发生在催化剂表面或相边界时，它并不作用于相邻的流体体积内部。相反，它改变了该表面处的物种通量。

在界面 (\Gamma) 处，边界条件可以表示为

[ -\mathbf{n}\cdot\mathbf{N}_i

R_{Si} ]

其中 (\mathbf{n}) 是外法线单位向量，而 (R_{Si}) 是表面生成或消耗速率。

对于多个表面反应，

[ R_{Si}

\sum_{j=1}^{N_s} \nu_{ij} r_{Sj} ]

典型单位为

[ r_{Sj}, R_{Si} ;:; \mathrm{mol,m^{-2},s^{-1}} ]

因此，反应出现在边界条件中，而不是作为主体流体方程中的体积源项。

界面通量跳跃

对于发生在两相界面上的反应，更一般的条件是通量不连续：

[ \mathbf{n}\cdot \left( \mathbf{N}{i,2}-\mathbf{N}{i,1} \right)

R_{Si} ]

这表示离开界面的物种通量与进入界面的物种通量之差等于在界面上生成或消耗的量。

这种形式对于气-液、气-固或液-固中试反应器非常有用，在这些反应器中，两相都有各自独立的输运方程需要求解。

为反应器模型转换基于催化剂的速率

基于每单位催化剂质量的速率

非均相催化动力学通常按每单位催化剂质量给出：

[ r_{i,\mathrm{cat}} ;:; \mathrm{mol,kg_{cat}^{-1},s^{-1}} ]

固定床反应器模型可能需要每单位反应器体积的源项。转换关系为

[ R_{Vi,\mathrm{eff}}

\rho_{b,\mathrm{cat}},r_{i,\mathrm{cat}} ]

其中 (\rho_{b,\mathrm{cat}}) 是每单位反应器体积的催化剂堆积量。

所得的有效源项单位为

[ R_{Vi,\mathrm{eff}} ;:; \mathrm{mol,m^{-3}_{reactor},s^{-1}} ]

这使得非均相反应可以插入一维反应器衡算中，同时保留原始的基于催化剂的动力学。

基于每单位催化剂表面积的速率

如果速率按每单位催化剂表面积表示，

[ r_{i,\mathrm{surf}} ;:; \mathrm{mol,m^{-2}_{cat},s^{-1}} ]

则相应的反应器体积源项为

[ R_{Vi,\mathrm{eff}}

a_s r_{i,\mathrm{surf}} ]

其中 (a_s) 是每单位反应器体积的催化剂表面积，单位为 (\mathrm{m^2_{cat},m^{-3}_{reactor}})。

这一区别很重要：表面速率不能直接添加到体积衡算中，除非乘以上述相关的界面面积密度。

中试规模建模选择

显式非均相模型

显式非均相模型分别对流体相和固相求解独立的守恒方程。流体方程包含输运和相间交换项，而固相或催化剂方程则考虑表面反应、颗粒内扩散或两者兼有。

例如，流体相衡算可以示意性地表示为

[ \frac{\partial C_{i,f}}{\partial t}

-\nabla\cdot\mathbf{N}{i,f} + a_s k{i,m}(C_{i,f}-C_{i,s}) ]

表面反应则通过单独的催化剂或界面条件来表示。

这种方法可以解析液膜阻力、流体与固体浓度之间的差异以及孔内扩散效应。

拟均相模型

拟均相固定床模型将流体和催化剂视为一个平均化的连续体。非均相反应被转换为有效的体积源项：

[ \frac{\partial C_i}{\partial t}

-\nabla\cdot\mathbf{N}i + a_s r{Si} ]

或者，基于催化剂质量基准，

[ \frac{\partial C_i}{\partial t}

-\nabla\cdot\mathbf{N}i + \rho{b,\mathrm{cat}}r_{i,\mathrm{cat}} ]

这种模型计算更简单，当局部流体-固体梯度不是设计问题的核心时通常适用。

选择表示方式

合适的形式取决于所需的详细程度。反应器尺寸计算可能使用平均体积源项，而催化剂设计或热管理可能需要显式的表面和颗粒内输运方程。

核心的建模问题是：催化剂表面和特定相位的梯度是被解析出来，还是被平均成一个反应器尺度的源项。

理解权衡

避免混用表面单位与体积单位

一个常见错误是将 (\mathrm{mol,m^{-2},s^{-1}}) 直接插入到要求源项单位为 (\mathrm{mol,m^{-3},s^{-1}}) 的方程中。必须先乘以单位反应器体积的界面面积，或将其作为边界通量施加。

使用基于每千克催化剂的速率时也会出现同样的问题。在将其代入体积物种衡算之前，必须先乘以单位反应器体积的催化剂质量。

不要混淆非均相化学与非均相建模

化学上的非均相反应可以通过拟均相模型来表示。在这种情况下，化学反应仍然发生在催化剂表面，但其净效应被平均为一个等效的体积反应速率。

相反，多相反应器模型并不自动需要表面反应项。即使化学反应发生在主体相中，也可能需要特定相位的方程来处理相间输运。

考虑催化剂内部表面

催化剂反应可能发生在内部孔表面上，而不仅仅是在外部颗粒边界上。如果这些表面在几何上没有解析，它们的贡献通常通过有效表面积、基于催化剂质量的速率、有效性因子或体积平均源项来表示。

将内部催化活性当作外部边界条件处理而不考虑可用的内部面积，可能会严重误判反应速率。

保持化学计量一致性

无论速率是体积基准、表面基准还是催化剂质量基准，所有物种源项都必须使用相同的反应进度速率和化学计量系数：

[ R_i=\sum_j \nu_{ij}r_j ]

改变归一化基准会改变单位和比例，但不会改变基本的化学计量关系。

如何应用到你的项目中

正确的实施取决于反应的物理位置以及可用动力学数据的基准。

如果你的主要关注点是均相反应器设计： 直接在物种守恒方程中使用 (R_{Vi})，单位为 (\mathrm{mol,m^{-3},s^{-1}})。
如果你的主要关注点是解析表面动力学： 将 (R_{Si}) 作为界面通量或边界条件施加，单位为 (\mathrm{mol,m^{-2},s^{-1}})。
如果你的主要关注点是固定床反应器放大： 使用 (\rho_{b,\mathrm{cat}}) 或 (a_s) 将基于催化剂质量或表面积的速率转换为有效体积源项。
如果你的主要关注点是催化剂或输运限制： 使用独立的流体和固体方程来解析液膜传递、温度差异和颗粒内扩散。
如果你的主要关注点是快速中试装置模拟： 使用具有明确定义的有效体积反应速率的拟均相模型。

只有当反应速率的物理位置、归一化基准和单位与其所应用的守恒方程或边界条件相匹配时，它在中试规模物种衡算中才是数学上正确的。

总结表：

反应或模型类型	数学表示	典型单位	中试规模用途
均相反应	体积源项，`Rᵥᵢ = Σνᵢⱼrᵥⱼ`	mol·m⁻³·s⁻¹	主体液相或气相反应
解析的非均相反应	表面通量或界面通量跳跃，`Rₛᵢ`	mol·m⁻²·s⁻¹	表面和相界面建模
基于催化剂质量基准	`Rᵥᵢ,eff = ρᵦ,cat rᵢ,cat`	mol·m⁻³ 反应器·s⁻¹	固定床反应器模型
基于催化剂表面积基准	`Rᵥᵢ,eff = aₛ rᵢ,surf`	mol·m⁻³ 反应器·s⁻¹	多孔催化剂和平均模型
拟均相模型	衡算中的有效体积源项	mol·m⁻³·s⁻¹	更快的反应器模拟和放大
显式多相模型	具有相间传递和表面条件的独立相衡算	取决于基准	液膜、孔道和热限制分析