全俊威 辛林斌 孫健

摘 要：間壁換熱結(jié)構(gòu)是隧道窯緩冷段的重要結(jié)構(gòu)，進一步了解其流動和傳熱能力具有重要意義。為了研究隧道窯間壁內(nèi)流動及傳熱特性，建立了隧道窯間壁的三維物理和數(shù)學(xué)模型。并采用計算流體力學(xué)（CFD）軟件對其進行了數(shù)學(xué)模擬，研究了分層間壁結(jié)構(gòu)、分層通道的大小和位置、內(nèi)壁溫度、出口抽力等因素對間壁換熱流動及換熱特性的影響。結(jié)果表明，分層結(jié)構(gòu)間壁可以有效提高換熱效率，其中M3-1的換熱性能最佳;在分層間壁層數(shù)一定的情況下，分層通道的大小和位置對于間壁整體的換熱性能影響不大;隨著風(fēng)機抽力的增大，間壁換熱效果更佳，但間壁的分層結(jié)構(gòu)對于換熱效果的影響降低，不同分層結(jié)構(gòu)的間壁熱風(fēng)出口溫差變小。

關(guān)鍵詞：隧道窯;間壁結(jié)構(gòu);傳熱

1 引 言

隧道窯是陶瓷生產(chǎn)中一種不可或缺的連續(xù)式燒成設(shè)備，在燒制一些體型較大的陶瓷，如衛(wèi)浴、大型花瓶等制品方面有著獨特的優(yōu)勢。間壁換熱結(jié)構(gòu)作為緩冷段一種重要調(diào)整手段，對于陶瓷制品冷卻，尤其是大型制品冷卻有著較為關(guān)鍵的作用。但是長期以來，研究人員和窯爐操控者對于間壁換熱的重視程度不高，在制品冷卻方面，更多的研究集中在制品的碼放方式及密度，抽熱風(fēng)口位置的布置等方面，對于隧道窯間壁換熱方面的研究較少[1～2]。

本文通過對某陶瓷隧道窯的間壁換熱的尺寸及部分熱工參量進行了測定，利用計算流體動力學(xué)軟件FLUENT，對換熱器內(nèi)部的流場和溫度場進行了數(shù)值模擬。根據(jù)模擬結(jié)果得到熱風(fēng)出口溫度，間壁內(nèi)流速等參量及計算區(qū)域內(nèi)的質(zhì)點速度矢量圖、溫度場云圖等。綜合分析，符合流體流動和傳熱規(guī)律，并在此基礎(chǔ)上提出了間壁結(jié)構(gòu)的優(yōu)化，有效提高了間壁換熱的換熱性能，同時可以在一定程度上提高窯爐熱效利用率[3～5]。

2 物理模型的建立

本文選取了陶瓷隧道窯緩冷段的一節(jié)窯墻作為研究對象，模型結(jié)構(gòu)的尺寸：長×寬×高=4250×1780×454 （mm），間壁的結(jié)構(gòu)尺寸=3850×1005×230 （mm）。如圖1-a所示，其中圖案填充部分為窯墻，空白部分為間壁換熱流動區(qū)域。如圖1-b所示為間壁的三維模型圖，窯墻由四種耐火材料砌筑而成，由內(nèi)及外所用耐火材料及耐火材料厚度分別是，輕質(zhì)高鋁磚114 mm，輕質(zhì)黏土磚230 mm，耐火纖維氈100 mm，硅鈣板10 mm。在間壁工作時，冷空氣由窯墻外壁上的冷空氣入口被吸入，在間壁內(nèi)完成熱交換后由抽熱風(fēng)機經(jīng)由熱空氣出口抽出。

