气液平衡模拟器 — 拉乌尔定律 返回
化学工程

气液平衡模拟器 — 拉乌尔定律

当二元混合物的液体和蒸气达到平衡时,各相的组成和温度如何确定?本工具通过拉乌尔定律和安托万方程实时计算沸点、露点、平衡组成和相对挥发度,并用T-x-y线图和y-x平衡线图直观展示蒸馏为什么能分离成分。

参数设置
二元系(轻质-重质)
自动设置安托万常数。成分1为轻质(易挥发)成分
液相轻质成分摩尔分率 x₁
沸点(气泡点)计算的输入组成
蒸气相轻质成分摩尔分率 y₁
露点(露点)计算的输入组成
全压 P
kPa
系统的操作压力。101.3 kPa 为标准大气压
计算结果
沸点(气泡点)T_b (℃)
平衡蒸气组成 y₁(与 x₁ 平衡)
露点(露点)T_d (℃)
平衡液组成 x₁(与 y₁ 平衡)
相对挥发度 α
气液浓度差 y₁ − x₁
沸点线图(T-x-y)— 操作点和连接线

下方是气泡点曲线,上方是露点曲线。两者之间是气液二相共存区。闪烁的点是当前操作点,水平线连接处于平衡状态的液和蒸气。

T-x-y 线图 — 沸点和露点曲线
y-x 平衡线图 — 平衡曲线和对角线
理论与主要公式

$$p_i=x_i\,P_i^{sat}(T),\qquad \sum_i x_i P_i^{sat}=P\ \ (\text{沸点})$$

拉乌尔定律。各成分的分压 p_i 是液相摩尔分率 x_i 与纯组分蒸气压 P_i^sat 的乘积。分压之和等于全压的温度就是沸点(气泡点)。

$$\log_{10}P^{sat}=A-\frac{B}{C+T},\qquad y_1=\frac{x_1 P_1^{sat}}{P}$$

纯组分蒸气压用安托万方程作为温度的函数求得(T:℃、P^sat:mmHg、常数A、B、C为成分固有)。平衡蒸气组成 y₁ 直接由拉乌尔定律得出。

$$\alpha=\frac{P_1^{sat}}{P_2^{sat}}$$

相对挥发度 α 是轻质和重质成分蒸气压的比值。蒸气总是富集易挥发的成分(轻质成分),α离1越远,蒸馏分离就越容易。

气液平衡与拉乌尔定律概述

🙋
「气液平衡」这个词在化学工程里经常出现…就是液体和蒸气混合后达到平衡的状态吗?
🎓
没错,你的理解是对的。把苯和甲苯的混合液放在密闭容器里,保持温度恒定,液面上方会积累蒸气。蒸发和凝聚的速度最终达到平衡,液体和蒸气的组成都不再变化——这就是气液平衡。问题是「液体是50:50时,上面的蒸气会是几比几呢」「什么温度开始沸腾」。拉乌尔定律就是来回答这些问题的。
🙋
拉乌尔定律好像是 p = x·P 这样的公式。那是什么意思呢?
🎓
简单说,就是「每个成分产生的蒸气压力(分压)等于该成分在液体中的比例(摩尔分率 x)乘以它作为纯液体时的蒸气压 P^sat」。公式是 p_i = x_i·P_i^sat。如果苯在液体中占一半,那么苯的分压就是纯苯蒸气压的一半。全压是各成分分压之和。对理想混合物,用这个公式就能解气液平衡。
🙋
纯组分的蒸气压 P^sat 会随温度变化吧。怎么计算呢?
🎓
那就要用安托万方程。形式是 log₁₀P^sat = A − B/(C+T),其中A、B、C是从实验数据确定的成分特定常数。输入温度T,就能算出纯组分蒸气压。这个工具用二分法找满足 Σx_i·P_i^sat = P 的温度T。那个温度就是沸点,也叫气泡点。你拖动 x₁ 滑块时,沸点会在苯的沸点80℃和甲苯的110℃之间平滑变化。
🙋
明白了。那为什么平衡的蒸气总是偏向轻质成分呢?我看结果卡上 y₁ 总是比 x₁ 大。
🎓
很好的问题。轻质成分(这里是苯)在同一温度下的纯组分蒸气压比甲苯高。根据拉乌尔定律,y₁ = x₁·P1^sat/P,P1^sat 越大,y₁ 相对于 x₁ 的增幅越大。也就是说蒸气必然富集「容易挥发的那个」。正是这种「每次气化都浓缩轻质成分」的现象才是蒸馏的核心。把它段段重复——气化、凝聚、再气化……就能从50:50的原料提出99%的苯。下面的y-x平衡线图里,平衡曲线远高于对角线 y=x,就是这个「每次气化的浓缩效果」。
🙋
还有个「相对挥发度 α」的数字。那是干什么用的?
🎓
可以把 α 看作分离「难度表」。α = P1^sat/P2^sat,就是两个成分蒸气压的比。α是2或3的话,平衡曲线离对角线很远,用少数几块塔板就能分离。α只有1.1左右的话,曲线紧贴对角线,需要几十块塔板才行。实际工作中,「α跌到1.05以下就别用普通蒸馏了,考虑萃取蒸馏或别的办法」。当α恰好等于1时,那就是共沸点,普通蒸馏永远分不了。试试切换不同的二元系,看 α 怎么变。

