蒸發(fā)器與熱力膨脹閥回路穩(wěn)定性研究

葉奇昉，劉杰，陳江平

(1-上海汽車集團(tuán)技術(shù)中心，上海 210804；2-上海交通大學(xué)，上海 200240）

符號(hào)說明：

C——比熱 L——長度 α——對(duì)流換熱系數(shù) t--——時(shí)間 D——直徑——質(zhì)量流量

Cd——流量系數(shù) m——質(zhì)量 ρ——密度 T——溫度

下標(biāo)

a——環(huán)境 i——進(jìn)口/內(nèi)表面 r——制冷劑 b——溫包 g——?dú)庀?sh——過熱

c——冷凝器/膜片 l——液相 tp——兩相 e——蒸發(fā)器 o——出口/外表面 w——壁面

0 引言

熱力膨脹閥作為節(jié)流機(jī)構(gòu)在制冷系統(tǒng)中有廣泛應(yīng)用，2012年國內(nèi)汽車空調(diào)系統(tǒng)上使用熱力膨脹閥就超過一千萬只。然而，熱力膨脹閥在使用過程中經(jīng)常要面對(duì)不同的動(dòng)力學(xué)問題。最典型的如：熱力膨脹閥系統(tǒng)工作不穩(wěn)定，具體表現(xiàn)為：膨脹閥頻繁開啟，制冷系統(tǒng)供液量出現(xiàn)周期性波動(dòng)，系統(tǒng)各處的溫度壓力產(chǎn)生振蕩。對(duì)這一現(xiàn)象，國內(nèi)外的研究人員進(jìn)行了大量的研究。Huelle[1]在對(duì)蒸發(fā)器完全蒸干點(diǎn)隨機(jī)振蕩研究結(jié)果的基礎(chǔ)上，對(duì)熱力膨脹閥控制的蒸發(fā)器的穩(wěn)定性進(jìn)行了研究，發(fā)現(xiàn)蒸發(fā)器振蕩現(xiàn)象容易發(fā)生在過熱度較低時(shí)。他認(rèn)為蒸發(fā)器最小穩(wěn)定過熱度（蒸發(fā)器由穩(wěn)定工作到產(chǎn)生振蕩時(shí)的最小過熱度）是蒸發(fā)器本身的一個(gè)特性，且與蒸發(fā)器的熱負(fù)荷有關(guān)。Broersen[2]建立了蒸發(fā)器-熱力膨脹閥系統(tǒng)的傳遞函數(shù)模型，利用Nyquist穩(wěn)定性判據(jù)詳細(xì)討論了蒸發(fā)器管壁時(shí)間常數(shù)對(duì)熱力膨脹閥控制的蒸發(fā)器過熱度系統(tǒng)的穩(wěn)定性裕度的影響。Ibrahim和Ghaleb[3]修正了Broersen的模型，研究制冷系統(tǒng)在不同載荷下的穩(wěn)定性，研究認(rèn)為當(dāng)系統(tǒng)的不穩(wěn)定幅度保持不變時(shí)，存在一個(gè)確定的臨界時(shí)間常數(shù)，當(dāng)實(shí)際的時(shí)間常數(shù)小于臨界常數(shù)時(shí)，系統(tǒng)是穩(wěn)定的。Ibrahim[4]研究蒸發(fā)器的水流量和入口溫度變化對(duì)制冷系統(tǒng)穩(wěn)定性的影響。Tahat等人[5]認(rèn)為蒸發(fā)器時(shí)間常數(shù)和換熱系數(shù)互相影響且對(duì)系統(tǒng)的穩(wěn)定性影響很大，其中換熱系數(shù)的提高對(duì)系統(tǒng)穩(wěn)定性更有利。

