跳转至

气体

气体通常是热学里最先详细研究的物态,因为它最容易建立简洁模型。气体分子之间平均距离较大、相互作用相对容易近似处理,所以很多温度、压强和状态方程的直观关系,会先在气体这里变得清晰。

学习目标

读完本页后,你应该能够:

为什么要先单独讲气体

本页将重点处理什么

正文扩写时,这一页会重点讲:

核心内容

理想气体的微观假设

  1. 质点(分子线度<<分子间平均距离)
  2. 遵从牛顿力学规律(经典粒子)
  3. 除碰撞瞬间,分子间、分子与器壁间无相互作用力
  4. 分子间、分子与器壁间的碰撞为完全弹性碰撞,即平动动能不向内部自由度(如转动、振动、激发等)转移

关于分子集体运动的统计假设

  1. 通过碰撞分子速度不断变化。
  2. 无外场时,平衡态分子按位置均匀分布。即 n=dNdV=NV
  3. 平衡态分子速度按方向的分布是均匀的(各向同性)。
vx=vy=vz=0,vx2=vy2=vz2=13v2,vα2=1Ni=1Nvαi2.
注意
  1. dV 应是宏观小微观大的体积元。例如:标准状态下,1cm3 的空气中有 2.7×1019 个分子。 若 dV=1015m3,则 dV 内有 1010 个分子!
  2. 涨落:各时刻的 dNdV 值相对于平均值 n 的差别。通常 dV 取足够大,使这一涨落比起平均值 n 可以小到忽略不计。

分子弹性碰撞

下面我们从两个分子碰撞的行为加入合理的统计假设得到气体分子热运动的微观统计行为与宏观热力学量T的关系。
现在我们仔细地考查一下两个分子碰撞中能量的交换。这两个分子可以都是气体中的,也可以一个是气体中的,另一个是器壁中的。设它们具有不同的质量 mM。令它们在碰撞前的速度分别为 v0V0,碰撞后的速度分别为 vV。把分子看成球形,取两球中心联线为直角坐标系的 x 轴,由于碰撞中垂直于 x 轴的速度分量不改变,问题就简化成只有 x 方向一维了。在此方向上分子 M 损失的动能为

(1)Δε=12M(V0x2Vx2)=12M(V0xVx)(V0x+Vx).

分子弹性碰撞

由于弹性碰撞中动能守恒,这也是分子 m 获得的动能。根据我们在动量与能量 - 弹性碰撞中得到的公式。把式中的符号适当地换过来,有

Vx=(Mm)V0x+2mv0xM+m,

从而