常见问题

拉乌尔定律是理想溶液的基本定律:各成分的分压 p_i 等于该成分的液相摩尔分率 x_i 与其纯液体蒸气压 P_i^sat 的乘积(p_i = x_i·P_i^sat)。全压 P 为各成分分压之和。本工具使用安托万方程将纯组分蒸气压表示为温度的函数,计算苯-甲苯等分子结构相似混合物的沸点、露点和平衡组成。对于化学性质相近的成分混合物,这个简单的理想模型可以给出实用精度。
沸点(气泡点)是加热某组成液体时第一个气泡出现的温度,由条件 Σx_i·P_i^sat = P 确定。露点(露点)是冷却某组成蒸气时第一滴液体出现的温度,由条件 Σy_i/P_i^sat = 1/P 确定。在相同的全压下,除非液相和蒸气相组成完全相同,否则沸点和露点不相等。在T-x-y线图中,下方的曲线是沸点,上方是露点,两者之间是气液二相共存区。
相对挥发度 α 是轻质和重质成分纯组分蒸气压的比值 α = P1^sat/P2^sat,表示蒸馏分离的难易程度。α越大,平衡蒸气与液体中轻质成分的差异越大,单次气化凝聚的组成变化越大,用更少的塔板就能达到高纯度分离。反之,α接近1的混合物(沸点接近的成分)分离困难,需要更多的塔板或更高的回流比。当α恰好等于1时为共沸点,普通蒸馏无法分离。
轻质成分在相同温度下具有更高的纯组分蒸气压。根据拉乌尔定律,平衡蒸气组成由 y_1 = x_1·P1^sat/P 给出,P1^sat 越大,y_1 相对于 x_1 的增幅越大。这意味着气相必然富集易挥发的成分。正是这种「每次气化都浓缩轻质成分」的现象构成了蒸馏的原理,通过多次气化和冷凝,可以从原料中获得远高于初始纯度的留出液。本工具的y-x平衡线图显示平衡曲线总是在对角线 y=x 上方,正说明了这种浓缩效应。

实际应用

蒸馏塔设计:石油精炼和石油化学工厂的蒸馏塔无法离开气液平衡计算进行设计。原油分馏成轻汽油、柴油、燃料油等,苯和甲苯的分离,都是基于相同的「每次气化都浓缩轻组分」的原理,经过多塔板的累积。设计时要逐段求解气液平衡,确定所需的理论塔板数和回流比。相对挥发度α直接决定分离难度,因此初期设计评估时首先要确认α。

