混凝燒杯實(shí)驗數(shù)值模擬研究

史 偉 ，劉宏遠(yuǎn) ，朱海濤 ，田孝禾 ，周 勇 ，孫海平 ，鐘葉華

（1.浙江工業(yè)大學(xué) 建筑工程學(xué)院，浙江杭州 310014；2.嘉源給排水有限公司，浙江嘉興 314000）

0 引言

燒杯實(shí)驗是水廠確定混凝劑投加量的一種簡便有效的方法。水廠一般直接采用購置的攪拌設(shè)備進(jìn)行實(shí)驗，不會去關(guān)注燒杯和槳葉形狀，而燒杯和槳葉形狀的不同在一定程度上會影響到燒杯實(shí)驗中的水力流態(tài)。因此，其實(shí)驗結(jié)果為生產(chǎn)運(yùn)行提供的可參考性值得商榷。本文基于CFD（Computational Fluid Dynamics）的 fluent6.3軟件模擬方形燒杯和圓形燒杯與彎葉槳和平槳組合情況下燒杯內(nèi)部的流態(tài)變化，通過數(shù)值模擬和相關(guān)分析，研究獲得較為合理的燒杯和槳葉形式，以便更有效地確定混凝劑投加量，增強(qiáng)混凝效果，從而為凈水廠的生產(chǎn)運(yùn)行提供更為可靠的參考數(shù)據(jù)[1]。

1 研究方法

在試驗研究中，混凝燒杯實(shí)驗采取數(shù)值模擬和現(xiàn)場實(shí)驗相結(jié)合的方式進(jìn)行。即：一方面通過基于CFD的fluent6.3軟件對不同組合情況下的燒杯（圓形燒杯或方形燒杯）和槳葉（平槳葉或彎葉槳）進(jìn)行數(shù)學(xué)建模，從而進(jìn)行混凝燒杯實(shí)驗的理論模擬研究，通過模型模擬的速度云圖和距離燒杯底部不同位置高度處的參數(shù)進(jìn)行綜合分析流態(tài)變化對混凝效果的影響；另一方面，通過現(xiàn)場實(shí)驗的方法進(jìn)行實(shí)測驗證，研究上述的燒杯和槳葉組合情況下對濁度及高錳酸鹽指數(shù)的去除效果。綜合上述兩者結(jié)論確定較佳的燒杯形狀和槳葉的組合方式，從而為水廠的生產(chǎn)實(shí)踐提供參考。

1.1 實(shí)驗水質(zhì)特點(diǎn)

嘉興地處杭嘉湖平原地帶，河流水系密布，水流緩慢，水體水質(zhì)常年處于IV-V類，甚至劣V類。燒杯實(shí)驗原水取自水廠水源南郊河，實(shí)驗期間具體原水水質(zhì)指標(biāo)見表1所列。鑒于水源水質(zhì)的實(shí)際特點(diǎn)，該項實(shí)驗中的現(xiàn)場實(shí)驗考察的實(shí)驗參數(shù)為濁度和高錳酸鹽指數(shù)。

表1 實(shí)驗期間原水水質(zhì)指標(biāo)一覽表

1.2 攪拌和燒杯設(shè)備

實(shí)驗攪拌裝置采用JJ-4六聯(lián)數(shù)顯電動攪拌器，它具有轉(zhuǎn)速穩(wěn)、噪聲小的特點(diǎn)。同時，設(shè)有燈光照明功能，可以觀察礬花的形成過程，有利于實(shí)驗信息的觀察和記錄。在混凝燒杯實(shí)驗進(jìn)行過程中，可以根據(jù)實(shí)驗要求控制操作界面調(diào)節(jié)不同的轉(zhuǎn)速（范圍為40 rad/min～1 200 rad/min）。攪拌裝置見圖1所示。

圖1 JJ-4六聯(lián)數(shù)顯電動攪拌器裝置實(shí)景

