孔板結(jié)構(gòu)換熱器傳熱與阻力性能的數(shù)值模擬

侯夏玲，周幗彥，涂善東

(華東理工大學(xué)機(jī)械與動(dòng)力工程學(xué)院承壓系統(tǒng)與安全教育部重點(diǎn)實(shí)驗(yàn)室，上海 200237)

孔板結(jié)構(gòu)換熱器傳熱與阻力性能的數(shù)值模擬

侯夏玲，周幗彥，涂善東

(華東理工大學(xué)機(jī)械與動(dòng)力工程學(xué)院承壓系統(tǒng)與安全教育部重點(diǎn)實(shí)驗(yàn)室，上海 200237)

基于周期性全截面模型及RNGk-ε湍流模型，運(yùn)用計(jì)算流體力學(xué)軟件FLUENT對不同孔板結(jié)構(gòu)換熱器殼程流體流動(dòng)以及傳熱性能進(jìn)行了數(shù)值模擬分析，并通過文獻(xiàn)試驗(yàn)數(shù)據(jù)驗(yàn)證了該數(shù)值模擬方法的可行性和準(zhǔn)確性。在此基礎(chǔ)上，對比分析了三葉孔、四葉孔、五葉孔、大圓孔、小圓孔等5種孔板結(jié)構(gòu)的傳熱與阻力性能，探討了支撐板等結(jié)構(gòu)參數(shù)對其傳熱與阻力性能的影響，進(jìn)一步采用場協(xié)同原理探討了孔板換熱器的強(qiáng)化傳熱機(jī)理。研究結(jié)果表明：采用RNGk-ε湍流模型以及周期性全截面模型可較為準(zhǔn)確地模擬孔板換熱器殼程流體流動(dòng)情況；5種模型中五葉孔換熱器的傳熱特性最好但阻力最大，小圓孔的傳熱效果最差但阻力最??；隨著支撐板間距以及開孔高度的增加，換熱器殼程的傳熱系數(shù)和壓力降均逐漸降低；在支撐板后，速度矢量與溫度梯度之間的夾角波動(dòng)幅度變化劇烈，起到了強(qiáng)化殼程傳熱的效果；其中五葉孔板的場協(xié)同角波動(dòng)幅度最大，強(qiáng)化傳熱效果最好。

孔板結(jié)構(gòu)；周期性；傳熱；計(jì)算流體力學(xué)；數(shù)值模擬

引 言

管殼式換熱器是應(yīng)用較為廣泛的一類換熱器，其中折流板換熱器具有結(jié)構(gòu)簡單且易于制造等特點(diǎn)，但同時(shí)有流動(dòng)阻力大、存在流動(dòng)死區(qū)、易結(jié)垢等缺點(diǎn)[1-3]。基于孔板的縱流殼程換熱器采用與傳統(tǒng)折流板換熱器不同的殼程結(jié)構(gòu)，使殼程流體由錯(cuò)流變?yōu)榭v向流動(dòng)，因而具有傳熱系數(shù)高、壓降小、流動(dòng)“死區(qū)”少、質(zhì)量輕、抗振動(dòng)、節(jié)省材料等諸多優(yōu)點(diǎn)，是目前被廣泛研究與應(yīng)用的一種新型管殼式換熱器[4-6]。

目前對于孔板式換熱器的研究主要包括孔板結(jié)構(gòu)和數(shù)值模擬兩個(gè)方面。最早的孔板式縱流殼程換熱器采用大圓孔支撐板[7]，殼程流體從圓孔中通過，呈縱向流動(dòng)，同時(shí)降低殼程壓降。隨后開發(fā)了小圓孔、梅花孔、三葉孔等其他形式的異形孔折流板[8-10]?？装迨娇v流殼程換熱器在支撐板處的流通面積較小，由于流通面積較小產(chǎn)生射流，提高了局部流體湍流程度，從而實(shí)現(xiàn)殼程強(qiáng)化傳熱[11-12]。El Maakoul等[13]的研究表明孔板式換熱器的傳熱效率優(yōu)于折流板以及螺旋折流板換熱器。1974年，Patankar等[14]首次提出了多孔介質(zhì)模型并將管殼式換熱器進(jìn)行簡化。Prithiviraj等[15-17]在多孔介質(zhì)模型基礎(chǔ)上結(jié)合分布阻力、體積多孔度以及表面滲透度等概念提出了一種三維CFD方法——HEATX。曾文良等[18]以多孔介質(zhì)、分布阻力模型為基礎(chǔ)，采用RNGk-ε湍流模型對大型換熱器殼程傳熱與流動(dòng)進(jìn)行了數(shù)值研究。運(yùn)用多孔介質(zhì)模型雖然可以減少數(shù)值模擬時(shí)間，但數(shù)值模擬過程所需要的分布阻力、分布熱源等多孔介質(zhì)參數(shù)需通過試驗(yàn)測試得到，這些都限制多孔介質(zhì)模型的應(yīng)用。董其伍等[19-20]提出了周期性“單元流道”模型和周期性全截面模型，這兩種模型都是對管殼式換熱器殼程數(shù)值模擬進(jìn)行簡化的計(jì)算方法。You等[21-22]運(yùn)用“單元流道”模型對殼程進(jìn)行模擬，并通過試驗(yàn)得到了殼程傳熱及準(zhǔn)則關(guān)系式。“單元流道”模型不能詳細(xì)準(zhǔn)確反映換熱器內(nèi)部流場，尤其是靠近進(jìn)出口和殼體壁面部分。周期性全截面模型可以真實(shí)地反映殼程流體流動(dòng)的形式、流場分布及整體的傳熱性能[23-27]。

