涂福炳，馬士偉，高晟揚(yáng)，劉慶赟，毛建豐

(中南大學(xué) 能源科學(xué)與工程學(xué)院，湖南 長沙，410083)

傳統(tǒng)的管殼式換熱器傳熱系數(shù)低、易結(jié)垢、易腐蝕、占地面積大且金屬耗量大，已不能滿足節(jié)能降耗和系統(tǒng)優(yōu)化的要求。而徑向熱管換熱器由于傳熱效率高、等溫性好、壓降小、結(jié)構(gòu)簡單、緊湊性好、抗露點(diǎn)腐蝕等優(yōu)點(diǎn)而廣受青睞[1]，已廣泛應(yīng)用于工業(yè)余熱的回收利用。Cotter[2]于 1965年提出熱管理論后，熱管的研究一直受到很大關(guān)注。Lee等[3]通過實(shí)驗(yàn)研究了影響熱管傳熱性能的各種內(nèi)外因素；Peretz等[4]從參變量的非線性特征出發(fā)，對熱管換熱器幾何結(jié)構(gòu)進(jìn)行了優(yōu)化；孫世梅等[5?6]借助計(jì)算流體力學(xué)軟件對徑向熱管換熱器內(nèi)流體的流動(dòng)和傳熱進(jìn)行了數(shù)值模擬研究；陳潔等[7]以經(jīng)濟(jì)效益為目標(biāo)函數(shù)對用于回收空調(diào)排風(fēng)熱量的熱管換熱器進(jìn)行了研究，得到了不同溫度下?lián)Q熱器的最優(yōu)結(jié)構(gòu)；許欣[8]通過引入多孔介質(zhì)模型和分布阻力方法，得到了徑向熱管換熱器在來流速度均勻情況下壓力分布特征以及壓降公式。在此，本文作者應(yīng)用數(shù)值模擬方法研究來流速度分布對換熱器流場、溫度場和整體性能的影響規(guī)律，以便為徑向熱管換熱器工程應(yīng)用提供參考。

1 模型的建立

1.1 物理模型

徑向熱管換熱器幾何結(jié)構(gòu)和殼程流場比較復(fù)雜，為了得到研究對象的殼程流場和壓力場，同時(shí)考慮到換熱器殼程煙氣流動(dòng)的實(shí)際情況，對研究對象進(jìn)行如下簡化：(1) 不考慮煙氣中煙塵顆粒的影響，將煙氣視為多質(zhì)單相流體；(2) 煙氣物理性質(zhì)穩(wěn)定，不發(fā)生相變和化學(xué)變化等；(3) 將煙氣視為各向同性、均勻連續(xù)性的不可壓縮牛頓流體；(4) 將煙氣流動(dòng)視為穩(wěn)態(tài)流動(dòng)；(5) 不考慮熱管污垢對換熱的影響。熱管換熱器結(jié)構(gòu)參數(shù)如表1 所示。

換熱器計(jì)算區(qū)域如圖1所示?？紤]到熱管的對稱性，只需對1?1和2?2這2個(gè)對稱面之間的區(qū)域進(jìn)行模擬。對于換熱器的網(wǎng)格劃分，以徑向熱管單元為中心將模型劃分成眾多的小體，然后，對這些小體分別進(jìn)行網(wǎng)格劃分；利用六面體網(wǎng)格來劃分徑向熱管單元及翅片部分，用三棱柱網(wǎng)格劃分換熱器殼程空間。圖2所示為計(jì)算區(qū)域局部三維網(wǎng)格圖。

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

數(shù)學(xué)模型包括質(zhì)量方程、動(dòng)量方程、能量方程以及氣體狀態(tài)方程。這些模型是較普遍的計(jì)算模型，具體方程式見文獻(xiàn)[9]。標(biāo)準(zhǔn)k?ε雙方程湍流模型可用于壁面邊界層、管流、無旋及弱旋的二維或三維流動(dòng)、通道流等，在工程實(shí)際及科學(xué)研究中得到了廣泛的應(yīng)用。本文運(yùn)用標(biāo)準(zhǔn)k?ε雙方程模型的輸運(yùn)方程模擬徑向熱管換熱器殼程煙氣的流動(dòng)。

表1 徑向熱管換熱器的結(jié)構(gòu)參數(shù)Table1 Structural parameters of radial heat pipe heat exchanger

