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 个;但闪蒸计算中进料组成和流量是已知输入,真正需要"再指定"的是温度与压力中的哪一个、以及另一个如何被约束——这正是区分闪蒸类型的关键。

记忆要点:等温闪蒸指定 T 和 p;绝热闪蒸指定 p 和"绝热"(即焓守恒,相当于用能量方程替代温度指定)。两者都是 C 个独立约束的不同组合。

三、两类闪蒸的差别与工业场景

等温闪蒸(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 使焓平衡残差为零。典型场景:节流阀后闪蒸、塔顶回流罐、反应器出口急冷分液、天然气露点控制。工程上最典型的例子是高压液体经阀门节流到低压,压力骤降导致部分汽化、温度自降,这正是"闪蒸"一词的来源。

闪蒸罐物料与相平衡关系示意 进料 F, z_i, T_F, p_F, H_F → 汽相 V, y_i, T, p, H_V + 液相 L, x_i, T, p, H_L 进料 F z_i, T_F, p_F 焓 H_F 闪蒸罐 T, p 平衡 y_i = K_i · x_i 汽相区 V 液相区 L β = V / F 汽相 V y_i, H_V 液相 L x_i, H_L 等温闪蒸 指定:T、p(进料已知) 求解:β、y_i、x_i 核心:Rachford-Rice 方程 Σ z_i(K_i−1)/(1+β(K_i−1)) = 0 场景:部分冷凝器、闪蒸罐、级间分液 绝热闪蒸 指定:p、绝热(Q=0) 求解:T、β、y_i、x_i 核心:焓平衡 + Rachford-Rice H_F = β·H_V(T,p) + (1−β)·H_L(T,p) 场景:节流阀后闪蒸、回流罐、急冷分液 相平衡三条件: T^V = T^L(热平衡) p^V = p^L(力学平衡) f_i^V = f_i^L(化学势/逸度平衡) 自由度:F = C − P + 2 = C(两相),故需 C 个独立约束;等温闪蒸用 T、p,绝热闪蒸用 p 与焓守恒。

四、在 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;绝热闪蒸的出口温度是否介于进料温度与露点温度之间。这些检查能快速暴露物性方法或输入错误。

常见误区:把"绝热"理解为"不指定温度就自动正确"。实际上 Flash2 默认可能带热负荷指定,务必确认 Q=0 或选择 Adiabatic 选项,否则算出的温度没有物理意义。