摘 " " "要： 換熱器內(nèi)折流板是影響其換熱效率的關(guān)鍵部件，對折流板進(jìn)行優(yōu)化設(shè)計可以有效提高換熱器的換熱效率，從而提高能源利用率，節(jié)約能源。對新型L型折流板進(jìn)行了研究，利用CFD原理，建立了L型折流板換熱器的幾何模型，并與普遍的單弓形折流板進(jìn)行對比分析。結(jié)果表明L型折流板橫板長度為80 mm時其綜合性能最佳，壓力損失較單弓形折流板降低了11%～13%，當(dāng)流速在1～3 m/s時，換熱效率提高了7.29%；當(dāng)流速大于3 m/s時，最大降幅為5.18%，且當(dāng)流速為1 m/s、3 m/s和5 m/s時較單弓形折流板換熱器綜合性能分別提高了12.3%、12.4%和13.7%，為換熱器工業(yè)設(shè)計提供了參考。

關(guān) "鍵 "詞：換熱器；折流板；優(yōu)化設(shè)計；換熱效率；綜合性能

中圖分類號：TQ013 " " 文獻(xiàn)標(biāo)識碼： A " " 文章編號： 1004-0935（2023）07-0983-06

目前，人類社會活動所需要的能源80%以上來自傳統(tǒng)化石能源[1]。據(jù)統(tǒng)計，中國對石油的需求量占世界的30.4%左右，到2024年，我國對煤炭的需求缺口將超過3.6億噸[2]。為了有效解決化石能源不斷減少及能源浪費而造成的世界“能源危機(jī)”問題，國內(nèi)外學(xué)者都致力于對節(jié)能高效新型技術(shù)的研究。優(yōu)化能源利用結(jié)構(gòu)，提高能源利用率，通過不斷地技術(shù)創(chuàng)新來解決迫在眉睫的能源問題至關(guān)重要[3]。換熱器是化工、石油、能源、海水淡化以及冶金等領(lǐng)域應(yīng)用最為廣泛的設(shè)備之一，在石油化工行業(yè)，換熱器的經(jīng)濟(jì)投入更是占整個行業(yè)42%以上[4]。以提高設(shè)備工作效率、降低燃料消耗、提高能源利用率、減小尺寸以及降低工業(yè)設(shè)備的投資成本為目的，對該設(shè)備進(jìn)行優(yōu)化設(shè)計，一直都是國內(nèi)外學(xué)者研究的方向。

迄今為止，已有換熱器種類諸多，其中管殼式換熱器在工業(yè)中被廣泛應(yīng)用[5]，其在世界范圍內(nèi)市場份額占36.9%左右，在石油化工領(lǐng)域市場份額占據(jù)70%左右[6]。管殼式換熱器主要由殼體、換熱管、折流板、接頭、管箱封頭等結(jié)構(gòu)組成。其中折流板是影響管殼式換熱器換熱效率的關(guān)鍵部件，折流板能夠使殼程流體沖刷管束次數(shù)增加，減少流動死區(qū)，提高流速從而提高工作效率，但折流板過多會導(dǎo)致?lián)Q熱器殼程壓力損失過大，導(dǎo)致設(shè)備壽命急劇減小。因此，對折流板進(jìn)行優(yōu)化設(shè)計從而有效提高換熱器的換熱效率，增長換熱器壽命，提高能源利用率，節(jié)約能源十分必要。

近年來，國內(nèi)外學(xué)者對換熱器的研究較為普遍，王志鵬等[7]提出了新型垂直布置百葉折流板管殼式換熱器，邵海磊等[8]提出新型的四分螺旋折流板換熱器，王英雙研發(fā)了一種新型的花格板管殼式換熱器，美國ABB公司研發(fā)了螺旋折流板換熱器，Mohebbi等[9]對帶有波紋的小型板式換熱器進(jìn)行研究，喻九陽等[10]對不同開孔數(shù)量的折流板進(jìn)行了實驗研究等。但是仍然存在設(shè)備工藝復(fù)雜、成本巨大、不易大規(guī)模生產(chǎn)使用、設(shè)備使用壽命短[11]、設(shè)備換熱工作效率低[12]、損失大、能源利用率不夠高等問題。針對該領(lǐng)域的現(xiàn)狀，本文基于傳統(tǒng)單弓形折流板換熱器，選取新型L型折流板進(jìn)行模擬研究。綜合運用Fluent軟件和Post軟件等對不同橫板長度的L型折流板管殼式換熱器進(jìn)行仿真模擬，通過壓力、溫度等與傳統(tǒng)單弓形折流板換熱器的數(shù)值相比較，最終確定最佳橫板長度，為換熱器工業(yè)設(shè)計提供了參考依據(jù)。

