不同波紋夾角的人字形板式熱交換器數(shù)值模擬

吳 私 李春蘭 王 森 李占英 朱海舟

(1. 蘭州蘭石能源裝備工程研究院有限公司；2. 蘭州蘭石換熱設(shè)備有限責(zé)任公司)

不同波紋夾角的人字形板式熱交換器數(shù)值模擬

吳 私*1李春蘭1王 森1李占英1朱海舟2

(1. 蘭州蘭石能源裝備工程研究院有限公司；2. 蘭州蘭石換熱設(shè)備有限責(zé)任公司)

建立了由波紋夾角60、120°人字形板片組成的3種換熱流道的主換熱區(qū)模型，利用計(jì)算流體力學(xué)軟件對(duì)3種流道內(nèi)流體的流動(dòng)和傳熱進(jìn)行了數(shù)值模擬，分析了流道內(nèi)的速度場(chǎng)、溫度場(chǎng)和壓力場(chǎng)。結(jié)果表明：人字形波紋的夾角是影響板片間流體流動(dòng)形態(tài)的因素之一，流道的換熱性能和壓力損失均隨人字形波紋夾角的增大而增大，并且波紋夾角的改變對(duì)流道壓力損失的影響比換熱性能更明顯。

板式熱交換器 人字形板片 數(shù)值模擬 傳熱系數(shù) 壓力降

由于板式熱交換器具有結(jié)構(gòu)緊湊、傳熱性能好等優(yōu)點(diǎn)，已成為食品、化工、石油、冶金及電力等眾多領(lǐng)域不可或缺的重要設(shè)備[1,2]。隨著近年來數(shù)值模擬計(jì)算軟件的應(yīng)用，出現(xiàn)了許多圍繞換熱器流體流動(dòng)和傳熱的數(shù)值模擬研究。徐志明等對(duì)不同速度工況下熱交換器內(nèi)流體的流動(dòng)和傳熱進(jìn)行了數(shù)值模擬[3]；蔡毅等運(yùn)用數(shù)值模擬軟件模擬了人字形波紋區(qū)域流道內(nèi)流體的流動(dòng)與傳熱特性[4]；崔立祺等利用數(shù)值模擬計(jì)算結(jié)果擬合了人字形波紋板各幾何參數(shù)與努塞爾數(shù)、壓力降之間的關(guān)系曲線[5，6]；李彥洲通過研究板式熱交換器的內(nèi)部流道特征，采用均勻試驗(yàn)優(yōu)化方法，得到最優(yōu)波紋傾角、法相節(jié)距和波紋高度數(shù)值[7]；趙元東等通過數(shù)值模擬，分析了熱混合板式熱交換器的流道模型換熱和阻力的情況，并結(jié)合場(chǎng)協(xié)同理論分析熱混合板式換熱器的綜合性能[8]；陳文超等采用雙流道形式對(duì)人字形板式熱交換器進(jìn)行了溫度場(chǎng)的數(shù)值模擬[9]。

在人字形板式熱交換器的設(shè)計(jì)中，不能只追求較高的傳熱系數(shù)，同時(shí)還要考慮壓力降，尋找傳熱與阻力的合理匹配才是工作的重點(diǎn)。通過波紋夾角60°人字形波紋板(W)與夾角120°人字形波紋板(H)的不同的組合，可以形成3種不同性能的流道，即H-H流道、W-W流道、H-W流道。筆者借助Fluent軟件對(duì)這3種流道模型的主換熱區(qū)進(jìn)行流場(chǎng)模擬，分析了速度場(chǎng)、溫度場(chǎng)和壓力場(chǎng)，探索傳熱、壓降合理的最佳組合。

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

1.1物理模型

流道由兩種板片構(gòu)成，具體數(shù)值見表1。

表1 板片的特征參數(shù)

