陶洪飛，滕曉靜，馬英杰，趙經(jīng)華，洪 明，馬 亮，付秋萍

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

過(guò)濾器是微灌系統(tǒng)中的重要組成部分，網(wǎng)式過(guò)濾器具有結(jié)構(gòu)簡(jiǎn)單、造價(jià)低廉、體積小等特點(diǎn)，一般作為最末級(jí)過(guò)濾器，其過(guò)濾性能好壞直接影響到下級(jí)輸配水管網(wǎng)的正常運(yùn)作和使用壽命。

國(guó)內(nèi)對(duì)過(guò)濾器的研究雖起步較晚，但近些年已有一些系統(tǒng)的研究。石河子大學(xué)的宗全利[1,2]、駱秀萍[3,4]、劉飛[5]等人針對(duì)自清洗網(wǎng)式過(guò)濾器的水力性能進(jìn)行了大量研究。山西省運(yùn)城市水利科學(xué)研究所研發(fā)了一種離心篩網(wǎng)一體式過(guò)濾器[6]，將篩網(wǎng)置于離心過(guò)濾器腹腔中，并對(duì)該過(guò)濾器的水力性能[7]、砂礫運(yùn)移規(guī)律[8]進(jìn)行研究。在國(guó)外，Capra[9]研究了不同類(lèi)型的過(guò)濾器對(duì)廢水的過(guò)濾能力，結(jié)果表明網(wǎng)式過(guò)濾器和疊片過(guò)濾器對(duì)過(guò)濾有機(jī)污物的能力較差。Yurdem H[10]、Duran-Ros M[11]等人對(duì)過(guò)濾器的水利性能和效率進(jìn)行了深入研究。

過(guò)濾器傳統(tǒng)的研發(fā)過(guò)程是需要先制造出網(wǎng)式過(guò)濾器的實(shí)體模型，再進(jìn)行水力性能和泥沙處理能力的研究，該研發(fā)方式不僅周期長(zhǎng)，而且成本高，也無(wú)法探尋過(guò)濾器內(nèi)部的流場(chǎng)狀況。隨著計(jì)算機(jī)普及以及計(jì)算軟件的發(fā)展，關(guān)于網(wǎng)式過(guò)濾器數(shù)值模擬的研究也開(kāi)始出現(xiàn)[12,13]。王新坤[14]用多孔介質(zhì)模型模擬濾網(wǎng)，將數(shù)值計(jì)算結(jié)果與試驗(yàn)結(jié)果相比較，發(fā)現(xiàn)吻合程度較高，證實(shí)了多孔介質(zhì)模型的可行性。目前，已有不少學(xué)者將多孔介質(zhì)模型應(yīng)用到其他型號(hào)的過(guò)濾器上，以探究過(guò)濾器內(nèi)部流場(chǎng)狀況并進(jìn)行結(jié)構(gòu)優(yōu)化[15-18]。相對(duì)于傳統(tǒng)過(guò)濾器的研發(fā)方式，將CFD技術(shù)應(yīng)用于過(guò)濾器的結(jié)構(gòu)優(yōu)化有著無(wú)可比擬的優(yōu)勢(shì)。

本文擬利用Fluent軟件模擬網(wǎng)式過(guò)濾器在過(guò)濾狀態(tài)下的清水流場(chǎng)，探究出水管與罐體之間的角度對(duì)過(guò)濾器內(nèi)部流場(chǎng)的影響，從而為網(wǎng)式過(guò)濾器的優(yōu)化設(shè)計(jì)提供理論依據(jù)。

1 全自動(dòng)網(wǎng)式過(guò)濾器的介紹