本文對孔板式縱流殼程換熱器建立周期性全截面模型，基于商業(yè)CFD軟件FLUENT 14.0及RNGk-ε湍流模型，提出一種三葉孔板換熱器殼程流動(dòng)和傳熱數(shù)值模擬方法，并通過試驗(yàn)驗(yàn)證該模型和數(shù)值模擬方法的可行性及準(zhǔn)確性。在此基礎(chǔ)上，對比分析三葉孔板、四葉孔板、五葉孔板、大圓孔板、小圓孔板等 5種異形孔支撐板結(jié)構(gòu)的傳熱與阻力性能，進(jìn)一步探討支撐板結(jié)構(gòu)及幾何參數(shù)對其傳熱與阻力性能的影響，并采用場協(xié)同原理探討孔板換熱器的強(qiáng)化傳熱機(jī)理，為孔板換熱器設(shè)計(jì)優(yōu)化提供理論依據(jù)。

1 數(shù)值模擬計(jì)算模型

對于縱流殼程換熱器性能的數(shù)值計(jì)算，整體模型、周期性全截面模型、單元流道模型已經(jīng)被廣泛應(yīng)用。綜合考慮3類計(jì)算模型的優(yōu)缺點(diǎn)，采用周期性全截面模型對5種孔板結(jié)構(gòu)進(jìn)行數(shù)值計(jì)算。

本文建立三葉孔、四葉孔、五葉孔、大圓孔、小圓孔等5種縱流殼程換熱器的模型，運(yùn)用周期性全截面模型對5種不同支撐板換熱器進(jìn)行數(shù)值模擬研究。用周期性全截面模型進(jìn)行數(shù)值模擬時(shí)，做了以下簡化假設(shè)：

（1）換熱器殼體與外界絕熱；

（2）殼程流體為湍流，且是穩(wěn)態(tài)流動(dòng)；

（3）換熱管壁面溫度恒定；

（4）忽略殼程各個(gè)流道之間的相互影響。

1.1 物理模型

圖1 整體結(jié)構(gòu)Fig.1 Heat exchangers with perforated-baffle structure

孔板換熱器整體模型如圖1所示。為了減小進(jìn)出口對流動(dòng)和傳熱的影響，選取了包含2塊支撐板在內(nèi)的幾何模型?？紤]到殼程流道模型具有明顯的對稱性，選取殼程中間段的一半進(jìn)行數(shù)值模擬?？装迨綋Q熱器的殼程周期性全截面模型如圖2所示。

圖2 孔板換熱器殼程周期性全截面模型Fig.2 Simplified periodic model of shell-side in perforated-baffle heat exchangers

朱凌云等[28]的數(shù)值模擬研究表明：六邊形結(jié)構(gòu)的導(dǎo)流筒換熱器殼程流體流動(dòng)“死區(qū)”比圓形導(dǎo)流筒換熱器小，流體流動(dòng)分布更均勻，因此本文均采用六邊形導(dǎo)流筒且尺寸相同。5種異型孔板換熱器的支撐板的間距、換熱管間距都相同，分別為 400 mm、19 mm。四葉孔由于結(jié)構(gòu)特殊采用正方形布管，換熱管數(shù)為29，其他異形孔板換熱器采用三角形布管，換熱管數(shù)為30。5種異型孔支撐板的開孔高度H一致，H的值為換熱管半徑與開孔半徑的差值，即H=|R-r|，5種異型孔支撐板的開孔結(jié)構(gòu)如圖3所示。

圖3 孔板換熱器開孔結(jié)構(gòu)示意圖Fig.3 Pore structure of heat exchanger with perforated-baffle

本文對不同支撐板結(jié)構(gòu)的5種周期性全截面模型進(jìn)行數(shù)值模擬分析。具體結(jié)構(gòu)參數(shù)及幾何尺寸見表1。

表1 孔板換熱器幾何尺寸Table 1 Geometric parameters of heat exchanger with perforated-baffle

1.2 控制方程和數(shù)值計(jì)算方法

