实际气体状态方程（范德瓦尔斯方程）

在工程热力学研究中，理想气体状态方程（$PV=nRT$）是一个极其重要的简化模型。它假设气体分子本身不占体积，且分子之间不存在相互作用力。然而，在实际的工程应用中，特别是在高压或低温环境下，气体分子表现出明显的非理想特性，此时理想气体方程的误差会迅速增大。

为了更准确地描述真实气体的行为，物理学家约翰·范德瓦尔斯（Johannes Diderik van der Waals）提出了修正后的状态方程，即范德瓦尔斯方程（Van der Waals Equation）。该方程通过引入两个关键的修正项，有效地描述了实际气体在偏离理想状态时的物理特性。
在理解范德瓦尔斯方程之前，必须明确理想气体模型在以下两个方面的假设缺陷：

1. 分子体积忽略不计：理想气体假设分子是“质点”，其占据的空间体积相对于容器总体积可以忽略。但在高压状态下，分子间的距离缩短，分子自身的体积（排斥体积）对气体占据的总空间产生了显著影响。
2. 分子间无相互作用力：理想气体假设分子间既没有吸引力也没有排斥力。但在实际情况中，分子间存在范德瓦尔斯力（主要是引力）。这种引力会减缓分子撞击容器壁的速度和频率，从而导致实际压力低于理想压力。

范德瓦尔斯方程的推导逻辑

范德瓦尔斯方程的核心思想是对理想气体方程的两个物理参数进行修正：压力（$P$）和体积（$V$）。

1. 体积修正（排斥体积）

由于实际气体分子本身具有一定的体积，它们在容器内运动时，实际上无法进入其他分子所占据的空间。因此，气体分子能够自由运动的“有效体积”应当小于容器的总体积 $V$。

对于 $n$ 摩尔气体，如果每个分子占据的排斥体积为 $b'$，那么修正后的有效体积 $V_{eff}$ 可以表示为：
$$V_{eff} = V - nb$$
其中，$b$ 是与气体种类相关的常数，代表单位摩尔气体所占有的排斥体积。

2. 压力修正（分子间引力）

在理想气体中，压力完全由分子撞击容器壁产生的冲量决定。但在实际气体中，由于分子间存在引力，靠近容器壁的分子会被后方的分子“向后拉”，这削弱了分子撞击壁面的力度。

这种引力导致的压力损失与分子密度的平方成正比（因为引力作用取决于分子间的距离，而密度平方反映了分子碰撞的频率和相互作用的强度）。因此，修正后的压力 $P_{real}$ 与理想压力 $P_{ideal}$ 的关系可以表示为：
$$P_{ideal} = P_{real} + \frac{an^2}{V^2}$$
或者写成实际压力：
$$P_{real} = P_{ideal} - \frac{an^2}{V^2}$$
其中，$a$ 是表征分子间引力强度的常数。

范德瓦尔斯方程的标准形式

综合上述两个修正项，我们将它们代入理想气体状态方程 $P_{ideal}V_{eff} = nRT$，得到范德瓦尔斯方程的标准形式：

$$\left( P + \frac{an^2}{V^2} \right) (V - nb) = nRT$$

若使用摩尔体积 $v = V/n$ 进行表示，方程可以简化为：
$$\left( P + \frac{a}{v^2} \right) (v - b) = RT$$

参数物理意义总结

• $a$ (吸引力系数)：反映了分子间相互作用力的强弱。$a$ 值越大，说明分子间的引力越强，气体越容易液化。
• $b$ (排斥体积系数)：反映了分子自身占据的空间大小。$b$ 值越大，说明分子体积越大。

实例计算

为了直观理解该方程的应用，我们来看一个关于二氧化碳（$\text{CO}_2$）的计算示例。

已知条件：
在温度 $T = 300\text{ K}$ 下，1 摩尔 $\text{CO}_2$ 气体的体积为 $V = 0.5\text{ L}$。
已知 $\text{CO}_2$ 的范德瓦尔斯常数为：

• $a = 3.59\text{ L}^2\cdot\text{atm/mol}^2$
• $b = 0.0427\text{ L/mol}$
• 气体常数 $R = 0.08206\text{ L}\cdot\text{atm/(mol}\cdot\text{K)}$

问题：

1. 使用理想气体方程计算压力 $P_{ideal}$。
2. 使用范德瓦尔斯方程计算实际压力 $P_{real}$。

计算过程：

1. 理想气体压力计算：
   $$P_{ideal} = \frac{nRT}{V} = \frac{1 \times 0.08206 \times 300}{0.5} = 49.236\text{ atm}$$

2. 范德瓦尔斯实际压力计算：
   首先计算修正项：

   • 压力修正项 $\frac{a}{v^2} = \frac{3.59}{0.5^2} = \frac{3.59}{0.25} = 14.36\text{ atm}$
   • 体积修正项 $v - b = 0.5 - 0.0427 = 0.4573\text{ L}$

   代入方程：
   $$\left( P_{real} + 14.36 \right) \times 0.4573 = 1 \times 0.08206 \times 300$$
   $$P_{real} + 14.36 = \frac{24.618}{0.4573} \approx 53.833$$
   $$P_{real} = 53.833 - 14.36 = 39.473\text{ atm}$$

结果分析：
通过对比可以发现，$P_{real} (39.47\text{ atm}) < P_{ideal} (49.24\text{ atm})$。这说明在当前条件下，分子间的引力起到了主导作用，显著降低了气体对容器壁的压力。

范德瓦尔斯方程的局限性与工程应用

尽管范德瓦尔斯方程比理想气体方程精确得多，但在极高压或接近临界点的复杂工况下，它仍然存在一定的误差。

• 局限性：
  • 它对临界参数（临界压力 $P_c$、临界温度 $T_c$）的预测精度有限。
  • 对于极性分子或复杂的大分子，其简单的二项修正不足以描述复杂的分子间相互作用。
• 工程应用：
  • 在热力学性质计算（如焓、熵的偏差计算）中，它是基础模型。
  • 它是后续更高级状态方程（如 Redlich-Kwong 方程、Peng-Robinson 方程）的理论基石。

在现代化工设计和石油工程中，工程师通常会根据具体的气体种类选择更精确的关联式，但理解范德瓦尔斯方程所蕴含的物理本质——体积排斥与引力修正，是掌握整个工程热力学理论的关键。