圖1 換熱器計(jì)算區(qū)域示意圖Fig.1 Schematic diagram of heat exchanger computational domain

圖2 熱管換熱器局部三維網(wǎng)格圖Fig.2 Local three-dimensional mesh of radial heat pipe heat exchanger

1.3 模型設(shè)置求解及邊界條件

采用穩(wěn)態(tài)、定常模型和PRESTO壓力插補(bǔ)格式；利用二階迎風(fēng)格式離散控制方程對流項(xiàng)；采用SIMPLE算法求解離散方程；采用標(biāo)準(zhǔn)壁面函數(shù)法處理近壁面煙氣流動(dòng)。

采用速度入口邊界條件，以現(xiàn)場運(yùn)行的徑向熱管換熱器為研究對象。由于煙氣進(jìn)入換熱器前在煙道中進(jìn)行了轉(zhuǎn)向，使得換熱器入口速度分布比較復(fù)雜。通過對其流場進(jìn)行模擬，可得到換熱器入口煙氣速度分布，速度擬合式為(不考慮煙氣沿?zé)峁茌S向流動(dòng))：

入口煙氣溫度為185 ℃，流量為128 203 m3/h。煙氣成分見表 2。入口煙氣湍流采用水力直徑和湍流強(qiáng)度。

表2 煙氣成分(體積分?jǐn)?shù))Table2 Components of flue gas %

出口邊界采用自由流出口邊界，近壁區(qū)采用壁面函數(shù)法處理。面1?1和2?2設(shè)置為對稱邊界；換熱器殼側(cè)壁面設(shè)置為無滑移、無滲透及絕熱邊界；熱管內(nèi)壁面設(shè)置為定溫；外壁面設(shè)定為無滑移，無滲透；選擇Coupled流固耦合條件。

采用等效導(dǎo)熱系數(shù)法對熱管工質(zhì)相變過程進(jìn)行簡化[10?11]。熱管等效導(dǎo)熱系數(shù)由下式[12]計(jì)算：

式中：leff為等效導(dǎo)熱系數(shù)，W/(m·K)；R為充液率，%，文中取R=43.2%。

2 模型驗(yàn)證

換熱器殼程不同點(diǎn)溫度測試結(jié)果與模擬結(jié)果對比如表3所示。

表3中4個(gè)測點(diǎn)均布于熱管換熱器的中心線上。從表3可以看出：數(shù)值模擬結(jié)果與測試結(jié)果有相同的變化趨勢，各測點(diǎn)溫度相對誤差小于 5%，說明所建模型合理，結(jié)果可靠。

表3 殼程溫度模擬結(jié)果與測試結(jié)果比較Table 3 Comparison between simulation results and test results of temperature in shell side

3 來流速度分布對換熱器性能影響的數(shù)值模擬

速度包含方向與大小2個(gè)特征，可以從入射角和均勻度2方面設(shè)計(jì)來流速度分布方案，模擬分析換熱器性能變化情況，以便為改善換熱器整體性能提供依據(jù)。

3.1 來流入射角方案

煙氣入射角φ定義為：在X?Y平面內(nèi)，煙氣速度與豎直方向所成的角度。方案設(shè)計(jì)以換熱器入口煙氣質(zhì)量流量不變?yōu)榛A(chǔ)，且換熱器入口處速度均勻，入射角方案如表4所示。

表4 煙氣入射方案Table 4 Schemes for smoke inclination angles

3.2 來流速度均勻度分布方案

在來流平均速度相等的基礎(chǔ)上，引入速度不均勻因子τ衡量來流速度的不均勻程度，其定義式為[13]

式中：n為控制單元數(shù)；S為換熱器迎風(fēng)面面積，m2；Si為換熱器第i個(gè)單元迎風(fēng)面面積，m2；vi為第i個(gè)單元進(jìn)口速度，m/s；為換熱器入口截面平均速度，m/s。

在入射角為0°的前提下，選取3種典型的階梯型即二階階梯型、三階階梯型和拋物階梯型速度分布進(jìn)行模擬計(jì)算，各階梯型速度分布如表5所示。從表5可見：方案序號越大，速度分布越不均勻，因子τ越大；各階梯型速度分布中的分速度在換熱器入口X軸方向上的寬度相等。

3.3 模擬結(jié)果分析

3.3.1 對流場的影響

