梁容真,田 偉,閻富生

(東北大學(xué)冶金學(xué)院,遼寧 沈陽110819)

輕燒鎂用途十分廣泛,作為提煉金屬鎂的基礎(chǔ)原料,主要用作耐火材料、化工原料和建材原料等[1]。過去反射窯、豎窯、回轉(zhuǎn)窯以及懸浮窯等爐窯是煅燒菱鎂礦制取輕燒鎂的主要設(shè)備[2]。然而,在低碳、低能耗、綠色環(huán)保的理念下以及對產(chǎn)品質(zhì)量要求的不斷提高,傳統(tǒng)爐窯均存在煅燒溫度不均勻、煅燒溫度與時(shí)間不易控制和能耗高等問題,不能滿足現(xiàn)階段輕燒鎂行業(yè)的高標(biāo)準(zhǔn)、高要求。旋流動(dòng)態(tài)煅燒系統(tǒng)成為目前國內(nèi)外比較先進(jìn)的煅燒工藝,尤其是制取高活性、高純度的輕燒鎂粉[3]。但是,在實(shí)際生產(chǎn)過程中,旋流動(dòng)態(tài)煅燒工藝存在動(dòng)力消耗太大、產(chǎn)品回收率不穩(wěn)定等問題。因此,為了保證整個(gè)輕燒鎂旋流動(dòng)態(tài)煅燒系統(tǒng)產(chǎn)品活性更高、系統(tǒng)更節(jié)能以及減少產(chǎn)品資源浪費(fèi),降低該系統(tǒng)的動(dòng)力消耗,對該工藝重要組成部分—旋風(fēng)分離器進(jìn)行研究,降低其壓力損失以及提高其分離效率成為該工藝的重中之重[4-5]。

蝸殼式旋風(fēng)分離器內(nèi)部是非常復(fù)雜的三維非對稱強(qiáng)烈旋轉(zhuǎn)湍流流動(dòng),采取實(shí)驗(yàn)的方法較難測量。雖然現(xiàn)在可以采用激光多普勒測速儀和PIV粒子成像測速儀等先進(jìn)測量儀器來測定旋風(fēng)分離器內(nèi)部復(fù)雜的氣流旋轉(zhuǎn)流動(dòng)特性,但是由于儀器本身的性能缺陷和較高成本的限制,實(shí)驗(yàn)測量仍然不是最理想的手段。同時(shí)蝸殼式旋風(fēng)分離器內(nèi)部是十分復(fù)雜的兩相流運(yùn)動(dòng),具體體現(xiàn)為固體顆粒在旋轉(zhuǎn)氣流中的分離過程,而固體顆粒的運(yùn)動(dòng)主要依賴于氣相的流動(dòng)。因此研究旋風(fēng)分離器內(nèi)氣固的運(yùn)動(dòng)規(guī)律對于提高分離器性能、優(yōu)化分離器的結(jié)構(gòu)具有指導(dǎo)意義。

基于Fluent軟件建立蝸殼式旋風(fēng)分離器物理模型,采用數(shù)值模擬的方法對蝸殼式旋風(fēng)分離器內(nèi)部流場進(jìn)行研究,分別研究其入口速度、入口顆粒粒徑和入口顆粒濃度對分離器內(nèi)分離效率和壓力損失的影響,從而為蝸殼式旋風(fēng)分離器內(nèi)氣固兩相研究提供理論支撐。

1 建立模型

1.1 分離器結(jié)構(gòu)和網(wǎng)格示意圖

選用THCX長錐體旋風(fēng)分離器,其結(jié)構(gòu)尺寸如圖1所示,其中具體尺寸參數(shù)如下:a=353 mm、b=94 mm、c=339 mm、d=250 mm、D=500 mm、e=200 mm、g=150 mm、h=359 mm、H=1759 mm 和r=47 mm。根據(jù)蝸殼式旋風(fēng)分離器結(jié)構(gòu),將其分為(180°蝸殼)入口段、圓柱段、錐體段、排料口和排氣管五個(gè)部分。采用Icem軟件對蝸殼式旋風(fēng)分離器進(jìn)行六面體網(wǎng)格劃分,生成三維計(jì)算域網(wǎng)格如圖2所示。經(jīng)過網(wǎng)格獨(dú)立性檢驗(yàn),最終確定網(wǎng)格數(shù)為23萬。