本文采用RNGk-ε湍流模型對三葉孔板換熱器殼程流體流動(dòng)及傳熱進(jìn)行模擬計(jì)算，換熱器殼程內(nèi)流體為單相連續(xù)流體，在不考慮體積力的情況下，應(yīng)滿足如下通用控制方程

式中，Γφ和Sφ分別為廣義擴(kuò)散系數(shù)和廣義源項(xiàng)。當(dāng)φ取1及U、T、k和ε等不同變量時(shí)，方程分別代表了連續(xù)性方程、動(dòng)量方程、能量方程、湍動(dòng)能k方程和湍流耗散率ε方程。

進(jìn)出口采用周期性邊界條件，選取水作為殼程流體介質(zhì)，入口溫度為 283 K，其物性參數(shù)由piecewise-linear分段線性插值給出。殼體壁面、支撐板面均設(shè)為絕熱面，換熱管壁面采用恒壁溫邊界條件，溫度為307 K。

數(shù)值模擬計(jì)算采用計(jì)算流體力學(xué)軟件FLUENT，求解質(zhì)量和動(dòng)量方程采用有限體積法，壓力與速度的耦合采用SIMPLE算法，壓力項(xiàng)采用標(biāo)準(zhǔn)格式離散，其余項(xiàng)采用二階迎風(fēng)格式。

1.3 網(wǎng)格劃分及無關(guān)性驗(yàn)證

采用商業(yè)軟件ANSYS中的Meshing對換熱器殼程整體模型進(jìn)行網(wǎng)格劃分，為了保證網(wǎng)格質(zhì)量，提高計(jì)算精度，對換熱器模型網(wǎng)格進(jìn)行了分塊劃分，在支撐板附近采用四面體網(wǎng)格，并進(jìn)行加密；在支撐板以外的區(qū)域，由于結(jié)構(gòu)比較規(guī)則，采用掃略的方法生成六面體網(wǎng)格。具體網(wǎng)格劃分情況如圖4所示。

為保證計(jì)算的準(zhǔn)確性，在數(shù)值計(jì)算前，同樣進(jìn)行了網(wǎng)格無關(guān)性分析。以三葉孔板的模型為例，采用 5種不同數(shù)量（3008753、3879544、4474431、5229371、6003791）的網(wǎng)格進(jìn)行了數(shù)值模擬，計(jì)算出來的殼程傳熱系數(shù)h和壓降Δp如表2所示。

當(dāng)殼程Reynolds數(shù)相同時(shí)，后3套網(wǎng)格計(jì)算出的傳熱系數(shù)h和壓降Δp分別為26076.3 W·(m2?K)-1、70.04 kPa，26086.5 W·(m2?K)-1、69.53 kPa，26096.7 W·(m2?K)-1、68.76 kPa，結(jié)果相差均在 1.5%以內(nèi)。綜合考慮計(jì)算精度和計(jì)算效率，選用第3套單元數(shù)為4.47×106網(wǎng)格模型進(jìn)行計(jì)算分析。

表2 不同網(wǎng)格數(shù)數(shù)值計(jì)算結(jié)果比較Table 2 Comparison of result calculated by different grid number

1.4 模型驗(yàn)證

目前使用最為廣泛的3類湍流模型：標(biāo)準(zhǔn)k-ε、RNGk-ε和realizablek-ε模型，都適用于湍流核心區(qū)的流動(dòng)，但這3類模型的使用情況并沒有明確的研究。對于本文研究的縱流殼程換熱器數(shù)值模擬所使用的湍流模型也沒有統(tǒng)一的標(biāo)準(zhǔn)。為了驗(yàn)證湍流模型的準(zhǔn)確性，建立周期性全截面模型，與文獻(xiàn)[29]試驗(yàn)數(shù)據(jù)進(jìn)行對比。試驗(yàn)用換熱器的結(jié)構(gòu)尺寸與數(shù)值模擬相同，試驗(yàn)過程中采用水-水換熱，其試驗(yàn)系統(tǒng)和測試工程與文獻(xiàn)[29]相同，測試結(jié)果經(jīng)后期處理得到換熱器殼程傳熱和阻力特性的試驗(yàn)結(jié)果。將數(shù)值計(jì)算結(jié)果與試驗(yàn)測試結(jié)果進(jìn)行對比，結(jié)果如圖5所示。

由圖5可知，采用不同湍流模型的殼程N(yùn)u和Δp變化趨勢與試驗(yàn)測試結(jié)果基本一致。由于試驗(yàn)測試系統(tǒng)本身存在一定誤差，為方便計(jì)算對物理模型也進(jìn)行了適當(dāng)簡化，使得數(shù)值計(jì)算結(jié)果與試驗(yàn)測試結(jié)果存在一定的誤差。由圖5的對比結(jié)果可知，RNGk-ε的計(jì)算結(jié)果與試驗(yàn)結(jié)果最為接近，其中殼程N(yùn)u最大誤差為17.3%，Δp最大誤差為13.42%。因此，本文采用RNGk-ε湍流模型以及周期性全截面模型是準(zhǔn)確的。

