楊 鑫,唐雪蓮,陳春林,張 千,劉向陽,趙天濤

(重慶理工大學 化學化工學院, 重慶 400054)

當前,我國養(yǎng)豬場規(guī)?；目焖侔l(fā)展為居民生活水平和國民生產總值提高做出了貢獻。然而,由于大多數養(yǎng)豬場沒有污水處理系統(tǒng),亂排放污染了環(huán)境。這些未經處理達標的污水會浸入地下水系統(tǒng)、流入河流,一定程度污染人們的生活用水[1-3],因此處理養(yǎng)豬場廢水使其達到國家排放標準成為一項急迫的任務。

近年來,養(yǎng)豬場廢水主要采用厭氧折流板反應器(anaerobic baffled reactor,ABR)[4],生物膜反應器(membrane bio-reactor,MBR)[5]等設備處理。生物轉盤(rotating biological contactor,RBC)因具有良好的凈化效果和低能耗特點,已成為處理養(yǎng)豬場廢水的一種主要反應器。生物轉盤是一種好氧生物膜處理裝置,被廣泛應用于有機物脫氮、脫磷。它由轉軸、圓盤、接觸反應槽和驅動裝置等組成[6-7]。接觸反應槽中裝有要處理的廢水,驅動裝置以一定的轉速帶動轉軸旋轉,嵌在轉軸上的盤片隨著轉軸轉動,盤片的一部分浸沒在接觸反應槽中,另一部分暴露在空氣中。轉盤在轉動過程中交替和空氣以及污水接觸,一段時間后,在轉盤盤片表面附著一層帶有大量微生物的生物膜,生物膜中的微生物吸附和降解廢水中的有機污染物,從而實現凈化污水的功能[8],因此生物膜在RBC處理廢水中起到最為關鍵的作用。因為生物膜的生長緩慢,采用傳統(tǒng)實驗方法對生物膜進行研究的整個實驗周期較長,同時流體的流動特性也會影響生物膜的生長以及新老更替,因此對生物轉盤流體的流動特性的研究很有必要。

隨著計算機技術的快速發(fā)展,計算流體動力學(computational fluid dynamics,CFD)技術在污水處理領域中應用越來越廣泛[9-10],因此可通過CFD模擬生物轉盤內部流場流動,研究流體的流動特性。該方法具有模擬仿真實驗周期短、成本低等優(yōu)點。

目前,國內外學者對生物轉盤的研究中,針對傳統(tǒng)生物轉盤啟動慢、掛膜性能差等缺點,馮迪等[11]、唐云鷺等[12]采用新型生物轉盤輔以生物接觸氧化技術處理有機污水,發(fā)現新型生物轉盤具有抗沖擊能力強、掛膜效率高、微生物量大等優(yōu)點,在最佳條件下處理城市生活污水時能夠達到污染物排放標準。 為了克服RBC堵塞和壓降問題,Ravi等[13]對RBC進行了改性處理,改性后的RBC無堵塞,無壓降,適用于有機化合物的處理。上述研究為改進生物轉盤結構提供了一定的指導意見,然而目前關于生物轉盤的氣液兩相流的研究較少。為此,鄧斌等[14]、喻江等[15]對盤面開窗類型的生物轉盤表面液膜更新和流動特性進行了數值模擬,發(fā)現圓盤表面液膜加速區(qū)液膜更新頻率最快,開窗形成的自由膜速率大,更新頻率快,不同厚度轉盤的液膜隨轉盤厚度的增加而增加。

劉向陽等[1]提出了一種新型生物膜反應器處理養(yǎng)豬場廢水。該反應器的主要裝置為轉盤,在旋轉盤片的面上開很多方孔,以異養(yǎng)硝化-好氧反硝化菌為生物強化劑,以聚丙烯腈基活性炭纖維作為填料的生物接觸氧化池(biological contact oxidation,BCO)為生物膜反應器,對養(yǎng)豬場廢水處理效果明顯,但未對反應器中的新型生物轉盤流場的流動特性進行CFD模擬和研究。為了解該新型生物膜反應器的流動特性,利用CFD數值模擬對其流場進行研究,重點探究所用生物轉盤的速度云圖和速度矢量圖、不同操作轉速下的速度圖和速度矢量圖,研究該生物轉盤的流動特性,為該新型生物轉盤操作參數優(yōu)化提供依據。

1 數值模擬

1.1 數值模擬物理模型

物理模型來源于文獻[1]采用的養(yǎng)豬場污水處理生物膜反應器實驗裝置。相對于傳統(tǒng)的實心生物轉盤盤片,該新型生物轉盤的盤片表面有很多方孔,如圖1所示。此生物轉盤由4個盤片構成1個生物轉盤組,生物轉盤的結構尺寸如表1所示,其三維物理模型如圖2所示。對該裝置在間歇操作狀態(tài)下的污水處理進行數值模擬分析。

圖1 生物轉盤盤片示意圖

表1 生物轉盤結構尺寸

圖2 生物轉盤三維物理模型結構示意圖

流場計算開始之前,對生物轉盤進行網格劃分以對計算區(qū)域進行離散。采用mesh進行非結構性網格劃分,如圖3所示。當模型網格個數分別為4 013 325、4 076 442、4 355 029,以轉速為6 r/min時為例,生物轉盤的網格劃分和速度的模擬效果如圖4所示,速度值隨位置的變化基本一致,得到的結果相差不到5%。因此,網格無關性驗證結果成立。

圖4 網格數對速度的影響曲線

1.2 數值模擬條件