圖1 蝸殼式旋風(fēng)分離器結(jié)構(gòu)圖

圖2 三維計(jì)算域網(wǎng)格

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

蝸殼式旋風(fēng)分離器氣相模型包含質(zhì)量、動(dòng)量守恒方程,控制方程為三維雷諾時(shí)均Navier-Stokes方程。

關(guān)于湍流模型的選擇,與各種k-ε模型相比,RSM模型完全拋棄了基于各向同性渦粘性的Boussinesq假設(shè),更加嚴(yán)格地考慮了流線型彎曲、旋渦、旋轉(zhuǎn)和張力快速變化,對于復(fù)雜流動(dòng)有更高的精度預(yù)測的潛力,尤其是對于旋風(fēng)分離器內(nèi)流場真實(shí)情況的預(yù)測。因此,采用RSM模型[6]。

1.3 邊界和初始條件

研究蝸殼式旋風(fēng)分離器內(nèi)氣固兩相流場。排氣管出口邊界設(shè)置為自由出流,流量權(quán)重為1;排料口邊界設(shè)置為自由出流,流量權(quán)重為0;壁面邊界條件:無滑移壁面邊界條件,壁面粗糙高度、壁面粗糙度分別設(shè)置為0和0.5,并且采用標(biāo)準(zhǔn)壁面函數(shù)法來處理邊界湍流。離散方程組采用求解壓力—速度耦合方程的半隱方法—SIMPLE求解,擴(kuò)散項(xiàng)離散采用中心差分,逐行迭代[7]。為提高求解精度,各控制方程中對流項(xiàng)的離散采用三階精度的QUICK格式。

2 模擬結(jié)果與討論

2.1 入口速度的影響

入口速度是操作參數(shù)中非常重要的一項(xiàng),它對旋風(fēng)分離器的分離性能有著很大影響。為了研究入口速度對蝸殼式旋風(fēng)分離器性能的影響,通過數(shù)值計(jì)算的方法,分別計(jì)算了入口速度為10 m/s、14 m/s、18 m/s、22 m/s、26 m/s的壓降和不同顆粒粒徑的分離效率。

2.1.1 入口速度對壓降的影響

目前,壓降的理論推導(dǎo)還沒有非常精確的方法,一般都根據(jù)實(shí)驗(yàn)數(shù)據(jù)來總結(jié)經(jīng)驗(yàn)公式,其通用的經(jīng)驗(yàn)公式可表示為:

式中:Vi為入口速度,m/s;ρg為氣體密度,kg/m3;ζi阻力系數(shù)。

旋風(fēng)分離器中的氣流是三維強(qiáng)旋轉(zhuǎn)湍流,雷諾數(shù)很高,因此,旋風(fēng)分離器的阻力系數(shù)與操作參數(shù)無關(guān),只與結(jié)構(gòu)參數(shù)相關(guān)。一般可以表示為:

式中:X,Z為旋風(fēng)分離器結(jié)構(gòu)參數(shù)的函數(shù);A0為旋風(fēng)分離器進(jìn)口截面積,m2;Dx為旋風(fēng)分離器的排氣管直徑,mm。

關(guān)于旋風(fēng)分離器壓降的經(jīng)驗(yàn)公式,其中被普遍認(rèn)可的主要有Coker公式,Cacal&Martinez-Benet公式,Shepherd&Lapple公式以及Dirgo公式[8]。這些公式都是由試驗(yàn)等方法歸納的旋風(fēng)分離器壓降的計(jì)算模型,它們往往都是些關(guān)于旋風(fēng)分離器結(jié)構(gòu)參數(shù)以及旋風(fēng)分離器內(nèi)氣固兩相流動(dòng)的數(shù)值計(jì)算總結(jié)的經(jīng)驗(yàn)、半經(jīng)驗(yàn)公式[9]。其中,這些相關(guān)公式可以表示為:

