陶洪飛，朱玲玲，馬英杰，滕曉靜，鄭文強

(新疆農(nóng)業(yè)大學(xué)水利與土木工程學(xué)院，烏魯木齊 830052)

2010年中央1號文件《關(guān)于加快水利改革發(fā)展的決定》提出加大水利建設(shè)，實施最嚴(yán)格的水資源管理制度，在新疆農(nóng)業(yè)用水占社會經(jīng)濟(jì)總用水的96%，而對GDP的貢獻(xiàn)率卻只有17.1%，而農(nóng)業(yè)灌溉僅有39.2%采用了高效節(jié)水灌溉，對于水資源短缺的新疆來說農(nóng)業(yè)節(jié)水潛力還很大[1]。根據(jù)《新疆高效節(jié)水建設(shè)方案》及《新疆農(nóng)業(yè)節(jié)水建設(shè)發(fā)展規(guī)劃》確定了每年完成農(nóng)業(yè)高效節(jié)水面積20萬hm2以上，至2020年新疆農(nóng)業(yè)高效節(jié)水面積累計達(dá)到286.667萬hm2以上的推廣目標(biāo)。目前，新疆已建的有關(guān)農(nóng)業(yè)高效節(jié)水工程，灌溉水源絕大部分為地下水，但許多地區(qū)的地下水超采嚴(yán)重，所以未來選擇地表水作為灌溉水源已成為新建微灌工程的最優(yōu)選擇[2]。新疆河水多為山溪性河流，泥沙含量高，易造成灌水器(如滴頭、微噴頭等)堵塞，目前，采用的是沉沙池和組合式過濾器的方法去除微灌用水中的泥沙，以防止灌水器堵塞，從而保障工程的正常運行[3,4]。過濾器是微灌系統(tǒng)中的關(guān)鍵性設(shè)備之一，其中網(wǎng)式過濾器相對于砂石過濾器和疊片過濾器而言，具有容易操作、更方便清洗及對無機污染物有較高的去除效率等優(yōu)勢[5]。目前，大部分學(xué)者[6-11]開展了有關(guān)網(wǎng)式過濾器系列產(chǎn)品的水力性能、去除效率及結(jié)構(gòu)優(yōu)化等方面的研究。WU[12]根據(jù)網(wǎng)式過濾器的結(jié)構(gòu)參數(shù)和大量的水頭損失試驗數(shù)據(jù)，得到了一個水頭損失計算的三維模型。

通過查閱相關(guān)文獻(xiàn)，發(fā)現(xiàn)對網(wǎng)式過濾器內(nèi)部流場的研究很少，分析原因是網(wǎng)式過濾器的內(nèi)部結(jié)構(gòu)較為復(fù)雜，很難通過現(xiàn)有的測試手段(如PIV或LDP)來測量速度場，但隨著計算機的廣泛應(yīng)用以及計算機容量的增大，仿真模擬流動已具有很高的可靠性和準(zhǔn)確度。故選擇正確的數(shù)學(xué)模型對網(wǎng)式過濾器中的流場進(jìn)行數(shù)值模擬不失為一種有效的方法。因此，筆者采用Fluent軟件開展了網(wǎng)式過濾器過濾過程中的清水流場模擬，并將計算結(jié)果同物理試驗結(jié)果對比，從而選擇最適合模擬網(wǎng)式過濾器內(nèi)部流場的數(shù)學(xué)模型，并在此基礎(chǔ)上分析了該過濾器的清水流場，計算結(jié)果為該過濾器的后續(xù)研究工作奠定了基礎(chǔ)，同時也為研究網(wǎng)式過濾器提供了一種新的研究手段。

1 網(wǎng)式過濾器的介紹

由圖1可知，網(wǎng)式過濾器由罐體、濾網(wǎng)、控制器等主要部分組成，其關(guān)鍵參數(shù)包括罐體長度1.031 m，直徑0.254 m；濾網(wǎng)長度0.961 m，直徑為0.2 m；進(jìn)、出口直徑均為0.2 m；排污口直徑0.05 m。

圖1 網(wǎng)式過濾器結(jié)構(gòu)示意圖(單位：m)Fig.1 Structural diagram of screen filter

