用CFD方法模擬膜生物反應(yīng)器內(nèi)部流場(chǎng)及布?xì)鈨?yōu)化

王 捷，吳 義，羅 南，張宏偉，肖志偉，環(huán)國(guó)蘭

（1.天津工業(yè)大學(xué)a.中空纖維膜材料與膜過(guò)程教育部重點(diǎn)實(shí)驗(yàn)室；b.環(huán)境與化學(xué)工程學(xué)院，天津 300387；2.中國(guó)科學(xué)院生態(tài)環(huán)境研究中心污染控制研究室，北京 100085；3.北京市交通運(yùn)輸環(huán)境保護(hù)中心，北京 100013；4.天津膜天膜科技股份有限公司，天津 300457）

用CFD方法模擬膜生物反應(yīng)器內(nèi)部流場(chǎng)及布?xì)鈨?yōu)化

王 捷1a，1b，吳 義1b，羅 南2，張宏偉1a，肖志偉3，環(huán)國(guó)蘭4

（1.天津工業(yè)大學(xué)a.中空纖維膜材料與膜過(guò)程教育部重點(diǎn)實(shí)驗(yàn)室；b.環(huán)境與化學(xué)工程學(xué)院，天津 300387；2.中國(guó)科學(xué)院生態(tài)環(huán)境研究中心污染控制研究室，北京 100085；3.北京市交通運(yùn)輸環(huán)境保護(hù)中心，北京 100013；4.天津膜天膜科技股份有限公司，天津 300457）

采用Eulerian多相流模型對(duì)膜生物反應(yīng)器進(jìn)行氣液兩相流數(shù)值模擬，對(duì)“對(duì)齊”和“對(duì)齊導(dǎo)流”兩種反應(yīng)器構(gòu)型內(nèi)部流場(chǎng)氣含率、速度場(chǎng)和膜面液體流速進(jìn)行了分析比較，并就布?xì)夥绞竭M(jìn)行了分析和優(yōu)化，同時(shí)借助縮小實(shí)驗(yàn)?zāi)Ｐ蛯?duì)模擬結(jié)果進(jìn)行了實(shí)驗(yàn)驗(yàn)證。結(jié)果表明：兩種反應(yīng)器構(gòu)型流場(chǎng)內(nèi)的氣泡呈現(xiàn)出匯集于膜組件中心位置，然后在膜組件頂端散開(kāi)的流動(dòng)狀態(tài)；導(dǎo)流作用對(duì)反應(yīng)器內(nèi)氣含率分布的影響不大，但對(duì)速度場(chǎng)分布特性影響顯著；體積縮小100倍的對(duì)齊導(dǎo)流模擬裝置中的氣液流動(dòng)狀態(tài)與CFD模擬結(jié)果基本一致；通過(guò)對(duì)布?xì)夥绞降膬?yōu)化模擬發(fā)現(xiàn)，不均勻布?xì)夥绞娇梢愿纳茪怏w分布狀況，提高反應(yīng)器內(nèi)氣含率和流場(chǎng)的湍流強(qiáng)度。

CFD；膜生物反應(yīng)器；流場(chǎng)模擬；曝氣；優(yōu)化

為了減緩膜面污染及維持高濃度的污泥生化降解能力，膜生物反應(yīng)器（Membrane Bioreactor／MBR）運(yùn)行過(guò)程中需要高強(qiáng)度的曝氣作用。因此與其它污水處理工藝相比，MBR運(yùn)行中需要很高的能耗，其中膜擦洗曝氣、生化工藝曝氣能耗占系統(tǒng)運(yùn)行能耗的76%［1］。因此了解各種布?xì)夥绞较翸BR內(nèi)部流場(chǎng)特性，對(duì)于合理設(shè)計(jì)MBR構(gòu)型及布?xì)夥绞揭赃_(dá)到節(jié)能降耗的目的至關(guān)重要。由于反應(yīng)器內(nèi)部流場(chǎng)復(fù)雜多變，因此僅通過(guò)實(shí)驗(yàn)手段獲得整體流場(chǎng)分布特性受到一定的局限性。近年來(lái)，隨著計(jì)算流體力學(xué)（Computational Fluid Dynamics／CFD）技術(shù)及其軟件的迅速發(fā)展，已在反應(yīng)器設(shè)計(jì)、優(yōu)化中得到了廣泛的應(yīng)用［2－4］，它的優(yōu)勢(shì)在于應(yīng)用數(shù)值軟件求解描述過(guò)程，可以實(shí)現(xiàn)過(guò)程設(shè)計(jì)、優(yōu)化以及放大。如今已成為解決工程問(wèn)題的一種新的手段［5－8］。

