基于IAPWS-IF97水和蒸汽體膨脹系數(shù)和壓縮系數(shù)

摘 要

摘要：根據(jù)國際水和水蒸氣性質(zhì)協(xié)會(huì)(IAPWS)1997年工業(yè)用計(jì)算標(biāo)準(zhǔn)(簡(jiǎn)稱IAPWS-IF97)中的狀態(tài)方程，推導(dǎo)出水和蒸汽的體膨脹系數(shù)和壓縮系數(shù)的計(jì)算公式。該公式具有較高的精度，滿足工
摘要：根據(jù)國際水和水蒸氣性質(zhì)協(xié)會(huì)(IAPWS)1997年工業(yè)用計(jì)算標(biāo)準(zhǔn)(簡(jiǎn)稱IAPWS-IF97)中的狀態(tài)方程，推導(dǎo)出水和蒸汽的體膨脹系數(shù)和壓縮系數(shù)的計(jì)算公式。該公式具有較高的精度，滿足工程計(jì)算的需要，提高了工程計(jì)算效率。
關(guān)鍵詞：蒸汽；體膨脹系數(shù)；壓縮系數(shù)；IAPWS-IF97標(biāo)準(zhǔn)
Volumetric Expansion Coefficient and Compressibility Coefficient of Water and Steam Based on IAPWS-IF97
LI Jia，SHI Xiumin
Abstract：Based on the state equation in IAPWS industrial formulation 1997(IAPWS-IF97)，the calculation tormula of volumetric expansion coefficient and compressibility coefficient of water and steam are deduced.With a high accuracy，the formulations meet engineering caleulation requirements and improve the engineering calculation efficiency.
Key words：steam；volumetric expansion coefficient；compressibility coefficient；IAPWS industrial formu]ation 1997(IAPWS-IF97)
   水和蒸汽的體膨脹系數(shù)及壓縮系數(shù)在傳熱學(xué)、工程熱力學(xué)和熱工水力計(jì)算中處于非常重要的地位。本文依據(jù)國際水和水蒸氣性質(zhì)協(xié)會(huì)(IAPWS)提供的1997年工業(yè)用計(jì)算標(biāo)準(zhǔn)(簡(jiǎn)稱IAPWS-IF97)中的狀態(tài)方程，推導(dǎo)出了水和蒸汽的體膨脹系數(shù)及壓縮系數(shù)的計(jì)算公式。IAPWS-IF97中的公式計(jì)算精度非常高，并且能夠極大限度地減少計(jì)算時(shí)間。因此，根據(jù)IAPWS-IF97推導(dǎo)出的計(jì)算公式計(jì)算所得的結(jié)果，比根據(jù)蒸汽表內(nèi)插值法計(jì)算出的結(jié)果具有更高的精度，而且提高了計(jì)算效率。
1 體膨脹系數(shù)和壓縮系數(shù)的有效區(qū)域
   IAPWS-IF97可以計(jì)算的有效區(qū)域由一個(gè)低溫區(qū)域和一個(gè)附加的高溫區(qū)域組成，其中低溫區(qū)域蒸汽溫度(T)和壓力(p)的范圍是273.15K≤T≤1073.15K，p≤100MPa；高溫區(qū)域蒸汽溫度和壓力的范圍是1073.15K≤T≤2273.15K，p≤10MPa。這也是本文推導(dǎo)出的計(jì)算公式所能計(jì)算的有效區(qū)域。本文中所有的壓力均為絕對(duì)壓力。
    IAPWS-IF97將整個(gè)有效區(qū)域劃分為5個(gè)子區(qū)域(見圖1)。
 

    圖1中，區(qū)域1為常規(guī)水區(qū)，區(qū)域2為常規(guī)蒸汽區(qū)，區(qū)域3為臨界水和臨界蒸汽區(qū)，區(qū)域4為飽和水和飽和蒸汽區(qū)(即蒸汽飽和曲線，曲線延伸至區(qū)域3中，圖中小圈表示包含端點(diǎn))，區(qū)域5為高溫區(qū)。區(qū)域2與區(qū)域3的分界由B23方程[1]確定。
2 體膨脹系數(shù)和壓縮系數(shù)的計(jì)算公式
根據(jù)工程熱力學(xué)的定義，體膨脹系數(shù)和壓縮系數(shù)分別是等壓情況下比體積隨溫度的變化率和等溫情況下比體積隨壓力的變化率，計(jì)算公式為：
 
式中αp——體膨脹系數(shù)，K-1
    v——比體積，m3/kg
    T——熱力學(xué)溫度，K
p——壓力，Pa
 