過濾器的工作原理描述如下：當(dāng)灌溉水從進(jìn)水管流入網(wǎng)式過濾器中，所有大于網(wǎng)孔孔徑尺寸的污物都會滯留在濾網(wǎng)上，而過濾水則流向出水管；隨著污物的累積，水頭損失逐漸增大，當(dāng)水頭損失達(dá)到預(yù)設(shè)壓差值時，控制器自動開啟排污閥進(jìn)行排污，排污的同時過濾依然進(jìn)行，滿足灌溉需水要求，當(dāng)達(dá)到預(yù)設(shè)的排污時間時，排污閥自動關(guān)閉，沖洗結(jié)束，繼續(xù)過濾。

2 網(wǎng)式過濾器的CFD模擬

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

采用Fluent軟件中的多孔介質(zhì)階躍(Porous Jump)模型來模擬濾網(wǎng)，多孔介質(zhì)相關(guān)參數(shù)表達(dá)式如下：

(3)

式中：C1為阻力系數(shù)，m-2；C2為慣性損失系數(shù)，m-1；α為滲透率，m2；D為濾網(wǎng)孔直徑，mm；ε為孔隙比(也叫“篩分面積百分比”)，%。

根據(jù)物理試驗中采用的濾網(wǎng)參數(shù)：濾網(wǎng)孔直徑D=0.18 mm，孔隙比ε=38%，絲徑d=0.112 mm。根據(jù)給定的參數(shù)值，可計算得到：C1=3.24×1010m-2，C2=219 703.2 m-1，α=3.09×10-11m2。數(shù)值計算時考慮過濾網(wǎng)的厚度，其值為0.001 m。

分別選用Fluent6.3提供的Standardk-ε、RNGk-ε及Realizablek-ε為數(shù)值計算的湍流模型，與多孔介質(zhì)階躍模型耦合求解。這3種模型的湍動能k和耗散率ε的方程見文獻(xiàn)[13]。

2.2 坐標(biāo)建立

圖2為網(wǎng)式過濾器的三維坐標(biāo)。從圖2中可知該過濾器的長度方向為X坐標(biāo)，寬度方向為Z坐標(biāo)，高度方向為Y坐標(biāo)。為分析該網(wǎng)式過濾器中的流場分布情況，分別截取X=0.51 m、Y=0 m及Z=0 m 3個典型代表面進(jìn)行分析研究。并定義遠(yuǎn)離出水口一側(cè)的濾網(wǎng)為上濾網(wǎng)，靠近出水口一側(cè)的濾網(wǎng)為下濾網(wǎng)；分別選取Y=±0.114、±0.1、±0.086斷面與Z=0斷面相交的6條線：Y1=0.114(上濾網(wǎng)外側(cè))、Y2=0.1(上濾網(wǎng))、Y3=0.086(上濾網(wǎng)內(nèi)側(cè))、Y4=-0.086(下濾網(wǎng)內(nèi)側(cè))、Y5=-0.1(下濾網(wǎng))、Y6=-0.114(下濾網(wǎng)外側(cè))來分析濾網(wǎng)及其內(nèi)、外側(cè)的流速。

圖2 三維坐標(biāo)Fig.2 Three-dimensional coordinate

2.3 網(wǎng)格劃分及計算方法

為減少數(shù)值模擬的計算工作量和保證計算精度，對過濾器的計算網(wǎng)格進(jìn)行了優(yōu)化。將過濾器的網(wǎng)格進(jìn)行分塊劃分，進(jìn)、出口因為結(jié)構(gòu)簡單，采用六面體結(jié)構(gòu)性網(wǎng)格，其余部分采用四面體非結(jié)構(gòu)網(wǎng)格，網(wǎng)格總數(shù)為123 252個。

對網(wǎng)式過濾器的數(shù)值計算采用非定常的壓力基隱式算法，研究的計算區(qū)域和控制方程的離散均采用有限體積法。采用Fluent6.3提供的SIMPLE算法來計算壓力與速度的耦合方程，差分格式采用精度更高的二階迎風(fēng)格式。殘差標(biāo)準(zhǔn)為1.0×10-3，迭代時間步長為0.000 1 s。

2.4 邊界條件

