張浩，何哲旺，武滔,2，吳國(guó)明，任滔，丁國(guó)良*

（1-廣東美的制冷設(shè)備有限公司，廣東佛山 528311；2-重慶大學(xué)城市建議與環(huán)境工程學(xué)院，重慶 400044；3-上海交通大學(xué)制冷與低溫工程研究所，上海 200240）

0 引言

由于高效緊湊化空調(diào)器的發(fā)展趨勢(shì)，空調(diào)的室外換熱器通常采用更小管徑的銅管來提升空調(diào)器的換熱性能[1-2]?？照{(diào)室外換熱器的管徑減小會(huì)導(dǎo)致管內(nèi)壓降增加[3-4]，不利于換熱器的性能進(jìn)一步提升。目前的空調(diào)器通常采用多流路室外換熱器來降低管內(nèi)壓降并提升空調(diào)器的性能[5]。

多流路室外換熱器最常見的問題是常規(guī)的流路設(shè)計(jì)無法保證空調(diào)器在不同的冷媒流動(dòng)特性下始終處于最佳的換熱狀態(tài)[6-7]。制冷劑的物性，比如流速、溫度、壓力等等，在換熱器的進(jìn)口階段和出口階段相差很大[8-9]。不同的制冷劑物性具有不同的壓降特性，從而需要不同的分流方案才能達(dá)到最佳的換熱。目前室外換熱器通常采用均一分流數(shù)目的流路[10-12]，無法保證換熱器內(nèi)所有的制冷劑均具有最佳的分流特性和換熱特性。

解決換熱器內(nèi)冷媒物性不同導(dǎo)致的分流問題的方法是在換熱器的不同位置處采用不同的分流方案[13-14]。當(dāng)制冷劑的干度很大時(shí)，制冷劑的流速較大[15-16]，導(dǎo)致壓降很大；此時(shí)制冷劑需要較大的分流數(shù)目來減少壓降，以便達(dá)到更好的換熱。當(dāng)制冷劑的干度很小時(shí)，制冷劑的流速很小，導(dǎo)致制冷劑的換熱系數(shù)很小；此時(shí)制冷劑需要更少的分流數(shù)目來增大換熱系數(shù)。

根據(jù)冷媒物性的變化設(shè)計(jì)新型流路的難點(diǎn)在于制冷劑物性和對(duì)應(yīng)的分流特性之間的理論關(guān)系未知。新型流路需要根據(jù)制冷劑的物性精確計(jì)算對(duì)應(yīng)的分流特性[17-18]。但是目前還沒有關(guān)于制冷劑物性和分流特性的相關(guān)理論研究。因此有必要開展制冷劑流動(dòng)特性和分流特性的理論研究。

本文的目的是提出一種能夠提升空調(diào)換熱的室外機(jī)分布式流路。

1 分布式流路模型開發(fā)原理

室外換熱器分布式流路模型開發(fā)思路是：針對(duì)換熱微元建立分流模型，求解使得微元換熱量最大的臨界質(zhì)流密度；并通過換熱器分流模型求解分布式流路，使得每個(gè)換熱微元的制冷劑質(zhì)流密度盡可能地接近極限質(zhì)流密度。

臨界質(zhì)流密度Gcri是使得換熱微元換熱量最大的質(zhì)流密度，由管內(nèi)制冷劑干度和制冷劑流量決定。不同換熱器位置處的制冷劑具有不同的干度和流量，如圖1(a)所示。在流路的入口階段，制冷劑干度小，液相制冷劑體積遠(yuǎn)大于氣相制冷劑體積，使得管內(nèi)制冷劑流速很小。此時(shí)管內(nèi)具有很小的壓降和很小的換熱系數(shù)，如圖1(b)所示。這說明換熱器需要通過減小分流數(shù)目，來提高質(zhì)流密度G和換熱系數(shù)。當(dāng)制冷劑流到室外機(jī)管程的中部時(shí)，氣相制冷劑逐漸增多，導(dǎo)致制冷劑流速增大，使得管內(nèi)制冷劑的壓降和換熱系數(shù)同時(shí)升高，并達(dá)到一個(gè)最佳的匹配，如圖1(c)所示。當(dāng)制冷劑處于室外機(jī)管程的出口階段，制冷劑的干度很大，氣相制冷劑體積遠(yuǎn)大于液相制冷劑體積，使得管內(nèi)制冷劑流速非常大，如圖1(d)。此時(shí)制冷劑的壓降非常大，嚴(yán)重影響了管內(nèi)換熱的換熱溫差。

整個(gè)換熱器的分布式流路通過求解換熱器分流模型得到。換熱器分流模型以一根U管為最小換熱單元，分別求解每一根U管的最佳分流方案，最終得到整個(gè)換熱器的分布式流路。

圖1 室外換熱器制熱工況下不同管程的換熱壓降特性

2 換熱微元分流模型

首先建立換熱微元換熱的控制方程，確定換熱微元的邊界條件和補(bǔ)充方程，然后求解換熱微元分流模型得到使換熱微元換熱量最大的臨界質(zhì)流密度Gcri和對(duì)應(yīng)的最佳分流數(shù)Nopt。

2.1 換熱微元分流模型控制方程

換熱微元分流的控制方程包括 1個(gè)換熱量方程和3個(gè)能量平衡方程。換熱量方程計(jì)算換熱微元內(nèi)的空氣與制冷劑之間的換熱量，用于評(píng)價(jià)室外換熱器的換熱性能。能量平衡方程是分別基于空氣側(cè)熱平衡、制冷劑側(cè)熱平衡以及空氣-制冷劑之間熱平衡獲得。

換熱量方程是通過換熱微元內(nèi)的制冷劑焓差變化得到。換熱微元是沿?fù)Q熱器管程方向上選取的一個(gè)ds長(zhǎng)度的微小換熱段，如圖2(a)所示。換熱微元的換熱量由制冷劑的進(jìn)出口焓差和制冷劑質(zhì)量流量計(jì)算得到，如公式1和圖2(b)所示。

