王小雨，張聿祥

(江漢大學(xué) 智能制造學(xué)院，湖北 武漢 430056)

能源問(wèn)題不僅關(guān)系到人民生活質(zhì)量和生態(tài)環(huán)境發(fā)展，還事關(guān)國(guó)家經(jīng)濟(jì)安全和戰(zhàn)略安全，其早已成為經(jīng)濟(jì)社會(huì)發(fā)展的重要問(wèn)題之一[1]。如何提高能源的利用已成為至關(guān)重要的問(wèn)題，人類對(duì)能源的利用過(guò)程絕大部分都涉及到熱能的傳遞和轉(zhuǎn)化。換熱器是進(jìn)行熱交換的設(shè)備，提高換熱器的換熱性能對(duì)節(jié)能具有重要意義。換熱管是管殼式換熱器的基本組成單元，是實(shí)現(xiàn)換熱過(guò)程的關(guān)鍵組件。為了提高換熱器的換熱性能，許多學(xué)者對(duì)換熱管內(nèi)流體的流動(dòng)和換熱性能展開(kāi)了研究，并形成了豐富的強(qiáng)化傳熱技術(shù)體系。根據(jù)是否需要提供除泵功之外的附加動(dòng)力，將強(qiáng)化傳熱技術(shù)分為主動(dòng)強(qiáng)化、被動(dòng)強(qiáng)化以及兩者相結(jié)合的復(fù)合強(qiáng)化[2]。主動(dòng)強(qiáng)化技術(shù)因需要增加設(shè)備，使系統(tǒng)變得更加復(fù)雜且投資成本增加，在工業(yè)中的使用受到了限制。被動(dòng)強(qiáng)化傳熱技術(shù)旨在修改傳熱表面的結(jié)構(gòu)或者在流體中增加擾流元件來(lái)擾動(dòng)流體和提高傳熱性能，因其系統(tǒng)簡(jiǎn)單、制造成本低、可靠性高等優(yōu)點(diǎn)在工業(yè)中具有更廣泛的應(yīng)用。

被動(dòng)強(qiáng)化傳熱技術(shù)主要從三個(gè)方面展開(kāi)，即管型改變[3-4]、管壁處理[5-7]和內(nèi)插擾流元件[8-13]。其中內(nèi)插擾流元件是一種常見(jiàn)的管內(nèi)被動(dòng)式強(qiáng)化傳熱技術(shù)，該技術(shù)通過(guò)在管內(nèi)插入不同結(jié)構(gòu)的擾流元件來(lái)擾動(dòng)流體，改變通道內(nèi)流體的流動(dòng)結(jié)構(gòu)，促使管內(nèi)形成縱向渦或二次流，從而減薄或破壞熱邊界層[14-17]。相比其他傳熱技術(shù)，該技術(shù)具有加工簡(jiǎn)單，制造成本低，可直接對(duì)舊設(shè)備進(jìn)行改進(jìn)，且易于維護(hù)等優(yōu)點(diǎn)。常見(jiàn)的擾流元件有扭帶[8]、螺旋片[9]、渦發(fā)生器[10]、多孔介質(zhì)[11]、小管插入物[12]、波浪帶插入物[13]等，它們各自有相應(yīng)的優(yōu)缺點(diǎn)，適用于不同的場(chǎng)合。

周璐瑤等[18-19]提出了一種縱向渦流強(qiáng)化管(內(nèi)插中心斜桿)的換熱管，研究了斜桿長(zhǎng)度、斜桿傾角和桿長(zhǎng)對(duì)換熱性能的影響。通過(guò)模擬發(fā)現(xiàn)這種類型的管插入件可以使流體形成具有多個(gè)縱向渦流的渦流結(jié)構(gòu)，并且可以在阻力增加不太大的前提下，有效提高換熱能力。通過(guò)內(nèi)插中心斜桿將管內(nèi)流體分為6個(gè)周向?qū)ΨQ的流動(dòng)區(qū)域，在每個(gè)區(qū)域中形成了一個(gè)縱向旋流的渦結(jié)構(gòu)，中心斜桿使管內(nèi)形成了多縱向旋流的流動(dòng)結(jié)構(gòu)。因此，內(nèi)插中心斜桿的換熱管能夠獲得優(yōu)異的強(qiáng)化傳熱性能。

