熱力學熵分析法改善雙級噴射制冷系統(tǒng)的研究

許浩，李香，徐昆，劉進陽，謝超許，王博韜, 王南, 王祥

?

熱力學熵分析法改善雙級噴射制冷系統(tǒng)的研究

許浩，李香，徐昆，劉進陽，謝超許，王博韜, 王南, 王祥

（廣西大學 化學化工學院，廣西石化資源加工與過程強化重點實驗室，廣西 南寧 530004）

利用熱力學第一定律和熱力學第二定律熵分析法對雙級噴射制冷系統(tǒng)進行研究，比較了多組制冷劑在不同發(fā)生溫度、蒸發(fā)溫度和冷凝溫度下，系統(tǒng)的噴射系數(shù)，性能系數(shù)及機械性能系數(shù)的變化趨勢。分析了系統(tǒng)各部件不可逆損失的分布情況，同時分析了工作參數(shù)的變化對系統(tǒng)各部分熵產(chǎn)的影響，進而確定使系統(tǒng)性能系數(shù)較佳同時滿足總熵產(chǎn)最小的運行參數(shù)，即系統(tǒng)工作性能最優(yōu)的工況。

雙級噴射式制冷；熵分析法；不可逆性；性能優(yōu)化

傳統(tǒng)制冷循環(huán)基于機械式蒸汽壓縮，結構緊湊、效率較佳而被廣泛應用，但其能耗較大，通過燃燒化石燃料來供電驅(qū)動，會產(chǎn)生大量空氣污染物[1]。隨著社會經(jīng)濟的發(fā)展，噴射制冷系統(tǒng)作為一種環(huán)保友好、有效節(jié)能的新途徑，可回收工業(yè)廢氣余熱等低品位熱源及利用太陽能、地熱能等新能源來驅(qū)動，提高了能源利用率，并減少了電力消耗。噴射器較早地應用于熱源驅(qū)動的制冷系統(tǒng)中，如利用太陽能集熱器提供的熱量加熱發(fā)生器中的制冷劑，進入噴射器內(nèi)不斷提供動力及壓差，使循環(huán)得以維持[2]。

國內(nèi)外對噴射式制冷系統(tǒng)的研究主要在噴射器幾何結構的改進[3]、制冷工質(zhì)的優(yōu)化[4]及運行狀態(tài)的設定[5]上，探究工況變化（壓力、溫度等）對噴射制冷過程噴射系數(shù)及性能系數(shù)COP的影響規(guī)律，從而改進噴射器的設計，使噴射制冷系統(tǒng)性能更為完善[6]。制冷效率較低及對運行條件的高靈敏度限制了噴射器的廣泛應用，但作為新型能源的利用方式，也逐漸被科研及工程人員所重視，發(fā)展前景和應用價值依然十分廣闊。僅從能量數(shù)量上的守恒將性能系數(shù)作為評判噴射制冷系統(tǒng)性能優(yōu)劣的標準，無法對能量品質(zhì)變化及傳遞方向給予評估，為了更全面地考察系統(tǒng)各環(huán)節(jié)的完善性，利用熵分析法對整個循環(huán)各部件的熵產(chǎn)及做功能力損失進行分析，找到不可逆程度最大的環(huán)節(jié)，從而更精準地反映系統(tǒng)的用能情況，指導噴射制冷系統(tǒng)的改善及優(yōu)化。

1 噴射器及兩級噴射制冷系統(tǒng)

1.1 系統(tǒng)熱力學模型