本文所研究的全自動(dòng)網(wǎng)式過(guò)濾器(下文簡(jiǎn)稱(chēng)“過(guò)濾器”)具有便于維護(hù)、工作流量大、水頭損失小、自動(dòng)反沖洗等優(yōu)點(diǎn)。該過(guò)濾器由罐體、不銹鋼濾網(wǎng)、自動(dòng)控制器等部件組成，圖1為過(guò)濾器的整體示意圖。過(guò)濾器主要尺寸如下：進(jìn)水管、出水管直徑均為0.2 m；排污管直徑0.05 m；濾網(wǎng)長(zhǎng)0.961 m，直徑0.2 m；罐體長(zhǎng)1.031 m，直徑0.257 m。為過(guò)濾器建立坐標(biāo)系，坐標(biāo)原點(diǎn)在進(jìn)水管中心處，濾網(wǎng)軸線、進(jìn)水管軸線、排污管軸線均與X軸重合，出水管軸線平行于Y軸，出水管與罐體之間的角度α=90°。取Z=0 m斷面為過(guò)濾器橫斷面。

圖1 全自動(dòng)網(wǎng)式過(guò)濾器示意圖Fig.1 Structural diagram of automatic screen filter

過(guò)濾器過(guò)濾時(shí)，排污閥處于關(guān)閉狀態(tài)，粒徑大于濾網(wǎng)孔徑的固體雜質(zhì)會(huì)被攔截在濾網(wǎng)內(nèi)部，而凈水在濾網(wǎng)外部。進(jìn)水管口、出水管口裝有壓力傳感器，當(dāng)進(jìn)、出水口壓降達(dá)到預(yù)設(shè)壓差時(shí)，控制器控制電磁閥打開(kāi)排污閥，此時(shí)固體雜質(zhì)從排污管排出，排污時(shí)間達(dá)到預(yù)設(shè)時(shí)間后，排污閥自動(dòng)關(guān)閉，直沖洗結(jié)束。

2 全自動(dòng)網(wǎng)式過(guò)濾器的數(shù)值模擬

2.1 三維建模

先用Gambit軟件對(duì)過(guò)濾器進(jìn)行前期建模及網(wǎng)格劃分，再用Fluent軟件對(duì)過(guò)濾器過(guò)濾狀態(tài)進(jìn)行數(shù)值模擬，最后利用Tecplot軟件進(jìn)行后期數(shù)據(jù)可視化處理。圖2是過(guò)濾器的三維圖，圖中深顏色的部分為濾網(wǎng)。

圖2 三維圖Fig.2 Three-dimension graph

2.2 網(wǎng)格劃分

由于過(guò)濾器模型為三維模型，故采用體網(wǎng)格劃分更加便捷。其中進(jìn)水管和排污管由于結(jié)構(gòu)簡(jiǎn)單，采用六面體結(jié)構(gòu)性網(wǎng)格；其余部分由于結(jié)構(gòu)復(fù)雜，同時(shí)為了避免在Fluent軟件中檢查時(shí)網(wǎng)格出現(xiàn)負(fù)體積，采用四面體非結(jié)構(gòu)網(wǎng)格。網(wǎng)格是后續(xù)計(jì)算的基礎(chǔ)，網(wǎng)格的數(shù)量關(guān)系到計(jì)算量的大小和計(jì)算結(jié)果的準(zhǔn)確性，為了探究網(wǎng)格節(jié)點(diǎn)間距(網(wǎng)格數(shù)量取決于網(wǎng)格節(jié)點(diǎn)間距)對(duì)計(jì)算結(jié)果的影響，針對(duì)本模型，劃分了三種不同數(shù)量的網(wǎng)格，具體節(jié)點(diǎn)間距、網(wǎng)格數(shù)量如表1所示。

表1 網(wǎng)格劃分Tab.1 Mesh generation

2.3 計(jì)算模型、邊界條件及求解方法

數(shù)值模采用非定常流二階隱式格式，以提高計(jì)算精度；湍流模型采用Realizablek-ε模型。

求解方法采用基于壓力-速度修正算法(SIMPLE)；差分格式采用精度更高的二階迎風(fēng)格式；殘差標(biāo)準(zhǔn)設(shè)置為1.0×10-3；為了盡快收斂，迭代時(shí)間步長(zhǎng)為1.0×10-4s。

