楊婧燁,陸冰清,陳江平
(上海交通大學(xué)制冷及低溫工程研究所,上海 200240)
能源緊缺和氣候環(huán)境的惡化使得能源的梯級(jí)有效利用受到廣泛關(guān)注。發(fā)動(dòng)機(jī)廢熱是中低溫余收的主要來(lái)源[1-3],燃油燃燒過(guò)程中60%~70%的能量會(huì)隨著尾氣排放至大氣造成資源浪費(fèi)和環(huán)境污染。發(fā)動(dòng)機(jī)廢熱主要分布在3個(gè)部分:發(fā)動(dòng)機(jī)高溫?zé)煔庥酂?、發(fā)動(dòng)機(jī)缸套水廢熱和發(fā)動(dòng)機(jī)冷卻液。小型有機(jī)朗肯循環(huán)作為應(yīng)用于中低溫廢熱回收的主要技術(shù)手段之一,其工質(zhì)選擇一直是國(guó)內(nèi)外汽車行業(yè)研究的熱點(diǎn)。
Shale[4]評(píng)估了多種工質(zhì)在中低溫余熱回收中的應(yīng)用,指出HFC型制冷劑如R134a和R245fa是適用于有機(jī)朗肯循環(huán)的低臨界溫度工質(zhì);Quolin[5]明確指出R245fa是適用于商用中低溫廢熱回收有機(jī)朗肯循環(huán)中的常見工質(zhì);Muhammad[6]建立了1kW的小型有機(jī)朗肯循環(huán),應(yīng)用R123作為工質(zhì)獲得4.66%的熱效率;Li[7]基于帶回?zé)峤Y(jié)構(gòu)的有機(jī)朗肯循環(huán)研究應(yīng)用R123作為工質(zhì)獲得約6kW輸出功。
隨著全球氣候變暖的影響,制冷劑選擇標(biāo)準(zhǔn)[8-9]不再局限于ODP(ozone depletion potential)為零,2017年 MAC Derivative提出未來(lái)GWP(global warming potential)高于150的制冷劑均應(yīng)被禁止使用。尋求環(huán)境友好型的清潔制冷劑的替代工作迫在眉睫。應(yīng)用于有機(jī)朗肯循環(huán)中的高GWP制冷劑也將面臨淘汰危機(jī)。一批新HCFO型低GWP制冷劑[10-12]因其物性與R134a和R245fa相近有望替代成為低溫有機(jī)朗肯循環(huán)工質(zhì),其中包括R1234yf,R1234ze,R1233zde和 R1336mzzZ。
為了研究R1233zde作為R245fa替代物應(yīng)用于低溫有機(jī)朗肯循環(huán)中的可行性,本文中針對(duì)有機(jī)朗肯循環(huán)建立了穩(wěn)態(tài)的系統(tǒng)仿真模型,包括4部件的獨(dú)立模型?;诜抡婺P头治霰容^了R12333zde和R245fa在給定工況下,各系統(tǒng)性能參數(shù)包括蒸發(fā)溫度、泵耗功和循環(huán)熱效率,說(shuō)明 R1233zde替代R245fa的可行性。
表1 工質(zhì)物性參數(shù)對(duì)比
圖1 T-S飽和曲線對(duì)比
圖2 飽和壓力對(duì)比
有機(jī)朗肯循環(huán)工質(zhì)的選擇遵循基本的制冷劑物性要求,無(wú)毒無(wú)腐蝕性且不含ODP。除此之外,根據(jù)有機(jī)朗肯的循環(huán)特性,工質(zhì)在性能參數(shù)上還必須滿足低臨界壓力、高臨界溫度的要求從而產(chǎn)生高的蒸發(fā)溫度以獲得更高的循環(huán)熱效率。常見應(yīng)用于低溫余熱回收的R245fa工質(zhì)雖然不含ODP,但其GWP高達(dá)1 030,在未來(lái)更加嚴(yán)格的排放法規(guī)的標(biāo)準(zhǔn)下面臨禁用的危機(jī)。R1233zde因其物性與R245fa十分相近被提出有望替代成為新余熱回收有機(jī)工質(zhì)。表1分別是R1233zde和R245fa的主要物性參數(shù)對(duì)比,圖1[13]為兩種工質(zhì)的飽和T-S曲線對(duì)比,圖2為兩種工質(zhì)相同飽和溫度下對(duì)應(yīng)的飽和壓力值的比較??梢钥闯?雖然兩者物性十分相似,但R1233zde的臨界壓力稍低于R245fa且臨界溫度略高,理論分析認(rèn)為R1233zde在相同的最大工作壓力條件下可以獲得更高的蒸發(fā)溫度,循環(huán)熱效率也會(huì)相應(yīng)提高。
系統(tǒng)仿真的目的是為了比較相同工況下R1233zde與R245fa主要性能參數(shù)表現(xiàn)。因此,針對(duì)有機(jī)朗肯循環(huán),本文中搭建了穩(wěn)態(tài)的半經(jīng)驗(yàn)系統(tǒng)模型,為了簡(jiǎn)化計(jì)算做出了必要的假設(shè)從而無(wú)法作為完全預(yù)測(cè)型的系統(tǒng)模型。圖3為靜態(tài)仿真模型的算法流程圖,算法主要包括3個(gè)收斂條件,對(duì)蒸發(fā)器、冷凝器和膨脹機(jī)進(jìn)行狀態(tài)參數(shù)的迭代計(jì)算。模型相關(guān)假設(shè)如下:
(1)工質(zhì)泵與膨脹機(jī)循環(huán)效率被視作常數(shù)計(jì)算,分別為0.75和0.8;
(2)忽略管路及部件的壓降;
(3)忽略部件與外部的傳熱損失等;
(4)過(guò)冷度作為模型的輸入條件之一。
圖3 算法流程圖
(1)熱交換器
熱交換器的獨(dú)立部件模型采用穩(wěn)態(tài)集中參數(shù)分布法。相區(qū)分布見圖4,圖中Tr和hr分別代表了制冷劑側(cè)溫度和焓值;下角標(biāo)l,tp和v分別表示液態(tài)、兩相和氣態(tài);下角標(biāo)sf表示第二側(cè)流體。
傳熱區(qū)域分為過(guò)冷、兩相和過(guò)熱區(qū)域。蒸發(fā)器與冷凝器的計(jì)算主要區(qū)別在于相區(qū)換熱關(guān)聯(lián)式的選擇,單相區(qū)均使用Muley(1998)。傳熱計(jì)算采用ENTU法計(jì)算換熱量和各相區(qū)換熱面積。表2列出了換熱器模型中所應(yīng)用的換熱關(guān)聯(lián)式。
圖4 相區(qū)分布
表2 各相區(qū)傳熱關(guān)聯(lián)式
(2)工質(zhì)泵及膨脹機(jī)
工質(zhì)泵采用的是Ylya Teleman[14]提出的簡(jiǎn)化模型,工質(zhì)泵的流量和功耗為
式中:vin為泵進(jìn)口的比體積;ηmp為泵的效率。
膨脹機(jī)的理想出口溫度為
式中k為該工質(zhì)的絕熱系數(shù)。
實(shí)際出口焓值由膨脹機(jī)的等熵效率修正:
式中:ηmech為膨脹劑的機(jī)械效率;Vs_cp為渦旋壓縮的工作容積;rv為設(shè)計(jì)體積比。
系統(tǒng)整體的循環(huán)效率為
模型參數(shù)的校核基于全工況實(shí)驗(yàn)數(shù)據(jù)所得。表3和表4分別為系統(tǒng)模型的輸入?yún)?shù)和傳熱關(guān)聯(lián)式的修正系數(shù)。圖5為仿真模型與實(shí)驗(yàn)流量值的對(duì)比,精度在3%左右。
表3 模型輸入?yún)?shù)
表4 傳熱關(guān)聯(lián)式修正參數(shù)
圖5 流量精度擬合
基于修正后的仿真模型對(duì)R1233zde和R245fa進(jìn)行分析計(jì)算。保證外部工況條件相同,控制熱源流量與冷源側(cè)不變,僅改變熱源溫度,從80℃連續(xù)變化至120℃,蒸發(fā)器出口過(guò)熱度控制在2℃。通過(guò)迭代計(jì)算得到兩種工質(zhì)的狀態(tài)點(diǎn)分布及性能參數(shù)指標(biāo)大小。圖6為工質(zhì)泵耗功的比較,隨著熱源溫度的上升,工質(zhì)泵耗功也逐漸增加,在相同的外部工況下,R1233zde系統(tǒng)所消耗的泵功明顯少于R245fa系統(tǒng),減少幅度在15.05%~17.02%之間。原因是在相同的運(yùn)行溫度下R1233zde對(duì)應(yīng)的工作壓力較低。工質(zhì)泵作為系統(tǒng)動(dòng)力的主要來(lái)源,其耗功大小也影響著系統(tǒng)性能的優(yōu)劣,減少泵耗功是系統(tǒng)優(yōu)化的方向之一。
圖6 工質(zhì)泵耗功對(duì)比
圖7 循環(huán)熱效率對(duì)比
圖7為各個(gè)熱源工況下R1233zde和R245fa系統(tǒng)的循環(huán)熱效率的對(duì)比。結(jié)果表明,R1233zde的循環(huán)效率略高于 R245fa。對(duì)應(yīng) 120℃熱源工況下,R1233zde的循環(huán)效率優(yōu)于R245fa約7.03%。原因是,相同工況下,R1233zde消耗的泵功相對(duì)更少,在最大工作壓力下,R1233zde可達(dá)到的蒸發(fā)溫度更高,相應(yīng)的循環(huán)熱效率也就越高。
圖8 100℃熱源下T-S狀態(tài)點(diǎn)分布對(duì)比
圖8列舉了100℃熱源工況下,R1233zde與R245fa在T-S曲線上各個(gè)狀態(tài)參數(shù)點(diǎn)的分布。由圖可知,在給定相同的工況條件下,R1233zde可以達(dá)到相對(duì)更高的蒸發(fā)溫度,在冷凝狀態(tài)相近的情況下,更高的蒸發(fā)溫度帶來(lái)更優(yōu)的熱效率。在Bala V.Dalta[15]的實(shí)驗(yàn)分析中,R1233zde的熱效率比R245fa高出了8.7%。綜上所述,R1233zde被認(rèn)為是性能略優(yōu)于R245fa的低GWP制冷劑,完全有能力替代R245fa成為中低溫有機(jī)朗肯余熱回收的主要有機(jī)工質(zhì)。
基于仿真結(jié)果認(rèn)為,低GWP的R1233zde在未來(lái)有望成為替代R245fa應(yīng)用于中低溫有機(jī)朗肯余熱回收技術(shù)中的有機(jī)工質(zhì)。
(1)R1233zde臨界壓力較R245fa低,臨界溫度相對(duì)高,可以在最大運(yùn)行壓力下比R245fa達(dá)到更高的蒸發(fā)溫度從而提升熱效率,120℃工況下大約提升7.03%。
(2)應(yīng)用R1233zde可以有效減少工質(zhì)泵耗功,優(yōu)化系統(tǒng)整體性能表現(xiàn)。
(3)因R1233zde與R245fa的物性十分相似,在零部件的設(shè)計(jì)上可以沿用R245fa系統(tǒng)部件而不存在硬件需求上的難度。
(4)R1233zde的GWP值僅為R245fa的1/150,ODP接近于0對(duì)環(huán)境無(wú)損害,安全等級(jí)A1且不可燃,符合新制冷劑替代的要求。