煙氣橫掠麻面管束顆粒沉積特性的數(shù)值分析

閆順林， 王皓軒

(華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003)

煙氣橫掠麻面管束顆粒沉積特性的數(shù)值分析

閆順林， 王皓軒

(華北電力大學(xué) 能源動(dòng)力與機(jī)械工程學(xué)院，河北 保定 071003)

為了減輕電廠鍋爐尾道內(nèi)煙氣對(duì)換熱管的磨損及積灰情況，特提出和建立了一種新型的管型——麻面管，并通過(guò)選擇適用的磨損計(jì)算和顆粒沉積模型，針對(duì)煙氣橫掠順列麻面管束進(jìn)行氣固兩相流動(dòng)的數(shù)值模擬，分析了煙氣物性參數(shù)和管束結(jié)構(gòu)參數(shù)對(duì)顆粒沉積特性的影響，并比較了相同條件下麻面管和光管的積灰情況。結(jié)果表明:等同工況下，麻面管表面顆粒沉積率和磨損量均比光管的低；較大煙速情況下，隨著顆粒平均粒徑的增大，磨損增加幅度較大；確定了當(dāng)s1/D=3.0，s2/D=1.8時(shí)，顆粒沉積率最低，為20.72%；對(duì)于一般的較小粒徑顆粒，e=0.75 mm的麻面管管束在減輕積灰上效果較其他幾種情況更好，另外，麻面管表面坑深度越大，飛灰顆粒對(duì)麻面管的磨損就越輕。

麻面管； 氣固兩相流； 顆粒沉積； 磨損； 數(shù)值模擬

0 引言

換熱管作為一種換熱元件，被廣泛應(yīng)用于鍋爐水冷壁、省煤器、凝汽器等換熱設(shè)備中。由于工業(yè)鍋爐中煙氣含塵量大，且煙氣流動(dòng)具有很大的隨機(jī)性，使得飛灰對(duì)換熱管的沖擊碰撞點(diǎn)也具有相對(duì)的隨機(jī)性，很容易造成煙道內(nèi)換熱設(shè)備的積灰和磨損,換熱器積灰問(wèn)題會(huì)導(dǎo)致其傳熱效率的下降和壓降損失的增加,同時(shí)也增加了煙氣流動(dòng)和設(shè)備運(yùn)行的不穩(wěn)定性,易引起安全隱患[1]。換熱管的積灰和磨損將嚴(yán)重影響到電廠設(shè)備運(yùn)行的安全性和經(jīng)濟(jì)性[2-3]。省煤器作為鍋爐煙道內(nèi)重要的換熱元件，研究省煤器管束在含塵煙氣中的積灰特性,提出有效的預(yù)防積灰措施, 對(duì)于各種設(shè)備的經(jīng)濟(jì)和安全運(yùn)行有著重要的意義。

雖然國(guó)內(nèi)外存在大量關(guān)于換熱管的磨損模擬研究，但僅有少部分學(xué)者從實(shí)驗(yàn)層面上研究了煙氣換熱器的積灰特性，很難解釋積灰機(jī)理和預(yù)測(cè)積灰情況。岑可法等[4]通過(guò)大量的實(shí)驗(yàn)研究，闡釋了煙氣中攜帶的飛灰顆粒對(duì)鍋爐換熱管的沖蝕磨損機(jī)理及顆粒沉積影響因素；董自翔等[5]采用格子Boltzmann與格子氣自動(dòng)機(jī)的耦合模型模擬分析了單排翅片管束外表面的顆粒沉積規(guī)律，得到了粒徑為5 μm和10 μm的顆粒在管道表面的沉積形態(tài)；Han等人[6]數(shù)值研究了管排換熱器表面的積灰特性, 并分析了顆粒粒徑、流速、管排形狀及布置方式等對(duì)積灰特性的影響；文獻(xiàn)[7-9]分別針對(duì)鰭片管、橢圓管、螺旋管等管型做了顆粒碰撞和磨損的數(shù)值研究，得出這些管型與光管相比，能有效均流煙氣流場(chǎng)、減輕磨損。