溶剂回收和精制:涂料、印刷、半导体制造等行业的废溶剂回收和再利用工艺中,蒸馏精制被广泛使用。甲苯、己烷这类烃类溶剂往往混含相似的组分,表现出接近理想的气液平衡,因此基于拉乌尔定律的计算方法有效。根据回收率和纯度目标,用这种平衡计算来决定塔的运行条件。

工艺模拟器的基础:Aspen Plus、Pro/II这样的商业工艺模拟软件,内部反复进行大量的气液平衡计算。理想溶液用拉乌尔定律,非理想系用活度系数模型(NRTL、UNIQUAC等)或状态方程来分别处理。本工具处理的理想模型是理解这些高级模型的出发点,也可用作模拟结果的粗略检验。

教学与概念理解:学习化学工程时,「会读」T-x-y线图和y-x平衡线图是基础中的基础。沸点曲线和露点曲线围成的二相区、连接线表示的液蒸气组成对应、平衡曲线与对角线的距离代表分离能力——不仅要从公式上理解,更要从图上直观感受。这种理解为后来学习McCabe-Thiele法的塔板数估算和蒸馏塔设计打下基础。

常见误区与注意事项

首先最大的误区是「拉乌尔定律对所有混合物都适用」。拉乌尔定律仅对理想溶液成立,即分子大小和种类相似、异种分子间作用力与同种分子间相近的混合物。本工具的苯-甲苯、正戊烷-正己烷确实接近理想,但乙醇-水这类涉及氢键的系统会严重偏离理想,甚至出现共沸(α=1,普通蒸馏无法分离)。非理想系需要引入活度系数 γ,改用修正的拉乌尔定律 p_i = γ_i·x_i·P_i^sat,本工具的结果不能直接套用。

其次,「混淆沸点和露点」。同一个数值,若作为液相组成 x₁ 读解,算出的温度是沸点;若作为蒸气组成 y₁ 读解,算出的温度是露点——两者完全不同。x₁=0.5 的液体加热到冒泡的温度(沸点),和y₁=0.5 的蒸气冷却到结露的温度(露点)是两码事。本工具为沸点和露点分别设了输入滑块(x₁ 和 y₁),正是为了区分。在T-x-y线图中,下方曲线(沸点)用液相组成的横轴读,上方曲线(露点)用蒸气组成的横轴读,千万别搞反。

最后,「提高压力就一定分离更容易」的想法是错的。提高压力会改变沸点和露点的温度,但相对挥发度α也会通过温度依赖性而改变。一般来说,升高压力(因而升高温度)时α会趋向1,分离反而变难。减压蒸馏的优势在于不仅能在低温下处理热敏物质,还能增大α。选择操作压力要综合考虑分离难度、运行温度、装置成本和冷却加热资源的可用性,不能单一追求高压。本工具的全压参数可调,你可以边改压力边观察α的变化,体会这种敏感性。

使用指南

  1. 用滑块设定液相组成 x₁ (0~1),输入全压 (1~101 kPa)
  2. 点击「计算沸点」,根据拉乌尔定律 P = x₁P₁ᵛᵃᵖ + (1-x₁)P₂ᵛᵃᵖ 计算出气泡点 T_b
  3. 由蒸气组成 y₁ = x₁P₁ᵛᵃᵖ / P_total 自动计算平衡蒸气组成,气液平衡曲线绘制在T-x-y线图上
  4. 从蒸气组成 y₁ 反算露点 T_d 和平衡液组成 x₁,显示相对挥发度 α = (y₁/y₂)/(x₁/x₂)

具体计算示例

苯-甲苯混合液在常压 101 kPa、x₁(苯) = 0.4 时:由安托万方程在80℃处得 P₁ᵛᵃᵖ = 135 kPa、P₂ᵛᵃᵖ = 47 kPa,气泡点 T_b ≈ 89℃,平衡蒸气组成 y₁ = 0.4×135/(101) ≈ 0.535。此时相对挥发度 α ≈ 2.6,气液浓度差 y₁ - x₁ = 0.135 能评估分离效率。

实务中的注意点