方程類型
Van der Waals 狀態方程
1873年J.D. Van der Waals提出了第一個適用於真實氣體的狀態方程,可寫為:
![立方型狀態方程](/img/d/654/wZwpmLyYTO2QjN3EjM3QTN1UTM1QDN5MjM5ADMwAjMwUzLxIzL4czLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
其中a,b是各種物質特有的常數,與其臨界參數有關。
Van der Waals方程是最簡單的立方型狀態方程。該方程能定量地描述流體的P-V-T關係,能夠同時描述氣、液兩相的性質。雖然精確度不高,但還是特別值得關注,因為它建立方程的推理方法對立方型方程、對比態原理以及後來與之類似的狀態方程的開發有著巨大的貢獻。
Redlich Kwong 方程
該方程由Redlich和Kwong於1949年提出,簡稱為RK方程,其形式如下:
![立方型狀態方程](/img/1/cb9/wZwpmLxITM0IzN4kjM3QTN1UTM1QDN5MjM5ADMwAjMwUzL5IzL2gzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
其中a,b為RK常數,與流體的特性有關。
RK方程與Van der Waals方程的區別僅在於壓力修正項的形式不同。
Soave Redlich Kwong方程
SRK方程,其形式為:
![立方型狀態方程](/img/f/9ab/wZwpmL2cjM2EDMwATM3QTN1UTM1QDN5MjM5ADMwAjMwUzLwEzL4gzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/1/ab8/wZwpmLwQDM1czN2gzM3QTN1UTM1QDN5MjM5ADMwAjMwUzL4MzLxgzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/d/1ae/wZwpmLxIzNxUzN2cDO2UzM1UTM1QDN5MjM5ADMwAjMwUzL3gzLxIzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/6/004/wZwpmL3ETM2UDN2IjM3QTN1UTM1QDN5MjM5ADMwAjMwUzLyIzL1UzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/c/440/wZwpmL0UjN3IzM0kjN0YzM1UTM1QDN5MjM5ADMwAjMwUzL5YzLxEzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
其中 是臨界參數有關的常數; 是對比溫度 的函式,其中的參數F可由純物質的飽和蒸氣壓和液體的密度數據關聯得到,對於非極性或弱極性物質,F可以通過下式進行計算:
![立方型狀態方程](/img/1/43f/wZwpmL0QDM4ADM1gTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL4EzLyAzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
其中ω為偏心因子。
與RK方程相比,SRK方程顯示出很大的優越性,特別是用它來計算純烴和烴類混合物系統的氣液平衡具有較高的精度,該方程在烴套用方面做出了很大的貢獻。
Peng Robinson 方程
RK和SRK方程有一個共同的不足,就是預測液相摩爾體積時精度相當差。為了彌補這一不足,Peng和Robinson 於1976年提出了如下形式的狀態方程:
![立方型狀態方程](/img/2/70c/wZwpmLzEjM5kTM3cjM3QTN1UTM1QDN5MjM5ADMwAjMwUzL3IzLwYzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/0/546/wZwpmL2ADM0UTN3cjM3QTN1UTM1QDN5MjM5ADMwAjMwUzL3IzLwczLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
。
![立方型狀態方程](/img/6/004/wZwpmL3ETM2UDN2IjM3QTN1UTM1QDN5MjM5ADMwAjMwUzLyIzL1UzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/c/440/wZwpmL0UjN3IzM0kjN0YzM1UTM1QDN5MjM5ADMwAjMwUzL5YzLxEzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
PR方程中的 是對比溫度 的函式,其中的參數F可由純物質的飽和蒸氣壓與飽和液體的密度數據關聯得到,對於非極性或弱極性物質,F亦可以通過下式計算:
![立方型狀態方程](/img/6/388/wZwpmLxATNzQzM0czM3QTN1UTM1QDN5MjM5ADMwAjMwUzL3MzLzUzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
。
![立方型狀態方程](/img/e/4ad/wZwpmLxYTNxczN5cTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL3EzLzMzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
PR方程的臨界壓縮因子 ,該值比RK方程的0.333有明顯改進,但仍偏離真實流體的數值。PR方程計算飽和蒸氣壓、飽和液體密度和氣液平衡的準確度均高於SRK方程。值得指出的是,SRK方程和PR方程在預測流體的蒸氣壓時顯示出突出的優勢,其主要原因在於溫度函式a(T)的表達式。
Patel Teja 方程
1982年,由Patel 和Teja導出了Patel Teja方程,其形式為:
![立方型狀態方程](/img/d/3b3/wZwpmL2MDM0YjMzgTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL4EzL1AzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/8/c3a/wZwpmL3QjM2MjM4cjM3QTN1UTM1QDN5MjM5ADMwAjMwUzL3IzL2czLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/b/b0b/wZwpmL4ITMxIDM3cTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL3EzL4YzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
其中 的計算方法如下:
![立方型狀態方程](/img/e/281/wZwpmLxAzN4QDN4cTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL3EzL0gzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
![立方型狀態方程](/img/5/1ac/wZwpmLyITN0UDO3czM3QTN1UTM1QDN5MjM5ADMwAjMwUzL3MzLzYzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/8/c6d/wZwpmL2MjNzUTOxkjM3QTN1UTM1QDN5MjM5ADMwAjMwUzL5IzL1AzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
是下列方程的最小正根:。
![立方型狀態方程](/img/8/a37/wZwpmLxQDM3cTNzcjM3QTN1UTM1QDN5MjM5ADMwAjMwUzL3IzL1IzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
這兩個參數是由純物質的飽和蒸氣壓和飽和密度數據關聯得到的,對非極性或弱極性物質也可由以下兩個式子求得:
![立方型狀態方程](/img/7/7e0/wZwpmLyMDO3UTM5kTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL5EzLyczLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/c/10f/wZwpmLzUjM4UTO0ATM3QTN1UTM1QDN5MjM5ADMwAjMwUzLwEzL4gzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
。
通用形式
一般式
立方型狀態方程可歸納為如下形式:
![立方型狀態方程](/img/7/1f9/wZwpmLyYjM1MjM1ATM3QTN1UTM1QDN5MjM5ADMwAjMwUzLwEzL4QzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
其中,對於不同的立方型狀態方程,m,n取不同的值,並且滿足如下關係:
![立方型狀態方程](/img/7/a0a/wZwpmL0YTM1cjN0ETM3QTN1UTM1QDN5MjM5ADMwAjMwUzLxEzL3gzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/d/80b/wZwpmL4QTO1kDOyAzM3QTN1UTM1QDN5MjM5ADMwAjMwUzLwMzL0gzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/3/7df/wZwpmLxYTMyIjN5EzM3QTN1UTM1QDN5MjM5ADMwAjMwUzLxMzLzYzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/0/f4b/wZwpmLxMTNxkzN2gTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL4EzL2UzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
其中 為關聯常數, 均為與臨界性質有關的參數。
各方程對應參數值
![立方型狀態方程](/img/2/0d5/wZwpmLxIjM1MDOxkjM3QTN1UTM1QDN5MjM5ADMwAjMwUzL5IzLzQzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
(1)Van der Waals方程:m=0,n=0, ;
![立方型狀態方程](/img/a/a00/wZwpmL4gDN5ETO2cDM3QTN1UTM1QDN5MjM5ADMwAjMwUzL3AzLxIzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
(2)Redlich Kwong 方程:m=b,n=0, ;
![立方型狀態方程](/img/7/885/wZwpmLyUTO1QzM1gTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL4EzL0gzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
(3)Soave Redlich Kwong方程:m=b,n=0, ;
![立方型狀態方程](/img/a/286/wZwpmL0EzN0EDNycjM3QTN1UTM1QDN5MjM5ADMwAjMwUzL3IzL0czLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![立方型狀態方程](/img/7/885/wZwpmLyUTO1QzM1gTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL4EzL0gzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
(4)Peng Robinson 方程:m=2b, , ;
![立方型狀態方程](/img/7/885/wZwpmLyUTO1QzM1gTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL4EzL0gzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
(5)Patel Teja 方程:m=b+c,n=-bc, 。
方程的求解
雖然立方型狀態方程可以用解析法求解三個體積根,但工程計算通常採用疊代法。下面具體介紹疊代過程,以Peng Robinson 方程為例進行討論。
已知:
![立方型狀態方程](/img/f/fe2/wZwpmLxcjN4IDMyAzM3QTN1UTM1QDN5MjM5ADMwAjMwUzLwMzL2czLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
為便於疊代,將上式恆等變形,得到:
![立方型狀態方程](/img/1/a71/wZwpmLyQTO2AzN4MzM3QTN1UTM1QDN5MjM5ADMwAjMwUzLzMzL0IzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
利用上式進行疊代計算時,體積根的初值設定方法為:
![立方型狀態方程](/img/c/f4a/wZwpmLwQTMxYDMwkTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL5EzLwEzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
(1)對於氣相,取初值 ,即以理想氣體作為初值;
![立方型狀態方程](/img/a/bdf/wZwpmL0cDN4MTOyQTM3QTN1UTM1QDN5MjM5ADMwAjMwUzL0EzLyAzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
(2)對於液相,取初值 。
計算出初值後代入上式右端,在疊代收斂後即可得到氣相摩爾體積或液相摩爾體積。
套用
立方型狀態方程的套用:
(1)用一個EOS即可精確地代表相當廣泛範圍內的實驗數據,藉此可精確計算所需的數據;
(2)EOS具有多功能性,除了pVT性質之外,還可用最少量的數據計算流體的其他熱力學函式、純物質的飽和蒸氣壓、混合物的氣一液相平衡、液一液相平衡.尤其是高壓下的相平衡計算;
(3)在相平衡計算中用一個EOS可進行兩相、三相的平衡數據計算,狀態方程中的混合規則與相互作用參數對各相使用同一形式或同一數值,計算過程簡捷、方便。