復(fù)合式波紋板傳熱與阻力特性數(shù)值研究

張 磊，杜 娟，白文剛，高 煒，吳帥帥

?

復(fù)合式波紋板傳熱與阻力特性數(shù)值研究

張 磊1，杜 娟2，白文剛1，高 煒1，吳帥帥1

（1.西安熱工研究院有限公司，陜西 西安 710054；2.西北農(nóng)林科技大學(xué)資源環(huán)境學(xué)院，陜西 楊凌 712100）

DN（double notched）型波紋板傳熱元件在回轉(zhuǎn)式空氣預(yù)熱器中應(yīng)用廣泛，為了揭示DN板波紋通道內(nèi)的微觀流動形態(tài)和強化傳熱機理，本文采用LBKE（Lam and Bremhorst k-ε model）湍流模型和計算流體動力學(xué)（CFD）數(shù)值模擬技術(shù)對DN板的傳熱、阻力特性展開研究，預(yù)測DN板的傳熱和阻力系數(shù)，分析DN板通道內(nèi)的流場和溫度場。模擬計算發(fā)現(xiàn)：DN板通道內(nèi)的流動可分為波紋通道流動和平直通道流動兩大類；DN板斜波紋段、凹槽段迎風(fēng)面努塞爾數(shù)（Nu）較大，背風(fēng)面Nu較小；Nu隨著雷諾數(shù)（Re）的增大而增大，DN板上的Nu分布隨著Re的增大而更趨均勻；增加DN板上斜波紋結(jié)構(gòu)比例可提升其傳熱能力，但阻力損失也會相應(yīng)增加；DN板在中等Re水平（Re≈5 000）下運行時效果最佳，此時既能保證DN板上Nu均勻分布又控制了阻力損失大小。

回轉(zhuǎn)式空氣預(yù)熱器；DN板；波紋板；數(shù)值模擬；傳熱；阻力系數(shù)

回轉(zhuǎn)式空氣預(yù)熱器（空預(yù)器）屬于再生式空預(yù)器的一種，具有結(jié)構(gòu)緊湊、安裝檢修方便、運行費用低、占地面積小等諸多優(yōu)點，是電站鍋爐不可缺少的組成部分[1-2]?；剞D(zhuǎn)式空預(yù)器利用鍋爐排煙的低品位余熱加熱入爐空氣，能夠改善爐膛內(nèi)著火條件和燃燒過程，并提高鍋爐熱效率。

傳熱芯體是回轉(zhuǎn)式空預(yù)器的核心部件，承擔(dān)著存儲、傳遞熱量的作用?；剞D(zhuǎn)式空預(yù)器的傳熱芯體一般由波紋板傳熱元件堆疊裝配而成。當(dāng)回轉(zhuǎn)式空預(yù)器運行時，圓盤狀的傳熱芯體連續(xù)旋轉(zhuǎn)，交替地與熱煙氣和冷空氣進行接觸換熱。當(dāng)煙氣與傳熱芯體接觸時，煙氣攜帶的熱量傳遞給芯體中的波紋板傳熱元件并蓄積起來；當(dāng)空氣與傳熱芯體接觸時，芯體中波紋板傳熱元件所蓄積的熱量又被傳遞給空氣。傳熱芯體旋轉(zhuǎn)時被周期性地加熱和冷卻，熱量就周期性地由煙氣傳遞給空氣[3]。

品質(zhì)優(yōu)良的回轉(zhuǎn)式空預(yù)器應(yīng)該具備以下特點：高效傳熱、低壓損、易吹灰和便于進行故障探測。為了滿足這些要求，需要對波紋板傳熱元件結(jié)構(gòu)進行細(xì)致的研究和設(shè)計[4]?；剞D(zhuǎn)式空預(yù)器波紋板傳熱元件一般由金屬平薄板壓制成波紋形狀，大量波紋板緊密堆疊，板片間隙構(gòu)成了工質(zhì)流動的波紋通道。波紋板通過特殊設(shè)計的壁面波紋形狀可誘發(fā)二次流動、增強流體混合，進而強化傳熱?；剞D(zhuǎn)式空預(yù)器中常用的典型波紋板型包括DN（double notched）板、CC（cross corrugated）板、CU（corrugated undulated）板和DU（double undulated）板等[5-9]。