圖5 不同湍流模型Nusselt數(shù)及壓降與試驗(yàn)結(jié)果對比Fig.5 Comparison of Nusselt number and pressure drop of different turbulence model with test results

2 殼程流體流動(dòng)性能對比分析

流體在異形孔板換熱器殼程都是作縱向流動(dòng)的，但由于支撐板上開孔形狀的差異，流體在流過支撐板后，形成的射流流動(dòng)不盡相同，由此對殼程傳熱的影響也不相同。分別截取支撐板間及其前后間隔20 mm距離的5個(gè)截面，對比分析5種異形孔板換熱器殼程內(nèi)流體的速度、溫度及壓力分布情況，如圖6所示。

圖6(a)為異形孔板換熱器5個(gè)截面的速度分布。從圖中可以看出：流體流動(dòng)在殼程呈現(xiàn)縱向流動(dòng)，流體流過支撐板時(shí)，由于流道變小導(dǎo)致流速增大，從而強(qiáng)化傳熱。支撐板前的流場分布較為均勻，在支撐板截面處，流體流過支撐板時(shí)流道減小，流速增大；由于小孔的射流作用，殼程局部流速顯著增大，且對支撐板后流場有較強(qiáng)的擾流作用。

圖6(b)為異形孔板換熱器 5個(gè)截面的溫度分布。從圖中可以看出：沿著流體流動(dòng)方向，流體溫度逐漸升高，在支撐板的強(qiáng)化后，溫度上升得較為明顯。其中五葉孔的溫度較高且分布較為均勻，五葉孔傳熱效果較好；四葉孔因特殊的結(jié)構(gòu)采取正方形布管，殼程的流通面積較大，溫度分布不均，傳熱效果較差。

圖6(c)為異形孔板換熱器5個(gè)截面的壓力分布。從圖中可以看出：沿著流體的流動(dòng)方向，壓降主要發(fā)生在支撐板處，這是由于流體流過支撐板時(shí)，流道變小、流速增大導(dǎo)致壓力減小，支撐板前后的壓力降增大；對比分析可知，流體經(jīng)過五葉孔支撐板比小圓孔支撐板的壓力降更大，而三葉孔的壓力降位于兩者之間。

圖6 5種換熱器截面流場分布情況Fig.6 Distribution of flow field for five heat exchangers cross section

3 殼程流體傳熱與阻力性能對比分析

5種異形孔板換熱器殼程傳熱及壓降的計(jì)算結(jié)果如圖7所示。從圖中可以看出，5種換熱器傳熱系數(shù)和壓降都隨Re增加而增大，其中五葉孔支撐板換熱器傳熱效果最好。在Re增大的過程中，三葉孔板換熱器殼程傳熱系數(shù)幾乎與四葉孔的相同。小圓孔的傳熱效果最差，這是由于小圓孔支撐板的開孔位置位于流道中間，流速大的流體對管壁的影響不大，而其余異形孔板的開孔位置都位于換熱管壁周圍，流體直接沖刷管壁，更有利于殼程的換熱。從圖7中還可以看出，由于五葉孔支撐板中流道較小，對流體阻礙作用較大，因此五葉孔支撐板換熱器的殼程壓力梯度最大，阻力性能最差。

從傳熱來看，在管壁周圍開孔的孔板結(jié)構(gòu)破壞了換熱管壁面的邊界層，更有利于流體換熱，其中五葉孔支撐板換熱器傳熱效果最好。小圓孔支撐板的開孔位置位于流道的中部，高速流體不對換熱管壁進(jìn)行直接沖刷，使得小圓孔的傳熱效果最差；從壓降來看，小圓孔的阻力最小，由于五葉孔支撐板的開孔面積小，使得五葉孔的阻力較大。

4 綜合性能對比分析

為了將傳熱與阻力性能進(jìn)行綜合比較，引入綜合性能評價(jià)因子——η=h/(Δp/L)，將傳熱性能和阻力性能結(jié)合在同一指標(biāo)內(nèi)綜合考慮。圖8給出了5種不同換熱器殼程綜合性能對比結(jié)果。從圖中可以看出：5種換熱器綜合性能都隨Re增加而減小，且彼此之間的差距逐漸減?。籖e＜15000時(shí)，三葉孔板換熱器綜合性能優(yōu)于其他孔板換熱器；而Re＞15000時(shí)，三葉孔板、大圓孔、小圓孔板換熱器綜合性能較為接近；由于五葉孔和四葉孔支撐板換熱器的壓降較大，使得五葉孔和四葉孔支撐板換熱器綜合性能較差。因此，在Re＜15000時(shí)優(yōu)先選用三葉孔板換熱器，而Re＞15000時(shí)，可優(yōu)先選擇三葉孔板或大圓孔板或小圓孔板換熱器。

