帶三角褶邊的百葉窗翅片傳熱與阻力特性數值模擬

王秋云，孫磊，張春路

（同濟大學機械與能源工程學院，上海　201804）

帶三角褶邊的百葉窗翅片傳熱與阻力特性數值模擬

王秋云，孫磊，張春路*

（同濟大學機械與能源工程學院，上海201804）

增加三角波紋褶邊，可以加強翅片邊緣強度，使翅片邊緣不輕易變形，但需要明確其對換熱器空氣側流動傳熱的影響。本文通過數值模擬，對比研究了帶三角波紋褶邊與不帶三角波紋褶邊的百葉窗翅片管換熱器空氣側換熱及壓降特性。首先與實驗關聯式比較，驗證了模型精度。其后數值模擬結果顯示，增加三角波紋褶邊后，百葉窗翅片的換熱因子及摩擦因子均有小幅增加，綜合換熱及壓降因子變化很小。本研究可為翅片管換熱器的設計提供參考。

百葉窗翅片；傳熱；壓降；數值模擬

0　引言

以百葉窗翅片管換熱器為代表的緊湊型翅片管換熱器已經廣泛應用在汽車、動力、空調制冷等多種工業領域。帶百葉窗翅片的換熱器能夠增加換熱面積、有效遏制換熱器空氣側熱邊界層的增長、減小空氣側換熱熱阻，從而提高換熱器換熱系數、增強換熱器的整體換熱能力［1］。翅片管換熱器通常包含多組平行管，管內流動的水、制冷劑等傳熱介質與直接橫掠管排的空氣進行熱交換。通常空氣側的傳熱熱阻比制冷劑側的傳熱熱阻大5到10倍［2］，因此很多研究工作都集中在翅片管外空氣側的傳熱特性［3-4］。Wang等［1-2，5］整理了91種扁管百葉窗翅片和49種圓管百葉窗翅片換熱器的實驗數據并分別擬合得到了適用于扁管和圓管百葉窗翅片傳熱及壓降的實驗關聯式。劉建等［6-7］對各類翅片換熱及壓降的實驗關聯式進行對比分析，認為Wang等在該領域已進行了長期和全面的研究，其采用的數據范圍較為廣泛，并且數據來源及處理方法可靠，推薦優先選用其相應的實驗關聯式。周宇等［8］選用Wang等的實驗關聯式驗證其二維計算模型，模擬百葉窗翅片的傳熱和阻力特性，發現百葉窗翅片傳熱性能最好的傾角是24°。寇磊等［9］選用Wang等的實驗關聯式驗證其三維計算模型，對百葉窗翅片空氣側傳熱和壓降特性進行模擬，發現百葉窗翅片傳熱最好的傾角為27°，同時阻力隨著百葉窗傾角和間距的增大而增大。

在工程實踐中，經常發現換熱器的翅片邊緣受外力影響而導致局部變形，阻塞了空氣流道，導致空氣側流動阻力上升、換熱性能下降。而通過在翅片邊緣增加三角波紋褶邊，可以加強翅片邊緣強度，使翅片邊緣不輕易變形。但是，增加三角波紋褶邊對換熱器空氣側流動傳熱會產生多少影響，目前尚缺乏相關研究。

本文采用數值模擬方法，使用FLUENT軟件對帶三角波紋褶邊百葉窗翅片和不帶褶邊百葉窗翅片的換熱特性進行三維數值模擬，以探究增加三角波紋褶邊對百葉窗翅片傳熱和壓降特性的影響。

1　物理模型及計算參數設定

百葉窗翅片管常用于家用或商用空調換熱器［10］。本文針對圓管百葉窗翅片進行研究。

研究模型分為帶三角波紋褶邊的圓管百葉窗翅片模型和不帶褶邊的圓管百葉窗翅片模型，以探究三角波紋對換熱的影響。在工程實踐中，同一三角波褶邊尺寸會匹配不同翅片間距。同樣，同一翅片間距也會匹配不同的褶邊尺寸。故本文研究同一褶邊尺寸變翅片間距以及同一翅片間距變褶邊尺寸的百葉窗翅片傳熱特性。表1給出的是百葉窗翅片固定參數的具體尺寸。表2給出的是褶邊尺寸與翅片間距的匹配參數。