(1)過濾器進(jìn)水管中運動的水流，可近似視為均勻流，取進(jìn)流方向與X軸正向一致，故設(shè)置進(jìn)水口的邊界條件為速度進(jìn)口：U=u，V=0，W=0。根據(jù)網(wǎng)式過濾器物理模型試驗情況，進(jìn)流量為352.1 m3/h，則根據(jù)流速定義可得該工況下u= 3.11 m/s。

(2)網(wǎng)式過濾器出水口的邊界條件設(shè)置成壓力出口，壓力大小為200 kPa。

(3)網(wǎng)式過濾器的罐體、管道內(nèi)壁及排污管均按固壁定律進(jìn)行處理，采用無滑移邊界條件，即U=V=W=0。

3 最優(yōu)數(shù)學(xué)模型的確定

3.1 水頭損失對比

對于過濾器而言，我們最為關(guān)注的是其水頭損失，因為其值的大小是評價過濾器性能的一個重要標(biāo)準(zhǔn)。表1為物理試驗結(jié)果與數(shù)值模擬結(jié)果的水頭損失對比情況，從表1中可知，物理試驗中網(wǎng)式過濾器的水頭損失為3 m，而紊流模型中的Standardk-ε模型、RNGk-ε模型、Realizablek-ε模型的計算水頭損失分別為3.91、3.18、3.12 m，與物理試驗結(jié)果的相對誤差分別為30.3%、6.0%、4.0%，因此，Realizablek-ε模型與多孔介質(zhì)階躍模型耦合的計算結(jié)果同物理試驗結(jié)果吻合最好。同時也驗證了Realizablek-ε模型可以模擬高應(yīng)變率及彎曲程度較大的流動[14,15]。

表1 結(jié)果對比Tab.1 Results comparison

3.2 計算耗時對比

網(wǎng)式過濾器的數(shù)值計算是在同一臺計算機上進(jìn)行，計算機內(nèi)存都占用127 MB，CPU都占用25%，而當(dāng)進(jìn)、出口流量平衡，即完全收斂時，從表2可知，RNGk-ε模型計算用時272 min> Standardk-ε模型計算用時271 min> Realizablek-ε模型計算用時269 mim，但總的來說，三者計算時間相差不大，最長與最短的計算相差僅為3 min。

表2 3種模型的計算時間對比Tab.2 Time comparison of three models

通過水頭損失和計算耗時對比，最終選擇Realizablek-ε模型與多孔介質(zhì)階躍模型來模擬網(wǎng)式過濾中的流場。

4 流場分析與討論

4.1 速度場對比

圖3為網(wǎng)式過濾器中不同斷面的速度矢量分布圖。從圖3(b)和圖3(c)中可知，水流從進(jìn)水管處均勻地流入網(wǎng)式過濾器中，當(dāng)水流進(jìn)入罐體時，水流的速度分布特性發(fā)生改變，一部分水流進(jìn)入上濾網(wǎng)，通過濾網(wǎng)孔進(jìn)行過濾；一部分水流進(jìn)入下濾網(wǎng)進(jìn)行過濾，當(dāng)突然遇到出水口邊界時，流速分布發(fā)生了巨大的變化，水流由沿X軸正向運動快速轉(zhuǎn)變?yōu)檠爻鏊苡覀?cè)斜向下運動，從而造成出水管左側(cè)的部分區(qū)域產(chǎn)生漩渦，水流運行33 cm左右的距離后，水流沿出水管均勻地向下運動；而X=0.61 m處以后的水流速度大部分呈“點狀”分布，說明流速很小。從圖3(a)中可知水流在出水管中的流速是關(guān)于Z=0截面對稱分布的，靠近管壁內(nèi)側(cè)的水流流速較大，在運動至距出水口10 cm處時，水流流速會減小。

圖3 不同斷面的速度矢量圖Fig.3 Velocity vector diagram of the different section

圖4為網(wǎng)式過濾器中不同斷面的速度云圖。從圖4(b)和圖4(c)中可以看出進(jìn)水管中的水流流速分布比較均勻，為3.0 m/s左右；當(dāng)水流運動至罐體時，過水?dāng)嗝婷娣e增大，水流流速減小，水流繼續(xù)運動一段路程，過水?dāng)嗝婷娣e沒變，但流速卻沿X方向逐漸減小，最終達(dá)到一個恒定值為0.5 m/s左右。這是因為受到排污口固體邊界條件的影響，當(dāng)水流運動至該處時，水流會發(fā)生反向流動，與沿X軸正向運動的水流發(fā)生碰撞、摩擦和混摻，從而消耗大量的能量，造成水流流速減小。由于出水口邊界條件影響，出水口中間處的水流流速依然較大，為5.0 m/s左右，而左、右側(cè)漩渦處的流速較小，為0.5 m/s左右。從圖4(a)可知出水管中的水流流速關(guān)于Z=0斷面對稱分布，水流在出水管Y=-0.12 m左右時的流速最大，5.0 m/s左右。

