梁忠偉，周俊輝，劉曉初

（廣州大學(xué)a.機(jī)械與電氣工程學(xué)院，b廣州市金屬材料強(qiáng)化研磨高性能加工重點(diǎn)實(shí)驗(yàn)室，c.廣東省水肥高效利用及太陽(yáng)能智能灌溉工程技術(shù)研究中心，廣東廣州　510006）

基于正交實(shí)驗(yàn)的空調(diào)換熱器流場(chǎng)模擬及性能分析

梁忠偉a，b，c，周俊輝a，劉曉初a，b，c

（廣州大學(xué)a.機(jī)械與電氣工程學(xué)院，b廣州市金屬材料強(qiáng)化研磨高性能加工重點(diǎn)實(shí)驗(yàn)室，c.廣東省水肥高效利用及太陽(yáng)能智能灌溉工程技術(shù)研究中心，廣東廣州510006）

文章運(yùn)用Fluent軟件對(duì)空調(diào)中換熱器氣態(tài)流場(chǎng)進(jìn)行數(shù)值模擬分析.在進(jìn)行計(jì)算網(wǎng)格劃分以及建立控制方程基礎(chǔ)上，對(duì)其溫度場(chǎng)、速度場(chǎng)、壓力場(chǎng)等關(guān)鍵性能要素進(jìn)行數(shù)值模擬及過(guò)程實(shí)驗(yàn).通過(guò)正交實(shí)驗(yàn)獲得復(fù)雜流場(chǎng)性能的動(dòng)態(tài)描述，有利于對(duì)其性能參數(shù)的數(shù)值分布及演變進(jìn)程提出量化機(jī)理方案.該研究為后續(xù)開(kāi)展的空調(diào)換熱器性能以及結(jié)構(gòu)參數(shù)優(yōu)化提供了理論基礎(chǔ)和數(shù)據(jù)準(zhǔn)備.

正交實(shí)驗(yàn)；空調(diào)換熱器；流場(chǎng)模擬；性能分析

1　換熱器流場(chǎng)研究的必要性及進(jìn)展

空調(diào)換熱器作為汽車(chē)空調(diào)系統(tǒng)的關(guān)鍵部件之一，其氣體流場(chǎng)動(dòng)態(tài)過(guò)程及工作性能對(duì)于空調(diào)系統(tǒng)結(jié)構(gòu)參數(shù)的優(yōu)化及壽命提高具有重要的指導(dǎo)意義，尤其是換熱工作過(guò)程更是決定了空調(diào)系統(tǒng)制冷量大小和整體性能優(yōu)劣.在換熱器空氣側(cè)性能研究方面，國(guó)內(nèi)外學(xué)者開(kāi)展了大量研究，并將主要研究方向集中在了空氣側(cè)結(jié)構(gòu)參數(shù)對(duì)換熱性能和管路壓降的影響上.文獻(xiàn)［1］通過(guò)實(shí)驗(yàn)發(fā)現(xiàn)，隨著換熱器中百葉窗角度的增大，換熱系數(shù)和壓降呈現(xiàn)不斷增長(zhǎng)的趨勢(shì)；與此同時(shí)，換熱系數(shù)和壓降會(huì)隨著百葉窗間距的增大而降低［2］，并隨著百葉窗長(zhǎng)度的增大而增加［3］；而換熱系數(shù)和壓降也將隨著翅片高度的增大而降低［4］.文獻(xiàn)［5］基于實(shí)驗(yàn)研究基礎(chǔ)上，提出了使得換熱系數(shù)最大化的百葉窗最優(yōu)角度.在目前所進(jìn)行的相關(guān)研究中，空氣側(cè)結(jié)構(gòu)參數(shù)基本涵蓋了百葉窗翅片的4個(gè)常規(guī)參數(shù)：翅片間距，翅片高度，翅片長(zhǎng)度和百葉窗角度，并有學(xué)者采用了控制變量法進(jìn)行最優(yōu)參數(shù)的迭代選擇［6－8］.另一方面，針對(duì)平行流換熱器的制冷劑性能提升，國(guó)外學(xué)者主要關(guān)注集流管內(nèi)流量分配均勻度等關(guān)鍵問(wèn)題［9－10］.經(jīng)過(guò)數(shù)據(jù)分析，業(yè)界普遍認(rèn)為集流管流量在扁管內(nèi)分配的不均勻性將會(huì)導(dǎo)致空調(diào)換熱性能的相應(yīng)弱化.有國(guó)內(nèi)學(xué)者采用了模擬實(shí)驗(yàn)和現(xiàn)場(chǎng)實(shí)驗(yàn)相結(jié)合的方法，對(duì)制冷劑中氣液兩相流均勻分布的溫度場(chǎng)進(jìn)行了探索，提出了獨(dú)到見(jiàn)解.同時(shí)也有學(xué)者采用了實(shí)驗(yàn)?zāi)M方法對(duì)控制扁管的結(jié)構(gòu)參數(shù)進(jìn)行了換熱器動(dòng)態(tài)性能分析，研究成果顯著提升了換熱器的工作效能［11－13］.盡管現(xiàn)有的大量研究為空調(diào)的性能提升和結(jié)構(gòu)優(yōu)化提供了理論支持和技術(shù)準(zhǔn)備，但是在氣體流場(chǎng)對(duì)于換熱器性能的影響機(jī)理、平行流內(nèi)部的復(fù)雜瞬時(shí)模擬，以至于多相相變流動(dòng)的實(shí)時(shí)改善優(yōu)化等多個(gè)方面，直至目前仍然無(wú)法獲得清晰可靠的理論解釋.這無(wú)疑為空調(diào)換熱器性能的進(jìn)一步顯著提升構(gòu)成了機(jī)理性障礙，因而在此方面需要進(jìn)一步深入研究，并提出有效可行的解決方法.

