邢少博，張金珠*，王振華，劉寧寧

(1. 石河子大學(xué)水利建筑工程學(xué)院，新疆 石河子 832000； 2. 現(xiàn)代節(jié)水灌溉兵團(tuán)重點(diǎn)實(shí)驗(yàn)室，新疆 石河子 832000)

滴灌系統(tǒng)中的滴灌灌水器通過(guò)不同結(jié)構(gòu)的流道將連續(xù)射流轉(zhuǎn)化為均勻穩(wěn)定的滴水狀出流[1]，發(fā)揮著穩(wěn)流降壓的作用，其流量受壓力變化影響的敏感程度稱為灌水器的水力性能[2].由于灌水器流道的尺寸較小且工作時(shí)間較長(zhǎng)，灌溉水質(zhì)復(fù)雜，易發(fā)生堵塞，因此灌水器的消能機(jī)理與抗堵塞性能成為當(dāng)前的研究熱點(diǎn).

目前，CFD(computational fluid dynamics, CFD)技術(shù)具有成本低、速度快且可視化流場(chǎng)易獲取的優(yōu)勢(shì)，故許多學(xué)者通過(guò)基于CFD的數(shù)值模擬與試驗(yàn)測(cè)試相結(jié)合的方法對(duì)灌水器的水力性能與消能機(jī)理進(jìn)行研究.唐學(xué)林等[3]通過(guò)兩相流數(shù)值模擬與PIV粒子測(cè)速試驗(yàn)相對(duì)比的方法，對(duì)迷宮流道內(nèi)的流場(chǎng)與顆粒運(yùn)動(dòng)情況進(jìn)行了計(jì)算，認(rèn)為數(shù)值模擬可以較為準(zhǔn)確地預(yù)測(cè)滴灌灌水器中的流場(chǎng)運(yùn)動(dòng).王新坤等[4]通過(guò)對(duì)高頻脈沖波下的迷宮流道灌水器進(jìn)行兩相流數(shù)值模擬，證明了高頻脈沖可以增強(qiáng)旋渦區(qū)的沖刷效應(yīng)以提高灌水器的抗堵塞能力.喻黎明等[5]研究了迷宮流道內(nèi)砂粒的運(yùn)動(dòng)規(guī)律及流體流場(chǎng)等，得出主流區(qū)流速較高、旋渦區(qū)流速較低、齒尖附近流速最大、靠近內(nèi)壁面和旋渦區(qū)中心流速最低的結(jié)論.門永南等[6]通過(guò)Fluent對(duì)矩形迷宮流道旋渦區(qū)進(jìn)行模擬，提出旋渦區(qū)可增強(qiáng)灌水器的消能效率與抗堵塞性能，但其零速區(qū)的存在也會(huì)導(dǎo)致灌水器堵塞.王新端[7]通過(guò)渾水試驗(yàn)對(duì)沉積泥沙的結(jié)構(gòu)特點(diǎn)以及粒徑分布進(jìn)行了分析，認(rèn)為灌水器結(jié)構(gòu)參數(shù)的變化對(duì)抗堵性能具有顯著影響.李云開等[8]提出了渦旋洗壁的概念，充分利用旋渦對(duì)壁面進(jìn)行清洗，使灌水器具備良好的自清洗和抗堵塞的能力.楚華麗等[9]通過(guò)多目標(biāo)優(yōu)化方法對(duì)卡門渦街形流道進(jìn)行優(yōu)化，認(rèn)為旋渦運(yùn)動(dòng)對(duì)水流耗能明顯，可促進(jìn)形成光滑無(wú)流動(dòng)死角的邊界，降低泥沙沉積.這些研究都肯定了CFD的數(shù)值模擬分析的可行性，并表明了旋渦區(qū)可以發(fā)揮抗堵塞及消能作用.