對百葉窗翅片換熱特性的數值模擬采用三維模型。如圖1所示，假定翅片間流道均勻，根據物理模型自身的對稱性和周期性［11］，選取框線區域內部區域為計算區域：以翅片中心為中心，以翅片間距大小為外側風道尺寸。為保證計算進口空氣流動狀態穩定以及出口充分發展，計算區域在進出口均進行延長，如圖2所示。

百葉窗翅片換熱特性的模擬采用商業軟件FLUENT［12］。空氣為不可壓縮粘性流體。翅片板材料為鋁箔。管壁為銅質標準管。計算模型選用帶能量方程的湍流兩方程RNG k-ε模型，壓力速度耦合采用SIMPLEC 算法，動量、湍流動能、湍流耗散率、能量均采用采用二階迎風差分格式，翅片采用內部導熱和表面對流換熱耦合的方式完成計算。

表1　百葉窗翅片固定參數的具體尺寸

表2　褶邊尺寸與翅片間距的匹配參數

圖1　計算單元

圖2　百葉窗翅片管換熱器計算區域

邊界條件設置如下。

1）入口邊界：以翅片外側空氣通道X負方向上端面截面作為空氣入口，采用速度進口邊界條件，同時給定空氣進口溫度（308 K）。

2）出口邊界：以翅片外側空氣通道x正方向上端面截面作為空氣出口，采用壓力出口邊界條件。

3）圓管外表面溫度恒定（330 K），采用固定溫度壁面條件。

4）翅片板及百葉窗板為固體壁面，傳熱耦合。

計算區域Z方向上下平面為周期性邊界，Y方向前后平面為對稱邊界。

2　計算結果與分析

2.1網格獨立性檢驗

對計算區域進行分區、劃分網格并進行網格獨立性驗證。取百葉窗翅片換熱系數（j因子）與摩擦系數（f因子）兩項作為驗證指標。如圖3和圖4所示，當網格數大于300萬時，計算得到的j因子數值不再變化。網格數為158萬時與網格數為107萬時的j因子計算數值相差為1.5%，網格數為319萬時與網格數為158萬時的j因子計算數值相差僅為0.2%，而f因子基本不變。可以認為，當網格數大于158萬時，計算結果已經不受網格數量影響。

圖3　不同網格數下換熱因子計算結果

圖4　不同網格數下摩擦因子的計算結果

2.2模型驗證

如圖3所示，模擬計算的j因子和f因子與Wang等［2］的實驗關聯式經驗值相比，變化趨勢基本相同，均是隨著迎面風速的增加而減小，且變化率逐漸變小。計算得到的j因子與經驗值平均誤差為1.3%，最大誤差為10.4%；f因子與經驗值平均誤差為4.1%，最大誤差為7.6%。由此可見，本文的計算模型精度可靠。

圖3　換熱/摩擦因子計算結果與實驗關聯式［2］的比較

2.3計算結果及分析

本文首先選擇目前已經投入生產使用的一種三角褶邊尺寸，固定褶邊尺寸，研究其在不同翅片間距、不同迎面風速時對百葉窗翅片換熱及壓降的影響。對翅片間距為1.6 mm、1.8 mm、2.0 mm、2.2 mm、2.4 mm、2.6 mm，迎面風速為1 m/s～6 m/s的百葉窗翅片換熱及壓降模擬。

2.3.1三角褶邊對百葉窗翅片管換熱器換熱及壓降的影響

圖4是不帶三角波紋褶邊與帶三角波紋褶邊的百葉窗翅片采用不同翅片間距和不同風速時的換熱因子和摩擦因子。其中（a）、（b）給出的是在不同迎面風速下的換熱因子j隨翅片間距的變化。可以看出，無論是否帶褶邊，相同翅片間距下，隨著迎面風速的增加，百葉窗翅片的流動換熱能力減弱，j因子變小。相同迎面風速下，隨著翅片間距增加，j因子有減小的趨勢，但減小的速率非常緩慢。可以認為相同迎面風速時，翅片間距的改變對百葉窗翅片換熱的影響較小。圖6（c）、（d）是摩擦因子f隨翅片間距的變化。可以看出，無論是否帶褶邊，相同翅片間距，隨著迎面風速的增加，f因子減小。當翅片間距大時，曲線較平緩，f因子減小的速率較小；翅片間距小時，曲線相對較陡，f因子減小的速率較快。相同迎面風速下，f因子隨翅片間距變大而減小，且隨著迎面風速的增加，變化程度逐漸減小。

圖4　不同形態翅片在不同翅片間距和迎面風速下的換熱因子和摩擦因子變化

