葉純杰,范俊生,潘紅良
(華東理工大學(xué) 機(jī)械與動(dòng)力工程學(xué)院,上海200237)
湍流射流沖擊表面時(shí)會(huì)在沖擊點(diǎn)附近產(chǎn)生較高的局部傳熱系數(shù)現(xiàn)象,被應(yīng)用于諸多工業(yè)場(chǎng)合.隨著工業(yè)自動(dòng)化的發(fā)展和對(duì)產(chǎn)品加工效率要求的提高,需在較為封閉的環(huán)境下對(duì)產(chǎn)品進(jìn)行自動(dòng)化加工和連續(xù)作業(yè),如帶鋼紊流酸洗以及乳化液漂洗,熱軋后帶鋼表面的快速降溫以及燃?xì)馔钙饺~片的冷卻等都是在封閉條件下對(duì)運(yùn)動(dòng)的受熱面進(jìn)行傳熱或傳質(zhì).
目前國(guó)內(nèi)外對(duì)半封閉沖擊射流下移動(dòng)平板表面?zhèn)鳠崽匦缘难芯枯^少,對(duì)傾斜射流下的移動(dòng)平板傳熱和傳質(zhì)特性尚未有研究.李少華等[1]采用數(shù)值模擬的方法對(duì)復(fù)合角度斜射流下不同排孔結(jié)構(gòu)對(duì)氣膜冷卻效率的影響進(jìn)行了研究,結(jié)果表明復(fù)合角度射流可以使靜止沖擊面受熱更為均勻;Sharif等[2]用RNGk-ε模型模擬了半封閉垂直入射流在不同雷諾數(shù)、射流高度和板速條件下平板表面?zhèn)鳠嵯禂?shù)和切應(yīng)力系數(shù)的分布;Senter等[3]對(duì)垂直入射至移動(dòng)平板的射流在不同入口雷諾數(shù)和不同平板速度與射流出口速度比下的湍動(dòng)特性做了對(duì)比和分析,研究得出當(dāng)平板運(yùn)動(dòng)速度大于射流出口速度時(shí),對(duì)平板表面的湍動(dòng)能分布影響較大;Chattopadhyay 等[4]采用大渦模擬湍流模型(LES)對(duì)三維半封閉垂直射流的流態(tài)和傳熱特性進(jìn)行了數(shù)值計(jì)算,由于其所采用的計(jì)算區(qū)域較小,因此得出的結(jié)論不具有普遍意義.
筆者采用雷諾應(yīng)力湍流模型(RSM)對(duì)二維半封閉狹縫斜射流沖擊移動(dòng)平板下的流場(chǎng)進(jìn)行了計(jì)算,研究了入射角度、平板移動(dòng)速度和運(yùn)動(dòng)方向?qū)ζ桨褰谕膭?dòng)能和板面局部努塞爾數(shù)分布的影響,對(duì)工業(yè)上相似工況下的傳熱傳質(zhì)具有一定參考價(jià)值.
二維半封閉射流模型如圖1所示,射流經(jīng)槽內(nèi)頂部正中間的狹縫噴嘴以一定角度入射至平板表面,頂面為絕熱固壁條件,底面為移動(dòng)平板,左右側(cè)為射流出口.入射噴嘴寬度D為10mm,計(jì)算區(qū)域長(zhǎng)度L=100D.Gori等[5]研究了半封閉空間內(nèi)射流入射高度對(duì)靜止平板表面?zhèn)鳠嵯禂?shù)的影響,得到最佳射流高度與射流入口直徑之比為8︰1,因此按最佳入射高度確定入射高度H為80mm.設(shè)定入口雷諾數(shù)Re為10 960,其他變量為平板移動(dòng)速度與入口流速比R(平板向右移動(dòng)為正方向)和入射角θ(入口射流中心線與x軸夾角).
圖2為二維模型的網(wǎng)格分布示意圖.分別對(duì)160×40、180×50、200×60、220×70、240×80和260×90的正四邊形結(jié)構(gòu)網(wǎng)格進(jìn)行網(wǎng)格數(shù)目影響計(jì)算結(jié)果的模擬.計(jì)算表明:當(dāng)網(wǎng)格分布達(dá)到260×90時(shí),平板運(yùn)動(dòng)方向上的局部努塞爾數(shù)分布基本上不再發(fā)生變化,因而對(duì)以下的模擬均采用260×90的網(wǎng)格.
圖1 二維半封閉射流模型圖Fig.1 Two-dimensional geometry of moving plate impinged by confined slot-jet
綜上所述,文中計(jì)算采用的噴嘴寬度D=10 mm,計(jì)算區(qū)域長(zhǎng)度L=1 000mm,入射高度H=80 mm,入口雷諾數(shù)Re=10 960,網(wǎng)格分布為260×90,平板移動(dòng)速度與入口流速比R和入射角度θ為變量.
Navier-Stokes方程能精確描述射流的流動(dòng)和傳熱問(wèn)題[6].系統(tǒng)控制方程如下:
連續(xù)方程
動(dòng)量方程
能量方程
式中:U,p和T分別為平均速度,壓力和溫度;u′和T′分別為脈動(dòng)速度和脈動(dòng)溫度;ρ,μ和Pr分別為平均流體密度、分子動(dòng)力黏性系數(shù)和普朗特?cái)?shù).本文中采用加熱空氣作為射流介質(zhì),其密度和黏度分別為1.225kg/m3和1.789 4×10-5kg/m·s.
Craft等[7]比較了k-ε模型和雷諾應(yīng)力方程對(duì)射流的模擬,發(fā)現(xiàn)雷諾應(yīng)力方程對(duì)壁面壓力應(yīng)變項(xiàng)的模擬精度更高.所以筆者選用雷諾應(yīng)力模型(RSM)來(lái)封閉流體的質(zhì)量方程、動(dòng)量方程和能量方程.RSM 在計(jì)算突擴(kuò)流動(dòng)分離區(qū)和計(jì)算湍流運(yùn)輸各向異性較強(qiáng)的流動(dòng)時(shí)優(yōu)于k-ε和RNGk-ε模型且計(jì)算量小于大渦模擬(LES),是一種非常有潛力的湍流模型.
將圖2所示的整個(gè)流場(chǎng)作為計(jì)算區(qū)域,采用有限體積法離散控制方程.其中,對(duì)流項(xiàng)采用延遲修正的Quick格式離散,其具有三階精度的截差,可產(chǎn)生比二階迎風(fēng)格式等更精確的計(jì)算結(jié)果;擴(kuò)散項(xiàng)采用二階截差的中心差分格式;壓力-速度耦合采用Simplic算法.求解方程組需給定如下邊界條件:
采用均勻的入口射速
式中:uin和vin分別為進(jìn)口處射流在x和y方向的速度分量;kin和εin分別為進(jìn)口處的湍動(dòng)能和湍動(dòng)能耗散項(xiàng);Tin為射流溫度;Uin為射流在入口的速度值.
出口選用局部單向化條件,即
頂部封閉板采用無(wú)滑移條件,即
對(duì)于底部移動(dòng)平板,壁面內(nèi)計(jì)算采用壁面函數(shù)法,壁面條件為:up=Up,vp=0.不考慮封閉板與外界的自然對(duì)流及輻射換熱,將封閉板面取為絕熱壁面,即?T/?y|y=0=0.平板面采用了等壁溫條件,取平板溫度Tp=300K.其中,uc和vc分別表示射流在封閉板處x和y方向的速度分量;up和vp則表示平板在x和y方向的速度分量;Up表示平板移動(dòng)的速度值.
由于平板移動(dòng)的斜射流方式缺少試驗(yàn)數(shù)據(jù),因此無(wú)法直接比較數(shù)值模擬與試驗(yàn)數(shù)據(jù)間的吻合程度.但在半封閉垂直射流對(duì)移動(dòng)平板流場(chǎng)分布的研究方面,Senter等[3]對(duì)入口雷諾數(shù)Re=10 600,入射高度入射寬度比H/D=8︰1,平板移動(dòng)速度(板速)和入口流速比R為0和0.5條件下的射流口軸心線速度和湍流強(qiáng)度分布做了試驗(yàn)分析,數(shù)據(jù)可靠.采用RSM 湍流模型對(duì)其進(jìn)行數(shù)值計(jì)算,得到無(wú)量綱化的射流軸心線速度(軸心線速度w與射流入口速度Uin之比)和湍流強(qiáng)度分布,與試驗(yàn)值進(jìn)行對(duì)比,分別見(jiàn)圖3和4.定義湍流強(qiáng)度為:
式中:w′和分別為射流軸心線上的湍流脈動(dòng)速度與平均速度.
從圖3和圖4中可以看出,RSM 模型對(duì)射流中心線的速度分布和移動(dòng)平板表面射流湍流強(qiáng)度分布也有很好的吻合性.因此,RSM 模型能夠很好地模擬射流沖擊移動(dòng)平板下的流動(dòng)和傳熱.
圖3 射流口軸心線速度分布模擬值和試驗(yàn)值對(duì)比Fig.3 Comparison of velocity distribution along the jet centerline between simulated and experimental data
圖4 射流口軸心線湍流強(qiáng)度分布的模擬值和試驗(yàn)值對(duì)比Fig.4 Comparison of turbulence intensity along the jet centerline between simulated and experimental data
由于上封閉板的作用,在低入射角下射流的卷吸作用使得射流在未到達(dá)底部移動(dòng)平板就向一側(cè)出口發(fā)展,因此將射入角大于70°的斜射流作為研究對(duì)象.Fitzgerald等[8]對(duì)半封閉平板射流的研究發(fā)現(xiàn),入口雷諾數(shù)的變化對(duì)流場(chǎng)結(jié)構(gòu)和近壁處湍流強(qiáng)度分布的影響較小,因此筆者只對(duì)入口雷諾數(shù)Re=10 960時(shí)的射流流場(chǎng)進(jìn)行分析.
近壁湍動(dòng)能作為表述平板近壁面速度脈動(dòng)的參數(shù),與壁面的傳熱和傳質(zhì)特性有著直接聯(lián)系.將湍動(dòng)能無(wú)量綱化處理得到k*=k/U2in.圖5為不同入射角和板速下近壁面處(y/D=7.98)湍動(dòng)能k*的分布.由于射流在平板沖擊點(diǎn)附近產(chǎn)生的速度梯度最大,速度邊界層最薄,因此在此處產(chǎn)生了湍動(dòng)能峰值.在平板靜止時(shí)(R=0),此峰值隨傾角的減小而減小.當(dāng)平板向右移動(dòng)時(shí),傾斜射流在沖擊點(diǎn)附近產(chǎn)生的湍動(dòng)能隨板速的增大先降低后大幅增高.這是因?yàn)榈退贂r(shí)平板運(yùn)動(dòng)方向與射流方向相同降低了沖擊點(diǎn)處的切應(yīng)力,而隨著平板速度的增大,平板運(yùn)動(dòng)產(chǎn)生的剪切應(yīng)力大大提高了近壁處的湍動(dòng)能.在x>|20D|的沖擊點(diǎn)以外區(qū)域,板速對(duì)湍動(dòng)能強(qiáng)度起主要作用,|R|=2時(shí)的近壁面湍流強(qiáng)度值約為平板靜止時(shí)的30倍.板速較高時(shí)在左側(cè)出口附近湍動(dòng)能出現(xiàn)峰值,這是由于向右移動(dòng)的平板在左側(cè)出口近壁處形成回流,產(chǎn)生了較高的湍動(dòng)能,向左移動(dòng)的平板在右出口位置處也出現(xiàn)了相似的峰值.
圖5 不同板速和入射角下近壁面的湍動(dòng)能分布Fig.5 Distribution of turbulent kinetic energy near the wall of moving plate at different plate velocities and jet angles
當(dāng)平板向左運(yùn)動(dòng)與入射角方向相反時(shí),傾斜射流在沖擊點(diǎn)處的湍動(dòng)能值較高.R=-1時(shí),入射角為70°和80°斜射流產(chǎn)生的湍動(dòng)能k*值最大,為0.03,大大高于垂直入射靜止平板時(shí)在沖擊點(diǎn)處的k*值,此時(shí)在沖擊點(diǎn)附近產(chǎn)生的湍動(dòng)能要高于平板靜止時(shí)的2 倍以上.對(duì)于R=-2條件下的斜射流,從圖5(a)和圖5(b)可以看出:此時(shí)湍動(dòng)能最大值的位置已偏向中心點(diǎn)左側(cè),說(shuō)明高速向左運(yùn)動(dòng)的平板使得射流沖擊點(diǎn)偏向流場(chǎng)左側(cè),因此R=-2 時(shí)在沖擊點(diǎn)的湍動(dòng)能峰值較R=-1時(shí)要小.
平板運(yùn)動(dòng)使得出口邊界處產(chǎn)生回流,為了減小出口回流對(duì)傳熱特性的影響,x方向選取為-25D~25D.局部努塞爾數(shù)是表征傳熱系數(shù)增大倍數(shù)的無(wú)量綱參數(shù),定義為:
不同板速和入射角下平板表面局部努塞爾數(shù)分布見(jiàn)圖6.平板靜止時(shí)射流在沖擊區(qū)域產(chǎn)生較高的溫度梯度,表征溫度梯度的局部努塞爾數(shù)值較高.平板運(yùn)動(dòng)引起射流的偏轉(zhuǎn)導(dǎo)致沖擊力下降,致使沖擊區(qū)域的速度和溫度邊界層增厚.因此,當(dāng)平板向右運(yùn)動(dòng)時(shí),局部努塞爾數(shù)在沖擊點(diǎn)處的峰值隨板速提高而降低,峰值對(duì)應(yīng)的位置點(diǎn)也向平板運(yùn)動(dòng)方向偏移,R=2時(shí)的峰值下降為平板靜止時(shí)的60%以下,而平板表面努塞爾數(shù)的整體分布隨板速的提高趨于均勻.Chattopadhyay 等[4]對(duì)垂直入射移動(dòng)平板的研究也得到了相同的結(jié)論.而對(duì)于沖擊點(diǎn)以外區(qū)域,板速的提高大大提高了局部努塞爾數(shù).當(dāng)R=-1時(shí),斜射流下在沖擊點(diǎn)處產(chǎn)生的局部努塞爾數(shù)峰值最高,約為R=1時(shí)的3倍,且略高于垂直入射靜止平板時(shí)的努塞爾數(shù)峰值.因此,對(duì)于連續(xù)作業(yè)的強(qiáng)化傳熱,采用一定板速與入口射速之比的斜射流能夠獲得較好的傳熱效果.
圖6 不同板速和入射角下平板表面局部努塞爾數(shù)分布Fig.6 Distribution of local Nusselt number on the surface of moving plate at different plate velocities and jet angles
沖擊點(diǎn)附近最大努塞爾數(shù)Numax隨板速的變化見(jiàn)圖7.傾斜射流下,當(dāng)板速與入射角方向相同時(shí),最大努塞爾數(shù)隨板速的提高先降低后升高.在板速較低時(shí),大入射角時(shí)產(chǎn)生的溫度邊界層較薄,最大努塞爾數(shù)值較大;而當(dāng)板速提高到R=2時(shí),板速的影響大于入射角的作用,此時(shí)各入射角下的最大努塞爾數(shù)趨于同一值.與最大湍動(dòng)能相同,當(dāng)平板運(yùn)動(dòng)方向與入射角方向相反且入射角為80°時(shí),射流在沖擊點(diǎn)處能獲得最大努塞爾數(shù).
平板表面平均努塞爾數(shù)分布見(jiàn)圖8.板速較低時(shí),垂直入射下的平均努塞爾數(shù)要稍高于斜射流下的平均值;隨著板速的提高,不同入射角下平均努塞爾數(shù)值的差距減小,此時(shí)板速使平板表面產(chǎn)生了較高的剪切流,大大提高了傳熱效率.Raju 等[9]對(duì)非封閉條件下射流對(duì)移動(dòng)平板表面?zhèn)鳠岬难芯康贸觯寒?dāng)平板移動(dòng)速度小于射流入口速度時(shí),板速對(duì)表面平均努塞爾數(shù)的分布影響不大;而板速大于入口射速時(shí),平均努塞爾數(shù)大大提高.從圖8可以看出:對(duì)于半封閉射流條件,當(dāng)板速超過(guò)入口射速時(shí),板速對(duì)平板表面?zhèn)鳠岬挠绊懼饾u顯著.而斜射流下R=-1時(shí)也能一定程度上保持較高的平均努塞爾數(shù),且受熱面越大,平板移動(dòng)產(chǎn)生的傳熱效果越顯著.
圖7 沖擊點(diǎn)附近最大努塞爾數(shù)隨板速的變化Fig.7 Variation of the maximum Nusselt number with plate velocity and jet angle
圖8 平板表面平均努塞爾數(shù)分布Fig.8 Variation of the average Nusselt number with plate velocity and jet angle
(1)板速的提高增加了平板近壁處的湍動(dòng)能,當(dāng)板速為入口射速的2倍時(shí),整個(gè)板面上的近壁湍動(dòng)能大幅提高;當(dāng)平板運(yùn)動(dòng)方向與入射角方向相反且R=-1時(shí)在沖擊點(diǎn)附近產(chǎn)生的湍動(dòng)能最高,為靜止時(shí)的2倍以上.
(2)當(dāng)板速與入射角方向相同時(shí),局部努塞爾數(shù)在沖擊點(diǎn)處的峰值隨板速升高而降低,而整體分布趨于均勻;沖擊點(diǎn)以外區(qū)域,板速的提高大大提高了局部努塞爾數(shù),當(dāng)R=-1且入射角為80°時(shí),斜射流下在沖擊點(diǎn)處產(chǎn)生的局部努塞爾數(shù)峰值最高.
(3)平板運(yùn)動(dòng)速度較低時(shí),垂直入射下的板面平均努塞爾數(shù)要高于斜射流下的平均值;隨著板速的提高,平均努塞爾數(shù)先降低后大幅升高,且不同入射角下平均努塞爾數(shù)的差距逐漸減??;此外,受熱面越大,平板移動(dòng)產(chǎn)生的傳熱效果越顯著.
(4)當(dāng)平板運(yùn)動(dòng)方向與傾斜射流方向相反且板速和射流速度大小相同時(shí),能夠在沖擊點(diǎn)附近以及整個(gè)板面上都獲得較好的紊動(dòng)和傳熱效果.
[1]李少華,彭濤,張玉,等.復(fù)合角度對(duì)稱(chēng)射流氣膜冷卻效率的數(shù)值研究[J].動(dòng)力工程,2009,29(3):236-240. LI Shaohua,PENG Tao,ZHANG Yu,etal.Numerical research on film-cooling efficiency of injection with opposite orientation compound angles[J].Journal of Power Engineering,2009,29(3):236-240.
[2]SHARIF M A R,BANERJEE A.Numerical analysis of heat transfer due to confined slot-jet impingement on a moving plate[J].Applied Thermal Engineering,2009,29(2):532-540.
[3]SENTER J,SOLLIEC C.Flow field analysis of a turbulent slot air jet impinging on a moving flat surface[J].International Journal of Heat and Fluid Flow,2007,28(4):708-719.
[4]CHATTOPADHYAY H,SAHA S K.Turbulent flow and heat transfer from a slot jet impinging on a moving plate[J].International Journal of Heat and Fluid Flow,2003,24(5):685-697.
[5]GORI F,BOSSI L.Optimal slot height in the jet cooling of a circular cylinder[J].Applied Thermal Engineering,2003,23(7):859-870.
[6]劉沛清.自由紊動(dòng)射流理論[M].北京:北京航空航天大學(xué)出版社,2008:32-33.
[7]CRAFT T J,GRAHAM L J W,LAUNDER B E.Impinging jet studies for turbulence model assessment:II.an examination of the performance of four turbulence models[J].International Journal of Heat Mass Transfer,1993,36(10):2685-2697.
[8]FITZGERALD J A,GARIMELLA S V.A study of the flow field of a confined and submerged impinging jet[J].International Journal of Heat and Mass Transfer,1997,41(8):1025-1034.
[9]RAJU K S,SCHLUNDER E.U.Heat transfer between an impinging jet and a continuously moving flat surface[J].Heat and Mass Transfer,1977,10(2):131-136.