MBR內(nèi)部流場(chǎng)數(shù)值模擬屬于氣液兩相流范疇［9－12］，氣液兩相流是常見(jiàn)的一種物質(zhì)流動(dòng)，廣泛應(yīng)用于工程和環(huán)境領(lǐng)域中［10－12］。在 MBR 運(yùn)行過(guò)程中，曝氣是重要的環(huán)節(jié)，曝氣的作用一方面是向反應(yīng)器內(nèi)充空氣，以保證好氧微生物所需的溶解氧；另一方面是吹掃膜組件，氣泡對(duì)膜面的沖刷剪切以減緩膜污染［13－15］。對(duì)于 MBR構(gòu)型及布?xì)夥绞降膬?yōu)化，由于工程實(shí)際裝置的復(fù)雜性及實(shí)驗(yàn)條件的局限性，有些問(wèn)題在實(shí)驗(yàn)室難以完成。而計(jì)算流體力學(xué)模擬可以實(shí)現(xiàn)及預(yù)見(jiàn)工程放大可能出現(xiàn)的一些問(wèn)題。Kroll等［16］運(yùn)用CFD對(duì)應(yīng)用最廣泛的MBR反應(yīng)池構(gòu)型進(jìn)行了二維數(shù)值模擬，對(duì)其進(jìn)行布?xì)夥绞絻?yōu)化，考察了特定曝氣強(qiáng)度和布?xì)夥绞綄?duì)MBR運(yùn)行特性的影響。

在上述研究的基礎(chǔ)上，運(yùn)用CFD數(shù)值模擬方法對(duì)曝氣引發(fā)的MBR內(nèi)流場(chǎng)的狀態(tài)進(jìn)行了分析，同時(shí)比較了曝氣管正對(duì)齊膜組件下方（簡(jiǎn)稱(chēng)“對(duì)齊”布?xì)猓┘啊皩?duì)齊”布?xì)鈺r(shí)在膜單元兩側(cè)加上曝氣導(dǎo)流板（簡(jiǎn)稱(chēng)“對(duì)齊”導(dǎo)流布?xì)猓?gòu)型的改進(jìn)給MBR內(nèi)部流場(chǎng)帶來(lái)的改變，由此提出了布?xì)夥绞降膬?yōu)化方案，以便為大型MBR工程設(shè)計(jì)提供相應(yīng)的理論基礎(chǔ)。

1 研究對(duì)象與方法

1.1 MBR兩種布?xì)鈽?gòu)型的CFD模型

研究以MBR中反應(yīng)器內(nèi)的流體為分析對(duì)象，以天津膜天膜科技股份有限公司的FP-AⅡ15型簾式膜組件為基礎(chǔ)模型。反應(yīng)器及膜組件結(jié)構(gòu)參數(shù)如下：反應(yīng)器尺寸為1 720 mm×1 050 mm×3 000 mm，膜組件尺寸為534 mm×450 mm×1 523 mm，單個(gè)組件膜面積20 m2，組件間距80 mm，膜組件距反應(yīng)池邊距300 mm，底部采用穿孔管曝氣，曝氣孔4 mm，每個(gè)膜組件下均勻排布十個(gè)曝氣孔，氣孔間距40 mm，排間距120 mm。圖1為兩種布?xì)夥绞较碌娜S反應(yīng)器簡(jiǎn)化模型。對(duì)“對(duì)齊導(dǎo)流”布?xì)馊鐖D1（a）及“對(duì)齊”布?xì)馊鐖D1（b）進(jìn)行三維數(shù)值模擬，分析其內(nèi)部流場(chǎng)特性。為模擬求解方便，假設(shè)反應(yīng)器內(nèi)均勻。由于本研究的重點(diǎn)是曝氣引發(fā)的反應(yīng)器內(nèi)流場(chǎng)狀態(tài)，且膜組件過(guò)濾引發(fā)的液體流動(dòng)速度遠(yuǎn)小于曝氣引發(fā)的氣、液混合流動(dòng)速度，因此將膜簡(jiǎn)化為不透水的無(wú)滑移壁面，模擬中不考慮膜過(guò)濾。

圖1 三維模型

1.2 反應(yīng)器的網(wǎng)格劃分及邊界條件的確定