{V0xVx=2m(V0xv0x)M+m,V0x+Vx=2(MV0x+mv0x)M+m.

代入(1)式,得

(2)Δε=2mM(M+m)2[MV0x2+mv0x2+(mM)v0xV0x].

以上是单次碰撞的结果,Δε 可正可负。我们假定

  1. Mm 各代表一类分子,宏观效果是两类分子之间大量碰撞的平均。
  2. 气体和器壁宏观上都是静止的,从而每类分子沿任何正负两个方向的概率都相等。

V0xv0x 的平均 V0x=v0x=0。再者,因 Mm 两类分子的热运动都是随机的,且在方向上彼此没有关联,故 v0xV0x 也等于 0。取(2)式的平均,得

(3)Δε=2mM(M+m)2(MV0x2mv0x2).

上式右端正比于两分子碰撞前后 x 方向平均动能之差。由于分子热运动是各向同性的,与 x 垂直的另外两个方向平均动能与 x 方向一样,故 x 方向的平均动能是各方向平均动能的 ⅓。故(3)式又可写为

(4)Δε=2mM3(M+m)2(MV2mv2)12MV212mv2.

M 类分子的平均动能大于 m 类分子,则 Δε>0,动能由 M 类分子传给 m 类分子;反之,若 M 类分子的平均动能小于 m 类分子,则 Δε<0,动能由 m 类分子传给 M 类分子;Δε=0 表示平均起来两类分子之间能量交换为0。

上节 - 热量及其本质所述,热量即热运动的能量,故 Δε 的正负决定了热量传递的方向,Δε=0Mm 两类分子之间达到热平衡。如此看来,按照热力学第零定律,分子的平均动能在宏观上具有温度的特征。我们暂时假定,分子的平均动能正比于热力学温度:

(1.26)ε=12mv2=KT,
  • 其中 ε 是分子的平均动能,K 是比例常量,T 是热力学温度。

因为热平衡状态温度相等的概念适用于任何物体,上式里的比例常量 K 应该是普适的。下面我们将进一步论证这个结论,并把这个普适常量定下来。

注意
  1. 温度T是描述热力学系统平衡态的宏观状态参量。
  2. 温度T是一个统计概念,只能用来描述大量分子的集体状态,单个分子谈论它的温度毫无意义。
  3. 温度所反映的运动是分子的无规则运动(分子热运动 εt ),和物体的整体运动无关。是相对于系统的质心参考系测量的。系统内所有分子的平动动能的总和就是系统的内动能。
  4. 分子热运动的平均转动动能和平均振动动能也都和温度有直接的联系。

理想气体压强公式

上面我们讨论了温度的微观意义,现在来看压强的微观意义。

什么是压强?

物体中通过内部假想截面ΔS相互作用的法向应力。

  • 应力:作用在单位面积上的力。

气体压强由哪部分贡献?
1. 分子力:ΔS 两侧附近的分子相互作用着; 2. 分子运动:大量分子携带着动量穿过 ΔS 的平均效果是均 匀而持续的压力。 - 对理想气体,只有 2 部分的贡献。

大量分子产生持续的平均冲力曲线

大量分子产生持续的平均冲力曲线

推导前置定义

符号定义与公式
m单个气体分子质量
N容器内气体分子总数量
V容器容积(气体总体积)
n总分子数密度:单位体积内的所有气体分子总数,公式为 n=NV
ni速度为vi的分子数密度:单位体积内速度大小方向等于vi的分子数量,公式为 ni=NiV
其中 Ni 是速度等于vi的分子总数

总分子数:

N=iNi

总分子数密度:

n=ini

示意图

冲量与压强的微观推导

一个分子对器壁上 ΔS 的冲量 —— 2pix

分子动量改变: pixpix

Δt 时间内所有动量为 pi 的分子对 ΔS 的冲量(考虑斜柱体积):

ΔIi=(nivixΔtΔS)(2pix)=2nipixvixΔtΔS

Δt 时间内所有分子对 ΔS 的冲量:

注意:vix>0vix<0 分子各占一半。

ΔI=vix>0ΔIi=12iΔIi=ΔtΔSinipixvix

压强:

=ΔFΔS=ΔIΔtΔS=inipixvix

求统计平均值

=inipixvix=nininpixvix=npxvx

其中:

pxvx=1ninipixvix=inipixvixini

由于动量和速度分布的各向同性,则:

pxvx=pyvy=pzvzpv=pxvx+pyvy+pzvz=3pxvx

故总压强公式可以写成如下形式:

=13npv(理想气体的压强公式)

不同情形下的压强公式

非相对论情形
p=mv,pv=mv2=2ε=13npv=13nmv2=23nεt

说明: 三个统计平均量 p,nεt 之间相互联系的一个统计规律,而不是一个力学规律。

极端相对论情形
p=m0v1v2/c2

vc 时:

p=ε2m02c4cε/cpvpcε

代入压强公式:

=13npv=13nε

在讨论光子气体时大有用处!

提示

压强不是只在器壁处才有,气体内部也是有压强的。器壁上气体的压强也与气体内部一样。

下面我们回到非相对论情形,讨论非相对论粒子组成的理想气体。

理想气体定律的推导

下面我们利用理想气体压强公式推导几条定律。

1. 理想气体物态方程

从理想气体压强公式出发:

p=23nε

其中: - n=N/V —— 分子数密度 - ε=αT —— 分子平均平动动能(与温度成正比)

代入得:

pV=23NαT

引入物质的量 ν 和阿伏伽德罗常数 NA

N=νNA

则:

pV=ν23NAαT
与实验定律对比

对比之前根据气体三大定律推导出来的物态方程 pV=νRT,可知:

23NAα=R

解得:

α=32RNA=32k

其中:

k=RNA=1.380649×1023 J/K

k 称为玻尔兹曼常数

分子的平均平动动能

因此:

εt=32kT

压强公式可写为:

p=nkT

重要结论:温度是分子平均平动动能的量度,温度越高,分子平均平动动能越大。

2. 方均根速率

方均根速率可以用来表征分子平均平动动能情况。

由:

εt=12mv2=32kT

可得方均根速率:

vrms=v2=3kTm=3RTM

其中: - m —— 单个分子质量 - M —— 摩尔质量(M=mNA

重要性质:同一温度下,质量大的分子其方均根速率小。

例题

求零摄氏度时,氢气和氧气分子的平均平动动能和方均根速率。 解: 温度 T=0C=273.15 K

平均平动动能

由于平均平动动能只与温度有关:

εt=32kT=32×1.38×1023×273.155.65×1021 J

氢气和氧气分子的平均平动动能相同

方均根速率

氢气MH2=2×103 kg/mol):

vrms, H2=3RTMH2=3×8.31×273.152×103=1.84×103 m/s

氧气MO2=32×103 kg/mol):

