房立偉,焦守一,房立帥,劉帥

(1.山東魯潤熱能科技有限公司,山東 濟(jì)南 250303;2.山東大學(xué) 能源與動(dòng)力工程學(xué)院,山東 濟(jì)南 250061)

在能源動(dòng)力、冶金工程、化工生產(chǎn)等重要的領(lǐng)域中,換熱過程都是必不可少的一部分,為了提高生產(chǎn)效率,人們對強(qiáng)化波紋板的傳熱效率的研究,一直都是極為重要的課題。隨著強(qiáng)化傳熱技術(shù)的研究和發(fā)展,換熱裝置也在逐步換代,在不斷地改進(jìn)和優(yōu)化中,換熱效率也在逐漸上升,又因?yàn)榘迨綋Q熱器的體積并不大,易于安裝,價(jià)格低廉,換熱效率高,流阻低,尤其是其具有良好的換熱特性,能夠靈活地適應(yīng)生產(chǎn)中的熱負(fù)荷變化,所以更受業(yè)內(nèi)青睞,被廣泛應(yīng)用,例如在相同大小的換熱面積下,板式換熱器的換熱量,比管式換熱器多一倍以上,而波紋板比平直板的換熱能力又大了數(shù)倍之多。而且,不同的板型對于換熱器的熱效率影響也很大,經(jīng)過長期的理論研究,還有不斷的實(shí)踐,人們已經(jīng)確定了板式換熱器的常規(guī)板型,也就是波紋板,在我們?nèi)粘Ｉa(chǎn)中有一些比較常用的波紋板型,如三角型波紋板、正弦型波紋板、以梯形交替的梯型波紋板、弧形波浪的弧型波紋板等,在這些造型迥異的波紋表面,最常用于生產(chǎn)加工的換熱裝置中的,是三角型波紋板和正弦型波紋板,因?yàn)閾?jù)大量的實(shí)驗(yàn)研究分析,這兩種形式的波紋通道,其換熱效率和經(jīng)濟(jì)性指標(biāo)是相對較高的。

自19世紀(jì)70年代板式換熱器問世以來,已有一百四十余年歷史[1]。在20世紀(jì)70年代,外國的板式換熱器就已經(jīng)進(jìn)入了發(fā)展的黃金階段,已經(jīng)研發(fā)出了很多不同種類的板式換熱器[2],而在國內(nèi),由于生產(chǎn)力水平的限制和當(dāng)時(shí)國際形勢的復(fù)雜變化,對板式換熱器的研究,比國外落后了幾十年,在20世紀(jì)70年代才開始進(jìn)行相關(guān)技術(shù)的研究、設(shè)計(jì),開始運(yùn)用于一些工業(yè)生產(chǎn)領(lǐng)域,第一臺(tái)板式換熱器直到1965年才被成功制造出來[3],又過了二十余年,才在國內(nèi)得到了較大的發(fā)展[4],并擴(kuò)大到民用建筑,同時(shí)提高了在層數(shù)較高的建筑設(shè)施中的使用,并開始廣泛地被應(yīng)用到了空調(diào)制冷等領(lǐng)域中[5-6]。

隨著現(xiàn)代科技的突破,更高精度設(shè)備的開發(fā),生產(chǎn)制造技術(shù)的跨越式發(fā)展,新世紀(jì)以來,人們通過大量的實(shí)驗(yàn)和仿真模擬,整理總結(jié)出了大量的實(shí)驗(yàn)數(shù)據(jù)和理論實(shí)踐成果,推動(dòng)了波紋板的改進(jìn),使板式換熱器換熱技術(shù)更加成熟,而目前作為板式以及板翅式換熱器的基礎(chǔ)性板型——波紋板,更是被廣泛地研究,研究表明,其各種不同的波紋通道中,正弦型和三角型的波紋板,具有較高的換熱效率,同時(shí)投入生產(chǎn)使用的成本相對較低。另外,大量對換熱阻力特性的研究,也為波紋板的優(yōu)化和改進(jìn)提供了理論與實(shí)踐數(shù)據(jù)支持。