已有針對(duì)熱力膨脹閥系統(tǒng)穩(wěn)定性的研究大多著眼于系統(tǒng)工況對(duì)系統(tǒng)穩(wěn)定性的影響。而針對(duì)熱力膨脹閥自身設(shè)計(jì)參數(shù)對(duì)穩(wěn)定性影響的研究很少。熱力膨脹閥溫包充注作為熱力膨脹閥最重要的設(shè)計(jì)參數(shù)之一對(duì)制冷系統(tǒng)性能有重要影響。當(dāng)前針對(duì)不同制冷系統(tǒng)開發(fā)的熱力膨脹閥溫包充注方式越來越多，而不同溫包充注的熱力膨脹閥系統(tǒng)穩(wěn)定性相差較大，溫包充注對(duì)系統(tǒng)穩(wěn)定性影響的相關(guān)研究也越來越凸現(xiàn)出其重要性。

本文針對(duì)蒸發(fā)器-蒸發(fā)器循環(huán)建立系統(tǒng)熱動(dòng)力學(xué)非線性模型，采用Nyquist穩(wěn)定性判據(jù)研究系統(tǒng)的穩(wěn)定性，獲得系統(tǒng)的不穩(wěn)定工作區(qū)域，分析了熱力膨脹閥溫包的不同充注方式、充注參數(shù)對(duì)系統(tǒng)穩(wěn)定性的影響。

1 蒸發(fā)器-熱力膨脹閥回路熱動(dòng)力學(xué)模型

1.1 蒸發(fā)器模型

目前換熱器模型采用的建模方法主要有集總參數(shù)法、分布參數(shù)法和移動(dòng)邊界法。移動(dòng)邊界法是目前在研究制冷系統(tǒng)控制中采用的最普遍的方法。它將熱交換器按照制冷劑所處的狀態(tài)分為兩相區(qū)和蒸汽區(qū)，并對(duì)每一個(gè)區(qū)域采用集總參數(shù)模型(如圖1所示)。制冷劑在制冷系統(tǒng)換熱器中呈兩相流動(dòng)，應(yīng)用質(zhì)量守恒方程、能量守恒方程可得到一組復(fù)雜的非線性偏微分方程表示式，在管長方向(界面移動(dòng))對(duì)微分方程組進(jìn)行積分后，可得到用狀態(tài)變量表示的狀態(tài)方程。

圖1 蒸發(fā)器模型示意圖

兩相區(qū)中對(duì)方程(1)、(2)積分，可得：

通過引入平均空泡系數(shù)可得：

式中，空泡系數(shù)采用滑移比為1的均相流模型：

因此方程(3)、(4)可寫為：

在過熱區(qū)，對(duì)方程(8)、(9)積分，可得：

兩相區(qū)與過熱區(qū)蒸發(fā)器壁面能量方程可分別寫為：

蒸發(fā)器換熱系數(shù)隨著管內(nèi)制冷劑的流速狀態(tài)的變化而變化并對(duì)制冷系統(tǒng)產(chǎn)生影響蒸發(fā)器制冷劑側(cè)和空氣側(cè)的換熱系數(shù)分別采用文獻(xiàn)[6-9]提出的公式進(jìn)行計(jì)算。

1.2 熱力膨脹閥模型

圖2為熱力膨脹閥示意圖，膨脹閥開度由溫包壓力Fbulb、蒸發(fā)壓力Fe、彈簧力Fspring共同決定。溫包壓力Fbulb由蒸發(fā)器出口溫度決定，F(xiàn)e由蒸發(fā)壓力決定，F(xiàn)spring由彈簧彈性系數(shù)決定。

圖2 熱力膨脹閥示意圖

當(dāng)熱力膨脹閥開度變化時(shí)，閥芯運(yùn)動(dòng)方程可表示為：

式中：

m——膨脹閥運(yùn)動(dòng)部件的質(zhì)量；

x——膨脹閥開度；

x0——彈簧預(yù)壓縮量；

Pb——溫包內(nèi)壓力；

k——彈簧彈性系數(shù)；