但內(nèi)插中心斜桿的數(shù)目和直徑等參數(shù)對(duì)換熱管的換熱性能的影響，尚未有學(xué)者進(jìn)行研究，因此，本文討論了斜桿數(shù)目、節(jié)距、直徑等參數(shù)對(duì)內(nèi)插中心斜桿換熱管換熱性能的影響，以期進(jìn)一步優(yōu)化內(nèi)插中心斜桿的幾何結(jié)構(gòu)，對(duì)內(nèi)插中心斜桿換熱管的換熱性能研究起到一定的指導(dǎo)意義。

1 光管的數(shù)值模擬及驗(yàn)證

圖1 光管的網(wǎng)格劃分Fig.1 Mesh generation of the smooth tube

1.1 光管的參數(shù)設(shè)置

本文所選管內(nèi)介質(zhì)為水，20 ℃的物性參數(shù)為：ρ=998.2 kg/m3，μ=1.003×10-3kg/(m·s)，CP=4 182 J/(kg·K)和λ= 0.6 W/(m·K)；光管管長(zhǎng)為L(zhǎng)=500 mm，管內(nèi)徑為di=20 mm，網(wǎng)格劃分如圖1所示。

1.2 控制方程與假設(shè)

換熱管的模擬要滿足3個(gè)基本方程：質(zhì)量守恒方程即連續(xù)性方程、能量守恒方程和動(dòng)量守恒方程。具體公式參考文獻(xiàn)[20]。

為了便于模型的建立和計(jì)算，做出以下簡(jiǎn)化：(1)忽略中心傾斜桿插入元件的熱傳導(dǎo)和換熱管壁厚；(2)流體的物理參數(shù)為固定值;(3)流體是一種連續(xù)不可壓縮、不考慮重力、沒(méi)有黏性耗散、各向同性牛頓流體；(4)熱輻射和重力的影響予以忽略。

1.3 邊界條件與計(jì)算方法

采用Fluent進(jìn)行計(jì)算時(shí)，選擇3D、穩(wěn)態(tài)和基于壓力的隱式求解器。管的內(nèi)壁采用2 000 W/m2的恒定熱通量密度，并定義為無(wú)滑壁面；中心斜桿擾流元件的表面為絕熱無(wú)滑壁面。在管道出口處采用自由出流邊界。湍流模型選擇標(biāo)準(zhǔn)的κ-ε模型，計(jì)算時(shí)采用SIMPLE算法，選擇基于最小二乘單元的梯度，壓力選擇二階格式，動(dòng)量和能量采用二階迎風(fēng)格式，湍流動(dòng)能和湍流耗散選擇一階迎風(fēng)格式。

1.4 網(wǎng)格無(wú)關(guān)性驗(yàn)證

使用Fluent軟件數(shù)值模擬計(jì)算出努塞爾數(shù)與網(wǎng)格數(shù)的關(guān)系。如圖2所示，在一定數(shù)量網(wǎng)格內(nèi)努塞爾數(shù)變化波動(dòng)較大，但隨著網(wǎng)格數(shù)量的增加努塞爾數(shù)趨于平穩(wěn)，模擬結(jié)果顯示在五十萬(wàn)的網(wǎng)格數(shù)后努塞爾數(shù)變化明顯較小，因此在后續(xù)的數(shù)值模擬中，出于電腦計(jì)算效率與能力問(wèn)題，最終選擇五十萬(wàn)以上的網(wǎng)格數(shù)量。

圖2 網(wǎng)格數(shù)量與努塞爾數(shù)的關(guān)系Fig.2 Relation between the mesh number and Nusselt number

1.5 光管的經(jīng)驗(yàn)公式計(jì)算值與模擬對(duì)比

為了驗(yàn)證數(shù)值計(jì)算結(jié)果的可靠性，將光管模擬得出的努塞爾數(shù)和壓降與經(jīng)驗(yàn)公式計(jì)算的值進(jìn)行比較。

因?yàn)楣艿赖闹睆皆O(shè)置得比較小，同時(shí)流體與壁表面之間的溫差也比較小，當(dāng)流體的μ/p值比較大時(shí)，自然對(duì)流對(duì)強(qiáng)制層流傳熱的影響可以忽略不計(jì)。此時(shí)，對(duì)流傳熱系數(shù)可以使用西德?tīng)?Sieder)和塔特(Tate)關(guān)聯(lián)式(1)計(jì)算：