目前，現(xiàn)有的關(guān)于換熱管積灰及磨損的大量研究成果，主要集中在對(duì)光管、膜式管束、H翅片管束等管型上。但是考慮到翅片等安裝上的難度，所以，本文特提出和建立了一種新型的管型——麻面管，并通過(guò)選擇適用的計(jì)算模型，針對(duì)煙氣橫掠順列麻面管束進(jìn)行氣固兩相流動(dòng)的數(shù)值模擬研究，為這種新式管型——麻面管管束的應(yīng)用和推廣提供參考價(jià)值。

1 物理模型和邊界條件

1.1 幾何模型

現(xiàn)有的研究已經(jīng)表明，錯(cuò)列布置的管束磨損要比順列布置管束嚴(yán)重[10-11]，因此，本文研究麻面管束磨損情況時(shí)采用5排順列布置管束。

為簡(jiǎn)化模型，本文建模選取一段長(zhǎng)度L=46 mm，直徑D=48 mm的管子，通過(guò)設(shè)置物理模型的對(duì)稱性和周期性邊界，進(jìn)行數(shù)值模擬計(jì)算。圖1為麻面管的結(jié)構(gòu)示意圖，管子表面排列的為深度d=1 mm、坑口直徑為4.47 mm的凹坑，凹坑采用錯(cuò)列緊密排列方式布置，凹坑設(shè)計(jì)來(lái)源于高爾夫球設(shè)計(jì)理念，坑口采用倒角分段平滑過(guò)渡。計(jì)算區(qū)域和邊界的設(shè)置如圖2所示，管軸向設(shè)置為周期性邊界，其中s1、s2分別表示管束的橫向、縱向管距，且s1/D=2.5,s2/D=1.8。

圖1 麻面管結(jié)構(gòu)示意圖

圖2 物理模型和邊界條件

1.2 邊界條件的設(shè)置

計(jì)算域采用ICEM進(jìn)行結(jié)構(gòu)化網(wǎng)格劃分，在網(wǎng)格劃分過(guò)程中采用了切塊方式，并將管道區(qū)域從計(jì)算域中獨(dú)立切分出來(lái)，再進(jìn)行O型剖分，且在麻面管管束表面附近區(qū)域采用Linear方式進(jìn)行由密及疏的加密處理，以簡(jiǎn)化計(jì)算量、提高準(zhǔn)確度。網(wǎng)格劃分過(guò)程中，通過(guò)對(duì)不同網(wǎng)格尺寸及數(shù)目下對(duì)計(jì)算結(jié)果的比較，驗(yàn)證了網(wǎng)格數(shù)量對(duì)計(jì)算結(jié)果無(wú)關(guān)性，最終確定網(wǎng)格數(shù)量約為68萬(wàn)。

氣相邊界條件:設(shè)定速度入口邊界條件，大小為6～10 m/s，入口煙氣流量為0.5 kg/m3,粘度系數(shù)u=3.328×10-5Pa·s煙氣溫度為常量690 K，壁面全部設(shè)定為無(wú)滑移邊界墻壁，出口為自由出流邊界條件。

顆粒相邊界條件：設(shè)定顆粒入口速度與氣相速度相同，并以面源的方式噴入，顆粒入口質(zhì)量流量為0.003 kg/s，密度為2 000 kg /m3，直徑為0～120 μm，符合Rosin-Rammler分布。顆粒與管束碰撞表面設(shè)為壁面反彈邊界條件，顆粒與壁面間的相互作用通過(guò)恢復(fù)系數(shù)來(lái)描述。

2 數(shù)值模型及計(jì)算方法

2.1 連續(xù)相控制方程

本文研究管束磨損的過(guò)程中，煙氣可視為不可壓縮流體，忽略溫度變化的影響。在控制方程的選擇上，整個(gè)過(guò)程都要滿足質(zhì)量守恒定律、動(dòng)量守恒定律以及能量守恒定律。目前Ansys Fluent15.0提供的湍流模型主要有標(biāo)準(zhǔn)k-ε模型，RNGk-ε模型，Realizablek-ε模型。由于各種模型都有各自的適用范圍[12]，考慮到煙氣沖刷管束的流動(dòng)中不會(huì)產(chǎn)生強(qiáng)旋流，本文研究中選擇應(yīng)用最為廣泛的標(biāo)準(zhǔn)k-ε模型和標(biāo)準(zhǔn)壁面函數(shù)模型。連續(xù)相湍流運(yùn)動(dòng)的基本方程如下：