能量平衡方程是在換熱微元被分成N個(gè)流路后針對(duì)每個(gè)流路微元建立的，如圖3所示。換熱微元被分成N個(gè)流路后，每個(gè)流路微元的長(zhǎng)度為ds/N，管內(nèi)面積為Ai/N，管外面積為Ao/N，進(jìn)口制冷劑質(zhì)流密度為Gin/N，進(jìn)口空氣流量為mair/N。空氣側(cè)溫度變化導(dǎo)致的能量變化是由空氣與銅管外壁面換熱導(dǎo)致的，如公式2所示；制冷劑側(cè)焓差的變化是由制冷劑與銅管內(nèi)壁面換熱導(dǎo)致的，如公式3所示；在微元換熱穩(wěn)定后，空氣側(cè)的能量變化等于制冷劑側(cè)的能量變化，如公式4所示。換熱微元的4個(gè)控制方程中含有 7個(gè)未知參數(shù)，即 Twall、Tref,out、Tair,out、dQ、N、hi、Href,out。因此需要添加補(bǔ)充方程，才能求解換熱微元的控制方程組。

式中：

D——管徑，m；

ds——管程，m；

Ao——管外換熱面積，m2；

Ai——管內(nèi)換熱面積，m2；

Tair,in——空氣進(jìn)風(fēng)溫度，K；

mair——空氣流量，kg/s；

ho——空氣側(cè)換熱系數(shù)，kJ/(kg?K)；

Tair,out——出風(fēng)溫度，K；

G——冷媒質(zhì)流密度，kg/(m2?s)；

x——干度，無量綱；

p——壓力，Pa；

T——溫度，K；

hi——管內(nèi)換熱系數(shù)，W/(m2?K)。

圖2 換熱微元內(nèi)的結(jié)構(gòu)參數(shù)和換熱參數(shù)

圖3 換熱微元被分成N個(gè)流路后的參數(shù)變化

2.2 換熱微元分流模型補(bǔ)充方程

補(bǔ)充方程包括管內(nèi)換熱系數(shù) hi計(jì)算方程和管內(nèi)壓降Δp計(jì)算方程。管內(nèi)換熱系數(shù)hi采用GUNGOR的經(jīng)驗(yàn)公式計(jì)算[19]，其平均計(jì)算精度經(jīng)過試驗(yàn)驗(yàn)證為17.6%以內(nèi)。管內(nèi)壓降采用胡海濤的經(jīng)驗(yàn)公式計(jì)算[20]，其平均計(jì)算精度經(jīng)過試驗(yàn)驗(yàn)證為10.8%以內(nèi)。

