閆順林,于興寶,張 莎
?
麻面式過熱器煙氣側(cè)傳熱與阻力特性數(shù)值研究
閆順林,于興寶,張 莎
(華北電力大學(xué)能源動力與機(jī)械工程學(xué)院,河北 保定 071003)
以麻面式過熱器交叉流動通道為研究對象,采用三維數(shù)值模擬的方法,對麻面管煙氣側(cè)換熱與阻力特性進(jìn)行研究。結(jié)果表明:麻面管外壁凹坑與管束布置形式改善了煙氣對流沖刷效果,破壞了氣體流動邊界層,縮小了傳熱惡化的尾流區(qū)面積;與光滑管束相比,順列與錯列麻面管束煙氣側(cè)的努塞爾數(shù)(Nu)分別提升了8%~63%、12%~48%,綜合換熱性能指標(biāo)CPEC分別為1.05~1.50、1.07~1.28;以煙氣側(cè)Nu、阻力系數(shù)f、CPEC為目標(biāo)函數(shù),優(yōu)化得到了最佳參數(shù)組合:順列管束橫向節(jié)距s1/Do=2.10,縱向節(jié)距s2/Do=1.75,對應(yīng)的Nu=42.63,CPEC=1.24,錯列管束s1/Do=2.55,s2/Do=2.50,對應(yīng)的Nu=47.02,CPEC=1.31;運(yùn)用多元線性回歸的方法,分別擬合出順列與錯列管束煙氣側(cè)的換熱準(zhǔn)則關(guān)聯(lián)式。本文研究結(jié)果可為麻面管的工程應(yīng)用提供參考。
麻面管;過熱器;強(qiáng)化換熱;換熱準(zhǔn)則關(guān)聯(lián)式;數(shù)值模擬
換熱管作為電廠中常用的換熱元件,廣泛應(yīng)用于過熱器、省煤器、凝汽器等換熱設(shè)備中[1]。電廠中換熱元件多為普通光滑管,熱阻較大。研究高效率、經(jīng)濟(jì)實(shí)用的新型強(qiáng)化換熱設(shè)備是減少電廠能耗、提高電廠經(jīng)濟(jì)效益的有效途徑之一。
強(qiáng)化換熱技術(shù)被廣泛應(yīng)用到實(shí)際工程中,表現(xiàn)出了良好的強(qiáng)化換熱能力[2-3]。無源強(qiáng)化換熱技術(shù)無需借助任何外力即可實(shí)現(xiàn)換熱強(qiáng)化[4],這項技術(shù)的研究與開發(fā)受到了國內(nèi)外學(xué)者的廣泛關(guān)注,Mahmood等人[5]利用大渦模擬技術(shù)對布置凹坑凸胞結(jié)構(gòu)的板片通道進(jìn)行了換熱研究,結(jié)果發(fā)現(xiàn)當(dāng)凹坑凸胞的深度增大時換熱性能提升,同時壓降損失隨之增加。Thianpong等人[6]通過實(shí)驗(yàn)研究了裝有扭曲帶渦旋發(fā)生器的凹坑管內(nèi)的流動與換熱特性,發(fā)現(xiàn)其換熱性能明顯增強(qiáng),且裝有扭曲帶的凹窩管的傳熱系數(shù)和阻力系數(shù)均隨著節(jié)距比和扭曲比的減小而增加。石磊等[7]對空氣冷凝器橢圓管外繞橢圓翅片的三排管外空氣側(cè)進(jìn)行數(shù)值模擬,結(jié)果表明隨著迎風(fēng)風(fēng)速的增大,三排管的流動阻力和傳熱系數(shù)隨之增加,并得出三排管外空氣流動阻力和平均傳熱系數(shù)的擬合計算式。Chen等人[8]以水為工質(zhì),利用實(shí)驗(yàn)研究了凹坑管換熱與流動特性,結(jié)果發(fā)現(xiàn)與光滑管相比,在恒定雷諾數(shù)和恒定泵功率下凹坑管傳熱分別增強(qiáng)了25%~137%、15%~84%。
電廠過熱器工作在高溫條件下,環(huán)境惡劣,而目前過熱器光滑管換熱性能較差、能源利用率低。對此,本文提出一種新式換熱表面的過熱器管型—麻面管,應(yīng)用流體力學(xué)分析軟件模擬麻面管煙氣側(cè)傳熱與阻力特性,并進(jìn)行煙氣側(cè)相關(guān)參數(shù)的優(yōu)化研究,旨在為這種新式換熱表面管型在強(qiáng)化換熱方面的應(yīng)用與推廣提供參考。
麻面式過熱器基管的內(nèi)徑i=43 mm,外 徑o=51 mm,管內(nèi)外壁面凸胞/凹坑投影直徑=6.0 mm,凸胞/凹坑深度=1.0 mm,同一截面圓周上凸胞/凹坑間呈40°等間距均勻分布,并且軸向呈交錯式布置,前后凸胞/凹坑間距離=16 mm,凹坑采用倒角分段平滑過渡。
將過熱器幾何模型進(jìn)行簡化,以多根麻面式過熱器管交叉流動通道為研究對象,管長,管束分別以順列與錯列方式排布,模擬研究煙氣側(cè)工質(zhì)流速、橫向節(jié)距1/o、縱向節(jié)距2/o對麻面管束交叉流動和換熱的影響。麻面管三維結(jié)構(gòu)、凸胞/凹坑的結(jié)構(gòu)示意、管束布置形式如圖1—圖3所示。
圖1 麻面管結(jié)構(gòu)示意
圖2 凸胞/凹坑的結(jié)構(gòu)示意
圖3 過熱器管錯列與順列布置形式
假設(shè)麻面管束兩側(cè)流體不可壓縮且穩(wěn)態(tài)流動,忽略重力、浮升力的影響,且滿足質(zhì)量、動量和能量守恒定律,則有[9]:
質(zhì)量守恒方程
動量守恒方程
流體域能量守恒方程
固體域能量守恒方程
式中:,=,,;為流體速度,m/s;為壓力,Pa;為密度,kg/m3;eff為有效熱導(dǎo)率,W/(m·K);為內(nèi)能,J/kg;為溫度,K;h為熱源項,J;J,為組分?jǐn)U散通量;τ為應(yīng)力張量;為固體導(dǎo)熱系數(shù),W/(m·K)。
在三維軟件SolidWorks中創(chuàng)建幾何模型,模擬過程在Fluent 17.0中完成。采用ANSYS中的Mesh功能進(jìn)行網(wǎng)格劃分,在管內(nèi)外壁面附近區(qū)域進(jìn)行網(wǎng)格縮小與加密處理,更好地展現(xiàn)麻面管內(nèi)外表面凸胞/凹坑結(jié)構(gòu),以提高近壁處邊界層內(nèi)數(shù)值求解的準(zhǔn)確度。網(wǎng)格劃分結(jié)果如圖4所示。
圖4 網(wǎng)格劃分示意
數(shù)值模擬采用比標(biāo)準(zhǔn)-方程更加適應(yīng)于低雷諾數(shù)流動、旋流、二次流、分離流等復(fù)雜流動的RNG-方程[10],壁面采用增強(qiáng)壁面函數(shù)法處理[11],壓力-速度耦合采用Simplec算法,壓力采用Standard格式,控制方程中其他各項采用二階迎風(fēng)格式,亞松弛因子采用默認(rèn)設(shè)置,模擬結(jié)果的能量殘差收斂精度達(dá)到10-6,其他殘差收斂值達(dá)到10-4。
邊界條件設(shè)置:管外工質(zhì)為煙氣,入口為速度入口,溫度為1 206 K,出口為壓力出口;管內(nèi)工質(zhì)為過熱蒸汽,入口為質(zhì)量流量入口,數(shù)值為0.2 kg/m3,溫度為752 K,出口為壓力出口;管內(nèi)外壁為無滑移邊界條件,其他壁面邊界條件如圖5所示;兩側(cè)流體與管壁交界面采用耦合壁面邊界條件,耦合交界面處滿足溫度及熱流密度連續(xù)(式(5))。
式中:T為溫度,K;q為熱流密度,W/m2;下標(biāo)s、f表示固體域、流體域;下標(biāo)w表示交界面。
采用等泵功率下的評價準(zhǔn)則來表征強(qiáng)化換熱管蒸汽側(cè)的綜合換熱性能。綜合換熱性能指標(biāo)PEC的定義式為[12]
式中,g分別為麻面管束與光滑管束煙氣側(cè)的平均努塞爾數(shù),、g分別為麻面管束與光滑管束煙氣側(cè)的平均阻力系數(shù)。
利用Design-Expert 8.0的設(shè)計、分析及優(yōu)化功能[13],以煙氣側(cè)、、PEC為多目標(biāo)函數(shù),麻面管束橫向節(jié)距1/0與縱向節(jié)距2/0為設(shè)計變量進(jìn)行參數(shù)優(yōu)化模擬研究。設(shè)計參數(shù)見表1。
計算中,換熱管煙氣側(cè)努賽爾數(shù)和阻力系數(shù)的數(shù)學(xué)表達(dá)式為:
式中:為換熱管平均對流換熱系數(shù),W/(m2·K);e為換熱管當(dāng)量直徑,m;為流體導(dǎo)熱系數(shù),W/(m·K);Δ為工質(zhì)進(jìn)出口壓力差,Pa;為過熱器管根數(shù);為煙氣密度,kg/m3;為煙氣平均流速,m/s。
表1 設(shè)計參數(shù)
Tab.1 The design parameters
以順列式麻面管束的結(jié)構(gòu)參數(shù)1/o=2.75、2/o=2.05為例,蒸汽側(cè)入口流量為0.2 kg/s,煙氣側(cè)入口速度為6.5 m/s,對模型進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證。依次劃分500萬、610萬、700萬、790萬網(wǎng)格,如圖6所示。隨著網(wǎng)格數(shù)增加,煙氣側(cè)和變化速率減緩,當(dāng)數(shù)量高于700萬時,其數(shù)值基本保持不變,網(wǎng)格數(shù)量及密度符合計算要求。
圖6 網(wǎng)格無關(guān)性驗(yàn)證
以光滑管作為基準(zhǔn)管進(jìn)行模型準(zhǔn)確性驗(yàn)證,模擬的參數(shù)設(shè)置按前所述。模擬結(jié)果與Zhukauskas公式[14]的計算數(shù)值進(jìn)行對比,煙氣側(cè)數(shù)和阻力系數(shù)的誤差在±10%以內(nèi),驗(yàn)證了模擬的準(zhǔn)確性。
以蒸汽側(cè)質(zhì)量流量0.2 kg/s煙氣側(cè)流速11.0 m/s,1/o=2.50,2/o=2.05為例進(jìn)行場態(tài)分析。圖7—圖9為煙氣側(cè)流場、速度場等值線圖及溫度場云圖。
從圖7、圖8可見:外壁凹坑使得煙氣流動狀態(tài)發(fā)生了改變,且速度分布較均勻;麻面管背風(fēng)側(cè)對換熱不利的尾流區(qū)[15]面積較光滑管明顯減小,最小速度出現(xiàn)在尾流渦區(qū);麻面管煙氣側(cè)平均流速高于光滑管,說明麻面管煙氣區(qū)湍流擾動加劇。
從圖9可見:凹坑的存在使得煙氣側(cè)高溫尾流區(qū)域面積縮小,越靠近凹坑,其附近流體域溫度越低,近壁處麻面管高溫域面積明顯小于光滑管,壁面的降溫速率增加;麻面管煙氣區(qū)平均溫度低于光滑管區(qū)域,說明凹坑的存在與管束的排布形式改善了煙氣對流沖刷效果,使得煙氣擾動加強(qiáng),破壞了煙氣流動邊界層,同時縮小了惡化換熱的尾流區(qū)域面積,實(shí)現(xiàn)了煙氣側(cè)的強(qiáng)化換熱。
圖10為順列管束時、1/o、2/o對煙氣側(cè)、的影響。由圖10可見:在不同、順列管束1/o、2/o下,與光滑管束相比,麻面式過熱器管束煙氣側(cè)明顯提升了傳熱效果,經(jīng)計算煙氣側(cè)提高了8%~63%;同時,隨著1/o減小、2/D增大,麻面管煙氣側(cè)增加明顯,同時阻力隨之增大,1/o<2.75、2/o>2.05時阻力增加明顯;隨著增加,隨之增大,阻力卻隨之減小且逐漸趨于平緩??梢?,適當(dāng)減小橫向節(jié)距、增大縱向節(jié)距均有助于增強(qiáng)順列麻面管束的傳熱效果。
圖11為錯列管束時、1/o、2/o對煙氣側(cè)、的影響。
由圖11可見:隨著錯列管束橫向節(jié)距1/o與縱向節(jié)距2/o的減小,麻面管的換熱性能得到增強(qiáng),相對而言管束橫向節(jié)距對煙氣側(cè)換熱影響較大,且阻力系數(shù)也隨之增加;隨著的增加,煙氣側(cè)隨之增大,阻力降低速率減慢并逐漸趨于平緩,適當(dāng)?shù)販p小錯列麻面管束橫向與縱向節(jié)距有助于增強(qiáng)傳熱效果。經(jīng)計算,在不同,錯列管束1/o、2/o下,與光滑管束相比,麻面管煙氣側(cè)提高了12%~48%,換熱性能顯著提升。
圖12為、1/o、2/o對順列與錯列管束煙氣側(cè)PEC的影響情況。由圖12可見:順列與錯列管束煙氣側(cè)的PEC分別為1.05~1.50、1.07~1.28;煙氣側(cè)PEC隨的增加而逐漸遞減并趨于平緩;相比于縱向節(jié)距,管束的橫向節(jié)距對煙氣側(cè)的綜合換熱性能影響較大。因此,可以通過改變管束橫向節(jié)距的大小來強(qiáng)化煙氣側(cè)的換熱性能。
圖12 Re、s1/Do、s2/Do對煙氣側(cè)CPEC的影響
水平煙道內(nèi)煙溫高,灰粒較軟,煙氣流速控制在10~15 m/s[16],管外煙速以11 m/s作為設(shè)計流速。以麻面管煙氣側(cè)、、PEC為多目標(biāo)函數(shù),利用軟件Design-Expert的優(yōu)化功能,取整得到了最佳參數(shù)組合為:順列管束1/o=2.10,2o=1.75,對應(yīng)的=42.63,PEC=1.24;錯列管束1/o=2.55,2o=2.50,對應(yīng)的=47.02,PEC=1.31。
通過對麻面式過熱器管煙氣側(cè)傳熱與阻力特性的數(shù)值分析,可知換熱管外徑o、煙氣側(cè)、橫向節(jié)距1/o縱向節(jié)距2/o對煙氣側(cè)影響重大。因此,以、1/o2/o作為煙氣側(cè)平均準(zhǔn)則關(guān)聯(lián)式的準(zhǔn)則變量,運(yùn)用多元線性回歸的方法分別擬合出順列與錯列管束煙氣側(cè)換熱的準(zhǔn)則關(guān)聯(lián)式,其數(shù)學(xué)表達(dá)見式(9)、式(10)。
順列管束:
式(9)適用范圍為=500~3 000,2.0≤1/o≤3.5,1.6≤2/o≤2.5。由式(9)計算所得的煙氣側(cè)平均數(shù)與模擬結(jié)果的標(biāo)準(zhǔn)誤差約為2.87%。
錯列管束:
式(10)適用范圍為=500~3 000,2.5≤1/o≤3.5,1.6≤2/o≤2.5。由式(10)計算所得的煙氣側(cè)平均數(shù)與模擬結(jié)果的標(biāo)準(zhǔn)誤差約為0.62%。
根據(jù)傳熱學(xué)理論[14],忽略換熱管內(nèi)外污垢的影響,以管外表面為基準(zhǔn)的總換熱系數(shù)數(shù)學(xué)表達(dá)式為
式中:i、o分別為管內(nèi)外表面換熱系數(shù),W/(m2·K);b為管壁導(dǎo)熱系數(shù),W/(m·K)。
以蒸汽側(cè)質(zhì)量流量0.2 kg/m3,順列與錯列管束煙氣側(cè)最佳參數(shù)組合為例,通過式(11)計算得出:順列麻面管束總換熱系數(shù)=73.58 W/(m2·K),是相同條件下光滑管束的1.51倍;錯列麻面管束總換熱系數(shù)=82.11 W/(m2·K),是相同條件下光滑管的1.38倍??梢姡槊媸綋Q熱管具有一定的強(qiáng)化換熱性能。
1)麻面管外壁的凹坑與管束的排布形式改善了煙氣對流沖刷效果,使得煙氣擾動加劇,破壞了氣體流動邊界層,同時縮小了惡化換熱的尾流區(qū)面積,實(shí)現(xiàn)了煙氣側(cè)的強(qiáng)化換熱。
2)與光滑管束相比,順列麻面管束煙氣側(cè)提高了8%~63%,PEC為1.05~1.50;錯列麻面管束煙氣側(cè)增加了12%~48%,PEC為1.07~1.28,適當(dāng)?shù)馗淖児苁臋M向與縱向節(jié)距有助于增強(qiáng)換熱。
3)順列麻面管束最佳參數(shù)組合為1/o=2.1,s/Do=1.75,對應(yīng)的=42.63,PEC=1.24;錯列麻面管束最佳參數(shù)組合為1/o=2.55,2o=2.50,對應(yīng)的=47.02,PEC=1.31。
4)運(yùn)用多元線性回歸的方法分別擬合出順列與錯列麻面管束煙氣側(cè)的換熱準(zhǔn)則關(guān)聯(lián)式,并給出其適用范圍,為麻面管的工程應(yīng)用提供參考。
[1] 閆順林, 王皓軒. 煙氣橫掠麻面管束顆粒沉積特性的數(shù)值分析[J]. 電力科學(xué)與工程, 2017, 33(11): 66-72.
YAN Shunlin, WANG Haoxuan. Numerical analysis on particle deposition characteristics of flue gas flowing across pitted tubes in aligned arrangement[J]. Electric Power Science and Engineering, 2017, 33(11): 66-72.
[2] CHENG Z D, HE Y L, CUI F Q. Numerical study of heat transfer enhancement by unilateral longitudinal vortex generators inside parabolic trough solar receivers[J]. International Journal of Heat and Mass Transfer, 2012, 55(21): 5631-5641.
[3] KUKULKA D J, SMITH R. Thermal-hydraulic performance of Vipertex 1EHT enhanced heat transfer tubes[J]. Applied Thermal Engineering, 2013, 61(1): 60-66.
[4] BERGLES A E. Some perspectives on enhanced heat transfer second generation heat transfer technology[J]. ASME Journal of Heat Transfer, 1988, 110(4b): 1082-1096.
[5] MAHMOOD G I, LIGRANI P M. Heat transfer in a dimpled channel: combined influences of aspect ratio, temperature ratio, Reynolds number and flow structure[J]. International Journal of Heat and Mass Transfer, 2002, 45(10): 2011-2020.
[6] THIANPONG C , EIAMSA-ARD P , WONGCHAREE K, et al. Compound heat transfer enhancement of a dimpled tube with a twisted tape swirl generator[J]. International Communications in Heat & Mass Transfer, 2009, 36(7): 698-704.
[7] 石磊, 石誠, 劉海峰, 等. 直接空冷凝汽器三排管外空氣流動和傳熱性能的數(shù)值研究[J]. 熱力發(fā)電, 2009, 38(5): 14-17. SHI Lei, SHI Cheng, LIU Haifeng, et al. Numerical investigation on air flow and heat transfer performance outside triple row tube used in direct air cooling conden- ser[J]. Thermal Power Generation, 2009, 38(5): 14-17.
[8] CHEN J, MüLLER-STEINHAGEN H, DUFFY G G. Heat transfer enhancement in dimpled tubes[J]. Applied Thermal Engineering, 2001, 21(5): 535-547.
[9] 王尊策, 何寶林, 韓建荒, 等. 高溫?fù)Q熱器多場耦合數(shù)值模擬研究[J]. 化工機(jī)械, 2013, 40(2): 188-192. WANG Zunce, HE Baolin, HAN Jianhuang, et al. Research on multifield coupling numerical simulation of high temperature heat exchangers[J]. Chemical Engi- neering & Machinery, 2013, 40(2): 188-192.
[10] 王福軍. 計算流體動力學(xué)分析—CFD軟件原理與應(yīng)用[M]. 北京: 清華大學(xué)出版社, 2004: 28. WANG Fujun. Computational fluid dynamics analysis: principles and applications of CFD software[M]. Beijing: Tsinghua University Press, 2004: 28.
[11] 李鵬飛, 徐敏義, 王飛飛. 精通CFD—工程仿真與案例實(shí)戰(zhàn)[M]. 北京: 人民郵電出版社, 2017: 6. LI Pengfei, XU Minyi, WANG Feifei. Proficient in CFD: engineering simulation and case study[M]. Beijing: Post and Telecom Press, 2017: 6.
[12] 張定才, 何雅玲, 劉啟斌, 等. 內(nèi)螺旋肋管流動與傳熱特性的實(shí)驗(yàn)研究[J]. 工程熱物理學(xué)報, 2006, 27(6): 1029-1031. ZHANG Dingcai, HE Yaling, LIU Qibin, et al. Experimental study of fluid flow and heat transfer in internal helical-rib tubes[J]. Journal of Engineering Thermophysics, 2006, 27(6): 1029-1031.
[13] 葛宜元. 試驗(yàn)設(shè)計方法與Design-Expert軟件應(yīng)用[M]. 哈爾濱: 哈爾濱工業(yè)大學(xué)出版社, 2014: 37. GE Yiyuan. Experimental design method and application of Design-Expert software[M]. Harbin: Harbin Institute of Technology Press, 2014: 37.
[14] 楊世銘, 陶文銓. 傳熱學(xué)[M]. 4版. 北京: 高等教育出版社, 2014: 66. YANG Shiming, TAO Wenquan. Heat transfer[M]. 4th ed. Beijing: Higher Education Press, 2014: 66.
[15] HEGGS P J, STONES P R. The effects of non-uniform heat transfer coefficients in the design of finned tube air-cooled heat exchangers[C]//Proceeding of International Heat Transfer Conference. Munich: International Institute of Refrigeration, 1982: 209-214.
[16] 周強(qiáng)泰. 鍋爐原理[M]. 3版. 北京: 中國電力出版社, 2013: 7. ZHOU Qiangtai. Boiler principle[M]. 3rd ed. Beijing: China Electric Power Press, 2013: 7.
Numerical study on heat transfer and resistance characteristics of gas side of superheater with new pitted tubes
YAN Shunlin, YU Xingbao, ZHANG Sha
(School of Energy, Power and Machinery Engineering, North China University of Electric Power, Baoding 071003,China)
Taking the cross flow path of superheater with pitted tubes as the research object, three-dimensional numerical simulation was performed to investigate the heat transfer and resistance characteristics of the flue gas side of the pitted tubes. The results show that, the added pits on external wall of the pitted tubes and the proper arrangement of tubes improve the flue gas convection and scour effect, damage the gas flow boundary layer, and reduce the area of heat flow deterioration region. Compared with the smooth tubes, at flue gas side of the aligned pitted tubes and staggered tubes, the Nusselt number (Nu) increases by 8%~63% and 12%~48%, and the overall heat transfer performance evaluation criterion (CPEC) is 1.05~1.5 and 1.07~1.28. Taking Nu, f and CPEC at the flue gas side as the objective functions, a group of optimal parameters are obtained. The optimal parameters of in-line bank are 2.1 for transversal pitch and 1.75 for longitudinal pitch under the condition of Nu=42.63, CPEC=1.24 while the optimal parameters of staggered bank are 2.55 for transversal pitch and 2.5 for longitudinal pitch under the condition of Nu=47.02, CPEC=1.31. The criterion formula of heat transfer at flue gas side of the aligned and staggered pitted tubes are obtained respectively by the method of multiple linear regression. The research results of this paper can provide references for the engineering application of pitted tubes.
pitted tube, superheater, enhanced heat transfer, heat transfer criterion formula, numerical simulation
TK124
A
10.19666/j.rlfd.201810177
閆順林, 于興寶, 張莎. 麻面式過熱器煙氣側(cè)傳熱與阻力特性數(shù)值研究[J]. 熱力發(fā)電, 2019, 48(5): 18-24. YAN Shunlin, YU Xingbao, ZHANG Sha. Numerical study on heat transfer and resistance characteristics of gas side of superheater with new pitted tubes[J]. Thermal Power Generation, 2019, 48(5): 18-24.
2018-10-04
閆順林(1959—),男,博士,教授,主要研究方向?yàn)闊崃ο到y(tǒng)節(jié)能理論及應(yīng)用。
于興寶(1993—),男,碩士研究生,主要研究方向?yàn)闊崃ο到y(tǒng)節(jié)能理論及應(yīng)用,hd_yxb@126.com。
(責(zé)任編輯 劉永強(qiáng))