圖2為混凝燒杯實(shí)驗中的燒杯和槳葉形狀圖。這兩種燒杯均是燒杯混凝實(shí)驗攪拌器配套使用的有機(jī)玻璃1.5L方形燒杯（簡稱方杯）和1 L圓形燒杯（簡稱圓杯）。圖2（c）、2（d）分別為自行加工的平槳，以及在攪拌器購買中提供的彎葉槳。燒杯實(shí)驗時槳葉下端距離燒杯底部10 mm,取樣口位置為燒杯液面下方25 mm處。燒杯和槳葉的具體參數(shù)見表2所列。

圖2 混凝燒杯實(shí)驗中的燒杯和槳葉形狀圖

表2 燒杯和槳葉參數(shù)具體信息表

1.3 混凝劑相關(guān)參數(shù)

實(shí)驗采用的混凝劑為液態(tài)聚合氯化鋁（PAC），密度1.24 g/mL，Al2O3含量在10.78%，鹽75.22%，不溶物含量0.08%。藥劑投加濃度為5%（以有效Al2O3百分比計）的溶液。

1.4 試驗操作流程

根據(jù)國標(biāo)水處理劑聚合氯化(GBl5892-2009)附錄A-混凝性能判定的試驗程序混凝控制指標(biāo)要求,并通過實(shí)驗現(xiàn)場燒杯混凝正交試驗的方法在原有攪拌工況條件[2]下進(jìn)行優(yōu)化。正交試驗研究表明，實(shí)驗快速攪拌轉(zhuǎn)速500 rad/min、快攪時間1 min,慢速攪拌轉(zhuǎn)速為60 rad/min、慢攪時間15 min，靜置沉淀時間20 min后取出上清液測定相關(guān)指標(biāo)，藥劑投加量為正交實(shí)驗確定的最優(yōu)投加量60 mg/L。

1.5 數(shù)值模擬fluent6.3

CFD數(shù)值模擬是從20世紀(jì)60年代中逐步發(fā)展起來的一種研究流體流動等物理現(xiàn)象的現(xiàn)代技術(shù)，通過計算機(jī)數(shù)值模擬可以得到流場內(nèi)部的各種細(xì)節(jié)[3,4]。近些年來，許多學(xué)者通過對混凝反應(yīng)器進(jìn)行數(shù)值模擬得出混凝效果受流體動力學(xué)條件的影響很大的結(jié)論[5]。本文基于CFD的fluent6.3對燒杯內(nèi)部三維穩(wěn)定流場進(jìn)行水力特性分析，能直觀清晰地觀察燒杯反應(yīng)器內(nèi)部流態(tài)變化，節(jié)省了粒子圖像測速法（PIV）監(jiān)測流場的高額費(fèi)用，同時為優(yōu)化水處理反應(yīng)器提供一種直觀易行的技術(shù)手段。

1.5.1 模型的建立

對不同燒杯與槳葉組合條件下的模型進(jìn)行數(shù)值模擬，采用多重參考系法（MRF）將燒杯內(nèi)部劃分為兩個區(qū)域，其中槳葉及其附近旋轉(zhuǎn)流動區(qū)域為旋轉(zhuǎn)運(yùn)動區(qū)，外部其它區(qū)域為靜止流體區(qū)域[6]。在網(wǎng)格上選擇采用四面體非結(jié)構(gòu)化體網(wǎng)格劃分形式,并對攪拌區(qū)域進(jìn)行網(wǎng)格加密處理?？紤]到計算機(jī)的運(yùn)算能力和速度,網(wǎng)格劃分大小為靜區(qū)域0.2,動區(qū)域0.1。

1.5.2 邊界參數(shù)的設(shè)定

在邊界條件的處理上，燒杯的壁面，攪拌槳葉及攪拌軸的外表面定義為WALL。內(nèi)部旋轉(zhuǎn)網(wǎng)格與外部網(wǎng)格交界面定義為INTERFACE，液面定義為SYMMETRY。

1.5.3 計算模型的方法

計算方法選用三維標(biāo)準(zhǔn)k-ε湍流模型，流動狀況為定常流，其速度壓力耦合方式采SIMPLE算法，差分格式為二階迎風(fēng)格式，收斂精度為10-4，其他參數(shù)按照默認(rèn)值設(shè)置。