邊界條件設(shè)置如下：①過(guò)濾器罐體按固壁處理。②本文僅模擬過(guò)濾狀態(tài)下的流場(chǎng)，此時(shí)排污口關(guān)閉，故排污管出口視為固壁。③操作壓力默認(rèn)為101 325 Pa，故以下討論的壓強(qiáng)均為相對(duì)壓強(qiáng)。④進(jìn)口、出口的條件同物理實(shí)驗(yàn)，即進(jìn)口為速度進(jìn)口，流速為1.769 m/s(此時(shí)流量為200 m3/h)。出水管出口為壓力出口，出口壓強(qiáng)為263 kPa。⑤濾網(wǎng)采用多孔階躍模型，即將濾網(wǎng)所在的面設(shè)為Porous Jump。有關(guān)多孔階躍模型的具體設(shè)置見(jiàn)2.4節(jié)。

2.4 多孔階躍模型

物理實(shí)驗(yàn)中，濾網(wǎng)孔徑為0.18 mm，如此微小的網(wǎng)孔，若按照實(shí)際狀況劃分網(wǎng)格，必定會(huì)十分麻煩，且計(jì)算量巨大。借鑒其他學(xué)者對(duì)多孔階躍模型的應(yīng)用[14,16,19]，可以通過(guò)孔隙率、滲透率等參數(shù)體現(xiàn)濾網(wǎng)的阻力特性。將濾網(wǎng)所在面定義為多孔階躍面，可以大大簡(jiǎn)化問(wèn)題。多孔階躍面(Porous Jump)是多孔介質(zhì)(Porous Media)模型的一維簡(jiǎn)化，用以模擬速度-壓降特性已知的薄膜。多孔介質(zhì)模型定義在流體區(qū)域上，而多孔階躍面定義在界面區(qū)上，相比多孔介質(zhì)模型，多孔階躍面更易收斂，更適合模擬濾網(wǎng)。

多孔介質(zhì)模型通過(guò)在動(dòng)量方程中增加黏性阻力項(xiàng)和慣性損失項(xiàng)這兩個(gè)源項(xiàng)來(lái)模擬計(jì)算域中多孔性材料對(duì)流體的流動(dòng)阻力。有關(guān)這兩個(gè)源項(xiàng)可以參考文獻(xiàn)[20]。ANSYS FLUENT參考手冊(cè)給出了各方向黏性阻力系數(shù)和慣性損失系數(shù)計(jì)算公式：

(3)

式中：C1為阻力系數(shù)，m-2；k為滲透率，m2；C2為慣性損失系數(shù)，m-1；D為濾網(wǎng)孔直徑，mm；ε為孔隙比，%。

為了便于比較，濾網(wǎng)有關(guān)參數(shù)設(shè)置與物理實(shí)驗(yàn)中所用濾網(wǎng)實(shí)際參數(shù)一致：濾網(wǎng)孔直徑D=0.18 mm；孔隙比ε=38 %；絲徑d=0.112 mm；濾網(wǎng)厚度取0.001 m。根據(jù)以上數(shù)據(jù)，可計(jì)算得到k=3.083 3×10-11m2，C1=3.24×1010m-2，C2=219 703.2 m-1。

2.5 模型驗(yàn)證

為了更好地比較，規(guī)定過(guò)濾器進(jìn)口壓強(qiáng)與出口壓強(qiáng)之差為該過(guò)濾器的壓降。將三種不同網(wǎng)格的計(jì)算結(jié)果同物理實(shí)驗(yàn)結(jié)果相比較。對(duì)比結(jié)果如表2所示。由表2可知，網(wǎng)格數(shù)量對(duì)計(jì)算結(jié)果有明顯的影響。網(wǎng)格數(shù)越多，即網(wǎng)格密度越大，計(jì)算結(jié)果越準(zhǔn)確。當(dāng)網(wǎng)格數(shù)量為69 189個(gè)時(shí)，計(jì)算得到的壓降與實(shí)際壓降誤差為12%，當(dāng)網(wǎng)格數(shù)量為123 887個(gè)時(shí)，相對(duì)誤差為6.77%。當(dāng)網(wǎng)格數(shù)量為168 681時(shí)，相對(duì)誤差僅為5.06%。當(dāng)網(wǎng)格數(shù)量多到一定程度時(shí)，再增加網(wǎng)格數(shù)量，對(duì)提高計(jì)算精確度幫助不大。對(duì)比工況2、3，網(wǎng)格個(gè)數(shù)多了44 794個(gè)，相對(duì)誤差卻只減小了1.71%，可以看出工況2的網(wǎng)格密度已經(jīng)足夠。由此得出結(jié)論：工況1的網(wǎng)格數(shù)量不夠，計(jì)算誤差太大，不宜采用。工況2、3計(jì)算誤差都在可接受范圍內(nèi)，但工況3網(wǎng)格數(shù)量太多，計(jì)算量太大，也不宜采用。故本模型采用工況2的網(wǎng)格節(jié)點(diǎn)間距較為合理。