2　換熱器空氣側(cè)流場(chǎng)分析

2.1控制方程與計(jì)算方法

空調(diào)層流模擬控制方程如下所示：

式中，μ，v，w為流體在XYZ方向的速度，單位m· s－1；μ為流體的粘性系數(shù)，單位kg/（m·s－1）；P為流體的壓力，單位Pa；T為流體的溫度，單位K；λ為流體的導(dǎo)熱系數(shù)，單位W/（m·k－1）；Cp為流體的比熱容，單位J/（kg·k－1）［14］.

2.2數(shù)值模型建立

圖1　百葉窗空氣側(cè)結(jié)構(gòu)示意圖Fig.1　Structure of air side shutter

圖2　翅片主要結(jié)構(gòu)Fig.2　Themain structural parameters of fins

圖1為平行流換熱器百葉窗空氣側(cè)結(jié)構(gòu)示意圖.高溫高壓的制冷劑（水蒸氣）經(jīng)過(guò)集流管流進(jìn)多孔扁管微通道，經(jīng)扁管節(jié)流的水蒸氣與管壁上下表面的百葉窗翅片進(jìn)行熱量交換；與此同時(shí)，空氣流經(jīng)百葉窗翅片通道與之進(jìn)行對(duì)流換熱，最終吸收翅片周?chē)諝鉄崃?，達(dá)到冷卻空氣的目的.圖2為平行流換熱器的結(jié)構(gòu)示意圖，百葉窗翅片表面有許多有角度的翅片，空氣通過(guò)這些開(kāi)窗翅片后，流動(dòng)方向會(huì)發(fā)生改變，邊界層也會(huì)不斷地破壞后重組，這樣會(huì)增加空氣的擾動(dòng)，從而使傳熱性能更好［15］.其結(jié)構(gòu)參數(shù)見(jiàn)表1.

表1　百葉窗翅片主要結(jié)構(gòu)參數(shù)Table 1　Main structure parameters of shutter fin

表2　計(jì)算網(wǎng)格敏感性驗(yàn)證Table 2　The sensitivity verification for computational grid