目前大量研究關(guān)注了優(yōu)化結(jié)構(gòu)參數(shù)以提高滴灌灌水器的抗堵塞性能與水力性能，但對(duì)旋渦區(qū)與抗堵塞及消能效果之間的具體關(guān)系缺少深入研究.依靠灌水器自身結(jié)構(gòu)產(chǎn)生的旋渦提高灌水器的抗堵塞與消能效果，從而增強(qiáng)流道“自清洗”能力[10]是灌水器發(fā)展的新思路.故文中以穿孔形流道為研究對(duì)象，通過(guò)不同壓力下的模擬從旋渦區(qū)流場(chǎng)、幾何特征、渦旋強(qiáng)度及壓力分布的角度，探究旋渦區(qū)的抗堵塞及消能作用機(jī)理.

1 材料與方法

1.1 穿孔形流道結(jié)構(gòu)參數(shù)

穿孔形流道的單元結(jié)構(gòu)如圖1所示.其中，流道深度為0.80 mm，寬度為2.76 mm,長(zhǎng)度為 60 mm，單元數(shù)為15；結(jié)構(gòu)1為流道腔，長(zhǎng)度為3.5 mm，結(jié)構(gòu)1中設(shè)有圓角，角3、角4、角5及角6的半徑分別為0.50，0.17，0.50，0.40 mm；結(jié)構(gòu)2為穿孔，長(zhǎng)度為0.50 mm，寬度為0.36 mm，中心軸距離為1.12 mm.

圖1 穿孔形流道單元結(jié)構(gòu)圖

1.2 流道模型

1.2.1 三維建模及網(wǎng)格無(wú)關(guān)性分析

采用SolidWorks及ICEM CFD對(duì)流道分別進(jìn)行三維建模與非結(jié)構(gòu)網(wǎng)格劃分.設(shè)置最大網(wǎng)格尺寸分別為0.20，0.19，0.18，0.17，0.15，0.14，0.13，0.12，0.11，0.10，0.09，0.08，0.07和0.05 mm，網(wǎng)格無(wú)關(guān)性分析結(jié)果顯示，網(wǎng)格尺寸減小，數(shù)量成倍增加，計(jì)算流量值先減小后趨于穩(wěn)定.隨網(wǎng)格數(shù)增加，計(jì)算時(shí)間成倍增加，綜合計(jì)算精度與效率的影響，將最大網(wǎng)格尺寸設(shè)置為0.07 mm，網(wǎng)格數(shù)為266.2萬(wàn).

1.2.2 流道數(shù)學(xué)模型

穿孔形流道內(nèi)部的水可近似認(rèn)為不可壓縮流動(dòng).由于標(biāo)準(zhǔn)k-ε模型計(jì)算較為穩(wěn)定，且穿孔形流道的流場(chǎng)中除近壁面外都為湍流，文中采用標(biāo)準(zhǔn)k-ε模型對(duì)流道內(nèi)部流場(chǎng)進(jìn)行計(jì)算[11].

1.2.3 邊界條件設(shè)置

進(jìn)水口設(shè)置為壓力入口，設(shè)定為0.050，0.075，0.100，0.125，0.150，0.175，0.200，0.225 MPa；出水口設(shè)置為壓力出口，壓力為0；所有壁面均采用無(wú)滑移邊界.顆粒在入口、出口處的邊界條件設(shè)為“逃逸”，邊壁的邊界條件設(shè)為“反彈”.顆粒在入口處釋放，顆粒密度為2 500 kg/m3，含沙率為1 g/L.針對(duì)目前常用的100目、120目和150目的過(guò)濾器，運(yùn)動(dòng)模擬采用的對(duì)應(yīng)的顆粒粒徑分別為0.150，0.125和0.100 mm.

1.2.4 計(jì)算方法設(shè)置

數(shù)值模擬計(jì)算采用非耦合隱式算法.其中，在Solve選項(xiàng)中將壓力項(xiàng)、動(dòng)量項(xiàng)、湍動(dòng)能項(xiàng)以及湍動(dòng)能耗散率項(xiàng)設(shè)置為精度更高的二階迎風(fēng)公式.為減小迭代計(jì)算誤差，將殘差收斂標(biāo)準(zhǔn)調(diào)整為1×10-4.為提高計(jì)算的穩(wěn)定性，壓力和速度耦合采用SIMPLE求解器，并將不同入口壓力下的模擬流量值與試驗(yàn)值進(jìn)行對(duì)比.