借助Gambit 2.3.16劃分反應(yīng)器網(wǎng)格結(jié)構(gòu)，網(wǎng)格決定著模型的準(zhǔn)確度。生成網(wǎng)格不但要考慮其細(xì)密程度來(lái)提高模擬過(guò)程的準(zhǔn)確程度，又要避免由于網(wǎng)格過(guò)于復(fù)雜帶來(lái)的計(jì)算復(fù)雜性。該研究采用四面體網(wǎng)格，分別畫(huà)了網(wǎng)格數(shù)為140萬(wàn)和180萬(wàn)2種，模擬結(jié)果相同，證明了網(wǎng)格的無(wú)關(guān)性，網(wǎng)格整體與局部圖如圖2。反應(yīng)器模擬活性污泥混合液，其混合液粘度及密度均為實(shí)測(cè)值，分別為0.007 03 kg·m－1·s－1和1 010 kg·m－3。從流體實(shí)際流動(dòng)情況出發(fā)，采用非定常流動(dòng)模式，簡(jiǎn)單且計(jì)算精度較高的標(biāo)準(zhǔn)k-ε兩方程湍流模型（該模型在一定程度上考慮了流動(dòng)場(chǎng)中各點(diǎn)的湍動(dòng)能傳遞和流動(dòng)的歷史作用）。求解器選用求解不可壓流場(chǎng)且需要內(nèi)存較小的壓力 速度耦和SIMPLE算法，穩(wěn)態(tài)3D分離顯式解算器，對(duì)流項(xiàng)選取一階迎風(fēng)離散格式。設(shè)計(jì)液面高度為2 600 mm，曝氣孔進(jìn)氣速度為0.995 2 m·s－1。設(shè)置入口邊界條件為速度入口邊界條件，出口邊界條件為壓力出口邊界條件，過(guò)濾器壁面及膜組件壁面為無(wú)滑移壁面邊界條件，反應(yīng)器為流場(chǎng)邊界條件。

圖2 三維模型網(wǎng)格劃分圖及局部網(wǎng)格

2 模擬分析

模擬過(guò)程中將污泥的固相和水的液相視作“假想的均一混合相”，則可將三相流問(wèn)題簡(jiǎn)化為“氣 液兩相流”問(wèn)題。通過(guò)對(duì)反應(yīng)器不同布?xì)夥绞降哪M，結(jié)果得到反應(yīng)器內(nèi)多個(gè)截面的氣液兩相速度場(chǎng)及速度流線(xiàn)等信息，以此來(lái)了解反應(yīng)器內(nèi)水力情況。通過(guò)實(shí)驗(yàn)室縮小模型與模擬結(jié)果的對(duì)比，驗(yàn)證模擬結(jié)果，并對(duì)改進(jìn)的布?xì)夥绞竭M(jìn)行了模擬計(jì)算。

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

3.1 對(duì)齊曝氣時(shí)反應(yīng)器內(nèi)流場(chǎng)特性

3.1.1 反應(yīng)器內(nèi)氣含率分布 圖3為模型的簡(jiǎn)化示意圖，實(shí)驗(yàn)選取了X＝0，Y＝0，Z＝0幾個(gè)截面從橫向和縱向分板反應(yīng)器內(nèi)流場(chǎng)的變化。圖4和5為反應(yīng)器不同截面及軸向上的氣含率分布情況，從圖4可以看出反應(yīng)器曝氣板底部的氣含率很低，從曝氣板到膜組件底部位置氣含率開(kāi)始顯著增加。在曝氣器上方的位置氣體沿膜組件中間位置比較集中的上升，然后在膜組件頂端散開(kāi)，這一模擬結(jié)果與Nicolas模擬結(jié)果一致。氣體在流場(chǎng)中的這種分布可能是由于從曝氣孔出來(lái)的小氣泡形態(tài)比較穩(wěn)定，對(duì)周?chē)鷼馀莸男螒B(tài)及速度變化沒(méi)有明顯的影響，組件之間的阻擋作用使大量的小氣泡在膜組件中心位置比較集中的緩慢上升，隨著氣泡的逐漸上升，到達(dá)組件頂部后，氣體沒(méi)有組件的阻擋作用，氣泡尺寸越來(lái)越大，較大的氣泡容易變形破裂或與周?chē)馀萑诤?，最終大量大氣泡在膜組件頂部破裂散開(kāi)，從膜組件頂端散開(kāi)的氣體在水流的作用下沿著壁面被帶到反應(yīng)器下部，在反應(yīng)器內(nèi)部形成循環(huán)流動(dòng)。從圖5也可以看出，氣體主要集中在反應(yīng)器中心位置，反應(yīng)器底部和反應(yīng)器兩側(cè)氣含率較低。

圖3 模型截面選取示意圖（原點(diǎn)坐標(biāo)位于體中心）

圖4 反應(yīng)器X＝0、Z＝0（如圖3所示）不同截面氣含率（%）分布圖