表2 結(jié)果對(duì)比Tab.2 Results comparison

3 不同出水管角度的數(shù)值模擬

3.1 原模型流場(chǎng)分析

這里僅分析工況2(網(wǎng)格數(shù)量為123 887個(gè))的計(jì)算結(jié)果。圖3是過(guò)濾器橫斷面的速度云圖、矢量圖和壓強(qiáng)云圖。從圖3(a)中可以看出水流均勻地從進(jìn)水管進(jìn)入濾網(wǎng)，一部分水流直接穿過(guò)濾網(wǎng)，流向出水管。而另一部分水流在濾網(wǎng)內(nèi)部繼續(xù)前行，與過(guò)濾器尾部折返回來(lái)的水流相互混摻，最后沿反方向流向出水管。進(jìn)水管流速為設(shè)定的1.769 m/s，濾網(wǎng)進(jìn)口段流速約為1.4 m/s，沿著濾網(wǎng)方向，流速逐漸變小，過(guò)濾器尾部流速很小，接近靜止。出水管段Y字型高流速區(qū)是整個(gè)過(guò)濾器流速最大的區(qū)域，而就在高流速區(qū)域左右兩側(cè)各有一個(gè)低流速的漩渦區(qū)域，這是因?yàn)槌鏊艽怪惫摅w，結(jié)構(gòu)突變，水流在這塊區(qū)域混摻和碰撞，從而消耗能量，流速減小。從速度矢量圖上可以看出，出水管左側(cè)流來(lái)的水流流速大于右側(cè)的水流流速，且左側(cè)漩渦區(qū)也大于右側(cè)的漩渦區(qū)。進(jìn)入出水管的水流流動(dòng)方向由近似水平逐漸轉(zhuǎn)變?yōu)榇怪?，在這個(gè)過(guò)程中，水流沖擊出水管右側(cè)邊壁，造成能量損失。

圖3 流場(chǎng)分析Fig.3 Flow field analysis

從圖3(b)中可以看出，進(jìn)水管、濾網(wǎng)、排污管內(nèi)部壓強(qiáng)分布均勻，且壓強(qiáng)較大。以濾網(wǎng)為分界面，濾網(wǎng)外部的壓強(qiáng)明顯小于濾網(wǎng)內(nèi)部，這正是由于多孔階躍面產(chǎn)生的壓降，反映到實(shí)際情況上，就是液體穿過(guò)濾網(wǎng)所產(chǎn)生的壓降。在濾網(wǎng)前部壓降約為6 kPa，濾網(wǎng)中部(即出水管上部)壓降約為8 kPa，濾網(wǎng)尾部壓降約為2 kPa。漩渦區(qū)域壓力低于周?chē)s為2 kPa。整個(gè)過(guò)濾器壓強(qiáng)損失約為10 kPa，可見(jiàn)絕大部分壓強(qiáng)損失在濾網(wǎng)附近。

3.2 不同出水管角度模型