Shepherd&Lapple公式:

式(1)-式(6)中:a,b分別為旋風(fēng)分離器進(jìn)口截面的高度和寬度,mm;Dx為排氣管直徑,mm;S為排氣管插入深度,mm;H為旋風(fēng)分離器總高度,mm;h為分離器圓柱段高度,mm;D為分離器筒體外徑,mm;Dd為排料口直徑,mm。

當(dāng)顆粒濃度為30 g/m3時(shí),通過改變邊界條件下的速度邊界條件,進(jìn)行數(shù)值模擬,計(jì)算出不同入口速度下蝸殼式旋風(fēng)分離器的壓降值和阻力系數(shù),如表1所示。為了便于分析研究,將表中模擬計(jì)算數(shù)據(jù)與經(jīng)驗(yàn)公式理論計(jì)算數(shù)據(jù)繪制成曲線圖,如圖3所示。

表1 不同入口速度下的壓降數(shù)值表Table 1 Statistics of pressure loss at different inlet speeds

圖3 入口速度和壓降的關(guān)系Fig.3 Relationship between inlet speed and pressure loss

由圖3可以看出,蝸殼式旋風(fēng)分離器的壓降隨著入口速度的增大而增大,基本上呈線性變化。與經(jīng)驗(yàn)公式計(jì)算結(jié)果相比,總體變化趨勢基本一致,只是壓降增加的幅度略有差別,這是因?yàn)樾L(fēng)分離器結(jié)構(gòu)不同導(dǎo)致的。其中,Dirgo公式的計(jì)算結(jié)果與本文的模擬結(jié)果最為接近。從能量變化來看,能量的損失會(huì)因分離器入口速度的增大而增大,這是由于旋風(fēng)分離器的磨損與氣體流速的四次方成正比[10],所以入口速度太大會(huì)增加旋風(fēng)分離器的壓力損失。因此,在旋風(fēng)分離器的分離效率和動(dòng)力損耗之間要有一個(gè)衡量指標(biāo),在保證所需分離效率的基礎(chǔ)上最大化的降低入口速度,降低能耗。

2.1.2 入口速度對分離效率的影響

當(dāng)顆粒濃度為10 g/m3,通過改變邊界條件下的速度邊界條件,進(jìn)行數(shù)值模擬來計(jì)算不同入口速度下的分離效率,將其數(shù)值計(jì)算結(jié)果繪成表2。

由表2所示,總體來看,旋風(fēng)分離器的分離效率隨著入口速度的增大而增大,并且可以發(fā)現(xiàn),入口速度的變化對分離效率的影響比較大。究其原因,切向速度決定著離心力的大小,從而對分離效率的影響起著決定性的作用。然而,當(dāng)顆粒粒徑為1μm時(shí),分離效率隨著入口風(fēng)速先增大再減小,入口風(fēng)速為26 m/s時(shí)的分離效率要略小于速度為22 m/s時(shí)的分離效率。這是因?yàn)?小顆粒受湍流以及顆粒碰撞反射等因素影響,很容易使已經(jīng)密集在內(nèi)壁上的顆粒重新卷吸起來;并且,分離器內(nèi)的向心徑向速度和上行軸向速度會(huì)隨著入口氣速的增大而增大,這使得小顆粒很容易被帶出排氣管而逃逸,使得分離效率大大降低;另外,還有許多二次渦流的因素也會(huì)降低小顆粒的分離效率。

表2 不同入口速度下的分離效率值Table 2 Separation efficiency at different inlet speeds

2.2 入口顆粒粒徑的影響

