緱偉隆 劉雨聲 馮嘉旺 李嘉偉 李萬勇 陳江平,
(1 江蘇中關(guān)村科技產(chǎn)業(yè)園節(jié)能環(huán)保研究有限公司 常州 213300;2 上海交通大學(xué) 制冷與低溫工程研究所 上海 200240)
散熱器在汽車?yán)鋮s系統(tǒng)中具有極其重要的作用,平行流換熱器因結(jié)構(gòu)輕便、換熱性能高效,是汽車熱泵空調(diào)中最常用的室外換熱器,主要由集流管、扁管、翅片等結(jié)構(gòu)組成。換熱器翅片種類很多,生產(chǎn)中常用的翅片形式主要有平直翅片、鋸齒形翅片及多孔翅片等,目前在車用空調(diào)中最常采用的是對稱式百葉窗翅片換熱器。百葉窗翅片有增大傳熱表面積、擾動流體邊界層的作用,以此來達(dá)到強(qiáng)化換熱的目的。因此對百葉窗翅片的研究具有重要意義。
目前,國內(nèi)外學(xué)者對百葉窗式換熱器傳熱性能開展了大量的研究工作。仿真計(jì)算方面,J.Y.Jang等[1-2]對不同角度的百葉窗翅片換熱器進(jìn)行研究,結(jié)果表明,變角度的百葉窗翅片換熱器綜合性能優(yōu)于均勻角度的百葉窗翅片換熱器。郭健忠等[3]提出8種變角度翅片結(jié)構(gòu),分析了5種入口風(fēng)速下8種變角度翅片結(jié)構(gòu)的傳熱因子與阻力因子,研究8種變角度翅片結(jié)構(gòu)的優(yōu)劣。P.Karthik等[4-5]對翅片間距、橫管間距、縱管間距和百葉窗間距幾何參數(shù)進(jìn)行模擬計(jì)算,結(jié)果表明,橫管間距及縱管間距對傳熱系數(shù)的影響可忽略不計(jì),百葉窗間距過小會導(dǎo)致空氣阻力過大,翅片間距的減小會導(dǎo)致阻力增大及傳熱系數(shù)提高。H.Huisseune等[6]對帶有和不帶有渦流發(fā)生器的換熱器進(jìn)行了三維仿真計(jì)算,M.E.Nakhchi等[7]分析了不同傾斜角度的多孔百葉窗換熱器內(nèi)流體流動的湍流特性和熱性能,均認(rèn)為強(qiáng)化傳熱的主要物理原因是管壁和區(qū)域間流體的良好混合。楊林等[8]提出一種新型斜針形百葉窗翅片,得到其窗翅的最佳布置方式,為新型百葉窗式散熱器的設(shè)計(jì)優(yōu)化提供理論依據(jù)。
實(shí)驗(yàn)研究方面,董軍啟等[9]提出了關(guān)于傳熱因子j和摩擦因子f的關(guān)聯(lián)式,關(guān)聯(lián)式的最大誤差小于12%。A.Okbaz等[10-11]對百葉窗翅片及波紋翅片換熱器行了傳熱和阻力特性實(shí)驗(yàn),研究了空氣流速、翅片間距以及管排數(shù)的影響。H.Huisseune等[12]用染料注入技術(shù)對不同開窗角度的百葉窗翅片圓管換熱器進(jìn)行了6種放大模型的可視化實(shí)驗(yàn),結(jié)果表明,低雷諾數(shù)下,流線沿著管壁外緣;高雷諾數(shù)下,管前發(fā)展出一個馬蹄形漩渦,這種漩渦會削弱傳熱效果。A.Nuntaphan等[13]研究了空氣自然對流下百葉窗翅片換熱器的整體傾角對換熱性能的影響,結(jié)果表明,在35°~40° 的傾角下,換熱器的傳熱性能提高顯著。B.Dogan等[14]對雙排和三排翅片的百葉窗翅片扁管換熱器進(jìn)行了實(shí)驗(yàn)研究,利用水和空氣的單相對流換熱實(shí)驗(yàn)研究了不同翅片排數(shù)的換熱器空氣側(cè)的性能差距。
綜上可知,國內(nèi)外學(xué)者對對稱式百葉窗翅片研究較多,但在非對稱式百葉窗領(lǐng)域的研究很少,因此本文在傳統(tǒng)的對稱式百葉窗翅片的基礎(chǔ)上,設(shè)計(jì)出6種非對稱式百葉窗翅片。通過分析各個翅片的換熱特性和阻力特性,來改進(jìn)傳統(tǒng)翅片結(jié)構(gòu),降低換熱器所受風(fēng)阻,充分利用整個換熱器流道,最大程度實(shí)現(xiàn)冷、熱流體熱量交換。
百葉窗換熱器整體結(jié)構(gòu)復(fù)雜,對其進(jìn)行完整的仿真計(jì)算較為困難,但翅片結(jié)構(gòu)呈周期性變化,因此取其中一個單元結(jié)構(gòu)進(jìn)行仿真計(jì)算,如圖1所示。百葉窗翅片的熱量交換主要集中在前半段靠近進(jìn)口方向,在后半段區(qū)域熱量交換較少,由于翅片結(jié)構(gòu)的對稱性,所受風(fēng)阻大小基本與前半段相同,因此換熱器整體所受風(fēng)阻較高。針對以上問題,在對稱式百葉窗翅片的結(jié)構(gòu)基礎(chǔ)上(幾何模型如圖2所示,該模型的尺寸如表1所示),提出幾種非對稱式百葉窗翅片結(jié)構(gòu),尺寸參數(shù)如表2所示,建立的模型如圖3所示。
圖1 百葉窗換熱器整體結(jié)構(gòu)Fig.1 Overall structure of louver heat exchanger
表1 對稱式百葉窗結(jié)構(gòu)參數(shù)Tab.1 Structural parameters of symmetrical shutters
圖2 對稱式百葉窗結(jié)構(gòu)Fig.2 Structure of symmetrical shutter
表2 新型非對稱式百葉窗翅片結(jié)構(gòu)參數(shù)Tab.2 Structure parameters of new asymmetric shutter fins
圖3 新型非對稱式百葉窗翅片結(jié)構(gòu)Fig.3 Structure of the new asymmetric shutter fin
取單個百葉窗翅片進(jìn)行數(shù)值計(jì)算,采用Starccm+軟件進(jìn)行網(wǎng)格劃分。對換熱器的扁管進(jìn)行簡化,在翅片的上下方建立0.1 mm×3 mm×15 mm的矩形薄板代替扁管結(jié)構(gòu)。流體域設(shè)置為5.4 mm×3 mm×21 mm的長方體,建立百葉窗翅片計(jì)算域模型,如圖4所示。
圖4 百葉窗翅片計(jì)算域模型Fig.4 Calculation domain model of louver fin
對上述百葉窗翅片模型的翅片和空氣的傳熱問題進(jìn)行整體數(shù)值求解,因流動結(jié)構(gòu)較為復(fù)雜,對不同區(qū)域進(jìn)行不同的網(wǎng)格劃分,對流體域和固體結(jié)構(gòu)均進(jìn)行非結(jié)構(gòu)化網(wǎng)格計(jì)算和多面體網(wǎng)格劃分。為了保證計(jì)算的準(zhǔn)確度,在與固體壁面相鄰的流體區(qū)域均劃分了網(wǎng)格較密的附面層網(wǎng)格;因翅片厚度較薄,對翅片和上下扁管部分進(jìn)行薄體網(wǎng)格劃分。截面網(wǎng)格示意圖如圖5所示。
圖5 截面網(wǎng)格Fig.5 Sectional mesh
采用合理的網(wǎng)格參數(shù)設(shè)置對計(jì)算精度及模擬結(jié)果的影響顯著。本文選取A1模型進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證,對不同網(wǎng)格數(shù)量的模型進(jìn)行數(shù)值計(jì)算,如圖6所示,當(dāng)網(wǎng)格數(shù)大于300萬時(shí),網(wǎng)格數(shù)量對換熱量的影響基本不變。綜合考慮實(shí)際計(jì)算能力和結(jié)果精度,選用劃分?jǐn)?shù)量約為350萬的計(jì)算模型。
圖6 網(wǎng)格數(shù)對換熱量的影響Fig.6 Influence of mesh number on heat transfer
2.1.1 邊界條件設(shè)置
圖4所示為百葉窗翅片計(jì)算域模型,具體邊界條件設(shè)置如下:1)進(jìn)口設(shè)置為速度進(jìn)口條件,根據(jù)不同的工況設(shè)置相應(yīng)的入口的速度幅值和進(jìn)口溫度;2)出口設(shè)置為壓力出口;3)流體域的左右兩側(cè)面設(shè)置為周期性邊界;4)固體部分包括翅片和上下扁管,扁管的壁面設(shè)置為恒定溫度,固體域材料為鋁,流體區(qū)域與固體區(qū)域之間相互耦合;5)其余未說明的壁面均設(shè)置為絕熱面。
2.1.2 數(shù)據(jù)處理
對百葉窗翅片換熱器在不同進(jìn)口氣流速度和不同的扁管溫度下進(jìn)行空氣側(cè)傳熱和阻力性能的分析,采用的數(shù)據(jù)處理公式如下:
(1)
(2)
式中:de為特征長度,m;ua為流體速度,m/s;ρa(bǔ)為流體密度,kg/m3;λa為空氣導(dǎo)熱系數(shù),W/(m·K);ha為空氣側(cè)的表面?zhèn)鳠嵯禂?shù),W/(m2·K);μ為動力黏度,Pa·s。
傳熱因子j的表達(dá)式[15-16]為:
(3)
式中:cpa為質(zhì)量定壓熱容,J/(kg·K);Pr為普朗特?cái)?shù),取0.7。
百葉窗翅片對流換熱時(shí)空氣側(cè)的對數(shù)平均溫差:
(4)
式中:TW為上下扁管溫度,℃;T0為出口溫度,℃;Ti為進(jìn)口溫度,℃。
百葉窗翅片對流換熱時(shí)空氣側(cè)的表面?zhèn)鳠嵯禂?shù):
(5)
式中:ma為流體流過百葉窗翅片時(shí)空氣的質(zhì)量流量,kg/s;A為百葉窗翅片對流換熱時(shí)總的傳熱面積,m2。
阻力因子f的計(jì)算式[15-16]為:
(6)
式中:Δp為空氣側(cè)壓降,Pa;Ac為最窄截面面積,m2。
本文首先對對稱式百葉窗翅片進(jìn)行了仿真計(jì)算,將其作為對照組(A1),分析冷卻空氣進(jìn)口速度分別為1.5、2.5、3.5、4.5 m/s時(shí)4種工況下的百葉窗換熱模型,可得到空氣流經(jīng)翅片區(qū)域時(shí)的速度、壓力、溫度分布情況,以及換熱量與風(fēng)阻大小,為非對稱百葉窗翅片性能機(jī)理研究提供更直接的方法。然后按照相同的分析方法,對6種非對稱式百葉窗翅片進(jìn)行仿真計(jì)算,分析其換熱性能和阻力性能的優(yōu)劣。
B1、B2組為開窗數(shù)量不對稱的百葉窗翅片,B5、B6組為開窗高度逐漸變化的非對稱式百葉窗翅片,通過分析可知:對于B1、B2組,將開窗數(shù)量少的一側(cè)作為進(jìn)口方向其綜合性能優(yōu)于將開窗數(shù)量多的一側(cè)作為進(jìn)口方向;對于B5、B6組,氣流沿翅片高度逐漸降低方向通入的綜合性能優(yōu)于氣流沿翅片高度逐漸增高方向通入。因此,在后文的敘述中,默認(rèn)B1、B2組開窗數(shù)量少的一側(cè)作為進(jìn)口方向,B5、B6組百葉窗高度高的一側(cè)作為進(jìn)口方向。
1)為了直觀顯示新型翅片相比于傳統(tǒng)對稱式百葉窗翅片的性能優(yōu)劣,如圖7所示,將6種非對稱式百葉窗翅片傳熱因子曲線與A1對稱式百葉窗翅片傳熱因子曲線進(jìn)行對比,可以明顯看出7條曲線的趨勢基本一致,傳熱因子j均隨氣流速度的增大而減小。B1、B2組傳熱性能相比A1組有所提升,在進(jìn)口速度為4.5 m/s時(shí),B1、B2傳熱因子j比A1翅片分別提高1.2%和1.0%。說明在結(jié)構(gòu)相同的條件下,適當(dāng)?shù)恼{(diào)整布置方式,有助于提高百葉窗翅片換熱性能。B組其余翅片相比A1組換熱性能均有所降低,但對于B3、B4組開窗角度可變的翅片,隨著氣流速度的增加,減弱效果會顯著降低。在空氣進(jìn)口速度分別為1.5、2.5、3.5、4.5 m/s時(shí),B3翅片傳熱因子分別降低28.5%、24.7%、21.6%、19.5%。因?yàn)殡S著進(jìn)口氣流速度的增大,百葉窗表面邊界層遭到破壞,導(dǎo)致邊界層厚度變薄,增強(qiáng)換熱效果,說明開窗角度可變的設(shè)計(jì)在氣流速度較大的條件下優(yōu)勢更顯著。
圖7 傳熱因子曲線Fig.7 Curve of heat transfer factor
2)分析各種非對稱式百葉窗翅片流動性能,以相同的方法對各組翅片的阻力因子進(jìn)行處理,結(jié)果如圖8所示。由圖8可知,7條曲線的趨勢基本一致,阻力因子f均隨氣流速度的增大而減小。B1、B2組非對稱式百葉窗翅片的阻力因子相比傳統(tǒng)對稱式百葉窗翅片A1有所增加,但增幅很小;B3~B6組的阻力因子相比A1組均顯著降低,在進(jìn)口速度為4.5 m/s時(shí),阻力因子f分別降低31.9%、20.4%、30.4%、16.0%。尤其是百葉窗角度可變的新型翅片可顯著降低風(fēng)阻。
圖8 阻力因子曲線Fig.8 Curve of resistance factor
2.3.1 傳熱分析
研究新型翅片結(jié)構(gòu)設(shè)計(jì)對流道內(nèi)流體流動的溫度分布影響,各翅片截面的溫度分布云圖如圖9所示。由圖9可知,換熱翅片上百葉窗的存在會改變流體流動的方向,增大了流體的湍流效應(yīng),破壞流體熱邊界層,提升不同流道內(nèi)在相同截面下流體溫度的均勻性,強(qiáng)化換熱器內(nèi)冷熱流體熱量交換效果。A1、B1、B2翅片的整體溫度分布基本一致,流體與換熱器翅片壁面的熱量交換主要發(fā)生在進(jìn)口段,換熱流體溫升較為顯著。B5、B6翅片的熱量交換也集中在氣流進(jìn)口段,但由于后半段翅片結(jié)構(gòu)的設(shè)計(jì),熱量交換不顯著,后半段翅片維持較高的溫度。B3、B4翅片在進(jìn)出口處沒有明顯的溫度梯度變化,而且整個翅片沿流動方向相比其他翅片換熱更加均勻,符合對新型翅片換熱器的要求。
圖9 翅片溫度分布云圖Fig.9 Temperature contours of fin
為了更好的說明非對稱式翅片結(jié)構(gòu)設(shè)計(jì)對換熱均勻性的影響,本文分析了不同翅片結(jié)構(gòu)設(shè)計(jì)對流道內(nèi)溫度場變化的影響。沿著翅片氣流方向從進(jìn)口到出口設(shè)置5個監(jiān)測面,分別為0 mm面(進(jìn)口面)、3.75 mm面、7.5 mm面、11.25 mm面和15 mm面(出口面),如圖10所示。測量不同截面的平均溫度,如圖11所示。對比三種新型翅片與傳統(tǒng)對稱式百葉窗翅片的溫度變化曲線可知,百葉窗開窗角度可變的B3翅片溫度曲線更加平直,說明溫度變化更均勻,在整個翅片上熱量交換不存在明顯集中部位;B5非對稱式百葉窗翅片前半段溫度增長較快,后半段溫度增長較慢,說明B5翅片換熱不均勻,熱量交換主要集中在進(jìn)口段。
圖10 溫度監(jiān)測截面的選取Fig.10 Selection of temperature monitoring section
圖11 監(jiān)測面平均溫度Fig.11 Average temperature of monitoring surface
2.3.2 流動分析
分析氣流在流道內(nèi)的流動情況,進(jìn)口速度為4.5 m/s時(shí)各翅片的壓力分布云圖如圖12所示。左側(cè)為進(jìn)口,右側(cè)為出口。由圖12可知,對于風(fēng)阻較大的A1、B1、B2翅片,阻力在進(jìn)口處較為集中,在翅片出口處也有顯著的壓力梯度變化。原因是在進(jìn)出口處,氣流遇到較大的轉(zhuǎn)折角,由于翅片破壞氣流邊界層的作用受翅片開窗角度θ的影響,θ越大,即轉(zhuǎn)折角越大,破壞氣流邊界層的作用越強(qiáng),所受風(fēng)阻也相應(yīng)增大。對于開窗角度變化的B3、B4翅片,進(jìn)出口處的阻力分布較為均勻,沒有明顯的風(fēng)阻集中現(xiàn)象,說明百葉窗角度可變的換熱器翅片設(shè)計(jì)對局部阻力的控制具有顯著優(yōu)勢。對于B5、B6翅片開窗高度逐漸降低的非對稱式百葉窗,風(fēng)阻主要集中在前半段,在后半段壓力降低不顯著,尤其是B5翅片,在后半段所受風(fēng)阻很小。因?yàn)閷τ贐5、B6開窗高度可變的翅片設(shè)計(jì),在后半段百葉窗高度逐漸降低,空氣與翅片的接觸面積減小,氣流流動更平順,風(fēng)阻也相應(yīng)降低。
圖12 翅片壓力分布云圖Fig.12 Pressure contours of fin
綜上可得,B組非對稱式百葉窗翅片B3~B6雖然降低了一定的傳熱能力,但可顯著降低風(fēng)阻,尤其是B3、B5翅片風(fēng)阻降低更顯著。在換熱均勻性方面,B3翅片相比于傳統(tǒng)翅片和其他新型翅片,換熱更均勻,不存在換熱集中在特定位置的問題。因此對于流阻要求較為苛刻且對換熱要求均勻的場景,B3翅片的設(shè)計(jì)具有顯著優(yōu)勢。
為驗(yàn)證數(shù)值模擬的可靠性進(jìn)行實(shí)驗(yàn)研究。通過3D打印建立百葉窗整體換熱器模型,如圖13所示,在風(fēng)量臺架上進(jìn)行風(fēng)阻測試。
圖13 整體換熱器實(shí)驗(yàn)?zāi)P虵ig.13 Experimental model of integral heat exchanger
為驗(yàn)證本文CFD仿真結(jié)果的可靠性,將所得的實(shí)驗(yàn)結(jié)果進(jìn)行處理,計(jì)算得到單個換熱器翅片的阻力因子,以同樣的方式繪制各個非對稱式百葉窗翅片的阻力因子曲線,并分析其流動性能。圖14所示為百葉窗翅片阻力因子實(shí)驗(yàn)數(shù)據(jù)和模擬結(jié)果的對比。由圖14可知,百葉窗翅片阻力因子曲線的實(shí)驗(yàn)結(jié)果和仿真結(jié)果的變化趨勢一致,模擬值與實(shí)驗(yàn)值的最大誤差為16%,除個別點(diǎn)誤差稍大外,其余計(jì)算點(diǎn)誤差均小于10%,在工程允許范圍內(nèi)。
圖14 阻力因子實(shí)驗(yàn)數(shù)據(jù)和模擬結(jié)果對比Fig.14 Comparison of experimental data and simulation results of resistance factor
為研究微通道平行流室外換熱器的翅片結(jié)構(gòu)優(yōu)化問題,基于換熱器換熱機(jī)理,本文提出6種非對稱式百葉窗翅片結(jié)構(gòu),對該新型翅片進(jìn)行傳熱特性和流動特性分析,得到如下結(jié)論:
1)提出的6種非對稱式百葉窗翅片結(jié)構(gòu),合理設(shè)置了百葉窗的高度、開窗角度以及開窗數(shù)量,對傳統(tǒng)對稱式百葉窗翅片換熱不均,翅片局部風(fēng)阻集中明顯的問題進(jìn)行了優(yōu)化。
2)采用三維建模對百葉窗翅片進(jìn)行數(shù)值模擬,得到百葉窗翅片截面處的溫度分布云圖和壓力分布云圖,研究其傳熱特性和流動特性,并對該新型翅片進(jìn)行實(shí)驗(yàn)驗(yàn)證,實(shí)驗(yàn)結(jié)果和仿真結(jié)果的誤差集中在10%以下,誤差在工程允許范圍內(nèi)。
3)通過MATLAB進(jìn)行數(shù)據(jù)計(jì)算,得到翅片傳熱因子和阻力因子曲線,與傳統(tǒng)對稱式百葉窗翅片性能曲線進(jìn)行對比分析,結(jié)果表明,在進(jìn)口速度為4.5 m/s時(shí),本文設(shè)計(jì)的B3翅片的阻力因子f相比傳統(tǒng)百葉窗翅片減小31.9%,且整個翅片換熱更加均勻,不存在換熱集中在特定位置的問題,適用于對流阻要求較為苛刻且對換熱要求均勻的場景。