多級(jí)空氣源熱泵系統(tǒng)仿真及試驗(yàn)研究*

姚陽(yáng)，孫杰，張華榮，宋恒，申凱

(國(guó)網(wǎng)江蘇省電力有限公司泰州供電分公司,江蘇泰州,225300)

0 引言

熱泵是一種利用外界做功將低位熱源轉(zhuǎn)化為高位熱源的高效制熱技術(shù),通過(guò)消耗少量電能得到大量熱能,具有綠色環(huán)保、節(jié)約等優(yōu)點(diǎn)符合綠色發(fā)展的要求[1-2]。因此,近年來(lái),熱泵技術(shù)引起越來(lái)越多的科學(xué)家關(guān)注空氣源熱泵式谷物烘干機(jī)在南方主要糧食產(chǎn)區(qū)也獲得越來(lái)越多的推廣應(yīng)用[3-4]。

巨永平等[5]對(duì)空氣源熱泵的經(jīng)濟(jì)性與節(jié)能效果進(jìn)行了分析,驗(yàn)證了空氣源熱泵具有明顯的環(huán)保性與節(jié)能性,證明了發(fā)展空氣源熱泵具有良好的經(jīng)濟(jì)價(jià)值。楊衛(wèi)芳等[6]基于試驗(yàn)臺(tái)研究了蒸發(fā)溫度、冷凝溫度等實(shí)驗(yàn)變量對(duì)熱泵性能的影響,對(duì)在不同條件下的熱泵性能做出了詳細(xì)的分析。蘇偉等[7]針對(duì)空氣源熱泵在低溫環(huán)境下容易結(jié)霜的問(wèn)題,提出了優(yōu)化換熱器結(jié)構(gòu),增加電磁場(chǎng)等方法進(jìn)行抑霜與除霜。Congedo等[8]調(diào)查了氣候變化對(duì)空氣源熱泵的影響,預(yù)測(cè)了未來(lái)各個(gè)地區(qū)空氣源熱泵的發(fā)展趨勢(shì)。Olympios等[9]運(yùn)用綜合熱網(wǎng)絡(luò)模型對(duì)空氣源熱泵系統(tǒng)進(jìn)行優(yōu)化,降低了運(yùn)營(yíng)成本,優(yōu)化了熱泵操作。由于熱泵系統(tǒng)的生產(chǎn)試驗(yàn)成本較高,很多學(xué)者通過(guò)仿真試驗(yàn),對(duì)系統(tǒng)的性能進(jìn)行研究。周富玉[10]建立了單、雙級(jí)壓縮空氣源熱泵系統(tǒng)的熱力學(xué)仿真模型,通過(guò)仿真,證明雙級(jí)壓縮模型性能強(qiáng)于單級(jí)壓縮,運(yùn)用雙級(jí)壓縮模型可以解決空氣源熱泵難以適應(yīng)環(huán)境的問(wèn)題。Shen等[11]建立空氣源熱泵仿真模型,通過(guò)模型進(jìn)行了建筑能源模擬,揭示空氣源熱泵節(jié)能和降低電力成本的潛力。劉業(yè)鳳等[12]基于KULI仿真建立CO2熱泵模塊與熱管理系統(tǒng)的整車(chē)模型,分析了余熱回收裝置對(duì)節(jié)能的重要作用。徐瑞成等[13]使用噴氣增焓技術(shù),解決了熱泵穩(wěn)定性差的問(wèn)題。牛建會(huì)等[14]對(duì)熱泵抑霜除霜技術(shù)進(jìn)行總結(jié),論述了空氣源熱泵除霜的重要性。張志剛等[15]對(duì)空氣源熱泵冬季性能進(jìn)行分析,尋找最好的節(jié)能效果。劉萌等[16]建立了熱泵仿真模型實(shí)現(xiàn)了對(duì)空氣源熱泵動(dòng)態(tài)調(diào)節(jié)控制。

