3. 拉乌尔定律与理想溶液气液相平衡

在气液相平衡计算中，最基础也最常用的一类模型是理想溶液模型。它的核心假设是：溶液中各组分分子间的相互作用与纯组分内部完全相同，混合过程既不放热也不吸热，体积具有加和性。在这一假设下，组分 i 在气相中的分压只取决于它的液相摩尔分数和该温度下的纯组分饱和蒸气压，这就是拉乌尔定律：

p_i = x_i · p_i^sat(T)

其中 p_i 为气相中组分 i 的分压，x_i 为液相摩尔分数，p_i^sat 为同温度下纯组分 i 的饱和蒸气压。对理想气体还叠加道尔顿分压定律 p_i = y_i·P，于是得到理想溶液的相平衡关系：

y_i · P = x_i · p_i^sat(T)   ⟹   y_i = (p_i^sat/P) · x_i

这个式子说明：在给定温度和总压下，气相组成与液相组成成正比关系，比例系数就是该组分的饱和蒸气压与总压之比。饱和蒸气压越高的组分（更易挥发），在相同液相组成下气相摩尔分数越大，这正是精馏分离的物理基础。

泡点方程与露点方程

对二元理想体系，两组分满足 x₁+x₂=1、y₁+y₂=1。把拉乌尔定律代入归一化条件，可得到两个最常用的方程。

泡点方程：已知液相组成 x，求体系在给定总压下开始沸腾的温度（泡点）及与之平衡的气相组成。由 P = x₁p₁^sat + x₂p₂^sat，由于 p_i^sat 是温度的函数，该式对 T 是隐式的，通常用试差或迭代求解：

给定 x1, x2=1-x1, 总压 P
猜一个温度 T
计算 p1sat(T), p2sat(T)   // 如用安托因方程
Pcalc = x1*p1sat + x2*p2sat
若 |Pcalc - P| > 容差, 调整 T 重算
收敛后: y1 = x1*p1sat/P,  y2 = 1 - y1

露点方程：已知气相组成 y，求给定总压下开始冷凝的温度（露点）及平衡液相组成。由 1/P = y₁/p₁^sat + y₂/p₂^sat 对 T 试差，收敛后 x₁ = y₁P/p₁^sat。同一总压下，同一组成的露点总高于泡点，两者之间就是气液两相共存区。

相对挥发度

工程上更常用相对挥发度 α 描述分离难易程度，定义为两组分挥发度之比。对理想体系，α 可简化为饱和蒸气压之比：

α₁₂ = (y₁/x₁) / (y₂/x₂) = p₁^sat / p₂^sat

结合 y₁+y₂=1，可导出只含 α 的相平衡方程：

y₁ = α x₁ / [1 + (α−1)x₁]

α 越大，相同液相组成下气相富集越明显，越容易分离；α 接近 1 时，两组分挥发度相近，普通精馏难以分开，需要特殊分离手段。α 随温度变化，但在一定温度范围内常取几何平均值近似为常数，这使平衡曲线绘制大为简化。

示例估算：苯—甲苯体系

苯—甲苯是典型的近似理想二元体系，常压（P=101.3 kPa）下可作为教学算例。设某温度下查得 p_苯^sat=120 kPa、p_甲苯^sat=40 kPa（示例假设值，仅用于演示算法）。

若液相组成 x_苯=0.4，则泡点压力 P = 0.4×120 + 0.6×40 = 72 kPa，对应气相组成 y_苯 = 0.4×120/72 = 0.667。可见气相中苯的摩尔分数明显高于液相，说明苯更易挥发。此时相对挥发度 α = 120/40 = 3.0，代入平衡方程 y = 3×0.4/[1+2×0.4] = 1.2/1.8 = 0.667，与分压法结果一致，可互相校核。

若已知气相组成 y_苯=0.667、总压 72 kPa，求露点液相组成：x_苯 = y_苯P/p_苯^sat = 0.667×72/120 = 0.400，正好回到原液相组成，验证了泡点与露点计算互为逆过程。

掌握上述三条主线——拉乌尔定律给出分压、泡点/露点方程完成组成与温度的互算、相对挥发度把平衡关系压缩为单一参数——就具备了理想二元体系相平衡数据估算的基本能力，也为后续非理想体系的修正与相图识读打下基础。