通常所說的分離效率是指旋風(fēng)分離器的總效率。然而,對于不同物性的顆粒,不同粒徑的顆粒,以及不同生產(chǎn)用途下分離器的總效率不能充分地反映分離器的分離性能。因此,分級分離效率,簡稱分級效率,即分離器對于不同顆粒粒徑的分離能力,它表示的是旋風(fēng)分離器對于特定粒徑顆粒的分離效率,與總效率相比更能充分地反映旋風(fēng)分離器的分離性能[9]。下面,主要研究不同顆粒粒徑下的分離效率變化。

2.2.1 入口顆粒粒徑對分離效率的影響

當(dāng)顆粒濃度為30 g/m3時(shí),以顆粒粒徑為變量,進(jìn)行數(shù)值模擬,計(jì)算出入口速度分別為18 m/s和22 m/s下蝸殼式旋風(fēng)分離器的分離效率與顆粒粒徑的關(guān)系,如表3所示。

表3 不同顆粒粒徑下的分離效率值Table 3 Separation efficiency under different size of particle

由表3可以看出,相同入口速度下,分離效率隨著顆粒粒徑的增大而增大,但增加的幅度有所降低,最終逐漸趨向穩(wěn)定,并且當(dāng)顆粒粒徑大于或等于11μm時(shí),分離器的分離效率均能達(dá)到100%。主要原因是粒徑較大的顆粒進(jìn)入分離器后受到的離心力較大,相比于小顆粒而言大顆粒在分離器內(nèi)隨氣流轉(zhuǎn)動(dòng)的圈數(shù)要少,所以能使大顆粒較早地碰到分離器器壁,從而較快地進(jìn)入排料口被捕集。小顆粒的分離效率低的主要原因是小顆粒所受的離心力小,其次小顆粒很容易受到徑向速度較大而引起的短路流的影響,直接被氣流帶出排氣管而逃逸;并且由于小顆粒對氣流有很大的跟隨性,較大一部分顆粒會(huì)跟隨旋轉(zhuǎn)氣流一直做旋轉(zhuǎn)運(yùn)動(dòng),最終被捕集,也會(huì)有另一小部分顆粒會(huì)在旋風(fēng)分離器內(nèi)做無限循環(huán)運(yùn)動(dòng),“上灰環(huán)”現(xiàn)象就是這種情況,這種情況被認(rèn)為是旋風(fēng)分離器無法分離該顆粒。

理論上每種旋風(fēng)分離器都有一個(gè)確定的臨界粒徑,大于臨界粒徑的顆粒會(huì)被捕集,而小于臨界粒徑的顆粒不會(huì)被捕集。但實(shí)際上,顆粒在分離器內(nèi)會(huì)受到許多不確定性因素的影響,比如顆粒間的碰撞、團(tuán)聚、破碎以及靜電和粒子吸引等因素,這使得顆粒在分離器內(nèi)的運(yùn)動(dòng)有很大的不確定性,因此,一些小于臨界粒徑的顆粒也會(huì)被捕集,而一些大于臨界粒徑的顆粒也會(huì)逃逸。

2.3 入口顆粒濃度的影響

入口顆粒濃度是工業(yè)生產(chǎn)過程中非常重要的操作參數(shù),它決定著生產(chǎn)產(chǎn)量的大小。為了研究入口速度對蝸殼式旋風(fēng)分離器性能的影響,通過數(shù)值計(jì)算的方法,分別計(jì)算了入口顆粒濃度為1 g/m3、10 g/m3、20 g/m3、30 g/m3、40 g/m3的壓 降 和 分 離效率。

2.3.1 入口顆粒濃度對壓降的影響

當(dāng)入口速度為18 m/s時(shí),變化入口顆粒濃度分別為通過改變離散相下射流的入口顆粒濃度,進(jìn)行數(shù)值模擬,計(jì)算出不同入口顆粒濃度下蝸殼式旋風(fēng)分離器的壓降。如表4所示。

表4 不同入口顆粒濃度下的壓降數(shù)值表Table 4 Statistics of pressure loss at differentinlet particle concentration

