1. 闪蒸与相平衡基础
闪蒸是化工流程模拟中最基础、也最常被低估的单元操作。它的本质可以一句话概括:一股进料在给定条件下被部分汽化(或部分冷凝),汽液两相在分离器中达到平衡后分别排出。与精馏不同,闪蒸不追求高纯度分离,而是回答一个更基本的问题——在指定温度、压力下,物料会分成多少汽、多少液,各自组成是什么。Aspen Plus 中的 Flash2、Flash3、Valve 等模块,以及几乎所有塔器的单级计算内核,都建立在同一套相平衡框架之上。
一、相平衡条件与气液平衡准则
当汽液两相接触足够充分、接触时间足够长,两相之间净的物质传递停止,系统即达到相平衡。热力学上,相平衡需同时满足三个条件:热平衡(两相温度相等,TV=TL)、力学平衡(两相压力相等,pV=pL)、化学势平衡(各组分在两相中的化学势相等,μiV=μiL)。工程上最常用的是化学势平衡的等价形式——逸度相等:
f_i^V = f_i^L
即 y_i · φ_i^V · p = x_i · γ_i · f_i^L(纯)
常用简化形式: y_i = K_i · x_i , K_i = γ_i · p_i^sat / (φ_i^V · p)
其中 Ki 为相平衡常数。对低压理想体系,γi≈1、φiV≈1,Ki 退化为 pisat/p,即拉乌尔定律;对高压轻烃体系,则需用状态方程法(如 PR、SRK)同时计算两相逸度系数。Aspen Plus 中这一选择体现在物性方法上:NRTL/UNIQUAC + 理想气体适用于低压极性体系,PR/SRK适用于中高压烃类,NRTL-RK 则是常见的混合策略。物性方法选错,后面所有闪蒸结果都不可信,这是初学者最容易踩的坑。
二、自由度分析:闪蒸到底需要几个独立变量
用吉布斯相律分析一个 C 组分、两相共存的闪蒸罐:F = C − P + 2 = C − 2 + 2 = C。也就是说,对 C 个组分,只需指定 C 个独立强度变量,系统状态即被确定。工程上通常把这 C 个变量组织为:进料组成(C−1 个摩尔分数)、温度、压力,共 C+1 个;但闪蒸计算中进料组成和流量是已知输入,真正需要"再指定"的是温度与压力中的哪一个、以及另一个如何被约束——这正是区分闪蒸类型的关键。
三、两类闪蒸的差别与工业场景
等温闪蒸(Isothermal Flash):给定进料、温度 T 和压力 p,求解汽化率 β 及两相组成。核心方程是 Rachford-Rice 方程:
Σ [ z_i (K_i − 1) / (1 + β (K_i − 1)) ] = 0 , β = V / F
该方程对 β 单调,用牛顿法或割线法迭代即可收敛,再配合 Ki 的内外循环更新。典型工业场景:部分冷凝器、闪蒸罐、脱气罐、压缩机级间分液。例如某轻烃进料在 30 ℃、0.8 MPa 下等温闪蒸,可快速判断气相分率与露点余量。
绝热闪蒸(Adiabatic Flash):给定进料、压力 p,并令闪蒸罐绝热(Q=0),温度由能量平衡自行确定。此时未知量多了一个 T,方程也多了一个焓平衡:
H_F = β · H_V(T,p) + (1 − β) · H_L(T,p)
求解时需在 T 与 β 之间做双层迭代:内层对给定 T 解 Rachford-Rice,外层调整 T 使焓平衡残差为零。典型场景:节流阀后闪蒸、塔顶回流罐、反应器出口急冷分液、天然气露点控制。工程上最典型的例子是高压液体经阀门节流到低压,压力骤降导致部分汽化、温度自降,这正是"闪蒸"一词的来源。
四、在 Aspen Plus 中如何落地
以 Flash2 模块为例,实际操作路径是:①在 Properties 中选定物性方法(如 NRTL-RK);②在 Setup 中输入组分;③连接进料流并给定流量、组成、T、p;④在 Flash2 的 Specification 页选择"Temperature"或"Pressure"作为第二指定变量——若选 Temperature 并给定值,即为等温闪蒸;若把 Pressure 指定、温度留空并保持默认绝热设置,则为绝热闪蒸;⑤运行后查看结果中的 Vapor Fraction、两相组成与热负荷 Q。若 Q 不为零而本应绝热,说明指定方式有误。
判断结果是否合理,可看三点:β 是否落在 0~1 之间(若为 0 或 1,说明处于泡点/露点之外,单相区);Σyi 与 Σxi 是否均为 1;绝热闪蒸的出口温度是否介于进料温度与露点温度之间。这些检查能快速暴露物性方法或输入错误。