吉布斯相律模拟器 返回
相平衡模拟器

吉布斯相律模拟器 — 相平衡的自由度

基于吉布斯相律 F = C − P + 2 − N,由组分数 C、相数 P、额外约束数 N 实时计算相平衡的自由度 F、1 个强度量固定下的 F、2 个强度量固定下的 F、可共存的最大相数 P_max。配合 T-P 相图概念图(水)与 F-P 直线图可视化不变点、单变点、二变点。

参数设置
组分数 C
相数 P
额外约束数 N
P_max 判定用约束数

默认值为 C=2、P=2、N=0、P_max 判定用约束数=0。计算式 F = C − P + 2 − N,F < 0 时判定为物理上不可能(指定相数 P 过大)。1 个强度量固定下 F = C − P + 1,2 个强度量固定下 F = C − P。

计算结果
自由度 F
1 强度量固定下的 F
2 强度量固定下的 F
最大相数 P_max
T-P 相图概念图(水)

横轴=温度 T/纵轴=log P/蓝色=固相(F=2)/绿色=液相(F=2)/橙色=气相(F=2)/白线=相边界(F=1)/黄点=三相点(F=0)/红点=临界点(相律外)

F-P 直线图(C 与 N 固定)

横轴=相数 P/纵轴=自由度 F/蓝色=F = C − P + 2 − N 直线/红区=F < 0(不可能)/黄标=当前 (P, F)

理论与主要公式

吉布斯相律确定了相平衡条件下可独立选择的强度变量数:

$$F = C - P + 2 - N$$

$F$ 为自由度(可独立变化的强度变量数),$C$ 为独立组分数,$P$ 为共存相数,$+2$ 对应温度 $T$ 与压力 $P$ 两个强度量,$N$ 为额外固定的约束数。$F=0$ 为不变点,$F=1$ 为单变点(曲线),$F=2$ 为二变点(区域)。可共存的最大相数由下式给出:

$$P_{\max} = C + 2 - N$$

纯物质($C=1$)无约束时 $P_{\max}=3$,这正是单组分系只有一个三相点的原因。恒压条件下($N=1$)的二组分系($C=2$)满足:

$$F = C - P + 1$$

这是 T-X 相图(金属合金相图)的理论基础:液相线、固相线上 $F=1$(共存曲线),两相区内 $F=2$。

吉布斯相律模拟器是什么

🙋
"吉布斯相律"我背的是 F = C − P + 2,但"自由度 F"具体是什么意思呢?是说我可以测量温度和压力的意思吗?
🎓
好问题。F 是「不破坏相平衡的前提下可独立变化的强度变量(温度、压力、组成等)的个数」。比如纯水 C=1,单液相 P=1 时 F = 1 − 1 + 2 = 2,T 和 P 可独立移动;液气两相共存 P=2 时 F = 1,给定 T 即唯一确定饱和蒸气压。本工具默认 C=2、P=2、N=0 时 F=2,是二变点区域。
🙋
那"三相点"就是 F=0 的不变点对吧?水的三相点是 0.01°C、611 Pa,为什么能定得这么精确?
🎓
正是如此,三相点是 F=0 的不变点。输入 C=1、P=3(固液气共存)即得 F = 1 − 3 + 2 = 0,T 和 P 唯一确定。SI 单位制(旧定义)中水的三相点 273.16 K = 0.01°C 曾被用于定义开尔文温标,原因正是相律保证了它作为「自然规律的唯一点」而非测量约定。
🙋
"额外约束数 N"这个滑块是什么呢?是不是和"恒压条件下少 1 个自由度"那种说法有关?
🎓
正是。设 N=1 相当于把压力固定在 1 atm,公式变为 F = C − P + 2 − 1 = C − P + 1。例如 Fe-C 合金(C=2)常压下两相共存(P=2)时 F = 1,给定温度即唯一确定组成。这就是金属相图(T-X 图)的理论基础。N=2 时 T 和 P 全固定 F = C − P。本工具的三个 F 统计卡片会随 N 滑块独立更新。
🙋
最后,F < 0 时显示"不可能",这是说相数太多了的意思吗?
🎓
正确。最大共存相数 P_max = C + 2 − N。纯物质(C=1、N=0)最多 3 相(即三相点),二组分系最多 4 相,三组分系最多 5 相。超出此数 F 即为负值,物理上不可能。Gibbs(1875 年)的伟大之处就在于:用这个简单的不等式一举解释了周期表上所有相图的形状。在本工具中把 C 和 P 推向极端值,F-P 图中即出现红色"不可能"区域。

常见问题

