查看“范德華方程式”的源代码

{{noteTA
|G1 = 物理學
}}

'''范德华方程'''（'''van der Waals equation'''）（一译'''范德-{}-瓦耳斯方程'''），简称'''范氏方程'''，是[[荷兰]][[物理学家]][[范德华]]于1873年提出的一种[[狀態方程|实际气体状态方程]]{{NoteTag|严格地说，范氏方程也适用于某些状态下的液体，但常用的还是在气体中，参见''适用范围''。}}。范氏方程是对[[理想气体状态方程]]的一种改进，特点在于将被[[理想氣體|理想气体]]模型所忽略的的气体分子自身大小和分子之间的相互作用力考虑进来，以便更好地描述气体的[[宏观]]物理性质。

== 方程的形式 ==
范德华方程具体形式为：
: <math>\left(p + \frac{a'}{v^2}\right)\left(v-b'\right) = kT</math>
式中 
*<math>p</math>为气体的[[压强]]
*<math>a'</math>为度量分子间引力的[[唯象模型|唯象]]参数
*<math>b'</math>为单个分子本身包含的体积
*<math>v</math>为每个分子平均占有的空间大小（即气体的体积除以总分子数量）;
*<math>k</math>为[[玻尔兹曼常数]]
*<math>T</math>为[[热力学温标|绝对温度]]

更常用的形式为：
: <math>\left(p + a \frac{n^2}{V^2}\right)\left(V-nb\right) = nRT</math>
在第二个方程裡
*<math>V</math>为总体积
*<math>n</math>为[[摩尔 (单位)|物质的量]]
*<math>a</math>为度量分子间引力的参数 <math>a=N_A^2 a'</math>
*<math>b</math>为1摩尔分子本身包含的体积之和 <math>b=N_A b'</math>,
*<math>R</math>为[[普适气体常数]]
* <math>N_A</math> 为[[阿伏加德罗常數]]. 

下表列出了部分气体的''a''，''b'' 的值
{| class="prettytable"
|- bgcolor="#efefef"
|气体种类 
| a [kPa (dm³/mol)²]
| b [dm³/mol]
|- 
| 氦气(He) 
| 3.45
| 0.024
|- 
| 氢气(H<sub>2</sub>)
| 24.32
| 0.027
|- 
| 氮气(N<sub>2</sub>)
| 141.86
| 0.039
|- 
| 氧气(O<sub>2</sub>)
| 137.80
| 0.032
|- 
| 二氧化碳(CO<sub>2</sub>)
| 364.77
| 0.043
|- 
| 水蒸气(H<sub>2</sub>O)
| 557.29
| 0.031
|}

在上述方程中必须严格区分总体平均性质和单个分子的性质。譬如，第一个方程中的<math>v</math>是每个分子平均占有空间的大小（可以理解成分子平均“势力范围”的大小），而<math>b'</math>则为单个分子本身“包含”的体积（若为单原子分子如[[稀有气体]]，<math>b'</math>就是[[原子半径]]内包含的体积）。

== 适用范围 ==
范氏方程对气-液临界温度以上流体性质的描写优于理想气体方程。对温度稍低于临界温度的液体和低压气体也有较合理的描述。

但是，当描述对象处于[[状态参量]]空间(P,V,T)中气液[[一级相变|相变]]区（即正在发生气液转变）时，对于固定的温度，气相的压强恒为所在温度下的[[饱和蒸气压]]，即不再随体积<math>V</math>（严格地说应该是单位质量气体占用的体积，即<U>比容</U>）变化而变化，所以这种情况下范氏方程不再适用。

== 方程的提出 ==
[[File:VdW-force-attractive.jpg|170px|right|thumb|水分子之间的范氏引力(凡得瓦力)（中国大陆的中学教科书称为“'''[[范德华力]]'''”）]]

[[File:VdW-V-Excluded.jpg|170px|left|thumb|一个[[双原子分子]]的排斥体积（图中黑色的部分）]]

下面以理想气体状态方程为基础，推导范氏方程。若把气体视为由体积无限小、相互之间无作用力的分子组成，这种模型便是理想气体模型，与其相对应的状态方程是：

:<math> p = \frac{kT}{v}</math>

若抛弃前一个的假设，把组成气体的分子视为有一定大小的<U>刚性球</U>（其半径称为[[范德华半径]]），用<math>b</math>表示这些“球”的体积，上面的方程便改写为：

:<math> p = \frac{kT}{v - b}</math>

在这里，每个分子的“占有体积”<math>v</math>被所谓“排斥体积”<math>v-b</math>代替，反映了分子在空间中不能重叠。若气体被压缩至体积接近分子体积之和（即分子间空隙<math>v-b</math>趋向于0），那么其压强将趋于[[无穷|无穷大]]。

下一步，我们考虑原子对之间的引力。引力的存在会使分子的平均[[亥姆霍兹自由能]]下降，减少量正比于流体的[[密度]]。但压强的大小满足热力学关系

:<math> p = - \left(\frac{\partial A^*}{\partial v}\right)_T </math>

式中''A*'' 为每个分子的亥姆霍兹自由能。由此得到，引力使压强减小的量正比于<math>\frac{1}{v^2}</math>。记该比例常数为<math>a</math>，可得

:<math> p = \frac{kT}{v-b}-\frac{a}{v^2}</math>

这便是范氏方程。

== 与理想气体方程模拟结果的比较 ==

=== 低压状况 ===
在气体[[压强]]不太高的情况下，以下事实成立：
* 排斥体积<math>b</math>的影响相对<math>V</math>而言极小，可以忽略；以[[二氧化碳]]（<chem>CO2</chem>）为例，在标准状况（0°C，1[[气压|标准大气压]]）下，1[[摩尔 (单位)|mol]]<chem>CO2</chem>体积<math>V</math>为 22414cm³，而相应的<math>b=43</math>cm³，比<math>V</math>小3个数量级；
* 分子间的距离足够大，<math>\frac{a}{V^2}</math>项完全可以视为0；譬如在一大气压下二氧化碳气体的<math>\frac{a}{V^2}</math>值只有7‰。
所以此时理想气体方程是范氏方程（也是对实际气体行为的）的一个良好近似。

[[File:VdW-Attr-Repul.jpg|thumb|left|250px|分别用理想气体方程和范德华方程模拟的二氧化碳气体70°C时的p-V等温线]]
=== 中高压状况 ===

随着气体压力的增加，范氏方程和理想气体方程结果的差别会变得十分明显（左图为<chem>CO2</chem>分别用理想气体方程和范德华方程模拟的p-V等温线，温度70&nbsp;°C）：

* 在压强为5000~15000kPa（50~150标准大气压）的<U>中压区</U>，由于体积被“压小”导致分子间距靠近，分子间的引力（表现为<math>\frac{a}{V^2}</math>项）变得不可忽略。<math>\frac{a}{V^2}</math>项的存在使得气体的压强比不考虑分子间引力的理想气体模型估计结果要小（所以左图的中压区里红线比蓝线要低）。
* 在压强为15000kPa以上的<U>高压区</U>，体积的急剧压缩致使<math>b</math>的影响不可忽略，于是范氏方程中的体积项<math>V-Nb</math>（或比容项<math>v-b</math>）将比理想气体方程中的体积项要小（或者说：对应相同体积/比容值的压强项会升高）。这一效应导致在高压区范氏气体的状态线重新赶上并超过理想气体线（见左图的左上角）。

== 用范氏方程描述气体的液化 ==
范氏方程适用于气体的[[液化]]过程。气体液化可能发生的'''最高'''温度称为'''临界温度'''，用<math>T_C</math>表示：
* 当温度<math>T>T_C</math>时，无论给气体施加多大的压强都无法将它液化；
* 当温度<math>T<T_C</math>时，气体可在压强大于一定值时液化，且这一压强随着温度<math>T</math>下降而下降；

[[File:VdW_liquefaction_CO2.jpg|thumb|right|250px|用范德华方程模拟的二氧化碳气体不同温度下压缩过程的p-V等温线，在临界温度以下时能看见明显的液化过程]]

右图所示为用范氏方程模拟的<chem>CO2</chem>在不同温度下的''p-V'' 等温线，从中可以明显看出范氏方程对液化过程的模拟（注意：若用理想气体状态方程作上述模拟，得到的只是一系列[[双曲线]]，因为在等温条件下理想气体状态方程就退化为[[玻意耳-马略特定律]]——''pV''=常数）。<chem>CO2</chem>气体的临界温度为<math>T_C=31</math>°C<math>=304</math>K。

* 70°C 时的曲线（右图中蓝线）形状仍与[[罗伯特·波义耳|玻意耳]]定律的结果（双曲线）类似，尽管位置要略低；
* 当温度下降到40&nbsp;°C，曲线（右图中右二的曲线）形状发生明显的变化，表现为两个拐点的出现。但此时二氧化碳仍然以气态存在；
* 温度进一步降至临界温度31°C（图中红线），若此时气体受压至体积小于某定值<math>V_C</math>（随温度变化而变化），则气体将发生液化。图中<math>V>V_C</math>时曲线对应气态<chem>CO2</chem>的<math>p</math>，<math>V</math>值，<math>V<V_C</math>时曲线对应液态<chem>CO2</chem>的''<math>p</math>''，<math>V</math>值；
* 图中13°C 和 21°C 对应的曲线只有两拐点以外的部分是与物理实际相符的。当气体被进一步压缩至比右拐点对应体积更小时，气体将进入液化区，在液化过程中实际气体的p-V线应是一段“平台”，而不是如图所示的“驼峰”型。但完全液化后，液态<chem>CO2</chem>的压强却仍能被图中曲线恰当地反映，此时曲线随体积的减小而剧烈上升，这一定程度上反映了液体的[[不可压缩性]]。另外，我们从图中能得到的另一个信息就是“液化平台”的长度随温度的下降而增加；

气体的临界状态参量<math>V_C</math>、<math>p_C</math>、<math>T_C</math>和范德华常数<math>a</math>、<math>b</math>之间存在下列数学关系：

:<math> {{V}_{C}} ={3Nb}</math>

:<math display="inline">{p}_{C} = {a}/{27}{b}^{2}</math>

:<math>{T}_{C} = {8a}/{27bR}</math>

我们可以利用这些关系通过测出气体的''<math>T_C</math>''和对应的''<math>p_C</math>''来得到<math>a</math>和<math>b</math>的值（由于测量上的困难，一般不使用''<math>V_C</math>''）。

== 其他热力学参量 ==
下面，我们不再考虑<math>v=\frac{V}{N}</math>（<math>N</math>为系统中的分子数），改为考虑总体体积<math>V</math>。

状态方程并不能告诉我们系统的所有热力学参量。我们可以照搬上面推导范氏方程的思路，从理想气体的亥姆霍兹自由能表达式出发，推得下面的结论：

:<math>A(T,V,N)=-NkT\left(1+\ln\left(\frac{(V-Nb)T^{\hat{c}_V}}{N\Phi}\right)\right)
-\frac{aN^2}{V}</math>

式中<math>A</math>为亥姆霍兹自由能，<math>\hat{c}_v</math>是无量纲的[[定容热容]]，<math>\Phi</math>是待定的[[熵]]常数。上述方程将<math>A</math>用它的[[自然变量]]<math>V</math>和<math>T</math>表示{{NoteTag|1 = 亥姆霍兹自由能''A'' 是一个[[特性函数]]，它作为温度''T'' 和体积''V'' 的[[全微分]]表达式为 d''A'' = -''S''d''T''-''p''d''V''，参见[[热力学势函数]]。}}，所以系统的所有热力学信息已全部知道。其力学状态方程就是前面导出的范氏方程

: <math>P = -\left(\frac{\partial A}{\partial V}\right)_T 
= \frac{NkT}{V-Nb}-\frac{aN^2}{V^2}</math>

系统的[[熵]](<math>S</math>)由下式决定

: <math>S = -\left(\frac{\partial A}{\partial T}\right)_V
=Nk\left[ \ln\left(\frac{(V-Nb)T^{\hat{c}_V}}{N\Phi}\right)+\hat{c}_V+1 \right]</math>

综合<math>A</math>和<math>S</math>的表达式，可由定义得到系统[[内能]]

: <math>U = A+TS = \hat{c}_V\,NkT-\frac{aN^2}{V}</math>

其他[[热力学函数|热力学势]]和[[化学势]]也可用类似的方程给出，但任何势函数若要用压强<math>P</math>表示都需要求解一个三阶多项式，使结果的形式变得很繁杂。所以，将[[焓]]和[[吉布斯能]]用它们相应的自然变量表示的结果都是复杂的（因为<math>P</math>是它们的自然变量之一）。

== 简化形式 ==
虽然在一般形式的范氏方程中，常数<math>a</math>和<math>b</math>因气体/流体种类而异，但我们可以通过改变方程的形式，得到一种适用于<U>所有</U>气体/流体的普适形式。 

按照下面的方式定义'''约减变量'''（亦称'''折合变量'''，就是把变量转换成其[[无量纲化|无量纲形式]]），其中下标<math>R</math>表示约减变量，下标<math>C</math>表示原变量的[[临界点 (热力学)|临界值]]： 
: <math>p_R=\frac{p}{p_C}</math>,
: <math>v_R=\frac{v}{v_C}</math>,
: <math>T_R=\frac{T}{T_C}</math>,

式中<math>p_C=\frac{a}{27b^2}</math>，<math>\displaystyle{v_C=3b}</math>，<math>kT_C=\frac{8a}{27b}</math>。

用约减变量代替原变量，范氏方程形式变为 

: <math>\left(p_R + \frac{3}{v_R^2}\right)(v_R - 1/3) = \frac{8}{3} T_R</math>

这就是范氏方程的<U>不变形式</U>，即这一形式不会因应用流体种类改变而改变。 

上述方程的不变性质亦称[[对应状态原理]]。

== 在可压缩流动中的应用 ==
在[[流体力学]]中，范氏方程可以作为可压缩流体（如液态[[高分子材料]]）的[[PVT状态方程]]。这种情况下，由于[[比容]]<math>V</math>变化不大，可将方程简化为：

<math>(p+A)(V-B)=CT\,</math>,

其中<math>p</math>为压强，<math>V</math>为比容，<math>T</math>为温度，<math>A</math>、<math>B</math>、<math>C</math>均为与对象相关的参数。

== 注释 ==
{{NoteFoot}}

== 參考文獻 ==
{{Reflist}}
* 赵凯华、罗蔚茵 著：《新概念物理教程·热学》（第二版），高等教育出版社，北京，2005，ISBN 7-04-017680-7
* 汪志诚 编：《热力学·统计物理》（第二版），高等教育出版社，北京，1993，ISBN 7-04-004360-2

== 外部链接 ==
{{Commons category|Van der Waals equation|范德华方程}}
* [https://web.archive.org/web/20060504071400/http://www.ac.wwu.edu/%7Evawter/PhysicsNet/Topics/Thermal/vdWaalEquatOfState.html 不同气体的范氏常数''a'' 和''b'' ]

== 参见 ==
{{Portal box|物理}}
* [[状态方程]]
* [[气体定律]]
* [[范德华常数表]]
* [[范德华力]]
* [[范德华半径]]

{{-}}
{{Statistical mechanics topics}}

[[Category:氣體定律|V]]
[[Category:流体力学]]
[[Category:物理化学]]
[[Category:方程]]
[[Category:热力学]]