由于這種過(guò)濾器具有實(shí)物的只有一種----出水管垂直罐體，即夾角α=90°。為研究其他夾角對(duì)過(guò)濾器內(nèi)部流場(chǎng)的影響，本文采用相同的數(shù)值計(jì)算方法，對(duì)不同出水管夾角的過(guò)濾器進(jìn)行數(shù)值模擬。在保證濾網(wǎng)面積、進(jìn)水管直徑、出水管直徑、流量均不變的前提下，提出以下方案：由于出水管中水流方向與水平方向呈斜交狀態(tài)，故可嘗試逆時(shí)針旋轉(zhuǎn)出水管，即改變圖1中的夾角α值，分別模擬α=75°、60°、45°、30°時(shí)的流場(chǎng)。

3.3 不同模型計(jì)算結(jié)果

表3為不同出水管角度的計(jì)算結(jié)果。由表3可以看出，出水管與罐體夾角α對(duì)壓降的影響非常小。當(dāng)α=30°時(shí)，壓降最小為0.997 m，當(dāng)α=60°時(shí)，壓降最大為1.051 m，兩者相差僅為0.054 m。出水管角度對(duì)壓降影響極小是因?yàn)闉V網(wǎng)面積、孔徑及濾網(wǎng)壓降特性均未變化，而壓降產(chǎn)生的主要部位是濾網(wǎng)，故夾角α的改變對(duì)壓強(qiáng)損失的影響可忽略不計(jì)。

表3 數(shù)值模擬結(jié)果Tab.3 Numerical simulation results

圖4為各個(gè)方案的速度云圖和矢量圖。從流場(chǎng)分布可以看出：①濾網(wǎng)內(nèi)有一塊低速回流區(qū)域(濾網(wǎng)尾部，靠近下側(cè)濾網(wǎng))，從濾網(wǎng)尾部一直延伸到出水管右側(cè)邊緣，在夾角α由75°變到30°的過(guò)程中，出水管與罐體相交的右側(cè)交點(diǎn)逐漸向右邊移動(dòng)(左側(cè)交點(diǎn)則向左邊移動(dòng))，導(dǎo)致回流區(qū)域也跟著向右側(cè)移動(dòng)，回流區(qū)域面積變小。在回流區(qū)域前部，水流流向?yàn)檎?，且流速分布均勻約為1.8 m/s，可見(jiàn)此處是濾網(wǎng)有效利用率最高的部位。在回流區(qū)域，流速較小，且流向各異，因此回流區(qū)域?yàn)V網(wǎng)的過(guò)濾效率較低。當(dāng)夾角α=30°時(shí)，回流區(qū)域最小，濾網(wǎng)有效利用率最高。除此之外，進(jìn)水管、排污管、濾網(wǎng)內(nèi)部速度分布特征基本與夾角α=90°時(shí)一致。②出水管內(nèi)流速分布有較大差異。隨著夾角α逐漸變小，出水管內(nèi)高流速區(qū)域面積先增大后減小，且位置也發(fā)生了較大變化。當(dāng)夾角α=90°時(shí)[圖3(a)]，高流速區(qū)域在出水管上部中間略微偏右的位置；當(dāng)夾角α=75°時(shí)[圖4(a)]，在出水管上部中間的位置；當(dāng)夾角α=60°時(shí)，在出水管中下部左側(cè)的位置；當(dāng)夾角α=30°時(shí)，在出水管左下部的位置，且形狀也不再是Y字型。③出水管內(nèi)漩渦區(qū)也隨著高流速區(qū)的改變而發(fā)生變化。當(dāng)夾角α=90°時(shí)[圖3(a)]，漩渦區(qū)在出水管左右兩側(cè)，左側(cè)大于右側(cè)。隨著夾角α變小，出水管內(nèi)左側(cè)漩渦漸漸消失，而右側(cè)漩渦卻在變大。由于夾角α變小，出水管左側(cè)水流能更平順地進(jìn)入出水管，因此左側(cè)漩渦消失。出水管右側(cè)水流為從過(guò)濾網(wǎng)尾部回流的水流，水流方向與出水管出流方向存在一個(gè)夾角，α越小該夾角越大，回流水需轉(zhuǎn)過(guò)的角度也越大，因此在出水管右側(cè)造成的漩渦也越大。值得注意的是，出水管左側(cè)水流流速大、流量大，右側(cè)水流流速小流量小，左側(cè)水流占據(jù)主流。因此將出水管逆時(shí)針旋轉(zhuǎn)，順應(yīng)主流是合理的。但是出水管右下部存在一個(gè)較大漩渦區(qū)同樣不可忽略，對(duì)于消除這個(gè)漩渦區(qū)可以嘗試將出水管向右平移，以進(jìn)一步減小出水管右側(cè)水流流量，達(dá)到減小漩渦的目的。

