劉春節(jié)，閔春華，齊承英，王丹，田麗亭，王進

（1.常州工學院機電工程學院，江蘇常州 213002；2.河北工業(yè)大學能源與環(huán)境工程學院，天津 300401）

渦量強度強化傳熱機理數(shù)值模擬

劉春節(jié)1，閔春華2，齊承英2，王丹2，田麗亭2，王進2

（1.常州工學院機電工程學院，江蘇常州 213002；2.河北工業(yè)大學能源與環(huán)境工程學院，天津 300401）

采用數(shù)值模擬的方法，研究了矩形通道布置矩形翼縱向渦發(fā)生器后二次流強度和渦量強度在流動方向上的變化規(guī)律，分析了渦量強度對強化換熱和阻力增加的貢獻．結果表明：渦量強度JABS在整個通道內(nèi)的變化規(guī)律與局部Nu相同，即通道入口的局部Nu較大．受縱向渦發(fā)生器的影響，局部Nu出現(xiàn)一極值后逐漸減??；與二次流強度JAxBS相比，渦量強度JABS更能反映流動與換熱之間的關系；得到了Nu和f的增加量與無量綱渦量強度之間的關系．

縱向渦發(fā)生器；渦量強度；強化傳熱；數(shù)值模擬

0 引言

通道中的二次流是一種常見的流動現(xiàn)象，是發(fā)生在垂直于主流方向上的一種伴隨流動[1]．這種伴隨運動的強度通常比主流運動小得多，但是它在工程運用中的作用卻不容忽視．研究表明，二次流在強化傳熱傳質(zhì)、血液動力學、化工工程、環(huán)境除污工程等傳輸過程有明顯的強化作用[2]．通道內(nèi)布置渦發(fā)生器后會形成二次流，縱向渦的大小、位置以及發(fā)生頻率會造成流體速度偏轉(zhuǎn)以及流場壓力畸變，詳細了解流場二次流強化換熱機理就有可能通過調(diào)整區(qū)段的幾何形狀或其它參數(shù)來控制二次流，使之達到較好的工作狀態(tài)[3]．

為了分析二次流的強化傳熱機理，文獻[4-5]中一種典型的處理方法是分析二次流的強度與換熱之間的關系．實際上，二次流只是縱向渦發(fā)生器誘導產(chǎn)生的渦量的1個分量，在強化換熱機理研究中，另2個方向往往被忽略．本文就是同時考慮到3個方向渦量，研究強化換熱機理．

1 物理模型及相關參數(shù)

1.1 物理模型

物理模型如圖1所示．矩形通道的長度(L)為300mm、寬度(W)為80mm、高度（H）為20 mm，矩形翼長（l）為20 mm，矩形翼的寬度（w）為5.5 mm，矩形翼的攻角（）為30°，矩形翼前端到通道入口的距離（S）為20mm，對稱位置的兩個矩形翼最短距離s為10mm．模擬的Re為4000、6 000、8000、10000和12000，考慮到模型的對稱性，取模型的一半進行模擬．

1.2 計算方法及邊界條件

圖1 物理模型Fig.1Physicalmodel

通道內(nèi)流體為湍流流動狀態(tài)，介質(zhì)為空氣，物性參數(shù)設定為常數(shù)．控制方程采用穩(wěn)態(tài)、三維離散格式，應用有限容積法離散控制方程，壓力與速度的耦合采用SIMPLE算法．連續(xù)性方程、動量方程和能量方程均采用二階迎風格式進行離散．采用Fluent軟件進行計算．

邊界條件設定為：入口速度與溫度恒定，出口為出口邊界邊界條件，固體邊界為無滑移邊界條件，通道底面為加熱面，給定為第二類邊界條件，其他邊界為絕熱邊界條件．

1.3 相關參數(shù)定義

2 模型驗證

為驗證模型與方法的可靠性，以光通道為研究對象，將模擬得到的Nu與理論公式進行比較，所選用于比較的公式為：Nu=0.023Re0.8Pr0.3[1+De/L0.7]．Re分別為4000和12000時，Nu的理論值分別為29.0和45.8，模擬得到的結果分別為29.1和43.2．可以看出，模擬結果與理論結果吻合較好，證明模型可靠．

圖2 主流方向上局部努賽爾數(shù)NuFig.2Distributionof localNusseltnumberalong the main flow direction

3 結果分析

3.1 局部Nu分布

圖2為沿著流動方向局部Nu的變化規(guī)律．可以看出，不同Re時，由于入口效應，通道入口段的Nu較大；Nu出現(xiàn)一極值后，Nu增加，這是受到矩形翼的作用，產(chǎn)生了縱向渦；之后，隨著渦量強度逐漸減弱，Nu逐漸減?。竺鎻臏u量強度分析縱向渦發(fā)生器的強化傳熱機理．

3.2 二次流強度和渦量強度的比較

