李文博 莫善軍
(中山大學(xué) 廣州 510000)
在制冷空調(diào)領(lǐng)域,獲取濕空氣的物性數(shù)據(jù)是系統(tǒng)實驗、熱力學(xué)計算、仿真計算以及新型材料(如:親水涂層)開發(fā)的過程中必須解決的一個問題。制冷系統(tǒng)仿真計算時,由于空氣側(cè)換熱的影響,濕度對制冷系統(tǒng)能力、能效影響非常大。因此濕空氣物性參數(shù)的計算也有較高的要求。本文通過對CoolProp軟件的使用方法進行介紹,以常用了Excel軟件和Visual Basic. NET語言為平臺,提供使用動態(tài)鏈接庫(DLL)技術(shù)將CoolProp直接引入到Excel和桌面程序中,該方法具有精確、方便的有點,可以滿足仿真計算的要求。
空調(diào)制造業(yè)正在向著智能化邁進,在制冷空調(diào)仿真系統(tǒng)中,濕空氣的物理性質(zhì)計算有著非常重要的應(yīng)用,由于濕空氣特性的原因,傳統(tǒng)方法都在使用查表法或者使用擬合公式,由于模擬精度存在問題,無法滿足設(shè)計要求。
當(dāng)前行業(yè)內(nèi)通用的計算方法主要有三種:
數(shù)據(jù)庫插值法:該方法主要使用前人測量繪制的圖表,將中間狀態(tài)進行插值計算,速度慢且精度較低。使用查表法進行計算的方法,行業(yè)內(nèi)通常會選用ASHRAE HANDBOOK FUNDAMENTALS CHAPTER1 PSYCHROMETRICS[1]中提供的濕空氣參數(shù)表格。
該方法一般行業(yè)內(nèi)簡化方法使用的是擬合公式,該方法速度快,但是由于擬合公式參數(shù)范圍問題,無法大范圍的使用,精度也無法保證。
公式法計算濕空氣物性,一般行業(yè)內(nèi)選擇使用ASHRAE HANDBOOK FUNDAMENTALS CHAPTER1 PSYCHROMETRICS[1]中提供的系列公式進行計算,該系列公式在其定義域內(nèi)有較高精度,但是如果作為仿真計算使用,精度尚無法達到預(yù)期的水平。
式(1)定義了冰點以下的水的飽和蒸氣壓,其中:
式(2)定義了0-100 ℃范圍內(nèi)的水的飽和蒸氣壓,其中:
該組公式參數(shù)較多,對于實際工程計算時造成了較大的麻煩,不利于工程應(yīng)用。
如果軟件作者使用了良好的計算方式,該方法可以保證精度和速度,但是由于第三方作者的軟件通常不開放源碼,無法判斷其計算方式,只能從計算結(jié)果來評估精度,另外第三方軟件也無法與自身的仿真系統(tǒng)相結(jié)合,使得計算效率下降。
CoolProp 是一個開源物性參數(shù)獲取工具,其本身就封裝了各種常用工質(zhì)的物性參數(shù)。CoolProp允許使用多種第三方程序語言進行編程調(diào)用,支持的編程語言如下:
可以提供完整功能支持的語言: Python (2.x, 3.x) ,Modelica, Octave, C#, VB.net, MathCAD, Java, Android,MATLAB
僅適用于部分支持的語言: Labview, EES, Microsoft Excel, LibreOffice, Javascript, PHP, FORTRAN, Maple,Mathematica, Scilab, Delphi & Lazarus, Julia。
CoolProp提供的導(dǎo)出函數(shù)格式非常簡單易用,理論上任何一種支持第三方dll的語言都可以使用CoolProp引用計算。CoolProp使用MIT授權(quán),可以免費的用于商用工具,而無需向作者付費,你甚至可以選擇將CoolProp的源碼修改后編譯,發(fā)布屬于自己的物性數(shù)據(jù)庫。
CoolProp在對濕空氣進行計算時,可以選擇設(shè)置濕空氣的相對濕度,這點比Refprop有很大的優(yōu)勢,Refprop中的空氣計算模型使用的是只有氮氣/氬氣/氧氣的混合物模型,其質(zhì)量比(75.57/1.27/23.16),摩爾比(0.7812/0.0092/0.2096),這使得Refprop無法滿足濕空氣的熱物性計算。
CoolProp的官方網(wǎng)址:http://www.CoolProp.org/,包含豐富的官方文檔,本文部分內(nèi)容也取自其官方文章或官方Github開源說明。
Windows環(huán)境下,CoolProp提供編譯并打包好的二進制文件可以直接進行安裝,訪問地址https://sourceforge.net/projects/CoolProp/files/CoolProp/6.2.1/Installers/,版本變化時,地址略有不同。
程序安裝過程非常簡單,雙擊安裝程序后點擊確定。其安裝包會根據(jù)當(dāng)前系統(tǒng)版本自動選擇默認(rèn)的dll版本,以更優(yōu)的方法適配。
無論選擇哪一種安裝模式,最終所有可能用到的文件都會被安裝到你的電腦中,如果是Windows10系統(tǒng),可以在路徑:C:Users用戶名AppDataRoamingCoolProp當(dāng)中可以找到這些文件。其中CoolProp.dll動態(tài)鏈接庫文件與用戶選擇的操作安裝版本一致,即:如果選擇了64位的安裝,則CoolProp.dll就是64位的。CoolProp_stdcall.dll是CoolProp的32位動態(tài)鏈接庫。
程序安裝后會自動安裝好Office Excel的加載項,可以提供默認(rèn)的函數(shù)使用。示例文件和計算演示會自動在桌面生成一個名為TestExcel.xlsx文件。
絕大多數(shù)情況下,CoolProp安裝后可以直接使用,無需任何修改和變更。
CoolProp內(nèi)置常用的濕空氣計算公式見表1。
CoolProp用于濕空氣計算的函數(shù)名為HAPropSI,其主要輸入輸出參數(shù)的定義見表2。
需要注意的是,CoolProp中使用的參數(shù),計算單位均為國際單位制(SI),應(yīng)用于北美地區(qū)或者其他非國際單位制地區(qū)時,需要注意單位換算。
使用CoolProp內(nèi)置函數(shù)計算了一些濕空氣物性參數(shù)點,與ASHRAE HANDBOOK FUNDAMENTALS CHAPTER1 PSYCHROMETRICS[1]中提供的濕空氣參數(shù)表格進行對比,從表3、表4的對比情況來看,二者的取值結(jié)果在含濕量和飽和空氣焓值計算上非常接近,最大偏差0.043 %,精度滿足仿真計算要求。
CoolProp允許用戶在使用時自定義任意組分比例的混合物,以下方法與直接計算R410A密度結(jié)果等價
圖1 CoolProp的安裝引導(dǎo)界面
表1 CoolProp濕空氣常用計算公式
表2 HAPropSI輸入輸出參數(shù)定義
表3 CoolProp計算的含濕量偏差
當(dāng)使用Refprop的計算接口進行計算時,計算結(jié)果與Refprop計算結(jié)果一致,注意僅限制冷劑類型,Refprop無法計算濕空氣。
CoolProp官方提供了一種VB.net的引用方法,其原理是基于C#的接口,使用vb.net二次打包的,需要引用非常多的文件,調(diào)用起來并不方便。
本文提供的是另一種方法。基于Excel版本進行修改,使其能夠在VB.net平臺上運行,兼容性高,容易從Excel/VBA向VB.net遷移,同時代碼量更少,更簡潔,提供更簡潔的引用方法,目錄結(jié)構(gòu)清晰分明。
CoolProp提供了多個版本的dll可供引用,由于用戶的操作系統(tǒng)存在64位和32位兩種情況,推薦引用32位版本的dll。
CoolProp默認(rèn)安裝路徑為C:Users用戶名AppDataRoamingCoolProp,其中CoolProp_stdcall.dll為32位版本的dll,CoolProp_x64.dll為64位版本的dll,如圖2。
使用VS創(chuàng)建VB.NET窗體項目,插入模塊命名為CoolProp,復(fù)制Excel加載項中的vba代碼到模塊中即可完成代碼引用,注意需要將CoolProp_stdcall.dll的文件路徑改為實際程序所在位置的真實路徑。
使用方法與Excel中相似,演示如下:
筆者將CoolProp整合了一套完整的濕空氣計算軟件,可以供同事們?nèi)粘9ぷ魇褂?,反饋良好,界面見圖3。
表4 CoolProp計算的濕空氣比焓值偏差
圖2 Clloprop安裝目錄內(nèi)的文件
圖3 整合的物性計算軟件界面
文章介紹了CoolProp的安裝和使用,為在應(yīng)用程序中調(diào)用CoolProp的子程序提供了方便。使用表格的形式給出了濕空氣不同物理性質(zhì)的計算對比結(jié)果,表明CoolProp的計算精度是滿足要求的。
CoolProp是一款優(yōu)秀的商業(yè)友好軟件,其計算結(jié)果精度高,速度快,通過對該軟件的使用,可以快速的實現(xiàn)濕空氣的物性計算,且入門難度極低,大幅度的降低了濕空氣仿真計算難度,提升了仿真計算速度和精度。