在計(jì)算網(wǎng)格建立過(guò)程中，模型網(wǎng)格尺寸越小，則理論上計(jì)算精度越高，模擬仿真效果及數(shù)值運(yùn)算結(jié)果也越逼近真實(shí)值，且計(jì)算時(shí)間將成幾何次方延長(zhǎng).但是同時(shí)由于受到計(jì)算能力限制，為了有效利用現(xiàn)有的計(jì)算能力，保證計(jì)算網(wǎng)格的可靠性及穩(wěn)定性，需要對(duì)其敏感性開(kāi)展驗(yàn)證，并據(jù)此選擇合適網(wǎng)格大小.在本研究中，針對(duì)可采用的7種網(wǎng)格大小方案分別取單元長(zhǎng)度為0.15，0.20，0.25，0.30，0.35，0.40及0.45mm開(kāi)展驗(yàn)證研究，以同一外界工作環(huán)境下標(biāo)記多個(gè)測(cè)試點(diǎn)所獲得的溫度、速度及壓力值瞬態(tài)值分布方差，以及達(dá)到穩(wěn)定值所需的收斂時(shí)間和收斂運(yùn)算步數(shù)作為計(jì)算收斂的判斷標(biāo)準(zhǔn)，從中綜合確定最優(yōu)的計(jì)算模型網(wǎng)格尺寸.而采用不同網(wǎng)格尺寸情況下，所獲得的溫度、速度及壓力值瞬態(tài)值分布方差見(jiàn)表2.由表2可見(jiàn)，當(dāng)網(wǎng)格尺寸大于0.15mm且逐漸增大時(shí)，溫度、速度及壓力值瞬態(tài)值分布方差均保持較低水平，至0.25mm時(shí)達(dá)到最低水平.與此同時(shí)，其所對(duì)應(yīng)的評(píng)價(jià)指標(biāo)收斂時(shí)間達(dá)到最短，而獲得穩(wěn)定值的收斂運(yùn)算步數(shù)也最少，反映出其參數(shù)并未隨著計(jì)算網(wǎng)格變化而產(chǎn)生顯著變化，表明在這區(qū)間內(nèi)計(jì)算模型精確性及仿真過(guò)程穩(wěn)定性不受網(wǎng)格尺寸大小的顯著影響.而當(dāng)網(wǎng)格尺寸小于0.15mm或者大于0.35mm時(shí)，各評(píng)價(jià)指標(biāo)的瞬時(shí)值方差顯著加大，而達(dá)到穩(wěn)定評(píng)價(jià)指標(biāo)的收斂時(shí)間及運(yùn)算步數(shù)也顯著增加，這顯示了在該網(wǎng)格尺寸區(qū)間內(nèi)，相關(guān)評(píng)價(jià)指標(biāo)的分布方差、收斂時(shí)間以及收斂運(yùn)算步數(shù)等將隨著網(wǎng)格大小而呈現(xiàn)出顯著變化，意味著此時(shí)計(jì)算網(wǎng)格的準(zhǔn)確性及可靠性對(duì)于網(wǎng)格尺寸變化將更為敏感，這對(duì)于保持模擬仿真過(guò)程的可靠性和準(zhǔn)確性尤為不利［16］.結(jié)合實(shí)驗(yàn)環(huán)境及計(jì)算條件，以上敏感性驗(yàn)證結(jié)果表明將網(wǎng)格尺寸選為0.25mm更有助于后續(xù)的流場(chǎng)模擬及數(shù)據(jù)分析，并獲得更優(yōu)的性能運(yùn)算效果，如表2中加粗字體所示.

圖3　平行流換熱器空氣側(cè)計(jì)算模型Fig.3　The computational air model of parallel flow heat exchanger

圖4　計(jì)算模型網(wǎng)格劃分圖Fig.4　The grid meshing structure of calculationmodel

圖5　換熱器邊界條件圖Fig.5　Boundary condition of exchanger

簡(jiǎn)化后的平行流換熱器Proe計(jì)算模型見(jiàn)圖3，主要由空氣流域和百葉窗翅片組成.通過(guò)敏感性驗(yàn)證后網(wǎng)格單元尺寸確定為0.25mm，總網(wǎng)格數(shù)為978 000個(gè).劃分后的網(wǎng)格模型見(jiàn)圖4.本文計(jì)算流域設(shè)置見(jiàn)圖5.計(jì)算區(qū)域包括空氣出入口和翅片接觸的空氣流域；入口采用速度入口邊界，入口溫度設(shè)為常溫300K；為加快計(jì)算收斂速度以提高效率，設(shè)出口為壓力出口，壓力為0Pa；并設(shè)置默認(rèn)壁面邊界和對(duì)稱(chēng)邊界；空氣流域的上下表面均設(shè)為默認(rèn)周期性邊界.

2.3數(shù)值模擬分析

通過(guò)改變空氣入口迎面風(fēng)速對(duì)模型進(jìn)行了從1.5m ·s－1到5.5m·s－1之間的5種工況計(jì)算.

（1）溫度場(chǎng)分布

圖6是空氣入口迎面風(fēng)速為1.5m·s－1、3.5m· s－1、5.5m·s－1的換熱器空氣側(cè)翅片對(duì)稱(chēng)面上的溫度分布云圖，可見(jiàn)：①在空氣入口處溫度的降低速度比出口段慢，因此，空氣側(cè)的熱交換主要集中在后半段；②溫度入口迎面風(fēng)速越低，出口處空氣溫度也越低，即空氣溫度降幅也越大（如圖6中的圓圈所示）.這是因?yàn)樵谄渌麠l件不變的情況下，迎面風(fēng)速的減小會(huì)導(dǎo)致空氣質(zhì)量流量的減小，單位空氣的換熱量則越大，因此空氣溫度降低幅度也就越大.③在前半段，3種工況下?lián)Q熱器翅片間的溫度基本都在350K以上（如圖6中的矩形方框所示），這顯示了空氣基本不從翅片間隙間流動(dòng).造成該現(xiàn)象的原因首先為翅片的間距太小，以至于對(duì)空氣造成過(guò)大的阻力；其次，翅片的角度過(guò)小也不利于空氣的流通.

圖6　換熱器翅片對(duì)稱(chēng)面溫度分布云圖Fig.6　The temperature distribution contour on the symmetry plane of heat exchanger fin

（2）速度場(chǎng)分布