谷物烘干機(jī)空氣源熱泵系統(tǒng)是一種近幾年發(fā)展起來(lái)的新型技術(shù)[17-18],國(guó)內(nèi)外對(duì)其研究比較少,關(guān)于多級(jí)谷物烘干機(jī)空氣源熱泵系統(tǒng)的仿真研究處于空白[19]。本文通過(guò)利用Simulink對(duì)系統(tǒng)的壓縮機(jī)、冷凝器、蒸發(fā)器、膨脹閥分別進(jìn)行仿真建模[20],然后根據(jù)熱泵的工作原理,建立整個(gè)熱泵系統(tǒng)的數(shù)學(xué)建模,最后通過(guò)上述模型進(jìn)行熱泵系統(tǒng)的性能的仿真研究,以期為空氣源熱泵系統(tǒng)的設(shè)計(jì)改進(jìn)和性能優(yōu)化提供理論參考。

1 空氣源熱泵谷物烘干機(jī)的組成

空氣源熱泵谷物烘干機(jī)的組成如圖1所示。

圖1 空氣源熱泵谷物烘干機(jī)組成圖

該機(jī)主要由循環(huán)式谷物烘干機(jī)(CPR-165型12噸谷物烘干機(jī)、額定風(fēng)量10 800 m3/h)、熱泵系統(tǒng)(自制,功率65 kW,額定制熱量90～150 kW)及控制系統(tǒng)三部分組成。其中,熱泵系統(tǒng)由4臺(tái)定頻壓縮機(jī)及1臺(tái)變頻壓縮機(jī),壓縮機(jī)型號(hào)見(jiàn)表1,共計(jì)5套制冷壓縮系統(tǒng)組成。

表1 壓縮機(jī)型號(hào)與相關(guān)參數(shù)Tab. 1 Compressor model and associated parameters

2 熱泵系統(tǒng)模型的建立

2.1 模型建立與假設(shè)

制冷壓縮系統(tǒng)主要由制冷壓縮機(jī)、冷凝器、蒸發(fā)器及膨脹閥構(gòu)成,通過(guò)分別建立壓縮機(jī)、冷凝器、蒸發(fā)器和膨脹閥的仿真模型,可以建立1組制冷壓縮系統(tǒng)的仿真模型,然后將5組制冷壓縮系統(tǒng)的仿真模型進(jìn)行組合,就可以構(gòu)成1套多級(jí)空氣源熱泵模型。在輸入冷凝溫度、蒸發(fā)溫度、過(guò)熱度、過(guò)冷度以及進(jìn)風(fēng)口風(fēng)量等參數(shù)后,便可以計(jì)算出熱泵機(jī)組的能效比、熱風(fēng)溫度、制熱量等參數(shù)。具體仿真流程如圖2所示。

圖2 多級(jí)熱泵仿真流程圖

為了簡(jiǎn)化試驗(yàn)步驟,提高試驗(yàn)的準(zhǔn)確性,本文做出如下假設(shè)。

假設(shè)1:假設(shè)壓縮機(jī)的壓縮過(guò)程為等熵過(guò)程。

假設(shè)2:假設(shè)膨脹閥節(jié)流過(guò)程為絕熱過(guò)程,節(jié)流前后焓值不變。

假設(shè)3:忽略熱泵工質(zhì)在管道內(nèi)流動(dòng)時(shí)的壓力損失。

假設(shè)4:忽略在冷凝器和蒸發(fā)器中的熱量損失。

2.2 壓縮機(jī)模型

為簡(jiǎn)化模型,本文只對(duì)壓縮機(jī)的制熱量與功率進(jìn)行建模討論,系統(tǒng)所用壓縮機(jī)參數(shù)與型號(hào)見(jiàn)表1。

由于壓縮機(jī)的壓縮過(guò)程是一個(gè)絕熱過(guò)程,壓縮過(guò)程中氣體的狀態(tài)方程和過(guò)程方程如式(1)～式(6)所示。

壓縮機(jī)的吸氣量