采用SolidWorks軟件對(duì)3種組合流道進(jìn)行了三維幾何建模，波紋沿流動(dòng)方向是周期性的，所以模型沿流動(dòng)方向截取長(zhǎng)100mm、寬50mm的區(qū)塊，幾何模型如圖1～3所示，圖1為由兩張H板調(diào)轉(zhuǎn)180°組成的換熱流道，圖2為由兩張W板調(diào)轉(zhuǎn)180°組成的換熱流道，圖3為由一張H板與另外一張W板調(diào)轉(zhuǎn)180°組成的換熱流道，圖3中上、下兩個(gè)換熱面分別和W板、H板相同。

圖1 H-H流道主換熱區(qū)周期性幾何模型

圖2 W-W流道主換熱區(qū)周期性幾何模型

圖3 H-W流道主換熱區(qū)周期性幾何模型

1.2相關(guān)假設(shè)

由于研究的傳熱過程沒有相變，流道內(nèi)的溫差較小，故進(jìn)行如下假設(shè)：流動(dòng)各物理量不隨時(shí)間變化，設(shè)為定常流動(dòng)；流體為不可壓縮的牛頓流體；重力和浮升力的影響忽略不計(jì)；忽略流體流動(dòng)時(shí)的粘性耗散作用所產(chǎn)生的熱效應(yīng)。

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

相關(guān)的控制方程如下：

(1)

(2)

(3)

湍流模型采用RNGk-ε模型

(4)

(5)

Gk——由平均速度梯度引起的湍動(dòng)能k的產(chǎn)生項(xiàng)；

p——壓力分量，Pa；

T——溫度，K；

Ui——i方向上的速度分量，m/s；

u、v、w——x、y、z方向上的速度分量，m/s；

α——流體熱擴(kuò)散率，m2/s；

αk、αε——湍動(dòng)能k和耗散率ε對(duì)應(yīng)的Prandtl，αk=αε=1.39；

μ——?jiǎng)恿φ扯?，Pa·s；

ρ——流體密度，kg/m3。

1.4網(wǎng)格劃分

筆者通過SolidWorks軟件構(gòu)建體模型，并通過ANSYS-Workbench軟件中的mesh工具進(jìn)行網(wǎng)格劃分。為提高計(jì)算精度、減少計(jì)算網(wǎng)格數(shù)，網(wǎng)格劃分采用六面體網(wǎng)格，網(wǎng)格尺寸為0.4mm，垂直板片方向設(shè)置20層網(wǎng)格，網(wǎng)格示意圖如圖4所示。

1.5邊界條件

進(jìn)出口采用周期性邊界條件，流道下側(cè)為入口，入口設(shè)置為速度入口，入口流速0.5m/s，溫度80℃；流道上側(cè)為出口，出口邊界采用壓力出口邊界，設(shè)定為0.1MPa。流道前后換熱面設(shè)置為恒熱流邊界條件，熱流密度為50kW/m2。兩側(cè)邊界設(shè)置為絕熱的對(duì)稱邊界條件。計(jì)算采用二階迎風(fēng)格式。因介質(zhì)流態(tài)為充分湍流，且存在周期性波動(dòng)，故為非穩(wěn)態(tài)計(jì)算。介質(zhì)為水，介質(zhì)物性如下：

密度ρ1 000kg/m3

比熱Cp4 182J/(kg·K)

導(dǎo)熱系數(shù)K0.6W/(m·K)

粘度μ1mPa·s

圖4 換熱板片主換熱區(qū)流道模型網(wǎng)格示意圖

2 數(shù)值模擬結(jié)果與分析

2.1板式熱交換器流道內(nèi)的速度場(chǎng)

人字形板式熱交換器的主要流動(dòng)形態(tài)有兩種[10]：十字交叉流和曲折流。流體先在一個(gè)板片上沿溝槽流動(dòng)，當(dāng)?shù)竭_(dá)板片的邊緣時(shí)，被反射折轉(zhuǎn)到通道另一個(gè)板片的溝槽中流向另一邊緣，這種形態(tài)稱為十字交叉流；流體仍沿溝槽流動(dòng)，但折返點(diǎn)不再出現(xiàn)在板片的左、右兩側(cè)，而是發(fā)生在波紋的觸點(diǎn)，流動(dòng)呈連續(xù)、并行的小波紋狀，這種形態(tài)稱為曲折流。