1.3 試驗(yàn)驗(yàn)證

采用高精度數(shù)控機(jī)床及透明無(wú)色有機(jī)玻璃對(duì)SolidWorks繪制的三維模型等比例加工流道樣機(jī)，并采用顯微鏡校核尺寸，再將雕刻好的流道與蓋板粘結(jié)組合.根據(jù)GB/T17187—2009《農(nóng)業(yè)灌溉設(shè)備-滴頭和滴灌管-技術(shù)規(guī)范和試驗(yàn)方法》布置試驗(yàn)系統(tǒng)，如圖2a所示，樣機(jī)加工5組，測(cè)試壓力分別為0.050，0.075，0.100，0.125，0.150，0.175，0.200，0.225 MPa，每片測(cè)試3 次，每次測(cè)試5 min，通過(guò)量筒容積法測(cè)試流量，并采用精度為0.001 g的電子天平校核，將實(shí)測(cè)值平均值與模擬值相互對(duì)照并繪制曲線圖，如圖2b所示，圖中Q為流量，p為壓力.

圖2 試驗(yàn)系統(tǒng)及結(jié)果對(duì)比曲線

模擬值與實(shí)測(cè)值的相對(duì)誤差為1.438%～2.299%，模擬可良好地反映流場(chǎng).穿孔形流道的流態(tài)指數(shù)為0.505，流量系數(shù)為0.302，相關(guān)系數(shù)大于0.999 8，水力性能優(yōu)良.

2 結(jié)果與分析

文中對(duì)0.05，0.10，0.15，0.20 MPa工作壓力下的穿孔形灌水器流道內(nèi)的流場(chǎng)分布、不同粒徑泥沙顆粒的運(yùn)動(dòng)情況、旋渦區(qū)的幾何特征、渦旋強(qiáng)度及壓力分布等進(jìn)行了分析.為降低入水口對(duì)流場(chǎng)的影響，根據(jù)流場(chǎng)在流道單元上的周期性，將第2單元作為主要研究對(duì)象.

2.1 流場(chǎng)分布

Z=0.40 mm截面為流道單元的中間截面，4種壓力下該截面上的速度云圖與流線圖相似.圖3a為0.10 MPa壓力下該截面的壓力云圖，圖3b為該截面的速度云圖與流線圖.流道的流場(chǎng)邊界光滑無(wú)流動(dòng)死角，可根據(jù)水流狀態(tài)將流場(chǎng)分為高速主流區(qū)與低速旋渦區(qū).主流區(qū)流速較高，旋渦區(qū)流速較低.文中針對(duì)旋渦區(qū)中的大旋渦進(jìn)行分析，以下大旋渦區(qū)稱為旋渦區(qū).

圖3 0.10 MPa壓力下Z=0.40 mm截面的壓力云圖、速度云圖與流線圖

表1為主流區(qū)和旋渦區(qū)的部分典型流速值.下文所稱主流區(qū)最大流速為流道腔中的最大流速.壓力增加3倍，主流區(qū)最大流速vmax1、最小流速vmin1與旋渦區(qū)邊界流速均增大了1倍，主流區(qū)的高速水流將顆粒沖入旋渦區(qū)，而旋渦區(qū)的離心力將顆粒甩出；旋渦區(qū)中心點(diǎn)流速量級(jí)較小且僅提高了0.86倍，導(dǎo)致顆粒在此處沉積.表1中vmax2為旋渦區(qū)邊界最大流速，vcent為旋渦區(qū)中心流速，△v1為主流區(qū)最大與最小流速差值，△v2為旋渦區(qū)最大與最小流速差值.

表1 不同壓力下主流區(qū)與旋渦區(qū)流速

圖4為穿孔形流道中粒徑dp為0.150, 0.125 和0.100 mm顆粒的運(yùn)動(dòng)模擬圖.在0.10 MPa工作壓力下，不同粒徑顆粒在流道中運(yùn)動(dòng)軌跡較為連貫，未在流道中發(fā)生堵塞與聚積；大部分顆粒在主流區(qū)中運(yùn)動(dòng)，小部分顆粒由主流區(qū)進(jìn)入旋渦區(qū)，隨后返回主流區(qū).

