范德华方程

**凡得瓦方程式**（**van der Waals equation**）（一譯**范德瓦耳斯方程式**），簡稱**范氏方程式**，是荷蘭物理學家范德華於 1873 年提出的一種實際氣體狀態方程式。范氏方程式是對理想氣體狀態方程式的一種改進，特點在於將被理想氣體模型所忽略的的氣體分子自身大小和分子之間的交互作用力考慮進來，以便更好地描述氣體的宏觀物理性質。

## 方程式的形式

凡得瓦方程式具體形式為：

　
: $\left(p+{\frac {a'}{v^{2}}}\right)\left(v-b'\right)=kT$

式中

 * $p$ 為氣體的壓力
* $a'$ 為度量分子間引力的唯象參數
* $b'$ 為單個分子本身包含的體積
* $v$ 為每個分子平均占有的空間大小（即氣體的體積除以總分子數量）;
* $k$ 為波茲曼常數
* $T$ 為絕對溫度

更常用的形式為：

　
: $\left(p+a{\frac {n^{2}}{V^{2}}}\right)\left(V-nb\right)=nRT$

在第二個方程式裡

 * $V$ 為總體積
* $n$ 為物質的量
* $a$ 為度量分子間引力的參數 $a=N_{A}^{2}a'$
* $b$ 為 1 莫耳分子本身包含的體積之和 $b=N_{A}b'$ ,
* $R$ 為普適氣體常數
* $N_{A}$ 為阿伏加德羅常數.

下表列出了部分氣體的 _a_，_b_ 的值

在上述方程式中必須嚴格區分總體平均性質和單個分子的性質。譬如，第一個方程式中的 $v$ 是每個分子平均占有空間的大小（可以理解成分子平均「勢力範圍」的大小），而 $b'$ 則為單個分子本身「包含」的體積（若為單原子分子如稀有氣體， $b'$ 就是原子半徑內包含的體積）。

## 適用範圍

范氏方程式對氣 - 液臨界溫度以上流體性質的描寫優於理想氣體方程式。對溫度稍低於臨界溫度的液體和低壓氣體也有較合理的描述。

但是，當描述對象處於狀態參量空間 (P,V,T) 中氣液相變區（即正在發生氣液轉變）時，對於固定的溫度，氣相的壓力恆為所在溫度下的飽和蒸氣壓，即不再隨體積 $V$ （嚴格地說應該是單位質量氣體占用的體積，即比容）變化而變化，所以這種情況下范氏方程式不再適用。

## 方程式的提出

下面以理想氣體狀態方程式為基礎，推導范氏方程式。若把氣體視為由體積無限小、相互之間無作用力的分子組成，這種模型便是理想氣體模型，與其相對應的狀態方程式是：

　
: $p={\frac {kT}{v}}$

若拋棄前一個的假設，把組成氣體的分子視為有一定大小的剛性球（其半徑稱為范德華半徑），用 $b$ 表示這些「球」的體積，上面的方程式便改寫為：

　
: $p={\frac {kT}{v-b}}$

在這裡，每個分子的「占有體積」 $v$ 被所謂「排斥體積」 $v-b$ 代替，反映了分子在空間中不能重疊。若氣體被壓縮至體積接近分子體積之和（即分子間空隙 $v-b$ 趨向於 0），那麼其壓力將趨於無窮大。

下一步，我們考慮原子對之間的引力。引力的存在會使分子的平均亥姆霍茲自由能下降，減少量正比於流體的密度。但壓力的大小滿足熱力學關係

　
: $p=-\left({\frac {\partial A^{*}}{\partial v}}\right)_{T}$

式中 _A\*_ 為每個分子的亥姆霍茲自由能。由此得到，引力使壓力減小的量正比於 ${\frac {1}{v^{2}}}$ 。記該比例常數為 $a$ ，可得

　
: $p={\frac {kT}{v-b}}-{\frac {a}{v^{2}}}$

這便是范氏方程式。

## 與理想氣體方程式模擬結果的比較

### 低壓狀況

在氣體壓力不太高的情況下，以下事實成立：

* 排斥體積 $b$ 的影響相對 $V$ 而言極小，可以忽略；以二氧化碳（ ${\ce {CO2}}$ ）為例，在標準狀況（0°C，1 標準大氣壓）下，1mol ${\ce {CO2}}$ 體積 $V$ 為 22414cm³，而相應的 $b=43$ cm³，比 $V$ 小 3 個數量級；
* 分子間的距離足夠大， ${\frac {a}{V^{2}}}$ 項完全可以視為 0；譬如在一大氣壓下二氧化碳氣體的 ${\frac {a}{V^{2}}}$ 值只有 7‰。

所以此時理想氣體方程式是范氏方程式（也是對實際氣體行為的）的一個良好近似。

### 中高壓狀況

隨著氣體壓力的增加，范氏方程式和理想氣體方程式結果的差別會變得十分明顯（左圖為 ${\ce {CO2}}$ 分別用理想氣體方程式和凡得瓦方程式模擬的 p-V 等溫線，溫度 70 °C）：

* 在壓力為 5000~15000kPa（50~150 標準大氣壓）的中壓區，由於體積被「壓小」導致分子間距靠近，分子間的引力（表現為 ${\frac {a}{V^{2}}}$ 項）變得不可忽略。 ${\frac {a}{V^{2}}}$ 項的存在使得氣體的壓力比不考慮分子間引力的理想氣體模型估計結果要小（所以左圖的中壓區里紅線比藍線要低）。
* 在壓力為 15000kPa 以上的高壓區，體積的急劇壓縮致使 $b$ 的影響不可忽略，於是范氏方程式中的體積項 $V-Nb$ （或比容項 $v-b$ ）將比理想氣體方程式中的體積項要小（或者說：對應相同體積 / 比容值的壓力項會升高）。這一效應導致在高壓區范氏氣體的狀態線重新趕上並超過理想氣體線（見左圖的左上角）。

## 用范氏方程式描述氣體的液化

范氏方程式適用於氣體的液化過程。氣體液化可能發生的**最高**溫度稱為**臨界溫度**，用 $T_{C}$ 表示：

* 當溫度 $T>T_{C}$ 時，無論給氣體施加多大的壓力都無法將它液化；
* 當溫度 $T<T_{C}$ 時，氣體可在壓力大於一定值時液化，且這一壓力隨著溫度 $T$ 下降而下降；

右圖所示為用范氏方程式模擬的 ${\ce {CO2}}$ 在不同溫度下的 _p-V_ 等溫線，從中可以明顯看出范氏方程式對液化過程的模擬（注意：若用理想氣體狀態方程式作上述模擬，得到的只是一系列雙曲線，因為在等溫條件下理想氣體狀態方程式就退化為玻意耳 - 馬略特定律 ——_pV_= 常數）。 ${\ce {CO2}}$ 氣體的臨界溫度為 $T_{C}=31$ °C $=304$ K。

* 70°C 時的曲線（右圖中藍線）形狀仍與玻意耳定律的結果（雙曲線）類似，儘管位置要略低；
* 當溫度下降到 40 °C，曲線（右圖中右二的曲線）形狀發生明顯的變化，表現為兩個拐點的出現。但此時二氧化碳仍然以氣態存在；
* 溫度進一步降至臨界溫度 31°C（圖中紅線），若此時氣體受壓至體積小於某定值 $V_{C}$ （隨溫度變化而變化），則氣體將發生液化。圖中 $V>V_{C}$ 時曲線對應氣態 ${\ce {CO2}}$ 的 $p$ ， $V$ 值， $V<V_{C}$ 時曲線對應液態 ${\ce {CO2}}$ 的 _$p$_ ， $V$ 值；
* 圖中 13°C 和 21°C 對應的曲線只有兩拐點以外的部分是與物理實際相符的。當氣體被進一步壓縮至比右拐點對應體積更小時，氣體將進入液化區，在液化過程中實際氣體的 p-V 線應是一段「平台」，而不是如圖所示的「駝峰」型。但完全液化後，液態 ${\ce {CO2}}$ 的壓力卻仍能被圖中曲線恰當地反映，此時曲線隨體積的減小而劇烈上升，這一定程度上反映了液體的不可壓縮性。另外，我們從圖中能得到的另一個資訊就是「液化平台」的長度隨溫度的下降而增加；

氣體的臨界狀態參量 $V_{C}$ 、 $p_{C}$ 、 $T_{C}$ 和范德華常數 $a$ 、 $b$ 之間存在下列數學關係：

　
: ${{V}_{C}}={3Nb}$

　
: ${p}_{C}={a}/{27}{b}^{2}$

　
: ${T}_{C}={8a}/{27bR}$

我們可以利用這些關係通過測出氣體的 _$T_{C}$_ 和對應的 _$p_{C}$_ 來得到 $a$ 和 $b$ 的值（由於測量上的困難，一般不使用 _$V_{C}$_ ）。

## 其他熱力學參量

下面，我們不再考慮 $v={\frac {V}{N}}$ （ $N$ 為系統中的分子數），改為考慮總體體積 $V$ 。

狀態方程式並不能告訴我們系統的所有熱力學參量。我們可以照搬上面推導范氏方程式的思路，從理想氣體的亥姆霍茲自由能表達式出發，推得下面的結論：

　
: $A(T,V,N)=-NkT\left(1+\ln \left({\frac {(V-Nb)T^{{\hat {c}}_{V}}}{N\Phi }}\right)\right)-{\frac {aN^{2}}{V}}$

式中 $A$ 為亥姆霍茲自由能， ${\hat {c}}_{v}$ 是無因次的定容熱容， $\Phi $ 是待定的熵常數。上述方程式將 $A$ 用它的自然變量 $V$ 和 $T$ 表示，所以系統的所有熱力學資訊已全部知道。其力學狀態方程式就是前面導出的范氏方程式

　
: $P=-\left({\frac {\partial A}{\partial V}}\right)_{T}={\frac {NkT}{V-Nb}}-{\frac {aN^{2}}{V^{2}}}$

系統的熵 ( $S$ ) 由下式決定

　
: $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]$

綜合 $A$ 和 $S$ 的表達式，可由定義得到系統內能

　
: $U=A+TS={\hat {c}}_{V}\,NkT-{\frac {aN^{2}}{V}}$

其他熱力學勢和化學勢也可用類似的方程式給出，但任何勢函數若要用壓力 $P$ 表示都需要求解一個三階多項式，使結果的形式變得很繁雜。所以，將焓和吉布斯能用它們相應的自然變量表示的結果都是複雜的（因為 $P$ 是它們的自然變量之一）。

## 簡化形式

雖然在一般形式的范氏方程式中，常數 $a$ 和 $b$ 因氣體 / 流體種類而異，但我們可以通過改變方程式的形式，得到一種適用於所有氣體 / 流體的普適形式。

按照下面的方式定義**約減變量**（亦稱**折合變量**，就是把變量轉換成其無因次形式），其中下標 $R$ 表示約減變量，下標 $C$ 表示原變量的臨界值：

　
: $p_{R}={\frac {p}{p_{C}}}$ ,

　
: $v_{R}={\frac {v}{v_{C}}}$ ,

　
: $T_{R}={\frac {T}{T_{C}}}$ ,

式中 $p_{C}={\frac {a}{27b^{2}}}$ ， $\displaystyle {v_{C}=3b}$ ， $kT_{C}={\frac {8a}{27b}}$ 。

用約減變量代替原變量，范氏方程式形式變為

　
: $\left(p_{R}+{\frac {3}{v_{R}^{2}}}\right)(v_{R}-1/3)={\frac {8}{3}}T_{R}$

這就是范氏方程式的不變形式，即這一形式不會因應用流體種類改變而改變。

上述方程式的不變性質亦稱對應狀態原理。

## 在可壓縮流動中的應用

在流體力學中，范氏方程式可以作為可壓縮流體（如液態高分子材料）的 PVT 狀態方程式。這種情況下，由於比容 $V$ 變化不大，可將方程式簡化為：

 $(p+A)(V-B)=CT\,$ ,

其中 $p$ 為壓力， $V$ 為比容， $T$ 為溫度， $A$ 、 $B$ 、 $C$ 均為與對象相關的參數。

## 注釋

## 參考文獻

* 趙凱華、羅蔚茵 著：《新概念物理教程・熱學》（第二版），高等教育出版社，北京，2005，ISBN 7-04-017680-7
* 汪志誠 編：《熱力學・統計物理》（第二版），高等教育出版社，北京，1993，ISBN 7-04-004360-2

## 外部連結

* 不同氣體的范氏常數 _a_ 和 _b_

## 參見

* 狀態方程式
* 氣體定律
* 范德華常數表
* 范德華力
* 范德華半徑