國外DN板的研發(fā)、制造技術(shù)經(jīng)過長期積累已趨成熟，但由于技術(shù)保密的原因，相關(guān)研究的公開報道十分稀少。由于無法借鑒國外的先進技術(shù)經(jīng)驗，加之DN板波紋通道具有結(jié)構(gòu)復(fù)雜、小尺度的特點，使得國內(nèi)DN板型相關(guān)研究工作進展緩慢。對此，本文采用有限容積法對DN板通道內(nèi)的流動、傳熱過程進行三維數(shù)值模擬研究，揭示DN板通道內(nèi)的微觀流動形態(tài)和強化傳熱機理。本文的研究有望為DN型復(fù)合式波紋板傳熱元件的自主研發(fā)提供技術(shù)支持。

1 DN型波紋板

DN板由多種波紋結(jié)構(gòu)復(fù)合而成，包括凹槽段、斜波紋段以及平直段3種結(jié)構(gòu)，其結(jié)構(gòu)示意如圖1所示。圖1中凹槽段波紋形線為近似正弦波，凹槽方向與主流方向平行。斜波紋段波紋形線為近似三角波，斜波紋與主流方向夾角為，斜波紋段位于相鄰凹槽段之間。平直段為平直板片，位于相鄰斜波紋段之間。復(fù)合式波紋板的特殊設(shè)計結(jié)構(gòu)使得DN板具有傳熱高效、故障探測靈敏、成本低等優(yōu)點，這些優(yōu)良特性使DN板成為回轉(zhuǎn)式空預(yù)器中應(yīng)用最廣的一種傳熱元件。

2 計算區(qū)域及網(wǎng)格劃分

2.1 計算區(qū)域

考慮到回轉(zhuǎn)式空預(yù)器DN板傳熱芯體結(jié)構(gòu)在、、方向上周期性重復(fù)變化的特點，為了節(jié)省計算資源和計算時間，確定了具有代表性的最小單元作為計算區(qū)域，該區(qū)域為方向上3塊DN板間的2個波紋通道，計算區(qū)域?qū)挾龋ǚ较颍┑扔谝粋€波長（140 mm），計算區(qū)域長度（方向）為145 mm，計算區(qū)域如圖1所示。

2.2 網(wǎng)格劃分

針對不規(guī)則形狀的DN板通道，本文采用GAMBIT軟件中的非結(jié)構(gòu)化Tet/Hybrid網(wǎng)格對計算區(qū)域進行離散化處理，并在壁面附近實施了局部加密。經(jīng)過網(wǎng)格無關(guān)性驗證，本文最終的計算網(wǎng)格數(shù)為1 396 578個。DN板網(wǎng)格劃分如圖2所示。

圖2 DN板網(wǎng)格劃分

3 數(shù)學(xué)模型及求解策略

3.1 數(shù)學(xué)模型

描述回轉(zhuǎn)式空預(yù)器DN板通道內(nèi)流動傳熱過程的雷諾時均控制方程包括質(zhì)量守恒、動量守恒以及能量守恒方程：

需要引入湍流模型才能實現(xiàn)上述控制方程組的封閉求解。Patel等人[10]采用不同的湍流模型預(yù)測通道內(nèi)受限湍流流動過程，發(fā)現(xiàn)Lam和Bremhorst[11]提出的低雷諾數(shù)-（LBKE）模型計算精度最高。Ciofalo等人[12]在數(shù)值研究CC通道內(nèi)湍流流動傳熱過程時發(fā)現(xiàn)，LBKE模型在所選用的5種湍流模型中計算偏差最小。考慮到DN與CC板結(jié)構(gòu)的相似性，本文借鑒上述研究成果采用LBKE湍流模型計算DN板通道內(nèi)湍流流動過程。LBKE湍流模型的表達式為[11]：

