帶有三角翼渦發(fā)生器的螺旋通道強化換熱研究

摘" " " 要： 采用數(shù)值模擬方法探究了帶有三角翼型渦發(fā)生器的半圓形截面螺旋通道內(nèi)流體流動與強化傳熱性能，考察了相鄰三角翼間隔角φ、長高比Г對強化傳熱性能的影響。結(jié)果表明：通道內(nèi)安裝三角翼的φ值越小，對流體的整體擾動作用越強，強化換熱效果越好，但同時流動阻力也越大；研究范圍內(nèi)，φ=180°布置的三角翼綜合強化傳熱效果最佳。同一雷諾數(shù)Re條件下，三角翼Г值越大綜合傳熱系數(shù)PEC值越高。φ=180°、Г=2.50時，綜合傳熱系數(shù)PEC在1.029～1.062之間。

關(guān)" 鍵" 詞：三角翼；螺旋通道；間隔角；強化換熱

中圖分類號：TQ021.3" " "文獻標(biāo)識碼： A" " "文章編號： 1004-0935（2023）05-0661-05

螺旋通道被廣泛應(yīng)用于石油化工、生物制藥、機械制造、制冷及能源動力等領(lǐng)域，具有制造簡單、結(jié)構(gòu)緊湊、高效換熱等優(yōu)點[1-2]。與其他通道相比，螺旋通道的主要關(guān)注點之一是二次流動，二次流是由離心力的作用引起的，有利于破壞或減薄邊界層，促進主流和壁面附近流體顆粒的混合，是強化傳熱傳質(zhì)的關(guān)鍵[3]。眾多科研學(xué)者探究了不同截面形狀螺旋通道內(nèi)流體流動和傳熱特性，其中半圓形截面螺旋通道可認(rèn)為是螺旋半管夾套的簡化模型[4-6]。螺旋半管夾套常常安裝在反應(yīng)釜外側(cè)用來傳遞反應(yīng)釜內(nèi)的反應(yīng)預(yù)熱以確保反應(yīng)釜內(nèi)溫度，從而大幅度降低化學(xué)反應(yīng)過程中的生產(chǎn)成本，產(chǎn)生巨大的經(jīng)濟效益。已有研究結(jié)果表明，高雷諾數(shù)下，螺旋通道中二次流的強化傳熱效果降低，導(dǎo)致螺旋通道的綜合強化傳熱性能降低。為了進一步提高螺旋通道的傳熱性能，各種強化傳熱方式已被學(xué)者應(yīng)用于螺旋通道以改善流場，實現(xiàn)強化換熱。

近些年來，采用固定安裝的渦發(fā)生器強化通道內(nèi)流體的換熱已被大量文獻報道[7-8]。其中，翼型渦發(fā)生器能夠產(chǎn)生作用范圍較大的縱向渦旋，有效減小通道內(nèi)的流動死區(qū)，在無需外加功率的條件下，能夠以增加少量壓力損失為代價來實現(xiàn)強化換" "熱[9-10]。翼型渦發(fā)生器產(chǎn)生的縱向渦旋會引起強烈的旋流運動，增加了流體混合，減薄或者破壞了熱邊界層的連續(xù)發(fā)展，使流體與換熱表面間的傳熱得以改善[11]。已有研究結(jié)果表明，渦流發(fā)生器的布置數(shù)量和結(jié)構(gòu)參數(shù)均會對換熱效果產(chǎn)生重要影響。SHAFIQUL[12]等在方形通道內(nèi)周期性布置渦流發(fā)生器，并對流體的換熱特性進行研究，結(jié)果表明傳熱性能明顯優(yōu)于光滑通道，但流動阻力同時增加。賈龍[13]等對比了渦發(fā)生器結(jié)構(gòu)參數(shù)對細(xì)微直通道換熱性能的影響，結(jié)果表明盡管長高比對流體換熱性能影響較小，但隨著長高比的增加流動阻力隨之" 增加。

本文通過三角單翼型渦流發(fā)生器復(fù)合半圓形截面螺旋通道達到強化換熱的目的，為深入分析安裝三角翼型渦發(fā)生器后流體的傳熱特性，主要考察了相鄰翼型渦發(fā)生器間隔角、結(jié)構(gòu)參數(shù)對螺旋通道內(nèi)傳熱效果的影響，為翼型渦發(fā)生器復(fù)合螺旋通道的強化傳熱方式在工程實際中的推廣和應(yīng)用提供理論參考。

1" 物理模型及幾何參數(shù)

物理模型如圖1（a）所示。半圓形截面螺旋通道螺旋半徑、螺距、半圓管半徑分別為：Rc=80 mm、P=40 mm、r0=10 mm，圈數(shù)為4圈。起始三角翼安裝在距螺旋通道入口θ=180°位置處。三角翼的結(jié)構(gòu)及布置方式如圖1（b）所示，翼的長和高取值分別為a=8 mm，b=4 mm，翼厚c=1 mm。三角翼斜邊迎流布置，定義長邊a與螺旋線夾角為攻角α，取值為α=30°。本研究中主要考察了間隔角φ與長高比Г對三角翼強化螺旋通道傳熱性能的影響，其中間隔角φ的取值分別為90°、180°和360°，長高比Г取值范圍為1.5～2.5。