從圖5可以看出，流體在H-H流道內(nèi)流動(dòng)，當(dāng)遇到波紋觸點(diǎn)時(shí)，一部分流體改變?cè)瓉淼牧鲃?dòng)方向而流入另一側(cè)流道，然后螺旋向前流動(dòng)，形成所謂曲折流。此流道特點(diǎn)是湍流程度大，換熱效果強(qiáng)，但壓力損失也大。

圖5 H-H流道主換熱區(qū)換熱面流體流線圖

從圖6可以看出，W-W流道波紋夾角明顯小于H-H流道，流體在W-W流道內(nèi)流動(dòng)時(shí)，受波紋觸點(diǎn)的影響很小，流體始終沿自身流道流動(dòng)，形式主要是十字交叉流。此流道特點(diǎn)是流體動(dòng)量損失和壓力損失小，同時(shí)湍流程度較小，換熱效果差。

圖6 W-W流道主換熱區(qū)換熱面流體流線圖

從圖7可以看出，H-W流道的波紋夾角介于H-H流道、W-W流道之間，觸點(diǎn)附近流體的曲折流特征相對(duì)明顯，相比H-H流道，該流道內(nèi)壓力損失變小，傳熱效果有所減弱；而相比W-W流道，壓力損失變大，傳熱效果有所增強(qiáng)。

a. H換熱面

b. W換熱面

相互交叉的人字形波紋通道內(nèi)流體的流動(dòng)形態(tài)變化是由流體沿波紋方向的動(dòng)量分量和垂直波紋溝槽的流動(dòng)阻力兩方面的共同作用決定的。由于波紋夾角改變了流體的動(dòng)力分量方向，也會(huì)對(duì)垂直波紋溝槽的流動(dòng)阻力有影響，因此除了波紋傾角、波間距和波高，人字形波紋的夾角也會(huì)影響到板片間流體的流動(dòng)形態(tài)。

2.2板式熱交換器流道內(nèi)的溫度場(chǎng)和壓力場(chǎng)

流體在溝槽中流動(dòng)時(shí)在波紋的交匯點(diǎn)上會(huì)受到相對(duì)板片溝槽中流體流動(dòng)的影響，兩股流體交匯時(shí)的這種切向力會(huì)使流體產(chǎn)生漩渦并增大湍流程度，漩渦可以使流體對(duì)壁面換熱的參與度大幅增加從而能強(qiáng)化熱傳遞，同時(shí)觸點(diǎn)對(duì)周圍流體產(chǎn)生擾動(dòng)，形成湍流，這都是人字形板式熱交換器在低雷諾數(shù)時(shí)就能使傳熱得到有效強(qiáng)化的主要原因；產(chǎn)生的漩渦會(huì)提高流體的阻力，同時(shí)觸點(diǎn)周圍流體的相互拖曳作用也是波紋通道中流體流動(dòng)阻力較大的主要原因。

由于H-H流道波紋夾角大，兩股流體在溝槽中流動(dòng)交匯時(shí)的切向力大，同時(shí)該流道單位面積內(nèi)觸點(diǎn)更多，所以湍流程度更強(qiáng)，強(qiáng)化傳熱更充分，從圖8可以看出，H-H流道換熱面溫度均勻分布，傳熱效果好，同時(shí)流體的動(dòng)量損耗大；從圖9可以看出，H-H流道進(jìn)出口壓力損失明顯，呈現(xiàn)明顯的遞減趨勢(shì)，壓力降增大。

圖8 H-H流道主換熱區(qū)換熱面溫度云圖

圖9 H-H流道主換熱區(qū)換熱面壓力云圖

由于W-W流道波紋夾角小，兩股流體在溝槽中流動(dòng)交匯時(shí)的切向力小，同時(shí)該流道單位面積內(nèi)觸點(diǎn)少，所以湍流程度較弱，強(qiáng)化傳熱不充分，從圖10可以看出，W-W流道換熱面溫度分布不均勻，觸點(diǎn)處與周圍流道溫差明顯，傳熱效果差，但該流道流體動(dòng)量損耗??；從圖11可以看出，W-W流道進(jìn)出口壓力比較接近，損失不明顯，遞減趨勢(shì)很弱，壓力分布均勻，壓降小。