圖7顯示了入口迎面風(fēng)速為1.5m·s－1和5.5m· s－1時(shí)，對(duì)稱(chēng)面上的局部速度流線(xiàn)圖，可見(jiàn)在入口速度為1.5m·s－1時(shí)空氣在翅片間隙分布的密度比速度為5.5m·s－1時(shí)大，且在轉(zhuǎn)向區(qū)都會(huì)改變流向，這說(shuō)明百葉窗結(jié)構(gòu)能夠增加空氣的擾動(dòng)，對(duì)其散熱性能是有利的.

（3）壓力場(chǎng)分布

圖8、圖9分別為空氣入口處迎風(fēng)速度為2.5m· s－1和4.5m·s－1時(shí)，對(duì)稱(chēng)面上的壓力等值線(xiàn)分布圖.由圖可知，①空氣入口迎風(fēng)速度為2.5m·s－1時(shí)，入口處的壓力為16.41Pa，中間轉(zhuǎn)向段為8.09Pa，出口處為－3.16Pa（如圖8）.即翅片前半段的壓降為8.32Pa，后半段的壓降為11.48Pa，總壓降為19.8Pa.此時(shí)，壓降主要集中在翅片的后半段.②空氣入口處速度為4.5m·s－1時(shí)，入口處的壓力為46.55Pa，中間轉(zhuǎn)向段為20.17Pa，出口處為－7.67Pa（如圖9）.前半段的壓降為26.38Pa，后半段為27.84Pa，總壓降為54.22Pa.前后壓降值更接近.因此，隨著風(fēng)速的增大，翅片壓降前后段的分布會(huì)更均勻，總壓降亦呈增大趨勢(shì).

圖7　入口風(fēng)速為1.5m·s－1和5.5m·s－1時(shí)對(duì)稱(chēng)面上的速度流線(xiàn)圖Fig.7　The speed flow map on the symmetry plane by inletwind speed of1.5m·s－1and 5.5m·s－1

圖8　入口迎面風(fēng)速為2.5m·s－1對(duì)稱(chēng)面上的壓力等值線(xiàn)分布圖Fig.8　The pressure contourmap on the symmetry plane by entrance wind speed of2.5m·s－1

圖9　入口迎面風(fēng)速為4.5m·s－1對(duì)稱(chēng)面上的壓力等值線(xiàn)分布圖Fig.9　The pressure contourmap on the symmetry plane by entrance wind speed of 4.5m·s－1

圖10為本文所模擬的5種不同工況下風(fēng)速對(duì)換熱器翅片換熱性能影響的趨勢(shì)圖.可見(jiàn)換熱器空氣側(cè)的換熱系數(shù)會(huì)隨著空氣入口的迎面風(fēng)速的增大而增大.而圖11為本文所模擬的5種不同工況下風(fēng)速對(duì)換熱器翅片壓降的影響趨勢(shì)圖.由圖可知，換熱器壓降隨著入口迎面風(fēng)速的增大而增大.

圖10　風(fēng)速對(duì)換熱器空氣側(cè)換熱系數(shù)影響Fig.10　Wind speed impacton heatexchanger in air side heat transfer coefficient

3　換熱器參數(shù)正交實(shí)驗(yàn)

為對(duì)平行流換熱器結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)選，需要從主要結(jié)構(gòu)參數(shù)中分析其對(duì)換熱器性能的影響；并選出最優(yōu)換熱器性能的結(jié)構(gòu)參數(shù)組合.結(jié)合實(shí)際正交實(shí)驗(yàn)要素，采用換熱器的換熱系數(shù)ha和摩擦因子f來(lái)表征其換熱性能.換熱系數(shù)越大，摩擦因子越小的結(jié)構(gòu)組合換熱器，綜合換熱性能越好［17－18］.

圖11　風(fēng)速對(duì)換熱器空氣側(cè)壓降的影響Fig.11　The influence ofwind speed on heat exchanger in air side pressure drop

3.1實(shí)驗(yàn)器材準(zhǔn)備