Islamoglu[7]研究了空氣為介質(zhì)的波紋通道,實(shí)驗(yàn)結(jié)果表明,隨著通道上下板間高度的增加,Nu和摩擦因子f增大。Kanaris[8]對以空氣為介質(zhì)的波紋微通道進(jìn)行了數(shù)值研究,分析了內(nèi)部的流場和應(yīng)力分布,得出了波紋通道提高了換熱系數(shù)的結(jié)論。Grijspeerdt[9]等對正弦波紋板進(jìn)行了兩維和三維的數(shù)值模擬,并根據(jù)可視化實(shí)驗(yàn)結(jié)果,對模擬結(jié)果進(jìn)行了定性驗(yàn)證,指出模擬結(jié)果可以確定近壁高溫區(qū),從而為板式換熱器的優(yōu)化設(shè)計(jì)提供了指導(dǎo)。陰繼翔[10]等通過數(shù)值模擬獲得了人字形波紋板的波紋間距、波紋深度等形狀參數(shù)與努塞爾數(shù)和壓降之間的關(guān)系曲線,得到了人字形波紋板片的最優(yōu)結(jié)構(gòu)參數(shù)。張冠敏[11]開發(fā)了一種復(fù)合波紋板式換熱器,并進(jìn)行了實(shí)驗(yàn)和數(shù)值模擬研究,結(jié)果表明此種類型的換熱器,在類似的換熱性能下,其流動(dòng)阻力更小。Fernandes[12]等建立了人字形正弦波紋板的一層流體模型,獲得了摩擦因子隨波紋高度比增大而增大的變化規(guī)律。馬學(xué)虎[13]指出波紋間距對傳熱效果的影響要高于波紋傾角和波紋高度。曲寧[5]提出了一種波浪形傳熱板片,并與人字形波紋板進(jìn)行了數(shù)值模擬,新板片具有相近的換熱性能,但壓降減小了50%以上。郭春生[14-15]等提出了一種不同波紋比例的波紋板式換熱器,并結(jié)合場協(xié)同理論進(jìn)行了換熱性能分析。

波紋板作為板式換熱器的核心部件,一直以來都是研究的重點(diǎn),何種參數(shù)的波紋板板型傳熱能力最強(qiáng),一直是人們不斷探索的課題,然而,由于影響波紋板傳熱能力的因素過多,不同因素的組合變換形式層出不窮,數(shù)據(jù)復(fù)雜多樣,對于波紋板的優(yōu)化改進(jìn)方向,依舊存在很大的研究空間,因此也一直是各國科學(xué)家、工程師研究的熱門課題,當(dāng)前尚存在以下問題:

1)在正弦型復(fù)合波紋板不同的結(jié)構(gòu)參數(shù)組合下,對氣氣換熱的流動(dòng)阻力特性、換熱因子、綜合換熱性能等方面的研究較少,仿真數(shù)據(jù)不夠全面。

2)在雷諾數(shù),努塞爾數(shù),進(jìn)出口速度等流動(dòng)參數(shù)的變化下,氣氣換熱的流動(dòng)阻力及換熱特性的對比數(shù)據(jù)較少。

1 數(shù)值模擬

波紋板作為介質(zhì)的流動(dòng)通道,承擔(dān)著介質(zhì)之間傳熱的中介的作用,在氣氣換熱的情況下,一定溫度的介質(zhì)從入口流入,通過波紋板的傳遞,從溫度基本恒定的外界介質(zhì)中吸收熱量,溫度升高,再從出口流出。不同的流動(dòng)參數(shù),如雷諾數(shù)、入口速度、入口壓力等,會(huì)導(dǎo)致?lián)Q熱能力的不同變化。因此,選定一組結(jié)構(gòu)參數(shù)作為基準(zhǔn),即選定一組確定的縱橫比、波長比、間距比等結(jié)構(gòu)參數(shù),建立波紋板模型,選取不同的入口雷諾數(shù),通過Fluent數(shù)據(jù)仿真模擬,對不同入口速度的介質(zhì)流過波紋板后,其溫度變化、摩擦因子、壓降損失、換熱性能等隨之變化的規(guī)律進(jìn)行分析。