(1)

壓降經(jīng)驗(yàn)公式如下

(2)

式中：Δpf為因克服阻力而引起的壓力降，Pa；hf為能量損失，J/kg；λ為導(dǎo)熱系數(shù)，W/(m·℃) ；u為流體流速，m/s。

在介質(zhì)流速為0.07、0.08、0.09、0.10 m/s時(shí)，經(jīng)驗(yàn)公式計(jì)算的努塞爾數(shù)和壓降與模擬值的比較分別如圖3和圖4所示。

圖3 光管努塞爾數(shù)的模擬值與經(jīng)驗(yàn)公式計(jì)算值對(duì)比Fig.3 Comparison between simulated and calculated values of empirical formula on the Nusselt number of the smooth tube

圖4 光管壓降的模擬值與經(jīng)驗(yàn)公式計(jì)算值對(duì)比Fig.4 Comparison between simulated and calculated values of empirical formula on the pressure drop of the smooth tube

2 內(nèi)插中心斜桿換熱管的模擬

內(nèi)插中心斜桿的換熱管是一種縱向渦流強(qiáng)化管，可以明顯改善管內(nèi)熱量的傳遞。為了研究不同結(jié)構(gòu)尺寸的中心斜桿換熱管對(duì)換熱性能的影響，本文選取了中心斜桿的三個(gè)幾何參數(shù)，即傾斜桿的數(shù)量(n=2、3、4、5、6)，節(jié)距(p=20、25、30、35、40 mm)和連桿的直徑(d=1.0、1.5、2.0、2.5、3.0 mm)，作為研究對(duì)象，分別探討了各參數(shù)對(duì)努塞爾數(shù)和壓降的影響。

2.1 中心斜桿幾何建模和網(wǎng)格劃分

本文設(shè)計(jì)的中心斜桿的幾何結(jié)構(gòu)如圖5所示。傾斜的中心連桿由直徑為2 mm的小連桿組成，這些連桿沿圓周對(duì)稱地傾斜，傾斜角為30°并以相等的間隔布置在中心2 mm直徑的連桿上。在模擬中，將換熱管的幾何參數(shù)選擇為管的長(zhǎng)度L=500 mm，管的內(nèi)徑di= 20 mm。

圖5 中心斜桿的幾何模型Fig.5 Mode of the central inclined rod

2.2 結(jié)果分析

通過(guò)Fluent后處理，得出內(nèi)插中心斜桿換熱管圓截面溫度分布云圖(圖6)，可以看出中心斜桿對(duì)于流體的混合熱交換有很大的幫助，因?yàn)榱黧w受到中心斜桿的阻礙后，被引導(dǎo)后有一個(gè)徑向的矢量，會(huì)與靠近管壁的高溫流體層進(jìn)行沖刷混合，從圖中比較容易看到高溫層明顯變薄。圖7為管壁溫度分布，可以看出中心斜桿的擾流對(duì)打破邊界層有一定作用。

圖6 內(nèi)插中心斜桿換熱管圓截面溫度分布云圖Fig.6 Temperature field of the central inclined rod tube in the circular section

圖7 內(nèi)插中心斜桿換熱管速度分布云圖Fig.7 Velocity field of the central inclined rod tube

圖8為內(nèi)插中心斜桿換熱管圓截面速度分布云圖，可以看出因中心斜桿的擾流產(chǎn)生的流體軌跡，十分符合文獻(xiàn)[14]所展示的狀態(tài)，靠近斜桿的區(qū)域速度相對(duì)較高，向其他區(qū)域有輻射。偏紅色區(qū)域就是因?yàn)橹行男睏U的結(jié)構(gòu)對(duì)流體產(chǎn)生了旋流，這對(duì)于冷、熱流體混合有很大幫助。從圖9可以看到在斜桿后還是存在有一定的運(yùn)動(dòng)死區(qū)。由此可以初步推斷如果過(guò)多地設(shè)置斜桿數(shù)量，或者是過(guò)大地設(shè)置中心連桿上的斜桿直徑會(huì)使這種運(yùn)動(dòng)死區(qū)增大，對(duì)于流體流動(dòng)以及冷熱流體混合都有負(fù)面影響。

圖8 內(nèi)插中心斜桿換熱管圓截面速度分布云圖Fig.8 Velocity field of the central inclined rod tube in the circular section