1 "數(shù)值模擬方法及數(shù)值計算方法

1.1 "數(shù)值模擬方法

首先要建立模型，主要為幾何模型，幾何模型是對描述現(xiàn)象或物體建立物理或數(shù)學(xué)模型。其次是要進(jìn)行網(wǎng)格劃分，將幾何模型劃分為多個小單元，這些小的單元存在著大量的數(shù)據(jù)信息，通過單元接觸的節(jié)點建立方程計算，從而形成有限元模型。接著要確定邊界條件，準(zhǔn)確的設(shè)置物性參數(shù)、邊界條件，選擇正確的計算方法，對計算的準(zhǔn)確性至關(guān)重要；最后要模擬計算并查看結(jié)果，結(jié)果可以通過圖像、表格、曲線等表示，同時也可導(dǎo)出平面的流場等信息。

1.2 "數(shù)值計算方法

1.2.1 "研究對象

本文研究對象為新型L型折流板，建立簡化后的幾何模型與合理的有限元模型，同時使用Fluent軟件模擬不同流速時單弓形、L型折流板管殼式換熱器的流場分布，得到殼程壓力、溫度場云圖[13]。利用Post軟件得出進(jìn)出口壓降、溫差，進(jìn)行分析。對不同橫板長度的L型折流板管殼式換熱器進(jìn)行仿真模擬，通過壓力、溫度較傳統(tǒng)單弓形折流板換熱器的數(shù)值相比較，評定最佳橫板長度。研究后發(fā)現(xiàn)，L型折流板換熱器結(jié)構(gòu)比傳統(tǒng)單弓形折流板結(jié)構(gòu)更穩(wěn)定，壓力損失更小[14]。橫板長度為80 mm時，壓力損失較單弓形折流板降低了11%～13%，當(dāng)流速在1 ～3 m/s之間時，換熱效率提高了7.29%；當(dāng)流速大于3 m/s時，最大降幅為5.18%。

1.2.2 "控制方程

本文中流體不可壓縮，無內(nèi)熱源，將標(biāo)準(zhǔn)k-ε模型簡化，簡化后的k-ε模型湍流方程如下：

（1）

（2）

物體各部分之間或相互接觸的兩個幾個物體，不發(fā)生宏觀位移，僅靠分子、原子、自由電子之間發(fā)生相對運動而產(chǎn)生的傳熱現(xiàn)象稱為導(dǎo)熱。根據(jù)傅里葉定律，單位時間內(nèi)通過平面的導(dǎo)熱能量與該處溫度梯度、平板面積成正比，即導(dǎo)熱方程：

（3）

流體由于宏觀位移，使得冷熱流體互相摻混時所導(dǎo)致的熱量交換過程稱為對流換熱過程，其基本關(guān)系式為牛頓冷卻公式：

^_D_Dd_______

物體之間通過電磁波來傳遞能量稱為輻射換熱。實際物體的斯特藩-玻爾茲曼定律如下：

（5）

動量守恒方程：

（6）

（7）

（8）

1.2.3 "網(wǎng)格劃分

本文通過ANSYS Workbench中的mesh模塊選取四面體非結(jié)構(gòu)化網(wǎng)格繪制[15]。網(wǎng)格劃分的好壞、多少對于計算準(zhǔn)確性、計算時間成本都至關(guān)重要。網(wǎng)格劃分可分為結(jié)構(gòu)化網(wǎng)格和非結(jié)構(gòu)化網(wǎng)格。非結(jié)構(gòu)化網(wǎng)格有三角形、四邊形或兼具兩種，也有六面體、四面體、三棱柱形以及它們的組合。對熱量變化明顯的地方應(yīng)該適當(dāng)提高網(wǎng)格密度，以提高計算精度。同時應(yīng)該在熱量變化平穩(wěn)的部分減少網(wǎng)格，以減少離散單元數(shù)，節(jié)約計算資源。同時要求不同網(wǎng)格區(qū)域結(jié)點相互連接。本文中L型折流板換熱器網(wǎng)格模型如圖1所示。