圖12（見(jiàn)封二）以圖示方式對(duì)于空調(diào)換熱器結(jié)構(gòu)進(jìn)行了說(shuō)明.在本實(shí)驗(yàn)中，結(jié)合具體的實(shí)驗(yàn)環(huán)境以及便于換熱過(guò)程數(shù)據(jù)參數(shù)的采集及監(jiān)控，本研究中采用真實(shí)尺寸的空調(diào)換熱器模型開(kāi)展流場(chǎng)優(yōu)化實(shí)驗(yàn).這能有效彌補(bǔ)真實(shí)換熱器工作過(guò)程中難以采集流場(chǎng)過(guò)程數(shù)據(jù)，難以對(duì)其換熱工作流場(chǎng)過(guò)程進(jìn)行實(shí)時(shí)監(jiān)控及參數(shù)標(biāo)定等多個(gè)方面的嚴(yán)重不足，更加有利于空調(diào)換熱器結(jié)構(gòu)的實(shí)驗(yàn)優(yōu)化及性能改善.其他實(shí)驗(yàn)部件/器材如圖12所示，圖12（a）中所示的發(fā)生器共有4個(gè)出口，分別為稀溶液入口、水蒸氣出口、濃溶液出口和抽真空出口［19］.在發(fā)生器內(nèi)部接有電加熱管，用來(lái)模擬汽車(chē)余熱對(duì)稀溶液進(jìn)行加熱.圖12（b）是冷凝器總成，由冷凝器芯體和冷凝風(fēng)機(jī)組成，用來(lái)對(duì)發(fā)生器出來(lái)的水蒸氣進(jìn)行散熱冷凝.圖12（c）所示的蒸發(fā)器總成主要由蒸發(fā)器芯體、蒸發(fā)風(fēng)機(jī)、膨脹閥及支撐外殼組成，用來(lái)蒸發(fā)水蒸氣吸熱降溫.圖12（d）、（e）、（f）均為實(shí)驗(yàn)所需器材，其中，圖12（d）所示真空泵用于系統(tǒng)真空度的保證，圖12（e）所示數(shù)顯溫度計(jì)用于檢測(cè)蒸發(fā)器出風(fēng)口溫度，以測(cè)試系統(tǒng)性能，圖12（f）所示蓄電池用于為系統(tǒng)提供電能.而表3所示為部分實(shí)驗(yàn)臺(tái)部件/器材.

3.2影響因素及水平確定

確定翅片寬度A、翅片間距B、翅片長(zhǎng)度C、百葉窗角度D為翅片的主要實(shí)驗(yàn)結(jié)構(gòu)參數(shù).各參數(shù)的取值如表4，翅片間距、翅片寬度和百葉窗角度均間隔2.5，翅片高度間隔1.每個(gè)因素分別取5個(gè)水平，見(jiàn)表5.由于本文所研究的平行流換熱器百葉窗結(jié)構(gòu)參數(shù)優(yōu)選問(wèn)題是4因素、5水平的問(wèn)題，因此選擇L25（56）的正交表（見(jiàn)表6），共有25種不同的計(jì)算模型.

3.3實(shí)驗(yàn)結(jié)果與分析

根據(jù)所選的正交實(shí)驗(yàn)表安排如下25組正交實(shí)驗(yàn)，共有25種不同結(jié)構(gòu)參數(shù)的百葉窗翅片參數(shù)組合（表6）.由于空調(diào)換熱器的迎面風(fēng)速一般為4.5m·s－1［20］，由實(shí)驗(yàn)所得的壓降與溫差計(jì)算其評(píng)價(jià)指標(biāo).表7為25組正交實(shí)驗(yàn)所求得評(píng)價(jià)指標(biāo)值.

為評(píng)價(jià)翅片各結(jié)構(gòu)參數(shù)對(duì)平行流換熱器空氣側(cè)翅片摩擦因子f和換熱系數(shù)的影響順序，在所選定的結(jié)構(gòu)參數(shù)中找出最優(yōu)組合［21－22］，根據(jù)表8的正交實(shí)驗(yàn)結(jié)果可以算出評(píng)價(jià)指標(biāo)摩擦因子和換熱系數(shù)的均值和極差（見(jiàn)表9）：由于極差表示各因素對(duì)評(píng)價(jià)指標(biāo)的影響程度，在此可以極差的大小判斷各因素對(duì)換熱器性能的影響先后順序［23－24］，極差越大的因素對(duì)評(píng)價(jià)指標(biāo)的影響也就越大，該因素也就越重要.

表3　實(shí)驗(yàn)臺(tái)部件/器材表Table 3　Component for equipment

表4　百葉窗結(jié)構(gòu)參數(shù)取值Table 4　The selected structure parameter for shutter fin

由表9極差R的大小和圖13（a，見(jiàn)封二）可見(jiàn)：翅片各結(jié)構(gòu)參數(shù)對(duì)空氣側(cè)摩擦因子的影響程度大小排序：翅片間距A（38.87％）＞百葉窗角度D（25.18％）＞翅片高度C（20.79％）＞翅片長(zhǎng)度B（15.15％），且翅片間距對(duì)摩擦因子的影響遠(yuǎn)大于其他3個(gè)因素，百葉窗角度和翅片高度對(duì)其影響程度相差不大，但都強(qiáng)于翅片高度.從圖13（b，見(jiàn)封二）中可見(jiàn)，各結(jié)構(gòu)參數(shù)對(duì)空氣側(cè)翅片換熱系數(shù)的影響程度的大小順序：翅片間距A（56.80％）＞翅片長(zhǎng)度B（18.05％）＞翅片高度C（13.37％）＞百葉窗角度D（11.79％），可以看出翅片間距對(duì)其影響程度仍然遠(yuǎn)強(qiáng)于其他因素.

