将组分流量作为未知数。 在化学工程单元操作和中试工厂物料平衡中，应将每个未知数设定为单独的组分质量流量、摩尔流量或组分量。避免使用“未知的总物流流量乘以未知的质量分数或摩尔分数”作为未知数，因为这种乘积会产生非线性方程，从而无法使用标准的线性矩阵方法直接求解。

最可靠的线性化方法是先求解组分流量，然后根据结果计算总流量和组成。这会产生形式为 Ax = b 的方程组，当方程独立且自由度为零时，可以使用标准的数值线性代数方法高效求解。

为什么未知数的选择决定了线性度

非线性物流表示法

组分流量通常写作：

[ \dot{m}_i = \dot{m}w_i ]

其中 (\dot{m}) 是总质量流量，(w_i) 是组分质量分数。

如果 (\dot{m}) 和 (w_i) 均为未知数，它们的乘积包含两个未知变量。由此产生的平衡方程是非线性的，即使底层的守恒定律是线性的。

线性组分流量表示法

直接将组分流量本身定义为未知数：

[ \dot{m}_i = x_i ]

此时组分平衡方程可直接写为：

[ \sum_j a_{ij}x_j = b_i ]

其中系数 (a_{ij}) 代表已知的物流关系、分配系数或化学计量系数。未知向量包含的是组分流量，而不是总流量与组成的乘积。

事后恢复物流属性

一旦计算出组分流量，通过以下公式确定总物流流量：

[ \dot{m}_{\text{total}} = \sum_i \dot{m}_i ]

然后使用以下公式计算组分质量分数：

[ w_i = \frac{\dot{m}i}{\dot{m}{\text{total}}} ]

同样的方法也适用于摩尔基准：

[ \dot{n}_{\text{total}} = \sum_i \dot{n}_i \qquad y_i = \frac{\dot{n}i}{\dot{n}{\text{total}}} ]

这实现了线性平衡计算与组成计算的分离。

如何构建平衡模型

建立计算基准

在定义未知数之前选择一个基准。对于连续稳态操作，通常是单位时间的流量；对于间歇操作，可以是一个批次或指定的投料量。

除非明确包含基准转换，否则在整个模型中应始终使用质量流量或摩尔流量。

绘制流程图并对物流编号

工艺流程图应标明每一股进入、流出、循环、分流和中间物流。与未知物流相关的每个组分流量都成为候选变量。

例如，一股双组分出口物流可以表示为：

[ x_1 = \dot{m}{A,\text{out}} \qquad x_2 = \dot{m}{B,\text{out}} ]

而不是由总出口流量和未知组成分数组成的对。

编写独立的组分平衡方程

对于非反应稳态单元：

[ \text{组分输入} = \text{组分输出} ]

对于反应单元，应包含已知的化学计量生产或消耗项。如果反应程度、速率表达式、平衡关系或其他量也是未知数且与其他未知数相乘，模型可能就不再保持线性。

组装矩阵形式

将每个方程整理，使未知项在左侧，已知项在右侧，然后写成：

[ A\mathbf{x}=\mathbf{b} ]

例如：

[ \begin{bmatrix} a_{11} & a_{12}\ a_{21} & a_{22} \end{bmatrix} \begin{bmatrix} x_1\ x_2 \end{bmatrix}

\begin{bmatrix} b_1\ b_2 \end{bmatrix} ]

此处，(A) 包含已知系数，(\mathbf{x}) 包含组分流量未知数，(\mathbf{b}) 包含已知的进料速率、规格或源项。

检查系统是否可解

执行自由度分析

计算：

[ \text{自由度}

\text{未知数数量}

\text{独立方程数量} ]

自由度为零的系统是精确指定的，如果系数矩阵满秩，则可能存在唯一解。

正值表示需要额外的测量值或规格。负值表示模型可能包含冗余方程或不一致的规格。

验证方程独立性

仅计算方程数量是不够的，因为某些平衡方程可能是线性相关的。在求解前检查 (A) 的秩。

对于方阵系统，使用直接求解器（例如下方的代码）获得唯一结果需要矩阵满秩：

import numpy as np

x = np.linalg.solve(A, b)

如果矩阵奇异或秩不足，请检查模型是否存在重复的平衡方程、遗漏的变量或未充分指定的物流。

求解后检查守恒性

将计算出的组分流量代回每个平衡方程。残差向量为：

[ \mathbf{r}=A\mathbf{x}-\mathbf{b} ]

较小的残差表明数值解满足所制定的方程。它们本身并不能证明工艺模型或测量值是正确的。

理解权衡

组分未知数增加了变量数量

使用组分流量可能会比使用“总流量加组成”的公式产生更多的未知数。这通常是一个可以接受的权衡，因为由此产生的方程保持线性，更容易求解和审计。

额外的变量也使组分守恒变得明确，这在中试工厂分析和实验室教学中非常有价值。

可能仍需要组成测量

当过程未充分指定时，组分流量未知数并不能消除对测量组成的需求。它们提供了更好的数学表示，但模型仍然需要足够的独立数据和操作规格。

已知的总流量或已知组成可以在适当的情况下作为线性方程包含在内。

某些工艺模型本质上是非线性的

物料守恒在组分流量上是线性的，但完整的单元操作模型可能并非如此。汽液平衡、压力-流量关系、反应动力学、传热关联式和物性计算都可能引入非线性项。

在可用信息支持的情况下使用线性平衡模型。不要仅仅为了使用矩阵求解器而强行将本质上非线性的工艺行为纳入线性公式。

计算出的负流量需要调查

负的组分流量通常预示着符号错误、遗漏了物流或工艺项、测量值不一致或工艺规格无效。它也可能表明假设的操作状态或模型结构与提供的数据不兼容。

在接受结果之前，请检查物流方向、平衡符号、单位、方程秩和测量输入。

为您的目标做出正确的选择

使用与模型目的相匹配的公式：

如果您的主要重点是稳健的线性求解器： 将各个组分的质量或摩尔流量定义为未知向量，并求解所得的 (A\mathbf{x}=\mathbf{b}) 系统。
如果您的主要重点是确定物流组成： 先求解组分流量，然后通过求和与归一化获得总流量和质量或摩尔分数。
如果您的主要重点是中试工厂数据协调： 将测量的物流和组成信息作为独立方程包含在内，然后检查矩阵秩和平衡残差。
如果您的主要重点是详细的工艺模拟： 在守恒平衡中保留组分流量变量，但在平衡、动力学或物性关系需要时使用非线性求解器。

通过选择组分流量作为主要未知数，您可以使守恒方程透明化，尽可能保持线性，并为数值求解器提供最清晰的求解问题。

汇总表：

建模要素	推荐方法	益处
主要未知数	各个组分的质量或摩尔流量	保持方程线性
矩阵模型	将平衡表示为 `A x = b`	支持标准线性求解器
总流量	求解后对组分流量求和	将计算与组成计算分离
组成	用每个组分流量除以总流量	提供质量或摩尔分数
可解性	检查自由度和矩阵秩	确认是否存在唯一解
验证	评估平衡残差并调查负流量	检测建模或数据问题
非线性行为	对平衡、动力学和物性关系使用非线性求解器	匹配完整的工艺模型

需要为您的教学实验室、中试工厂或工艺研究项目提供可靠的物料平衡模型吗？联系 LABPARK，讨论用于化学工程、生物工艺与生物技术、以及环境与水处理应用的教育和职业单元操作中试工厂。LABPARK 通过实用的中试工厂系统、基于组分的建模支持以及适用于实验室教学、工艺开发和放大研究的设备，帮助大学、研究机构和企业构建更清晰的培训和测试工作流程。