圖3和圖4分別所示為不同Re數(shù)下，渦量強度JABS和二次流強度JAxBS沿主流方向的變化規(guī)律．可以看出，對于不同的Re，在入口段，JAxBS較小，這是因為入口段沒有縱向渦發(fā)生器的作用．流體流過矩形翼后，產(chǎn)生了二次流，使JxABS突然增加，隨后由于粘性耗散的作用，JAxBS逐漸減?。畬τ跍u量強度，當Re較高時，入口段的JABS較大，且Re越大，JABS的入口效應越明顯．這是因為受到固體壁面的作用，在固體壁面附近的速度梯度較大．之后的變化規(guī)律與JAxBS相同．比較JAxBS和JABS與Nu的變化規(guī)律，可以看出，JABS與Nu的變化規(guī)律比較一致，而JAxBS與Nu的變化規(guī)律在入口段的差別較大．這表明，JAxBS更能反映對流換熱的強度，而文獻中廣泛使用的利用二次流強度表征對流換熱強弱的做法有一定局限：二次流強度只能表征縱向渦發(fā)生器后的換熱強度，而不能表征縱向渦發(fā)生器前面的換熱強度．本文后面將討論JABS對流動與換熱的影響．

圖3 主流方向渦量強度JABS的分布Fig.3Distribution ofvorticity intensity along themain flow direction

圖4 主流方向二次流強度分布Fig.4 Distribution ofsecondary flow intensityalong themain flow direction

3.3 渦量強度對強化換熱和阻力增加的貢獻的分離

為分析渦量強度對換熱的影響，將渦量強度對換熱和阻力損失的貢獻分離出來．與光通道相比，布置縱向渦發(fā)生器通道的Nu與阻力因子f的增加量的定義分別為：Nu=Nu Nu0，f=f f0，式中下角標“0”表示光通道．

4 結論

模擬了矩形通道布置矩形翼縱向渦發(fā)生器后二次流強度和渦量強度的分布規(guī)律，分析了渦量強度對強化換熱和阻力增加的貢獻的分離，得到如下主要結論：

1）通道入口的局部Nu較大，受縱向渦發(fā)生器的影響，局部Nu出現(xiàn)一極值后逐漸減?。?/p>

2）與二次流強度JAxBS相比，渦量強度JABS在整個通道內(nèi)的變化規(guī)律與局部Nu比較一致，能反映流動與換熱之間的關系．

[1]郭小勇，趙創(chuàng)要，王良璧，等．螺旋管中二次流強度的數(shù)值研究[J]．甘肅科學學報．2011，23（3）：87-92．

[2]GuptaR，Wanchoo RK，JafarA liTRM．Laminar flow in helical coils:aparametric study[J]．Ind Eng Chem Res,2011，50（2）：1150-1157．

[3]樊洪明，李先庭，何鐘怡，等．方形截面彎管二次流數(shù)值模擬[J]．熱能動力工程，2002，17（5）：510-513．

[4]Song KW，Wang LB．Theeffectivenessofsecondary flow produced by vortex generatorsmounted on both surfacesof the fin toenhanceheat transfer in a flattubebank fin heatexchanger[J]．ASMEJournalof Heat Transfer，2013，135（4）：041902．

[5]Wang LC，Su M，HuW L，Lin ZM，etal．The characteristic temperature in the definitionofheattransfer coefficienton the fin sidesurface in tube bank fin heatexchanger[J]．Numerical Heat Transfer,Part A:Applications，2011，60（10）：848-866．

[責任編輯 田豐]

Numericalsimulation forheat transferenhancementmechanism of vorticity intensity

LIU Chun-jie1，M INChun-hua2，QICheng-ying2，WANG Dan2，TIAN Li-ting2，WANG Jin2

(1.School of Mechanical and Electronic Engineering,Changzhou Institute of Technology,Jiangsu Changzhou 213002,China; 2.School of Energy and Environmental Engineering,HebeiUniversity of Technology,Tianjin 300401,China)

The secondary flow intensity and vorticity intensity distributionsalong the flow direction in a rectangular channelarranged w ith rectangular longitudinal vortex generatorsw ere numerically studied.The contribution of vorticity intensity on heat transfer enhancementand pressured drop isanalyzed.The resultshows that the distribution of vorticity intensity in thewhole channel issim ilar to thatof localNusseltnumber,i.e.the localNusseltnumber in the inletsection of thechannelis relatively higher,and then itdecreasesgraduallyafterapeak.Comparedw ith secondary flow intensity, vorticity intensity ismore suitable for revealing themechanism of theheattransferenhancementand fluid flow.The relationship of incrementsof Nusseltnumberand drag factorand the dimensionlessvorticity intensitywasobtained.

longitudinal vortex generator;vorticity intensity;heat transferenhancement;numericalsimulation

TK 124

A

1007-2373(2014)03-0065-04

2014-03-12

國家自然科學基金（51106041）

劉春節(jié)（1972-），男（漢族），副教授，博士．