圖5 反應(yīng)器內(nèi)部X＝0、Z＝0軸線(xiàn)上氣含率分布

3.1.2 反應(yīng)器內(nèi)液速分布 為進(jìn)一步分析流體在膜面的流動(dòng)特性，研究考察了曝氣槽中平行膜面以及垂直膜面幾個(gè)典型截面上流體的速度等值線(xiàn)及矢量場(chǎng)分布情況。圖6為各截面液體速度云圖，圖7為各截面速度矢量圖。速度云圖及流線(xiàn)圖中的顏色深淺代表速度絕對(duì)值的大小，箭頭指向代表流體的運(yùn)動(dòng)方向。從圖6中可以看出，膜單元中間，及膜單元頂部存在大量的高流速區(qū)域，而反應(yīng)器內(nèi)底部、靠近反應(yīng)器壁面及膜單元外圍的膜組件存在一定范圍的水力死區(qū)。曝氣帶動(dòng)的混合液流動(dòng)呈現(xiàn)出了不均勻的狀態(tài)，這會(huì)致使膜面污染的不均勻性，外側(cè)膜組件受到曝氣清洗的作用較弱，導(dǎo)致污染較嚴(yán)重。

圖6 “對(duì)齊”布?xì)夥磻?yīng)器各截面流體速度分布云圖

圖7 平行膜面以及垂直膜面各剖面流體速度的矢量分布圖

分析圖7可知，流體在反應(yīng)器內(nèi)呈現(xiàn)為循環(huán)運(yùn)動(dòng)，漩渦的形成有利于氣液兩相的充分混合。氣泡的軌跡傾向于流向膜組件中心位置，然后在膜組件頂端散開(kāi)。從平行膜面截面速度矢量圖7（a）可以看出，氣泡帶動(dòng)液體在膜組件之間形成的一束高流速流體比較集中的向上流動(dòng)，到達(dá)組件上部時(shí)由于大氣泡的相互作用形成大量的渦旋，大量渦旋的擾動(dòng)使流體在組件之間形成內(nèi)循環(huán)流動(dòng)，這有利于提高膜面的沖刷作用。但膜組件下部?jī)蓚?cè)存在一定范圍的低流速區(qū)域，這些區(qū)域渦旋不明顯，流體擾動(dòng)不強(qiáng)烈，流體對(duì)膜面的沖刷作用較弱。

結(jié)合圖6和7可以看出，流體在反應(yīng)器內(nèi)的流動(dòng)呈現(xiàn)一定的不均勻性，組件中間的升流區(qū)流速較大而兩側(cè)的流速較小，這使得膜污染的分布也呈現(xiàn)一種不均勻性。同時(shí)流體在膜單元中間位置形成流速較大的升流區(qū)，到達(dá)膜組件頂部擴(kuò)散，流體往反應(yīng)器兩側(cè)向下流動(dòng)，再與膜單元外圍膜組件間隔區(qū)域形成循環(huán)流動(dòng)。

3.2 導(dǎo)流作用對(duì)反應(yīng)器內(nèi)部流場(chǎng)的影響

3.2.1 導(dǎo)流板對(duì)氣含率分布的影響 圖8、圖9為導(dǎo)流結(jié)構(gòu)的反應(yīng)器不同截面及軸向上的氣含率分布情況。對(duì)比圖4和8可以看出導(dǎo)流板的增加，使氣體分布有所改變，但變化不大。分析圖9數(shù)據(jù)也可以看出，反應(yīng)器軸向上氣含率分布基本一致，導(dǎo)流布?xì)鈺r(shí)氣含率略高，但總的來(lái)說(shuō)導(dǎo)流作用對(duì)反應(yīng)器內(nèi)氣含率分布的影響不是很大。

圖8 反應(yīng)器不同截面氣含率（%）分布圖（截面選取圖4）

圖9 反應(yīng)器內(nèi)部X＝0、Z＝0軸線(xiàn)上氣含率分布

3.2.2 導(dǎo)流板對(duì)流場(chǎng)液速分布的影響 從圖10可以看出，在導(dǎo)流板的作用下，反應(yīng)器內(nèi)速度場(chǎng)分布趨于均勻，高流速區(qū)域明顯增加，底部流體流速有所提升，沒(méi)有出現(xiàn)大范圍的低流速區(qū)。對(duì)比圖6和10的反應(yīng)器內(nèi)流速分布可以看出，導(dǎo)流板在一定程度上提升了膜單元外側(cè)組件周?chē)幕旌弦毫魉?，能夠使外?cè)的組件受到更為均勻的水力沖刷作用，對(duì)保障膜單元中各組件的污染均勻性具有一定作用。