表5　百葉窗結(jié)構(gòu)因素水平Table 5　The factor levels for shutter fin structure

表6　L25（56）正交表Table 6　Orthogonal table for L25（56）

表7　平行流換熱器百葉窗結(jié)構(gòu)參數(shù)實(shí)驗(yàn)方案Table 7　The parameter experimental scheme for parallel flow heat exchanger shutter

表8　正交實(shí)驗(yàn)結(jié)果Table 8　Results for orthogonal experiment

表9　正交實(shí)驗(yàn)結(jié)果極差分析表Table 9　Analysis for the orthogonal experimental results

3.4平行流換熱器百葉窗結(jié)構(gòu)參數(shù)的優(yōu)選

圖14、15為各因素水平對(duì)應(yīng)的摩擦因子和換熱系數(shù)的均值變化趨勢(shì)分布，表明了摩擦因子和換熱系數(shù)隨各因素的變化情況［25］.從中可見(jiàn)，①翅片間距對(duì)摩擦因子的影響最大，隨著翅片間距的增大，摩擦因子呈現(xiàn)逐漸增大后趨于平穩(wěn)的趨勢(shì)；②隨著翅片長(zhǎng)度的增大，摩擦因子是逐漸減小的；③摩擦因子隨著翅片參數(shù)的增大呈現(xiàn)先減小后增大的趨勢(shì)，隨著百葉窗角度的增大整體呈現(xiàn)增大的趨勢(shì).

圖14　因素水平對(duì)摩擦因子的影響Fig.14　The influences of factor levels on friction coefficient

圖15　因素水平對(duì)換熱系數(shù)的影響Fig.15　The influences of factor levels on heat exchange coefficients

由圖15可知，①翅片間距對(duì)換熱系數(shù)的影響程度也是最大的，隨著翅片間距的增大，換熱系數(shù)基本呈現(xiàn)直線(xiàn)上升的狀態(tài)；②隨著翅片長(zhǎng)度的增大，換熱系數(shù)亦緩慢增大；③翅片高度和百葉窗角度對(duì)換熱系數(shù)的影響趨勢(shì)基本一致，隨著翅片高度和百葉窗角度的增大，換熱系數(shù)呈現(xiàn)先增大后減小的趨勢(shì).

由于摩擦因子與換熱器性能的大小呈反比，而換熱系數(shù)與之成正比［23－25］.綜合上述分析得出最優(yōu)水平組合為組合25（翅片間距A＝2，翅片長(zhǎng)度B＝30，翅片高度C＝9，百葉窗角度D＝25°）.

4　結(jié)　論

本文運(yùn)用Fluent軟件對(duì)平行流換熱器空氣側(cè)進(jìn)行流場(chǎng)數(shù)值模擬分析.首先以平行流換熱器迎面風(fēng)速對(duì)空氣側(cè)換熱性能影響為參考值，將其與迎風(fēng)面積5個(gè)工況下的壓降與換熱系數(shù)實(shí)驗(yàn)值進(jìn)行對(duì)比與誤差分析.再通過(guò)正交實(shí)驗(yàn)對(duì)平行流換熱器的結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)選，以期保證空氣側(cè)的換熱性能達(dá)到最優(yōu)，可得到以下結(jié)論：

（1）選取翅片間距、翅片長(zhǎng)度、翅片高度和百葉窗角度為影響換熱性能的主要因素，每個(gè)因素選取5個(gè)水平，采用正交實(shí)驗(yàn)法得出25種不同的實(shí)驗(yàn)組合，建立了25種不同的計(jì)算模型，采用fluent軟件對(duì)其評(píng)價(jià)指標(biāo)摩擦因子和換熱系數(shù)進(jìn)行分析研究.

（2）以摩擦因子和換熱系數(shù)為換熱性能的評(píng)價(jià)指標(biāo)，對(duì)正交實(shí)驗(yàn)采用極差評(píng)價(jià)法進(jìn)行評(píng)價(jià)，得出影響摩擦因子的結(jié)構(gòu)參數(shù)的順序?yàn)槌崞g距A（38.87％）＞百葉窗角度D（25.18％）＞翅片高度C（20.79％）＞翅片長(zhǎng)度B（15.15％），影響換熱系數(shù)的結(jié)構(gòu)參數(shù)的順序：翅片間距A（56.80％）＞翅片長(zhǎng)度B（18.05％）＞翅片高度C（13.37％）＞百葉窗角度D（11.79％），最優(yōu)的結(jié)構(gòu)參數(shù)組合為組合25，結(jié)構(gòu)參數(shù)分別為（A＝2mm，B＝30mm，C＝9mm，D＝25°）.

［1］莫羚.多元平行流冷凝器熱力性能的數(shù)值模擬［D］.長(zhǎng)沙：中南大學(xué)，2007.

MO L.Flat tube numerical simulation of the parallel flow in condenser［D］.Changsha：Central South University，2007.