2" 研究方法

2.1" 數(shù)值模擬

運用Fluent軟件對安裝多個三角翼的螺旋通道進行數(shù)值模擬計算。基于Realizable k-ε湍流模型，采用SIMPLEC算法，二階迎風(fēng)格式離散，壁面處理采用標(biāo)準(zhǔn)壁面函數(shù)。入口邊界條件為速度入口，入口溫度Tin=293 K，出口為壓力出口。在計算過程中半圓形直壁面為換熱壁面，定義壁面溫度 Tw=373 K，外壁面為絕熱。收斂條件為連續(xù)性方程和動量方程殘差小于1.0×10-5，能量方程殘差小于1.0×10-6。

對計算域網(wǎng)格采用分塊劃分，三角翼周圍局域為非結(jié)構(gòu)化網(wǎng)格，并針對此區(qū)域進行局部加密，無三角翼區(qū)域采用結(jié)構(gòu)化網(wǎng)格，并進行網(wǎng)格獨立性驗證以確保計算結(jié)果的準(zhǔn)確性。

2.2" 數(shù)值模擬結(jié)果的驗證

為驗證數(shù)值模擬的正確性，采用文獻[14]中的實驗裝置對光滑半圓形截面螺旋通道內(nèi)換熱特性進行實驗研究。實驗中螺旋通道曲率半徑、螺距、半圓管半徑分別為：Rc=80 mm、P=40 mm、r0=10 mm，圈數(shù)為4。圖2給出了對流換熱系數(shù)hm和壓差Δp的模擬值與實驗值的對比結(jié)果。

由圖2可以看出，模擬值與實驗值的趨勢吻合較好，hm相對偏差絕對值3.54%～6.79%，Δp相對偏差絕對值為3.35%～6.12%。由此可見，本研究采用的數(shù)值模擬方法是可靠的。

3" 結(jié)果與討論

3.1" 間隔角φ對換熱性能的影響

為了分析相鄰三角翼間隔角φ對換熱和流動阻力的影響，圖3和圖4分別給出了不同φ值下努塞爾數(shù)Nu和阻力系數(shù)f的變化曲線。由圖3可以看出，隨著φ值的減小，Nu值逐漸增加。這是由于隨著φ值的減小，螺旋通道內(nèi)三角翼的數(shù)量越多，產(chǎn)生縱向渦旋的數(shù)量增加，促進了流體混合程度，使換熱效果增強。

由圖4可以看出，與單一螺旋通道相比，在通道內(nèi)安裝三角翼導(dǎo)致阻力系數(shù)f顯著增加，隨著φ值的減小，f值逐漸增加，導(dǎo)致此現(xiàn)象的原因是，一方面流體繞過三角翼后產(chǎn)生縱向渦旋會造成一定的能量損失，另一方面渦旋隨著流動演變過程中，減薄了換熱邊界層的同時，導(dǎo)致了近壁處流體速度梯度的增加，進而使摩擦阻力變大。而小φ值下，三角翼數(shù)量的增多會造成更多的流動阻力損失。

從圖3和圖4中可以看出，φ值越小，Nu和f值同時增加，為了綜合評價三角翼間隔角φ對通道綜合強化換熱效果的影響，定義綜合強化傳熱因子PEC如下。

圖5給出了不同φ值下的PEC值，可以看出研究范圍內(nèi)PEC值均大于1，說明三角翼可以在半圓形螺旋通道內(nèi)獲得良好的強化傳熱效果。同時可以看出φ=180°的三角翼的綜合強化換熱效果最好。分析結(jié)果表明，研究范圍φ=180°時，PEC在" " " 1.050 9～1.018 4之間。

3.2" 長高比Γ對換熱性能的影響

為了考察三角翼長高比Γ對半圓形通道換熱效果的影響，圖6、圖7分別給出了努塞爾數(shù)Nu、流動阻力系數(shù)f隨Re的變化曲線。

從圖6、圖7中可以看出，隨著Γ值增加，Nu和f值同時增加，同單一螺旋通道相比，Nu值增加了3.79%～10.56%，f值增加了8.83%～13.34%。這是由于隨著Γ值越大，流體繞過三角翼時形成的二次渦旋強度越強，范圍越大，促進強化換熱效果。但同時通道內(nèi)由三角翼引起的形狀阻力也隨之增大，導(dǎo)致通道內(nèi)流動阻力的增加。圖8給出了三角翼不同Г下的PEC值，從圖中可以看出，所有Г值下的PEC值都高于1，說明在螺旋通道中安裝三角翼是一種有效的強化換熱措施。同一雷諾數(shù)下，三角翼Г越大PEC值越高，即綜合強化換熱效果最佳。