噴射器廣泛應用于石油化工、航空動力、空調(diào)制冷等領域，相比于機械增壓型設備(壓縮機、膨脹機、渦流管等)，降低了膨脹過程的節(jié)流損失，提高了引射流體壓力而不直接消耗機械能[7]。如圖1所示，噴射器由噴嘴、吸收室、混合室及擴壓室組成，屬于非容積型節(jié)流機構[8]，其結構簡單、運行穩(wěn)定、安裝維護方便、適用于兩相流工況，可靠性較高。當蒸發(fā)溫度較低或冷凝溫度較高時，噴射器壓縮比過大，單級噴射制冷系統(tǒng)的性能會出現(xiàn)惡化，采用如圖2所示串聯(lián)增壓式的雙級噴射制冷系統(tǒng)適用性更強，包括兩個能量子循環(huán)1-g2-c2-4-5-7-1、1-g1-c1-3-e2-c2-4-5-7-1（卡諾循環(huán)，為整個系統(tǒng)供能、為升壓提供動力）和制冷子循環(huán)2-e1-c1-3-e2-c2-4-5-6-2（逆向卡諾循環(huán)，完成蒸發(fā)制冷負荷）。本文沒有限定發(fā)生器的熱源形式，可根據(jù)應用場合的不同，利用太陽能、汽車尾氣、工業(yè)廢氣余熱等來加熱傳熱工質(zhì)。

圖1 噴射器結構示意圖

圖2 雙級噴射式制冷循環(huán)示意圖

1.2 系統(tǒng)工作原理

本文所述雙級噴射制冷系統(tǒng)對應的壓-焓圖與溫-熵圖分別如圖3、圖4所示。制冷劑在發(fā)生器中吸熱汽化形成高溫高壓蒸氣(狀態(tài)1)，分別經(jīng)噴射器Ⅰ、Ⅱ的縮放型噴嘴絕熱膨脹，形成壓力降低、流速升至超音速的工作流體，在噴嘴出口附近產(chǎn)生真空攜帶作用以抽吸引射流體；從蒸發(fā)器流出的低溫低壓蒸氣(狀態(tài)2)作為第一級的引射流體(狀態(tài)e1)進入噴射器引射室，與第一級的工作流體(狀態(tài)g1)在混合室充分混合(狀態(tài)c1)，經(jīng)擴壓室進一步減速升壓到一定背壓，最終形成一股速度均勻、壓力居中的流體(狀態(tài)3)；第二級的引射流體(狀態(tài)e2)恰好為噴射器Ⅰ出口的混合流體，與第二級的工作蒸氣(狀態(tài)g2)充分混合(狀態(tài)c2)；從擴壓室出來的混合流體(狀態(tài)4)進入冷凝器(風冷式)凝結放熱；流出的冷凝液(狀態(tài)5)分為兩股，一股通過膨脹閥節(jié)流降壓后(狀態(tài)6)回流至蒸發(fā)器吸熱汽化以蒸發(fā)制冷；另一股則通過循環(huán)泵增壓后(狀態(tài)7)抽送回發(fā)生器中, 被重新加熱成高溫高壓氣態(tài)。

圖3 雙級噴射式制冷系統(tǒng)壓-焓圖

圖4 雙級噴射式制冷循環(huán)溫-熵圖

2 系統(tǒng)能量守恒計算模型

噴射系數(shù)定義為噴射器內(nèi)單位質(zhì)量工作流體所能抽吸引射流體的量。確定了工作流體和引射流體進口狀態(tài)（壓力、溫度）以及噴射器的出口壓力，可計算出系統(tǒng)最大的噴射系數(shù)[9]；在建立系統(tǒng)數(shù)學模型時，假設工質(zhì)在發(fā)生器、蒸發(fā)器、冷凝器的流動損失忽略不計，且管路內(nèi)的流動過程絕熱。

噴射器Ⅰ、Ⅱ的噴射系數(shù)可由下式求得：

兩級噴射器的總噴射系數(shù)：

由噴射器質(zhì)量和能量守恒方程有：

傳熱介質(zhì)通過集熱器時將熱流量提供給系統(tǒng)，發(fā)生器和管路散失到環(huán)境中的熱量忽略不計，則發(fā)生輸送給制冷子循環(huán)的熱量:

式中：g—兩級噴射器工作流體的質(zhì)量流量之和，kg·s-1；