為了探究分層間壁結(jié)構(gòu)對于間壁內(nèi)流體流動及換熱特性的影響，在圖1間壁的基礎(chǔ)上，保持間壁的體積，窯墻材料及進出口面積，位置不變的情況，設(shè)計了兩種不同的分層間壁結(jié)構(gòu)，即圖2-a、2-b和圖3-a、3-b所示。同時為了探究分層通道的大小及位置對于分層間壁內(nèi)流體流動及換熱特性的影響，在圖2-a、2-b和圖3-a、3-b所示分層間壁的基礎(chǔ)上，提出了圖2-c、3-c和3-d三種間壁結(jié)構(gòu)。為了方便后續(xù)對于模型的論述及分析，后文將用“MX-X”代表不同間壁模型，“M1”代表模型一的間壁結(jié)構(gòu)，整個間壁不分層，即圖1所示;“M2”代表模型二的間壁結(jié)構(gòu)，如圖2所示，整個間壁分為兩層，其中圖2-a和2-b所示的結(jié)構(gòu)為M2-1，圖2-c所示結(jié)構(gòu)為M2-2;“M3”代表模型三的間壁結(jié)構(gòu)，如圖3所示，整個間壁分為三層，其中圖3-a和3-b所示的結(jié)構(gòu)為M3-1，圖3-c所示結(jié)構(gòu)為M2-2，圖3-d所示結(jié)構(gòu)為M3-3。

3 間壁換熱的數(shù)學(xué)描述

3.1? 控制方程

進行數(shù)值模擬時，簡化計算模型并做出假設(shè)：

（1）對于穩(wěn)定運行的隧道窯，窯內(nèi)溫度場和速度場不隨時間變化，研究其穩(wěn)定狀態(tài)下間壁內(nèi)溫度場和速度場分布。

（2）流體為不可壓縮牛頓流體，耐火材料的導(dǎo)熱系數(shù)及比熱容為定值或與溫度成一次線性關(guān)系、無內(nèi)熱源，具體物性參數(shù)如表1所示。

（3）在工程問題中利用空氣對流進行強制換熱時均為湍流流動，在本研究中空氣在間壁內(nèi)的流動也為三維穩(wěn)態(tài)湍流運動。

迪卡爾直角坐標(biāo)系三維定常不可壓湍流流動控制方程如下[7～8]：

其中，i，j=1，2，3，分別表示表示坐標(biāo)軸x，y，z方向;Fi 為體積力;T 為溫度;ρ，μ，λ，Cp 分別為流體密度、運動黏度、導(dǎo)熱系數(shù)、比熱。

本文計算模擬采用Realizable k-ε模型進行模擬，也使用了壁面函數(shù)。為使模型的控制方程容易進行分析計算，即建立其通用形式，通用方程的表達式如下[6]。

式中：? ? 是通用變量?？纱韚、v、w、T等求解變量;? ? 是廣義擴散系數(shù); S是廣義源項。

Realizable k-ε模型中的k 和ε 的方程如下：

3.2? 邊界條件

在計算中，除窯墻內(nèi)外壁面外，其余窯墻外壁面均視為絕熱壁面。

3.3? 網(wǎng)格劃分和模型的求解

3.3.1 網(wǎng)格的劃分

采用控制容積法將控制方程離散化，用ANSYS ICEM軟件繪制網(wǎng)格。圖4為計算模型整體網(wǎng)格，采用六面體結(jié)構(gòu)化網(wǎng)格。為了使得計算結(jié)果更加準(zhǔn)確，對流體進口，出口及整個流體域進行加密處理。如圖4-b所示，由于窯墻材料厚度跨度較大，窯墻最外層材料僅有10 mm，為了保證網(wǎng)格質(zhì)量和計算結(jié)果的準(zhǔn)確性，在Z軸方向網(wǎng)格劃分時對該部分也進行了加密處理，加密結(jié)果如圖4-c所示。

3.3.2 模型的求解