劃分網(wǎng)格結(jié)束后，需檢查網(wǎng)格質(zhì)量，確保滿足計算要求。ANSYS mesh中提供了多種網(wǎng)格質(zhì)量檢測項目，如正交質(zhì)量、正交斜歪率、最大長寬比等，本文選取Element Quality，該值在0～1之間，數(shù)值越大網(wǎng)格越優(yōu)良，本文模型平均網(wǎng)格質(zhì)量都在0.82以上。因此網(wǎng)格質(zhì)量滿足計算要求。本文中L型折流板換熱器網(wǎng)格質(zhì)量如圖2所示。

1.2.4 "邊界條件

在繪制不同折流板換熱器幾何模型與有限元模型后，在ANSYS Fluent 2020軟件進(jìn)行計算之前，需要在軟件中確定邊界條件。

計算時，液態(tài)水由殼程入口流入換熱器，機(jī)油從管程入口流入換熱器，邊界條件如下，

速度入口：殼程入口速度分別為1 m/s、3 m/s和5 m/s，管程入口速度為0.5 m/s。

入口溫度：液態(tài)水為20 ℃，機(jī)油為100 ℃。

壓力出口：管程、殼程壓力出口邊界靜壓為零，參考?xì)鈮簽?01 kPa。

1.2.5 "網(wǎng)格無關(guān)性驗證

本文選取網(wǎng)格數(shù)量為422 044、706 901、1 349 714、1 709 204、2 021 597、2 591 680的6種情況進(jìn)行無關(guān)性驗證，取兩流體出口溫度為參考因素。對比發(fā)現(xiàn)，1 349 714個網(wǎng)格與1 709 204個網(wǎng)格的結(jié)果僅相差2.56%，可采取第三種方案。為了節(jié)約計算資源，最終采用網(wǎng)格數(shù)量為135 W。

2 "實驗結(jié)果分析

2.1 "壓力分布模擬結(jié)果分析

從圖4得到不同橫板長度的L型折流板換熱器壓力場分布相似。殼程流體入口壓力最高，沿流道逐漸降低。流體沖刷管束時，向流處壓力較大，背流處壓力較小。流體沖刷折流板導(dǎo)致折流板附近壓力較高。

通過對比，折流板橫板長度過長，橫板與豎板間的流道過窄，導(dǎo)致此處流體壓力過大，對管束與折流板沖刷嚴(yán)重，不利于結(jié)構(gòu)穩(wěn)定，會降低壽命，嚴(yán)重時甚至流體無法順利流出，導(dǎo)致殼體破裂，產(chǎn)生極大的安全隱患。折流板橫板長度過小，折流板橫板對流體流場影響不大，壓力分布與單弓形折流板相似。

通過Post導(dǎo)出數(shù)據(jù)可以發(fā)現(xiàn)，殼程壓降隨著橫板長度變成先降低后升高。其中，在80～120 mm附近殼程壓降最低，三種流速下80～120 mm橫板長度的折流板均比單弓形折流板壓降低了11%～13%。由此可以看出，L型折流板承受能力較傳統(tǒng)單弓形折流板更好。由于本文初始的L型折流板結(jié)構(gòu)就較好，故80 mm、100 mm與120 mm橫板長度L型折流板壓降差別不大。橫板長度較小時，由于L型折流板橫板能很好的引流減小流體對管束、折流板的橫向沖刷，故隨著橫板長度增加，壓降逐漸減小。當(dāng)橫板長度較長時，由于折流板橫板與豎板缺口較小，流體不易通過，故殼程壓降升高，若缺口過小，嚴(yán)重時會導(dǎo)致?lián)Q熱器殼程破損。

2.2 "溫度分布模擬結(jié)果分析

從圖8可以看出，不同橫板長度的L型折流板換熱器溫度場分布相似。殼程流體入口流體最低，沿流道逐漸升高。流體沖刷管束時，向流處溫度較低，背流處出現(xiàn)滯留區(qū)，流體溫度較高。兩折流板中間位置由于流體沖刷嚴(yán)重，換熱明顯，故溫度較高，折流板附近流體與冷流體混合，溫度較低。