1—工質(zhì)流出發(fā)生器時的比焓值，J·kg-1；

7—工質(zhì)流入發(fā)生器時的比焓值，J·kg-1。

蒸發(fā)器的制冷量：

式中：e1—噴射器Ⅰ中引射流體的質(zhì)量流量，kg·s-1；

2—工質(zhì)流出蒸發(fā)器時的比焓值，J·kg-1；

6—工質(zhì)流入蒸發(fā)器時的比焓值，J·kg-1。

冷凝器的散熱量：

式中：c2—噴射器Ⅱ中混合流體的質(zhì)量流量，kg·s-1；

4—工質(zhì)流入冷凝器時的比焓值，J·kg-1；

5—工質(zhì)流出冷凝器時的比焓值，J·kg-1。

循環(huán)泵的耗功量：

假定循環(huán)工質(zhì)在膨脹閥節(jié)流前后的焓值不變，即：

噴射制冷子循環(huán)系統(tǒng)的性能系數(shù)：

由于噴射制冷系統(tǒng)由太陽能、廢氣余熱等免費提供熱源，通常也可用機械性能系數(shù)來評價[8]：

整個雙級噴射制冷系統(tǒng)的總性能系數(shù)：

現(xiàn)采用R600a、R290、R141b、R134a、R123、R152a作為制冷劑，工質(zhì)物性在各狀態(tài)點下的熱力學參數(shù)可通過調(diào)用NIST Refprop 8.0數(shù)據(jù)庫求得，運用Matlab軟件編制程序計算，探討不同制冷劑溫度變化對系統(tǒng)性能的影響，分別如圖4-圖12所示。設定發(fā)生溫度Tg、蒸發(fā)溫度Te、冷凝溫度Tc的變化范圍分別為71~85 ℃，-15~10 ℃，45~65 ℃，發(fā)生器和蒸發(fā)器中的過熱度分別取為5 ℃和3 ℃。發(fā)生壓力、蒸發(fā)壓力、冷凝壓力為各制冷工質(zhì)對應溫度下的飽和壓力[10]。設定噴射制冷系統(tǒng)的制冷量為100 kW，循環(huán)泵的絕熱效率為85%。

圖5 不同制冷劑的性能系數(shù)隨發(fā)生溫度變化的影響

分析圖4-圖12曲線特性可得出以下規(guī)律：

圖6 不同制冷劑的機械性能系數(shù)隨發(fā)生溫度變化的影響

圖7 不同制冷劑的噴射系數(shù)隨蒸發(fā)溫度變化的影響

圖8 不同制冷劑的性能系數(shù)隨蒸發(fā)溫度變化的影響

增加發(fā)生溫度g，工作流體通過縮放型噴嘴時能提供更多由焓值轉(zhuǎn)化而來的動能，相應抽吸引射流體的能力增強。但系統(tǒng)的機械性能系數(shù)COPm隨發(fā)生溫度的升高有所降低，如圖6所示，因為在保持發(fā)生器吸熱量一定時，工質(zhì)泵耗功隨發(fā)生溫度的升高而增加了。

圖9 不同制冷劑的機械性能系數(shù)隨蒸發(fā)溫度變化的影響

圖10 不同制冷劑的噴射系數(shù)隨冷凝溫度變化的影響

圖11 不同制冷劑的性能系數(shù)隨冷凝溫度變化的影響

圖12 不同制冷劑的機械性能系數(shù)隨冷凝溫度變化的影響

3 雙級噴射制冷系統(tǒng)的熵分析