圖4 0.10 MPa壓力下不同粒徑顆粒運(yùn)動(dòng)模擬圖

根據(jù)表1及圖3b可知，主流區(qū)流速較高，不易形成顆粒沉積.部分顆粒進(jìn)入旋渦區(qū)后，又在邊界的離心力作用下返回主流區(qū)進(jìn)入下一流道單元，未在流道中聚積.不同壓力下，主流區(qū)帶動(dòng)旋渦區(qū)運(yùn)動(dòng)，旋渦區(qū)與主流區(qū)間發(fā)生摩擦.同時(shí)，旋渦區(qū)內(nèi)部流速不同，流層間又相互摩擦，水流能量下降，旋渦區(qū)發(fā)揮了消能作用.圖3a顯示，隨著主流區(qū)與旋渦區(qū)的相互作用，水流能量下降較為顯著，消能效果明顯.

2.2 不同工作壓力下的旋渦區(qū)幾何特征

圖5為在4種壓力下，Z=0.40 mm截面的旋渦區(qū)流線與速度矢量圖.由圖可知，壓力增加，旋渦區(qū)形態(tài)相對(duì)穩(wěn)定，其內(nèi)部各質(zhì)點(diǎn)速度矢量的運(yùn)動(dòng)趨勢(shì)為脫離旋渦區(qū)轉(zhuǎn)向主流區(qū).顆粒在進(jìn)入旋渦區(qū)時(shí)會(huì)受到水流的阻擋，少量進(jìn)入的顆粒在旋渦區(qū)受水流裹挾發(fā)生脫離，減少了進(jìn)入旋渦區(qū)的顆粒，旋渦區(qū)可發(fā)揮抗堵塞作用.

圖5 不同壓力下旋渦速度矢量與流線圖

將圖5中的最外側(cè)流線作為旋渦區(qū)邊界，采用樣條曲線對(duì)邊界進(jìn)行擬合繪制出旋渦區(qū)的形狀.以旋渦區(qū)最外側(cè)流線起點(diǎn)作為旋渦區(qū)中心點(diǎn)，樣條曲線包含區(qū)域?yàn)樾郎u區(qū)面積.

表2為不同壓力下旋渦區(qū)面積s與中心點(diǎn)位置.由表2可得，隨著壓力變化，旋渦區(qū)面積增加0.222 mm2，增幅較小，為2.35%，中心點(diǎn)向壁面及水流方向偏移，旋渦中流動(dòng)滯止區(qū)的位置不斷變化，表明旋渦區(qū)中心點(diǎn)不會(huì)產(chǎn)生大量顆粒聚積，不會(huì)引起堵塞.

表2 不同壓力下旋渦區(qū)面積與中心點(diǎn)位置

2.3 不同工作壓力下的旋渦區(qū)渦旋強(qiáng)度

旋渦區(qū)內(nèi)質(zhì)點(diǎn)運(yùn)動(dòng)形成的渦旋為物理渦，為表征渦運(yùn)動(dòng)的強(qiáng)弱，引入量綱為一的參數(shù)ω[6]，即渦旋強(qiáng)度為

ω=v/d，

(1)

式中：d為旋渦區(qū)等效圓半徑；v為旋渦區(qū)邊界平均切向速度.

根據(jù)面積與半徑的關(guān)系式公式，得到等效半徑r為

(2)

式中：s為根據(jù)旋渦區(qū)邊界線計(jì)算所得的旋渦區(qū)面積.

邊界平均切向速度為邊界上100個(gè)數(shù)據(jù)點(diǎn)的速度平均值.不同壓力下的邊界平均切向速度、等效半徑及渦旋強(qiáng)度見表3.由表3可得，壓力增大，等效半徑、切向速度及離心力增大，旋渦區(qū)內(nèi)的顆粒脫離趨勢(shì)增強(qiáng)，減緩了顆粒聚積引起的堵塞.渦旋強(qiáng)度增強(qiáng)，旋渦對(duì)邊壁的沖刷及與主流區(qū)的摩擦增強(qiáng)，顆粒在邊壁附近的聚積減緩，此外，旋渦受到邊壁的反作用而消耗能量，旋渦區(qū)發(fā)揮了抗堵塞作用及消能作用.