圖9 內(nèi)插中心斜桿換熱管速度分布云圖Fig.9 Velocity field of the inserted central inclined rod tube

3 內(nèi)插中心斜桿換熱管幾何結(jié)構(gòu)的研究

3.1 斜桿數(shù)量對(duì)內(nèi)插中心斜桿換熱管性能的影響

對(duì)改變斜桿數(shù)量后換熱管的換熱性能進(jìn)行分析，得到了努塞爾數(shù)和壓降隨斜桿數(shù)量的變化關(guān)系，如圖10所示。從圖中可以看出，在2到5個(gè)斜桿時(shí)，換熱系數(shù)是不斷增大的，這是由于通過(guò)斜桿進(jìn)行擾流，使得流體在徑向空間發(fā)生冷熱流體的混合，更好地進(jìn)行熱交換；但是當(dāng)斜桿數(shù)量增加到6時(shí)，這時(shí)的努塞爾數(shù)開(kāi)始下降，結(jié)合其壓降的激增，可以推斷是因?yàn)樽钃跷镞^(guò)多，而形成了較多的流動(dòng)死區(qū)，導(dǎo)致?lián)Q熱系數(shù)無(wú)法進(jìn)一步提高，甚至開(kāi)始下降。因此為了提高換熱系數(shù)不能一味地增加其斜桿數(shù)量，應(yīng)選擇適當(dāng)?shù)臄?shù)量。從圖10中也可以看出，隨著斜桿數(shù)量的增加，其進(jìn)出口端面的壓降也一直在增加，說(shuō)明增加斜桿數(shù)量雖然可以一定程度增大換熱系數(shù)，但其帶來(lái)阻力的負(fù)面影響也在隨之增大。

內(nèi)插中心斜桿換熱管的換熱性能主要受到努塞爾數(shù)和壓降的影響，為了綜合評(píng)價(jià)兩者對(duì)換熱管綜合換熱性能的影響，采用了文獻(xiàn)中廣泛認(rèn)可的性能評(píng)價(jià)指標(biāo)IPEC[21-22]來(lái)研究二者的影響，計(jì)算公式如式(3)所示：

(3)

其中，Nu為內(nèi)插中心斜桿換熱管努塞爾數(shù)，Nus為光管努塞爾數(shù)，f為內(nèi)插中心斜桿換熱管阻力系數(shù)，fs為光管阻力系數(shù)。

圖11為斜桿數(shù)量對(duì)內(nèi)插中心斜桿換熱管綜合換熱性能的影響，由圖11可知，隨著斜桿數(shù)量的增加，綜合換熱性能下降，但是內(nèi)插中心斜桿換熱管的綜合換熱性能始終大于光管的換熱性能(IPEC>1)。本文斜桿數(shù)量為3時(shí)，內(nèi)插中心斜桿換熱管的綜合換熱性能較好。

圖10 斜桿數(shù)量對(duì)內(nèi)插中心斜桿換熱管換熱性能的影響Fig.10 Variation in heat transfer performance with the number of inclined rods

圖11 斜桿數(shù)量對(duì)內(nèi)插中心斜桿換熱管的綜合換熱性能的影響Fig.11 Variation in IPEC value with the number of inclined rods

3.2 斜桿節(jié)距對(duì)內(nèi)插中心斜桿換熱性能的影響

對(duì)改變斜桿節(jié)距后換熱管的換熱性能進(jìn)行分析，得到了努塞爾數(shù)和壓降隨斜桿數(shù)量的變化關(guān)系。從圖12中可以看出，隨著斜桿節(jié)距的增加，其努塞爾數(shù)和壓降隨之減小。這是因?yàn)樾睏U節(jié)距增加后，斜桿對(duì)流體的擾流作用降低，而流體的擾流會(huì)有助于溫度的傳遞。因?yàn)樾睏U節(jié)距的增加，使得流動(dòng)阻力減小，故而壓降降低。

圖13為斜桿節(jié)距對(duì)內(nèi)插中心斜桿換熱管綜合換熱性能的影響，由圖13可知，隨著斜桿節(jié)距的增加，綜合換熱性能下降，但是內(nèi)插中心斜桿換熱管的綜合換熱性能始終大于光管的換熱性能(IPEC>1)。本文斜桿節(jié)距為20 mm時(shí)，內(nèi)插中心斜桿換熱管的綜合換熱性能較好。