2 模擬結(jié)果分析和討論

2.1 槳葉形狀對原水混凝效果的影響

數(shù)值模擬實(shí)驗主要研究轉(zhuǎn)速為60 rad/min的絮凝階段，進(jìn)入該階段前混凝劑和水體在高速旋轉(zhuǎn)下已經(jīng)充分混合均勻，大部分膠體已經(jīng)脫穩(wěn)，能否相互碰撞形成大的顆粒對沉降非常關(guān)鍵[7]。

2.1.1 速度流場的分析

圖3（a）和圖3（b）分別為方杯內(nèi)部采用平槳葉和彎葉槳時軸向截面的速度流場分布云圖。彩圖中的紅色和藍(lán)色對應(yīng)燒杯內(nèi)相對高速和低速的水體流速分布情況。由圖3可知水體在槳葉附近區(qū)域產(chǎn)生強(qiáng)烈的徑向運(yùn)動,這是因為槳葉附近的液體流速較快并且具有射流特征,這部分相對高速的流體卷吸周圍低速流體從而一起運(yùn)動，遇到燒杯杯壁分別向上方和下方流動,在槳葉的上下方都呈現(xiàn)出一個循環(huán)。

圖3 方杯不同槳葉下的軸向截面后處理速度云圖

為了比較在方杯內(nèi)距離燒杯底部不同位置高度處的兩種槳葉對水體絮凝反應(yīng)的影響，對方杯內(nèi)橫截面的速度流場進(jìn)行了速度云圖的分析。圖4為平槳葉和彎葉槳距離方杯底部高度分別為10 mm、40 mm、70 mm處橫截面的速度流場分布云圖。和平槳不同的是，在偏離槳葉附近區(qū)域時彎葉槳內(nèi)部速度下降的更快，在彎葉槳下端出現(xiàn)大量“死區(qū)”，非常容易發(fā)生堆積現(xiàn)象，影響水體絮凝效果。相比于彎葉槳，平槳很大程度的避免了這一問題，在軸向和橫向位置上流場分布都比較均勻，且速度大小適宜，有利于成熟絮凝顆粒的形成，在一定程度上又避免了大顆粒的破碎。

圖4 方杯內(nèi)距離底部不同位置高度下的速度流場分布云圖

以橫截面與自由液面高度的比值作為橫坐標(biāo)，速度等相關(guān)參數(shù)作為縱坐標(biāo)來研究反應(yīng)器內(nèi)部的整體絮凝效果[8]。圖5為方杯內(nèi)兩種槳葉的平均速度變化趨勢圖，在燒杯內(nèi)彎葉槳對燒杯內(nèi)部水體的擾動較小，這也直接減少了顆粒間的碰撞機(jī)會，同時顆粒之間的碰撞強(qiáng)度不夠會影響絮體的密實(shí)程度，在沉淀階段很可能因為尺度不夠而不能很快地下沉，并最終影響了出水水質(zhì)。

圖5 方杯內(nèi)兩種槳葉平均速度變化趨勢曲線圖

2.1.2 平均湍動能和有效能耗的分析

平均湍動能、有效能耗作為數(shù)值模擬軟件的常用的分析指標(biāo)在很多研究中得到應(yīng)用，其中湍流動能K主要用來衡量流體湍流的發(fā)展或衰退，有效能耗ε常和湍動能一起用來研究攪拌效果的好壞，因為有效能耗是總能耗中顆粒碰撞絮凝能得到的那一部分有用的能量，它主要反應(yīng)的是渦旋的尺度。按照相關(guān)絮凝理論，湍動能和有效能耗越大，絮凝效果越好。

圖6為方杯內(nèi)兩種槳葉下的平均湍動能變化趨勢圖，圖7為方杯內(nèi)兩種槳葉下的有效能耗變化趨勢圖。由圖6和圖7可知，在槳葉附近區(qū)域是產(chǎn)生有效碰撞的主要區(qū)域，平槳葉在方杯中的平均湍動能和平均有效能耗都高于彎葉槳，表明平槳整體效果優(yōu)于彎葉槳。