PV=RgT

(1)

式中:P——壓強(qiáng),Pa;

V——體積,m3;

Rg——?dú)怏w常量;

T——絕對(duì)溫度,℃。

壓縮機(jī)的容積比

(2)

式中:Vth——壓縮機(jī)吸氣容積,m3;

V2——壓縮機(jī)壓縮結(jié)束時(shí)的壓縮腔體積,m3;

v1——壓縮機(jī)吸氣結(jié)束時(shí)的比容,m3/kg;

v2——壓縮機(jī)壓縮結(jié)束時(shí)的比容,m3/kg。

壓縮機(jī)的吸氣量

(3)

fv——壓縮機(jī)的容積效率,取0.96;

ρ1——壓縮機(jī)吸氣孔出的制冷劑密度,kg/m3;

n——壓縮機(jī)的轉(zhuǎn)速,r/min。

壓縮機(jī)的壓力比

(4)

式中:P1——吸氣結(jié)束時(shí)出口的壓力,N;

P2——壓縮結(jié)束后的壓力,N;

k——制冷劑的等熵指數(shù)。

制冷劑在這一過(guò)程所做的功

(5)

式中:v——制冷劑流速,m/s;

dp——制冷劑通過(guò)的管徑長(zhǎng)度,m;

γ——制冷劑傳熱系數(shù);

p1——壓縮機(jī)吸氣結(jié)束時(shí)的制冷劑壓力,kPa;

p2——壓縮機(jī)壓縮結(jié)束時(shí)的制冷劑壓力,kPa。

壓縮機(jī)的質(zhì)量流量

(6)

式中:nv——壓縮機(jī)轉(zhuǎn)速,r/min;

λ——阻力系數(shù)。

根據(jù)上述計(jì)算公式,利用Simulink中的數(shù)學(xué)運(yùn)算、離散系統(tǒng)及函數(shù)調(diào)用等模塊,建立壓縮機(jī)的模型如圖3所示。其中,Te表示蒸發(fā)溫度,℃;Tc表示冷凝溫度,℃;mr表示制冷劑質(zhì)量流量,kg/h;V1表示壓縮機(jī)的工作容積,m3。

圖3 壓縮機(jī)建模

3.3 冷凝器模型

冷凝器在工作過(guò)程中可以分為過(guò)冷、兩相與過(guò)熱區(qū)三部分,每個(gè)階段的溫度變化規(guī)律如圖4所示。其中,tr,in表示制冷劑過(guò)熱區(qū)進(jìn)口溫度,℃;tr,1表示制冷劑兩相區(qū)進(jìn)口溫度,℃;tr,2表示制冷劑兩相區(qū)出口溫度,℃;tr,out表示制冷劑過(guò)冷區(qū)出口溫度,℃;ta,in表示空氣過(guò)冷區(qū)進(jìn)口溫度,℃;ta,1表示空氣兩相區(qū)進(jìn)口溫度,℃;ta,2表示空氣兩相區(qū)出口溫度,℃;ta,out表示空氣過(guò)熱區(qū)出口溫度,℃。

圖4 冷凝器工作原理

根據(jù)傳熱過(guò)程的規(guī)律與管翅式冷凝器的換熱關(guān)系,可以對(duì)換熱過(guò)程進(jìn)行計(jì)算,首先對(duì)過(guò)冷段進(jìn)行計(jì)算。

1) 過(guò)冷段

工質(zhì)側(cè)能量

Qr,sc=mout,com(hr,2-hr,out)

(8)

式中:mout,com——制冷劑進(jìn)口質(zhì)量流量,kg/h;

hr,2——制冷劑過(guò)冷區(qū)進(jìn)口焓值,J/kg,也是制冷劑兩相區(qū)出口焓值;

hr,out——制冷劑過(guò)冷區(qū)進(jìn)口焓值,J/kg。

空氣側(cè)能量

Qa,sc=ma(ha,1-ha,in)

(9)