圖7 5種孔板換熱器殼程傳熱與阻力性能Fig.7 Heat transfer coefficient and pressure drop on shell-side of five perforated-baffle heat exchangers

圖8 5種孔板換熱器綜合性能隨Reynolds數(shù)的變化關(guān)系Fig.8 Change of comprehensive properties with Reynolds number of five perforated-baffle heat exchangers

5 結(jié)構(gòu)參數(shù)敏感性分析

5.1 支撐板間距的影響

建立不同板間距的幾何模型，在Re=30143時(shí)，支撐板間距對換熱器殼程傳熱系數(shù)和壓力降的變化關(guān)系如圖9所示。從圖中可以看出，在Re相同時(shí)，換熱器殼程的傳熱系數(shù)和壓力降都隨著支撐板間距的增加而逐漸降低，且板間距的變化對傳熱系數(shù)和壓力降的影響逐漸減小。支撐板間距相同時(shí)，五葉孔支撐板的傳熱效果最好但壓降最大，小圓孔的傳熱效果最差，大圓孔的壓降最小。

圖9 支撐板間距對傳熱及壓降的影響Fig.9 Effects of baffle spacing on performance of heat transfer and pressure drop

5.2 開孔高度的影響

建立不同開孔高度的幾何模型，在Re=30143時(shí)，開孔高度對換熱器殼程傳熱系數(shù)和壓力降的影響如圖10所示。從圖中可以看出，隨著開孔高度的增大，換熱器殼程的傳熱系數(shù)和壓力降都逐漸減小。這是由于在Re不變的情況下，較小的開孔高度有利于提高流體流過開孔時(shí)的射流流速，可以更好地沖刷換熱管壁，減薄流動(dòng)邊界層，從而加強(qiáng)了局部區(qū)域內(nèi)流體的換熱。然而，較小的開孔高度對流體流動(dòng)阻力也相應(yīng)增大，流體流過支撐板時(shí)的局部動(dòng)量損失也增大。開孔高度相同時(shí)，五葉孔支撐板的傳熱效果最好但壓降最大，小圓孔的傳熱效果最差，大圓孔的壓降最小。

圖10 開孔高度對傳熱及壓降的影響Fig.10 Influence of hole-height on performance of heat transfer and pressure drop

6 強(qiáng)化傳熱機(jī)理分析

過增元[30]提出了對流傳熱性能不僅與溫度場和速度場分布有關(guān)，還與速度矢量與溫度梯度之間的夾角有關(guān)，這個(gè)角叫協(xié)同角。將能量方程寫出點(diǎn)積形式

其中θ是速度矢量和溫度梯度的夾角，即協(xié)同角。

從式(2)可以看出，速度場與溫度場的協(xié)同作用能夠強(qiáng)化傳熱，在相同的速度和溫度邊界條件下，它們的協(xié)同程度越好，則換熱強(qiáng)度就越高，該理論即為場協(xié)同理論。

為了分析不同孔板結(jié)構(gòu)換熱器殼程流體流動(dòng)及傳熱場協(xié)同性，本文對換熱器殼程的速度場及溫度梯度場的協(xié)同性進(jìn)行計(jì)算，通過編譯 UDF程序計(jì)算了兩者之間的夾角。中間段協(xié)同角與溫度之間的耦合關(guān)系見圖11。

圖11 協(xié)同角沿縱向位置變化關(guān)系Fig.11 Variation of synergy angle along with z position

從圖中可以看出：在沒有支撐板段，速度矢量與溫度梯度的夾角在某個(gè)范圍內(nèi)變化幅度很??；在支撐板前，協(xié)同角有一定幅度的變化；在支撐板尾部，協(xié)同角波動(dòng)幅度變大，后逐漸趨于穩(wěn)定。其中五葉孔的協(xié)同角波動(dòng)幅度最為劇烈，小圓孔的協(xié)同角波動(dòng)幅度較小。由于不同支撐孔板的特殊結(jié)構(gòu)，協(xié)同角在支撐板后的變化幅度不同，因此不同程度地強(qiáng)化了殼程傳熱。

7 結(jié) 論

本文建立了包括三葉孔、四葉孔、五葉孔、大圓孔、小圓孔等5種異形孔支撐板換熱器殼程周期性全截面模型，并采用計(jì)算流體力學(xué)軟件對其傳熱與阻力性能進(jìn)行了數(shù)值模擬，對比分析了其殼程流場和結(jié)構(gòu)參數(shù)對其熱力性能的影響關(guān)系，并進(jìn)一步探討了孔板支撐結(jié)構(gòu)的殼程強(qiáng)化傳熱機(jī)理。主要得出以下結(jié)論。