連續(xù)性方程

(1)

式中：u為煙氣的平均流速，m/s；i=1、2、3，代表速度沿x、y、z軸的分量。

動(dòng)量守恒方程

(2)

式中：P為煙氣平均壓力，Pa；ρ為煙氣密度，kg/m3；μ為煙氣動(dòng)力粘度，Pa·s；t為時(shí)間，s；j=l、2、3，為沿x、y、z軸的分量。

能量守恒方程

(3)

式中：λ為煙氣的導(dǎo)熱系數(shù)，W/(m2·K)；T為煙氣溫度，K；Cp為煙氣比熱容，J/(kg·K)。

湍動(dòng)能方程k

(4)

湍動(dòng)能耗散率方程ε

(5)

式中：Gk為湍能體積生成率；

表1 系數(shù)取值參考

2.2 顆粒運(yùn)動(dòng)方程

本文中把顆粒相作為稀疏相，但顆粒密度遠(yuǎn)遠(yuǎn)大于煙氣密度，沒(méi)有較大的溫度梯度和壓力梯度存在，同時(shí)又忽略顆粒間的碰撞和旋轉(zhuǎn)，所以認(rèn)為顆粒主要受重力、曳力和慣性力的作用。因此，單個(gè)顆粒的運(yùn)動(dòng)方程可以寫(xiě)為：

(6)

式中：m為單個(gè)顆粒的質(zhì)量；FD為單位質(zhì)量的顆粒所受曳力；Fx為x方向的其他作用力。

(7)

式中：μ是動(dòng)力粘度，Pa·s;ρp是顆粒密度，kg/m3;dp是顆粒直徑，m;Re是相對(duì)雷諾數(shù)。

曳力系數(shù)是由Schiller-Naumann曳力法則來(lái)決定[13]，見(jiàn)公式(8)：

(8)

2.3 顆粒沉積及磨損計(jì)算模型

研究飛灰顆粒對(duì)麻面管束的磨損情況，本文選擇Edwards等提出的磨損計(jì)算模型，并采用隨機(jī)軌道模型[14]來(lái)考慮煙氣湍流脈動(dòng)對(duì)顆粒運(yùn)動(dòng)軌跡的影響，磨損方程為

(9)

式中：Rr為磨損速率，kg/(m2·s)；mp為顆粒流的質(zhì)量流率，kg/s；C(dp)為顆粒粒徑函數(shù)；α為顆粒與管壁的碰撞角度；b(up)為顆粒相對(duì)速度的函數(shù)；Af為顆粒撞擊壁面的面積，m2；n為碰撞面積Af的顆粒數(shù)目。C、f和b依據(jù)壁面邊界條件定義，與材料物性無(wú)關(guān)。

顆粒沉積率計(jì)算表達(dá)式為：

(10)

式中：Ndeposition為沉積的飛灰顆粒數(shù)；N為入射的顆?？倲?shù)。

3 計(jì)算結(jié)果與分析

首先，針對(duì)麻面管束模型進(jìn)行了數(shù)值模擬，定性分析了氣固兩相流外掠麻面管束的流動(dòng)特性及麻面管表面的磨損情況。然后定量分析了煙氣物性參數(shù)及管束結(jié)構(gòu)參數(shù)對(duì)積灰特性的影響。模擬所得結(jié)果與一些模擬及實(shí)驗(yàn)研究[6-8]做了相關(guān)對(duì)比，驗(yàn)證了結(jié)果的可靠性。

3.1 流動(dòng)與磨損特性分析

在分析含灰煙氣橫掠麻面管管束過(guò)程中顆粒的運(yùn)動(dòng)軌跡時(shí)，進(jìn)口煙氣速度選取8 m/s,飛灰顆粒平均粒徑為25 μm，并且模擬了相同條件下煙氣橫掠光管管束的情況，以便進(jìn)行對(duì)比分析。