式中:ma——空氣進(jìn)口質(zhì)量流量,kg/h;

ha,1——空氣過(guò)冷區(qū)出口焓值,J/kg,也是空氣兩相區(qū)進(jìn)口焓值;

ha,in——空氣過(guò)冷區(qū)進(jìn)口焓值,J/kg。

工質(zhì)端的平均溫度

(10)

空氣端平均溫度

(11)

2) 兩相段

工質(zhì)側(cè)能量

Qr,tp=mout,com(hr,1-hr,2)

(12)

式中:hr,1——制冷劑兩相區(qū)進(jìn)口焓值,J/kg,也是制冷劑冷卻段出口焓值。

空氣側(cè)能量

Qa,tp=ma(ha,2-ha,1)

(13)

式中:ha,2——空氣兩相區(qū)出口焓值,J/kg,也是空氣冷卻段進(jìn)口焓值。

工質(zhì)端的平均溫度

(14)

3) 冷卻段

工質(zhì)側(cè)能量

Qr,sh=mout,com(hr,in-hr,1)

(15)

式中:hr,in——制冷劑冷卻段進(jìn)口焓值,J/kg。

空氣側(cè)能量

Qa,sh=ma(ha,out-ha,2)

(16)

式中:ha,out——空氣冷卻段出口焓值,J/kg。

工質(zhì)端的平均溫度

(17)

空氣端平均溫度

(18)

4) 換熱系數(shù)

流動(dòng)中的空氣雷諾數(shù)

(19)

式中:DF——流體流速;m/s;

VF——特征尺寸,m;

vAC——流體運(yùn)動(dòng)粘度,m2/s。

空氣側(cè)換熱系數(shù)

(20)

式中:C1、C2、n、m——方程參數(shù),查表可得;

λAC——進(jìn)出口平均熱導(dǎo)率,W/(m·℃);

DE——當(dāng)量直徑,m;

L——翅片寬度,m。

工質(zhì)側(cè)凝結(jié)換熱系數(shù)

αCR=0.555Bc0.25DI-0.25(Tc-Tb)-0.25

(21)

式中:Bc——工質(zhì)特性參數(shù),查表可得;

DI——管直徑,m;

Tc——工質(zhì)冷凝溫度,℃;

Tb——管內(nèi)表面溫度,℃。

冷凝器傳熱公式

αCRAI(Tc-Tb)=ηEOαACAOF(Tb-TAC)

(22)

式中:ηEO——傳熱效率;

AI——冷凝器表面積,m2;

AOF——冷凝器翅片表面積,m2;

TAC——管外表面溫度,℃;

傳熱系數(shù)

(23)

式中:RTCO——外側(cè)傳熱表面的污垢熱阻。

根據(jù)上述公式,使用Simulink中的模塊對(duì)過(guò)冷區(qū)、兩相區(qū)、冷卻區(qū)進(jìn)行模型的建立,通過(guò)仿真得出出口溫度、冷凝器制熱量等參數(shù),模型如圖5所示。其中,Tc表示環(huán)境溫度,℃;Sc表示制冷劑進(jìn)口的熵,J/(kg·℃);tm,1表示冷卻段制冷劑出口溫度,℃;tm,2表示過(guò)冷段進(jìn)口溫度,℃;h3表示冷卻段出口焓值,J/kg;h6表示過(guò)冷段出口焓值,J/kg,也是膨脹閥進(jìn)口處的焓值。

圖5 冷凝器建模

3.4 蒸發(fā)器模型

蒸發(fā)器模型與冷凝器模型大致相同,把蒸發(fā)器分為過(guò)熱區(qū)與兩相區(qū),其工作時(shí)的具體過(guò)程如圖6所示。

圖6 蒸發(fā)器工作原理