2.3 換熱微元分流模型求解結(jié)果

4個(gè)控制方程和2個(gè)補(bǔ)充方程組成的換熱微元分流方程組中包括7個(gè)未知參數(shù)，即Twall、Tref,out、Tair,out、dQ、N、hi、Href,out。換熱微元分流方程組存在無數(shù)個(gè)解，即每一個(gè)分流數(shù)N均對(duì)應(yīng)一個(gè)微元換熱量 dQ。在無數(shù)個(gè)解當(dāng)中，存在一個(gè)分流數(shù) Nopt會(huì)產(chǎn)生最大的微元換熱量 dQmax。此時(shí)微元內(nèi)的質(zhì)流密度即為臨界質(zhì)流密度。

臨界質(zhì)流密度Gcri的求解方法是：

1）首先假設(shè)一個(gè)分流數(shù)N值，帶入 Δp補(bǔ)充方程，求解得制冷劑出口壓力 pref,out和制冷劑出口溫度Tref,out；

2）再將分流數(shù)N和制冷劑出口溫度Tref,out帶入控制方程和 hi補(bǔ)充方程，求解空氣出風(fēng)溫度Tair,out，制冷劑出口焓Hout和微元總換熱量dQ；

3）循環(huán)計(jì)算，直到總換熱量 dQ 達(dá)到最大，并輸出臨界質(zhì)流密度Gcri和最佳分流數(shù)Nopt，如圖4所示。

臨界質(zhì)流密度Gcri的求解結(jié)果如圖5所示。臨界質(zhì)流密度Gcri隨著冷媒干度增大而減小，表明換熱微元的最佳分流數(shù)目隨著制冷劑的干度增大而增大。

圖4 換熱微元最佳分流求解思路

圖5 換熱微元臨界質(zhì)流密度求解結(jié)果

3 換熱器分流模型

室外換熱器分布式流路通過求解換熱器分流模型得到。換熱器分流模型是以單位U管長(zhǎng)度的換熱量為目標(biāo)函數(shù)，求解使得單位U管長(zhǎng)度換熱量最大的分流方案，如公式5所示。

式中：

Q——換熱量，W；

ds——計(jì)算單元的長(zhǎng)度，m；

Ni——第i輪的分流數(shù)目；

L——U管長(zhǎng)度，m。

分布式分流模型進(jìn)行若干輪分流計(jì)算，并且每一輪各個(gè)分流支路的制冷劑均流過1個(gè)U管后進(jìn)入?yún)R總管等待下一輪的分流計(jì)算，如圖6(a)所示。在每一輪分流之前，分布式流路模型均進(jìn)行一次最佳分流數(shù)的求解，使得室外換熱器的任何一個(gè)U管內(nèi)的制冷劑質(zhì)流密度盡可能地接近臨界質(zhì)流密度Gcri，如圖6(b)所示。

圖6 室外換熱器分布式流路設(shè)計(jì)方法

圖7 室外換熱器分布式流路求解思路

換熱器分流模型的求解，如圖7所示，主要包括以下5個(gè)步驟：

1）分流前，進(jìn)行最佳分流數(shù)目求解，獲得使單位U管長(zhǎng)度的換熱量最大的分流數(shù)目，即最佳分流數(shù)目Ni；

2）求解最佳分流數(shù)目Ni后，制冷劑按照最佳分流數(shù)目分別流入Ni個(gè)U管；

3）在流過1個(gè)U管長(zhǎng)度后，所有制冷劑進(jìn)入?yún)R聚管等待下一次的分流；

4）計(jì)算下一次分流的最佳分流數(shù)目Ni+1，如果剩余U管數(shù)目大于Ni+1，則進(jìn)入下一次分流；

5）重復(fù)步驟1）至4），直到剩余U管數(shù)目不足，即剩余U管數(shù)目小于下一次分流數(shù)目Ni+1，則將剩余的U管連接到上一次分流中，室外換熱器分布式流路設(shè)計(jì)結(jié)束。

4 案例分析

為了更清楚地展示室外換熱器分布式流路的設(shè)計(jì)方法，本文提供一個(gè)案例的設(shè)計(jì)過程。