圖10 W-W流道主換熱區(qū)換熱面溫度云圖

圖11 W-W流道主換熱區(qū)換熱面壓力云圖

由于H-W流道波紋夾角介于H-H流道、W-W流道之間，所以兩股流體在溝槽中流動(dòng)交匯時(shí)的切向力和流道單位面積內(nèi)觸點(diǎn)數(shù)都介于H-H流道、W-W流道之間，湍流程度、流體動(dòng)量損耗也處于中間階段，從圖12可以看出，H-W流道換熱面的溫度均勻分布程度介于H-H流道、W-W流道之間，相比W-W流道傳熱效果較好，相比H-H流道傳熱效果較差；從圖13可以看出，H-W流道進(jìn)出口壓降的遞減趨勢(shì)介于H-H流道、W-W流道之間，與之前的速度場(chǎng)分析相吻合。

3 模擬數(shù)據(jù)對(duì)比

為了更具體地研究3種流道模型的傳熱和流動(dòng)阻力差異，為開發(fā)設(shè)計(jì)提供更多參考依據(jù)，選取熱交換器的板間流速在0.3～1.1m/s的工況進(jìn)行研究。

圖14、15為3種流道換熱面的傳熱系數(shù)和努塞爾數(shù)對(duì)比，從圖中可以看出，隨著流體雷諾數(shù)的增大，傳熱系數(shù)和努塞爾數(shù)也隨之增大，但趨勢(shì)逐漸趨于平緩，相比于H-W流道和H-H流道的傳熱系數(shù)提高了近10%，W-W流道的傳熱系數(shù)降低了11%左右，各流道的努塞爾數(shù)變化與傳熱系數(shù)相當(dāng)；圖16為3種流道換熱面的壓力梯度對(duì)比，從圖中可以看出，隨著流體雷諾數(shù)的增大，壓力梯度也隨之增加，且增加的趨勢(shì)越來越明顯，相比于H-W流道，H-H流道的壓力梯度提高了一倍之多，W-W流道的壓力梯度降低了80%左右；圖17為3種流道換熱面的歐拉數(shù)對(duì)比，歐拉數(shù)代表流體的動(dòng)量損失率，從圖中可以看出，H-H流道的動(dòng)量損失率最大，W-W流道的動(dòng)量損失率最小，并且隨著流體雷諾數(shù)的增大，3種流道換熱面的歐拉數(shù)均呈遞減趨勢(shì)。綜上所述，H-H流道特點(diǎn)是高傳熱系數(shù)、高阻力降，H-W流道特點(diǎn)是中等傳熱系數(shù)、中等阻力降，W-W流道特點(diǎn)是低傳熱系數(shù)、低阻力降，H-W流道換熱面的傳熱系數(shù)與壓力降均介于H-H流道、W-W流道之間，與之前的數(shù)值模擬分析相吻合。

a. H換熱面 b. W換熱面

a. H換熱面 b. W換熱面

圖14 3種流道換熱面的傳熱系數(shù)對(duì)比

圖15 3種流道換熱面的努塞爾數(shù)對(duì)比

圖16 3種流道換熱面的壓力梯度對(duì)比

圖17 3種流道換熱面的歐拉數(shù)對(duì)比

4 結(jié)論

4.1隨著人字形板片波紋夾角的增大，板間流體的流動(dòng)形式出現(xiàn)從“十字交叉流”向“曲折流”形式轉(zhuǎn)變的趨勢(shì)，所以除了波紋的深度、傾角、節(jié)距3項(xiàng)基本參數(shù)外，人字形波紋的夾角也是影響板片間流體流動(dòng)形態(tài)的因素之一。

