1. 课程导论与模型定位
在化工过程模拟、精馏与萃取工艺设计、共沸物分离以及溶剂筛选等工作中,只要体系偏离理想溶液,就离不开活度系数模型。NRTL(Non-Random Two-Liquid,非随机双液体模型)与 UNIQUAC(Universal Quasi-Chemical,通用准化学模型)是工程上使用频率最高的两个局部组成型活度系数模型。它们既不像 Wilson 模型那样无法描述液液分层,也不像 UNIFAC 那样完全依赖基团贡献而缺少二元交互参数的灵活性,因此在极性体系、部分互溶体系以及汽液平衡(VLE)与液液平衡(LLE)计算中占据核心位置。本课程不打算逐条罗列公式推导,而是以“工况对照”为主线,帮助你在面对具体体系时判断:该选哪一个、参数从哪来、算出来的结果能不能信。
一、活度系数模型谱系中的位置
理解 NRTL 与 UNIQUAC 的定位,先要看清整个模型谱系的分层逻辑。最底层是理想溶液假设(γ=1),只适用于分子大小相近、相互作用相似的非极性体系,如苯—甲苯。往上一层是正规溶液与 Scatchard-Hildebrand 类模型,用溶解度参数描述非理想性,但难以处理强极性与缔合。再往上是局部组成模型族,包括 Wilson、NRTL、UNIQUAC,它们承认“局部组成不等于整体组成”这一核心物理图像,因而能描述强非理想性。最上层是基团贡献法(UNIFAC、modified UNIFAC),通过基团拆分预测缺少实验数据体系的活度系数。
NRTL 与 UNIQUAC 同属局部组成模型,但二者与 Wilson 的关键区别在于:Wilson 模型在二元部分互溶体系中无法给出液液分层,而 NRTL 与 UNIQUAC 都能同时描述 VLE 与 LLE。这一能力差异直接决定了它们在萃取、共沸精馏、水—有机溶剂体系中的不可替代性。相对 UNIFAC,二者又属于“有参数模型”——需要由实验数据回归二元交互参数,精度更高但依赖数据基础。
二、共同假设与核心差异
两者共享同一物理出发点:液体混合物在微观尺度上并非均匀,围绕某一分子的局部组成与整体组成不同。这一“局部组成”概念使它们能刻画氢键、极性相互作用引起的强非理想性。同时,二者都建立在过量吉布斯自由能 gE 的表达之上,通过 gE/RT 对组成的偏导得到活度系数,因此都具备热力学一致性,可由二元参数推算多元体系。
差异主要体现在三处。第一,结构描述方式:UNIQUAC 把 gE 拆成组合项(combinatorial)与剩余项(residual)。组合项只依赖纯组分结构参数 r(体积参数)和 q(表面积参数),描述分子大小与形状差异带来的熵效应;剩余项才包含二元交互参数 τij。NRTL 则没有独立的组合项,其非理想性全部由 τij 与 αij 承担。第二,参数体系:NRTL 每个二元对需要 τ12、τ21 和 α12 三个参数,α 常取 0.2~0.47 的经验值;UNIQUAC 每个二元对只需 τ12、τ21 两个交互参数,r、q 由分子结构查得或由基团贡献估算。第三,适用倾向:NRTL 对强极性、部分互溶体系(如水—醇、水—酯)拟合灵活,α 参数提供了额外的形状调节能力;UNIQUAC 因组合项显式处理分子大小差异,在聚合物溶液、分子尺寸悬殊的体系以及需要外推至多元时表现稳健,且参数更少、回归更稳定。
三、以“工况对照”为主线的学习地图
本课程后续 29 章将围绕一条主线展开:给定体系与工况,如何判断模型适用性并落地计算。学习地图可概括为四个递进层次。
第一层是参数层。你需要掌握二元交互参数的来源与质量判断:是实验回归、文献查取还是 UNIFAC 估算;参数的温度适用范围;以及 αij 取值的敏感性。示例性地,若某水—乙醇体系在 101.3 kPa 下回归得到 τ12=0.85、τ21=−0.42、α=0.30(此处数值为教学假设,非实测),则应检查该参数组在 60~90 ℃ 区间内预测的 VLE 是否与实验点吻合,而不是直接套用到减压工况。
第二层是工况层。同一模型在 VLE、LLE、VLLE 以及高压、稀溶液、多组分条件下的表现不同。例如 NRTL 在液液分层体系中需保证 α 与 τ 组合能产生两个液相根;UNIQUAC 在稀溶液端对 r、q 的微小误差较敏感。课程会用对照表与判据帮助你在具体工况下做取舍。
第三层是计算与验证层。包括活度系数计算、泡露点迭代、参数回归目标函数设置,以及用实验数据做残差分析。你会看到如何用少量二元数据回归 τ 参数,并用另一组独立数据做验证,而不是用同一组数据既回归又验证。
第四层是工程决策层。在精馏塔设计、萃取剂筛选、共沸分离方案比较中,模型选择直接影响塔板数、回流比与溶剂比。课程将给出典型体系的对照结论与常见误用排障路径。
| 对照维度 | NRTL | UNIQUAC |
|---|---|---|
| 二元参数个数 | τ12、τ21、α12(3 个) | τ12、τ21(2 个)+ 纯组分 r、q |
| 结构描述 | 无独立组合项,靠 α 调节 | 组合项显式处理分子大小与形状 |
| 典型优势工况 | 强极性、部分互溶、LLE/VLLE | 分子尺寸悬殊、聚合物溶液、多元外推 |
| 主要风险点 | α 取值不当导致分层预测失真 | r、q 估算误差在稀溶液端放大 |