vrms, O2=3RTMO2=3×8.31×273.1532×103=461 m/s

发现:结果与声波在空气中传播的速率同数量级。

结论
  1. 同一温度下,不同气体的分子平均平动动能相等,为 32kT
  2. 分子质量越大,方均根速率越小
  3. 氢气分子的方均根速率约为氧气的 32/2=4 倍。

3. 玻尔兹曼常数 k 与普适气体常数 R 的物理意义

根据分子平均平动动能公式:

εt=32kT

玻尔兹曼常数 k 可以看成是温度与能量之间的单位换算比率。它建立了微观粒子热运动能量(单位:焦耳 J)与宏观温度(单位:开尔文 K)之间的桥梁。

而根据摩尔内能公式:

Etmol=NAεt=32RT

普适气体常数 R 则是温度与每摩尔物质所含能量之间的单位换算比率。它连接的是宏观尺度下的“摩尔量”与温度。

关键关系
R=NAk 即:1 摩尔粒子的总“热容系数”等于单个粒子的热容系数乘以阿伏伽德罗常数。

4. 各种能量单位之间的换算关系

在统计物理和热力学中,不同领域习惯使用不同的能量单位。下表提供了常用能量单位间的换算因子(基于 T=1K 时的能量值或等效能量)。

各种能量单位之间的换算关系
eVKaJ*kJ/molkcal/mol
eV11.16×1041.60×1019.6×1012.3×101
K0.86×10411.38×1058.3×1031.98×103
aJ6.257.25×10416.0×1021.43×102
kJ/mol1.04×1021.20×1021.66×10312.4×101
kcal/mol4.4×1025.0×1026.9×1034.21

注:1aJ=1018J (attojoule,阿焦耳)

使用说明
  • 表格中的数值表示:1 单位行对应的能量 = 列标题单位的多少倍
  • 例如:第一行第二列 “1.16×104” 表示:
    $$ 1\,\text{eV} = 1.16 \times 10^4\,\text{K} \cdot k $$
    即:1 电子伏特对应的“温度当量”约为 11600 K(通过 E=kT 转换)。
  • 同理,第三行第一列 “6.25” 表示:
    $$ 1\,\text{aJ} = 6.25\,\text{eV} $$
    因为 1eV=1.602×1019J=0.1602aJ,所以倒数即为 6.25。
  • 这些换算是基于以下基本常数:
  • k=1.380649×1023J/K
  • NA=6.02214076×1023mol1
  • 1eV=1.602176634×1019J
  • 1kcal=4184J

阿伏伽德罗定律

  • 原本是一条经验的定律;
  • 从实验事实得知标准状态下的摩尔体积 ω0=22.41410 L/mol 是个与物质无关的普适常数。根据 R=p0ω0/T0,可知气体常量 R 也是普适的。
  • 根据这章节的微观理论,从微观上可以解释阿伏伽德罗定律。
ε=αT热平衡状态温度相等,适用于任何物体,α 是普适的。α=32RNA 普适R 普适ω0 普适

道尔顿分压定律

混合气体的压强等于各组分的分压强之和,即

p=pα+pβ+

由于混合气体中各组分处于热平衡,它们的温度相同,因而各组分的分压强

pα=23nαεα=nαkT,pβ=23nβεβ=nβkT,

因而总压强

p=(nα+nβ+)kT=nkT
例题

题目:已知空气中几种主要组分的分压百分比是氮(N2) 78%,氧(O2) 21%,氩(Ar) 1%,求它们的质量百分比和空气在标准状态下的密度。

:因温度相同,

pαp=nαkTnkT=nαn=vαNA/VvNA/V=vαv

各组分的质量比

vN2MN2mol:vO2MO2mol:vArMArmol=0.78×28:0.21×32:0.01×40 mN2:mO2:mAr=75.4%:23.2%:1.4%

