gPROMS二次开发教程07物性与热力学——把真物性接进你的模型版本声明块工具/软件gPROMS 桌面建模环境gPROMS ModelBuildergPROMS ModelBuilder为桌面建模工作台检索期官方发布锚点gPROMS Process 2022.1.0适用版本以官方发布说明为准语言/环境gPROMS 语言equation-oriented modelling实体MODELS段PARAMETER/VARIABLE本文目标把物性从模型里的几个常数还原成方程组的一半学会判断物性接入是否可靠、物性参数是否可引用一句话结论在 equation-oriented 模型里物性不是一个装饰性的查询表而是被联立进方程组的一批约束——一个闪蒸模型F*z(i) L*x(i) V*y(i)配y(i) K(i)*x(i)里K(i)随P、T和组成变化物性一旦失真物料衡算会照样守恒结论却整体偏移也因此热力学一致性比单点物性精确更优先。〇、本篇要解决的认知问题Q1物性在模型里的位置到底在哪它为什么不是常数而是方程组的一部分Q2模型与物性之间的抽象层物性包 / 属性调用是如何工作的为什么这层抽象能让你换物性而不改模型Q3相平衡与热力学一致性为什么会直接影响反应与分离模型的结论不一致有哪些可观察的症状Q4物性参数与文献数据要遵循什么来源与引用纪律为什么查到一个数就填进去是高风险动作Q5物性错误会导致哪些典型求解失败到第 11 篇做参数估计时物性参数该怎么处理一、机制解析1.1 物性在模型里的位置它是方程的一部分不是常数初学者的心智模型通常是“物性 一串常数填进PARAMETER段即可”。这在教学例子里勉强成立在真实模型里会立刻失真物性类别在模型中的真实形态是否随工况变化纯组分基础物性分子量、临界参数、生成焓近似常数进PARAMETER基本不变或仅随压力轻微变化密度、黏度、导热系数、热容温度乃至组成、压力的函数显著变化相平衡常数、逸度系数、活度系数P、T、全组成的函数强变化是分离模型的核心非线性来源反应热、溶解热温度的函数影响能量衡算与热点位置【错误心智模型】 【正确心智模型】 模型 物性常数 可解模型 模型方程 ┐ 物性错了 ⇒ 数值偏一点 物性约束 ├→ 联立求解 ⇒ 唯一解 求解器 ┘ 物性错了 ⇒ 非线性行为变了 ⇒ 解的位置与稳定性都可能变关键推论物性失真不会让模型报错只会让它的非线性结构悄悄改变——唯一稳态可能变成多重稳态收敛的迭代可能开始振荡。这类失败在现象上与数值设置不对极像因此排查时必须先排除物性。经验法则拿到莫名不收敛的模型先冻结物性再看把强非线性物性临时换成工作点常数若模型立刻稳定收敛问题就在物性侧。1.2 抽象层模型只认属性调用不认物性细节gPROMS 平台在模型与物性之间隔了一层抽象你的方程调用属性物性系统负责按当前的 T、P、组成算出这个属性的值。这层抽象带来三个直接好处┌──────────────────────────────────────────────┐ │ 你的模型方程 │ │ 守恒方程 / 本构方程 / 相平衡关系 │ └───────────────┬──────────────────────────────┘ │ 属性调用向量T, P, 组成 → 属性值 ┌───────────────▼──────────────────────────────┐ │ 物性 / 热力学 抽象层 │ │ 物性方法选择 · 参数集 · 混合规则 │ └───────────────┬──────────────────────────────┘ ┌───────────────▼──────────────────────────────┐ │ 物性系统的具体实现平台提供 / 外部接入 │ └──────────────────────────────────────────────┘ 收益 1换物性方法 ≠ 改模型。模型里只是调用属性换了底层物性方法 方程一字不改这让同一模型跑多种热力学假设成为低成本实验。 收益 2参数集集中管理。物性参数如二元交互参数集中在一处维护 不必散落在各个模型里避免同一个物性在不同模型里取值不同。 收益 3可追溯。属性调用点就是模型对物性的全部依赖做证据审查时 能一条条列出来核对铁律 2证据分级。铁律 4商业闭源软件的诚实边界属性的具体名称、物性方法的选择方式、参数集的声明与覆盖语法、以及外部物性与高级热力学包的接入方式全部属于 gPROMS 语言与平台细节一律以随附官方《gPROMS ModelBuilder Guide》与平台文档为准。本文只讲抽象层的职责划分与方程如何消费物性结果不点任何未在本系列事实底账中登记的物性产品或方法名。1.3 相平衡与热力学一致性为什么它决定结论相平衡是分离与反应模型的核心非线性来源。以气液平衡为例模型里有三类约束同时存在约束类别表达作用物料衡算F*z(i) L*x(i) V*y(i)组分守恒无论物性对错这一条永远成立相平衡关系y(i) K(i)*x(i)K由物性系统给出决定两相如何分配组分归一化Σx(i) 1、Σy(i) 1保证摩尔分数自洽这三类约束的耦合方式决定了一个重要事实物料衡算永远会守恒因此物性错误不会被守恒方程拦住——它只会把解推向一个守恒但物理上错的位置。这就是物性错误最难查的根本原因。而热力学一致性指的是不同物性属性之间必须相互自洽例如相平衡关系与焓、熵的推导链条不能互相矛盾。一致性被破坏后典型症状有三种能算出负浓度或摩尔分数越界归一化被强行满足但某个x(i)落到 0 以下能量衡算与相平衡互相打架按相平衡算出的组成代入焓模型后能量衡算残差始终不降同一个体系在不同初值下收敛到不同解非线性结构被人为扭曲出现伪多重稳态。最佳实践“物性一致优先于单点物性精确”。一个精确的纯组分常数若与体系其余物性推导链不匹配往往比一套自洽但稍粗的物性更危险。工程做法是整套方法与参数集一起选定、一起记录、一起做回归。1.4 接入路线平台内与平台外物性来源可以分两条路线两条路线的风险点完全不同路线形态主要风险纪律要点平台内物性系统模型调用平台提供的属性参数集在平台内维护选型错误方法与体系不匹配、参数缺失、工况超出适用范围记录方法 参数集 适用范围三件套外部物性接入通过平台支持的外部连接机制引入外部物性/热力学计算连接前提许可模块、外部软件、环境与路径、数据往返与单位约定按铁律 7 列出全部前置条件并声明详细参数以官方文档为准铁律 7外部连接必声明前提凡涉及外部连接第 14 篇将系统展开gO:Python/gO:MATLAB/AdaptorFO/ gPROMS Foreign Process Interface 等机制必须列出前置条件并声明参数细节以官方文档为准。物性接入是外部连接的一个高频场景连接方式的具体支持范围与配置细节以官方文档为准本文不做点名式推荐。1.5 物性参数的来源与引用纪律物性参数是模型里最容易被凭印象填的部分因此本系列把它列为铁律级问题铁律 10数值与参数不臆造动力学常数、物性参数、经验系数一律取自可引用文献或官方示例并标注来源不得凭记忆填数。执行层面要求每个物性参数都能回答三个问题1值是多少2从哪来3在什么范围内有效第三问最常被忽略——很多参数只在有限温区内回归得到把模型跑到温区之外“不收敛其实是外推失效”。1.6 物性引发的典型求解失败症状物性侧根因处置方向迭代不收敛、残差震荡物性强非线性 初值工作点过远冻结物性二分法定位给接近工作点的初值分步求解第 10 篇摩尔分数越界负值/大于 1相平衡关系与归一化冲突、参数不适用检查物性方法适用性检查温度/组成是否超出参数回归范围能量衡算残差不降焓与相平衡来自不一致的推导链统一物性方法避免混用不同来源的焓与相平衡参数求解时间突然暴涨物性计算被高频调用、或非线性陡峭简化分析阶段可先用常物性正式计算再开全物性结果对物性极端敏感体系接近临界区/共沸区缩小工况范围或在结论中显式声明物性假设的边界1.7 与第 11 篇参数估计的分工参数估计时最大的陷阱是把物性参数和动力学参数混在一起估两者在数据上的可辨识性往往很差速率常数调大一点 相平衡常数调小一点可能给出几乎相同的曲线估出来的参数组即使残差很小也可能毫无物理意义。本系列的处置策略是第一步物性参数 → 用独立的物性数据纯组分性质、二元相平衡数据单独确定 不放进动力学回归的决策变量集合。 第二步动力学参数 → 在物性参数固定的前提下用反应/过程数据回归第 11 篇。 第三步一致性检查 → 把两步的结果放回同一模型核对是否仍能复现全部数据。经验法则一次只估一类参数。参数估计的自由度是稀缺资源铁律 8 的另一种表现把它花在物性上动力学就失去了可辨识性。物性侧通常有独立数据可用不该占用过程数据的解释力。二、完整模型与逐行剖析代码 2-1气液平衡闪蒸模型结构示意% --------------------------------------------------------------------------- % 说明本块为结构示意。MODELS / PARAMETER / VARIABLE 为已确证关键字 % 属性调用物性系统 → 模型与相平衡属性的确切写法、物性方法的选择语法 % 以官方《gPROMS ModelBuilder Guide》为准此处仅示意方程如何消费物性。 % 数值为教学示例值正式使用须替换为可引用文献或官方示例数据铁律 10。 % 物理背景单级等温闪蒸。进料 F组成 z在给定 T、P 下分成液相 L 与气相 V。 % --------------------------------------------------------------------------- MODELS PARAMETER N 2 % 组分数示意用二元体系 F 1.0 % 进料摩尔流量示意单位 mol/s z(i) % 进料摩尔分数向量参数Σz 1 T 350.0 % 闪蒸温度示例值 P 1.0e5 % 闪蒸压力示例 Pa VARIABLE x(i) % 液相摩尔分数 y(i) % 气相摩尔分数 K(i) % 相平衡常数由物性系统按 T、P、组成给出物性侧量 L % 液相摩尔流量 V % 气相摩尔流量 %% 物性侧相平衡常数 K(i) %% K(i) 由物性系统按当前 T、P 与组成计算其内部走逸度系数或活度系数路径。 %% 属性调用与物性方法选择的语法以官方《gPROMS ModelBuilder Guide》为准。 %% 本处只声明K(i) 是物性侧的量并把它当作模型的输入属性消费。 %% —— 这一行是科研复现的关键模型对物性的依赖必须显式写出来。 %% 物料衡算每组分一条 F*z(i) L*x(i) V*y(i) % E1 组分 i 的摩尔衡算进料 液相 气相 %% 相平衡关系 y(i) K(i)*x(i) % E2 相平衡气相 相平衡常数 × 液相 %% 归一化 SUM(x(i)) 1.0 % E3 液相摩尔分数归一SUM 为求和算子示意写法 SUM(y(i)) 1.0 % E4 气相摩尔分数归一 %% 总物料衡算 F L V % E5 总摩尔衡算逐行剖析z(i)写在PARAMETER段进料组成是已知量建模前确定。注意它是向量参数其维度必须与x(i)、y(i)一致第 06 篇的维度一致性纪律。K(i)写在VARIABLE段并在注释中明确它是物性侧的量。这不是形式主义读模型的人必须一眼看出这个模型依赖哪些物性属性把依赖藏起来模型就变成不可审查的黑箱。E1是逐组分物料衡算E5是总物料衡算。两者不是独立的把 E1 对全部组分求和即得 E5在归一化与总流量定义下。写全两条可以让物理意义更清楚但自由度核算时必须意识到这一点。E3/E4两条归一化这两条之中只有一条是独立约束下面自由度核对里会证明。这是闪蒸模型的经典结构性陷阱。全式采用残差式第 05 篇的规范F*z(i) - L*x(i) - V*y(i) 0的写法把进料减两股出料原样保留便于逐项核对是否漏项。自由度核对铁律 8——这是本篇最值得带走的一段项数量说明未知数x(i)N 2液相组成未知数y(i)N 2气相组成未知数K(i)N 2相平衡常数由物性系统给出但在方程组内是待定的非线性量未知数L、V2两相流量未知数合计3N 2 8方程 E1物料衡算N 2方程 E2相平衡N 2方程 E3/E4归一化2其中仅 1 条独立见下方证明方程 E5总物料衡算1有效方程合计2 2 1 1 6对 N 2 时E1 两条 E2 两条 归一化一条 总衡算一条归一化冗余证明对 N 2可逐项推广 把 E1 对 i 求和 F*Σz(i) L*Σx(i) V*Σy(i) 代入 Σz(i) 1进料归一、Σx(i) 1E3 F L V*Σy(i) 再代入 E5F L V Σy(i) 1 ← 即 E4 结论E4 可由 E1(全部) E3 E5 导出因此 E3 与 E4 中只有一个带来新信息。这意味着什么若把 E3、E4 都当作独立约束硬写进去模型看起来方程比未知数多一条实际自由度却是正确的——因为求解器面对的是线性相关的方程组不会报错但会给出奇怪的诊断信息如方程冗余或残差不降。初学者常在此处把一条归一化删掉结果模型真正变成欠定。正确处置是保留两条提高可读性与鲁棒性但在模型头里写明两条归一化中仅一条独立。铁律 8建模先算自由度自由度核算必须识别线性相关/可导出的方程。方程数不等于独立约束数闪蒸归一化、总衡算与分衡算之间的依赖关系是最常见的数不对但解得了的场合。代码 2-2物性消费清单——把模型对物性的依赖显式登记结构示意% % 物性消费清单模型对物性的全部依赖逐条登记 % 用途换物性方法时按此清单核对做证据审查时按此清单追溯写论文时按此 % 清单在物性假设一节逐条声明。 % % 序 | 模型中的量 | 消费方式 | 依赖变量 | 参数来源 | 适用范围 % ---------------------------------------------------------------------------------- % 1 | K(i) 相平衡常数 | 属性调用 | T, P, x(i), y(i)| 参数集 A | 273~450 K % 2 | 液相密度 rho_L | 属性调用 | T, P, x(i) | 参数集 A | 常压附近 % 3 | 气相黏度 mu_V | 属性调用 | T, P, y(i) | 参数集 B | 300~600 K % 4 | 反应焓 dH_rxn | 属性调用生成焓差| T | 参数集 C | 标准态附近 % % --------------------------------------------------------------------------- % 换物性方法时的核对动作三步 % 第 1 步清单第 4 列的依赖变量是否与新方法一致如新方法还依赖压力 % 第 2 步清单第 5 列的参数来源是否齐全缺失参数必须先补齐再试算。 % 第 3 步清单第 6 列的适用范围是否覆盖全部工况超范围必须改方法或缩工况。 % % --------------------------------------------------------------------------- % 回归基线铁律 9改完必验 % 物性方法/参数集一经变更必须用固定算例复跑并对比下列量 % 基线量液相组成 x(i)、气相组成 y(i)、两相流量 L 与 V、目标转化率/回收率 % 对比阈值由你的工程精度要求决定但变更物性后结果一致必须留下记录 % 逐行剖析这张清单把物性依赖从散落在方程里的若干调用变成一张可审查的表。它的三个用途分别对应三类真实场景换物性方法时的核对表、证据审查时的追溯链铁律 2官方源优于第三方、论文写作时的物性假设小节素材。最后一段的回归基线要求则是把铁律 9 落到物性这个最容易被悄悄改掉的位置上——物性变更往往被当成无关紧要的调整而它恰恰是能让所有历史结论失效的那类变更。代码 2-3含相平衡的分馏/吸收类模型的物性接入骨架结构示意% --------------------------------------------------------------------------- % 说明本块为结构示意示范反应 相平衡 能量衡算三者耦合时的物性消费 % 骨架。属性调用、物性方法选择与参数集声明的语法以官方《gPROMS % ModelBuilder Guide》为准。数值为教学示例值铁律 10。 % 物理背景单个平衡级理论板进料与两股出料含一股热负荷 Q。 % --------------------------------------------------------------------------- MODELS PARAMETER N 2 % 组分数 z(i) % 进料摩尔分数Σz 1 F 1.0 % 进料摩尔流量 P 1.0e5 % 级压力 Q 0.0 % 级热负荷示意 W0 表示绝热级 VARIABLE x(i) % 液相摩尔分数出液相 y(i) % 气相摩尔分数出气相 K(i) % 相平衡常数物性侧 h_L % 液相摩尔焓物性侧 h_V % 气相摩尔焓物性侧 h_F % 进料摩尔焓物性侧 T % 级温度由相平衡与能量衡算共同决定 L % 液相流量 V % 气相流量 %% 1) 物性侧量全部由物性系统按 (T, P, 组成) 给出 %% K(i)、h_L、h_V、h_F 的调用语法以官方《gPROMS ModelBuilder Guide》为准 %% 它们是本模型对物性的全部依赖已在物性消费清单中登记。 %% 2) 相平衡把物性结果消费进方程 y(i) K(i)*x(i) % E1 相平衡连接气相与液相组成 %% 3) 逐组分物料衡算 F*z(i) L*x(i) V*y(i) % E2 组分 i进料 出液相 出气相 %% 4) 能量衡算焓必须与相平衡来自一致的热力学链条 F*h_F Q L*h_L V*h_V % E3 级能量衡算残差式 %% 5) 归一化 SUM(x(i)) 1.0 % E4 液相归一 SUM(y(i)) 1.0 % E5 气相归一与 E4 中仅一条独立见模型头 %% 6) 总物料衡算 F L V % E6 总摩尔衡算逐行剖析E3是本块的教学重点四个焓值全部来自物性侧。若相平衡方法产生K(i)与焓模型来自不同的推导链或参数集E3与E1、E2就会互相矛盾——残差在数学上可以下降物理上得到的温度与组成却是焓或相平衡其中一方说了假话的结果。这就是 §1.3一致性优先于单点精确的具体后果。T在本模型中不是参数而是变量闪蒸温度由能量衡算与相平衡共同决定。这解释了为什么物性参数必须带适用范围——温度是被解出来的你无法事先保证它落在参数回归的温区内。Q 0.0表示绝热级把Q改为非零值即转成带热负荷的级这正是精馏塔逐板计算的单元级原型。结构实验只需改一行第 06 篇参数化思路的延续。注释里再次出现以官方《gPROMS ModelBuilder Guide》为准——这不是啰嗦而是本系列的双层证据制在每一处涉及闭源语言细节的位置都必须显式落地。三、常见报错与排查报错 3-1相平衡不收敛 / 残差长时间抖动。现象迭代过程残差震荡或长时间不下降尤其出现在高温高压或接近临界区。根因物性强非线性相平衡关系在工作点附近陡峭或初值远离工作点。解法先做冻结物性二分法把所有强非线性物性临时替换为工作点常数判断问题在物性侧还是求解器侧其次给T、x(i)、y(i)设接近预期的初值最后按第 10 篇的初始化思路分步求解官方能力名Initialization procedures in gPROMS。报错 3-2摩尔分数越界或出现负值。现象求解成功但x(i)或y(i)出现负值或大于 1。根因相平衡关系与实际工况不匹配物性方法选错、参数超出回归范围归一化被强行满足而某一相组分被挤成负值。解法核对物性方法的适用体系与温压范围检查参数集是否缺项在物性消费清单里补登适用范围并据此缩小工况。报错 3-3能量衡算残差始终不降。现象物料衡算与相平衡看起来都满足唯独能量衡算残差居高不下。根因焓与相平衡来自不一致的热力学链条例如焓用一套参数、相平衡用另一套。解法把焓与相平衡统一到同一物性方法与同一参数集若必须混用则退回等温/绝热简化版本先跑通再逐步放开。报错 3-4报方程冗余或结构性异常但模型其实可解。现象诊断信息提示方程数多于未知数或提示某些方程线性相关。根因闪蒸/平衡级模型中归一化方程与总物料衡算、逐组分衡算之间存在可导出关系本篇代码 2-1 已给证明。解法在模型头里写明两条归一化中仅一条独立不要盲目删掉归一化——删掉会让模型真正欠定正确做法是保留并把冗余关系记录在案。报错 3-5物性参数缺失或超出适用范围外推失效。现象报错提示某物性参数不可用或模型在某个工况下突然失败。根因该组分/二元对的参数未在参数集中提供或解出的温度超出参数回归温区。解法按物性消费清单逐条核对参数完整性对超范围工况改为适用范围更宽的方法或在结论中显式声明物性假设的边界并禁止用超范围结果做工程判断。报错 3-6求解时间异常增长。现象加了物性之后求解时间成倍增加。根因物性计算被高频调用且属性本身非线性强。解法分析阶段先用常物性做趋势判断正式计算再打开全物性把物性简化假设与正式计算假设分别记录铁律 9两者结果必须对比留痕。四、动手练习练习 1自由度与冗余为代码 2-1 做完整自由度核对要求写出归一化冗余的三行证明对 N 2。判定证明步骤完整能说明为什么不能简单删掉 E4模型求解成功且Σx、Σy均等于 1数值误差范围内。练习 2物性冻结二分法把代码 2-1 中的K(i)临时替换为两个常数取一个合理工作点值重新求解并与全物性版本对比x(i)、y(i)、L、V。判定记录两组结果的差异幅度能据此判断该体系对相平衡的敏感程度并把差异写入工程记录铁律 9。练习 3物性消费清单落地为你的课题里的模型按代码 2-2 格式写一份物性消费清单。判定至少 5 行每行都有依赖变量、参数来源、适用范围三列且非空换物性方法的三步核对动作能逐条执行。练习 4范围越界观察把代码 2-1 的T从示例值逐步提高到清单中的适用范围之外观察求解行为变化。判定能记录到至少一种异常不收敛、越界、时间暴涨能说明外推失效与物性方法选错在现象上的区别。练习 5思考题无标准答案为什么热力学一致性比单点物性精确更优先请用残差可以下降但结论错误这一现象给出两条理由。五、小结与下一篇预告本篇把物性从模型里的常数还原为方程组里的一半物性是被联立进方程的一批约束它的失真不会破坏物料衡算的守恒性只会把解推到守恒但错误的位置模型与物性之间的抽象层属性调用让你能换物性而几乎不改模型前提是把模型对物性的全部依赖显式登记相平衡与焓必须来自一致的热力学链条一致性优先于单点精确物性参数必须回答值多少、从哪来、在什么范围内有效三问铁律 10闪蒸类模型里两条归一化仅一条独立是必须写进模型头的结构性事实铁律 8。最后参数估计时一次只估一类参数——物性参数用独立物性数据先定动力学参数留给过程数据第 11 篇。第 08 篇进入动态世界第一个动态模型——从LUMPED到DISTRIBUTED。我们讲清动态与稳态的本质差别时间自变量、初值、微分方程、集中参数LUMPED与分布参数DISTRIBUTED模型的分野对应官方能力名Implementing a lumped model in gPROMS与Implementing a distributed model in gPROMS、空间离散与网格/阶数选择对精度与耗时的影响以及动态仿真对初值的敏感性——那是从会算稳态到会算过程的分水岭。本篇认知问题回显FAQQ1物性在模型里的位置到底在哪为什么它不是常数A只有纯组分基础物性分子量、临界参数、生成焓近似常数可放进 PARAMETER 段密度、黏度、热容随温度乃至组成与压力变化相平衡常数、逸度系数、活度系数则是 P、T 与全组成的强非线性函数反应焓也随温度变化。因此物性不是查询表而是被联立进方程组的约束。它失真的后果不是报错而是非线性结构悄悄改变例如唯一稳态变多重稳态、收敛变振荡所以排查时要先排除物性。Q2模型与物性之间的抽象层是如何工作的为什么能换物性而不改模型A抽象层让模型只做属性调用即把当前的 T、P、组成作为输入、取回属性值由物性系统负责具体计算。这带来三个好处换物性方法时模型一字不改同一模型跑多种热力学假设成为低成本实验物性参数集集中维护避免同一物性在不同模型里取值不同属性调用点即模型对物性的全部依赖便于逐条追溯审查。属性的具体名称、物性方法选择、参数集覆盖语法与外部物性接入方式均以官方《gPROMS ModelBuilder Guide》为准。Q3相平衡与热力学一致性为什么会直接影响结论不一致有哪些症状A因为物料衡算无论物性对错都永远成立物性错误不会被守恒方程拦住只会把解推向守恒但物理上错的位置。热力学一致性要求相平衡关系与焓、熵等属性相互自洽。不一致的三种典型症状是摩尔分数越界甚至为负归一化被强行满足而某组分被挤成负值能量衡算残差始终不降因为焓与相平衡来自不同推导链同一体系在不同初值下收敛到不同解出现伪多重稳态。Q4物性参数与文献数据要遵循什么纪律A按铁律 10物性参数必须取自可引用文献或官方示例并标注来源不得凭记忆填数。执行上每个参数都要回答三问值是多少、从哪来、在什么范围内有效。第三问最常被忽略因为动态模型里的温度往往是被解出来的变量无法事先保证落在参数回归温区内把模型跑到温区之外时所谓不收敛其实是外推失效。此外所有物性依赖应按清单登记并在物性方法或参数集变更后按铁律 9 用固定算例复跑对比基线。Q5物性错误会造成哪些典型失败参数估计时物性参数怎么处理A典型失败包括相平衡不收敛或残差抖动、摩尔分数越界、能量衡算残差不降、诊断提示方程冗余而模型实际可解、物性参数缺失或超出适用范围、求解时间异常增长。参数估计的处置是三步先用独立的物性数据纯组分性质、二元相平衡数据单独确定物性参数不放进动力学回归的决策变量集合再在物性固定的前提下用过程数据回归动力学参数最后把两步结果放回同一模型核对能否复现全部数据。原则是一次只估一类参数。本文涉及的语言细节以所购版本随附的官方 gPROMS 文档为准。 
长河君 · 2026/10/10 2:44:05 · 阅读 1.2k gPROMS二次开发教程（07）：物性与热力学——把真物性接进你的模型
gPROMS二次开发教程07物性与热力学——把真物性接进你的模型版本声明块工具/软件gPROMS 桌面建模环境gPROMS ModelBuildergPROMS ModelBuilder为桌面建模工作台检索期官方发布锚点gPROMS Process 2022.1.0适用版本以官方发布说明为准语言/环境gPROMS 语言equation-oriented modelling实体MODELS段PARAMETER/VARIABLE本文目标把物性从模型里的几个常数还原成方程组的一半学会判断物性接入是否可靠、物性参数是否可引用一句话结论在 equation-oriented 模型里物性不是一个装饰性的查询表而是被联立进方程组的一批约束——一个闪蒸模型F*z(i) L*x(i) V*y(i)配y(i) K(i)*x(i)里K(i)随P、T和组成变化物性一旦失真物料衡算会照样守恒结论却整体偏移也因此热力学一致性比单点物性精确更优先。〇、本篇要解决的认知问题Q1物性在模型里的位置到底在哪它为什么不是常数而是方程组的一部分Q2模型与物性之间的抽象层物性包 / 属性调用是如何工作的为什么这层抽象能让你换物性而不改模型Q3相平衡与热力学一致性为什么会直接影响反应与分离模型的结论不一致有哪些可观察的症状Q4物性参数与文献数据要遵循什么来源与引用纪律为什么查到一个数就填进去是高风险动作Q5物性错误会导致哪些典型求解失败到第 11 篇做参数估计时物性参数该怎么处理一、机制解析1.1 物性在模型里的位置它是方程的一部分不是常数初学者的心智模型通常是“物性 一串常数填进PARAMETER段即可”。这在教学例子里勉强成立在真实模型里会立刻失真物性类别在模型中的真实形态是否随工况变化纯组分基础物性分子量、临界参数、生成焓近似常数进PARAMETER基本不变或仅随压力轻微变化密度、黏度、导热系数、热容温度乃至组成、压力的函数显著变化相平衡常数、逸度系数、活度系数P、T、全组成的函数强变化是分离模型的核心非线性来源反应热、溶解热温度的函数影响能量衡算与热点位置【错误心智模型】 【正确心智模型】 模型 物性常数 可解模型 模型方程 ┐ 物性错了 ⇒ 数值偏一点 物性约束 ├→ 联立求解 ⇒ 唯一解 求解器 ┘ 物性错了 ⇒ 非线性行为变了 ⇒ 解的位置与稳定性都可能变关键推论物性失真不会让模型报错只会让它的非线性结构悄悄改变——唯一稳态可能变成多重稳态收敛的迭代可能开始振荡。这类失败在现象上与数值设置不对极像因此排查时必须先排除物性。经验法则拿到莫名不收敛的模型先冻结物性再看把强非线性物性临时换成工作点常数若模型立刻稳定收敛问题就在物性侧。1.2 抽象层模型只认属性调用不认物性细节gPROMS 平台在模型与物性之间隔了一层抽象你的方程调用属性物性系统负责按当前的 T、P、组成算出这个属性的值。这层抽象带来三个直接好处┌──────────────────────────────────────────────┐ │ 你的模型方程 │ │ 守恒方程 / 本构方程 / 相平衡关系 │ └───────────────┬──────────────────────────────┘ │ 属性调用向量T, P, 组成 → 属性值 ┌───────────────▼──────────────────────────────┐ │ 物性 / 热力学 抽象层 │ │ 物性方法选择 · 参数集 · 混合规则 │ └───────────────┬──────────────────────────────┘ ┌───────────────▼──────────────────────────────┐ │ 物性系统的具体实现平台提供 / 外部接入 │ └──────────────────────────────────────────────┘ 收益 1换物性方法 ≠ 改模型。模型里只是调用属性换了底层物性方法 方程一字不改这让同一模型跑多种热力学假设成为低成本实验。 收益 2参数集集中管理。物性参数如二元交互参数集中在一处维护 不必散落在各个模型里避免同一个物性在不同模型里取值不同。 收益 3可追溯。属性调用点就是模型对物性的全部依赖做证据审查时 能一条条列出来核对铁律 2证据分级。铁律 4商业闭源软件的诚实边界属性的具体名称、物性方法的选择方式、参数集的声明与覆盖语法、以及外部物性与高级热力学包的接入方式全部属于 gPROMS 语言与平台细节一律以随附官方《gPROMS ModelBuilder Guide》与平台文档为准。本文只讲抽象层的职责划分与方程如何消费物性结果不点任何未在本系列事实底账中登记的物性产品或方法名。1.3 相平衡与热力学一致性为什么它决定结论相平衡是分离与反应模型的核心非线性来源。以气液平衡为例模型里有三类约束同时存在约束类别表达作用物料衡算F*z(i) L*x(i) V*y(i)组分守恒无论物性对错这一条永远成立相平衡关系y(i) K(i)*x(i)K由物性系统给出决定两相如何分配组分归一化Σx(i) 1、Σy(i) 1保证摩尔分数自洽这三类约束的耦合方式决定了一个重要事实物料衡算永远会守恒因此物性错误不会被守恒方程拦住——它只会把解推向一个守恒但物理上错的位置。这就是物性错误最难查的根本原因。而热力学一致性指的是不同物性属性之间必须相互自洽例如相平衡关系与焓、熵的推导链条不能互相矛盾。一致性被破坏后典型症状有三种能算出负浓度或摩尔分数越界归一化被强行满足但某个x(i)落到 0 以下能量衡算与相平衡互相打架按相平衡算出的组成代入焓模型后能量衡算残差始终不降同一个体系在不同初值下收敛到不同解非线性结构被人为扭曲出现伪多重稳态。最佳实践“物性一致优先于单点物性精确”。一个精确的纯组分常数若与体系其余物性推导链不匹配往往比一套自洽但稍粗的物性更危险。工程做法是整套方法与参数集一起选定、一起记录、一起做回归。1.4 接入路线平台内与平台外物性来源可以分两条路线两条路线的风险点完全不同路线形态主要风险纪律要点平台内物性系统模型调用平台提供的属性参数集在平台内维护选型错误方法与体系不匹配、参数缺失、工况超出适用范围记录方法 参数集 适用范围三件套外部物性接入通过平台支持的外部连接机制引入外部物性/热力学计算连接前提许可模块、外部软件、环境与路径、数据往返与单位约定按铁律 7 列出全部前置条件并声明详细参数以官方文档为准铁律 7外部连接必声明前提凡涉及外部连接第 14 篇将系统展开gO:Python/gO:MATLAB/AdaptorFO/ gPROMS Foreign Process Interface 等机制必须列出前置条件并声明参数细节以官方文档为准。物性接入是外部连接的一个高频场景连接方式的具体支持范围与配置细节以官方文档为准本文不做点名式推荐。1.5 物性参数的来源与引用纪律物性参数是模型里最容易被凭印象填的部分因此本系列把它列为铁律级问题铁律 10数值与参数不臆造动力学常数、物性参数、经验系数一律取自可引用文献或官方示例并标注来源不得凭记忆填数。执行层面要求每个物性参数都能回答三个问题1值是多少2从哪来3在什么范围内有效第三问最常被忽略——很多参数只在有限温区内回归得到把模型跑到温区之外“不收敛其实是外推失效”。1.6 物性引发的典型求解失败症状物性侧根因处置方向迭代不收敛、残差震荡物性强非线性 初值工作点过远冻结物性二分法定位给接近工作点的初值分步求解第 10 篇摩尔分数越界负值/大于 1相平衡关系与归一化冲突、参数不适用检查物性方法适用性检查温度/组成是否超出参数回归范围能量衡算残差不降焓与相平衡来自不一致的推导链统一物性方法避免混用不同来源的焓与相平衡参数求解时间突然暴涨物性计算被高频调用、或非线性陡峭简化分析阶段可先用常物性正式计算再开全物性结果对物性极端敏感体系接近临界区/共沸区缩小工况范围或在结论中显式声明物性假设的边界1.7 与第 11 篇参数估计的分工参数估计时最大的陷阱是把物性参数和动力学参数混在一起估两者在数据上的可辨识性往往很差速率常数调大一点 相平衡常数调小一点可能给出几乎相同的曲线估出来的参数组即使残差很小也可能毫无物理意义。本系列的处置策略是第一步物性参数 → 用独立的物性数据纯组分性质、二元相平衡数据单独确定 不放进动力学回归的决策变量集合。 第二步动力学参数 → 在物性参数固定的前提下用反应/过程数据回归第 11 篇。 第三步一致性检查 → 把两步的结果放回同一模型核对是否仍能复现全部数据。经验法则一次只估一类参数。参数估计的自由度是稀缺资源铁律 8 的另一种表现把它花在物性上动力学就失去了可辨识性。物性侧通常有独立数据可用不该占用过程数据的解释力。二、完整模型与逐行剖析代码 2-1气液平衡闪蒸模型结构示意% --------------------------------------------------------------------------- % 说明本块为结构示意。MODELS / PARAMETER / VARIABLE 为已确证关键字 % 属性调用物性系统 → 模型与相平衡属性的确切写法、物性方法的选择语法 % 以官方《gPROMS ModelBuilder Guide》为准此处仅示意方程如何消费物性。 % 数值为教学示例值正式使用须替换为可引用文献或官方示例数据铁律 10。 % 物理背景单级等温闪蒸。进料 F组成 z在给定 T、P 下分成液相 L 与气相 V。 % --------------------------------------------------------------------------- MODELS PARAMETER N 2 % 组分数示意用二元体系 F 1.0 % 进料摩尔流量示意单位 mol/s z(i) % 进料摩尔分数向量参数Σz 1 T 350.0 % 闪蒸温度示例值 P 1.0e5 % 闪蒸压力示例 Pa VARIABLE x(i) % 液相摩尔分数 y(i) % 气相摩尔分数 K(i) % 相平衡常数由物性系统按 T、P、组成给出物性侧量 L % 液相摩尔流量 V % 气相摩尔流量 %% 物性侧相平衡常数 K(i) %% K(i) 由物性系统按当前 T、P 与组成计算其内部走逸度系数或活度系数路径。 %% 属性调用与物性方法选择的语法以官方《gPROMS ModelBuilder Guide》为准。 %% 本处只声明K(i) 是物性侧的量并把它当作模型的输入属性消费。 %% —— 这一行是科研复现的关键模型对物性的依赖必须显式写出来。 %% 物料衡算每组分一条 F*z(i) L*x(i) V*y(i) % E1 组分 i 的摩尔衡算进料 液相 气相 %% 相平衡关系 y(i) K(i)*x(i) % E2 相平衡气相 相平衡常数 × 液相 %% 归一化 SUM(x(i)) 1.0 % E3 液相摩尔分数归一SUM 为求和算子示意写法 SUM(y(i)) 1.0 % E4 气相摩尔分数归一 %% 总物料衡算 F L V % E5 总摩尔衡算逐行剖析z(i)写在PARAMETER段进料组成是已知量建模前确定。注意它是向量参数其维度必须与x(i)、y(i)一致第 06 篇的维度一致性纪律。K(i)写在VARIABLE段并在注释中明确它是物性侧的量。这不是形式主义读模型的人必须一眼看出这个模型依赖哪些物性属性把依赖藏起来模型就变成不可审查的黑箱。E1是逐组分物料衡算E5是总物料衡算。两者不是独立的把 E1 对全部组分求和即得 E5在归一化与总流量定义下。写全两条可以让物理意义更清楚但自由度核算时必须意识到这一点。E3/E4两条归一化这两条之中只有一条是独立约束下面自由度核对里会证明。这是闪蒸模型的经典结构性陷阱。全式采用残差式第 05 篇的规范F*z(i) - L*x(i) - V*y(i) 0的写法把进料减两股出料原样保留便于逐项核对是否漏项。自由度核对铁律 8——这是本篇最值得带走的一段项数量说明未知数x(i)N 2液相组成未知数y(i)N 2气相组成未知数K(i)N 2相平衡常数由物性系统给出但在方程组内是待定的非线性量未知数L、V2两相流量未知数合计3N 2 8方程 E1物料衡算N 2方程 E2相平衡N 2方程 E3/E4归一化2其中仅 1 条独立见下方证明方程 E5总物料衡算1有效方程合计2 2 1 1 6对 N 2 时E1 两条 E2 两条 归一化一条 总衡算一条归一化冗余证明对 N 2可逐项推广 把 E1 对 i 求和 F*Σz(i) L*Σx(i) V*Σy(i) 代入 Σz(i) 1进料归一、Σx(i) 1E3 F L V*Σy(i) 再代入 E5F L V Σy(i) 1 ← 即 E4 结论E4 可由 E1(全部) E3 E5 导出因此 E3 与 E4 中只有一个带来新信息。这意味着什么若把 E3、E4 都当作独立约束硬写进去模型看起来方程比未知数多一条实际自由度却是正确的——因为求解器面对的是线性相关的方程组不会报错但会给出奇怪的诊断信息如方程冗余或残差不降。初学者常在此处把一条归一化删掉结果模型真正变成欠定。正确处置是保留两条提高可读性与鲁棒性但在模型头里写明两条归一化中仅一条独立。铁律 8建模先算自由度自由度核算必须识别线性相关/可导出的方程。方程数不等于独立约束数闪蒸归一化、总衡算与分衡算之间的依赖关系是最常见的数不对但解得了的场合。代码 2-2物性消费清单——把模型对物性的依赖显式登记结构示意% % 物性消费清单模型对物性的全部依赖逐条登记 % 用途换物性方法时按此清单核对做证据审查时按此清单追溯写论文时按此 % 清单在物性假设一节逐条声明。 % % 序 | 模型中的量 | 消费方式 | 依赖变量 | 参数来源 | 适用范围 % ---------------------------------------------------------------------------------- % 1 | K(i) 相平衡常数 | 属性调用 | T, P, x(i), y(i)| 参数集 A | 273~450 K % 2 | 液相密度 rho_L | 属性调用 | T, P, x(i) | 参数集 A | 常压附近 % 3 | 气相黏度 mu_V | 属性调用 | T, P, y(i) | 参数集 B | 300~600 K % 4 | 反应焓 dH_rxn | 属性调用生成焓差| T | 参数集 C | 标准态附近 % % --------------------------------------------------------------------------- % 换物性方法时的核对动作三步 % 第 1 步清单第 4 列的依赖变量是否与新方法一致如新方法还依赖压力 % 第 2 步清单第 5 列的参数来源是否齐全缺失参数必须先补齐再试算。 % 第 3 步清单第 6 列的适用范围是否覆盖全部工况超范围必须改方法或缩工况。 % % --------------------------------------------------------------------------- % 回归基线铁律 9改完必验 % 物性方法/参数集一经变更必须用固定算例复跑并对比下列量 % 基线量液相组成 x(i)、气相组成 y(i)、两相流量 L 与 V、目标转化率/回收率 % 对比阈值由你的工程精度要求决定但变更物性后结果一致必须留下记录 % 逐行剖析这张清单把物性依赖从散落在方程里的若干调用变成一张可审查的表。它的三个用途分别对应三类真实场景换物性方法时的核对表、证据审查时的追溯链铁律 2官方源优于第三方、论文写作时的物性假设小节素材。最后一段的回归基线要求则是把铁律 9 落到物性这个最容易被悄悄改掉的位置上——物性变更往往被当成无关紧要的调整而它恰恰是能让所有历史结论失效的那类变更。代码 2-3含相平衡的分馏/吸收类模型的物性接入骨架结构示意% --------------------------------------------------------------------------- % 说明本块为结构示意示范反应 相平衡 能量衡算三者耦合时的物性消费 % 骨架。属性调用、物性方法选择与参数集声明的语法以官方《gPROMS % ModelBuilder Guide》为准。数值为教学示例值铁律 10。 % 物理背景单个平衡级理论板进料与两股出料含一股热负荷 Q。 % --------------------------------------------------------------------------- MODELS PARAMETER N 2 % 组分数 z(i) % 进料摩尔分数Σz 1 F 1.0 % 进料摩尔流量 P 1.0e5 % 级压力 Q 0.0 % 级热负荷示意 W0 表示绝热级 VARIABLE x(i) % 液相摩尔分数出液相 y(i) % 气相摩尔分数出气相 K(i) % 相平衡常数物性侧 h_L % 液相摩尔焓物性侧 h_V % 气相摩尔焓物性侧 h_F % 进料摩尔焓物性侧 T % 级温度由相平衡与能量衡算共同决定 L % 液相流量 V % 气相流量 %% 1) 物性侧量全部由物性系统按 (T, P, 组成) 给出 %% K(i)、h_L、h_V、h_F 的调用语法以官方《gPROMS ModelBuilder Guide》为准 %% 它们是本模型对物性的全部依赖已在物性消费清单中登记。 %% 2) 相平衡把物性结果消费进方程 y(i) K(i)*x(i) % E1 相平衡连接气相与液相组成 %% 3) 逐组分物料衡算 F*z(i) L*x(i) V*y(i) % E2 组分 i进料 出液相 出气相 %% 4) 能量衡算焓必须与相平衡来自一致的热力学链条 F*h_F Q L*h_L V*h_V % E3 级能量衡算残差式 %% 5) 归一化 SUM(x(i)) 1.0 % E4 液相归一 SUM(y(i)) 1.0 % E5 气相归一与 E4 中仅一条独立见模型头 %% 6) 总物料衡算 F L V % E6 总摩尔衡算逐行剖析E3是本块的教学重点四个焓值全部来自物性侧。若相平衡方法产生K(i)与焓模型来自不同的推导链或参数集E3与E1、E2就会互相矛盾——残差在数学上可以下降物理上得到的温度与组成却是焓或相平衡其中一方说了假话的结果。这就是 §1.3一致性优先于单点精确的具体后果。T在本模型中不是参数而是变量闪蒸温度由能量衡算与相平衡共同决定。这解释了为什么物性参数必须带适用范围——温度是被解出来的你无法事先保证它落在参数回归的温区内。Q 0.0表示绝热级把Q改为非零值即转成带热负荷的级这正是精馏塔逐板计算的单元级原型。结构实验只需改一行第 06 篇参数化思路的延续。注释里再次出现以官方《gPROMS ModelBuilder Guide》为准——这不是啰嗦而是本系列的双层证据制在每一处涉及闭源语言细节的位置都必须显式落地。三、常见报错与排查报错 3-1相平衡不收敛 / 残差长时间抖动。现象迭代过程残差震荡或长时间不下降尤其出现在高温高压或接近临界区。根因物性强非线性相平衡关系在工作点附近陡峭或初值远离工作点。解法先做冻结物性二分法把所有强非线性物性临时替换为工作点常数判断问题在物性侧还是求解器侧其次给T、x(i)、y(i)设接近预期的初值最后按第 10 篇的初始化思路分步求解官方能力名Initialization procedures in gPROMS。报错 3-2摩尔分数越界或出现负值。现象求解成功但x(i)或y(i)出现负值或大于 1。根因相平衡关系与实际工况不匹配物性方法选错、参数超出回归范围归一化被强行满足而某一相组分被挤成负值。解法核对物性方法的适用体系与温压范围检查参数集是否缺项在物性消费清单里补登适用范围并据此缩小工况。报错 3-3能量衡算残差始终不降。现象物料衡算与相平衡看起来都满足唯独能量衡算残差居高不下。根因焓与相平衡来自不一致的热力学链条例如焓用一套参数、相平衡用另一套。解法把焓与相平衡统一到同一物性方法与同一参数集若必须混用则退回等温/绝热简化版本先跑通再逐步放开。报错 3-4报方程冗余或结构性异常但模型其实可解。现象诊断信息提示方程数多于未知数或提示某些方程线性相关。根因闪蒸/平衡级模型中归一化方程与总物料衡算、逐组分衡算之间存在可导出关系本篇代码 2-1 已给证明。解法在模型头里写明两条归一化中仅一条独立不要盲目删掉归一化——删掉会让模型真正欠定正确做法是保留并把冗余关系记录在案。报错 3-5物性参数缺失或超出适用范围外推失效。现象报错提示某物性参数不可用或模型在某个工况下突然失败。根因该组分/二元对的参数未在参数集中提供或解出的温度超出参数回归温区。解法按物性消费清单逐条核对参数完整性对超范围工况改为适用范围更宽的方法或在结论中显式声明物性假设的边界并禁止用超范围结果做工程判断。报错 3-6求解时间异常增长。现象加了物性之后求解时间成倍增加。根因物性计算被高频调用且属性本身非线性强。解法分析阶段先用常物性做趋势判断正式计算再打开全物性把物性简化假设与正式计算假设分别记录铁律 9两者结果必须对比留痕。四、动手练习练习 1自由度与冗余为代码 2-1 做完整自由度核对要求写出归一化冗余的三行证明对 N 2。判定证明步骤完整能说明为什么不能简单删掉 E4模型求解成功且Σx、Σy均等于 1数值误差范围内。练习 2物性冻结二分法把代码 2-1 中的K(i)临时替换为两个常数取一个合理工作点值重新求解并与全物性版本对比x(i)、y(i)、L、V。判定记录两组结果的差异幅度能据此判断该体系对相平衡的敏感程度并把差异写入工程记录铁律 9。练习 3物性消费清单落地为你的课题里的模型按代码 2-2 格式写一份物性消费清单。判定至少 5 行每行都有依赖变量、参数来源、适用范围三列且非空换物性方法的三步核对动作能逐条执行。练习 4范围越界观察把代码 2-1 的T从示例值逐步提高到清单中的适用范围之外观察求解行为变化。判定能记录到至少一种异常不收敛、越界、时间暴涨能说明外推失效与物性方法选错在现象上的区别。练习 5思考题无标准答案为什么热力学一致性比单点物性精确更优先请用残差可以下降但结论错误这一现象给出两条理由。五、小结与下一篇预告本篇把物性从模型里的常数还原为方程组里的一半物性是被联立进方程的一批约束它的失真不会破坏物料衡算的守恒性只会把解推到守恒但错误的位置模型与物性之间的抽象层属性调用让你能换物性而几乎不改模型前提是把模型对物性的全部依赖显式登记相平衡与焓必须来自一致的热力学链条一致性优先于单点精确物性参数必须回答值多少、从哪来、在什么范围内有效三问铁律 10闪蒸类模型里两条归一化仅一条独立是必须写进模型头的结构性事实铁律 8。最后参数估计时一次只估一类参数——物性参数用独立物性数据先定动力学参数留给过程数据第 11 篇。第 08 篇进入动态世界第一个动态模型——从LUMPED到DISTRIBUTED。我们讲清动态与稳态的本质差别时间自变量、初值、微分方程、集中参数LUMPED与分布参数DISTRIBUTED模型的分野对应官方能力名Implementing a lumped model in gPROMS与Implementing a distributed model in gPROMS、空间离散与网格/阶数选择对精度与耗时的影响以及动态仿真对初值的敏感性——那是从会算稳态到会算过程的分水岭。本篇认知问题回显FAQQ1物性在模型里的位置到底在哪为什么它不是常数A只有纯组分基础物性分子量、临界参数、生成焓近似常数可放进 PARAMETER 段密度、黏度、热容随温度乃至组成与压力变化相平衡常数、逸度系数、活度系数则是 P、T 与全组成的强非线性函数反应焓也随温度变化。因此物性不是查询表而是被联立进方程组的约束。它失真的后果不是报错而是非线性结构悄悄改变例如唯一稳态变多重稳态、收敛变振荡所以排查时要先排除物性。Q2模型与物性之间的抽象层是如何工作的为什么能换物性而不改模型A抽象层让模型只做属性调用即把当前的 T、P、组成作为输入、取回属性值由物性系统负责具体计算。这带来三个好处换物性方法时模型一字不改同一模型跑多种热力学假设成为低成本实验物性参数集集中维护避免同一物性在不同模型里取值不同属性调用点即模型对物性的全部依赖便于逐条追溯审查。属性的具体名称、物性方法选择、参数集覆盖语法与外部物性接入方式均以官方《gPROMS ModelBuilder Guide》为准。Q3相平衡与热力学一致性为什么会直接影响结论不一致有哪些症状A因为物料衡算无论物性对错都永远成立物性错误不会被守恒方程拦住只会把解推向守恒但物理上错的位置。热力学一致性要求相平衡关系与焓、熵等属性相互自洽。不一致的三种典型症状是摩尔分数越界甚至为负归一化被强行满足而某组分被挤成负值能量衡算残差始终不降因为焓与相平衡来自不同推导链同一体系在不同初值下收敛到不同解出现伪多重稳态。Q4物性参数与文献数据要遵循什么纪律A按铁律 10物性参数必须取自可引用文献或官方示例并标注来源不得凭记忆填数。执行上每个参数都要回答三问值是多少、从哪来、在什么范围内有效。第三问最常被忽略因为动态模型里的温度往往是被解出来的变量无法事先保证落在参数回归温区内把模型跑到温区之外时所谓不收敛其实是外推失效。此外所有物性依赖应按清单登记并在物性方法或参数集变更后按铁律 9 用固定算例复跑对比基线。Q5物性错误会造成哪些典型失败参数估计时物性参数怎么处理A典型失败包括相平衡不收敛或残差抖动、摩尔分数越界、能量衡算残差不降、诊断提示方程冗余而模型实际可解、物性参数缺失或超出适用范围、求解时间异常增长。参数估计的处置是三步先用独立的物性数据纯组分性质、二元相平衡数据单独确定物性参数不放进动力学回归的决策变量集合再在物性固定的前提下用过程数据回归动力学参数最后把两步结果放回同一模型核对能否复现全部数据。原则是一次只估一类参数。本文涉及的语言细节以所购版本随附的官方 gPROMS 文档为准。