此時(shí),tr,in表示制冷劑兩相進(jìn)口溫度,℃;tr,1表示制冷劑兩相區(qū)出口溫度,℃;tr,out表示制冷劑過(guò)熱區(qū)出口溫度,℃;ta,in表示空氣過(guò)熱區(qū)進(jìn)口溫度,℃;ta,1表示空氣兩相區(qū)進(jìn)口溫度,℃;ta,out表示空氣兩相區(qū)出口溫度,℃。

兩相區(qū)與過(guò)熱區(qū)公式與上述基本相同,不再贅述,僅介紹換熱系數(shù)。

液相表面換熱系數(shù)

(24)

式中:λL——進(jìn)出口平均熱導(dǎo)率,W/(m·℃);

D1——管內(nèi)徑,m;

ReL——飽和制冷劑雷諾數(shù);

PγL——普朗常數(shù)。

流動(dòng)中的空氣雷諾數(shù)

(25)

管內(nèi)的蒸發(fā)表面換熱系數(shù)

αER=αL[C1C2C0(25FrL)5+C3B0C4F0]

(26)

式中:FrL——內(nèi)外傳熱面的中間傳熱面的面積,m2;

C0——對(duì)流特征系數(shù);

B0——沸騰特征系數(shù);

F0——取決于工質(zhì)的系數(shù);

C3——方程參數(shù)。

通過(guò)研究冷凝器在過(guò)冷與兩相兩個(gè)工作狀態(tài),建立Simulink模型,仿真得出冷凝器出口溫度、冷凝器規(guī)格等物理量,模型如圖7所示。其中,Sh表示蒸發(fā)器制冷劑進(jìn)口焓值,J/kg;Te表示蒸發(fā)器制冷劑進(jìn)口溫度,℃;tz,1表示制冷劑過(guò)熱段進(jìn)口溫度,℃。

圖7 蒸發(fā)器仿真模型

3.5 膨脹閥模型

制冷劑通過(guò)膨脹閥時(shí)由高溫高壓的液體轉(zhuǎn)化為低溫低壓的蒸汽,其主要作用就是調(diào)節(jié)進(jìn)入蒸發(fā)器制冷劑的壓力、流量、物質(zhì)狀態(tài)等。在這個(gè)過(guò)程中,制冷劑與外界沒(méi)有沒(méi)有能量交換,故為等焓過(guò)程,故膨脹閥能量方程為

h6=h7

(27)

式中:h6——膨脹閥進(jìn)口處的焓值,kg/kJ;

h7——膨脹閥出口處的焓值,kg/kJ。

3.6 單級(jí)和多級(jí)熱泵系統(tǒng)模型

制冷劑依次經(jīng)過(guò)壓縮機(jī)、冷凝器、膨脹閥和蒸發(fā)器,在蒸發(fā)器吸熱,然后至冷凝器放熱,完成從環(huán)境吸熱對(duì)空氣加熱的過(guò)程。將上述壓縮機(jī)等模型串聯(lián)起來(lái),即構(gòu)成一個(gè)單級(jí)熱泵系統(tǒng)模型。

通過(guò)Simulink構(gòu)建多級(jí)熱泵模型如圖8所示。單極熱泵系統(tǒng)模型如圖9所示。

圖8 多級(jí)熱泵模型

圖9 單極熱泵系統(tǒng)模型

多級(jí)熱泵系統(tǒng)由5個(gè)制冷壓縮系統(tǒng)組成,采取逐級(jí)加熱,即五個(gè)熱泵機(jī)組相互串聯(lián),第五級(jí)冷凝器端的溫度便是熱泵的制熱溫度,出口的空氣溫度為熱風(fēng)溫度。其中,Q表示輸出的總熱量,kJ;h表示熱風(fēng)中的焓值,J/kg;COP表示熱泵性能系數(shù);W表示壓縮機(jī)做的總功,kW;T2表示熱風(fēng)溫度,℃;Sc2表示過(guò)熱度;Sh2表示熱空氣所具有的熵,J/(kg·℃)。

4 仿真試驗(yàn)及結(jié)果分析