圖3～4分別為煙氣繞流光管管束和麻面管管束的顆粒運(yùn)動(dòng)軌跡線，從圖中可以看到，粒徑較小的顆粒受到煙氣渦流作用明顯，在管束間呈波紋式交替運(yùn)動(dòng)，這種運(yùn)動(dòng)形式在麻面管束中體現(xiàn)更為明顯。另外，煙氣橫掠麻面管管束中，小粒徑顆粒的運(yùn)動(dòng)軌跡線前后變化不大，而光管管束中，小粒徑顆粒的軌跡線由密變疏，部分軌跡線消失，這說(shuō)明有一部分顆粒在該處沉積，通過(guò)兩種管束的對(duì)比，很好地說(shuō)明了麻面管管束在減輕積灰方面具有優(yōu)勢(shì)。

圖3 煙氣繞流光管管束的各粒徑顆粒運(yùn)動(dòng)軌跡

圖4 煙氣繞流麻面管管束的各粒徑顆粒運(yùn)動(dòng)軌跡

圖5表示的是相同條件下兩種管束的磨損情況。圖中可以看出，光管磨損情況比麻面管的嚴(yán)重，而且光管最大磨損發(fā)生的位置也比麻面管的靠前。這是因?yàn)槁槊婀鼙砻娴目咏Y(jié)構(gòu)，對(duì)于煙氣流動(dòng)有很好的導(dǎo)流作用，緩解了部分顆粒對(duì)管束壁面的沖擊，從而有效減輕磨損。

圖5 磨損分布云圖

3.2 煙氣物性參數(shù)對(duì)沉積率的影響

分析煙氣物性參數(shù)對(duì)顆粒沉積的影響時(shí)，主要從煙氣速度和顆粒粒徑兩個(gè)方面考慮，模擬分析了不同進(jìn)口煙速、不同粒徑下，麻面管外表面顆粒的沉積情況。

由圖6可看出，對(duì)于三種不同進(jìn)口煙速繞流麻面管管束，由v=6 m/s增加到v=8 m/s時(shí)，煙氣繞流麻面管過(guò)程中，管束前后的低速區(qū)域明顯減小，形成的渦結(jié)構(gòu)也相對(duì)較小，這有利于降低飛灰顆粒在管束間沉積的可能性，減輕積灰；而當(dāng)v=10 m/s時(shí)，圖中可以看到前排管束間的低速區(qū)域增大的較多，這對(duì)于減輕積灰不利。

圖6 不同進(jìn)口煙速的流場(chǎng)速度云圖

圖7 粒徑大小對(duì)沉積率的影響

3.3 粒徑對(duì)磨損特性影響

論文針對(duì)不同平均粒徑的顆粒對(duì)于磨損的影響情況進(jìn)行了模擬分析，圖8表示的是不同煙速下顆粒平均粒徑的大小對(duì)于磨損的影響曲線。由圖8可知，同一煙速下，隨著顆粒平均粒徑的增大，磨損也隨之增大。當(dāng)粒徑為15 μm時(shí)，顆粒主要受到煙氣的粘性作用，慣性小，隨著煙氣一起繞管外壁流動(dòng)，與壁面的碰撞幾率很小，磨損較??；當(dāng)粒徑為100 μm時(shí)，顆粒運(yùn)動(dòng)主要受到自身的慣性作用，受到煙氣流體的影響較小，在慣性作用下，大粒徑顆粒很容易穿過(guò)近壁區(qū)的邊界層與壁面發(fā)生相撞，導(dǎo)致碰撞次數(shù)增加，磨損加劇。

另外，對(duì)于較大煙速情況下，隨著顆粒平均粒徑的增大，磨損增加幅度較大。煙速為6 m/s下，粒徑為15 μm的顆粒造成的磨損量為1.66×10-8kg/(m2·s)，粒徑為100 μm時(shí)，磨損量為8.17×10-8kg/(m2·s)，磨損增加6.51×10-8kg/(m2·s)；煙速為10 m/s下，粒徑為15 μm，磨損量為7.75×10-8kg/(m2·s)，粒徑為100 μm時(shí)，磨損量為47.9×10-8kg/(m2·s)，磨損增加40.15×10-8kg/(m2·s)，后者增大的磨損量約是前者的6倍。

3.4 橫縱管距對(duì)沉積特性影響

為了分析麻面管管束間距、麻面管表面坑結(jié)構(gòu)參數(shù)對(duì)顆粒沉積率的影響情況，針對(duì)不同橫/縱管距、不同凹坑深度下，分別進(jìn)行了煙氣橫掠管束的模擬研究，該模擬過(guò)程中，煙速為8 m/s,飛灰顆粒平均粒徑為25 μm。