通過對比，折流板橫板長度過長，橫板與豎板間的流道過窄，導(dǎo)致此處流體滯留，溫度較高，對管束與折流板沖刷劇烈，換熱更為明顯，但不利于結(jié)構(gòu)穩(wěn)定，會產(chǎn)生較大的熱應(yīng)力，從而降低使用壽命。折流板橫板長度過小，折流板橫板對流體流場影響不大，溫度分布與單弓形折流板相似。

通過對比三種流速下的殼程出口溫度隨橫板長度變化曲線可以發(fā)現(xiàn)，當(dāng)橫板長度過長時，殼程流體出口溫度會升高，這是因為橫板長度過長時，橫板與豎板間缺口較小，隨著流速增加，缺口處壓強(qiáng)增大，流體也不能很好的流出，故此處流體停留時間較長，流體溫度升高。由于折流板橫板過厚與管束交匯造成的熱損失和折流板橫板長度增加對換熱效率提高會有一個綜合的最佳值，約在60～80 mm附近。當(dāng)流速為1 m/s時，80 mm的L型折流板比單弓形折流板殼程流體出口溫度提高了7.29%，在不同流速下，殼程流體出口攝氏溫度較單弓形折流板最多降低5.18%。

2.3 "速度分布模擬結(jié)果分析

折流板橫板能使得流體在其阻擋下進(jìn)行回轉(zhuǎn)流動，流體在橫掠管束后產(chǎn)生回轉(zhuǎn)流動。雖然殼程出口流速變化較為明顯，但由Post導(dǎo)出出口流量與入口流量持平，確保了質(zhì)量守恒。在橫板與另一塊折流板豎板間夾縫處，由于流道尺寸發(fā)驟變，該處流體流速較高，以橫板長度為200 mm最為明顯，這不利于結(jié)構(gòu)穩(wěn)定。在橫板較短的折流板背部易形成滯留區(qū)，不利于換熱。故適中的折流板橫板長度有利于穩(wěn)定換熱器結(jié)構(gòu)，提升換熱性能。

2.4 "綜合性能分析

考慮換熱器的溫度變化與壓力損失，本文根據(jù)Z因子[16]進(jìn)行換熱器綜合性能評價。Z因子評價標(biāo)準(zhǔn)如下：

（9）

式中：Δt，Δt0—優(yōu)化前后換熱器殼程進(jìn)出口溫差，K。

ΔP/ΔP0—優(yōu)化前后換熱器殼程壓力損失，Pa。

由于邊界條件設(shè)定各種情況下流體入口溫度與出口壓力相同，故可以對Z因子進(jìn)行簡化：

（10）

式中：T0—殼程出口溫度，K；

Pi—殼程入口壓力，Pa。

通過Z因子評價換熱器綜合性能，得殼程入口流速分別為1 m/s、3 m/s和5 m/s時，80 mm橫板長度的L型折流板比單弓形折流板換熱器綜合性能分別提高了12.3%、12.4%和13.7%，比其他橫板長度都高。故選擇80 mm橫板長度為最佳長度。

4 "結(jié) 論

本文對1 m/s、3 m/s和5 m/s殼程入口流速下的橫板長度為20 mm、40 mm、80 mm、120 mm、160 mm和200 mm的L型折流板換熱器的壓力場、溫度場、流場進(jìn)行分析，并與單弓形折流板換熱器內(nèi)部流場特性進(jìn)行對比，得出如下結(jié)論：

1）L型折流板橫板較短時，換熱器內(nèi)壓力、溫度和流場分布與單弓形折流板相似；橫板長度過長時，折流板橫板豎板間的缺口較小，流體通過困難，換熱器內(nèi)壓力損失較大，溫差較大，對管束與殼體的影響十分明顯，不利于結(jié)構(gòu)穩(wěn)定。

2）從壓力損失上看，在1 m/s、3 m/s和5 m/s三種流速下，L型折流板橫板長度為80 mm時，殼程壓力損失最小，比單弓形折流板低了11%～13%。

3）從換熱效率來看，流速為1 m/s時，80 mm時L型折流板比單弓形折流板殼程流體出口溫度提高了7.29%，在3 m/s和5 m/s兩種流速下，L型折流板殼程流體出口溫度較比弓形折流板的最大降幅為5.18%。

4）根據(jù)Z因子，L型折流板換熱器橫板長度為80 mm時綜合性能最高，比單弓形折流板換熱器提高了12.3%。

參考文獻(xiàn)：

［1］李虎. 低溫余熱發(fā)電系統(tǒng)的試驗分析與應(yīng)用研究[D]. 天津大學(xué)，2012.