c——運(yùn)動(dòng)阻尼系數(shù)。

熱力膨脹閥的制冷劑流量可表示為：

式中：

Cd——流量系數(shù)，由實(shí)驗(yàn)獲得[10]；

ρc——膨脹閥進(jìn)口制冷劑密度；

A(x) ——膨脹閥流通面積。

當(dāng)熱力膨脹閥工作時(shí)，溫包內(nèi)工質(zhì)進(jìn)行強(qiáng)迫對(duì)流換熱，其熱交換過程可以表示為：

當(dāng)熱力膨脹閥工作時(shí)，溫包膜片的換熱過程可表示為：

式(8)～式(17)為蒸發(fā)器-熱力膨脹閥系統(tǒng)的非線性數(shù)學(xué)模型，用Nyquist穩(wěn)定性判據(jù)即可判斷系統(tǒng)在不同的條件下工作狀態(tài)是否穩(wěn)定。

2 結(jié)果分析

圖3給出了平行充注的熱力膨脹閥在不同工作條件下的工作狀態(tài)。由圖可知，當(dāng)蒸發(fā)溫度大于277 K后，系統(tǒng)存在不穩(wěn)定工作區(qū)域。系統(tǒng)必須存在一定的過熱度，否則將工作不穩(wěn)定。隨著蒸發(fā)溫度的升高，所需的過熱度逐漸增加。

圖3 平行充注的熱力膨脹閥系統(tǒng)穩(wěn)定性

當(dāng)熱力膨脹閥溫包內(nèi)采取交叉充注時(shí)，系統(tǒng)穩(wěn)定性如圖4。由圖可知，采用交叉充注的熱力膨脹閥存在“交叉點(diǎn)”，交叉點(diǎn)壓力大約為0.347 MPa。當(dāng)蒸發(fā)壓力低于“交叉點(diǎn)”壓力時(shí)，溫包內(nèi)工質(zhì)壓力高于相同溫度下的蒸發(fā)壓力，即使不存在過熱度膨脹閥也存在微小開啟度；當(dāng)蒸發(fā)壓力高于“交叉點(diǎn)”壓力時(shí)溫包內(nèi)工質(zhì)壓力低于相同溫度下的蒸發(fā)壓力，蒸發(fā)器出口必須存在一定過熱度膨脹閥才能開啟。蒸發(fā)壓力越高，所需的開啟過熱度也越高。

圖4 交叉充注的熱力膨脹閥系統(tǒng)穩(wěn)定性

交叉充注的熱力膨脹閥也存在不穩(wěn)定工作區(qū)。當(dāng)蒸發(fā)壓力低于“交叉點(diǎn)”壓力時(shí)，熱力膨脹閥不需要過熱度即可開啟，有一定質(zhì)量的制冷劑通過膨脹閥，有利于系統(tǒng)的穩(wěn)定工作。不穩(wěn)定工作區(qū)主要存在于蒸發(fā)壓力大于交叉點(diǎn)壓力的情況下。交叉充注的熱力膨脹閥系統(tǒng)不穩(wěn)定工作區(qū)域遠(yuǎn)比平行充注的熱力膨脹閥系統(tǒng)要小。

改變溫包內(nèi)充注工質(zhì)，使得“交叉點(diǎn)”壓力升高到0.377 MPa，系統(tǒng)工作狀態(tài)如圖5所示。對(duì)比圖4可知當(dāng)“交叉點(diǎn)”壓力升高，系統(tǒng)不穩(wěn)定工作區(qū)域減小。對(duì)于交叉充注的熱力膨脹閥制冷系統(tǒng)“交叉點(diǎn)”壓力越高，對(duì)系統(tǒng)穩(wěn)定工作越有利。但是，提高“交叉點(diǎn)”壓力將增大蒸發(fā)器出口兩相的可能性，容易造成壓縮機(jī)液擊。