使用建模軟件,根據(jù)不同的結(jié)構(gòu)參數(shù),分別進(jìn)行模型建立,把入口設(shè)為inlet,出口設(shè)為outlet,上下壁面設(shè)為up-down-wall,前后壁面設(shè)為front-back-wall,并根據(jù)模擬要求劃分網(wǎng)格,網(wǎng)格物理偏好改為CFD,在網(wǎng)格插入多區(qū)域方法,進(jìn)行源面設(shè)置,設(shè)置網(wǎng)格單元尺寸為0.15 mm,生成后,將膨脹應(yīng)用到特定網(wǎng)格邊界,對上下及前后壁面應(yīng)用平滑過渡的邊界條件,生成后模型如圖1,因?yàn)椴y板流道較為復(fù)雜,劃分網(wǎng)格時(shí)逐漸減小網(wǎng)格尺寸,細(xì)化調(diào)整網(wǎng)格,保證網(wǎng)格質(zhì)量。

圖1 基準(zhǔn)波紋板模型示意圖及網(wǎng)格情況

為檢驗(yàn)網(wǎng)格無關(guān)性,可不斷調(diào)整網(wǎng)格尺寸,當(dāng)努塞爾數(shù)在網(wǎng)格尺寸變化到一定程度時(shí),能夠穩(wěn)定在一定的范圍之內(nèi),不再有明顯的大范圍波動(dòng),就說明此時(shí)的網(wǎng)格密度符合數(shù)值模擬的精度要求。波紋板的網(wǎng)格劃分如圖1所示。

1.1 控制方程

根據(jù)波紋板的結(jié)構(gòu)變化,以及實(shí)際情況中的流動(dòng)換熱特性,本章進(jìn)行的數(shù)值模擬可以根據(jù)以下假設(shè):

1)在波紋板中的流動(dòng)和傳熱過程是穩(wěn)定的;

2)本章研究的是不可壓縮的、連續(xù)定常流動(dòng)的牛頓流體,對于此類流體,可以應(yīng)用無滑移的Navier-Stokes方程進(jìn)行求解;

3)介質(zhì)低速流動(dòng),忽略流體流動(dòng)中的浮升力,以及黏度性耗散產(chǎn)生的熱效應(yīng);

4)物性參數(shù)為常數(shù);

5)忽略壓力和溫度對流體實(shí)驗(yàn)定參數(shù),即比熱及黏度性系數(shù)等的影響;

6)波紋板上下壁面溫度恒定,設(shè)置為400 K;

7)假設(shè)入口流量是均勻的。

流體的連續(xù)性方程:

(1)

動(dòng)量方程:

(2)

能量方程:

(3)

1.2 湍流方程

波紋板中發(fā)生湍流,即流體內(nèi)部產(chǎn)生小漩渦,流層被破壞,每一點(diǎn)的無規(guī)則流體團(tuán),流動(dòng)狀態(tài)都在不斷改變,速度大小和方向隨機(jī)變化,部分流層間滑動(dòng)混合以復(fù)雜多變的無規(guī)律形式流動(dòng)。湍流會(huì)增大阻力,影響實(shí)驗(yàn)數(shù)值的精確測量,因此,在Fluent中使用Realizablek～ε模型,減小湍流對數(shù)值模擬的影響。

Realizablek～ε模型的湍動(dòng)能方程為:

(4)

Realizablek～ε模型的湍動(dòng)耗散率方程為:

(5)

在Fluent設(shè)置中,介質(zhì)入口邊界選擇速度入口,其中把入口溫度設(shè)置為300 K,設(shè)置入口雷諾數(shù)分別為2 000,4 000,6 000,8 000,10 000,入口的湍流強(qiáng)度為5%,湍流黏度度比為10。在出口邊界條件,選擇抑制回流,能在出口出現(xiàn)回流時(shí),有更好的收斂速度。對于波紋板的上下壁面,采用恒定的溫度邊界條件,溫度設(shè)定為400 K,無滑移速度。求解方案選用SIMPLEC算法,求解數(shù)據(jù)更加精確,對動(dòng)量、湍流動(dòng)能、湍流耗散率及能量,均采用SOU格式求解。

2 結(jié)果與討論

2.1 雷諾數(shù)對換熱性能的影響