（1）建立了周期性全截面模型，對比分析標(biāo)準(zhǔn)k-ε、RNGk-ε和 realizablek-ε模型數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果，驗(yàn)證了RNGk-ε湍流模型以及周期性全截面模型的可行性和準(zhǔn)確性。

（2）在管壁周圍開孔的孔板結(jié)構(gòu)破壞了換熱管壁面的邊界層，更有利于流體換熱，其中五葉孔支撐板換熱器傳熱效果最好，小圓孔的傳熱效果最差；由于五葉孔支撐板的開孔面積較小，使得五葉孔的流動(dòng)阻力較大。

（3）不同開孔結(jié)構(gòu)的換熱器綜合性能都隨Reynolds數(shù)增加而減小，且彼此之間的差距也逐漸降低；Re＜15000時(shí)，三葉孔板換熱器綜合性能優(yōu)于其他孔板換熱器；Re＞15000時(shí)，三葉孔板換熱器綜合性能與大圓孔、小圓孔板換熱器較為接近。五葉孔和四葉孔支撐板換熱器綜合性能較差。

（4）在Re=30143下，殼程傳熱系數(shù)和壓力梯度都隨著支撐板間距的增大而逐漸減小，且支撐板間距對壓降的影響大于其對傳熱的影響。

（5）在Re=30143下，隨著開孔高度的增大，換熱器殼程的傳熱系數(shù)和壓力降都逐漸減小。較小的開孔高度有利于提高流體流過開孔時(shí)的射流流速，可以更好地沖刷換熱管壁，減薄流動(dòng)邊界層，從而強(qiáng)化了局部區(qū)域內(nèi)流體的換熱。

（6）由于不同支撐孔板的特殊結(jié)構(gòu)，改變了支撐板后速度矢量與溫度梯度之間的夾角，使得協(xié)同角在支撐板后波動(dòng)幅度較大，對殼程流體起到強(qiáng)化傳熱的效果。

[1]方書起,祝春進(jìn),吳勇,等.強(qiáng)化傳熱技術(shù)與新型高效換熱器研究進(jìn)展[J].化工機(jī)械,2004,31(4)：249-253.FANG S Q,ZHU C J,WU Y,et al.Research progress of strengthened heat transfer technique and new type high efficient heat exchangers[J].Chemical Engineering & Machinery,2004,31(4)：249-253.

[2]錢才富,高宏宇,孫海陽.曲面弓形折流板換熱器殼程流體流動(dòng)與傳熱[J].化工學(xué)報(bào),2011,62(5)：1233-1238.QIAN C F,GAO H Y,SUN H Y.Shell-side fluid flow and heat transfer in curved baffle heat exchanger[J].CIESC Journal,2011,62(5)：1233-1238.

[3]杜文靜,王紅福,程林.螺旋折流板換熱器殼程流動(dòng)特性[J].化工學(xué)報(bào),2014,65(8)：2970-2975.DU W J,WANG H F,CHENG L.Shell side fluid flow characteristics in helical baffled heat exchanger[J].CIESC Journal,2014,65(8)：2970-2975.

[4]董其伍,杜慶飛,劉敏珊,等.三葉孔板換熱器殼程流體流動(dòng)和傳熱特性數(shù)值研究[J].化工設(shè)備與管道,2012,49(2)：21-23.DONG Q W,DU Q F,LIU M S,et al.Numerical research on flow and heat transfer characteristics in shell-side of heat exchanger with trefoil-tube-support[J].Process Equipment & Piping,2012,49(2)：21-23.

[5]BUTTERWORTH D.A model for heat transfer during three-dimensional flow in tube bundles[C]// 6th International Heat Transfer Conference.Toronto,1978：219-224.

[6]姚璃,韓東,羅朝嘉,等.梅花形孔板換熱器殼程流動(dòng)與換熱特性模擬[J].能源化工,2017,38(1)：85-90.YAO L,HAN D,LUO C J,et al.Numerical simulation on shell-side flow and heat transfer of cinquefoil orifice baffle heat exchanger[J].Energy Chemical Industry,2017,38(1)：85-90.

[7]GREEN S J.Thermal hydraulic and corrosion aspects of PWR steam generator problems[J].Heat Transfer Engineering,1988,9(1)：19-22.

[8]WANG Y S,LIU Z C,HUANG S Y,et al.Experimental investigation of shell-and-tube heat exchanger with a new type of baffles[J].Heat& Mass Transfer,2011,47(7)：833-839.