圖4 不同斷面的速度云圖Fig.4 Velocity contours diagram of the different section

4.2 壓強場對比

圖5為網(wǎng)式過濾器中不同斷面的壓強場。從圖5(a)可知出水口正上方罐體中的水流壓強由中心向四周逐漸減小，與速度流場分布規(guī)律相反，中心處壓強最大，可達(dá)230 kPa左右，濾網(wǎng)附近最小的壓強為210 kPa左右，兩者相差約20 kPa。從圖5(b)可知，X=0～0.2 m處，壓強分布比較均勻，X=0.2～0.75 m處，上、下濾網(wǎng)的內(nèi)外壓差較大，這是因為水流進(jìn)入罐體后會以一個相對較大的流速分別流向上、下濾網(wǎng)，并經(jīng)過濾網(wǎng)孔完成過濾，從而造成濾網(wǎng)的內(nèi)外壓差；X>0.75 m，上、下濾網(wǎng)的內(nèi)外壓差較小，壓強分布較均勻。從圖5(c)可知進(jìn)水管處的水流壓強明顯大于出水管處的水流壓強。由于受到罐體中上、下濾網(wǎng)及出水口邊界條件的影響，水流從進(jìn)水管流至出水管過程中會發(fā)生局部水頭損失和沿程水頭損失。故當(dāng)位置水頭和流速水頭一定時，進(jìn)水管處的壓強大于出水管處的水流壓強，與實際液體恒定總流的能量方程相符合。

圖5 不同斷面的壓強場Fig.5 Pressure flow of the different section

5 流速討論

圖6為水流流速沿X軸的變化曲線。從圖6(a)中可得出2點：①上濾網(wǎng)及上濾網(wǎng)外側(cè)的水流流速沿X軸的變化包含兩個階段：流速迅速增大階段。X=0.2～0.25 m時，上濾網(wǎng)及上濾網(wǎng)外側(cè)的水流流速沿X軸迅速增加，但上濾網(wǎng)的水流流速增加幅度明顯大于上濾網(wǎng)外側(cè)，上濾網(wǎng)的流速最大值達(dá)2.2 m/s左右，約為上濾網(wǎng)外側(cè)流速最大值的兩倍。流速逐漸減小階段。X=0.25～1.23 m時，上濾網(wǎng)及上濾網(wǎng)外側(cè)的水流流速逐漸減小，雖然減小幅度不同，但減小趨勢大體相同。②上側(cè)濾網(wǎng)內(nèi)側(cè)的水流流速沿著X軸逐漸減小，且水流流速比上濾網(wǎng)外側(cè)的水流流速大很多，但兩者的流速差卻沿著X軸逐漸減小，兩者的最大流速差約為2 m/s；當(dāng)X=1.05～1.191 m時，上濾網(wǎng)內(nèi)、外側(cè)的水流流速大體相同，流速差幾乎為0。從圖6(b)中可知，下濾網(wǎng)及其內(nèi)外側(cè)的水流流速沿X軸方向的變化規(guī)律，與上濾網(wǎng)及其內(nèi)外側(cè)的有所不同，主要是因為出水口邊界條件的影響。

圖6 水流流速沿X軸的變化曲線Fig.6 Change curve of water velocity along X axis

綜上所述，濾網(wǎng)內(nèi)、外側(cè)的水流流速沿X軸方向的分布很不均勻，從而影響網(wǎng)式過濾器的過濾效率，特別是網(wǎng)式過濾器尾部附近的濾網(wǎng)，內(nèi)、外側(cè)流速相差較小，水的滲透率就較小，導(dǎo)致過濾效率較低，如何提高網(wǎng)式過濾器的過濾效率值得進(jìn)一步研究，可通過改變出水管出水口位置解決水流流速分布不均勻的問題。

6 結(jié) 語