步驟1：輸入室外換熱器結(jié)構(gòu)參數(shù)和運(yùn)行參數(shù)，如表1所示。計(jì)算管內(nèi)面積Ai，管外面積Ao，管外換熱系數(shù)ho，制冷劑進(jìn)口壓力prefer,in，制冷劑進(jìn)口干度xin，制冷劑進(jìn)口焓Hin，制冷劑質(zhì)流密度Gin。

步驟2：根據(jù)制冷劑進(jìn)口參數(shù)求解第1輪分流的最佳分流數(shù)目 N1。計(jì)算結(jié)果表明最佳分流數(shù)目N1為2。故第1輪分流方案是：制冷劑分成兩個(gè)流路，分別流進(jìn)1個(gè)U管，最后進(jìn)入?yún)R聚管。

步驟3：計(jì)算第1輪分流出口制冷劑狀態(tài)參數(shù)，并根據(jù)制冷劑狀態(tài)參數(shù)計(jì)算第2輪分流的最佳分流數(shù)目N2。第2輪分流的最佳分流數(shù)目N2為3。故第2輪分流方案是：制冷劑分成3個(gè)流路，分別流進(jìn)1個(gè)U管，最后進(jìn)入?yún)R聚管。

步驟4：計(jì)算第2輪分流出口制冷劑狀態(tài)參數(shù)，并根據(jù)制冷劑狀態(tài)參數(shù)計(jì)算第3輪分流的最佳分流數(shù)目N3。第2輪分流的最佳分流數(shù)目N3為4。故第3輪分流方案是：制冷劑分成4個(gè)流路，分別流進(jìn)1個(gè)U管，最后進(jìn)入?yún)R聚管。

步驟5：剩余U管數(shù)目為零，分布式流路設(shè)計(jì)結(jié)束。室外換熱器的最終流路設(shè)計(jì)方案是2-3-4流路，如圖8(a)所示。

為了比較分布式流路和常規(guī)流路的換熱特性，本文將2-3-4分布式流路的換熱量與2種常規(guī)流路進(jìn)行比較。2種常規(guī)流路分別是2流路和3流路，分別如圖8(b)和8(c)所示。

室外換熱器流路與理論最佳流路之間的性能差異可以通過分析制冷劑沿程質(zhì)流密度與臨界質(zhì)流密度之間的差異來獲得。理論上最佳的流路設(shè)計(jì)應(yīng)該滿足管程任意位置處的質(zhì)流密度等于臨界質(zhì)流密度Gcri這個(gè)條件。但室外換熱器的最小流路長(zhǎng)度是1個(gè)U管，使得實(shí)際流路的制冷劑質(zhì)流密度呈階梯狀，如圖9(a)所示。當(dāng)實(shí)際流路的階梯狀的質(zhì)流密度越接近臨界質(zhì)流密度 Gcri，表明該流路分流越有利于換熱。2-3-4流路的質(zhì)流密度比2流路和3流路更接近臨界質(zhì)流密度Gcri，表明2-3-4流路的性能比2流路和3流路更好。計(jì)算結(jié)果表明采用2-3-4分布式流路的室外換熱器的額定制熱量和低溫制熱量分別比采用流路數(shù)目均一的 3流路的高 8.8%和5.6%，如圖9(b)所示。

表1 室外換熱器輸入?yún)?shù)

圖8 室外換熱器流路設(shè)計(jì)結(jié)果

圖9 不同流路室外換熱器性能分析

5 結(jié)論

1）采用分布式流路的室外換熱器的額定制熱量和低溫制熱量分別比采用流路數(shù)目均一的3流路提高8.8%和5.6%。

2）本文提出的分布式流路設(shè)計(jì)方法是：針對(duì)室外換熱器微元建立分流模型，求解使得微元換熱量最大的臨界質(zhì)流密度；并通過室外換熱器分布式分流模型求解分布式流路，使得室外換熱器的制冷劑質(zhì)流密度盡量接近極限質(zhì)流密度。

3）臨界質(zhì)流密度可以作為流路設(shè)計(jì)的評(píng)價(jià)標(biāo)準(zhǔn)。制冷劑在流路中的沿程質(zhì)流密度越接近臨界質(zhì)流密度，表明該流路越有利于換熱。