圖2和圖3為計(jì)算所得基準(zhǔn)波紋板的努塞爾數(shù)Nu和換熱因子j隨雷諾數(shù)Re的變化趨勢。由圖可知,Nu隨著Re的增大而增大,因此Re越大,對流換熱越強(qiáng)。對于充分發(fā)展的湍流流動(dòng),流速的增大會(huì)增大湍流度,進(jìn)而強(qiáng)化近壁面流體與核心區(qū)域流體的摻混傳熱傳質(zhì)過程。這是高Re數(shù)下,換熱系數(shù)較高的主要原因。相比其他類型的板式換熱器,正弦波紋通道的換熱性能較好,具有很好的應(yīng)用前景。對于本文的研究結(jié)果,Nu基本與Re呈現(xiàn)線性增大的規(guī)律。主要的原因在于:充分發(fā)展的湍流狀態(tài)下,換熱系數(shù)基本與Re的0.8～1次方成比例,相關(guān)的公式可以由圓管內(nèi)部充分發(fā)展湍流狀態(tài)下的DB公式獲得。雖然正弦波紋板片與圓管內(nèi)部的強(qiáng)制對流換熱特性有一定的差異,但其本質(zhì)上都屬于受限空間內(nèi)的強(qiáng)迫對流,湍流對于換熱特性的影響規(guī)律基本類似。基于此規(guī)律,也有很多的波紋板式換熱器的預(yù)測公式,可以指導(dǎo)工業(yè)設(shè)計(jì)和生產(chǎn)以及結(jié)構(gòu)參數(shù)的優(yōu)化。

圖2 基準(zhǔn)波紋板的努塞爾數(shù)Nu隨雷諾數(shù)Re的變化

圖3 基準(zhǔn)波紋板的換熱因子j隨雷諾數(shù)Re的變化

而換熱因子j的變化趨勢不同,j隨著Re的增大而減小,且隨著Re增大,j的變化量越來越小?？芍猂e對j的影響隨著Re的增大而減弱,可以認(rèn)為當(dāng)Re足夠大,即流速足夠大的時(shí)候,在流道達(dá)到完全湍流狀態(tài),j的變化趨于穩(wěn)定,換熱能力趨于穩(wěn)定。換熱因子j是換熱器設(shè)計(jì)中的重要參數(shù),從其定義可以得出,j與Re的倒數(shù)成正比,也就對應(yīng)著j隨著Re的增大而逐漸降低。在低Re工況,由于Nu隨Re的增大趨勢較為明顯,此時(shí)j隨Re增大而降低的趨勢更加顯著。在高Re工況,由于Nu隨Re的增大趨勢較為平緩,j隨Re增大而降低的趨勢也逐漸趨于平穩(wěn)。換熱因子j是一個(gè)綜合的無量綱參數(shù)組合,主要由Re、Nu和Pr決定,其值的大小也更能反映換熱器的綜合性能,比單純的Nu所包含的影響因素更為全面,推薦在相關(guān)的波紋板式換熱器優(yōu)化設(shè)計(jì)過程中,考慮采用j因子作為優(yōu)化指標(biāo)。

圖4為速度矢量云圖,可以看出正弦波紋板式換熱器內(nèi)部的速度矢量變化規(guī)律。由圖可得,工質(zhì)沿著波紋表面呈現(xiàn)曲折流動(dòng),主要的原因在于正弦波紋板片對于氣流的限制作用。在波紋板片的限制作用下,流體只能沿著波紋方向向前流動(dòng)。與傳統(tǒng)的人字形波紋板片相比,正弦波紋板片的流場變化更加光順,流體主要沿著主流方向呈現(xiàn)蛇形彎曲流動(dòng)。由于正弦波紋板式換熱器不存在波紋傾角,氣流基本沿著通道的長度方向流動(dòng),在寬度方向的折返流動(dòng)較少。此種流型,雖然對換熱特性來說,缺少縱向旋渦的擾動(dòng),換熱性能有一定的損失。但在壓降方面卻有重要的優(yōu)勢,由于流體順著波紋流動(dòng),導(dǎo)致流程相對較短。同時(shí),無縱向旋渦的擾動(dòng),也進(jìn)一步降低了正弦波紋板式換熱器的壓降損失。導(dǎo)致了本文研究的正弦波紋板式換熱器具有較好的綜合性能,即傳熱較好,壓降不大的特性,更適合用于氣氣換熱裝置中需要降低壓降的場景。尤其是針對鋼鐵行業(yè)的煤氣預(yù)熱器,在實(shí)際生產(chǎn)過程中,要求板式換熱器的壓降損失必須非常小,不能阻礙煤氣進(jìn)入煉鋼爐膛,以免出現(xiàn)停機(jī)的現(xiàn)象。對于此種換熱需求,正弦波紋板式換熱器可以很好的滿足。

