实际气体与理想气体的能量特性
实际气体与理想气体的能量特性
在工程热力学中，气体的能量特性是分析热机、制冷循环以及过程设计的基础。理想气体模型因其数学简洁而被广泛采用，但在高压、低温或强相互作用的工况下，实际气体的行为偏离理想模型，需要引入修正。本文系统阐述理想气体与实际气体的内部能、焓、比热以及状态方程的差异，并通过典型气体的数值示例说明如何进行能量计算。
理想气体的能量特性
1. 基本假设
- 分子之间无相互作用（仅碰撞）。
- 分子体积相对于容器体积可忽略不计。
- 运动遵循经典统计力学。
2. 内部能 (U) 与温度的关系
[
U = n C_{V,m} T
]
其中 (C_{V,m}) 为理想气体的摩尔定容比热，仅为常数（单原子气体 (C_{V,m}= \frac{3}{2}R)，双原子气体约为 (\frac{5}{2}R)）。
3. 焓 (H) 与温度的关系
[
H = U + pV = n C_{p,m} T
]
(C_{p,m}=C_{V,m}+R) 同样为常数。
4. 比热的温度依赖
在理想气体模型中，(C_{V,m}) 与 (C_{p,m}) 与温度无关，这为工程计算提供了极大的便利。
实际气体的能量特性
1. 影响因素
- 分子间相互作用（吸引力、排斥力）。
- 分子体积（有限体积导致可压缩性下降）。
- 温度、压力 对相互作用的显著影响。
2. 内部能的表达
实际气体的内部能不再仅是温度的函数，而是 (U = U(T,,p)) 或 (U = U(T,,v))。常用的热力学关系式为
[
\left(\frac{\partial U}{\partial v}\right)_T = T\left(\frac{\partial p}{\partial T}\right)_v - p
]
该式表明在等温条件下，体积变化会导致内部能的变化。
3. 焓的表达
[
H = U + p v = H(T,,p)
]
同样需要考虑压强对分子间势能的贡献。
4. 比热的温度与压强依赖
实际气体的定容比热 (C_V) 与定压比热 (C_p) 均随 (T) 与 (p) 变化。常用的计算方法是：
[
C_p - C_V = T\frac{\alpha^2}{\beta_T} , R
]
其中
(\alpha = \frac{1}{v}\left(\frac{\partial v}{\partial T}\right)_p) 为膨胀系数，
(\beta_T = -\frac{1}{v}\left(\frac{\partial v}{\partial p}\right)_T) 为等温压缩系数。
常用模型与修正
| 模型 | 形式 | 关键参数 | 适用范围 |
|---|---|---|---|
| 范德瓦尔斯方程 | ((p + a/v^2)(v - b) = RT) | (a)（吸引常数），(b)（排斥常数） | 中等压强、低至中温 |
| 雷德利克-库恩方程 | (\displaystyle p = \frac{RT}{v-b} - \frac{a}{v(v+b)}) | 同上，加入 (b) 的二次项 | 高压气体 |
| Virial（Virial）方程 | (\displaystyle \frac{p v}{RT}=1 + \frac{B(T)}{v}+ \frac{C(T)}{v^2}+ \dots) | (B(T), C(T)) 为温度函数的系数 | 低至中等压强，适合实验拟合 |
| 对应状态法 | (\displaystyle Z = \frac{p v}{RT}=f!\left(\frac{p_r}{T_r}\right)) | (Z) 为压缩因子，(p_r, T_r) 为相对压强/温度 | 当有可靠的对应图（如亨利图）时 |
计算内部能的常用积分式
[
U(T,p) = U^{\text{ideal}}(T) + \int_{0}^{p}!\left[ v - T\left(\frac{\partial v}{\partial T}\right)_p \right] ,dp
]
该式利用实际气体的 压缩因子 (Z)：
[
v = Z \frac{RT}{p}
]
把 (Z) 的经验表达式代入，即可得到修正项。
示例计算：氮气在 10 MPa、300 K 下的焓
1. 选取模型
氮气在 10 MPa（约 100 bar）属于高压工况，常采用 范德瓦尔斯方程 或 对应状态法。这里使用范德瓦尔斯参数：
- (a = 1.390 \ \text{Pa·m}^6\text{/mol}^2)
- (b = 3.91 \times 10^{-5}\ \text{m}^3\text{/mol})
2. 计算实际比容
[
(p + \frac{a}{v^2})(v - b) = RT
]
令 (p = 10\ \text{MPa}=10^7\ \text{Pa})，(R = 8.314\ \text{J/(mol·K)})，(T = 300\ \text{K})。通过数值迭代（Newton‑Raphson）得到
[
v \approx 2.68 \times 10^{-3}\ \text{m}^3\text{/mol}
]
对应压缩因子
[
Z = \frac{p v}{RT} \approx \frac{10^7 \times 2.68\times10^{-3}}{8.314 \times 300} \approx 1.08
]
3. 计算内部能修正
先求 (\displaystyle \left(\frac{\partial v}{\partial T}\right)_p)。对范德瓦尔斯方程求偏导可得
[
\left(\frac{\partial v}{\partial T}\right)_p = \frac{R}{p + a/v^2}
]
代入数值：
[
\left(\frac{\partial v}{\partial T}\right)_p \approx \frac{8.314}{10^7 + 1.390/ (2.68\times10^{-3})^2}
\approx 8.2 \times 10^{-7}\ \text{m}^3\text{/mol·K}
]
修正积分（近似为常数乘以压强）：
[
\Delta U = \int_{0}^{p}!\bigl[ v - T(\partial v/\partial T)_p \bigr],dp
\approx \bigl[ v - T(\partial v/\partial T)_p \bigr] p
]
[
\Delta U \approx \bigl[2.68\times10^{-3} - 300 \times 8.2\times10^{-7}\bigr]\times10^7
\approx 2.68\times10^{4} - 2.46\times10^{3}
\approx 2.43\times10^{4}\ \text{J/mol}
]
4. 求焓
理想气体焓（以 300 K 为基准）：
[
H^{\text{ideal}} = n C_{p,m}^{\text{ideal}} T = 1 \times \frac{7}{2}R \times 300 \approx 8.73\times10^{3}\ \text{J/mol}
]
实际焓：
[
H = H^{\text{ideal}} + \Delta U + p v - R T
]
其中 (p v = Z R T = 1.08 \times 8.314 \times 300 \approx 2.69\times10^{3}\ \text{J/mol})。
[
H \approx 8.73\times10^{3} + 2.43\times10^{4} + 2.69\times10^{3} - 2.49\times10^{3}
\approx 3.64\times10^{4}\ \text{J/mol}
]
结果：在 10 MPa、300 K 条件下，氮气的焓约为 36.4 kJ/mol，比理想气体的 8.7 kJ/mol 高出约 27.7 kJ/mol，体现了高压下分子间相互作用对能量的显著影响。
小结
- 理想气体：内部能、焓仅随温度线性变化，比热为常数，计算简便。
- 实际气体：内部能、焓受压强和体积的共同影响，需要通过状态方程（如范德瓦尔斯、Virial）引入压缩因子 (Z) 进行修正。
- 比热：实际气体的 (C_V) 与 (C_p) 随 (T) 与 (p) 变化，可利用热膨胀系数 (\alpha) 与等温压缩系数 (\beta_T) 进行推导。
- 工程应用：在设计高压容器、压缩机或低温制冷系统时，必须采用实际气体模型，否则会低估能量需求或安全裕度。
掌握上述理论与计算方法，能够在工程热力学分析中准确评估实际气体的能量特性，为系统优化提供可靠依据。