圖4 斷面速度場(chǎng)Fig.4 Velocity field of the section

圖5為各個(gè)方案壓強(qiáng)云圖。與α=90°時(shí)的壓強(qiáng)分布規(guī)律相比較，當(dāng)夾角α改變時(shí)，進(jìn)水管、排污管、罐體內(nèi)壓強(qiáng)分布仍然比較均勻，最大壓差約為2 kPa。主要的壓降還是在濾網(wǎng)上。但是隨著夾角α的改變，出水管內(nèi)壓強(qiáng)分布規(guī)律發(fā)生了明顯的變化：夾角α越小，出水管內(nèi)壓強(qiáng)分布越均勻。當(dāng)夾角α=75°時(shí)，出水管上方(即出水管與罐體相交處)，左右兩側(cè)即出現(xiàn)漩渦低壓區(qū)，且壓強(qiáng)分布不均勻。當(dāng)夾角α=45°時(shí)，出水管左側(cè)低壓區(qū)基本消失，雖然右側(cè)低壓區(qū)還存在，但是出水管內(nèi)總體壓強(qiáng)分布趨于均勻。當(dāng)夾角α=30°時(shí)，出水管上部壓強(qiáng)分布十分均勻，這對(duì)濾網(wǎng)均勻受壓，不至于局部破壞來(lái)說(shuō)十分有利。整個(gè)出水管內(nèi)，壓強(qiáng)沿水流方向逐漸減小，壓降有明顯的過(guò)渡過(guò)程。

圖5 斷面壓強(qiáng)場(chǎng)Fig.5 Pressure field of the section

4 結(jié) 語(yǔ)

(1)根據(jù)對(duì)過(guò)濾器壓強(qiáng)場(chǎng)的分析可知，壓降主要的部位為濾網(wǎng)。因此保持濾網(wǎng)有關(guān)特性不變，而改變出水管與罐體的夾角，對(duì)過(guò)濾器整體壓降影響不大，從焊接罐體與出水管的工藝來(lái)說(shuō)，選擇α=90°最為合適。

(2)改變出水管與罐體的夾角α，對(duì)進(jìn)水管、濾網(wǎng)內(nèi)部、排污管流場(chǎng)的影響比較小，而對(duì)出水管的流場(chǎng)有明顯的影響，主要表現(xiàn)在：隨著夾角α的減小，出水管上部流速、壓強(qiáng)分布均趨于均勻；出水管左側(cè)漩渦區(qū)域逐漸消失，右側(cè)漩渦區(qū)域有所增大。

(3)由于改變出水管角度對(duì)過(guò)濾器內(nèi)部流場(chǎng)影響不大，若要對(duì)過(guò)濾器進(jìn)行結(jié)構(gòu)優(yōu)化，改善流場(chǎng)分布，則需要從優(yōu)化濾網(wǎng)部分著手。