圖8 平均粒徑對(duì)磨損的影響

圖9和圖10分別為橫縱向管距對(duì)顆粒沉積率的影響關(guān)系曲線。可以看出，隨著橫向管距的增大，顆粒沉積率降低的幅度較大，這是因?yàn)殡S著橫向管距的增大，相鄰管道間煙氣的流通面積也增大，更多的飛灰顆粒從管間流過(guò)，使得管束間背風(fēng)側(cè)內(nèi)的顆粒沉積減少。隨著縱向管距的增大，顆粒沉積率整體呈現(xiàn)先降低再增大趨勢(shì)，s2/D=1.5～1.8時(shí)，沉積率降低的幅度很小，s2/D=1.8～2.0時(shí)，沉積率增大趨勢(shì)逐漸明晰，以s2/D=1.8時(shí)沉積率最低，從兩個(gè)圖比較可以得出，在進(jìn)行的模擬所有情況中，當(dāng)s1/D=3.0，s2/D=1.8時(shí)，顆粒沉積率最低，為20.72%。

圖9 橫向管距對(duì)沉積率的影響

究其原因，增大縱向管距，煙氣內(nèi)顆粒在流經(jīng)前排管子后有充足的空間繼續(xù)擴(kuò)展運(yùn)動(dòng)，降低了沉積的可能性，然而隨著縱向管距的繼續(xù)增大，顆粒沉積率有所上增，原因是前后空間距離的加大，煙氣繞流前排管束后，在管束間形成的渦流運(yùn)動(dòng)，不足以擾動(dòng)管束間的整個(gè)空間，使得管束間的軸向或徑向運(yùn)動(dòng)的煙氣對(duì)麻面管表面的沖刷能力減弱，導(dǎo)致顆粒沉積增多。

圖10 縱向管距對(duì)沉積率的影響

3.5 麻面坑深度對(duì)沉積特性影響

圖12表示s1/d=2.5，s2/d=1.8時(shí)坑深對(duì)麻面管管束磨損的關(guān)系曲線。從圖中可以看出，對(duì)于同一平均粒徑的飛灰顆粒，其對(duì)麻面管管束造成的磨損量隨坑深的增加而減小，而且對(duì)于d=100 μm的大粒徑顆粒，隨著坑深的增加，磨損量降低的幅度較大，效果更明顯。這是因?yàn)殡S著坑深的增加，煙氣流經(jīng)管束的流通面積增大，煙氣流動(dòng)速度相對(duì)減小，另外，坑深增加，管束表面的坑結(jié)構(gòu)對(duì)顆粒的導(dǎo)流作用增強(qiáng)，使顆粒產(chǎn)生“避中”效應(yīng)，不易靠近壁面，從而降低顆粒對(duì)管壁的碰撞沖擊幾率，使磨損減輕。

圖11 麻面管表面坑深度對(duì)沉積率的影響

圖12 坑深度對(duì)磨損的影響

4 結(jié)論

(1)煙氣橫掠麻面管管束時(shí)，麻面管表面的坑結(jié)構(gòu)對(duì)于煙氣具有很好的導(dǎo)流作用，麻面管表面顆粒沉積率和磨損量均比光管的低。

(2)同一煙速下，隨著顆粒平均粒徑的增大，磨損也隨之增大；較大煙速情況下，隨著顆粒平均粒徑的增大，磨損增加幅度較大，如平均粒徑為15 μm、100 μm的顆粒分別在煙速為6 m/s、10 m/s下造成的磨損量差值，后者約是前者的6倍。

(3)增大橫向管距，相鄰管道間煙氣通流面積增大，更多的飛灰顆粒從中間流過(guò)，積灰減少；縱向管距增大，顆粒沉積率整體呈現(xiàn)先降低再增大趨勢(shì)，在進(jìn)行所有情況分析中，確定了當(dāng)s1/D=3.0，s2/D=1.8時(shí)，顆粒沉積率最低，為20.72%。

[1]王宇，陳偉雄，于飛，等. 氣固兩相流外掠H型翅片管度損特性數(shù)值模擬[J]. 工程熱物理學(xué)報(bào)，2017，38(1)：129-133.