利用Fluent6.3軟件對網(wǎng)式過濾器的內(nèi)部流場進(jìn)行數(shù)值模擬，主要得到以下結(jié)論。

(1)Realizablek-ε模型與多孔介質(zhì)階躍模型能更好地模擬該網(wǎng)式過濾器的內(nèi)部流場，這為今后該過濾器的結(jié)構(gòu)優(yōu)化、沖洗過程及水沙兩相流流場的數(shù)值模擬奠定了基礎(chǔ)。

(2)得到了水流從進(jìn)水管進(jìn)入網(wǎng)式過濾器并流出出水管的運動狀態(tài)，指出排污口固體邊界和出水口邊界條件對該過濾器的速度流場和壓強場分布規(guī)律影響很大。

(3)濾網(wǎng)內(nèi)外側(cè)的水流流速沿X軸的分布不均勻，但由于出水口固壁邊界條件條件的影響，上濾網(wǎng)水流的流速沿X軸的分布規(guī)律不同于下濾網(wǎng)，降低了網(wǎng)式過濾器的過濾效率；同時濾網(wǎng)內(nèi)外側(cè)較大的壓強差，容易導(dǎo)致濾網(wǎng)變形，直至損壞，故減小流速和壓強的不均勻分布有待進(jìn)一步研究。

[1] 蘇 薈.新疆農(nóng)業(yè)高效節(jié)水灌溉技術(shù)選擇研究[D]. 新疆石河子：石河子大學(xué)，2013.

[2] 孫 娟.河水滴灌重力沉沙過濾池的設(shè)計與應(yīng)用推廣[J].節(jié)水灌溉，2014,(1):60-64.

[3] 陶洪飛，邱秀云，李 巧，等.不同鰓片間距下的分離鰓內(nèi)部流場三維數(shù)值模擬[J].農(nóng)業(yè)機械學(xué)報，2014，45(6):183-189.

[4] 陶洪飛，邱秀云，趙麗娜，等.水沙分離鰓內(nèi)部流場的數(shù)值模擬[J].農(nóng)業(yè)工程學(xué)報，2013，29(17)：38-46.

[5] Tajrishy MA,Hills DJ, Tchobanoglous G. Pretreatment of secondary effluent for drip irrigation[J]. Journal of Irrigation and Drainage,1994,120(4): 716-731.

[6] 劉煥芳，劉 飛，谷趁趁，等.自清洗網(wǎng)式過濾器水力性能試驗[J]. 排灌機械工程學(xué)報，2012，(2)：203-208.

[7] Avner Adin, Giora Alon. Mechanisms and process parameters of filter screens[J]. Irrigation and Drainage Engineering, 1996,293(4):112-115.

[8] Thokal R T, Raghavendra A G, Suresh N M, et al. Effect of sand particle size and concentration on performance of screen filter in trickle irrigation[J]. An-nals of Arid Zone, 2004,43(1):65-71.

[9] 宗全利，劉 飛，劉煥芳，等.滴灌用自清洗網(wǎng)式過濾器排污壓差計算方法[J].農(nóng)業(yè)機械學(xué)報，2012，43(11)：107-112.

[10] 劉煥芳，鄭鐵剛，劉 飛，等．自吸網(wǎng)式過濾器過濾時間與自清洗時間變化規(guī)律分析[J]．農(nóng)業(yè)機械學(xué)報，2010，41(7)：80-83.

[11] 劉 飛，劉煥芳，宗全利，等.自清洗網(wǎng)式過濾器水頭損失和排污時間研究[J].農(nóng)業(yè)機械報，2013，44(5)：127-134.

[12] Wenyong WU, Wei CHEN, Honglu LIU, et al. A new model for head loss assessment of screen filters developed with dimensional analysis in drip irrigation systems [J]. Irrigation and Drainage, 2014 (63):523-531.

[13] 白兆亮，李 琳. 有壓管道中孔板相對間距對局部阻力系數(shù)的影響及其機理研究[J].水電能源科學(xué)，2015，33(1)：177-182.

[14] 于 勇. FLUENT入門與進(jìn)階教程[M].北京：北京理工大學(xué)出版社，2010.

[15] 江 帆，黃 鵬. FLUENT高級應(yīng)用與實例分析[M].北京：清華大學(xué)出版社，2010.