圖5對比了帶三角波紋褶邊與不帶三角波紋褶邊的百葉窗翅片換熱及壓降的變化情況。圖5（a）表示相應翅片間距下，帶三角波紋褶邊的百葉窗翅片換熱因子j1相對于不帶三角波紋褶邊的百葉窗翅片換熱因子j0的增量。除了低風速、小翅片間距情況下影響較大以外，其它情況下的影響都不大，整體上帶褶邊的百葉窗翅片的j因子有0.36%左右的增加。圖5（b）表示相應翅片間距下，帶三角波紋褶邊的百葉窗翅片摩擦因子f1相對于不帶三角波紋褶邊的百葉窗翅片摩擦因子f0的增量。沒有突出的變化，整體上帶褶邊的百葉窗翅片f因子有0.55%左右的增加。

此外，統計84個帶三角波紋褶邊的百葉窗翅片換熱及壓降模擬算例中翅片進出口褶邊段與翅片中間段各占總壓降的比例，可以發現，中間段壓降均占總壓降的95%以上。進口褶邊段壓降均占總壓降4%以下。由此可見，影響帶褶邊的百葉窗翅片流動壓降的主要部分仍然是翅片中間段。三角波紋褶邊因其尺寸較小對翅片整體壓降影響不大。

采用j/f因子綜合考慮換熱器翅片流動換熱情況，其意義在于判斷通過使用強化換熱技術后，相同流量下換熱能力的增加是否大于阻力的增加。對比帶三角波紋褶邊與不帶三角波紋褶邊的百葉窗翅片換熱壓降綜合評價因子j/f，增加褶邊后，j/f因子略有減小，大多為不帶褶邊j/f因子的0.998左右。如圖6所示，j1/f1為帶三角波紋褶邊的百葉窗翅片綜合換熱因子，j0/f0為不帶三角波紋褶邊的百葉窗翅片綜合換熱因子。以翅片間距2 mm為例說明該綜合性能評價指標變化的趨勢。可以認為，增加褶邊后，綜合了壓降的換熱效果與不帶褶邊時幾乎相同。

圖5　不同翅片間距下帶褶邊比不帶褶邊的百葉窗翅片換熱/摩擦因子增加百分比

圖6　翅片間距為2 mm時不同迎面風速下百葉窗翅片的綜合性能因子之比

2.3.2三角褶邊尺寸對百葉窗翅片換熱及壓降的影響

如上文所述，采用該特定尺寸三角波紋褶邊，在增加百葉窗翅片強度的同時，對其換熱及壓降性能影響不大。分析三角波紋褶邊的幾何尺寸，我們發現：三角波紋褶邊的高度、寬度、深度、以及其與百葉窗翅片上百葉窗板的高度等尺寸的匹配性都可能影響百葉窗翅片的流動換熱及壓降性能。在此我們僅研究三角波紋褶邊的高度對百葉窗翅片換熱及壓降特性的影響。

選擇工程實踐中常見的翅片間距2.0 mm的百葉窗翅片管換熱器作為研究對象。對褶邊高度為0.5 mm、0.8 mm、1.0 mm、1.5 mm時，迎面風速為2～5 m/s分別進行換熱及壓降的數值模擬。下圖均以迎面風速為2 m/s為例進行說明：

對于定翅片間距，改變褶邊高度的算例，當褶邊高度增加，換熱因子增加，同時壓降增加，但因相同流量下阻力的增加大于換熱能力的增加導致綜合換熱因子j/f下降。但是，這個下降的幅度是很小的。以翅片間距為2.0 mm，迎面風速為2 m/s時的情形為例，如下圖所示，褶邊高度減小，j/f因子不斷減小，但1.5 mm時的j/f因子只比0.5 mm時的減小了0.67%。

圖7　2.0 mm翅片間距不同褶邊高度下綜合傳熱因子

3　結論

本文對帶三角波紋褶邊的百葉窗翅片與不帶三角波紋褶邊的百葉窗翅片換熱及壓降特性進行了研究，通過大量數值模擬計算與分析，獲得以下主要結論。

1）相同翅片間距、相同迎面風速下，帶褶邊的百葉窗翅片與不帶褶邊的百葉窗翅片相比，f因子有0.55%左右的增加，j因子有0.36%左右的增加，綜合換熱因子j/f略有減小，大多為不帶褶邊j/f因子的0.998左右。三角波紋褶邊的增加對百葉窗翅片的換熱和壓降的影響不大，工程中可以忽略。