選用不可壓縮粘性流體的三維、穩(wěn)態(tài)、湍流流動與傳熱模型作為間壁內(nèi)流體流動和傳熱的物理模型。數(shù)學(xué)模型中的控制方程包括連續(xù)性方程、動量方程、能量方程。利用FLUENT軟件對間壁內(nèi)流場及溫度場進行數(shù)值模擬，模擬采用3D雙精度求解器，湍流模型選用Realizable k-ε雙方程模型，壓力修正選用SIMPLE算法，除湍流能量和湍流耗散率選用一階迎風(fēng)式外，其余各項均采用二階迎風(fēng)式。當(dāng)連續(xù)性方程、動量方程和能量方程變量殘差均小于且不發(fā)生變化時，認為計算收斂。

3.3.3 邊界條件及物性參數(shù)的設(shè)置

本文選取了隧道窯緩冷段的一段窯墻進行研究，因此除了窯墻內(nèi)外壁面外，窯墻側(cè)面全部設(shè)定為絕熱邊界層。同時，為了使得模擬更加準(zhǔn)確，將熱空氣出口設(shè)置為壓力出口邊界條件，冷空氣進口設(shè)置壓力入口邊界條件。模型中所涉及的材料物性參數(shù)如表1所示。

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

為了探究分層通道的大小和位置對于同一種分層間壁換熱流動及傳熱特性的影響，本文對M2和M3兩種分層間壁的四個結(jié)構(gòu)在內(nèi)壁溫度為773K的工況下進行了數(shù)值計算，表2是不同結(jié)構(gòu)模型的熱風(fēng)出口溫度。對于M2模型而言，M2-1和M2-2溫差始終約為1 K，分析其原因是當(dāng)間壁層數(shù)為兩層時，開口的大小和位置對于間壁內(nèi)流體流動的擾流效果并不明顯，因此換熱效果基本一致;對于M3結(jié)構(gòu)而言，M3-1、M3-2和M3-3三者溫差約在3～4K，其中M3-3模型換熱效果最佳，M3-1和M3-2換熱效果基本一致。綜合考慮換熱效果和窯爐的施工砌筑，M2-1和M3-1結(jié)構(gòu)較為簡單，且換熱效果也較好，是同一分層結(jié)構(gòu)中的最優(yōu)結(jié)構(gòu)。

4.1? 溫度場分析

如圖5所示，是M1、M2-1和M3-1模型在抽力為100 Pa，內(nèi)壁溫度為773 K 工況下，間壁中心的溫度云圖。圖5-a可以看出M1模型的間壁中心溫度分布較為均勻，在靠近出口上部的溫度較高約在330 K，其他大部分區(qū)域溫度約為315 K，溫度最低的位置在靠近冷空氣入口位置，間壁的平均溫度;從圖5-b可以看出M2-1模型的換熱效果較好，間壁中心大部分區(qū)域溫度約為320 K，比M1模型中心溫度高出約5 K。間壁的高溫區(qū)域在上層間壁的最左側(cè)，間壁內(nèi)最大溫差約為25 K，間壁的平均溫度;圖5-c和5-d可以看出M3-1模型的溫度分布均勻，上中下三層沒有明顯溫差，僅在上層靠近出口的位置溫度較高，比其他位置高出約10 K，間壁的平均溫度。由此可以得出，M3-1結(jié)構(gòu)的間壁，換熱性能最好，換熱均溫性也最好。

4.2? 流場分析

如圖6所示，是M1、M2-1和M3-1模型在抽力為100 Pa，內(nèi)壁溫度為773 K工況下，間壁中心的流場圖。從圖6-a可以看出M1的流場整體流速較小，空氣擾動也不明顯，僅在流場的中上部，出現(xiàn)兩個速度較低的小型渦流;從圖6-b可以看出M2-1流場較M1而言更為復(fù)雜，流場中出現(xiàn)了多個渦旋區(qū)，空氣擾動較為明顯。同時由于分層結(jié)構(gòu)的影響，流體在通過分層通道時，流速明顯增大，進一步加強了空氣擾動;從圖6-c和6-d可以看出M3-1流場的流速整體分布較為均勻，沒有出現(xiàn)較為明顯的渦流區(qū)和低速區(qū)，流速相較與M1、M2-1也明顯較大。從流場的流速大小和速度分布的綜合情況可以看出，M3-1的流場更加穩(wěn)定且流速更大，進一步證明M3結(jié)構(gòu)性能優(yōu)于其它兩種結(jié)構(gòu)，這與溫度場分析結(jié)論一致。