圖6 方杯內(nèi)兩種槳葉湍動能變化趨勢曲線圖

圖7 方杯內(nèi)兩種槳葉有效耗散率變化趨勢曲線圖

2.2 燒杯形狀對原水混凝效果的影響

2.2.1 速度流場的分析

圖8為方杯內(nèi)部采用圓杯平槳時軸向面后處理的速度云圖。圖9為圓杯內(nèi)不同位置高度處的速度云圖。和圖4比較發(fā)現(xiàn)，在低慢轉(zhuǎn)速下，圓杯反應(yīng)器由于四壁光滑，顆粒之間的速度梯度小，出現(xiàn)了“打漩”現(xiàn)象從而影響了絮凝反應(yīng)的充分發(fā)生。當(dāng)然，方杯四周尖角處有小范圍的“死區(qū)”，也會對絮凝過程產(chǎn)生一定的影響，但是方杯具有自己獨(dú)特的優(yōu)勢，四周杯壁具有阻擋流體的作用，在旋轉(zhuǎn)平槳葉帶動下產(chǎn)生共同旋轉(zhuǎn)的機(jī)會大大減弱，增加了不同轉(zhuǎn)速下顆粒之間的碰撞機(jī)會，形成大的顆粒后在沉淀階段更有利于沉降。

圖8 圓杯平槳軸向截面后處理速度云圖

圖9 圓杯內(nèi)距離底部不同位置高度下的速度流場分布云圖

2.2.2 平均湍動能和有效能耗的分析

由數(shù)值模擬方法所得的平槳條件下兩種燒杯內(nèi)的平均湍動能和有效能耗的變化分別見圖10和圖11所示。由圖10和圖11可知，方杯內(nèi)部的平均湍動能和有效能耗稍高于圓杯，也就是說方杯更能有效地利用槳葉傳遞過來的能量。

圖10 兩種燒杯湍動能變化趨勢曲線圖

圖11 兩種燒杯有效耗散率變化趨勢曲線圖

3 現(xiàn)場實(shí)測實(shí)驗結(jié)果與分析

3.1 濁度去除效果分析

圖12為在方杯內(nèi)采用不同槳葉形狀和水廠高密度沉淀池對原水濁度去除的效果，圖13為平槳葉條件下不同燒杯和水廠高密度沉淀池對原水濁度去除的效果。研究發(fā)現(xiàn)，混凝燒杯實(shí)驗在確定的攪拌工況條件和相同的PAC投加量進(jìn)行過程中，隨著原水濁度的小幅度波動，兩種燒杯進(jìn)行混凝燒杯實(shí)驗都能有效地降低濁度。一方面方杯平槳條件下的剩余濁度和水廠沉淀池的出水濁度基本吻合，能有效地代表生產(chǎn)運(yùn)行的實(shí)際情況。另一方面，同樣的方杯條件下，采用平槳葉時的剩余濁度低于采用彎葉槳時的剩余濁度，而同樣的平槳葉下，方杯中剩余濁度略低于圓杯中的剩余濁度。試驗結(jié)果表明，在該項實(shí)驗條件下方杯-平槳葉組合模擬運(yùn)行實(shí)際情況更理想，去除原水濁度效果最好，而槳葉類型對出水濁度的影響略高于燒杯形狀對其的影響。

圖12 方杯不同形狀槳葉和水廠沉淀池對濁度的去除柱狀圖

圖13 不同形狀燒杯和水廠沉淀池對濁度的去除柱狀圖

3.2 高錳酸鹽指數(shù)去除效果分析