[2]樊泉桂,閻維平，閆順林,等.鍋爐原理[M].北京：中國(guó)電力出版，2011:178-180.

[3]李鈞, 閻維平, 邢德山,等. 某電廠670t/h鍋爐省煤器磨損原因分析及改造設(shè)計(jì)[J]. 熱力發(fā)電, 2007, 36(4):100-101.

[4]岑可法,樊建人,池作和,等.鍋爐和熱交換器的積灰、結(jié)渣、磨損和腐蝕的防止原理與計(jì)算[M].北京: 科學(xué)出版社,1994:4-5.

[5]董自翔，何雅玲，李冬，等.應(yīng)用格子Boltzmann方法研究顆粒橫掠管束的沉積過(guò)程[J].工程熱物理學(xué)報(bào)，2014,35(9):1802-1806.

[6]HAN H, HE Y L, TAO W Q, et al. A parameter study of tube bundle heat exchangers for fouling rate reduction[J]. International Journal of Heat & Mass Transfer, 2014, 72(72):210-221.

[7]鐘磊, 趙伶玲, 儲(chǔ)墨. 鰭片管磨損特性的數(shù)值研究[J]. 動(dòng)力工程學(xué)報(bào), 2016,36(5):389-394.

[8]王迎慧, 孫寧, 歸柯庭. 煙氣橫掠螺旋槽管束磨損特性的數(shù)值模擬[J].東南大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,44(3):585-590.

[9]李飛,史月濤,孫奉仲. 橢圓H型翅片管管束阻力特性的試驗(yàn)研究[J].中國(guó)電機(jī)工程學(xué)報(bào),2013,(32):34-39.

[10]王則力,羅坤,樊建人.顆粒與順列管束磨損的數(shù)值模擬[J].工程熱物理學(xué)報(bào),2008,29(9):1518-1520.

[11]陳述國(guó).省煤器磨損影響因素分析與數(shù)值模擬[J].電力學(xué)報(bào),2012,27(1): 90-93.

[12]PATWARDHAN A W.CFD modeling ofjet mixed tanks[J]. Chemical Engineering Science,2002,57(8):1307-1318.

[13]Moser R D, Kim J, Mansour N N. Direct numerical simulation of turbulet channel flow up to re[J]. Physics of Fluids,1999,11(7):943-945.

[14]歐益宏,杜揚(yáng),肖杰,等.粒子在柱形旋流器中的隨機(jī)軌道數(shù)值模擬[J].石油礦場(chǎng)機(jī)械,2006,35(1):27-29.

Numerical Analysis on Particle Deposition Characteristics of Flue Gas Flowing Across Pitted Tubes in Aligned Arrangement

YAN Shunlin, WANG Haoxuan

(School of Energy Power and Mechanical Engineering，North China Electric Power University，Baoding 071003，China)

In order to reduce the erosion and relieve the ash accumulation of heat exchange tube in boiler onomichi of power plant, the pitted tube, a new type of tube, is put forward and established in this paper. By choosing a suitable wear calculation and a particle deposition model, this paper carries on the numerical simulation of gas-solid flow on the basis of flue gas flowing cross pitted tubes in aligned arrangement; and the effects of the smoke’s physical parameters and tubes’s structural parameters on the pitted tubes’s particle deposition characteristics are analyzed, and the rough and smooth tube’s deposition under the same conditions are compared. The results show that the particle deposition rate and the wear of pitted tube’s surface are lower than those of the light tube under equivalent conditions. With the increase of the average particle size, the wear rate increases greatly under the condition of large smoke velocity; whens1/D=3.0,s2/D=1.8, the particle deposition rate is the lowest, and it is 20.72%. For the average small particle size, pitted tube withe=0.75 mm is better than other kinds when it comes to the effect of relieving the ash accumulation. In addition, the greater the pipe surface pits’ depth is, the lighter the wear is.

pitted tube; gas-solid flow; particle deposition; wear; numerical simulation

2017-08-28。

10.3969/j.ISSN.1672-0792.2017.11.012

TK 223.3

A

1672-0792(2017)11-0066-07

閆順林(1959-)，男，教授，研究方向?yàn)闊崃ο到y(tǒng)節(jié)能理論及應(yīng)用。