[9]李雅俠,婁巖,戰(zhàn)洪仁,等.梅花孔板縱向流換熱器殼程流動(dòng)與傳熱的三維數(shù)值模擬[J].過程工程學(xué)報(bào),2015,15(4)：639-645.LI Y X,LOU Y,ZHAN H R,et al..Three-dimensional numerical simulation of fluid flow and heat transfer in the shell side of quincunx orifice baffle longitudinal flow type heat exchanger[J].The Chinese Journal of Process Engineering,2015,15(4)：639-645.

[10]羅朝嘉,韓東,牟春燕,等.孔板結(jié)構(gòu)對梅花形孔板換熱器殼程傳熱影響的模擬研究[J].能源化工,2016,37(1)：55-59.LUO C J,HAN D,MOU C Y,et al.Numerical simulation study on the effects of orifice structure on the shell-side heat transfer of cinquefoil orifice baffle heat exchanger[J].Energy Chemical Industry,2016,37(1)：55-59.

[11]朱凌云,周幗彥,朱冬生.縱流殼程換熱器傳熱性能研究進(jìn)展[J].化工機(jī)械,2014,41(3)：269-272.ZHU L Y,ZHOU G Y,ZHU D S.Progress in investigating heat transfer property of heat exchangers with longitudinal flow in shell side[J].Chemical Engineering & Machinery,2014,41(3)：269-272.

[12]孫海陽.大小孔折流板換熱器殼程場協(xié)同分析[J].機(jī)械工程與自動(dòng)化,2014,(4)：40-41.SUN H Y.Field synergy analysis LASH baffle heat exchanger in shell-side[J]. Mechanical Engineering & Automation,2014,(4)：40-41.

[13]EL MAAKOUL A,LAKNIZI A,SAADEDDINE S,et al.Numerical comparison of shell-side performance for shell and tube heat exchangers with trefoil-hole,helical and segmental baffles[J].Applied Thermal Engineering,2016,109：175-185.

[14]PATANKAR S V,SPALDING D B.A calculation procedure for the transient and steady state behavior of shell-and-tube heat exchanger[M]// AFGAN A A,SCHLUNDER E U.Heat Exchanger：Design and Theory Source Book.Washington：McGraw-Hill,1974.

[15]PRITHIVIRAJ M,ANDREWS M J.Three-dimensional numerical simulation of shell-and-tube heat exchangers(Ⅰ)：Foundation and fluid mechanics[J].Numerical Heat Transfer,Part A：Applications,1998,33(8)：799-816.

[16]PRITHIVIRAJ M,ANDREWS M J.Three-dimensional numerical simulation of shell-and-tube heat exchangers(Ⅱ)：Heat transfer[J].Numerical Heat Transfer,Part A：Applications,1988,33(8)：817-828.

[17]PRITHIVIRAJ M,ANDREWS M J.Comparison of a three-dimensional numerical model with existing methods for prediction of flow in shell-and-tube heat exchangers[J].Heat Transfer Engineering,1999,20(2)：15-19.

[18]曾文良,張復(fù)興,王劍秋,等.軸流管殼式換熱器殼側(cè)流體優(yōu)化分布數(shù)值研究[J].衡陽師范學(xué)院學(xué)報(bào),2009,30(3)：65-69.ZENG W L,ZHANG F X,WANG J Q,et al.Numerical investigation of fluid maldistribution of shell-side of shell-and-tube heat exchanger with longitudinal flow[J].Journal of Hengyang Normal University,2009,30(3)：65-69.

[19]董其伍,白彩鵬,劉敏珊,等.三葉孔板換熱器殼程流體流動(dòng)和傳熱研究[J].工程熱物理學(xué)報(bào),2013,34(1)：115-117.DONG Q W,BAI C P,LIU M S,et al.Research on flow and heat transfer in shell-side of the trefoil-baffle support type heat exchanger[J].Journal of Engineering Thermophysics,2013,34(1)：115-117.

[20]DONG Q W,WANG Y Q,LIU M S.Numerical and experimental investigation of shell side characteristics for ROD baffle heat exchanger[J].Applied Thermal Engineering,2008,28(7)：651-660.

[21]YOU Y H,CHEN Y Q,XIE M Q,et al.Numerical simulation and performance improvement for a small size shell-and-tube heat exchanger with trefoil-hole baffles[J].Applied Thermal Engineering,2015,89(5)：220-228.

[22]YOU Y H,FAN A W,LAI X J,et al.Experimental and numerical investigations of shell-side thermo-hydraulic performances for shell-and-tube heat exchanger with trefoil-hole baffles[J].Applied Thermal Engineering,2013,50(1)：950-956.

[23]周幗彥,史燕華,郭震,等.基于單元流道模型的三葉孔板換熱器熱力性能研究[J].壓力容器,2016,33(4)：10-17.ZHOU G Y,SHI Y H,GUO Z,et al.A study on the thermal performance of heat exchanger with trefoil-baffles by using unit duct model[J].Pressure Vessel Technology,2016,33(4)：10-17.