式中：t為湍流黏性系數(shù)，kg/(m·s)；為湍流動能，J/kg；為湍動能耗散率；為密度，kg/m3；為時間，s；u、u為速度分量，m/s；x、x為直角坐標(biāo)分量，m；為流體黏性系數(shù)，kg/(m·s)；m、m、k、s、1、1、2、2為模型系數(shù)。

3.2 求解策略

本文利用Fluent軟件求解控制方程組并完成計算數(shù)據(jù)的后處理工作。求解策略為：

1）計算中忽略DN板片導(dǎo)熱熱阻、不考慮沾污的影響；

2）工質(zhì)為空氣，其物性參數(shù)均擬合為關(guān)于溫度的二次多項式；

3）邊界條件為速度入口、壓力出口、左右側(cè)面周期性邊界、壁面無滑移、恒壁溫；

4）壓力-速度耦合解法為SIMPLEC，對流項離散格式為QUICK，擴散項采用中心差分格式。

4 計算結(jié)果及討論

4.1 模型驗證

為了說明本文數(shù)值模擬方法的準(zhǔn)確性，首先將模擬結(jié)果與文獻[13]中的DN板通道實測數(shù)據(jù)進行了比較。圖3匯總比較了數(shù)值計算和瞬態(tài)實驗得到的DN板通道傳熱系數(shù)因子和阻力系數(shù)。由圖3可見，本文數(shù)值模擬得到的和值計算偏差分別為 –9.25%~11.38%、–17.46%~13.99%，計算結(jié)果和實驗數(shù)據(jù)較為接近，驗證了本文數(shù)值模擬方法的可靠性和準(zhǔn)確性。

圖3 DN板通道模擬與實驗j和f值

4.2 速度場和溫度場

選擇DN板通道內(nèi)4個代表性截面-、-、-和-（圖1c)），其中截面-、-垂直于軸，截面-垂直于斜波紋，截面-垂直于主流方向。截面-、-的速度矢量分布和溫度場模擬結(jié)果如圖4—圖6所示。為了更清晰地展示流動形態(tài)，速度矢量圖為局部放大后的結(jié)果。

圖4 N-N截面速度矢量分布

由圖4可見：-截面內(nèi)不同對應(yīng)的速度場均呈現(xiàn)出從左向右的整體流動趨勢；未發(fā)現(xiàn)如Ciofalo等人[12]、Zhang等人[14]在CC波紋通道中發(fā)現(xiàn)的渦流形態(tài)。這是由于DN板和CC板的結(jié)構(gòu)差異造成的，雖然DN板上也設(shè)有類似CC板的交錯斜波紋，但DN板上的斜波紋被尺寸更大的凹槽隔開而未能與相鄰的交錯斜波紋形成緊密有效接觸，大幅減弱了DN板通道內(nèi)交錯流體間的相互作用，從而不足以誘發(fā)渦流形態(tài)。由于缺少渦流形態(tài)，DN板通道內(nèi)的流體混合、邊界層擾動效果比CC板通道顯著減弱，導(dǎo)致DN板的傳熱系數(shù)、阻力系數(shù)均小于CC板，這與文獻[13]中的實驗結(jié)果一致。

圖5 O-O截面速度矢量分布

圖6 O-O截面溫度場示意

由圖5可見：-截面凹槽段二次流動幾乎不可見，凹槽段內(nèi)流動與平直通道內(nèi)流動相近，流股間的相互剪切作用很??；斜波紋段二次流動劇烈，斜波紋段中交錯流體間的相互剪切作用促進了相鄰板片斜波紋通道內(nèi)流體的相互遷移、混合，進而強化了傳熱。