表3 不同壓力下旋渦區(qū)渦旋強(qiáng)度

2.4 不同工作壓力下的旋渦區(qū)內(nèi)部壓力

在4種工作壓力下，取流道第2單元中Z=0.40 mm截面與過(guò)旋渦區(qū)中心點(diǎn)橫截面的交線，分析旋渦區(qū)內(nèi)壓力p′沿該交線的變化，并繪制曲線圖如圖6.該交線與流道邊壁交點(diǎn)為A，與旋渦區(qū)下邊界的交點(diǎn)為B.

表4為不同壓力下旋渦區(qū)內(nèi)壓力差，表中，pb為邊界壓力，pc為中心點(diǎn)壓力，pw為邊壁壓力，△p1為邊界到中心點(diǎn)壓差，△p2為邊壁到中心點(diǎn)壓差，△p3為邊界與邊壁壓差.結(jié)合圖6及表4可得，壓力增大會(huì)顯著影響旋渦區(qū)壓力分布，但對(duì)壓力變化影響較小.不同壓力下，壓力最低及最高點(diǎn)分別出現(xiàn)在中心點(diǎn)O與邊界點(diǎn)B，BO，AO及AB間壓力差增加.結(jié)合圖5得，BO為AO單位距離壓力差的2倍以上；壓力增大3倍，各壓力差增加5倍左右，壓力差對(duì)壓力變化較為敏感.壓力差可減緩顆粒在邊壁的聚積，也增強(qiáng)了顆粒進(jìn)入旋渦區(qū)的趨勢(shì)并引起堵塞.但進(jìn)入旋渦區(qū)的顆粒又受到離心力及水流裹挾作用.根據(jù)圖4得，旋渦區(qū)的離心力及水流裹挾的作用強(qiáng)于壓力差的作用，顆粒不會(huì)受壓力差的影響而在旋渦區(qū)大量聚積，旋渦區(qū)可以發(fā)揮抗堵塞作用.

表4 不同壓力下旋渦區(qū)內(nèi)壓力差

3 結(jié) 論

以穿孔形灌水器為研究對(duì)象，結(jié)合流道內(nèi)不同粒徑泥沙顆粒的運(yùn)動(dòng)情況對(duì)4種壓力下流道的流場(chǎng)分布、旋渦區(qū)幾何特征、渦旋強(qiáng)度及壓力分布進(jìn)行了分析.

1) 旋渦區(qū)可發(fā)揮抗堵塞作用.少量進(jìn)入旋渦區(qū)的顆粒，在脫離旋渦區(qū)的水流裹挾及強(qiáng)于旋渦區(qū)內(nèi)部壓力的離心力作用下返回主流區(qū)，不會(huì)引起堵塞；旋渦的洗壁作用隨壓力增大而增強(qiáng)，邊壁顆粒沉積被抑制；旋渦區(qū)中心的流動(dòng)滯止區(qū)的位置隨壓力變化，不會(huì)產(chǎn)生顆粒聚積.

2) 旋渦區(qū)可發(fā)揮消能作用.旋渦區(qū)與主流區(qū)流速差較大，彼此摻混，旋渦區(qū)由主流區(qū)驅(qū)動(dòng)，消耗能量；旋渦區(qū)內(nèi)各流速層流速差異較大，摻混、摩擦消耗能量；旋渦區(qū)渦旋強(qiáng)度較高，與邊壁沖刷摩擦消耗能量.

3) 不同壓力下的旋渦區(qū)形態(tài)穩(wěn)定，可穩(wěn)定存在并發(fā)揮消能與抗堵塞作用.通過(guò)適當(dāng)利用灌水器自身結(jié)構(gòu)產(chǎn)生的旋渦區(qū)可提高灌水器的抗堵塞與消能性能.