沈 聰,董振洲,王佳音,楊景軒,郝曉剛
(太原理工大學(xué) 化學(xué)化工學(xué)院,太原 030024)
高效率和低壓降一直是旋風(fēng)分離技術(shù)開(kāi)發(fā)追求的目標(biāo)。隨著持續(xù)不斷的改進(jìn),現(xiàn)代高效旋風(fēng)分離器已能基本除凈粒徑大于10 μm的顆粒,但對(duì)于5 μm以下的超細(xì)顆粒的捕集能力仍需進(jìn)一步挖掘。目前,國(guó)內(nèi)外學(xué)者通過(guò)實(shí)驗(yàn)和數(shù)值模擬的方法對(duì)旋風(fēng)分離器內(nèi)顆粒運(yùn)動(dòng)進(jìn)行研究,發(fā)現(xiàn)排塵口顆粒返混夾帶現(xiàn)象是超細(xì)顆粒逃逸的主要因素之一。因此,深入地認(rèn)識(shí)顆粒的返混夾帶現(xiàn)象對(duì)旋風(fēng)分離器結(jié)構(gòu)參數(shù)的設(shè)計(jì)優(yōu)化具有十分重要的意義。
黃學(xué)東[1]通過(guò)實(shí)驗(yàn)方法研究顆粒返混夾帶問(wèn)題,發(fā)現(xiàn)顆粒返混夾帶主要影響內(nèi)旋流區(qū)域的顆粒濃度,內(nèi)旋流對(duì)返混顆粒具有一定的二次分離作用。韓恒標(biāo)等[2]通過(guò)實(shí)驗(yàn)方法研究了PSC型旋風(fēng)管排塵錐內(nèi)顆粒濃度分布,發(fā)現(xiàn)開(kāi)縫的排塵結(jié)構(gòu)具有二次分離作用。萬(wàn)古軍等[3]、薛小虎等[4]、宋健斐等[5]從數(shù)值計(jì)算的角度去研究顆粒返混問(wèn)題,發(fā)現(xiàn)排塵口附近存在明顯的顆粒返混。孫濤[6]通過(guò)實(shí)驗(yàn)回歸分析得到灰斗返混量的計(jì)算公式,計(jì)算得到的總效率與實(shí)驗(yàn)得到的總效率誤差在0.79%以內(nèi)。上述研究揭示了影響顆粒返混的因素,但并未說(shuō)明顆粒返混的形成原因。祝華騰等[7]對(duì)不同結(jié)構(gòu)的旋風(fēng)分離器二次渦的數(shù)值模擬和分析發(fā)現(xiàn)排塵口處的局部旋渦流動(dòng)會(huì)對(duì)顆粒返混產(chǎn)生影響。元少昀等[8]、吳小林等[9]通過(guò)實(shí)驗(yàn)研究發(fā)現(xiàn)旋進(jìn)渦核現(xiàn)象能夠引起這種局部流動(dòng)。旋風(fēng)分離器全空間內(nèi)都存在旋進(jìn)渦核現(xiàn)象,但渦核的運(yùn)動(dòng)頻率會(huì)隨軸向位置而變化,并在排塵口達(dá)到最大[9-11]。因而,旋進(jìn)渦核會(huì)將部分已被壁面捕集的顆粒重新帶入內(nèi)旋流,造成大量顆粒的返混逃逸[12]。KOSAKI et al[13]通過(guò)實(shí)驗(yàn)及CFD模擬的方法發(fā)現(xiàn)在分離器內(nèi)部加入穩(wěn)渦桿或在排塵口下方加入圓錐能夠穩(wěn)定旋風(fēng)分離器內(nèi)部流場(chǎng),從而提高分離效率。吳小林等[14-15]、YOSHIDA et al[16-18]通過(guò)對(duì)相似內(nèi)構(gòu)件的實(shí)驗(yàn)與模擬研究,也得到類似觀點(diǎn)。雖然眾多學(xué)者發(fā)現(xiàn)了旋進(jìn)渦核的旋轉(zhuǎn)特性與顆粒返混存在聯(lián)系,但是至今還鮮少有對(duì)于兩者之間定量關(guān)系的研究。因此,深入研究渦核的旋轉(zhuǎn)特性與顆粒返混的定量關(guān)系成為本文重點(diǎn)。
為了研究渦核的旋轉(zhuǎn)特性與顆粒返混的定量關(guān)系需要排除旋風(fēng)分離器的基本結(jié)構(gòu)參數(shù)對(duì)渦核旋轉(zhuǎn)特性的影響。通過(guò)KOSAKI et al[13]的研究發(fā)現(xiàn),內(nèi)構(gòu)件能夠有效地改變渦核的旋轉(zhuǎn)特性。所以筆者選用圓錐作為穩(wěn)渦內(nèi)構(gòu)件,通過(guò)對(duì)比在不同入口氣速下有無(wú)圓錐的兩個(gè)尺寸旋風(fēng)分離器的內(nèi)流場(chǎng),來(lái)揭示渦核的旋轉(zhuǎn)特性與顆粒返混的定量關(guān)系,并對(duì)分離效率和壓降進(jìn)行分析。
求解控制方程時(shí)壓力速度耦合項(xiàng)采用SIMPLEC(semi-implicit method for pressure linked equations consistent)算法,壓力梯度項(xiàng)采用PRESTO(pressure staggering option)方法進(jìn)行處理,各對(duì)流項(xiàng)采用QUICK(quadratic upwind interpolation of convective kinematics)差分格式。由于旋進(jìn)渦核現(xiàn)象是隨時(shí)間變化的,采用非穩(wěn)態(tài)模擬,時(shí)間步長(zhǎng)選為2×10-4s.
入口氣體為常溫下的空氣,入口邊界條件為速度入口,分別給出垂直進(jìn)入入口截面的時(shí)均速度值,速度大小分別為15,20,25 m/s,入口處設(shè)置相應(yīng)的水力直徑和湍流強(qiáng)度。出口邊界條件設(shè)置為充分發(fā)展出口(outflow邊界條件),為保證充分發(fā)展條件的成立,計(jì)算中將旋風(fēng)分離器出口段加長(zhǎng)。氣相流場(chǎng)在壁面采用無(wú)滑移邊界條件。
兩相流場(chǎng)模擬計(jì)算時(shí),選用密度為2 700 kg/m3,粒徑為1 μm的顆粒,在入口處的射流源采用面源,即顆粒由入口截面上的每一個(gè)網(wǎng)格的中心位置射入。顆粒的入口速度設(shè)定與氣相入口速度相同;入口質(zhì)量濃度取10 g/m3.
圖1展示了旋風(fēng)分離器的基本結(jié)構(gòu)型式。坐標(biāo)系原點(diǎn)在頂板圓截面圓心處,旋風(fēng)分離器的中心軸為Z軸,且方向向下為正,氣體進(jìn)入旋風(fēng)分離器的方向?yàn)閅軸。Type-A旋風(fēng)分離器是標(biāo)準(zhǔn)平頂蝸殼式旋風(fēng)分離器,其內(nèi)徑為186 mm,筒體長(zhǎng)180 mm,錐體長(zhǎng)560 mm,矩形入口高度為110 mm,入口寬度為45 mm,排氣管直徑為90 mm,排氣管下端與入口底端平齊。Type-B旋風(fēng)分離器的結(jié)構(gòu)參數(shù)與Type-A旋風(fēng)分離器相同,區(qū)別僅是在排塵口處增加了圓錐。采用六面體結(jié)構(gòu)化網(wǎng)格對(duì)計(jì)算區(qū)域進(jìn)行離散,為了方便網(wǎng)格劃分,將圓錐頂角去除,具體尺寸見(jiàn)圖1.
為了驗(yàn)證計(jì)算模型的準(zhǔn)確性,與高翠芝等[19]文中筒體直徑186 mm模型的試驗(yàn)結(jié)果進(jìn)行了對(duì)比,模擬模型與其試驗(yàn)采用的旋風(fēng)分離器模型完全相同,具體尺寸如上文所述,入口氣速為20 m/s,進(jìn)行模擬。對(duì)網(wǎng)格數(shù)目的無(wú)關(guān)性進(jìn)行驗(yàn)證,通過(guò)對(duì)比筒體Z=100 mm處的速度沿徑向的分布圖,結(jié)果如圖2,綜合考慮精度和計(jì)算負(fù)荷,選取了網(wǎng)格數(shù)目為17萬(wàn)的網(wǎng)格為最終的模型,計(jì)算的總網(wǎng)格數(shù)為170 284.
圖1 旋風(fēng)分離器幾何模型
圖2 網(wǎng)格無(wú)關(guān)性驗(yàn)證
高翠芝等[19]采用點(diǎn)動(dòng)態(tài)壓力信號(hào)采集系統(tǒng)測(cè)量分離器各測(cè)點(diǎn)的壓力,并利用快速傅立葉變換(Fast Fourier Transform,簡(jiǎn)稱FFT)將測(cè)量到的原始時(shí)域信號(hào)轉(zhuǎn)換成易于分析的頻域信號(hào)(信號(hào)的頻譜),發(fā)現(xiàn)在分離器筒體及錐體段,分離器外旋流區(qū)壓力信號(hào)沒(méi)有明顯的主頻,而內(nèi)旋流壓力信號(hào)具有150 Hz的波動(dòng)頻率,因此,判定PVC的頻率為150 Hz.
本文使用相近的手段分析渦核的運(yùn)動(dòng)頻率。首先在排塵口截面建立一個(gè)虛擬面并監(jiān)測(cè)它的面最小靜壓隨時(shí)間的變化情況,這里認(rèn)為面最小靜壓點(diǎn)即為渦核中心點(diǎn),然后將這種變化記錄下來(lái)。其次,對(duì)記錄下來(lái)的數(shù)據(jù)進(jìn)行快速傅立葉變換,得到頻譜關(guān)系(如圖3)。其中有一個(gè)主峰,對(duì)應(yīng)的頻率為152 Hz.模擬得到壓力頻譜特性與實(shí)驗(yàn)的基本相同,僅是峰值頻率略低于實(shí)驗(yàn)觀測(cè)。
此外,王璐等[20]提出利用監(jiān)測(cè)面最小壓力方法可以描述出排塵口橫截面處渦核旋轉(zhuǎn)一周的時(shí)間。本文利用王璐所提出的方法得到的渦核旋轉(zhuǎn)一周的時(shí)間為0.006 8 s,與之對(duì)應(yīng)的旋轉(zhuǎn)頻率為147 Hz,這與FFT得到的基本頻率152 Hz非常接近,見(jiàn)圖3.
圖3 模擬驗(yàn)證結(jié)果圖
本文通過(guò)對(duì)比在不同入口氣速下,排塵口處有無(wú)圓錐的旋風(fēng)分離器分離空間的切向速度來(lái)研究圓錐的加入對(duì)切向速度的影響。圖4為Type-A與Type-B旋風(fēng)分離器內(nèi)流場(chǎng)的切向速度在不同軸向位置的分布圖。由圖可知,在整個(gè)分離空間內(nèi),Type-A與Type-B的切向速度分布曲線沿徑向基本保持一致,為明顯的雙駝峰分布,各個(gè)位置處的速度值也基本相同。但是,在接近排塵口的地方,由于圓錐的約束,Type-B中切向速度分布的對(duì)稱性明顯優(yōu)于Type-A,見(jiàn)圖4(f).總體來(lái)說(shuō),圓錐對(duì)于分離空間中氣流的切向速度影響較小。
圖5為Type-A與Type-B在不同氣速下軸向速度分布圖,從中可以發(fā)現(xiàn)圓錐對(duì)于軸向速度影響主要存在于上行流中,特別是幾何中心附近上行流速減小的現(xiàn)象(以下簡(jiǎn)稱為滯流現(xiàn)象),而對(duì)下行氣流的軸向速度影響不明顯。觀察不同的軸向位置可以發(fā)現(xiàn)在筒體的幾何中心附近,Type-B的軸向速度比Type-A的大;但是觀察圖5(f)可以發(fā)現(xiàn),在靠近排塵口附近的區(qū)域,Type-A與Type-B的軸向速度分布有較大的差異,Type-A軸向速度的滯留層消失,而Type-B仍存在較明顯的滯留層,且呈現(xiàn)較好的軸對(duì)稱性。圓錐一方面減輕了分離空間上部的滯流程度,另一方面也擴(kuò)大了發(fā)生滯流現(xiàn)象的區(qū)域??傮w來(lái)看,圓錐并未明顯改善分離空間內(nèi)的軸向速度分布。
目前,國(guó)內(nèi)外學(xué)者還很少研究圓錐對(duì)短路流的影響。筆者通過(guò)分析對(duì)比Type-A與Type-B的下行流量的變化來(lái)分析研究圓錐對(duì)短路流的影響。圖6為不同氣速下Type-A與Type-B的下行流量圖,從圖中可以發(fā)現(xiàn),在軸向位置Z=110 mm到Z=135 mm之間,兩個(gè)分離器內(nèi)的下行流量都會(huì)突然減小,而這個(gè)位置正是排氣管下放氣流短路流出分離器的位置。圓錐的加入并不會(huì)改變短路流區(qū)的大小。統(tǒng)計(jì)這兩個(gè)位置的下行流量之差即得到短路流量。在三個(gè)入口氣速下,Type-A的短路流量分別是0.016 69 m3/s、0.021 06 m3/s和0.026 83 m3/s,Type-B的短路流量分別是0.015 74 m3/s、0.021 33 m3/s和0.026 22 m3/s.對(duì)應(yīng)入口氣速下,Type-B的短路流量與Type-A相差最多5.7%,變化很小。
圖4 Type-A與Type-B的切向速度分布圖
圖5 Type-A與Type-B的軸向速度分布圖
對(duì)旋風(fēng)分離器Type-A與Type-B進(jìn)行純氣相模擬,待其流場(chǎng)穩(wěn)定后,得到渦核擺動(dòng)的云圖。圖7為20 m/s下,Type-A與Type-B的純氣相條件下的靜壓分布圖??梢钥闯鯰ype-A中旋風(fēng)分離器的旋渦尾端已經(jīng)進(jìn)入到了料腿中部,并且產(chǎn)生了很明顯的擺尾現(xiàn)象;而Type-B的旋進(jìn)渦核由于受到了圓錐的阻礙作用,無(wú)法進(jìn)入料腿中,并且渦核尾端被牢牢固定到了圓錐上方,看不到明顯的渦核擺尾現(xiàn)象。
在純氣相流場(chǎng)的基礎(chǔ)上對(duì)兩種旋風(fēng)分離器內(nèi)的氣固兩相流動(dòng)進(jìn)行模擬。計(jì)算過(guò)程中隨時(shí)監(jiān)測(cè)排氣管出口的顆粒濃度,取濃度穩(wěn)定之后2 s內(nèi)的流場(chǎng)數(shù)據(jù)分析氣固兩相流場(chǎng)的非穩(wěn)態(tài)性及顆粒捕集情況。
圖6 Type-A與Type-B的下行流量圖
圖7 純氣相條件下的Y=0截面靜壓分布圖
圖8為Type-A與Type-B在氣固兩相下的靜壓云圖。對(duì)比圖7與圖8可以發(fā)現(xiàn),Type-A的兩幅靜壓云圖在插入管與分離空間差異較小,但是在排塵口與料腿處有較為明顯的差異。在純氣相條件下,Type-A中的旋進(jìn)渦核已經(jīng)進(jìn)入到料腿中,但是加入顆粒后可以明顯發(fā)現(xiàn)它的自然旋風(fēng)長(zhǎng)變短了。
圖8 氣固兩相共存穩(wěn)定的Y=0截面的壓力云圖
同時(shí),觀察Type-B可以明顯發(fā)現(xiàn),加入圓錐后通入顆?;静粫?huì)對(duì)自然旋風(fēng)長(zhǎng)造成影響,其旋進(jìn)渦核依舊是被牢牢地固定在圓錐的上方。圓錐對(duì)于自然旋風(fēng)長(zhǎng)的影響很大。
圖9 氣固兩相下排塵口截面渦核旋轉(zhuǎn)特性圖
上文所開(kāi)展的流場(chǎng)分析表明,圓錐的加入主要改變了旋進(jìn)渦核的運(yùn)動(dòng)特性,而對(duì)各個(gè)分離效率理論模型所關(guān)注的外旋流區(qū)域的影響很有限。因此,有圓錐和無(wú)圓錐兩種結(jié)構(gòu)的捕集效率和顆粒逃逸情況的差別主要?dú)w因于旋進(jìn)渦核特性的差別。由于短路流也會(huì)攜帶部分顆粒逃逸,因此,在旋風(fēng)分離器的排氣管正下方設(shè)置了如圖10所示的監(jiān)測(cè)面,分別統(tǒng)計(jì)短路流攜帶的顆粒(短路逃逸率)和上行氣流夾帶逃逸的顆粒(返混夾帶率)。結(jié)果如表1所示。
圖10 顆粒逃逸的類型:短路流逃逸和返混夾帶逃逸
表1 捕集效率、短路流逃逸率與返混夾帶率
從中可以看出在各個(gè)入口氣速下,Type-A的效率均遜于Type-B,加入圓錐后效率至少提高了20.9%。從顆粒逃逸情況看,圓錐的加入降低了短路逃逸的顆粒量,但影響非常有限,不足5%,而返混夾帶的顆粒量則降低了大約20%,是提升捕集效率的關(guān)鍵,與渦核在排塵口處的運(yùn)動(dòng)特性密不可分。分離空間所捕集的顆粒匯聚于排塵口,而此處又是整個(gè)分離空間直徑最小的部分,也是被捕集顆粒最靠近高速上行氣流的地方。渦核內(nèi)的壓力是低于邊壁處的,當(dāng)其靠近邊壁,在壓差作用下邊壁處的氣流匯入,部分已分離的小顆粒受氣流裹挾而進(jìn)入上行流。如圖9所示,Type-A中,渦核遠(yuǎn)離幾何中心而非??拷叡冢w粒受氣流裹挾而返混的幾率大;Type-B中,渦核基本被固定在了圓錐上方,偏心幅度被極大的削弱,氣流裹挾顆粒進(jìn)入上行流的幾率降低,另外,渦核旋轉(zhuǎn)頻率的降低意味著相同時(shí)間內(nèi)渦核與壁面的接觸減少,也有利于削弱顆粒因氣流夾帶而返混。
壓降是旋風(fēng)分離器的另一性能指標(biāo),小壓降意味著能量的節(jié)約。本文通過(guò)讀取入口截面中心點(diǎn)和出口截面中心點(diǎn)的靜壓值,以兩者的差作為旋風(fēng)分離器的壓降,結(jié)果如圖11.從圖中可以看出,隨著氣速的增大,壓強(qiáng)呈線性增大的趨勢(shì),但是Type-A和Type-B的壓降差值并不大,最大不超過(guò)200 Pa,不足壓降的0.5%.因此,圓錐結(jié)構(gòu)的加入在增加旋風(fēng)分離器效率的同時(shí)并不會(huì)消耗額外的能量。
圖11 兩種結(jié)構(gòu)的壓降對(duì)比
本文采用數(shù)值模擬方法研究了排塵口加入小圓錐對(duì)旋風(fēng)分離器內(nèi)流場(chǎng)與顆粒返混的影響,得到如下結(jié)論:
1) 排塵口處安放的圓錐可以改善其附近氣流切向和軸向速度沿徑向分布的對(duì)稱性,但對(duì)分離空間大部分區(qū)域內(nèi)氣流的分布影響很小。
2) 在旋風(fēng)分離器的排氣管下方存在的氣流短路逃逸現(xiàn)象不會(huì)因圓錐的加入而明顯變化,相應(yīng)的,顆粒隨短路氣流逃逸情況并沒(méi)有明顯改善。
3) 圓錐能夠很好地抑制渦核偏離分離器幾何中心的現(xiàn)象,降低已分離顆粒在氣流裹挾下進(jìn)入高速上行流的幾率。相比未安裝圓錐的結(jié)構(gòu),分離效率提升20%以上,同時(shí)壓降沒(méi)有明顯增加。