由圖6可見：隨著增大，DN板通道近壁區(qū)的熱邊界層逐漸變薄，傳熱強度逐漸增大；不同對應(yīng)的工質(zhì)溫度分布特征相近，均為斜波紋段溫度高，凹槽段溫度低。如前文所述，斜波紋段傳熱能力比凹槽段更優(yōu)，工質(zhì)在斜波紋段從DN板壁吸收的熱量更多，故而溫度也更高。

-和-截面上的速度矢量分布、溫度場示意如圖7、圖8所示。圖7、圖8中計算工況對應(yīng)的工質(zhì)主流方向為從左向右。由圖7、圖8可見：-截面流動具有典型的波紋通道流動特征，工質(zhì)流經(jīng)該截面時對迎風(fēng)面施加周期性沖刷，所以此處的邊界層較薄，傳熱強度較高；-截面流動則與平直通道內(nèi)流動相近，對邊界層的擾動效果較差；-截面熱邊界層比-截面的薄。

圖7 L-L和M-M截面速度矢量分布

圖8 L-L和M-M截面溫度示意

4.3 努塞爾數(shù)（Nu）分布

圖9展示了計算域底部DN板上的努塞爾數(shù)（）分布情況，圖9中顏色由深到淺對應(yīng)值由低到高。由圖9可見：斜波紋和凹槽迎風(fēng)面上值較高，強化傳熱效果明顯，這是流體對壁面的沖刷作用造成的；斜波紋和凹槽背風(fēng)面以及DN板間接觸部位的值則較低；DN板上分布隨的增大而趨于均勻。Gaiser和Kottke[15]在對CC板的實驗研究中也發(fā)現(xiàn)了相似的分布現(xiàn)象。

圖9 DN板上Nu分布

為了深入分析對分布的影響，計算了沿DN板寬度方向（圖1中軸）的面積積分平均值分布情況，結(jié)果如圖10所示。圖10中將數(shù)值計算區(qū)域沿軸等分為20份，每份的面積積分平均努賽爾數(shù)為Nu，20個Nu中的最大值即為max。由圖10可見：值隨增大而增大，且趨于均勻；隨增大，DN板通道強化傳熱中心從凹槽段向斜波紋段轉(zhuǎn)移。較?。ā? 381）時，曲線中部（=9~13，凹槽段）出現(xiàn)波峰，凹槽段傳熱系數(shù)高于兩側(cè)斜波紋段，此時DN板通道中為層流流態(tài)，凹槽段流動阻力比斜波紋段小，流量較大，所以凹槽段傳熱系數(shù)更高。隨著的增大（≥3 288），斜波紋段強化傳熱效果逐漸顯現(xiàn)，其傳熱系數(shù)逐漸接近凹槽段，當(dāng)增至3 288時斜波紋段平均傳熱系數(shù)幾乎與凹槽段傳熱系數(shù)峰值持平，當(dāng)進一步增至5 252時甚至超過了后者?；谏鲜鼋Y(jié)果，建議將DN板運行工況控制在合適的雷諾數(shù)范圍內(nèi)，這樣既限制了流動損失又使DN板上的傳熱系數(shù)均勻分布，提高了DN板的利用率和使用壽命。

圖10 Re對Nu分布的影響

5 結(jié) 論

1）本文采用有限容積法和LBKE湍流模型，對回轉(zhuǎn)式空預(yù)器DN板通道內(nèi)空氣的流動傳熱過程進行了三維數(shù)值模擬研究。研究發(fā)現(xiàn)LBKE湍流模型可準(zhǔn)確反映DN板通道內(nèi)流動傳熱過程。

2）DN板通道內(nèi)交錯流體間的相互作用大為減弱，不足以誘發(fā)渦流形態(tài)，DN板通道流場中并未發(fā)現(xiàn)渦流結(jié)構(gòu)。DN板通道內(nèi)流動包含波紋通道流動和平直通道流動，其中前者的強化傳熱效果明顯，對整體傳熱貢獻更大。