圖3所示為來流速度均勻且入射角φ為0°時(shí)換熱器殼程局部速度矢量圖。由圖3可知：換熱器殼程流場較均勻，煙氣流無明顯偏移。前排管束的尾跡與單圓柱繞流的尾跡非常相似，在管束后形成穩(wěn)定的楔形雙渦旋回流區(qū)。沿?zé)煔饬飨颍捎诠苁母蓴_，煙氣流湍流度增加，熱管背風(fēng)側(cè)產(chǎn)生渦脫落現(xiàn)象。

圖4所示為來流速度均勻且入射角φ為45°時(shí)換熱器殼程局部速度矢量圖。由圖4可知：由于煙氣流斜向進(jìn)入換熱器，熱管背風(fēng)側(cè)回流區(qū)發(fā)生偏移，下游管束被上游管束尾跡覆蓋，出現(xiàn)類似于順排管束的流動(dòng)狀態(tài)，減少了煙氣與管束的碰撞以及煙氣壓力損失，換熱效率有所下降；此外，在換熱器入口左側(cè)有一范圍較小的煙氣滯流區(qū)，且隨來流速度的增大，滯流區(qū)范圍增大。

表5 不同階梯型速度分布Table 5 Different step velocity distributions m/s

圖3 φ=0°, Z=1.75 mm且來流速度均勻時(shí)局部速度矢量圖Fig.3 Local velocity vectors with uniform velocity inflow when φ=0° and Z=1.75 mm

圖4 φ=45°, Z=1.75 mm且來流速度均勻時(shí)局部速度矢量圖Fig.4 Local velocity vector with uniform velocity inflow when φ=45° and Z=1.75 mm

圖5所示為方案4的三階階梯型速度分布中換熱器殼程速度矢量圖。由圖5可知：由于來流速度分布不均勻，換熱器殼程出現(xiàn)明顯的高速區(qū)與低速區(qū)，且不均勻度越大，高速區(qū)與低速區(qū)越明顯；在不同速度交界處，煙氣由高速區(qū)向低速區(qū)偏移。煙氣流的偏移使熱管背風(fēng)側(cè)回流區(qū)發(fā)生偏移，同樣出現(xiàn)類似于順排管束的流動(dòng)狀態(tài)，且速度不均勻程度越大，煙氣產(chǎn)生偏移范圍越大；沿殼程流向，煙氣偏移逐漸減弱，速度逐漸均勻；此外，煙氣經(jīng)過管束間隙時(shí)速度越大，管束的磨損、振動(dòng)越嚴(yán)重，故應(yīng)避免換熱器入口局部速度過大的情況出現(xiàn)。

3.3.2 對溫度場的影響

圖6所示為來流速度均勻且入射角φ為0°時(shí)換熱器殼程溫度分布云圖。由圖6可知：當(dāng)來流速度均勻時(shí)，不同位置的煙氣經(jīng)過每排管束的溫降基本相同，煙氣在同一橫截面上溫度較均勻；靠近換熱器殼壁的煙氣溫度較高，這是換熱器殼壁設(shè)置為絕熱邊界的緣故。

圖5 φ=0°, Z=1.75 mm且來流速度為三階階梯型(方案4)時(shí)速度分布云圖Fig.5 Velocity profile with third-order type of velocity inflow for Scheme 4 when φ=0° and Z=1.75 mm

圖6 φ=0°, Z=1.75 mm且來流速度均勻時(shí)溫度分布云圖Fig.6 Temperature profile with uniform velocity inflow when φ=0° and Z=1.75 mm

圖7所示為方案4的拋物階梯型速度分布中換熱器殼程溫度分布云圖。由圖7可知：由于來流速度分布不均勻，換熱器溫度場均勻性變差，煙氣在來流煙氣速度小的區(qū)域停留時(shí)間較長，與管束換熱時(shí)間也較長，出現(xiàn)低溫區(qū)。在同一階梯型分布中，因子τ越大，溫度場均勻性越差，低溫區(qū)范圍也越大，此時(shí)，若低溫區(qū)溫度低于煙氣露點(diǎn)，將使換熱器產(chǎn)生局部腐蝕。

圖7 φ=0°, Z=1.75 mm且來流速度為拋物階梯型(方案4)時(shí)溫度分布云圖Fig.7 Temperature profile with parabolic type of velocity inflow for Scheme 4 when φ=0° and Z=1.75 mm