圖10 “對(duì)齊”導(dǎo)流布?xì)庀路磻?yīng)器各截面流體速度分布云圖（截面選取同圖6）

由圖11可知，加導(dǎo)流板之后，流體在反應(yīng)器內(nèi)仍為循環(huán)運(yùn)動(dòng)，與“對(duì)齊”布?xì)獠煌氖?，?dǎo)流板的設(shè)置使流場(chǎng)內(nèi)的循環(huán)流動(dòng)明顯增強(qiáng)。流場(chǎng)內(nèi)液相質(zhì)點(diǎn)形成的旋渦非常明顯，導(dǎo)流板的設(shè)置增強(qiáng)了反應(yīng)器內(nèi)流體的循環(huán)流動(dòng)，形成了非常明顯的升流區(qū)和降流區(qū)，這有利于流體濃度及速度的分布均勻。通過(guò)圖12中這兩種反應(yīng)器構(gòu)型不同截面及膜表面平均速度的對(duì)比也可以發(fā)現(xiàn)，導(dǎo)流結(jié)構(gòu)能提高反應(yīng)器內(nèi)部的液體流速，改變流場(chǎng)分布情況，使得膜表面的流體速度顯著增加，且分布更均勻，有利于更均勻的對(duì)膜組件沖刷與水力清洗。

圖11 平行膜面以及垂直膜面各剖面流體速度的矢量分布

圖12 對(duì)齊與對(duì)齊導(dǎo)流布?xì)夥磻?yīng)器各截面與各膜面平均速度對(duì)比

3.3 布?xì)夥绞浇Y(jié)構(gòu)優(yōu)化模擬

通過(guò)對(duì)“對(duì)齊”與“對(duì)齊導(dǎo)流”這兩種反應(yīng)器構(gòu)型的模擬發(fā)現(xiàn)，均勻排布?xì)饪讜r(shí)，膜單元流場(chǎng)中氣泡軌跡流向膜單元中心位置，然后在膜單元頂端散開(kāi)。集中上升的氣泡不利于氧的傳質(zhì)和膜組件的清洗。為了改善反應(yīng)器內(nèi)部的氣含率分布情況，研究對(duì)曝氣孔的排布進(jìn)行了改進(jìn)。對(duì)布?xì)饪孜恢眠M(jìn)行了調(diào)整，將10個(gè)曝氣孔分成2組，2組5個(gè)曝氣孔均勻排布，2組之間間距100 mm。為了求解方便，實(shí)驗(yàn)對(duì)單個(gè)膜組件的小型反應(yīng)器進(jìn)行了模擬，實(shí)驗(yàn)?zāi)Ｐ腿鐖D13。其網(wǎng)格劃分、計(jì)算模型選擇與全尺寸模型模擬相同。

圖13 不均勻布?xì)馊S模擬模型

改進(jìn)曝氣形式的內(nèi)部流場(chǎng)氣含率分布如圖14。對(duì)比圖4可以看出，這一布?xì)夥绞绞狗磻?yīng)器內(nèi)形成2束氣流，這2束氣流在膜組件上部匯合并散開(kāi)。不均勻布?xì)夥绞礁纳屏藲怏w分布狀況，使氣體分布更均勻，提高了反應(yīng)器內(nèi)氣含率。

圖14 反應(yīng)器中心位置不同截面氣含率