圖4 基準(zhǔn)波紋板速度矢量圖

2.2 雷諾數(shù)對阻力性能的影響

圖5為進(jìn)出口壓降ΔP隨Re的變化規(guī)律?？芍隨著Re的增大而增大,且隨著Re的增大,ΔP呈增大趨勢。板式換熱器的壓降損失是其優(yōu)化設(shè)計(jì)的重要參數(shù),也是決定其適用工況的重要指標(biāo)。對于氣氣換熱器,其應(yīng)用領(lǐng)域常為生產(chǎn)過程中的余熱回收。由于換熱器與整個(gè)生產(chǎn)流程串聯(lián),必須要求經(jīng)過板式換熱器的壓降損失要足夠小,不能影響整個(gè)生產(chǎn)流程。對于本文所研究的正弦波紋板式換熱器,由于優(yōu)秀的板片形式,其壓降損失較小。由圖6可得,在本文所研究的工況范圍內(nèi)(2 000≤Re≤10 000),正弦波紋板式換熱器的壓降損失在50～450。根據(jù)實(shí)際的換熱需求,選定合理的流速,可以保證正弦波紋板式換熱器滿足相關(guān)的工業(yè)生產(chǎn)需求。與傳統(tǒng)的人字形波紋板式換熱器相比,正弦波紋板式換熱器對應(yīng)的壓降損失更小,更適合作為氣氣換熱器,應(yīng)用于生產(chǎn)過程中氣態(tài)余熱介質(zhì)的回收,提升整個(gè)生產(chǎn)流程的經(jīng)濟(jì)性,降低能耗,實(shí)現(xiàn)節(jié)能減排,助力“雙碳”戰(zhàn)略目標(biāo)的實(shí)現(xiàn)。

圖5 基準(zhǔn)波紋板的進(jìn)出口壓降ΔP隨雷諾數(shù)Re的變化

圖6 基準(zhǔn)波紋板的摩擦因子f隨雷諾數(shù)Re的變化

圖6為摩擦因子f隨Re的變化規(guī)律,可知f隨著Re的增大而減小,即Re越大,摩擦阻力越小,且隨著Re的增大,f變化趨勢減小。f的變化趨勢與j因子類似,均為隨著Re的增大而逐漸降低。與j因子的定義類似,f也是一個(gè)組合的無量綱參數(shù),低Re下,Re對于壓降的增大作用較為明顯,導(dǎo)致f較高,隨著Re的增大,流速對于f的削弱作用更為明顯,導(dǎo)致f逐漸降低。

3 結(jié)論

利用數(shù)值模擬的方法,研究了正弦波紋通道內(nèi)氣氣流動(dòng)換熱特性,主要結(jié)論如下:

1)在一定條件內(nèi),雷諾數(shù)增大,努塞爾數(shù)和對流換熱系數(shù)不斷增大,對于充分發(fā)展的湍流流動(dòng),流速的增大會(huì)增大湍流度,進(jìn)而強(qiáng)化近壁面流體與核心區(qū)域流體的摻混傳熱傳質(zhì)過程。相比其他類型的板式換熱器,正弦波紋通道換熱性能較好,具有很好的應(yīng)用前景。

2)工質(zhì)沿著波紋表面呈現(xiàn)曲折流動(dòng),主要的原因在于正弦波紋板片對于氣流的限制作用。在波紋板片的限制作用下,流體只能沿著波紋方向向前流動(dòng)。與傳統(tǒng)的人字形波紋板片相比,正弦波紋板片的流場變化更加光順,流體主要沿著主流方向呈現(xiàn)蛇形彎曲流動(dòng)。導(dǎo)致了本文研究的正弦波紋板式換熱器具有較好的綜合性能,即傳熱較好,壓降不大的特性,更適合用于氣氣換熱裝置中需要降低壓降的場景。

3)進(jìn)出口壓降增大,且曲線變陡;進(jìn)出口溫差逐漸減小,趨勢變小;對流換熱系數(shù)變大;波紋板的換熱因子和摩擦因子都減小,最后趨于穩(wěn)定。