[24]WANG D,WANG K,WANG Y Q,et al.Numerical and experimental investigation of the shell side characteristics of the trefoil-hole baffle heat exchanger[J].Heat Transfer Research,2017,48(1)：81-95.

[25]魏志國,李華峰,柯漢兵,等.三葉孔板強(qiáng)化換熱性能及機(jī)理分析[J].化工進(jìn)展,2017,36(2)：465-472.WEI Z G,LI H F,KE H B,et al.Analysis on performance and mechanism of heat transfer enhancement of trefoil-hole baffle[J].Chemical Industry and Engineering Progress,2017,36(2)：465-472.

[26]董其伍,杜慶飛,劉敏珊,等.三葉孔板換熱器殼程流體流動(dòng)和傳熱特性數(shù)值研究[J].化工設(shè)備與管道,2012,49(2)：21-23.DONG Q W,DU Q F,LIU M S,et al.Numerical research on flow and heat transfer characteristics in shell-side of heat exchanger with trefoil-tube-support[J].Process Equipment & Piping,2012,49(2)：21-23.

[27]朱凌云,郎紅方,周幗彥,等.三葉孔板換熱器殼程流動(dòng)及傳熱數(shù)值模擬[J].化工學(xué)報(bào),2014,65(3)：829-835.ZHU L Y,LANG H F,ZHOU G Y,et al.Numerical simulation on shell side fluid flow and heat transfer in heat exchanger with trefoil-baffles[J].CIESC Journal,2014,65(3)：829-835.

[28]朱凌云,楊錦春,周幗彥,等.三葉孔板換熱器熱力性能及其影響因素分析[J].化工進(jìn)展,2014,33(12)：3183-3188.ZHU L Y,YANG J C,ZHOU G Y,et al.Investigation of the thermal performance of heat exchangers with trefoil-baffles[J].Chemical Industry and Engineering Progress,2014,33(12)：3183-3188.

[29]郭震.三葉孔板換熱器殼程傳熱與阻力性能的試驗(yàn)研究[D].上海：華東理工大學(xué),2014.GUO Z.Experimental investigation of the shell-side heat transfer and friction characteristics of shell-and-tube heat exchangers with trefoil-hole baffles[D].Shanghai：East China University of Science and Technology,2014.

[30]過增元.換熱器中的場協(xié)同原則及其應(yīng)用[J].機(jī)械工程學(xué)報(bào),2003,39(12)：1-9 GUO Z Y.Principle of field coordination in heat exchangers and its application[J].Chinese Journal of Mechanical Engineering,2003,39(12)：1-9.

date:2017-06-09.

ZHOU Guoyan,zhougy@ecust.edu.cn

supported by the National Natural Science Foundation of China(51675181) and Shanghai Pujiang Program(14PJD015).

Numerical study on heat transfer and resistance performance of heat exchangers with perforated-baffle structure

HOU Xialing,ZHOU Guoyan,TU Shandong
(Key Laboratory of Pressure Systems and Safety,Ministry of Education,School of Mechanical and Power Engineering,East China University of Science and Technology,Shanghai200237,China)

Based on simplified periodic model and RNGk-εturbulence model,the flow characteristics and heat transfer properties of shell-side fluid in perforated-baffle heat exchangers were numerically simulated by using CFD software FLUENT.The feasibility and accuracy of numerical simulation method was verified by experiment.The heat transfer and resistance performance of trefoil-baffle,four-leaf-hole baffle,five-leaf-hole baffle,big-hole baffle and small-hole baffle are analyzed.Then the influences of structural parameters on heat transfer and resistance performance were discussed.Furthermore,heat transfer enhancement mechanism of perforated-baffle structure heat exchangers was investigated based on the field synergy principle.The results showed that the fluid flow in the shell-side of the perforated-baffle structure heat exchanger can be accurately simulated by using the RNGk-εturbulence model and the simplified periodic model.The heat transfer coefficient of five-leaf-hole baffle is the best while the resistance the greatest.The heat transfer coefficient of the small-hole baffle is the worst,while the resistance is the least.The heat transfer coefficient and pressure drop decrease when baffle pitch and hole-height increase.Behind the support plate,the angle between the velocity vector and the temperature gradient varies violently,which enhances the heat transfer in the shell-side.The field synergy angle of five-leaf-hole baffle has the largest fluctuation range,and the enhancement of heat transfer effect is the best.

perforated-baffle structure; periodicity; heat transfer; computational fluid dynamics; numerical simulation

TK 124

A

0438—1157（2017）12—4517—09

10.11949/j.issn.0438-1157.20170744

2017-06-09收到初稿，2017-08-10收到修改稿。

聯(lián)系人：周幗彥。

侯夏玲（1992—），女，碩士研究生。

國家自然科學(xué)基金項(xiàng)目（51675181）；上海市浦江人才計(jì)劃（14PJD015）。