［2］YANG C Y，WEBB R L.Friction pressure drop of R-12in small hydraulic diameter extruded aluminum tubes with and withoutmicro-fins［J］.Internat JHeat Mass Transf，1996，39（4）：801-809.

［3］HORUZ I.Experimental investigation of a vapor absorption refrigeration system［J］.Internat JRefrig，2004，2：10-16.

［4］HORUZ I.Vapor absorption refrigeration in road transport vehicles［J］.Energ Engin，2009，125（2）：48-58.

［5］ANKIT S.Theoretical investigation ofwaste heat recovery from an IC engine using vapor absorption refrigeration system and thermoelectric converter［J］.Heat Transf，2014，2（15）：40-46.

［6］周東一，石楚平，袁文華，等.利用汽車(chē)發(fā)動(dòng)機(jī)余熱的溴化鋰吸收式制冷研究［J］.節(jié)能，2008（6）：23-26.

ZHOU D Y，SHICP，YUANW H，etal.Research on LiBr absorption refrigeration system using automotive exhausted heat ［J］.Energ Conserv，2008（6）：23-26.

［7］周東一，石楚平，肖飚，等.汽車(chē)發(fā)動(dòng)機(jī)余熱溴化鋰吸收式制冷系統(tǒng)研究［J］.制冷空調(diào)與電力機(jī)械，2008，29（122）：25-28.

ZHOU D Y，SHICP，XIAO B，etal.Research on LiBr absorption refrigeration system using the exhausted heatofautomotive engines［J］.Refrig Air Condit Electr Pow Mach，2008，29（122）：25-28.

［8］李曉科，紀(jì)威.汽車(chē)余熱制冷研究［J］.現(xiàn)代機(jī)械，2007（6）：35-37.

LIX K，JIW.A study on LiBr absorption refrigeration with automotive exhaust heat［J］.Modern Machin，2007（6）：35-37.

［9］李小華，唐景立.利用汽車(chē)余熱的吸收式制冷系統(tǒng)的研究［J］.制冷與空調(diào)，2013，27（2）：130-136.

LIX H，TANG JL.Study on LiBr absorption refrigeration system with automotive engines'exhausted heat［J］.Refrig Air Condit，2013，27（2）：130-136.

［10］劉合心，舒水明，胡中元.實(shí)用型載貨汽車(chē)空調(diào)器［J］.流體機(jī)械，2007，35（3）：76-79.

LIU H X，SHU SM，HU Z Y.Application on air-conditioner in truck［J］.Fluid Machin，2007，35（3）：76-79.

［11］龍劍.溴化鋰吸收式制冷系統(tǒng)溶液熱交換器的傳熱性能研究［J］.制冷與空調(diào)，2006（3）：68-70.

LONG J.Research on heat transfer performance of solution heatexchanger of lithium bromide absorption on refrigeration system［J］.Refrig Air-condit，2006（3）：68-70.

［12］LIANG ZW，LIU X C，YE B Y，et al.Performance investigation of fitting algorithms in surfacemicro-topography grinding processes based on multi-dimensional fuzzy relation set［J］.Internat JAdvan Manuf Tech，2013，67（7）：2779-2798.

［13］LIANG ZW，XIE B H，LIAO S P，et al.Concentration degree prediction of AWJ grinding effectiveness based on turbulence characteristics and the improved ANFIS［J］.Internat JAdvan Manuf Tech，2015，80（5）：887-905.

［14］LIANG ZW，YE B Y，WANG Y J.Four-dimensional fuzzy relation investigation in turbulence kinetic energy distribution，surface clustermodeling［J］.Arabian JSci Engin，2014，39（1）：2339-2351.

［15］谷波，方繼華，趙鵬程.空氣側(cè)結(jié)構(gòu)對(duì)多元微通道平行流冷凝器傳熱流動(dòng)性能的影響［J］.制冷學(xué)報(bào)，2014，35（6）：9-15.

GU B，F(xiàn)ANG JH，ZHAO PC.Effects of air-side structure on heat transfer and flow performance ofmultiplemicro-channel parallel flow condenser［J］.JRefrig，2014，35（6）：9-15.

［16］侯獻(xiàn)軍，云祥，劉志恩.基于Fluent的汽車(chē)空調(diào)冷凝器翅片寬度實(shí)驗(yàn)分析［J］.武漢理工大學(xué)學(xué)報(bào)，2013，35（6）：838-841.

HOU X J，YUN X，LIU Z E.Simulation analysis of automobile air conditioning condenser fin width based on Fluent［J］.J WUT，2013，35（6）：838-841.

［17］涂小蘋(píng)，梁祥飛.微通道冷凝器空氣側(cè)性能的實(shí)驗(yàn)研究［J］.制冷與空調(diào)，2013，13（9）：48-53.

