化工分离过程仿真系列教程（5）：MESH 逐板方程组结构——每板四类方程消四类未知量
化工分离过程仿真系列教程（5）：MESH 逐板方程组结构——每板四类方程消四类未知量
版本声明块
- 工具基线：aspenONE V15；RadFrac（Rigorous Multistage Separation——正文首现给全称，保留官方全大写）。方程族口径以权威教学教程为证据（旁证——标注口径），逐字书写以官方文档为准（铁律 9 精神）。
- 所有数值标注为示例值（非官方工况数据）；判定线阈值标注为工程约定（必须标注）。
- 本文目标：把铁律 2 落成逐板账本——简捷法假设在这里显形为"被省略的方程"
一句话结论：每块理论板有四类未知量（该板液相组成、汽相组成、温度、两相流量）与四类方程——MESH：物料衡算（Material）、相平衡（Equilibrium，含 K 值）、加和式（Summation，Σx=1/Σy=1）、热量衡算（Heat）；每板四类方程消四类未知量即逐板账本闭合（铁律 2）。简捷法的"恒摩尔流"与"恒 α"在逐板账本里显形为被省略的热量方程与 K 值变化——假设不是消失，是被替换成常数。
〇、本篇要解决的认知问题
- MESH 四类方程各自消哪个未知量——逐板消元映射的完整格式是什么？
- 简捷法两条假设（CMO、恒 α）在逐板账本里显形的位置——"被省略的方程"与"被替换成常数的方程"分别指哪类？
- 规格层（回流比、板数、进料规格）如何补边界而不重复写物理——与系列账本纪律（铁律 3/4）的对应？
- 逐板欠定与超定的可观察信号是什么——板级局部缺口的现场形态？
- 判定线怎么立——逐板剖面哪些量可与手工/简捷对照（工程容差 <5%——必须标注为约定）？
一、机制解析
1.1 逐板账本与消元映射（价值锚点：MESH 不是四个公式，是每板一套消元）
教学教程口径的单板账本（旁证——逐字以官方文档为准；示例值口径——必须标注）：
| # | 未知量（每板 j） | 提供方程 | 消元映射 |
|---|---|---|---|
| 1 | x_j（液相组成，n 组分→n−1 独立量+加和式） | M：物料衡算（进入=离开） | M_j → x_j |
| 2 | y_j（汽相组成） | E：相平衡 y_ij = K_ij·x_ij | E_j → y_j |
| 3 | T_j（板温——K 值的载体） | H：热量衡算（进出焓+潜热） | H_j → T_j |
| 4 | L_j、V_j（两相流量，n−1 独立量） | S：加和式 Σx=1、Σy=1（与组分物料联立消流量） | S_j → L_j/V_j |
账本闭合口径：每板"四类方程消四类未知量"——组分维按 n 列展开后仍是逐条消元（第 04 篇 gPROMS 数组口径同构）。判定线（练习用）：对任一板能写出 M/E/H/S 四行与消元映射；缺任一行即该板局部欠定（1.4）。
1.2 简捷假设在逐板账本里显形（"被省略/被替换"的准确位置）
- CMO 显形：逐板热量衡算 H_j 被替换为"V、L 沿塔常数"——即 H 类方程退出消元、流量升为全局常数。失效现场（第 03 篇操作线折线）在此账本上可指认：H 类方程回归后流量逐板变化。
- 恒 α 显形：K_ij 被替换为常数（或平均）——E 类方程的 K 值不再随组成/温度响应。失效现场（塔顶/塔底 K 比值偏差 >10%——工程约定）在此账本上可指认：逐板 K 剖面非常数。
经验法则：假设不是消失，是被替换成常数——交付文档写"CMO 成立"必须带 H 类方程的适用性证据（潜热差异量级——工程容差 <5% 口径）；写"恒 α 成立"必须带逐板 K 剖面的常数性证据（<10%——工程约定）。
1.3 规格层补边界，不重复写物理（铁律 3/4 的分离段落点）
- 规格层给的是边界条件：回流比 R、理论板数 N、进料位置与热状态、产品规格——它们进入账本作为"某些量升为已知"；
- 分界线判定：对任一规格能指出"它把哪个量从解改为已知、释放了哪类方程的消元对象"——指不出即重复写物理风险（例如已给全塔板数又逐板规定组成——超定）。
1.4 逐板欠定/超定的可观察信号（板级局部形态）
| 形态 | 现场 | 可观察信号 |
|---|---|---|
| 板级欠定 | 某板缺 H 或 E 类方程（如效率模型只给 y 不给 T 闭合路径——第 08 篇接口） | 求解报该板约束不足；或该板 T/x 给出不可复现任意值（自由变量现场——与系列铁律 2 同构） |
| 板级超定 | 把某板组成当规格同时 M/E 决定它 | 残差无法低于容差（第 04 篇形态 B 的板级版） |
| 假闭合 | 组分物料与加和式写成等价对（线性相关） | 能进求解但残差不降——指出等价对即达标（铁律 2 防它） |
判定线：三形态各制造一次并在日志指认；"板级局部缺口"必须能还原到具体板号与方程类别。
1.5 判定线（逐板剖面的可观察量——工程容差 <5%，必须标注为约定）
① 逐板 T/x/y 剖面与手工递推（近理想体系小塔）对照 <5%（工程容差）； ② 简捷法结果（工作点在经验区间内时）与严格解的关键量（分割、Rmin 对应操作线斜率）对照 <5%——第 18 篇一致化流程在此立线； ③ K 值剖面常数性检查（恒 α 适用域证据——<10%——工程约定）。二、完整代码与逐行剖析
2.1 单板账本（概念级——方程族口径以教学教程为旁证；逐字以官方文档为准；示例值口径——必须标注）
Stage j (以官方文档为准的书写形态——此处给账本口径): M_j: L_{j-1} x_{i,j-1} + V_{j+1} y_{i,j+1} = L_j x_{i,j} + V_j y_{i,j} ← 消 x_j（组分维逐条） E_j: y_{i,j} = K_{i}(T_j, x_j) * x_{i,j} ← 消 y_j（K 随 T/x 响应——恒 α 时退化为常数——1.2） S_j: Σ_i x_{i,j} = 1 ; Σ_i y_{i,j} = 1 ← 与 M 联立消 L_j/V_j H_j: L_{j-1} h_{L,j-1} + V_{j+1} H_{V,j+1} = L_j h_{L,j} + V_j H_{V,j} ← 消 T_j（CMO 时此方程被"流量常数"替换——1.2） 规格层补边界（不重复写物理——1.3）: R、N、进料位置与 q、产品规格 → 把对应量升为已知逐行要点：四行即每板账本的四列消元映射（1.1 表）；"简捷法=删掉 H_j+把 K 当常数"是同一账本的两种退化形态——假设链与账本在此显形（铁律 1/2 的交汇点）。
2.2 一次可复现的板级欠定（判定线）
# 现场：某板只给 Murphree 效率定义（y 与平衡 y* 的关系——第 08 篇口径）但无该板能量闭合路径 # 信号：该板 T 给出不可复现的任意值（自由变量——同一模型两次运行数值不同）；其余板正常——板级局部欠定（1.4） # 处置：补该板热量路径或改全局效率口径并标注假设；交付文档按消元映射逐板列（第 19 篇三件套收口）三、常见报错与排查
- 板级局部欠定—— 现象：个别板量不可复现、整体报约束不足。根因：某类方程缺失（效率/能量路径——1.4）。解法：按 M/E/H/S 逐板列消元映射；缺行即 FAIL 并可还原到板号+类别。
- 把某板组成当规格—— 现象：残差无法低于容差。根因：M/E 已决定该量又被规定（形态 B 板级版——铁律 3）。解法：规格与解二选一；指出释放的消元对象。
- 加和式与组分物料写成等价对—— 现象：假闭合、残差不降。根因：线性相关（1.4 假闭合）。解法：查等价对方程——保留加和式的独立角色（S 类是流量消元的提供方，不是 M 的重复）。
- 假设链与账本脱节—— 现象："CMO 成立"无证据、K 剖面非常数却按恒 α 引用。根因：被替换成常数的方程未给适用性证据（1.2）。解法：交付文档逐条给证据（潜热差异量级、K 剖面常数性——工程约定必须标注）。
四、动手练习
- 单板账本复现：按 2.1 写一个五组分塔的 M/E/H/S 消元映射表。判定标准：每板四类方程消四类未知量；制造一次板级欠定并还原到板号+类别——写出"假设是被替换成常数，不是消失"的一条证据（K 剖面检查 <10%——工程约定）。
- 规格层分界演练：给一个"已给全塔板数又逐板规定组成"的超定场景。判定标准：指出哪条规格把哪个量升为已知、哪类消元被释放；修复后按 1.5 判定线①对照 <5%（工程容差——必须标注为约定）。
- 思考题（无标准答案）：为什么加和式（S）不能被组分物料衡算吸收？验证要点清单：① 写出 S 类方程在流量消元中的独立角色；② 举一个"S 与 M 写成等价对导致假闭合"的场景；③ 给出交付文档逐板消元映射表在二十板塔的列设计（第 19 篇收口呼应）。
五、小结与下一篇预告
本篇把铁律 2 落成逐板账本：M/E/H/S 四类方程消四类未知量、简捷假设显形为被省略/被替换的方程、规格层只补边界不重复写物理。第 06 篇进入 RadFrac 规格账本——塔级"输入还是输出二选一"的完整判定线（铁律 3 在工具层的落点）；第 07 篇把逐板求解的两阶段（初值剖面与主迭代）与第 04 篇 gPROMS 初始化口径同构收拢。
本篇认知问题回显（FAQ）
Q1：MESH 四类方程各自消哪个未知量？
A：每板四类未知量为液相组成 x、汽相组成 y、板温 T、两相流量 L/V；提供方程为 M（物料衡算→x）、E（相平衡 y=Kx→y）、H（热量衡算→T）、S（加和式 Σx=1/Σy=1，与组分物料联立消流量）。每板四类方程消四类未知量即逐板账本闭合（铁律 2）——组分维按 n 列展开后仍是逐条消元。
Q2：简捷法的两条假设在逐板账本里显形的位置是什么？
A：CMO 显形为 H 类方程被"流量沿塔常数"替换；恒 α 显形为 E 类方程的 K 值被替换成常数（不再随组成/温度响应）。假设不是消失，是被替换成常数——交付文档写"成立"必须带适用性证据：潜热差异量级（工程容差 <5%）与逐板 K 剖面常数性（<10%，均必须标注为工程约定）。
Q3：规格层如何补边界而不重复写物理？
A：规格（回流比、板数、进料位置与 q、产品规格）把对应量升为已知——判定线是对任一规格能指出"它把哪个量从解改为已知、释放了哪类方程的消元对象"；指不出即重复写物理风险（已给全塔板数又逐板规定组成即超定）。
Q4：逐板欠定与超定的可观察信号是什么？
A：板级欠定——某板缺 H 或 E 闭合路径时该板量给出不可复现的任意值（自由变量现场）或求解报该板约束不足；板级超定——把某板组成当规格同时 M/E 决定它，残差无法低于容差；假闭合——加和式与组分物料写成等价对（线性相关），能进求解但残差不降。三形态必须可还原到具体板号与方程类别。
Q5：逐板判定的判定线是什么？
A：①逐板 T/x/y 剖面与手工递推对照 <5%（工程容差——必须标注为约定，近理想体系小塔）；②简捷结果（工作点在经验区间内时）与严格解关键量对照 <5%；③K 剖面常数性检查（恒 α 适用域证据——<10%——工程约定）。第三方按逐板消元映射表复跑得到同一剖面才算可追责（铁律 10）。