4" 結(jié) 論

通過數(shù)值模擬方法，研究了帶有三角翼型渦發(fā)生器的螺旋通道內(nèi)流體的強化傳熱特性，考察了三角翼間隔角以及長高比對強化傳熱效果的影響。得到如下結(jié)論：

1）相鄰三角翼間隔角φ越小，螺旋通道內(nèi)努塞爾數(shù)Nu和阻力系數(shù)f越大，φ=180°時具有最佳的換熱效果。

2）隨著三角翼長高比Γ的增加，更有利于產(chǎn)生強度更大且作用范圍更大的縱向渦旋，能夠有效促進通道內(nèi)流體的充分混合，達到高效的強化換熱效果，但同時導(dǎo)致螺旋通道內(nèi)流體流動阻力增大，研究范圍內(nèi)Г=2.5時三角翼具有最佳的綜合換熱" 性能。

參考文獻：

［1］ZHANG L Y， CUI X， LU Z， et al. A novel spiral channel with the growing waviness on the sidewalls for compact high-efficiency heat exchanger[J]. Applied Energy，2021，299： 117332.

［2］李雅俠，史曉航，張靜，等．對翼幾何參數(shù)對半圓螺旋通道內(nèi)流體流動特性的影響[J].過程工程學(xué)報，2021，21（6）：671-679．

［3］RENNIE T J， RAGHAVAN V G S. Numerical studies of a double-pipe helical heat exchanger[J]. Applied Thermal Engineering， 2006， 26（11-12）：1266-1273.

［4］WANG Y F， LIU P， SHAN F， et al. Effect of longitudinal vortex generator on the heat transfer enhancement of a circular tube[J]. Applied Thermal Engineering， 2019， 148：1018-1028.

［5］李雅俠，許鵬宇，韓澤民，等．矩形截面高寬比對射流強化螺旋通道傳熱性能的影響[J].過程工程學(xué)報，2022，22（5）：612-621．

［6］李凡，陸高鋒，馬光柏，等．縱向渦強化圓管內(nèi)換熱的數(shù)值模擬及性能分析[J].化工學(xué)報，2021，72（S1）：120-126．

［7］吳佩霖，林清宇，艾鑫，等．渦發(fā)生器對矩形截面螺旋細(xì)通道內(nèi)液體傳熱與熵產(chǎn)的影響[J].輕工機械，2020，38（1）：28-34．

［8］AWAIS M， BHUIYAN A A. Heat transfer enhancement using different types of vortex generators （VGs）： A review on experimental and numerical activities[J]. Thermal Science and Engineering Progress， 2018， 5：524-545.

［9］李雅俠，吳凱，史曉航，等．矩形翼強化半圓螺旋通道內(nèi)流體流動與換熱[J].沈陽化工大學(xué)學(xué)報，2021，35（1）：54-60．

［10］LI Y X， WANG X， ZHANG J， et al. Comparison and analysis of the arrangement of delta winglet pair vortex generators in a half coiled jacket for heat transfer enhancement[J]. International Journal of Heat and Mass Transfer， 2019， 129：287-298.

［11］倪亞琴．渦流發(fā)生器研制及其對邊界層的影響研究[J].空氣動力學(xué)學(xué)報，1995（1）：110-116．

［12］ISLAM M S， HINO R， HAGA K， et al. Experimental study on heat transfer augmentation for high heat flux removal in rib-roughened narrow channels[J]. Journal of nuclear science and technology， 1998， 35（9）： 671-678.

［13］賈龍，張劍飛．縱向渦發(fā)生器幾何參數(shù)對微通道流動傳熱特性的影響研究[J].熱能動力工程，2019，34（10）：97-105．

［14］李雅俠，張騰，張春梅，等．翼型渦發(fā)生器對半圓形螺旋通道的換熱強化機理[J].化工學(xué)報，2016，67（5）：1814-1821．

Abstract：" The fluid flow and enhanced heat transfer performance in a helical channel with a semicircular section with a triangular airfoil vortex generator were investigated by numerical simulation method. The effects of the interval angle φ and aspect ratio Г of adjacent triangular winglet on the enhanced heat transfer performance were investigated. The results showed that the smaller the φ value of the triangular winglet installed in the channel， the stronger the overall disturbance effect on the fluid， the better the heat transfer enhancement effect， but the greater the flow resistance; Within the research range， the triangular winglet arranged when φ=180° had the best comprehensive heat transfer enhancement effect. Under the same Reynolds number Re， the larger the value of Г of the triangular winglet， the higher the PEC value of the comprehensive heat transfer coefficient. When φ=180° and Г=2.50， the comprehensive heat transfer coefficient PEC was between 1.029 and 1.062.

Key words： Triangular winglet; Helical channel; Interval angle; Enhanced heat exchange