式中κT——壓縮系數(shù)，Pa-1
    在IAPWS-IF97中，區(qū)域1、2和5給出的基本方程是Gibbs比自由能方程g(p，T)(對(duì)應(yīng)的無量綱形式的Gibbs比自由能方程為γ(π，τ))，區(qū)域3給出的基本方程是Helmholtz比自由能方程f(ρ，T)(對(duì)應(yīng)的無量綱形式的Helmholtz比自由能方程為φ(δ，τ))，區(qū)域4給出的基本方程是飽和壓力方程pS(T)和反推的飽和溫度方程TS(p)，同時(shí)給出了基本方程的微分計(jì)算方程和熱力學(xué)參數(shù)關(guān)系方程。根據(jù)不同區(qū)域的方程，得到的體膨脹系數(shù)和壓縮系數(shù)計(jì)算公式也各不相同。
2.1 區(qū)域1中體膨脹系數(shù)和壓縮系數(shù)的計(jì)算公式
    ①體膨脹系數(shù)的計(jì)算公式
在區(qū)域1中，γ(π，τ)以折算壓力π和折算溫度τ為獨(dú)立變量。比體積與基本方程之間的關(guān)系表達(dá)式為：
 
式中R——水的比氣體常數(shù)，kJ/(kg·K)，取0.461526kJ/(kg·K)
    p0——基準(zhǔn)壓力，Pa
    T0——基準(zhǔn)溫度，K
    區(qū)域1中，p0=16.53MPa，T0=1386K。
根據(jù)體膨脹系數(shù)的定義，體膨脹系數(shù)的計(jì)算公式推導(dǎo)如下：
 
式中γπ——γ(π，τ)對(duì)折算壓力的偏導(dǎo)數(shù)
    γπτ——γ(π，τ)對(duì)折算壓力和折算溫度的二階偏導(dǎo)數(shù)
將式(3)代入式(1)，得到體膨脹系數(shù)的計(jì)算公式：
 
式中ni，Ii，Ji——計(jì)算參數(shù)
    本文各區(qū)域中的ni，Ii，Ji均參見文獻(xiàn)[1]中相應(yīng)區(qū)域的取值。
    ② 壓縮系數(shù)的計(jì)算公式
根據(jù)壓縮系數(shù)的定義，壓縮系數(shù)的計(jì)算公式推導(dǎo)如下：
 
式中γππ——γ(π，τ)對(duì)折算壓力的二階偏導(dǎo)數(shù)
將式(4)代入式(2)，得到壓縮系數(shù)的計(jì)算公式：
 
    區(qū)域1的計(jì)算壓力、溫度范圍為：
    pS(T)≤p≤100MPa
    273.15K≤T≤623.15K
式中pS(T)——溫度為T時(shí)的飽和壓力，MPa
2.2 區(qū)域2中體膨脹系數(shù)和壓縮系數(shù)的計(jì)算公式
   ① 體膨脹系數(shù)的計(jì)算公式
在區(qū)域2中，γ(π，τ)由理想氣體部分和剩余部分的無量綱Gibbs比自由能組成，以折算壓力π和折算溫度τ為獨(dú)立變量。比體積與基本方程的關(guān)系表達(dá)式如下：
 
區(qū)域2中p0=1MPa，T0=540K。
根據(jù)體膨脹系數(shù)的定義，體膨脹系數(shù)的計(jì)算公式推導(dǎo)如下：
 
將式(5)代入式(1)，得到體膨脹系數(shù)的計(jì)算公式：
 
② 壓縮系數(shù)的計(jì)算公式
根據(jù)壓縮系數(shù)的定義，壓縮系數(shù)的計(jì)算公式推導(dǎo)如下：
 
將式(6)代入式(2)，得到壓縮系數(shù)的計(jì)算公式：
 
    區(qū)域2中的計(jì)算壓力、溫度范圍為：
  
式中p(T)B23——溫度為T時(shí)，由B23方程所確定的區(qū)域2、3的邊界壓力，MPa
2.3 區(qū)域3中體膨脹系數(shù)和壓縮系數(shù)的計(jì)算公式
    ① 體膨脹系數(shù)的計(jì)算公式
    在區(qū)域3中，φ(δ，τ)以折算密度δ和折算溫度τ為獨(dú)立變量。壓力的關(guān)系表達(dá)式如下：
 
    區(qū)域3中，ρ0=322kg/m3，T0=540K。
由于在區(qū)域3中給出的基本方程不是直接以p，T為獨(dú)立變量，因此，要使用隱函數(shù)求導(dǎo)的方法，將p看成是相對(duì)于T獨(dú)立的變量，然后將壓力的關(guān)系方程兩邊同時(shí)對(duì)T求偏導(dǎo)數(shù)，即：
 