4.2相比于H-W流道，H-H流道的膜傳熱系數(shù)提高了近10%，W-W流道的膜傳熱系數(shù)降低了11%左右，各流道的努塞爾數(shù)變化與傳熱系數(shù)相當(dāng)。

4.3相比于H-W流道，H-H流道的壓力損失提高了一倍之多，W-W流道的壓力損失降低了80%左右，相對(duì)于傳熱系數(shù)，人字形波紋夾角的改變對(duì)流道壓力損失的影響更明顯。

4.4隨著流體雷諾數(shù)的增大，各流道換熱面的傳熱系數(shù)和努塞爾數(shù)也隨之增大，但趨勢(shì)逐漸趨于平緩，增長(zhǎng)率不斷減小；各流道的壓力梯度隨之增加，且增加的趨勢(shì)越來越明顯，增長(zhǎng)率不斷加大；3種流道的流體動(dòng)量損失均隨雷諾數(shù)的增大而呈遞減趨勢(shì)。

[1] 秦叔經(jīng)，葉文邦.換熱器[M].北京：化學(xué)工業(yè)出版社，2003.

[2] 余建祖.換熱器原理與設(shè)計(jì)[M].北京：北京航空航天大學(xué)出版社，2006.

[3] 徐志明，王月明，張仲彬.板式換熱器性能的數(shù)值模擬[J].動(dòng)力工程學(xué)報(bào)，2011，31(3)：198～202.

[4] 蔡毅，賈志剛，周文學(xué)，等.人字形波紋板式換熱器性能數(shù)值模擬的研究[J].計(jì)算機(jī)與應(yīng)用化學(xué)，2009，26(1)：105～108.

[5] 崔立祺.基于FLUENT的板式換熱器三維數(shù)值模擬[D].杭州：浙江大學(xué)，2008.

[6] 張仲彬，董鵬飛，王月明，等.板式換熱器結(jié)構(gòu)優(yōu)化三維數(shù)值分析[J].化工機(jī)械，2013，40(2)：206～210.

[7] 李彥洲.板式換熱器板片換熱和流動(dòng)特性的研究[D].長(zhǎng)春：長(zhǎng)春工業(yè)大學(xué)，2014.

[8] 趙元東，張井志，張冠敏，等.熱混合板式換熱器換熱阻力性能[J].壓力容器，2013，30(8)：46～53.

[9] 陳文超，張鎖龍，梁欣.人字形板式換熱器雙流道模型的溫度場(chǎng)數(shù)值模擬[J].化工機(jī)械，2010，37(4)：465～468.

[10] 欒志堅(jiān)，張冠敏，張俊龍，等.波紋幾何參數(shù)對(duì)人字形板式換熱器內(nèi)流動(dòng)形態(tài)的影響機(jī)理[J].山東大學(xué)學(xué)報(bào)，2007，37(2)：34～37.

NumericalSimulationofChevron-typePlateHeatExchangerwithDifferentWaveAngles

WU Si1, LI Chun-lan1, WANG Sen1, LI Zhan-ying1, ZHU Hai-zhou2

(1.LanzhouLSEnergyEquipmentEngineeringInstituteCo.,Ltd.,Lanzhou730314,China; 2.LanzhouLSHeatExchangerEquipmentCo.,Ltd.,Lanzhou730314,China)

The main heat-exchanging zone’s model of three heat flow channels which constituted by chevron-type plates which boasting of 60°and 120° wave angles was built; and making use of CFD software to simulate both flow and heat transfer of the fluid within these three channels and to analyze the velocity field, temperature field and the pressure field there was implemented to show that, the inclined angle of chevron-type plate is one of factors which influencing the fluid’s flow pattern; and both heat transfer performance and pressure loss within channels become increased with the increase of chevron-type wave angles; and compared to the heat transfer performance, the influence of the change of wave angles on the flow channel’s pressure loss is more obvious.

plate heat exchanger, chevron-type plate, numerical simulation, heat transfer coefficient, pressure drop

*吳 私，男，1988年11月生，助理工程師。甘肅省蘭州市，730314。

TQ051.5

A

0254-6094(2016)06-0770-06

2015-11-24)