改變溫包內(nèi)充注工質(zhì)，使得溫包內(nèi)工質(zhì)的壓力-溫度曲線斜率降低，系統(tǒng)工作狀態(tài)如圖6所示。對(duì)比圖4可知：當(dāng)壓力-溫度曲線斜率降低，系統(tǒng)所需的膨脹閥開啟過熱度增加，而系統(tǒng)不存在不穩(wěn)定工作區(qū)域。對(duì)于交叉充注的熱力膨脹閥制冷系統(tǒng)，壓力-溫度曲線斜率越低，對(duì)系統(tǒng)穩(wěn)定工作越有利。

圖5 提高交叉點(diǎn)壓力的熱力膨脹閥系統(tǒng)穩(wěn)定性

圖6 降低溫包內(nèi)壓力溫度曲線斜率的熱力膨脹閥系統(tǒng)穩(wěn)定性

3 結(jié)論

本文對(duì)蒸發(fā)器-熱力膨脹閥回路的穩(wěn)定性進(jìn)行了研究，分析了熱力膨脹閥溫包充注特性對(duì)系統(tǒng)穩(wěn)定性的影響，得到如下結(jié)論：

1) 相同工作條件下，交叉充注的熱力膨閥系統(tǒng)不穩(wěn)定工作區(qū)域遠(yuǎn)小于平行充注的系統(tǒng)；

2) 采用交叉充注的熱力膨脹閥“交叉點(diǎn)”壓力對(duì)系統(tǒng)穩(wěn)定性有較大影響，交叉點(diǎn)壓力越高，系統(tǒng)的不穩(wěn)定工作區(qū)域越小，系統(tǒng)穩(wěn)定性越好；

3) 采用交叉充注的熱力膨脹閥，溫包內(nèi)充注工質(zhì)的壓力-溫度曲線斜率越低，系統(tǒng)不穩(wěn)定工作區(qū)域越小，系統(tǒng)穩(wěn)定性越好。

[1] HUELLE Z R. The MSS Line—a new approach to hunting problem[J]. ASHRAE Journal, 1972(10): 43-46.

[2] BROERSEN P M T. Control with a thermostatic expansion valve[J]. International Journal of Refrigeration,1982, 5(4): 209-212.

[3] IBRAHIM G A. Theoretical investigation into instability of a refrigeration system with an evaporator controlled by a thermostatic expansion valve[J]. The Canadian Journal of Chemical Engineering, 1998, 76(4): 722-727.

[4] IBRAHIM G A. Effect of sudden changes in evaporator external parameters on a refrigeration system with an evaporator controlled by a thermostatic expansion valve[J]. International Journal of Refrigeration, 2001,24(6): 566-576.

[5] TAHAT M A, IBRAHIM G A, PROBERT S D.Performance instability of a refrigerator with its evaporator controlled by a thermostatic expansion valve[J]. Applied Energy, 2001, 70(3): 233-249.

[6] DITTUS F W, BOELTER L M K. Heat transfer in single phase for plate finned tube heat exchangers[C]//Proceedings of the 5th International Heat Transfer Conference. San Francisco, 1980: 2275-2280.

[7] KANDILKAR S G. A model for correlating flow boiling heat transfer in augmented tubes and compace evaporators[J]. Journal of Heat Transfer, 1991, 113(4):966-972.

[8] MANGLIK R M, BERGLES A E. Heat transfer and pressure drop correlations for the rectangular offset strip fin compact heat exchanger[J]. Experimental Thermal and Fluid Science, 1995(10): 171-180.

[9] ROBERTSON J M, LOVEGROVE P C. Boling heat transfer with Freon 11 (R11) in brazed aluminum,plate-fin heat exchangers[J]. ASME Journal of Heat Transfer, 1983(105): 605-620.

[10] 葉奇昉. 制冷系統(tǒng)典型自控閥件動(dòng)力學(xué)研究[D]. 上海:上海交通大學(xué), 2010.