以上分析基于熱力學第一定律的能量守恒角度，但忽略了能量品質(zhì)因素。雙級噴射制冷系統(tǒng)包括節(jié)流、熱交換等絕熱過程，從能量平衡角度來考慮這些過程并無能量損失，但節(jié)流前后工質(zhì)的壓力降低了，能量品質(zhì)隨之下降，因此對系統(tǒng)內(nèi)部損失的分析還不夠準確。熵由過程的不可逆性產(chǎn)生，用來判定熱現(xiàn)象過程進行的方向、限度和條件，其大小可作為熱力過程不可逆性的量度，其變化指明了過程進行的方向。熵增表明不可用功的增加，也可用來衡量過程的不可逆程度及做功能力的損失[11]。制冷技術的關鍵在于從被冷卻物體中抽取熱量，相當于抽取熵，制冷系統(tǒng)的熵產(chǎn)與性能損失有關，工質(zhì)流過各部件引起的熵產(chǎn)等于系統(tǒng)內(nèi)工質(zhì)的熵增及進出系統(tǒng)熱量熵變之和[12]。通過熵分析可以優(yōu)化噴射器結構，提高設備傳能及傳質(zhì)效率，確定出最佳工作參數(shù)，從而更全面地改善制冷系統(tǒng)的性能。

不可逆過程開口系熵方程的一般表達式[13]為：

積分可得：

—系統(tǒng)與外界的換熱量，kJ；

—熱源溫度，K；

—制冷工質(zhì)的比熵，J/(kg·K)；

g—系統(tǒng)的熵產(chǎn)，J/(kg·K)；

—質(zhì)量流量，kg·s-1；

其中下標：—第個部件；—第個部件；

in—流進；out—流出；

cv—控制體體積。

將系統(tǒng)內(nèi)的流動視為穩(wěn)態(tài)，忽略各部件散失到環(huán)境中的熱量，由式(19)可得熵產(chǎn)的計算公式：

發(fā)生器中的熵產(chǎn)：

式中：ac—集熱器中傳熱工質(zhì)的質(zhì)量流量，kg·s-1；

in—傳熱工質(zhì)在發(fā)生器進口的比熵，J/(kg·K) ；

out—傳熱工質(zhì)在發(fā)生器出口的比熵，J/(kg·K)。

噴射器的熵產(chǎn)：

式中：r—室內(nèi)溫度，取為20 ℃。

式中：0—環(huán)境溫度，取為30 ℃。

由Gouy-Stodla方程式[14]可得到各過程做功能力的損失：

系統(tǒng)總做功能力損失[11]：

運用Matlab軟件編制程序計算，由上文分析可知以R152a作為制冷劑時的整體性能較佳，由于篇幅有限，現(xiàn)只選取R152a在某一工況下（發(fā)生溫度g=75 ℃, 蒸發(fā)溫度e=-10 ℃, 冷凝溫度c=60 ℃）為例進行系統(tǒng)各部件的熵產(chǎn)分析，如表1所示。

表1 雙級噴射制冷系統(tǒng)中各部件熵產(chǎn)及做功能力損失分布情況

由表可知循環(huán)中各部件的最大做功能力損失發(fā)生在兩級噴射器上，占總熵產(chǎn)的44.765%，由于工作流體在噴嘴出口附近的超音速流動及混合流體在擴壓室背壓影響下會產(chǎn)生激波，蒸氣穿過激波時受到激波阻力，造成較大的不可逆損失；此外噴射系數(shù)較低時，高速工作流體產(chǎn)生真空度卷吸引射流體的能力下降，工作流體在噴射器內(nèi)的摩擦增加，熵增變大，流體的能量損失也就增加[15]。其次發(fā)生在冷凝器與發(fā)生器，由于非理想絕熱與不等溫傳熱引起，然后是蒸發(fā)器，它們共占總熵產(chǎn)和總做功能力損失的93.977%，提高這些關鍵部件的換熱性能可有效降低整個雙級噴射制冷系統(tǒng)的不可逆損失，進而優(yōu)化系統(tǒng)性能。

下面對該系統(tǒng)各部件在不同工況下相應的熵產(chǎn)變化進行了分析，如圖13-18所示。

圖13 各部件的熵產(chǎn)隨發(fā)生溫度變化的影響