從圖15、16中對(duì)比“對(duì)齊”均勻布?xì)夂透倪M(jìn)布?xì)鈨煞N布?xì)夥绞浇っ嫣幍牧黧w速度及速度矢量分布可知：均勻布?xì)鈺r(shí)，流場(chǎng)內(nèi)高流速區(qū)域集中在膜組件中心位置與膜組件上部，膜組件兩側(cè)、反應(yīng)器壁面及底部有大范圍的低流速區(qū)域；不均勻布?xì)鈺r(shí)，流場(chǎng)內(nèi)速度分布均勻，反應(yīng)器內(nèi)膜組件及壁面附近均有較大的流速。反應(yīng)器底部也有較大的流速。通過(guò)對(duì)圖16中各特征區(qū)域局部放大的流線(xiàn)軌跡（如圖17所示）進(jìn)行分析，可以對(duì)上述2種速度分布做出解釋。分析圖17（a）、17（c）發(fā)現(xiàn)：均勻布?xì)鈺r(shí)組件上端有較大的流速，循環(huán)流動(dòng)不明顯，這是由于均勻曝氣時(shí)小氣泡集中上升，流體速度主要集中在這一區(qū)域；不均勻布?xì)鈺r(shí)，氣體分成兩束上升，分布較廣，組件上部區(qū)域的速度要小于均勻布?xì)夥绞剑h(huán)流非常明顯。通過(guò)對(duì)比圖17（b）、17（d）可以發(fā)現(xiàn)：均勻布?xì)鈺r(shí)，流體速度主要集中在膜組件中心及組件上部，沒(méi)有在反應(yīng)器兩側(cè)形成循環(huán)流動(dòng)，使得反應(yīng)器底部也幾乎沒(méi)有循環(huán)流動(dòng)，流體在這一區(qū)域速度較低；不均勻布?xì)鈺r(shí)，在2束氣流的作用下，液體在膜組件兩側(cè)形成了2個(gè)方向相反的渦流。渦流提高了反應(yīng)器內(nèi)流場(chǎng)的湍流強(qiáng)度，使得反應(yīng)器壁面附近流體速度變大，如圖15（b）所示。由于漩渦將流體帶至反應(yīng)器底部，使底部壁面附近流體產(chǎn)生了一定的流速，這能夠有效地改善反應(yīng)器底部低流速區(qū)。

圖15 “對(duì)齊”布?xì)庀虏煌瑲饪着挪冀っ妫? mm）流體速度分布云

16 “對(duì)齊”布?xì)庀虏煌瑲饪着挪冀っ媪黧w速度矢量圖

3.4 實(shí)驗(yàn)與模擬結(jié)果對(duì)比

影響反應(yīng)器內(nèi)部流場(chǎng)的因素很多，譬如：反應(yīng)器構(gòu)型（有無(wú)導(dǎo)流板），曝氣孔的排布方式（均勻和非均勻布置曝氣孔），曝氣孔的大小等。而本研究的重點(diǎn)是從宏觀(guān)上考察一定曝氣強(qiáng)度下反應(yīng)器不同構(gòu)型及不同布?xì)夥绞綄?duì)反應(yīng)器內(nèi)部流場(chǎng)的影響，并且結(jié)合小型實(shí)驗(yàn)驗(yàn)證數(shù)值模擬結(jié)果準(zhǔn)確性。研究對(duì)體積縮小100倍的對(duì)齊導(dǎo)流模擬裝置進(jìn)行了流動(dòng)實(shí)驗(yàn)，圖18（a）為對(duì)齊導(dǎo)流曝氣小型實(shí)驗(yàn)?zāi)Ｐ?，模型尺寸?.72 mm×1.05 mm×3.00 mm，曝氣孔直徑為1.50 mm，實(shí)驗(yàn)各項(xiàng)參數(shù)的獲得均通過(guò)佛洛德數(shù)相似性準(zhǔn)則計(jì)算得到。為了便于觀(guān)察反應(yīng)器內(nèi)部流場(chǎng)運(yùn)動(dòng)特性，選取密度與流體相近的壓實(shí)海綿顆粒進(jìn)行示蹤，并通過(guò)高清相機(jī)（Canon，LEGRIA HF200）連續(xù)采集反應(yīng)器內(nèi)流體流動(dòng)狀態(tài)的圖片與數(shù)值模擬結(jié)果進(jìn)行對(duì)比分析。

17 “對(duì)齊”布?xì)獠煌瑲饪着挪枷率噶繄D各局部速度流線(xiàn)

從圖18（d）中可看出壓實(shí)海綿示蹤顆粒在反應(yīng)器內(nèi)部做旋轉(zhuǎn)運(yùn)動(dòng)，這與模擬結(jié)果圖11中的的速度矢量圖中液體運(yùn)動(dòng)方向一致。圖18（b）、（c）中的觀(guān)測(cè)到從曝氣孔出來(lái)的大量小氣泡在組件及導(dǎo)流板的作用下在膜組件中心位置比較集中的緩慢上升。隨著氣泡的逐漸上升，氣泡尺寸越來(lái)越大，大氣泡上升速度相對(duì)較快并且容易變形破裂或與周?chē)馀萑诤?，最終大量大氣泡在膜組件頂部破裂散開(kāi)，從膜組件頂端散開(kāi)的氣體在水流的作用下沿著壁面被帶到反應(yīng)器下部，在反應(yīng)器內(nèi)部形成循環(huán)流動(dòng)。從圖18（c）可看出延曝氣孔往上流體紊動(dòng)程度明顯增大。這可能是由于組件及導(dǎo)流板對(duì)氣泡的束縛作用減弱，且較大的氣泡對(duì)周?chē)鷼馀莸男螒B(tài)及速度變化產(chǎn)生的影響較大，對(duì)于其上方的氣泡有頂托作用，對(duì)位于的氣泡其尾渦會(huì)加速氣泡的上升，而對(duì)位于其側(cè)面稍遠(yuǎn)的氣泡則有抑制作用甚至?xí)蛊鋫?cè)面的氣泡小幅度下降。而組件底端距曝氣口較近的地方，剛從曝氣孔出來(lái)的大量小氣泡受組件及導(dǎo)流板的阻滯作用，上升及擴(kuò)散速度相對(duì)于大氣泡較慢，形態(tài)比大氣泡穩(wěn)定，但小氣泡對(duì)周?chē)鷼馀莸男螒B(tài)及速度變化沒(méi)有明顯的影響。氣泡軌跡與模擬結(jié)果相同，由此說(shuō)明整個(gè)模擬過(guò)程中參數(shù)、計(jì)算模型及計(jì)算公式的選擇都比較合理，能使模擬結(jié)果和實(shí)際情況很好的吻合。