3.3.3 對換熱器整體性能的影響

在評價(jià)換熱器整體性能時(shí)，換熱量與壓降往往是要兼顧的，有多種將換熱量與壓降結(jié)合的評價(jià)方法，如單位壓降換熱系數(shù)法、熵分析法、熱經(jīng)濟(jì)學(xué)分析法等[14]。實(shí)踐證明：將單位壓降換熱系數(shù)α/?p作為衡量換熱器綜合性能的指標(biāo)能獲得理想效果[15]；α/?p越大，表明換熱器整體整體性能越好。

圖8 速度均勻時(shí)單位壓降換熱系數(shù)α/?p 隨入射角變化情況Fig.8 Relationship between α/?p and inclination angle when velocity inflow is uniform

圖 8所示為α/?p隨入射角的變化情況。由圖 8可知：隨著煙氣入射角的增大，α/?p逐漸增大，并在入射角約為 45°時(shí)達(dá)到最大，此時(shí)，換熱器整體性能最好，提高幅度為2.31%；此后，α/?p隨入射角的進(jìn)一步增大急劇減小，換熱器性能惡化。

當(dāng)入射角為0°時(shí)，α/?p隨τ因子的變化情況如圖9所示。從圖9可以看出：在二、三階階梯型下，隨因子τ的增大，α/?p逐漸減小且下降幅度逐漸增大，換熱器整體性能惡化；在拋物階梯型下，α/?p隨τ的增大先逐漸增大，并在τ為20%左右時(shí)達(dá)到最大；此后，隨τ的繼續(xù)增大，α/?p逐漸減小。經(jīng)比較可知：二階階梯型速度分布下?lián)Q熱器性能惡化最嚴(yán)重，α/?p最大降幅為11.88%；其次是三階階梯型，最大降幅為5.14%；而拋物階梯型在τ較小時(shí)提高了換熱器整體性能。因此，當(dāng)τ較小時(shí)，多股擾流能提高換熱器整體性能；而少股擾流將削弱換熱器整體性能，且τ越大，整體性能越差。

圖9 入射角為0°時(shí)單位壓降換熱系數(shù)α/?p隨不均勻因子τ變化情況Fig.9 Relationship between α/?p and τ when inclination angle is 0°

4 結(jié)論

(1) 當(dāng)來流速度均勻時(shí)，入射角增大使熱管背風(fēng)側(cè)回流區(qū)發(fā)生偏移，下游管束被上游管束尾跡覆蓋，減少了煙氣與管束的碰撞，降低了煙氣壓力損失，換熱效率有所降低；當(dāng)來流速度分布不均勻時(shí)，換熱器殼程煙氣出現(xiàn)明顯的高速區(qū)與低速區(qū)，煙氣產(chǎn)生偏移，且不均勻度越大，高速區(qū)與低速區(qū)越明顯，煙氣偏移范圍也越大。

(2) 當(dāng)來流速度均勻時(shí)，同一橫截面上煙氣溫度分布較均勻；當(dāng)來流速度不均勻時(shí)，換熱器溫度場均勻性變差，來流速度小的區(qū)域煙氣停留時(shí)間較長，與管束換熱時(shí)間也較長，出現(xiàn)低溫區(qū)，且因子τ越大，溫度場均勻性越差，低溫區(qū)范圍也越大。

(3) 單位壓降換熱系數(shù)α/?p隨煙氣入射角的增大而逐漸增大，并在入射角為 45°時(shí)達(dá)到最大，此時(shí)換熱器整體性能最好。此后，α/?p隨入射角的增大急劇減小，換熱器性能惡化明顯。

(4) 在τ較小時(shí)，多股擾流能提高換熱器整體性能；而少股擾流將降低換熱器整體性能，且煙氣流不均勻度越大，整體性能惡化越嚴(yán)重。

[1] 毛建豐, 涂福炳, 許欣, 等. 同軸徑向熱管換熱器殼程流場數(shù)值模擬[J]. 熱科學(xué)與技術(shù), 2010, 9(1): 31?35.MAO Jianfeng, TU Fubing, XU Xin, et al. Numerical simulation on flow field of jacketed heat pipe heat exchanger[J]. Journal of Thermal Science and Technology, 2010, 9(1): 31?35.

[2] Cotter T P. Theory of heat pipes[R]. Los Alamos: Los Alamos Scientific Lab, 1965: 1?37.