解方程(7)可得：
 
將式(8)代入式(1)，得到體膨脹系數(shù)的計(jì)算公式：
 
式中φδτ——φ(δ，τ)對(duì)折算密度和折算溫度的二階偏導(dǎo)數(shù)
    φδδ——φ(δ，τ)對(duì)折算密度的二階偏導(dǎo)數(shù)
    ② 壓縮系數(shù)的計(jì)算公式
與體膨脹系數(shù)的計(jì)算方法類似，將壓力的關(guān)系方程兩邊同時(shí)對(duì)p求偏導(dǎo)數(shù)，即：
 
解方程(9)可得：
 
將式(10)代入式(2)，得到壓縮系數(shù)計(jì)算公式：
 
    上式中，φδ，φδδ的計(jì)算方法與本區(qū)域體膨脹系數(shù)計(jì)算公式中的計(jì)算方法相同。
    區(qū)域3中的計(jì)算壓力、溫度范圍為：
    p(T)B23≤p≤100MPa
    623.15K＜T≤T(p)B23
式中p(T)B23——溫度為T時(shí)，由B23方程確定的壓力，MPa
    T(p)B23——壓力為p時(shí)，由B23方程確定的溫度，K
2.4 區(qū)域4中體膨脹系數(shù)和壓縮系數(shù)的計(jì)算方法
    在區(qū)域4中，IAPWS-IF97給出的是飽和壓力方程pS(T)和飽和溫度方程TS(P)。因此，當(dāng)分別求該區(qū)域中的飽和蒸汽和飽和水的體膨脹系數(shù)和壓縮系數(shù)時(shí)，可根據(jù)飽和壓力方程和飽和溫度方程，在已知p或T時(shí)，求出飽和時(shí)的另1個(gè)參數(shù)，確定計(jì)算參數(shù)所處的實(shí)際區(qū)域。
    當(dāng)溫度T≤623.15K時(shí)，飽和水的體膨脹系數(shù)和壓縮系數(shù)采用區(qū)域1中的計(jì)算公式進(jìn)行計(jì)算，飽和蒸汽的體膨脹系數(shù)和壓縮系數(shù)采用區(qū)域2中的計(jì)算公式進(jìn)行計(jì)算。
    當(dāng)溫度T＞623.15K時(shí)，處于區(qū)域3中，此時(shí)在計(jì)算飽和水和飽和蒸汽的體膨脹系數(shù)和壓縮系數(shù)時(shí)，應(yīng)分別計(jì)算出飽和水和飽和蒸汽的比體積，然后根據(jù)區(qū)域3中的計(jì)算公式進(jìn)行計(jì)算。
2.5 區(qū)域5中體膨脹系數(shù)和壓縮系數(shù)的計(jì)算方法
    在區(qū)域5中，γ(π，τ)由理想氣體部分和剩余部分的無量綱Gibbs比自由能組成，γ(π，τ)與比體積的關(guān)系表達(dá)式與區(qū)域2中的形式相同，因此，可以得到區(qū)域5中體膨脹系數(shù)和壓縮系數(shù)的計(jì)算公式。
體膨脹系數(shù)計(jì)算公式：
 
壓縮系數(shù)計(jì)算公式：
 
    區(qū)域5中，p0=1MPa，T0=1000K。
    區(qū)域5中的計(jì)算壓力和溫度范圍覆蓋了整個(gè)高溫區(qū)域，計(jì)算范圍為：
    0MPa＜p≤10MPa
    1073.15K＜T≤2273.15K
3 結(jié)論
    根據(jù)IAPWS-IF97給出的基本方程，推導(dǎo)出了水和蒸汽的體膨脹系數(shù)和壓縮系數(shù)計(jì)算公式。由于IAPWS-IF97中的公式精確度較高(比體積相對(duì)誤差為±0.01%～±0.3%)，并通過與實(shí)際數(shù)據(jù)進(jìn)行對(duì)比，得知本文所推導(dǎo)出的計(jì)算公式計(jì)算結(jié)果相對(duì)誤差在±0.5%之內(nèi)，完全能夠滿足工程計(jì)算的需要，同時(shí)能夠顯著提高工程計(jì)算效率。
參考文獻(xiàn)：
[1] WAGNER W，COOPER J R，DITFMANN A.et al.The IAPWS industrial formulation 1997 for the thermodynamic properties of water and steam[J].Eng.Gas Turbines Power，2000(122)：150-182.
 
(本文作者：李佳 史秀敏 中國石油大學(xué)(華東)儲(chǔ)運(yùn)與建筑工程學(xué)院 山東青島 266555)