圖18 流動(dòng)模擬實(shí)驗(yàn)

4 結(jié) 論

1）通過(guò)模擬氣孔均勻排布的布?xì)夥绞桨l(fā)現(xiàn)反應(yīng)器內(nèi)部氣泡軌跡傾向于流向膜組件中心位置，然后在膜組件頂端散開(kāi)，從膜組件頂端散開(kāi)的氣體在水流的作用下，沿著壁面被帶到反應(yīng)器下部，在整個(gè)反應(yīng)器內(nèi)形成循環(huán)。

2）對(duì)比“對(duì)齊”與“對(duì)齊導(dǎo)流”反應(yīng)器內(nèi)流體速度分布場(chǎng)發(fā)現(xiàn)，導(dǎo)流板的導(dǎo)流作用，對(duì)流場(chǎng)特性影響顯著，能很好的改善反應(yīng)器內(nèi)部速度分布不均的情況。導(dǎo)流板的設(shè)置增加了流場(chǎng)局部區(qū)域內(nèi)的紊流強(qiáng)度，強(qiáng)化了反應(yīng)器內(nèi)流體的循環(huán)流動(dòng)，形成了非常明顯的升流區(qū)和降流區(qū)，有利于氣體在反應(yīng)器內(nèi)的擴(kuò)散和流體濃度分布均勻。

3）通過(guò)對(duì)不均勻布?xì)獾母倪M(jìn)結(jié)構(gòu)的模擬發(fā)現(xiàn)，氣泡軌跡仍為沿著膜組件中心集中，但改進(jìn)結(jié)構(gòu)改善了氣體分布狀況，提高了反應(yīng)器內(nèi)氣含率和流場(chǎng)的湍流強(qiáng)度，液體在膜組件兩側(cè)形成了2個(gè)方向相反的渦流，使反應(yīng)器內(nèi)流體速度分布更均勻，能改善膜面的不均勻污染。

4）對(duì)體積縮小100倍的對(duì)齊導(dǎo)流小型模型裝置進(jìn)行了流動(dòng)實(shí)驗(yàn)，模擬結(jié)果與實(shí)驗(yàn)結(jié)果基本一致，即模擬過(guò)程中參數(shù)、計(jì)算模型及計(jì)算公式的選擇都比較合理。

［1］Verrecht B，Maere T，Nopens I，et al.The cost of a large-scale hollow fiber MBR ［J］.Water Research，2010，44（10）：5274-5283.

［2］張景騫，王艷紅，曹越.膜生物反應(yīng)器處理技術(shù)的研究現(xiàn)狀［J］.環(huán)境科學(xué)與管理，2010，35（1）：66-69.

Zhang J Q，Wang Y H，Cao Y.Research situation of membrane bioreactor processing technology［J］.Environmental Science and Management，2010，35（1）：66-69.

［3］Brannock M，Wang Y，Shane C.Mixing characterisation of full-scale membrane bioreactors：CFD modelling with experimental validation ［J］.Water Research，2010，44（10）：3181-3191.

［4］Wang Y，Brannock M，Shane C，et al. CFD simulations of membrane filtration zone in a submerged hollow fiber membrane bioreactor using a porous media approach［J］.Journal of Membrane Science，2010，363（1／2）：57-66.

［5］Darcovich K，Dal-Cin M M，Ballèvreb S，et al.CFD-assisted thin channel membrane characterization module design［J］.Journal of Membrane Science，1997，124（2）：181-193.

［6］Brans G，van der Sman R G M，Schro C G P H，et al.Optimization of the membrane and pore design for micro-machined membranes ［J］.Journal of Membrane Science，2006，278（1／2）：239-250.