再归一,可得

ρ=vN2MN2mol+vO2MO2mol+vArMArmolv22.4 L/mol=0.78×28g+0.21×32g+0.01×40g22.4×103 cm3

标准状态下空气的密度为

=1.29×103 g/cm3=1.29 kg/m3

实际气体

在推导理想气体物态方程时,我们几乎把分子力完全忽略了,但在实际气体中它还是有影响的。不过在气态中分子力的效应毕竟比较小,我们把它当作对理想气体模型的修正来处理。

分子力的效应表现在体积压强两个方面。先看它对体积的影响。如 热量及其本质 - 分子力 所述,分子力由短程的排斥力和较为长程的吸引力组成,对体积的影响主要在于前者。

  • 体积

形象地说,排斥力的作用相当于分子有一个有一定半径的硬壳,在硬壳占有的体积内,其他分子是不能侵入的,所以应从理想气体物态方程中的体积V里把这部分体积扣除。所有分子占有的体积正比于气体中分子总数,亦即正比于气体的摩尔数ν。于是我们把式中的V换成V=Vνb,这里b是1 mol气体的所有分子所占有的体积,V代表每个分子真正能够在其中自由活动的有效体积。

“硬壳”

  • 压强

压强来源:
- 分子运动 —— 动理压强 pk
- 分子力(相互作用势,影响主要在于较为长程的吸引力部分)—— 内压强 pU

对于处在气体内部的一个分子,周围分子给它的吸引力平均说来抵消了,从而对它的自由飞行不产生影响,只要把体积的修正考虑进去,下式对描写pk仍旧是有效的,即

pk=23nε=νRTV=νRTVνb.

现在来看内压强 pU,它代表气体的内聚力,正比于施力者和受力者的粒子数密度。

pUn2ν2V2pU=ν2aV2

其中 a 是反映分子之间吸引力的常量(吸引势)。

则气体对器壁的实际压强 p 为:

p=pk+pU=νRTVνbν2aV2

这里的p是气体内部某一点的压强,作用在器壁上的压强怎么样?器壁对气体分子的吸引力会对压强产生什么影响?你不妨自己思考看看。

整理后得到范德瓦尔斯物态方程(1910年诺贝尔物理奖):

(p+ν2aV2)(Vνb)=νRT

问题:那么, ab 怎么定出来?
哈,那当然是由实验来确定啦!

范德瓦耳斯修正量 a,b 的实验值

气体a/[atm(Lmol1)2]b/(Lmol1)
He0.03410.0234
H20.2470.0265
O21.3690.0315
N21.3610.0385
CO23.6430.0427
H2O5.5070.0304

注意:使用的时候注意转换成国际单位制。

对应态原理

三个临界参量之间还有一个简单的关系,K 叫做临界系数

K=νRTKpKVK=83=2.667

这个无量纲的比值 K 叫做临界系数。利用临界参量可以把范德瓦尔斯方程无量纲化

定义对比参数:

τ=T/TK,ω=V/VK,π=p/pK

代入原方程可得对应态方程

(π+3ω2)(3ω1)=8τπ=8τ3ω13ω2

此式是一个用对比参数表示的状态方程式。特点是式中不包含反映个别物质特性的常数。

范德瓦尔斯提出的两参数对应态原理:在相同的对比状态下,所有的物质表现出相同的性质。即对于不同的气体,当具有相同的对比温度和对比压力时,则具有相同的对比体积(或压缩因子)。

系统地处理非理想气体物态方程的方法是位力展开(virial expansion),即将压强与温度的关系按密度的幂次展开:

(1.36)p=νRTV[1+νVB(T)+(νV)2C(T)+],

式中的B(T)和C(T)分别称为第二和第三位力系数,它们都只是温度的函数。
位力系数的表达式可由统计物理理论导出,其数值也可由实验来测定。将范德瓦耳斯方程按密度(ν/V)的幂次展开,

p=νRTV(1νbV)1ν2aV2=νRTV[1+νbV+(νbV)2+]ν2aV2,

与范德瓦尔斯方程对比即可得范德瓦耳斯气体的位力系数:

B(T)=baRT,C(T)=b2.

B(T)在低温下是负的,高温时变正,在转折点

TB=aRb

处为0。这时物态方程与理想气体的玻意耳定律偏离最小,转折温度TB称为玻意耳温度。



评论