在給水處理中高錳酸鹽指數(shù)是衡量天然有機(jī)物和人工合成有機(jī)化合物的重要依據(jù)，也是該項實(shí)驗現(xiàn)場所在水廠的重要考核指標(biāo)。方杯內(nèi)不同形狀槳葉和水廠生產(chǎn)高密度沉淀池對原水高錳酸鹽指數(shù)去除的影響見圖14所示。靜沉取樣后對水樣分析發(fā)現(xiàn)，方杯平槳組合和水廠高密度沉淀池對高錳酸鹽指數(shù)的去除非常接近，同時方杯內(nèi)部采用平槳葉的出水高錳酸鹽指數(shù)低于彎葉槳，平槳葉對高錳酸鹽指數(shù)的去除量比彎葉槳對高錳酸鹽指數(shù)的去除量高0.4～0.6 mg/L。圖15為平槳條件下不同燒杯和水廠高密度沉淀池對高錳酸鹽指數(shù)的去除效果影響，可以發(fā)現(xiàn)，在平槳條件下，方杯對高錳酸鹽指數(shù)的去除效果略好于圓杯。表明,在該項實(shí)驗條件下，方杯-平槳組合模擬實(shí)際運(yùn)行高密度沉淀池吻合度高，去除原水高錳酸鹽指數(shù)效果更好，同時槳葉類型對出水高錳酸鹽指數(shù)的影響高于燒杯形狀，因此較為良好的混凝條件更有利于高錳酸鹽指數(shù)的去除。

圖14 方杯不同形狀槳葉和水廠沉淀池對原水高錳酸鹽指數(shù)的去除柱狀圖

圖15 不同形狀燒杯和水廠沉淀池對原水高錳酸鹽指數(shù)的去除柱狀圖

4 結(jié)論

對嘉興某自來水廠水源水進(jìn)行混凝燒杯實(shí)驗的數(shù)值模擬表明，在該項實(shí)驗工況下，方杯平槳組合進(jìn)行混凝燒杯實(shí)驗可以作為確定水廠運(yùn)行高密度沉淀池混凝劑投加量的參考依據(jù)。方杯-平槳葉比其他燒杯槳葉組合方式更適合混凝燒杯實(shí)驗，其中槳葉的選擇是主要因素?，F(xiàn)場的實(shí)驗數(shù)據(jù)表明使用方杯和平槳葉時對原水濁度和高錳酸鹽指數(shù)的去除效果好于其他燒杯槳葉的組合情況，因此建議水廠采用方杯-平槳葉進(jìn)行混凝燒杯實(shí)驗。

混凝燒杯實(shí)驗的數(shù)值模擬和現(xiàn)場實(shí)驗的實(shí)測結(jié)果基本一致，這表明CFD數(shù)值模擬軟件fluent6.3從一定程度上能解決“黑箱”內(nèi)部流態(tài)問題，為同類問題的深入研究提供一種研究思路和方向。

[1]邱振華.杯罐實(shí)驗的優(yōu)化設(shè)計及在飲用水生產(chǎn)中的應(yīng)用[J].凈水技術(shù),2003,22(1):28-31.

[2]朱海濤,周曉龍,劉宏遠(yuǎn),孫海平,許衛(wèi)初,張剛.中置式高密度沉淀池的改造與優(yōu)化運(yùn)行[J].中國給水排水,2014，(4):95-99.

[3]范蘢,施漢昌,徐農(nóng),施慧明.污水處理反應(yīng)器的計算流體力學(xué)[M].北京：中國建筑工業(yè)出版社,2012.

[4]Bertrand J,Poux M,Aubin J.European Federation of Chemical Engineering Working Party on Mixing:Development,State of the Art and Future of the Scientific Field[J].Chemical Engineering Research and Design,2004,82(12):1575-1579.

[5]楊小林,楊開明,王華,李太富.基于FLUENT的折板絮凝池三維渦旋水流數(shù)值模擬[J].給水排水,2009，(8):21-23.

[6]王福軍.計算流體動力學(xué)分析—CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社,2005:13-22.

[7]Jarvis P,Jefferson B,Gregory J,Parsons SA.A Review of Floc Strength and Breakage[J].Water Research,2005,39:3121-3137.

[8]尉鵬翔.攪拌反應(yīng)器內(nèi)三維流場對絮凝效果影響研究[J].人民黃河,2011，(5):43-45.