［7］員文權(quán)，楊慶峰.計(jì)算流體動(dòng)力學(xué)在反滲透膜分離過(guò)程中的應(yīng)用［J］.化工進(jìn)展，2008，27（9）：1357-1363.

Yuan W Q，Yang Q F.Application of computational fluid dynamics in reverse osmosis membrane separation process［J］.Chemical Industry and Engineering Progress，2008，27（9）：1357-1363.

［8］Mercier-Bonin M，Langne C，F(xiàn)onade C.Influence of a gas／liquid two-phase flow on the ultrafiltration and microfiltration performance：case of a ceramic flat sheet membrane［J］.Journal of Membrane Science，2000，180（9）：93-102.

［9］Taha T，Cui Z F.CFD modeling of gas sparged ultrafiltration in tubular membranes ［J］.Journal of Membrane Science，2002，2（10）：13-27.

［10］Phattaranawik J，F(xiàn)ane A G，Pasquier A C S，et al.Membrane bioreactor with bubble-size transformer：design and fouling control［J］.AICh E Journal，2007，53（4）：243-248.

［11］Bodmer S，Imark C，Kneubühl M.Biogenic amines in foods：histaine and food processing ［J］.Inflammation Research，1999，48（2）：296-300.

［12］Latorre-Moratalla M L，Bover-Cid S，Aymerich T，et al.Aminogenesis control in fermented sausages manufactured with pressurized meat batter and starter culture［J］.Meat Science，2007，75（7）：460-469.

［13］Docom G，Puech F P，Cabassud C.Air sparging with flat sheet nanofiltration a link between wall shear stress and flux enhancement［J］.Desalination，2002，14（5）：97-102.

［14］Martinelli L，Guigui C，Line A.Characterisation of hydrodynamics induced by air injection related to membrane fouling behaviour［J］.Desalination，2010，250（2）：587-591.

［15］Fazzolari A， Gauger N R，Brezillon J. Efficient aerodynamic shape optimization in MDO context［J］.Journal of Computationl and Applied Mathematics，2007，20（3）：548-560.

［16］Kroll N，Gauger N R，Brezillon J，et al.Flow simulation and shape optimization for aircraft design［J］.Journal of Computationl and Applied Mathematics，2007，20（3）：397-411.

（編輯 郭 飛）

CFD Simulation and Optimization of Internal Hydrodynamic Characteristics and Aerating Distribution in Membrane Bioreactor

Wang Jie1a，1b，Wu Yi1b，Luo Lan2，Zhang Hongwei1a，Xiao Zhiwei3，Huan Guolan4

（1a.State Key Laboratory of Hollow Fiber Membrane Materials and Processes，Tianjin Polytechnic University，1b.School of Environmental and Chemical Engineering，Tianjin Polytechnic University，Tianjin 300387，P.R.China；2.Department of Water Pollution Control Technology，Research Center for Eco-Environmental Science，Chinese Academy of Scionce，Beijing 100085，P.R.China；3.Environmental Protection Center Ministry of Transport，Beijing 100013，P.R.China；4.Tianjin MOTIMO Membrane Thechnology CO.，LTD，Tianjin，300457，P.R.China）

The simulation of Eulerian multiphase flow model with two-phase（air ＆ water）in membrane bioreactor were aimed at obtaining and comparing gas hold-up distribution，velocity field and membrane surface rate among the two reactor internal configuration of“alignment”and “alignment”diversion.Furthermore，the gas distribution mode was analyzed as well as optimized.The mini-model experiment was conducted to validate the simulation result.It is shown that：the bubbles track of internal field tended to flow to the membrane assembly center position，and then dispersed on the top of membrane modules in the two reactor configurations；the influence of diversion on gas holdup distribution was in apparent while notable on the velocity field distribution.The simulation results of an experimental model shrunk 100 times was proved to be consistent with the experimental results.Through an optimization simulation of gas distribution mode，the non-uniform gas distribution mode was proved that it can improve the range of gas distribution and the turbulence intensity.

CFD；membrane bioreactor；flow field simulation；aeration；optimization

X 703.1

A

1674-4764（2014）02-0119-07

10.11835／j.issn.1674-4764.2014.02.019

2012-11-27

國(guó)家自然科學(xué)青年基金（51378349、51108314）；國(guó)家水體污染控制與治理科技重大專(zhuān)項(xiàng)（2011ZX07317-001）

王 捷（1979-），男，副教授，博士，主要從事膜法水處理技術(shù)的研究，（E-mail）wangjie＠tjpu.edu.cn。