圖12 斜桿節(jié)距對(duì)內(nèi)插中心斜桿換熱管換熱性能的影響Fig.12 Variation in heat transfer performance with the pitch of inclined rods

圖13 斜桿節(jié)距對(duì)內(nèi)插中心斜桿換熱管的綜合換熱性能的影響Fig.13 Variation in IPEC value with the pitch of inclined rods

3.3 斜桿直徑對(duì)內(nèi)插中心斜桿換熱性能的影響

對(duì)改變斜桿直徑后換熱管的換熱性能進(jìn)行分析，得到了努塞爾數(shù)和壓降隨斜桿數(shù)量的變化關(guān)系，如圖14所示。從圖14可以看出，增大中心斜桿的直徑對(duì)于傳熱有一定幫助，努塞爾數(shù)隨著斜桿直徑的增大也逐漸變大，特別是從1.0 mm到1.5 mm時(shí)其傳熱性能有一個(gè)突變的激增，而從1.5 mm到3.0 mm時(shí)變化較為平緩。又因?yàn)殡S著斜桿直徑的增大，其帶來(lái)的壓降也隨之增加。故斜桿直徑不宜過(guò)大或過(guò)小。

圖15為斜桿直徑對(duì)內(nèi)插中心斜桿換熱管綜合換熱性能的影響，由圖15可知，隨著斜桿直徑的增加，綜合換熱性能先上升后下降。在斜桿直徑為1.0 mm時(shí)，內(nèi)插中心斜桿換熱管的綜合換熱性能小于光管的換熱性能(IPEC<1)。但當(dāng)斜桿直徑大于等于1.5 mm時(shí)，內(nèi)插中心斜桿換熱管的綜合換熱性能大于光管的換熱性能(IPEC>1)。本文斜桿直徑為2.0 mm時(shí)，內(nèi)插中心斜桿換熱管的綜合換熱性能較好。

圖14 斜桿直徑對(duì)內(nèi)插中心斜桿換熱管換熱性能的影響Fig.14 Variation in heat transfer performance with the diameter of inclined rods

圖15 斜桿直徑對(duì)內(nèi)插中心斜桿換熱管的綜合換熱性能的影響Fig.15 Variation in IPEC value with the diameter of inclined rods

4 結(jié)論

本文主要探討了斜桿幾何機(jī)構(gòu)對(duì)內(nèi)插中心斜桿換熱管換熱性能的影響，根據(jù)模擬結(jié)果得出以下結(jié)論：

(1)通過(guò)綜合換熱性能的比較，發(fā)現(xiàn)除斜桿直徑為1.0 mm時(shí)，內(nèi)插中心斜桿換熱管的綜合換熱性能略低于光管(IPEC<1)，其余參數(shù)下的內(nèi)插中心斜桿換熱管的綜合換熱性能都優(yōu)于光管(IPEC>1)，且綜合換熱性能較光管提高了50%以上。

(2)在斜桿數(shù)量為2到5時(shí)，其努塞爾數(shù)是不斷增大的，但是當(dāng)其斜桿數(shù)量增加到6時(shí)，這時(shí)的努塞爾數(shù)開(kāi)始下降。本文通過(guò)對(duì)內(nèi)插中心換熱管綜合性能的分析得出斜桿數(shù)量為3時(shí)，內(nèi)插中心斜桿換熱管的綜合換熱性能較好。

(3)隨著斜桿節(jié)距的增加，其努塞爾數(shù)隨之減小，流體阻力減小，壓降降低。對(duì)內(nèi)插中心斜桿換熱管斜桿節(jié)距的綜合換熱性能分析得出，斜桿節(jié)距為20 mm時(shí)，內(nèi)插中心斜桿換熱管的綜合換熱性能較好。

(4)增大中心斜桿的直徑對(duì)傳熱有一定幫助，努塞爾數(shù)隨著斜桿直徑的增大也逐漸變大，特別是從1.0 mm到1.5 mm時(shí)其傳熱性能有一個(gè)突變的激增，而從1.5 mm到3.0 mm時(shí)變化較為平緩。故斜桿直徑不宜過(guò)大或過(guò)小，本文對(duì)內(nèi)插中心斜桿換熱管斜桿直徑進(jìn)行綜合換熱性能分析得出斜桿直徑為2.0 mm時(shí)，內(nèi)插中心斜桿換熱管的綜合換熱性能較好。