3）DN板通道斜波紋和凹槽的迎風(fēng)面值較高，斜波紋和凹槽的背風(fēng)面及DN板間接觸部位值較低。隨增大而增大，分布隨增大而更趨均勻。

4）建議DN板在約為5 000的條件下運行，這樣既能保證DN板上均勻分布又控制了阻力損失。

［1］ 應(yīng)靜良, 李永華. 電站鍋爐空氣預(yù)熱器[M]. 北京: 中國電力出版社, 2002: 9-10. YING Jingliang, LI Yonghua. Air preheater for power station boiler[M]. Beijing: China Electric Power Press, 2002: 9-10.

［2］ 張續(xù)墩. 回轉(zhuǎn)式空預(yù)器換熱板數(shù)值模擬及優(yōu)化[D]. 北京: 華北電力大學(xué), 2016: 1-2. ZHANG Xudun. Numerical simulation and optimization on rotary air preheater[D]. Beijing: North China Electric Power University, 2016: 1-2.

［3］ 王和平, 程紹兵. 熱電聯(lián)產(chǎn)鍋爐機組設(shè)備及經(jīng)濟運行[M]. 北京: 中國石化出版社, 2014: 216-217. WANG Heping, CHENG Shaobing. Cogeneration boiler unit equipment and economic operation[M]. Beijing: China Petrochemical Press, 2014: 216-217.

［4］ STASIEK J, COLLINS M W, CIOFALO M, et al. Investigation of flow and heat transfer in corrugated passages-I. Experimental results[J]. International Journal of Heat and Mass Transfer, 1996, 39(1): 149-164.

［5］ CIOFALO M, COLLINS M W, STASIEK J A. Flow and heat transfer predictions in flow passages of air preheaters: assessment of alternative modeling approaches[M] //SUNDéN B, FAGHRI M. Computer Simulations in Compact Heat Exchangers. Southampton, United Kingdom,Boston, MA, USA: Computational Mechanics Publications, 1998: 169-225.

［6］ 劉建民, 陳國慶, 黃啟龍, 等. 燃煤脫硝機組空氣預(yù)熱器蓄熱片表面飛灰沉積板結(jié)機理研究[J]．中國電機工程學(xué)報, 2016, 36(增刊1): 132-139. LIU Jianmin, CHEN Guoqing, HUANG Qilong, et al. Study on mechanism of fly ash deposition and hardening on the air preheater regenerative piece surface of the coal-fired and denitration unit[J]. Proceedings of the CSEE, 2016, 36(Suppl.1): 132-139.

［7］ 王利民, 樊強, 鄧?yán)? 等. DN板通道傳熱阻力特性的實驗和數(shù)值研究[J]. 工程熱物理學(xué)報, 2015, 36(10): 2250-2254.WANG Limin, FAN Qiang, DENG Lei, et al. Experimentaland numerical studies of thermal-hydraulic performance of DN plate channel[J]. Journal of Engineering Thermophysics, 2015, 36(10): 2250-2254.

［8］ 危日光, 高建強, 張磊, 等. 回轉(zhuǎn)式空氣預(yù)熱器搪瓷蓄熱元件傳熱與流動數(shù)值模擬研究[J]. 華北電力大學(xué)學(xué)報, 2016, 43(2): 87-92.WEI Riguang, GAO Jianqiang, ZHANG Lei, et al. Numerical simulation on heat transfer and flow of enamel heat storage elements of rotary air preheater[J]. Journal of North China Electric Power University, 2016, 43(2): 87-92.

［9］ 黃風(fēng)良, 孫志堅, 李鵬程, 等. 帶擾流孔波紋板的傳熱和阻力特性[J]. 浙江大學(xué)學(xué)報, 2015, 49(7): 1242-1248.HUANG Fengliang, SUN Zhijian, LI Pengcheng, et al. Heat transfer and resistance characteristics of corrugated plate with spoiler holes[J]. Journal of Zhejiang University, 2015, 49(7): 1242-1248.