4.3? 換熱性能分析

如圖7，8所示，是M1、M2-1和M3-1模型在不同工況下，熱風(fēng)出口溫度和間壁內(nèi)平均努賽爾數(shù)隨風(fēng)機抽力的變化情況。隨著風(fēng)機抽力的增大，平均努塞爾增數(shù)大，間壁的整體換熱效果更佳，但是由于流量的增大，導(dǎo)致熱風(fēng)出口溫略有下降;同時三種模型之間的平均努塞爾數(shù)和熱風(fēng)出口溫度之間的差值不斷變小，這是因為隨著風(fēng)機抽力的不斷增大間壁內(nèi)的紊流程度不斷增強，間壁的分層結(jié)構(gòu)對流體擾動的影響逐漸變小。

從熱風(fēng)出口溫度來看，在相同的壁溫和風(fēng)機抽力的工況下，M3-1結(jié)構(gòu)的間壁熱風(fēng)出口溫度高于M1和M2-1結(jié)構(gòu)的間壁熱風(fēng)出口溫度;對于M2結(jié)構(gòu)的間壁，內(nèi)壁溫度為973 K和773 K時，其換熱性能較M1更佳;當(dāng)內(nèi)壁溫度為573 K時，其換熱效果與M1基本一致。

5 結(jié)? 論

通過對隧道窯間壁換熱內(nèi)部溫度場和流場數(shù)值模擬結(jié)果的分析可以得出以下結(jié)論：

（1）分層結(jié)構(gòu)間壁可以有效增強間壁換熱效果，其中三層結(jié)構(gòu)的M3-1間壁模型換熱效果最優(yōu)，熱風(fēng)出口溫度較M1結(jié)構(gòu)間壁提高了約20K。

（2）在分層結(jié)構(gòu)一定的情況下，改變分層通道的位置和大小對于間壁換熱的整體換熱效果影響不大，熱風(fēng)出口溫差約在3K。

（3）風(fēng)機抽力對換熱效果有著顯著影響，當(dāng)風(fēng)機抽力從80Pa增大至16Pa時，熱風(fēng)出口溫差由最大近20K減小至約1K，此時間壁的分層結(jié)構(gòu)對于間壁換熱效果的影響較小。

參考文獻

[1]. 胡國林， 周露亮， 陳功備. 陶瓷工業(yè)窯爐[M]. 武漢：武漢理工大學(xué)出版社，2001.

[2] 許寧，王雪峰，張柱銀，王君，黃建軍，肖慧，郅曉.隧道窯冷卻帶流-固耦合散熱數(shù)值模擬[J].中國陶瓷工業(yè)，2016，23（02）：26-30.

[3] 李杰，胡國林，孫健，朱永紅.某隧道窯煙氣換熱器換熱性能的研究[J].陶瓷學(xué)報，2013，34（04）：482-486.

[4] 陳靜.隧道窯冷卻帶抽氣系統(tǒng)的改進[J].磚瓦，2018（12）：51-52.

[5] 王洪利，黃杰，劉建雄.陶瓷廠隧道窯熱平衡測試與節(jié)能分析[J].中國陶瓷，2013，49（02）：35-38.

[6] 《ANSYS ICEM CFD工程實例詳解》胡坤，李振北. ANSYS ICEM CFD工程實例詳解[M].北京：人民郵電出版社， 2014.

[7] 《數(shù)值傳熱學(xué)》陶文銓編著.數(shù)值傳熱學(xué)[M]. 西安交通大學(xué)出版社 ， 2001

[8] 楊世銘， 陶文栓. 傳熱學(xué)[M]. 北京： 高等教育出版社，2006.