第二章:热力学基础回顾——吉布斯自由能、焓、熵、化学势

各位同学,咱们今天聊点“硬核”的。做CALPHAD,说白了就是在跟热力学打交道。你如果不懂吉布斯自由能,那后面的相图计算基本就是瞎蒙。我当年刚入行时,也觉得这些概念太抽象,直到有一次算一个三元体系的相平衡,死活算不出来,最后发现是化学势的符号搞反了……嗯,从那以后,我再也不敢轻视这些“基础”了。

2.1 吉布斯自由能:相平衡的“总指挥”

先问大家一个问题:为什么有的合金在高温下是液体,冷却后就变成了固体?

答案很简单——系统在寻找能量最低的状态。这个“能量”在热力学里,就是吉布斯自由能(G)

它的定义式是:

G = H - T·S

其中:

  • H 是焓(可以理解为系统的总热量)
  • T 是绝对温度
  • S 是熵(系统的混乱程度)

我个人习惯把吉布斯自由能想象成一个“裁判”。它决定了在某个温度、压力下,哪个相最稳定。哪个相的G值最低,哪个相就“赢”了。

核心判据:在恒温恒压下,系统总是朝着吉布斯自由能减小的方向进行,直到达到最小值——这就是相平衡的判据。

我在做Fe-Cr-Ni体系的计算时,就遇到过这种情况。同一个成分,在800℃时是奥氏体稳定,到了600℃就变成了铁素体。为什么?因为温度变了,G-T·S这一项的影响权重变了。说白了,就是熵在高温下“说了算”。

2.2 焓与熵:一对“欢喜冤家”

咱们把G拆开来看。

焓(H),代表系统的内能加上压力体积功。在固态相变中,它主要反映原子间的键合强度。你想想看,一个相如果原子排列紧密、键能大,它的焓值就低,就更稳定。

熵(S),代表混乱度。温度越高,原子振动越剧烈,混乱度越大,熵的贡献就越重要。

这两者之间有个“博弈”:

  • 低温时,焓占主导——系统倾向于形成有序、低能量的结构
  • 高温时,熵占主导——系统倾向于形成无序、高混乱度的结构

举个例子:纯铁在912℃发生α-Fe(体心立方)到γ-Fe(面心立方)的转变。α-Fe的焓更低,但γ-Fe的熵更高。到了912℃以上,T·S这一项超过了焓的差异,γ-Fe就变成了稳定相。

避坑指南:我曾经在计算一个高温合金的相变温度时,忽略了熵的温度依赖性,结果算出来的转变温度差了50℃。后来我改用更精确的熵模型,才把结果拉回来。所以,千万别把熵当成常数!

2.3 化学势:相平衡的“微观推手”

化学势(μ)这个概念,很多初学者觉得难。其实说白了,它就是某个组元在某个相中的“逃跑倾向”

定义式:

μ_i = (∂G / ∂n_i)_{T,P,n_j≠i}

翻译成人话:在恒温恒压下,往系统里加入一丁点i组元,吉布斯自由能的变化量,就是i组元的化学势。

为什么化学势重要?因为相平衡的条件就是:每个组元在所有相中的化学势相等

举个例子:你在一个两相区(比如液相+固相)里,如果A组元在液相中的化学势比固相中高,那A就会从液相往固相跑,直到两边化学势相等——这就是扩散的驱动力。

记住这个公式:

μ_i^α = μ_i^β = μ_i^γ = ...

这就是CALPHAD计算中求解相平衡的核心方程。

我在做焊接熔池凝固模拟时,就靠这个原理判断元素在固液界面的分配行为。如果某个元素的化学势在液相中更高,它就会偏析到液相中,导致成分不均匀。

2.4 吉布斯自由能最小化原理:CALPHAD的“灵魂”

好了,咱们把前面三个概念串起来。

对于一个多组元、多相系统,CALPHAD要解决的核心问题是:在给定的温度、压力和总成分下,系统应该由哪些相组成?每个相的含量和成分是多少?

答案就是——寻找吉布斯自由能最小的状态

具体做法是:

  1. 建立每个相的吉布斯自由能模型(比如用Redlich-Kister多项式描述溶液相)
  2. 给定总成分和温度
  3. 用数学优化方法(比如牛顿-拉夫森法)求解使总G最小的相分数和相成分

下面我用一个简单的二元体系示意图来说明这个逻辑:

吉布斯自由能最小化原理示意图 成分 (X_B) 吉布斯自由能 G 液相 (L) 固相 (S) X_B^S X_B^L 两相平衡时,公切线连接两个自由能曲线,切点对应平衡成分 液相自由能曲线 固相自由能曲线 公切线(平衡条件)

这张图里,红色虚线是液相的自由能曲线,蓝色实线是固相的。在给定温度下,系统会找到一条公切线,切点对应的成分就是两相平衡时的成分。这就是CALPHAD计算的基本逻辑。

注意:实际CALPHAD计算中,自由能模型往往是非线性的,可能存在多个局部极小值。我曾经在计算一个含金属间化合物的体系时,优化算法掉进了局部极小值,算出来的相图完全不对。后来我改用全局优化算法(比如模拟退火)才找到真正的全局最小值。

2.5 小结:这些概念怎么用?

咱们做个快速回顾:

概念 物理意义 在CALPHAD中的作用
吉布斯自由能 G 系统稳定性的总度量 相平衡判据:G最小化
焓 H 原子间键合强度 决定低温下的相稳定性
熵 S 系统混乱度 决定高温下的相稳定性
化学势 μ 组元的“逃跑倾向” 相平衡条件:各相化学势相等

我个人习惯,每次开始一个新体系的计算前,都会先手动画一下这几个概念的草图。别小看这个习惯,它能帮你避免很多低级错误。比如有一次,我算一个Al-Cu体系的共晶反应,发现液相线形状不对,回头一查,原来是化学势的参考态选错了。

好了,这一章的内容就到这儿。记住:热力学基础是CALPHAD的根,根扎得深,树才能长得高。


专注资料整理