[1] 宗全利,劉 飛,劉煥芳,等.大田滴灌自清洗網(wǎng)式過(guò)濾器水頭損失試驗(yàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2012，16:86-92．

[2] 宗全利,劉煥芳,鄭鐵剛,等. 微灌用網(wǎng)式新型自清洗過(guò)濾器的設(shè)計(jì)與試驗(yàn)研究[J]. 灌溉排水學(xué)報(bào),2010,29(1):78-82.

[3] 駱秀萍,劉煥芳,宗全利,等. 自清洗網(wǎng)式過(guò)濾器排污流量的計(jì)算[J]. 排灌機(jī)械工程學(xué)報(bào),2012,30(5):588-591,608.

[4] 駱秀萍,劉煥芳,宗全利,等. 微灌自清洗網(wǎng)式過(guò)濾器水頭損失的試驗(yàn)研究[J]. 石河子大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,29(1):98-102.

[5] 劉 飛,劉煥芳,宗全利,等. 新型自清洗網(wǎng)式過(guò)濾器結(jié)構(gòu)優(yōu)化研究[J]. 中國(guó)農(nóng)村水利水電,2010,(10):18-21.

[6] 孫新忠. 離心篩網(wǎng)一體式微灌式過(guò)濾器的試驗(yàn)研究[J]. 排灌機(jī)械,2006,24(3):20-23.

[7] 李正平. 離心篩網(wǎng)一體式微灌過(guò)濾器水力性能試驗(yàn)[J]. 水利與建筑工程學(xué)報(bào),2012,10(6):156-159.

[8] 李振成,孫新忠. 離心篩網(wǎng)一體式微灌過(guò)濾系統(tǒng)砂粒運(yùn)移規(guī)律的試驗(yàn)研究[J]. 節(jié)水灌溉,2013,(2):14-16.

[9] Capra, Scicolone.Assessing dripper clogging and filtering performance using municipal wastewater[J].Irrigation and Drainage, 2005(54):871-879.

[10] Yurdem H, Demir V, Degirmencioglu A. Development of a mathematical model to predict clean water head losses in hydrocyclone filters in drip irrigation systems using dimensional analysis[J]. Biosystems Engineering, 2010,105(4):495-506.

[11] Duran-Ros M, Puig-Bargues J, Arbat G, et al.Performance and backwashing efficiency of disc and screen filters in microirrgation systems[J]. Biosystems Engineering, 2009,103(1):35-42.

[12] 巴 鵬,房元燦,譚效武. 基于CFD技術(shù)的管道過(guò)濾器內(nèi)部流場(chǎng)模擬及其結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J]. 潤(rùn)滑與密封,2011,36(4):98-101.

[13] 王棟蕾,宗全利,劉建軍. 微灌用自清洗網(wǎng)式過(guò)濾器自清洗結(jié)構(gòu)流場(chǎng)分析與優(yōu)化研究[J]. 節(jié)水灌溉,2011,(12):5-8,12.

[14] 王新坤,高世凱,夏立平,等. 微灌用網(wǎng)式過(guò)濾器數(shù)值模擬與結(jié)構(gòu)優(yōu)化[J]. 排灌機(jī)械工程學(xué)報(bào),2013,31(8):719-723.

[15] 阿力甫江·阿不里米提,虎膽·吐馬爾白,馬合木江·艾合買(mǎi)提,等. 直沖洗魚(yú)雷網(wǎng)式過(guò)濾器內(nèi)流場(chǎng)的數(shù)值模擬[J].節(jié)水灌溉,2014,(10):6-10.

[16] 阿力甫江·阿不里米提,陶洪飛,馬英杰,等. 不同出水口位置下全自動(dòng)魚(yú)雷網(wǎng)式過(guò)濾器內(nèi)部流場(chǎng)的數(shù)值模擬[J]. 灌溉排水學(xué)報(bào),2015,34(12):47-51.

[17] 李 浩,韓啟彪,黃修橋,等. 基于多孔介質(zhì)模型下微灌網(wǎng)式過(guò)濾器CFD湍流模型選擇及流場(chǎng)分析[J]. 灌溉排水學(xué)報(bào),2016,35(4):14-19.

[18] 王忠義,任翱宇,王紀(jì)達(dá),等. 管道過(guò)濾器流場(chǎng)數(shù)值模擬與實(shí)驗(yàn)[J]. 華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,43(1):75-79.

[19] 姚卓飛. 高溫陶瓷過(guò)濾器的設(shè)計(jì)與數(shù)值模擬研究[D]. 北京：華北電力大學(xué),2015.

[20] 劉 洋. 地下輸油管道泄漏污染物遷移過(guò)程仿真研究[D]. 黑龍江大慶：東北石油大學(xué),2015.