基于盤片表面上連續(xù)更新的液膜,生物轉盤內部流場實際上是氣-液兩相傳質。因需要得到盤片自由表面流動的情況,所以采用Hirt等[16]提出的流體體積(volume of fluid,VOF)模型進行計算。VOF模型是固定在歐拉網格下的表面跟蹤方法,通過求解單個動量方程和處理通過該區(qū)域的每種流體的體積分數來模擬2種或3種不能混合的流體。在VOF模型中,不同的流體組分共用一套動量方程,通過引進相體積分數這一變量,實現對每個計算單元相界面的追蹤。跟蹤相之間的界面是通過求解一相或者多相的容積比率的連續(xù)方程來完成的。對第q相有:

(1)

其中:Sα為質量源項,在默認情形下方程(1)右端源項為零,但當給每一項指定常數或用戶定義的質量源,則右端不為零。

采用ANSYS FLUENT對新型生物轉盤進行數值模擬計算?；趬毫Φ乃矐B(tài)求解器,采用實驗室的養(yǎng)豬場污水和空氣作為模擬介質,空氣為連續(xù)相,污水為第二相(污水參數見表2),采用間歇操作,水位保持在7 cm深度,即轉盤盤片浸沒在水中一半。

表2 空氣和污水的物性參數

計算流場采用的是多重參考系法(multiple reference frame,MRF),將計算區(qū)域劃分為圓盤區(qū)域和其他區(qū)域,前者放入動參考系中,后者放入靜止參考系中進行模擬計算。使用標準k-ε湍流模式來模擬仿真湍流狀態(tài)下生物轉盤的流體流動。

2 計算結果與討論

轉盤表面液膜更新與轉盤速度分布密切相關[14,17],為了探討空氣與廢水在新型生物轉盤中的流動情況,分別采用2種轉速對該生物轉盤的流動情況進行模擬分析。

2.1 轉速為3 r/min下的模擬結果

在轉速為3 r/min,即線速度為2.35 m/min時,最后得到的混合相的速度云圖如圖5,速度矢量圖如圖6。由于新型生物轉盤呈順時針方向旋轉,上方為空氣,下方為污水。

圖5 盤片表面速度云圖(z=40 mm)

圖6 盤片表面速度矢量圖 (z=40 mm)

從圖5與圖6中可發(fā)現,在轉速為3 r/min時,空氣的流動速度遠大于下方污水的速度,其原因是空氣的黏度和密度遠小于污水的黏度和密度。上方空氣和下方污水會形成2個不同區(qū)域的循環(huán),空氣和污水沒有充分接觸,氣液兩相傳質不很明顯。在轉速為3 r/min時方孔內沒有明顯觀察到氣液兩相傳質效果。

2.2 轉速為6 r/min下模擬的結果

為改變轉速為3 r/min時氣液兩相傳質不明顯的問題,將轉速提高為6 r/min,即線速度為4.71 m/min。得到混合相的速度云圖如圖7,盤片表面方孔局部放大云圖如圖8,速度矢量圖如圖9,盤片表面局部放大矢量圖如圖10。

通過圖7與圖9可發(fā)現,盤片周圍以及盤片上方附近流速變大。在生物轉盤旋轉過程中,盤片的轉動會加快盤片周圍的水和空氣的流動,氣液兩相不斷接觸混合?？焖倭鲃拥目諝饽苁股锬ぶ械奈⑸镂崭嗟难鯕庥糜谘趸饔梅纸馑械挠袡C物。

圖7 盤片表面速度云圖(z=40 mm)

圖8 盤片表面方孔局部放大云圖(z=40 mm)

圖9 盤片表面速度矢量圖(z=40 mm)

圖10 盤片表面圓孔局部放大矢量圖(z=40 mm)

2.3 盤片表面y方向分布

因該生物轉盤的形狀以及尺寸大小,在轉速為6 r/min,即線速度為4.71 m/min時,在y=0面所截取的速度云圖如圖11,速度矢量圖如圖12。

圖11 盤片表面速度云圖(y=0 mm)

圖12 盤片表面速度矢量圖(y=0 mm)

圖11與圖12分別為y=0 mm位置處盤片表面速度云圖和速度矢量圖。這2個圖體現了y方向所看到的生物膜的分布情況以及氣液兩相的流動情況。相較于傳統(tǒng)的實心轉盤中氣液兩相只有徑向流動而無軸向流動,本文所研究的新型生物轉盤盤片因為方孔的存在是有軸向流動,這更加有利于氣液兩相傳質過程,利于生物膜的新老更替,生物轉盤能更有效地凈化污水。

根據以上結果及分析,新型生物轉盤的操作轉速不同,所得到的流場結果也不同。本文對比6 r/min和3 r/min時轉速發(fā)現,在轉速為6 r/min時,空氣和污水會在一起形成一個循壞,空氣和污水都會進入方孔,兩相之間有所接觸。這是因為轉盤的快速轉動會帶動上方的空氣進入槽中水體,轉盤上升時會有少量夾帶的液體污水,因此氣-液兩相接觸面增大,氣液傳質效果更好。

分析圖5—9發(fā)現,無論是在何種轉速,所用方槽的左右下角均存在死區(qū),所以后續(xù)需要對實驗用方槽進行改進。

3 結論

1) 該新型生物轉盤能夠增大轉盤盤片表面生物膜的分布區(qū)域,使生物膜分布表面變大,而不是只分布在一個方向上。

2) 該新型生物轉盤改進了生物轉盤對污水的凈化區(qū)域,提高了生物轉盤的利用率,能達到更好的污水凈化作用。

3) 為避免死區(qū)出現,污水水槽可改進為梯形槽。相比于轉速為3 r/min時,轉速為6 r/min時能夠增加氣-液兩相的接觸面,達到更好的傳質效果,可為生物轉盤的操作參數優(yōu)化提供參考。