在江蘇地區(qū),水稻收獲季節(jié)環(huán)境溫度多在0 ℃～20 ℃ 之間,為此,對(duì)熱泵系統(tǒng)在環(huán)境溫度為0 ℃、5 ℃、10 ℃、15 ℃及20 ℃時(shí),熱風(fēng)溫度、制熱量和COP的仿真結(jié)果與試驗(yàn)數(shù)據(jù)結(jié)果進(jìn)行對(duì)比分析,通過(guò)平均絕對(duì)誤差MAD和相對(duì)平均誤差MRE,對(duì)模型的性能進(jìn)行評(píng)價(jià)。對(duì)熱風(fēng)溫度、制熱量及COP的仿真與試驗(yàn)結(jié)果分別如圖10～圖12所示。

圖10 熱風(fēng)溫度對(duì)比圖

圖11 制熱量對(duì)比圖

圖12 COP仿真值與試驗(yàn)值

由圖10可知,在環(huán)境溫度為0 ℃時(shí),熱泵系統(tǒng)的實(shí)際出風(fēng)溫度可以達(dá)到55 ℃,在環(huán)境溫度為20 ℃時(shí),熱泵系統(tǒng)的實(shí)際出風(fēng)溫度最高可以達(dá)到80 ℃,完全能夠滿(mǎn)足稻谷烘干的工藝要求。

總體來(lái)說(shuō),仿真值總是高于試驗(yàn)結(jié)果,但變化趨勢(shì)一致。隨著環(huán)境溫度的增高,熱風(fēng)溫度和COP都隨之增大。制熱量雖有明顯波動(dòng),但也呈增加的趨勢(shì)。誤差分析顯示,在環(huán)境溫度為0 ℃～20 ℃范圍內(nèi),系統(tǒng)輸出熱風(fēng)溫度的平均絕對(duì)誤差MAD=8.16,平均相對(duì)誤差MRE=12.0%;制熱量平均絕對(duì)誤差MAD=12.04,平均相對(duì)誤差MRE=9.5%;COP平均絕對(duì)誤差MAD=0.516,平均相對(duì)誤差MRE=14.02%。導(dǎo)致誤差的原因,可能是因?yàn)橄到y(tǒng)實(shí)際工作中,其工作溫度、轉(zhuǎn)速、排氣容積、壓力等數(shù)值會(huì)不斷發(fā)生變化,但仿真計(jì)算時(shí),采用固定不變的理論值,因此兩者會(huì)產(chǎn)生一定的誤差。另外,在仿真過(guò)程中,忽略了熱泵系統(tǒng)的熱量損失與內(nèi)部消耗,也是導(dǎo)致仿真數(shù)值大于試驗(yàn)數(shù)值的一個(gè)重要因素。但是,試驗(yàn)結(jié)果也顯示出,盡管仿真與實(shí)際值之間存在一定的誤差,但平均相對(duì)誤差均小于15%,處于工程上可以接受的范圍內(nèi)。說(shuō)明所建仿真模型可行,可以用于熱泵系統(tǒng)的模擬計(jì)算和試驗(yàn)研究。

5 結(jié)論

1) 環(huán)境溫度對(duì)熱泵系統(tǒng)的性能有明顯的影響,隨著環(huán)境溫度的增加,熱泵系統(tǒng)的熱風(fēng)溫度、制熱量、COP均呈增加的趨勢(shì)。

2) 在環(huán)境溫度為0 ℃～20 ℃范圍內(nèi),多級(jí)熱泵系統(tǒng)的出風(fēng)溫度可在55 ℃～80 ℃變化,完全能夠滿(mǎn)足糧食烘干的工藝要求。

3) 盡管存在一定的誤差,但可以利用所建立的Simulink模型對(duì)系統(tǒng)的熱風(fēng)溫度、制熱量和COP值進(jìn)行預(yù)測(cè),為熱泵系統(tǒng)的設(shè)計(jì)和試驗(yàn)研究提供了一種有效的技術(shù)手段。