分析圖13-圖18曲線特性可得出以下規(guī)律：

（1）由圖14可知，系統(tǒng)的總熵產(chǎn)隨發(fā)生溫度g的升高而增大，且在等溫差梯度的蒸發(fā)溫度下隨發(fā)生溫度升高而增大的趨勢不變；由圖16可知，系統(tǒng)總熵產(chǎn)隨蒸發(fā)溫度e的升高而減小，且在等溫差梯度下發(fā)生溫度越高，總熵產(chǎn)下降越快；由圖18可知，總熵產(chǎn)隨冷凝溫度c的升高而增大，且在等溫差梯度下發(fā)生溫度越低，總熵產(chǎn)增加趨勢越緩慢。

圖14 不同蒸發(fā)溫度下總熵產(chǎn)與發(fā)生溫度的關系

圖15 各部件的熵產(chǎn)隨蒸發(fā)溫度變化的影響

圖16 不同發(fā)生溫度下總熵產(chǎn)與蒸發(fā)溫度的關系

（2）由圖13可知，隨著發(fā)生溫度g增大，兩級噴射器和發(fā)生器的熵產(chǎn)持續(xù)增加，冷凝器的熵產(chǎn)有所下降，蒸發(fā)器的熵產(chǎn)幾乎不變。提高發(fā)生溫度可使發(fā)生壓力增加，工作流體流速增快，卷吸引射流體的能力增強，噴射系數(shù)增大。但由于噴射器結構固定，當工況變化時，只能通過改變工作流體和引射流體進入噴射器的壓力，進而影響兩股流體在噴射器內(nèi)的能量和動量交換，當溫度變化達極限值時，繼續(xù)增加反而會使系統(tǒng)性能減小。而發(fā)生器熵產(chǎn)增加則因為其傳熱損失增大，效率下降。

圖17 各部件的熵產(chǎn)隨冷凝溫度變化的影響

圖18 不同發(fā)生溫度下總熵產(chǎn)與冷凝溫度的關系

（3）由圖15可知，提高蒸發(fā)溫度e，兩級噴射器的熵產(chǎn)下降趨勢最明顯，總熵產(chǎn)降低，因為與室溫溫差減小，不等溫傳熱造成的不可逆損失減少。

（4）由圖17可知，增大冷凝溫度c，兩級噴射器熵產(chǎn)增加，發(fā)生器熵產(chǎn)減小，冷凝器熵產(chǎn)增加，系統(tǒng)的總熵產(chǎn)不斷增加，因為噴射系數(shù)和系統(tǒng)的性能系數(shù)隨冷凝溫度的增大而下降的幅度較為明顯。

4 結論

根據(jù)實驗條件的不同選擇合適的制冷劑可以優(yōu)化系統(tǒng)性能。通過熵分析可知，噴射器熵產(chǎn)最大，是影響制冷系統(tǒng)動態(tài)性能的主要因素。在實際系統(tǒng)設計時，可對制冷系統(tǒng)采用雙級或多級噴射器隨工況變化進行調(diào)節(jié)，使系統(tǒng)運行性能保持最佳。

雙級噴射制冷系統(tǒng)的噴射系數(shù)、性能系數(shù)隨著發(fā)生溫度和蒸發(fā)溫度的增大而不斷增加，卻隨著冷凝溫度的增大而顯著降低；系統(tǒng)的總熵產(chǎn)隨著發(fā)生溫度和冷凝溫度的增大而不斷增加，卻隨著蒸發(fā)溫度的增大而明顯降低。為了提高系統(tǒng)的整體性能，即在滿足系統(tǒng)運行效率較佳的同時使系統(tǒng)總熵產(chǎn)盡可能小，可適當提高系統(tǒng)的蒸發(fā)溫度，降低冷凝溫度；并在其他條件一定時，選擇適當?shù)陌l(fā)生溫度。