2）統計進出口褶邊段與翅片中間段各占總壓降的比例，中間段均占總壓降的95%以上。進口褶邊段壓降均占總壓降4%以下。增加了三角波紋褶邊的百葉窗翅片壓降主要來源仍為翅片中段。褶邊因其尺寸較小對翅片整體壓降影響不大。

3）定翅片間距，當褶邊高度增加，換熱效果增加，同時壓降增加，j/f因子下降，但其變化幅度均較小。

綜上所述，在百葉窗翅片邊緣增加三角褶邊對與百葉窗翅片本身的換熱及壓降影響程度不大。從工程實踐的需要出發，在翅片邊緣增加三角波紋褶邊以增加其強度的做法是可以提倡的。

［1］CHANG Y J，WANG C C.A generalized heat transfer correlation for Iouver fin geometry［J］.International Journal of Heat and Mass Transfer，1997，40（3）∶533-44.

［2］WANG C C，LEE C J，CHANG C T，et al.Heat transfer and friction correlation for compact louvered fin-and-tube heat exchangers［J］.International Journal of Heat and Mass Transfer，1998，42（11）∶1945-56.

［3］SUGA K，AOKI H.Numerical study on heat transfer and pressure drop in multilouvered fins［J］.Journal of Enhanced Heat Transfer，1995，2（3）∶231-238.

［4］KIM M H，BULLARD C W.Air-side thermal hydraulic performance of multi-louvered fin aluminum heat exchangers［J］.International Journal of Refrigeration，2002，25（3）∶390-400.

［5］CHANG Y J，HSU K C，LIN Y T，et al.A generalized friction correlation for louver fin geometry［J］.International Journal of Heat and Mass Transfer，2000，43（12）∶2237-2243.

［6］劉建，魏文建，丁國良，等.翅片管式換熱器換熱與壓降特性的實驗研究進展——關聯式［J］.制冷學報，2003，24（4）∶21-27.

［7］劉建，魏文建，丁國良，等.翅片管式換熱器換熱與壓降特性的實驗研究進展——實驗研究［J］.制冷學報，2003，24（3）∶25-30.

［8］周宇，李惠珍，陶文銓.百葉窗翅片的傳熱與阻力特性的數值模擬［C］.中國動力工程學會第三屆青年學術年會論文集，2005.

［9］寇磊，廖勝明，劉玉涵.百葉窗翅片傳熱特性的數值模擬［J］.建筑熱能通風空調，2009，28（1）∶6-9.

［10］WANG C C，CHI K Y，CHANG Y J，et al.An experimental study of heat transfer and friction characteristics of typical louver fin-and-tube heat exchangers［J］.International Journal of Heat and Mass Transfer，1998，41（4）∶817-22.

［11］張智，涂旺榮，韓蔚，等.空調用冷凝器中的空氣流動與傳熱分析［J］.制冷技術，2002，22（1）∶8-13.

［12］劉睿盈，吳靜怡，孫曉琳.直膨式太陽能熱泵熱水器集熱/蒸發器流道結構分析與實驗［J］.制冷技術，2014，34（2）∶1-6.

Numerical Simulation on Heat Transfer and Pressure Drop Characteristics of Louvered Fins with Triangle Corrugated Ruffles

WANG Qiu-yun，SUN Lei，ZHANG Chun-lu*
（School of Mechanical Engineering，Tongji University，Shanghai 201804，China）

The louvered fin can be strengthened and protected from deformation，if its edges were added triangle corrugated ruffles.However，the impact of the ruffles on the heat transfer and pressure drop of the heat exchanger should be clear.Heat transfer and pressure drop of louvered fins（with or without ruffles）were numerically studied with FLUENT.The precisions of the numerical models were firstly validated by experimental correlations.The simulation results showed that both the heat transfer factor j and the friction factor f of louvered fins with the triangle corrugated ruffles have small increase in comparison with those of the louvered fins without ruffles.The comprehensive index of heat transfer and pressure drop changes little.The study provides a valuable reference for the design of the fin-and-tube heat exchanger.

Louvered fin；Heat transfer；Pressure drop；Numerical simulation

10.3969/j.issn.2095-4468.2015.03.105

*張春路（1971-），男，教授，博士。研究方向：制冷空調系統仿真、優化與控制。聯系地址：同濟大學機械與能源工程學院。郵編：201804。聯系電話：13671825133。E-mail∶chunlu.zhang@gmail.com。

本論文選自2014 年第八屆全國制冷空調新技術研討會。