［10］ PATEL V C, RODI W, SCHEUERER G. Turbulence models for near-wall and low Reynolds number flows: a review[J]. AIAA Journal, 1985, 23(9): 1308-1319.

［11］ LAM C K G, BREMHORST K. A modified form of the-model for predicting wall turbulence[J]. ASME Journal of Fluids Engineering, 1981, 103(3): 456-460.

［12］ CIOFALO M, STASIEK J, COLLINS M W. Investigation of flow and heat transfer in corrugated passages-II. numerical simulations[J]. International Journal of Heat and Mass Transfer, 1996, 39(1): 165-192.

［13］ 張磊, 聶鵬, 劉文娟, 等. 復(fù)合式波紋板傳熱與阻力特性單吹瞬態(tài)實驗研究[J]. 中國電機工程學(xué)報, 2017, 37(22): 6675-6680. ZHANG Lei, NIE Peng, LIU Wenjuan, et al. Single-blow transient test on the thermal-hydraulic performance of the composite corrugated plate[J]. Proceedings of the CSEE, 2017, 37(22): 6675-6680.

［14］ ZHANG L, CHE D F. Turbulence models for fluid flow and heat transfer between cross-corrugated plates[J]. Numerical Heat Transfer, Part A: Applications, 2011, 60(5): 410-440.

［15］ GAISER G, KOTTKE V. Flow phenomena and local heat and mass transfer in corrugated passages[J]. Chemical Engineering & Technology, 1989, 12(1): 400-405

Numerical investigation on heat transfer and resistance characteristics of composite corrugated plate

ZHANG Lei1, DU Juan2, BAI Wengang1, GAO Wei1, WU Shuaishuai1

(1. Xi’an Thermal Power Research Institute Co., Ltd., Xi’an 710054, China; 2. College of Natural Resources and Environment, Northwest A&F University, Yangling 712100, China)

The double notched (DN) type corrugated plate heat transfer element is commonly used in rotary air preheaters. To reveal the micro-flow morphology and strengthening heat transfer mechanism in the corrugated channel of the DN plate, the LBKE (Lam and Bremhorst k-ε model) turbulence model and CFD numerical simulation technology are applied to investigate the heat transfer and flow resistance characteristics of the DN plate, including predicting the heat transfer and flow resistance coefficient, and analyzing the flow field and temperature field inside the channel of the DN plate. The results show that, two types of flow are observed in the predicted velocity fields, wavy channel flow and pipe flow. The Nusselt number Nu in upwind surface of the DN plate is larger while that at the leeward side is lower. As the Reynolds number Re increases, the Nu increases and its distribution on the DN plate becomes more uniform. Increasing the proportion of the corrugated structure can enhance the heat transfer capacity of the DN plate, but will also increase the resistance loss simultaneously. It suggests that the DN plate operate at moderate Re values (≈5 000), to obtain the best effect and ensure the uniform distribution of Nu and control the resistance loss.

rotary air preheater, DN plate, corrugated plate, numerical simulation, heat transfer, friction factor

TK124

A

10.19666/j.rlfd.201810199

張磊, 杜娟, 白文剛, 等. 復(fù)合式波紋板傳熱與阻力特性數(shù)值研究[J]. 熱力發(fā)電, 2019, 48(5): 44-48. ZHANG Lei, DU Juan, BAI Wengang, et al. Numerical investigation on heat transfer and resistance characteristics of composite corrugated plate[J]. Thermal Power Generation, 2019, 48(5): 44-48.

2018-10-31

張磊(1982—)，男，博士，高級工程師，主要研究方向為強化傳熱、換熱器優(yōu)化設(shè)計，zhangleix@tpri.com.cn。

（責(zé)任編輯 楊嘉蕾）