[3] Lee Y, Mital U. A two-phase closed thermosyphon[J].International Journal of Heat and Mass Transfer, 1972, 15(9):1695?1707.

[4] Peretz R, Horbaniuc B. Optimal heat pipe heat exchanger design[J]. Journal of Heat Recovery Systems, 1984, 4(1): 9?24.

[5] 孫世梅, 張紅. 熱管換熱器流動(dòng)與傳熱的CFD模擬及試驗(yàn)[J].南京工業(yè)大學(xué)學(xué)報(bào), 2004, 26(2): 62?66.SUN Shimei, ZHANG Hong. Analysis of CFD simulation with experiment of heat transfer and pressure drop for heat pipe heat exchanger[J]. Journal of Nanjing University of Technology, 2004,26(2): 62?66.

[6] 孫世梅, 張紅. 熱管換熱器傳熱性能及溫度場數(shù)值模擬[J].化工學(xué)報(bào), 2004, 55(3): 472?475.SUN Shimei, ZHANG Hong. Numerical simulation of thermal performance and temperature field in heat pipe heat exchanger[J].Journal of Chemical Industry and Engineering, 2004, 55(3):472?475.

[7] 陳潔, 曹家樅. 空調(diào)排風(fēng)能量回收用熱管換熱器的優(yōu)化設(shè)計(jì)[J]. 制冷與空調(diào), 2007, 21(3): 44?47.CHEN Jie, CAO Jiazong. Optimum design for heat pipe heat exchanger in air-conditioning ventilation system[J].Refrigeration and Air-Conditioning, 2007, 21(3): 44?47.

[8] 許欣. 夾套式熱管傳熱特性試驗(yàn)研究及其換熱器殼程壓力場數(shù)值模擬[D]. 長沙: 中南大學(xué)能源科學(xué)與工程學(xué)院, 2009:1?55.XU Xin. Study on heat transfer characteristics of the jacketed heat pipe and numerical simulation on flow field of the exchanger[D]. Changsha: Central South University. School of Energy Science and Engineering, 2009: 1?55.

[9] 周萍, 周乃君, 蔣愛華, 等. 傳遞過程原理及其數(shù)值仿真[M].長沙: 中南大學(xué)出版社, 2006: 299?301.ZHOU Ping, ZHOU Naijun, JIANG Aihua, et al. Principle and numerical stimulation of transport process[M]. Changsha:Central South University Press, 2006: 299?301.

[10] Zhang C, Sousa A C M, Venart J E S. The numerical and experimental study of a power plant condenser[J]. Heat Transfer,1993, 115(2): 435?445.

[11] Jacobi A M, Shall R K. Heat transfer surface enhancement through the use of longitudinal vortices: A review of recent progress[J]. Experimental Thermal and Fluid Science, 1995,11(3): 295?309.

[12] 張嶺. 夾套式熱管換熱器多場數(shù)值模擬與優(yōu)化研究[D]. 長沙:中南大學(xué)能源科學(xué)與工程學(xué)院, 2010: 1?72.ZHANG Ling. Numerical simulation and optimum research on jacket-type heat pipe heat exchanger[D]. Changsha: Central South University. School of Energy Science and Engineering,2010: 1?72.

[13] 鄧斌, 陶文銓, 余汪洋. 來流不均勻性對冷凝器性能影響的數(shù)值模擬研究[J]. 制冷與空調(diào), 2007, 21(3): 1?6 DENG Bin, TAO Wenquan, YU Wangyang. Numerical investigation of condensers with non-uniform air distribution[J].Refrigeration and Air-Conditioning, 2007, 21(3): 1?6.

[14] 史美中, 王中錚. 熱交換器原理與設(shè)計(jì)[M]. 南京: 東南大學(xué)出版社, 2005: 275?278.SHI Meizhong, WANG Zhongzheng. Principia and design of heat transfer device[M]. Nanjing: Southeast University Press,2005: 275?278.

[15] 王玉琴, 宋天民, 張國福. 橢圓管管殼式換熱器殼程傳熱性能的試驗(yàn)研究[J]. 壓力容器, 2008, 25(10): 11?15.WANG Yuqin, SONG Tianmin, ZHANG Guofu. Experiment research on heat transfer performance in shell of the shell-and-tube heat exchanger with elliptic tube[J]. Pressure Vessel Technology, 2008, 25(10): 11?15.