TU X P，LIANG X F.Experimental study on air-side performance of micro-channel condenser［J］.Refrig Air-condit，2013，13（9）：48-53.

［18］吳瑋，陳文單，付永健，等.電動(dòng)汽車(chē)熱泵空調(diào)冷凝器空氣側(cè)性能模擬與結(jié)構(gòu)優(yōu)化［J］.制冷與空調(diào)，2013，13（1）：88-91.

WUW，CHENW D，F(xiàn)U Y J，etal.Simulation and structure optimization of the condenser of heat pump air conditioner for electric vehicle［J］Refrig Air-condit，2013，13（1）：88-91.

［19］CAVALLINIA，COLD D，DORETTIL，etal.Two-phase frictional pressure gradientof R236ea，R134a and R410A inside multiportmini-channels［J］.Exper Therm Fluid Sci，2005，29：861-870.

［20］鄧敏鋒.多元平行流蒸發(fā)器空氣側(cè)百葉窗翅片的數(shù)值模擬［J］.建筑熱能通風(fēng)空調(diào)，2012，31（1）：69-73.

DENG M F.Numerical simulation of louvered fins ofmulti-unit parallel flow type evaporator air side［J］.Build Energ Environ，2012，31（1）：69-73.

［21］NAE H K，SUNG PH.Distribution of air-water annular flow in a header of a parallel flow heat exchanger［J］.Internat J HeatMass Transf，2008，51：977-992.

［22］方繼華，谷波，田鎮(zhèn)，等.制冷劑側(cè)結(jié)構(gòu)多元微通道平行流冷凝器傳熱與流動(dòng)性能的影響［J］.上海交通大學(xué)學(xué)報(bào)，2014，48（9）：1315-1322.

FANG JH，GU B，TIAN Z，etal.Effects of refrigerant-side structure on heat transfer and flow performance ofmultimicrochannel parallel flow condenser［J］.JShanghai Jiaotong Univ，2014，48（9）：1315-1322.

［23］陸平，陳江平，陳芝久.微通道平行流氣冷器流量分配的數(shù)值模擬［J］.應(yīng)用科學(xué)學(xué)報(bào)，2007，25（3）：317-322.

LU P，CHEN JP，CHEN Z J.Computational fluid dynamics simulation of flow distribution in a gas coolerwith micro-channel parallel flow［J］.JAppl Sci，2007，25（3）：317-322.

［24］向立平，曹小林，歐陽(yáng)琴.汽車(chē)空調(diào)平行流式冷凝器空氣側(cè)性能研究［J］.流體機(jī)械，2007，35（10）：667-670.

XIANG L P，CAO X L，OUYANG Q.Research on air-side performance of parallel flow condenser［J］.Fluid Machin，2007，35（10）：667-670.

［25］LEE J K，LEE S Y.Distribution of two-phase annular flow at header channel junctions［J］.Exper Thermal Fluid Sci，2004，28：217-222.

Flow field simulation and performance analysis of heat exchanger in air conditioning system based on orthogonal experiment

LIANG Zhong-weia，b，c，ZHOU Jun-huia，LIU Xiao-chua，b，c

（a.School of Mechanical and Electrical Engineering，b.The Guangzhou Key Laboratory for Strengthen Grinding and High Performance Machining of Metal Material，c.The Guangdong Engineering Research Centre for Manure Utilization and Intelligent Irrigation of Solar Energy，Guangzhou University，Guangzhou 510006，China）

In this research，F(xiàn)luent softwarewas employed tomake a detailed simulation and performance analysis for the air flow field of heat exchanger in air conditioning system.Based on the computational grid meshing and control functions，a detailed flow simulation and performance analysis of heat exchanger has been achieved in such representative domains as temperature field，speed field，pressure field and the flow type of header pipe，thereafter a quantitative description for air flow performance can be obtained，and an orthogonal-experiment-based mechanism demonstration concentrating on the spatial distribution and evolution progresses of performance parameters can be ensured also.This paper provides theoretical basis and data preparation for the following performance improvement of heat exchanger，and the combinatorial optimization of structural parameters in air conditioning system aswell.

orthogonal experiment；heat exchanger of air conditioning system；flow simulation；performance analysis

1671-4229（2015）06-0069-10

TH 4

A

2015－09－10；

2015－10－21

國(guó)家自然科學(xué)基金資助項(xiàng)目（51575116，51205073）；國(guó)家星火計(jì)劃資助項(xiàng)目（2013GA780063）；廣東省水利科技創(chuàng)新資助項(xiàng)目（2012－11）；廣東省高校科技創(chuàng)新重點(diǎn)資助項(xiàng)目（2013KJCX0142）.

梁忠偉（1978－），男，副教授，博士，碩士生導(dǎo)師.E-mail：lzwstalin＠126.com

【責(zé)任編輯：陳鋼】