吉布斯相律给出了相平衡条件下「可独立变化的强度变量数」即自由度 F 的计算式:F = C − P + 2 − N,其中 C 为独立组分数,P 为共存相数,+2 对应温度 T 与压力 P 两个强度量,N 为额外的固定约束数(例如恒压实验取 1)。本工具默认值 C=2、P=2、N=0 时 F = 2 − 2 + 2 − 0 = 2,为可独立变化 T 与 P 的二变点区域。
F = 0 称为不变点,温度与压力均唯一确定,在相图上是「点」。水的三相点(0.01°C、611.657 Pa)即此类型,三相同时共存。F = 1 称为单变点,给定温度即唯一确定压力,对应「线」上的平衡(饱和蒸气压曲线、熔化曲线、升华曲线)。F = 2 称为二变点,温度与压力可独立变化,是「区域」上单相存在的状态。
N 表示额外固定的强度变量数。实验中将压力固定在 1 atm 时取 N = 1,将温度和压力同时固定时取 N = 2,自由度分别减少 1 和 2。1 atm 固定下的二组分系(C=2)有 F = C − P + 1,最大相数为 C + 1。最大共存相数 P_max 由 F = 0 推得 P_max = C + 2 − N:C = 2 且无约束时最多 4 相,C = 1 的纯物质 P_max = 3 相(即三相点)。
对单组分系(C=1)相律化为 F = 3 − P:P=1 是二变点(区域),P=2 是单变点(线),P=3 是不变点(点)。三相点(0.01°C、611.657 Pa)是固液气三相共存的 F=0 点,在相图上唯一确定。临界点(373.946°C、22.064 MPa)是液相与气相区别消失的奇点,属于相律框架以外(连续过渡到超临界流体)。本工具将 C=1、P=3 输入即可得 F=0,确认三相点条件。

实际应用

金属合金相图设计(Fe-C 系、Al-Cu 系):常压(N=1)下二组分系(C=2)满足 F = 3 − P。Fe-C 合金(铸铁、钢)相图中奥氏体+渗碳体+液相的三相共存点(共晶 1148°C、4.3wt%C;共析 727°C、0.76wt%C)作为 F=0 的不变点出现,跨越此处组织发生剧变。本工具输入 C=2、N=1、P=3 即得 F=0,揭示共晶/共析点为何在相图上呈现一个点。热处理工艺(淬火、回火、时效硬化)的理论基础。

蒸馏与萃取塔设计(Aspen Plus、ASPEN HYSYS):化工过程中三组分系(C=3,如甲醇-乙醇-水)的常压蒸馏塔设计直接受相律约束。F = 3 − P + 1 = 4 − P,气液两相(P=2)时 F=2,温度与一种组分的组成独立,其余两组分受约束。共沸点(azeotrope)是这种自由度的退化点,相律强制存在「蒸馏边界(distillation boundary)」无法跨越。本工具切换到 C=3 即可判断何时必须采用萃取蒸馏或变压蒸馏。

地球科学与岩浆结晶:地壳硅酸盐矿物(SiO₂、Al₂O₃、CaO、Na₂O、MgO 等)的多组分系相平衡直接套用相律。N. L. Bowen 反应序列描述玄武岩浆冷却时依次结晶橄榄石→辉石→角闪石→云母的过程,由四组分 Mg-Fe-Si-O 系相律 F = 6 − P 解释。地球深部超高压下的矿物相变也用同一框架讨论。本工具中 C=4、P=2〜5 扫描即可观察结晶路径上自由度的递减。

食品科学与冷冻干燥:食品工业的冷冻干燥(freeze drying)工艺要求压力降至水的三相点(0.01°C、611 Pa)以下,使冰直接升华为水蒸气。本工具输入 C=1、P=2(固+气)即得 F=1,给定货架温度即唯一确定升华压。咖啡、疫苗、宇航食品的品质都由这个简单的相律约束支配,操作点必须远离三相点之下侧,才能「不破坏形状地脱水」。

常见误解与注意事项

最常见的误解是把 「C 解读为系统中存在的物种数」。严格地讲,C 是「独立组分数」——所列物种数减去化学反应、化学计量与电中性条件后的独立约束数。例如用 H₂、O₂、H₂O 三种描述水时,平衡式 2H₂ + O₂ ⇌ 2H₂O 提供 1 个约束,C = 3 − 1 = 2;若进一步固定化学计量初始组成 C = 1。本工具以「独立组分数」处理 C,复杂反应系应先整理化学计量再输入。

第二个常见混淆是 「相律与热力学第二定律无关」。事实上相律是化学势 μ_i 在各相中相等的连续性条件与独立变量数清点的直接结果,本质上是「变量数的问题」而非第二定律的核心内容。Gibbs(1875 年)在其原始论文《论非均匀物质的平衡》(On the Equilibrium of Heterogeneous Substances)中作为吉布斯能极小化的必要条件提出。本工具仅做简单的四则运算,没有显式使用任何热力学势函数。

最后是把 「相律应用于超临界流体或二级相变」的过度推广限制。临界点处液相与气相区别连续消失,吉布斯相律严格不适用(前提是一级相变)。同样液晶向列-近晶相变、超导转变、合金有序-无序转变等二级相变,均不能简单套用相律,需要 Landau 对称性破缺理论或临界指数讨论。本工具是限于一级相变的教学工具,超临界流体与临界现象的分析请参考其他框架。