［1］Omidvar A, Ghazikhani M. Entropy analysis of a solar-driven variable geometry ejector using computational fluid dynamics [J]. Energy Conversion and Management, 2016, 119: 435-443.

［2］張玲玲,陶樂仁. 噴射式制冷的發(fā)展研究現(xiàn)狀[J].制冷與空調(diào), 2012, 26(5): 504 -510.

［3］Zhu Y, Cai W, Wen C, et al. Numerical investigation of geometry parameters for design of high performance ejectors[J]. Applied Thermal Engineering, 2009, 29(5–6): 898-905.

［4］王菲, 沈勝強.不同制冷劑噴射制冷性能計算分析[J]. 化工學報,2010, 61(2): 275-280.

［5］Rogdakis E D, Alexis G K. Design and parametric investigation of an ejector in an air-conditioning system[J]. Applied Thermal Engineering, 2000, 20(2): 213-226.

［6］李慶普，陶樂仁，劉效德. 太陽能噴射式制冷系統(tǒng)性能分析[J]. 制冷技術, 2015, 35(6): 67-71.

［7］Huang B J, Chang J M. Empirical Correlation for Ejector Design[J]. International Journal of Refrigeration, 1999, 22（5）: 379-388.

［8］丁寶存, 楊洪源, 張琨.噴射器理論發(fā)展的論述[J]. 區(qū)域供熱, 2013 (1): 65-67, 74.

［9］索科洛夫 E R，津格爾 H M. 噴射器[M].北京: 科學出版社，1977.

［10］劉娜，楊啟容，吳榮華. 太陽能噴射式制冷系統(tǒng)熱力學分析[J]. 石油化工設備, 2014, 43(5): 23-28.

［11］季一新，陶樂仁，王金鋒. 太陽能噴射式制冷循環(huán)的熱力學熵分

析[J]. 制冷技術, 2008, 36(12): 52-55, 65.

［12］邱金友, 余曉明, 倪鉒棟,等.一種變蒸發(fā)溫度冷水機組的設計及節(jié)能熵分析[J]. 制冷技術,2013, 41(4): 65-68.

［13］嚴濟慈. 熱力學第一定律和第二定律[M]. 北京：人民教育出版社，1978.

［14］Jones J B, Dugan R E. Engineering thermodynamics[M]. New Jersey Prentice Hall Inc,1996.

［15］趙良舉,曾丹芬,袁鵬,等. 汽液兩相混合物的加速與激波的熱力學分析[J]. 工程熱物理學報, 2001, 22(3): 284-286.

Research on Improving the Two-stage Ejector Refrigeration System Based on the Thermodynamic Entropy Analysis Method

,,,,,,,

（Key Laboratory of Guangxi Petrochemical Resource Processing and Process Intensification, School of Chemistry and Chemical Engineering, Guangxi University, Guangxi Nanning 530004, China）

The first law of thermodynamics and entropy analysis in the second law of thermodynamics were used to study the two-stage ejector refrigeration system. The variation trends of the entrainment ratio, the performance coefficient and the mechanical performance coefficient of the system using different several kinds of refrigerants were studied and compared at different generation temperature, evaporation temperature and condensation temperature. The irreversible loss distribution of each component in the system was analyzed, the influence of working parameters on each part of the system entropy production was also analyzed, then the best operating parameters were determined to ensure better system performance coefficient and satisfy the total entropy production minimum simultaneously.

Two-stage ejector refrigeration; Entropy analytic method; Irreversibility; Optimum performance

TQ 013.1

A

1671-0460（2017）10-2066-07

2017-02-16

許浩（1991-），男，湖北仙桃人，碩士研究生，廣西大學化學化工學院化工過程機械專業(yè)，研究方向：化工設備性能分析及優(yōu)化設計。 E-mail：249453257@qq.com。

李香（1987-），女，助教，碩士，研究方向：化學工程與化工設備的傳遞過程。E-mail：gxlixiang@163.com。