［2］梁紅娟. 淺談我國的能源現(xiàn)狀及能源對策[J]. 甘肅科技，2019，35 （15）： 6-8.

［3］薛欽源，張祝愷. 后疫情時代中國能源發(fā)展形勢與應(yīng)對措施[J]. 中國經(jīng)貿(mào)導(dǎo)刊（中），2021 （02）： 85-87.

［4］楊世銘，陶文銓. 傳熱學(xué)[M]. 北京：高等教育出版社，1998.

［5］LI H D，KOTTKE V. Effect of the leakage on pressure drop and local heat transfer in shell-and-tube heat exchangers for staggered tube arrangement[J]. International Journal of Heat and Mass Transfer，1998，41 （2） ：425-433.

［6］ARTHUR E，BERGLES. ExHFT for fourth generation heat transfer technology[J]. Experimental Thermal and Fluid Science，2002，26 （2）：335-344.

［7］梁帥，黃雪琴，余海洋.不同螺旋角折流板換熱器流動與傳熱性能研究[J]. 化工設(shè)備與管道，2019，56 （05）： 32-38.

［8］邵海磊，劉彤，許偉峰，等. 折流板結(jié)構(gòu)對換熱器殼程性能影響的對比研究[J]. 化工設(shè)備與管道，2016，53 （01）： 24-27.

［9］SAEED M， FARZAD V. Numerical investigation of small plate heat exchangers performance having different surface profiles[J]. Applied Thermal Engineering，2021，188：116616.

[10] 喻九陽，王澤武，馮興奎. 單弓形折流板開孔試驗研究[J]. 石油化工設(shè)備，2004 （02）： 4-6.

[11] WEN J，YANG H Z， WANG S M，et al. Experimental investigation on performance comparison for shell-and-tube heat exchangers with different baffles[J]. International Journal of Heat and Mass Transfer，2015，84：990-997.

[12] 王英雙. 縱流管殼式換熱器流動與傳熱性能的理論與實驗研究[D].華中科技大學(xué)，2011.

[13] 江竹，秦健，張輝. 管殼式換熱器折流板對換熱性能影響的數(shù)值模擬[J]. 山東化工，2020，49 （02）：159-161.

[14] 田朋. 管殼式換熱器結(jié)構(gòu)設(shè)計研究[J]. 中國設(shè)備工程，2021（23）：2.

[15] 卜英南，陳祥彬，張瑩瑩，等. 折流板結(jié)構(gòu)對換熱器殼程流動和傳熱性能的影響[J]. 遼寧石油化工大學(xué)學(xué)報，2018，38 （05） ：72-76.

[16] 張少琴，楊維陽，姚河省，等. 關(guān)于復(fù)合材料板斷裂分析的z因子討論[J]. 太原重型機(jī)械學(xué)院學(xué)報，1999（03）：6.

Optimal Design of L-type Baffle in Heat Exchangers

WAN Pan1， YU Wen-zhu2， MA Ying2， ZHANG Tao2， ZHU Jin-yuan2

（1. Wuhan Maritime Communication Research Institute （WMCRI）， Wuhan Hubei 430205， China;

2. School of Energy and Environment， Shenyang Aerospace University， Shenyang Liaoning 110136， China）

Abstract： "The baffle plate in the heat exchanger is the key component affecting its heat exchange efficiency. The optimal design of the baffle plate can effectively improve the heat exchange efficiency of the heat exchanger， so as to improve the energy utilization and save energy. In this paper， the new L-shaped baffle was studied. Using the principle of CFD， the geometric model of L-shaped baffle heat exchanger was established， and it was compared with the general single bow baffle. The results showed that when the length of L-shaped baffle was 80 mm， its comprehensive performance was the best， the pressure loss was 11% ～13% lower than that of single bow baffle， and the heat exchange efficiency was increased by 7.29% when the flow rate was 1～3 m·s-1; When the flow rate was greater than 3 m·s-1， the maximum decrease was 5.18%， and when the flow rate was 1 m·s-1， 3 m·s-1 and 5 m·s-1， the comprehensive performance of the heat exchanger was improved by 12.3%， 12.4% and 13.7% respectively， which could provide some reference for the industrial design of the heat exchanger.

Key words： Heat exchanger; Baffle plate; Optimization design; Heat exchange efficiency; Comprehensive performance