4. 等温闪蒸模拟实战

等温闪蒸是化工流程模拟中最基础也最常用的单元操作模型之一。它的物理含义非常明确:一股进料在给定温度 T 和压力 P 下进入闪蒸罐,在罐内达到气液两相平衡后分离为气相产品 V 和液相产品 L。Aspen Plus 中承担这一计算的是 Flash2 模块,其自由度规定为:进料状态(流量与组成)已知,再指定温度与压力两个独立变量,气化率、两相组成和相分率全部由相平衡方程求解得到。本节以正戊烷/正己烷二元混合物为例,完整走一遍建模、参数设置、灵敏度分析与结果解读的流程。

一、物系与进料设定

教学示例采用正戊烷(n-C5)与正己烷(n-C6)各 50 mol% 的二元混合物,总进料量 100 kmol/h,进料温度 60 ℃、压力 1 atm。物性方法选择 Raoult 定律(IDEAL 或 Raoult 选项),因为该体系在低压下接近理想溶液,用 Raoult 定律可以手算校核,便于理解 Flash2 内部到底在算什么。若体系含极性组分或操作压力较高,则应改用 NRTL、PR 等活度系数或状态方程方法,这一点在后续章节展开。

在 Aspen Plus 中依次建立:组分定义(C5H12、C6H14)→ 物性方法(IDEAL)→ 模拟环境放置 Flash2 模块(命名为 FLASH1)→ 连接进料流 FEED、气相产品 VAP、液相产品 LIQ。注意 Flash2 是两出口模块,必须同时连接气相与液相出口,否则模块无法完成相分离计算。

二、Flash2 参数设置要点

双击 FLASH1 打开输入面板,在 Specification 页中,闪蒸规格(Flash specifications)选择 Temperature 与 Pressure 两个变量,分别填入 60 ℃ 和 1 atm。这是等温闪蒸的标准设置;若改为 Pressure 与 Vapor fraction,则变成等压气化率指定型闪蒸,属于另一种常见工况。有效相态(Valid phases)保持默认的 Vapor-Liquid,若体系可能出现液液分层,才需要改为 Vapor-Liquid-Liquid。

收敛页一般无需改动,Flash2 采用序贯求解,单级闪蒸通常 3~5 次迭代即可收敛。需要提醒的是:当指定温度落在泡点与露点之间时,闪蒸有唯一的两相解;若指定温度低于泡点,计算结果会退化为全液相(气相流量为 0),高于露点时则退化为全气相。这不是计算错误,而是热力学上确实不存在两相区,解读结果时必须先判断这一点。

三、灵敏度分析:压力与温度对气化率的影响

气化率(Vapor fraction)是等温闪蒸最核心的输出指标。用 Aspen Plus 的 Sensitivity 工具建立两组分析:第一组固定温度 60 ℃,让压力在 0.5~3 atm 之间变化;第二组固定压力 1 atm,让温度在 40~90 ℃ 之间变化。被测量变量(Flowsheet variable)选择 FLASH1 的 Vapor fraction,以及 VAP 流股的总摩尔流量。

从热力学上可以预判结果趋势:对同一进料,压力升高会抑制气化,气化率下降;温度升高会促进气化,气化率上升。在 1 atm 下,该二元混合物按 Raoult 定律估算的泡点约 55 ℃、露点约 66 ℃(示例估算值,实际以模拟结果为准),因此 60 ℃ 恰好落在两相区内,气化率约为 0.4~0.5 量级。当温度升到 70 ℃ 以上时,气化率迅速趋近 1;降到 50 ℃ 以下时趋近 0。灵敏度曲线在泡点与露点附近斜率最大,这正是相变区间的典型特征。

等温闪蒸 Flash2 计算流程与气化率响应 FEED 进料 100 kmol/h C5/C6 各 50% FLASH2 T = 60 ℃ P = 1 atm V-L 平衡 VAP 气相 y_C5 富集 LIQ 液相 x_C6 富集 自由度:进料 + T + P → 求 V、L、y、x 气化率 vs 温度 40 90℃ 1.0 0 60℃ 灵敏度分析:压力对气化率的影响(固定 T = 60 ℃) 0.8 0.4 0 0.72 0.5 atm 0.46 1.0 atm 0.28 1.5 atm 0.15 2.0 atm 0.05 3.0 atm 压力升高 → 气化率单调下降(示例估算趋势)

四、结果解读与相图绘制

模拟收敛后,打开 FLASH1 的结果页,重点看三组数据:一是 Vapor fraction,即气化率;二是 VAP 与 LIQ 两股流股的摩尔流量,二者之和应等于进料 100 kmol/h,这是物料衡算的快速校验;三是两相的组成,气相中正戊烷摩尔分数应明显高于液相,因为正戊烷挥发度更大。以 60 ℃、1 atm 的示例工况估算,气相中 C5 约 0.62、液相中 C5 约 0.38(示例估算值),符合轻组分向气相富集的规律。

相图的绘制可以借助灵敏度分析的数据直接完成。将温度扫描得到的泡点、露点数据整理成 T-x-y 图:横坐标为 C5 的液相摩尔分数 x 与气相摩尔分数 y,纵坐标为温度。上曲线为露点线(T-y),下曲线为泡点线(T-x),两线之间的区域即两相区。60 ℃ 的水平线与两条曲线的交点,分别对应平衡液相组成与平衡气相组成,而闪蒸点就落在这条水平线(等温线)上——这正是"等温闪蒸"名称的几何含义。若改用压力扫描,则可绘制 P-x-y 图,分析方法完全对称。

判断结果是否合理的三个快速检查:① 物料衡算闭合(V + L = F);② 气化率介于 0 与 1 之间,若等于 0 或 1 说明指定条件已越出两相区;③ 轻组分在气相中的摩尔分数大于其在液相中的摩尔分数。三条同时满足,等温闪蒸结果基本可信。
常见误区:把 Flash2 的"温度-压力"规格与"压力-气化率"规格混用。前者求解的是平衡温度或平衡组成,后者求解的是达到指定气化率所需的温度。两者物理含义不同,切换规格后必须重新检查结果页中的温度或气化率是否落在合理区间,不能直接沿用旧结论。

完成上述步骤后,读者应能独立完成一个等温闪蒸的完整闭环:从物系与物性方法选择,到 Flash2 的温度压力规格设置,再到灵敏度扫描与 T-x-y 相图绘制,最后用物料衡算与组成富集规律校验结果。这套流程是后续绝热闪蒸、多级精馏与严格相平衡计算的基础,务必把每一步的物理含义与数值结果对应起来理解,而不是只记住点哪个按钮。