在本文研究的低濃度入口下,蝸殼式旋風(fēng)分離器的壓降隨入口顆粒濃度的增大而減小,但是總壓降的數(shù)值減小的幅度逐漸變緩,這是因?yàn)樵诜蛛x器排氣管出口位置沒有安裝穩(wěn)流裝置或者回復(fù)裝置,導(dǎo)致該位置旋轉(zhuǎn)氣流的動(dòng)壓較大,使得總壓較大,因此總壓降就較小。

2.3.2 入口顆粒濃度對分離效率的影響

除了入口速度之外,入口氣流顆粒濃度對分離器分離效率也有較大影響。下面探索不同顆粒濃度與分離效率的關(guān)系,確定入口速度為18 m/s,變化入口顆粒濃度分別為1 g/m3、10 g/m3、20 g/m3、30 g/m3、40 g/m3,通過改變離散相下射流的入口顆粒濃度,進(jìn)行數(shù)值模擬,計(jì)算出不同入口顆粒濃度下蝸殼式旋風(fēng)分離器的分離效率。如表5所示。

表5 不同入口顆粒濃度下的分離效率值Table 5 Separation efficiency under different inlet particle concentration

如表5入口氣流顆粒濃度與分離效率的關(guān)系所示,旋風(fēng)分離器的總分離效率和粒級分離效率隨著氣流顆粒濃度的增大而增大;大顆粒增加的幅度較小,而小顆粒增加的幅度較大。而且,入口氣流濃度越大,小顆粒分離效率增加的越多,主要是濃度較高的氣流對小顆粒有著更加明顯的攜帶作用,所以更多的小顆粒得到捕集。然而,當(dāng)濃度增加到一定值時(shí),如表5所示,9μm的分離效率在濃度大于20 g/m3時(shí)和5μm的分離效率在濃度大于30 g/m3時(shí),它們各自的分離效率逐漸趨于穩(wěn)定;大顆粒趨于穩(wěn)定的濃度較于小顆粒而言要低,這是因?yàn)闅饬鲗Υ箢w粒的攜帶作用沒有對小顆粒的攜帶作用大。

雖然入口氣流顆粒濃度的增加會(huì)提高分離效率,但是,在旋風(fēng)分離器的實(shí)際工業(yè)應(yīng)用中,顆粒對分離器壁面的磨損也會(huì)隨著顆粒濃度的增加而加劇,使得旋風(fēng)分離器的維修周期和使用壽命變短,而顆粒也會(huì)因碰撞而粉碎變細(xì),更加不利于分離。因此,在大多數(shù)情況下,人們不會(huì)指望只經(jīng)過一個(gè)分離器的分離而達(dá)到分離目的,而是經(jīng)過幾次分離,逐漸降低小顆粒群的濃度和粒度,最終達(dá)到分離效率指標(biāo);同樣,在輕燒鎂旋流動(dòng)態(tài)煅燒系統(tǒng)中的回收系統(tǒng)和冷卻系統(tǒng)都是兩級回收,目的是為了達(dá)到更高的分離效率。

3 結(jié)語

通過相間耦合的隨機(jī)軌道模型在拉格朗日坐標(biāo)下模擬蝸殼式旋風(fēng)分離器內(nèi)的兩相流運(yùn)動(dòng)?？疾炝朔蛛x器入口操作參數(shù)對分離性能的影響。主要模擬計(jì)算結(jié)果可以總結(jié)出以下結(jié)論:

(1)蝸殼式旋風(fēng)分離器的壓降隨入口風(fēng)速的增大而增大,兩者基本上呈指數(shù)變化關(guān)系。增加蝸殼式旋風(fēng)分離器的入口風(fēng)速能夠提高分離效率,但是壓力損失也逐漸增大。

(2)相同入口速度下,分離效率隨著顆粒粒徑的增大而增大,但增加的幅度有所降低,最終逐漸趨向穩(wěn)定,并且當(dāng)顆粒粒徑大于或等于11μm時(shí),分離器的分離效率均能達(dá)到100%。

(3)在一定范圍內(nèi),適當(dāng)?shù)卦黾尤肟陬w粒濃度既能降低壓力損失又能提高分離效率。