不同動態(tài)水壓模式下迷宮流道內顆粒物運動特性研究

鄭 超 吳普特, 張 林, 朱德蘭 趙 笑 安伯達

(1.西北農林科技大學水土保持研究所， 陜西楊凌 712100；2.西北農林科技大學中國旱區(qū)節(jié)水農業(yè)研究院， 陜西楊凌 712100)

不同動態(tài)水壓模式下迷宮流道內顆粒物運動特性研究

鄭 超1吳普特1,2張 林1,2朱德蘭2趙 笑2安伯達2

(1.西北農林科技大學水土保持研究所， 陜西楊凌 712100；2.西北農林科技大學中國旱區(qū)節(jié)水農業(yè)研究院， 陜西楊凌 712100)

采用PTV技術，觀測三角函數(shù)、三角、臺階、矩形波形動態(tài)水壓模式下迷宮流道內顆粒物的運動情況，分析不同動態(tài)水壓模式下迷宮流道內單個顆粒物運動軌跡、運動速度及流體流場，揭示動態(tài)水壓抗堵塞機理。結果表明：相對于三角、臺階及矩形波形動態(tài)水壓模式，三角函數(shù)波形動態(tài)水壓模式下顆粒物滯留滯止區(qū)內的時間最短，顆粒物的沉積概率最小，水流對顆粒物的輸移能力最大；三角函數(shù)波形動態(tài)水壓模式下流道內產生的水流波動效應不斷沖擊流道滯止區(qū)內存在的低速漩渦，加劇了滯止區(qū)水流紊動，進入滯止區(qū)的顆粒物隨水流不斷運動，并對滯留甚至沉積在滯止區(qū)內的顆粒物發(fā)生強烈沖擊，使得顆粒物離底懸浮，返回主流區(qū)，并迅速通過流道，增大了水流對顆粒物的輸移能力，增強了流道的抗堵塞性能。

動態(tài)水壓； 迷宮流道； 顆粒物； 運動速度； 運動軌跡

引言

圖2 迷宮流道灌水器及有機玻璃試驗樣件Fig.2 Schematics of labyrinth channel emitter and specimen

迷宮流道灌水器內流道結構復雜、尺寸微小(1 mm左右)，即使是水質良好且有較完備過濾措施的滴灌系統(tǒng)，仍有微小的固體顆粒物進入灌水器流道造成灌水器堵塞[1-4]。解決迷宮流道灌水器堵塞問題的關鍵是使灌水器流道具有較高的抗堵塞性能[5]。因此，眾多學者常通過優(yōu)選流道結構參數(shù)及優(yōu)化設計流道邊界，來保證流道內具有良好的流動狀態(tài)，提升流道內水流對顆粒物的輸移能力，進而達到提高流道抗堵塞能力的目的[6-12]。魏正英等[6]分析了迷宮流道內部流動場的情況及其堵塞機理，提出了迷宮流道主航道抗堵優(yōu)化設計方法。蘆剛等[7]通過兩相流數(shù)值模擬，分析了灌水器內水流的流場、固相物的運動軌跡及密度分布，通過短周期抗堵塞試驗，對3種渦體流道進行了短周期堵塞測試，并結合快速原形/快速制模(RP/RT)技術，提出了一種灌水器快速低成本開發(fā)的方法。喻黎明等[8]分析了梯形流道內含沙量分布及水沙流速分布，并以某一含沙量分布線作為流道邊界，對流道進行優(yōu)化，由此獲得抗堵塞性能較好的流道模型。但由于流道邊界條件及灌溉水質狀況復雜，使得灌水器物理堵塞問題仍未完全被解決。近年來，為提高滴灌系統(tǒng)及灌水器的抗堵塞性能，有學者嘗試提出動態(tài)水壓滴灌系統(tǒng)，這種方法有別于已有的脈沖式滴灌方式[13-14]，它通過動態(tài)的工作水壓加劇系統(tǒng)和灌水器內水流紊動，提高水流對顆粒物的輸移能力，從而達到改善系統(tǒng)和灌水器抗堵性能的目的[15-17]。王聰[16]分析了水壓模式、波動周期及振幅等波動參數(shù)對灌水器流量、均勻度和堵塞率的影響規(guī)律，得出隨著動態(tài)水壓波幅的增大，4種水壓模式產生的波動效應逐漸減弱，波幅為2～3 m時波動效應效果最好。劉潔等[17]通過正交壓力流量測試試驗，分析了基礎水壓、波幅及波動周期等波動參數(shù)對灌水器水力性能的影響，得出基礎水壓對灌水器流量影響最大。這些研究主要集中于灌水器宏觀水力特性方面，較少涉及灌水器流道內顆粒物運動情況及流場分析，動態(tài)水壓抗堵塞機理尚不清楚。

本文通過在水源中添加一定濃度的固體顆粒物，并采用高速攝像儀分別對4種典型的動態(tài)水壓模式(三角函數(shù)、三角、臺階及矩形波形)下灌水器流道內顆粒物運動情況進行連續(xù)拍攝，以獲得顆粒物的運動軌跡和速度等信息，進而對比分析不同動態(tài)水壓模式下迷宮流道內的顆粒物運動特性，揭示動態(tài)水壓抗堵塞機理，并優(yōu)選能最有效提高流道抗堵性能的動態(tài)水壓模式。

1 材料與方法

1.1 試驗材料

本研究選取常用的鋸齒型迷宮流道結構形式(以楊凌秦川節(jié)水灌溉設備工程有限公司生產的內鑲貼片式滴灌帶中灌水器迷宮流道為原型，其額定流量為2 L/h)進行水流運動特性研究，流道結構及參數(shù)詳見圖1。為了便于觀測流道內水流的運動情況，利用數(shù)控激光加工機床雕刻技術(雕刻精度為0.01 mm，表面粗糙度Ra≤0.25 μm)將平角齒形流道按1∶1的比例雕刻在透明的1 mm厚的有機玻璃板上，采用2塊3 mm厚的透光性好的有機玻璃板將流道密封，并采用8個螺栓加以固定，制成試驗樣件，如圖2所示。流道共13個單元，總長22.1 mm，流道寬度0.802 3 mm，齒高1 mm，流道深度1 mm，流道進出口平直流道加長，以便流體充分發(fā)展。在流道進口前4 mm處刻有測壓孔，并通過直徑2 mm的不銹鋼管及透明軟管連接壓力傳感器，以此達到實時監(jiān)測流道進口壓力變化的目的。

圖1 迷宮流道結構及參數(shù)(單位：mm)Fig.1 Form and structure parameters of labyrinth channel

1.2 試驗裝置與方法

試驗在西北農林科技大學中國旱區(qū)節(jié)水農業(yè)研究院灌溉水力學試驗廳進行，試驗裝置如圖3所示，主要由水箱、攪拌機(電動機轉速為900 r/min)、50WBZS15-22型不銹鋼自吸式耐腐蝕微型電泵、篩網過濾器(篩網孔徑為0.125 mm)、變頻柜、壓力表(量程40 m，精度0.5%)、壓力傳感器(西安新敏CYB型，量程40 m，精度0.1%)、數(shù)據(jù)采集器和PTV可視化試驗平臺(由不銹鋼支架、連續(xù)光源、高速攝像儀及VS-M0910型放大鏡組成)等組成。試驗采用的高速攝像儀每秒可采集0～200 000幀圖像，能清晰地捕捉到粒子在微小迷宮流道內的運動情況，利用后處理軟件Movias Pro Viewer 1.63可獲得粒子運動軌跡及瞬時速度等運動特征參數(shù)。試驗不同的動態(tài)水壓模式通過變頻柜控制加壓水泵電機的轉速來實現(xiàn)。采用水流跟隨性較好、密度接近于水的Nylon粒子(粒徑為4 μm，密度為1.04 g/cm3)作為固體懸浮顆粒物。

圖3 迷宮流道顆粒物運動特性試驗裝置Fig.3 Schematic of test system for particles movement characteristics of labyrinth channel1.不銹鋼離心泵 2.閥門 3.過濾器 4.水箱 5.攪拌機 6.壓力表 7.高速攝像儀 8.變頻柜 9.壓力采集器 10.塑料量杯 11.試驗樣件 12.連續(xù)光源 13.壓力傳感器

試驗時，在水箱中加入濃度為0.01 g/L的固體顆粒物，開啟攪拌器進行攪拌，使得顆粒物均勻分布于水箱的水體中；開啟變頻柜并設置變頻器運行參數(shù)，然后開啟水泵。調節(jié)完畢后，待壓力穩(wěn)定，通過連接在流道進口處的壓力傳感器測量流道進口壓力并由計算機記錄壓力傳感器瞬時采集(1次/s)的數(shù)據(jù)，用稱量法(測量2次，每次測量時間為20 min，取平均值)測量并記錄流道末端出流量。開啟連續(xù)光源及高速攝像儀，調整并設定高速攝像儀的焦距、數(shù)字式位移、快門速度及每秒采集幀數(shù)等參數(shù)，最后開始采集并保存高速攝像儀拍攝的圖像。采集圖像時，設定每秒采集6 000幀圖像，快門速度、數(shù)字移動及觸發(fā)模式等參數(shù)采用系統(tǒng)默認值。考慮到高速攝像儀觀測區(qū)域有限，無法觀測整個迷宮流道。因此，本研究選取位于迷宮流道中后部且流體發(fā)展充分的第7和第8流道單元作為主要的PTV觀測區(qū)域，其觀測區(qū)域長度為5.63 mm。

參照文獻[18]，對恒定水壓下流道壓力流量關系進行測定，得到流道的壓力流量計算公式為q=1.153H0.534(q為流道出流量，L/h；H為流道進口平均工作壓力水頭，m)。為了分析不同動態(tài)水壓模式下迷宮流道出流量的變化，對三角函數(shù)波形、三角波形、矩形波形及臺階波形4種動態(tài)水壓模式(不同動態(tài)水壓模式下實測工作壓力水頭變化見圖4)下波動周期T為30 s、振幅為3 m及流道進口基礎壓力水頭(流道進口平均工作壓力水頭)為4、6、8、10、12 m時迷宮流道出流量進行測定，結果表明相比恒定水壓，4種動態(tài)水壓模式下流道內流量未出現(xiàn)明顯下降，動態(tài)水壓對滴灌系統(tǒng)的供水能力影響較小。低壓灌溉是當代滴灌發(fā)展的趨勢，滴灌系統(tǒng)工作壓力的降低雖然降低了滴灌系統(tǒng)的實施及運行成本，但在低壓條件下滴灌灌水器流道更加容易發(fā)生堵塞。因此，在PTV觀測試驗時選定三角函數(shù)波形、三角波形、矩形波形及臺階波形4種動態(tài)水壓模式，選定較低的工作壓力水頭(4 m)作為流道進口基礎壓力水頭，波動周期T和振幅分別設為30 s和3 m，共4個試驗處理。

2 結果與分析

2.1 不同動態(tài)水壓模式下顆粒物運動速度

圖4 動態(tài)水壓模式示意圖Fig.4 Schematic diagrams of different dynamic water pressure modes

圖5 不同動態(tài)水壓模式下PTV觀測顆粒物運動軌跡及瞬時運動速度Fig.5 Trajectory and velocity of PTV observation about particles under different dynamic water pressure modes

圖5為不同動態(tài)水壓模式下工作壓力水頭為4 m左右及顆粒運動時間均為0.012 s左右時PTV觀測的第7和第8流道單元內單個顆粒物運動軌跡及瞬時運動速度。如圖5所示，4種動態(tài)水壓模式下顆粒物均進人了滯止區(qū)(流道可分為主流區(qū)和滯止區(qū)[19])，且顆粒物在滯止區(qū)的瞬時運動速度均較小，絕大部分在0.2 m/s(沙粒直徑為10 μm左右時的止動比速[20])以下，說明4種動態(tài)水壓模式下進入流道的顆粒物易隨水流進入滯止區(qū)，滯止區(qū)水流流速較低，顆粒物長時間滯留滯止區(qū)易發(fā)生沉積。三角函數(shù)波形動態(tài)水壓模式下顆粒物往返主流區(qū)和滯止區(qū)3次，平均每個流道單元內發(fā)生1次，且每次滯留滯止區(qū)的時間較短；三角及臺階波形動態(tài)水壓模式下發(fā)生2次往返運動，且每次滯留在滯止區(qū)的時間較長；而矩形波形動態(tài)水壓模式下僅發(fā)生1次往返運動，且每次滯留在滯止區(qū)的時間最長。由于三角函數(shù)波形動態(tài)水壓模式下流道內水流流速始終保持大幅度地連續(xù)上下波動，水流紊動加劇，并產生強烈的水流波動效應，對滯止區(qū)水流不斷發(fā)生沖擊，滯止區(qū)與主流區(qū)水流和能量交換更為頻繁，進入流道的顆粒物極易隨水流進入滯止區(qū)，但在水流波動效應的作用下能迅速地返回主流區(qū)，并最終通過整個流道；三角波形動態(tài)水壓模式下流道內工作水壓的瞬時波動幅度低于三角函數(shù)波形動態(tài)水壓模式下工作水壓的瞬時波動幅度，其水流紊動強度及產生的水流波動效應相對較小，進入滯止區(qū)的顆粒物隨水流不斷運動，在水流波動效應的作用下返回主流區(qū)，并最終通過流道；臺階及矩形波形動態(tài)水壓模式下流道內水流雖然流速波動也較大，但在波動期間較長時間流道內水流流速穩(wěn)定在一定值，使得水流波動效應減弱，且對滯止區(qū)水流沖擊的作用減弱，使得進入滯止區(qū)的顆粒物較長時間滯留在滯止區(qū)內，增大了顆粒物沉積的可能。

2.2 不同動態(tài)水壓模式下顆粒物運動軌跡線

由于示蹤粒子的布朗運動是由擴散引起的，單個顆粒物的運動情況存在較大的隨機性，難以代表流場內所有顆粒物的運動規(guī)律，但所有顆粒物的布朗運動屬于正態(tài)分布的隨機運動，因此可以采用概率統(tǒng)計的方法，分析流場內所有顆粒物的運動規(guī)律[21]。本研究以顆粒物在流道流場內運動圈數(shù)作為統(tǒng)計指標，對大量顆粒物的運動軌跡進行統(tǒng)計分析，進而分析整個流道流場內顆粒物的運動情況。采用后處理軟件Movias Pro Viewer 1.63對基礎壓力水頭為4 m時不同動態(tài)水壓模式下采集的顆粒物運動動態(tài)圖像進行處理，分別從圖像中隨機選取300個進入滯止區(qū)的顆粒物，描繪每個粒子的運動軌跡線，并對顆粒物運動圈數(shù)按1、2、3、4、5(含5以上)5個水平進行統(tǒng)計。圖6為不同動態(tài)水壓模式下顆粒物運動圈數(shù)累積概率曲線。如圖6所示，當顆粒物運動圈數(shù)為1時，三角函數(shù)波形動態(tài)水壓模式下顆粒物運動圈數(shù)的累積概率最大，累積概率為35%；三角波形及臺階波形的次之，累積概率均為31%左右；矩形波形的最小，累積概率為29%左右；當顆粒物運動圈數(shù)為3時，三角函數(shù)波形動態(tài)水壓模式下顆粒物運動圈數(shù)累積概率最大，累積概率為86%；三角波形及臺階波形的次之，累積概率分別為77.8%和77%；矩形波形的最小，累積概率為75.8%。當顆粒物運動圈數(shù)為3～5時，三角函數(shù)累積概率曲線斜率最為平緩，三角及臺階波形次之，矩形波形最大。由此可以看出，對于不同動態(tài)水壓模式，粒子運動圈數(shù)從少到多的次序為：三角函數(shù)、三角、臺階、矩形波形，說明三角函數(shù)波形動態(tài)水壓模式下水流對顆粒物的輸移能力可能最大。

圖6 不同動態(tài)水壓模式下流道內顆粒物運動 圈數(shù)累積概率曲線Fig.6 Probability cumulative curves for particles circle number under dynamic water pressure modes

圖7 不同運動圈數(shù)下顆粒物平均通過時間Fig.7 Histogram of particles average transit time in different particles circle numbers

圖7為不同動態(tài)水壓模式下顆粒物通過PTV觀測區(qū)域的平均時間。由圖7可知，隨著顆粒物運動圈數(shù)的增加，4種動態(tài)水壓模式下顆粒物平均通過時間均不斷增大；不同運動圈數(shù)時，三角函數(shù)波形動態(tài)水壓模式下顆粒物平均通過時間最短，三角及臺階波形動態(tài)水壓模式下的次之，矩形波形動態(tài)水壓模式下的最長。由此可以看出，對于不同動態(tài)水壓模式，顆粒物平均通過時間從小到大的次序為：三角函數(shù)、三角、臺階、矩形波形。進一步說明三角函數(shù)波形動態(tài)水壓模式下顆粒物在流道內停留時間最短，顆粒物在流道內的沉積概率最低，水流對顆粒物的輸移能力最大。

2.3 迷宮流道滯止區(qū)局部流場分析

從上述4種動態(tài)水壓模式下單個顆粒物運動情況的對比分析及多個顆粒物運動情況的統(tǒng)計分析可知，三角函數(shù)波形動態(tài)水壓模式能最有效提高流道抗堵塞性能。為了進一步揭示三角函數(shù)波形動態(tài)水壓模式下流道高抗堵塞性能的原因，對其流道內流場進行詳細分析。

對于鋸齒型迷宮流道，其齒尖后部流場較為復雜，且存在低速旋渦區(qū)，此處水流流速較小，甚至可能低于堵塞物質的起動速度，是流道中極易發(fā)生堵塞的區(qū)域[22-23]。因此，為了更加有效地分析迷宮流道流場細節(jié)，選取第7流道單元齒尖后部頂角局部區(qū)域作為流道流場分析的特征區(qū)域B(圖8)。采用最近鄰法對進入B區(qū)的顆粒物的運動情況進行匹配，獲得B區(qū)顆粒物速度矢量圖[24]，并分析進入B區(qū)近流道頂部0～0.1 mm處區(qū)域A(屬于流道近壁區(qū)域，是顆粒物易于發(fā)生沉積的區(qū)域[5])內的顆粒物平均運動速度變化。

圖8 PTV觀測區(qū)域Fig.8 Schematic of PTV observation region

圖10 迷宮流道頂角區(qū)域速度矢量圖Fig.10 Velocity vectors of regional flow fields of micro-channels vertex angle

圖9為三角函數(shù)波形動態(tài)水壓模式下基礎壓力水頭為4 m時一個完整周期T內近壁區(qū)內顆粒物平均運動速度變化曲線。隨著時間變化，三角函數(shù)波形動態(tài)水壓模式下流道內工作水壓和近壁區(qū)內顆粒物的平均運動速度均呈先不斷增大后不斷減小的趨勢，且其變化趨勢基本一致；三角函數(shù)波形動態(tài)水壓模式下近壁區(qū)內顆粒物的運動速度仍較低，但仍有一定程度的上下波動，從0.03 m/s左右不斷增至0.159 m/s左右，而后不斷下降至0.03 m/s左右。說明三角函數(shù)波形動態(tài)水壓模式下近壁區(qū)水流受流道內水流波動效應的影響，水流紊動加劇，進入近壁區(qū)的顆粒物易隨水流不斷運動，降低了顆粒物沉積的可能，甚至部分顆粒物在水流波動效應的作用下能夠迅速重返主流區(qū)，并最終通過流道。

圖9 近壁區(qū)內顆粒物平均運動速度變化曲線Fig.9 Variation curves of particles average velocity in near-wall region

圖10為三角函數(shù)波形動態(tài)水壓模式下基礎壓力水頭為4 m時一個完整周期T(進口水壓先從Hmin變化至Hmax，再從Hmax變化至Hmin)內迷宮流道齒尖后部頂角區(qū)域速度矢量圖。如圖10所示，流道內工作水壓波動初期，流道內工作壓力水頭最低，迷宮流道齒尖后部頂角區(qū)域存在較大的漩渦區(qū)，漩渦區(qū)水流流速及顆粒物運動速度均較低；隨著時間的變化，流道內工作壓力水頭不斷增大，迷宮流道齒尖后部頂角區(qū)域內低速漩渦區(qū)遭到沖擊，在接近主流區(qū)的邊緣位置顆粒物運動方向劇烈變化，并形成了為數(shù)較少的新的小漩渦；當工作壓力水頭增大到最大值時，迷宮流道齒尖后部頂角區(qū)域內低速漩渦區(qū)遭到強烈沖擊，整個漩渦區(qū)內顆粒物運動方向變化更為劇烈，原有的漩渦結構遭到嚴重破壞，并形成了為數(shù)較多的新的小漩渦；隨著時間的繼續(xù)變化，流道內工作壓力水頭不斷降低，迷宮流道齒尖后部頂角區(qū)域內遭到破壞的低速漩渦區(qū)逐漸恢復，僅在接近主流區(qū)的邊緣位置仍然存在為數(shù)較少的小漩渦；當流道內工作壓力水頭降低至最低值時，迷宮流道齒尖后部頂角區(qū)域內遭到破壞的低速漩渦區(qū)已經完全恢復。由此可以看出，隨著時間的變化，流道內工作壓力水頭不斷發(fā)生變化，流道內水流紊動劇烈，形成的水流波動效應使得流道滯止區(qū)低速漩渦不斷遭到沖擊并發(fā)生破壞，而后又逐漸恢復；在流道工作水壓較低時水流的波動效應使得進入滯止區(qū)的顆粒物隨水流不斷運動，不易在滯止區(qū)內產生沉積；在流道內工作壓力水頭較大時水流的波動效應使得進入滯止區(qū)的顆粒物能隨水流迅速通過流道，并對滯留甚至沉積在滯止區(qū)內的顆粒物發(fā)生強烈沖擊，使得顆粒物離底懸浮，返回主流區(qū)，并迅速通過流道，極大地增大了水流對顆粒物的輸移能力，從而增強流道的抗堵塞性能。

3 結論

(1)對于不同動態(tài)水壓模式，顆粒物運動圈數(shù)從少到多的次序為：三角函數(shù)、三角、臺階、矩形波形；顆粒物平均通過時間從小到大的次序為：三角函數(shù)、三角、臺階、矩形波形，說明三角函數(shù)波形動態(tài)水壓模式下顆粒物在流道內停留時間最短，顆粒物在流道內的沉積概率最低，水流對顆粒物的輸移能力最大。

(2)三角函數(shù)波形動態(tài)水壓模式下流道內水流流速始終保持大幅度地連續(xù)上下波動，水流紊動加劇，并產生強烈的水流波動效應，不斷沖擊流道滯止區(qū)的水流，使得流道滯止區(qū)內的低速漩渦不斷遭到沖擊并發(fā)生破壞，而后又逐漸恢復；使得進入流道滯止區(qū)的顆粒物在流道工作壓力水頭較低時隨水流不斷運動，不易在流道內產生沉積，在流道內工作壓力水頭較大時能隨水流迅速通過流道；使得滯留甚至沉積在滯止區(qū)內的顆粒物離底懸浮，返回主流區(qū)，并迅速通過流道，極大地增大了水流對顆粒物的輸移能力，從而增強流道的抗堵塞性能。

1 HUANG S, DONG W, WEI Q, et al. Advanced methods to develop drip emitters with new channel types [J]. Applied Engineering in Agriculture, 2006, 22(2): 243-249.

2 TAYLOR H D, BASTORS R K X, PERSON H W, et al. Drip irrigation with waste stabilization pond effluents: solving the problem of emitter fouling[J]. Water Science Technology, 1995, 31(12): 417-424.

3 王文娥, 王福軍, 嚴海軍. 迷宮滴頭CFD分析方法研究[J]. 農業(yè)機械學報, 2006, 37(10): 70-73. WANG Wen’e, WANG Fujun, YAN Haijun. Study on CFD method for flow simulation in labyrinth in the emitter[J]. Transactions of the Chinese Society for Agricultural Machinery, 2006, 37(10): 70-73. (in Chinese)

4 李云開, 楊培嶺, 任樹梅, 等. 分形流道設計及幾何參數(shù)對滴頭水力性能的影響[J]. 機械工程學報, 2007, 43(7): 109-114. LI Yunkai, YANG Peiling, REN Shumei, et al. Effects of fractal flow path designing and its parameters on emitter hydraulic performance[J]. Chinese Journal of Mechanical Engineering, 2007, 43(7): 109-114. (in Chinese)

5 馮吉, 孫昊蘇, 李云開. 滴灌灌水器內顆粒物運動特性的數(shù)字粒子圖像測速[J]. 農業(yè)工程學報, 2013, 29(13): 90-97. FENG Ji, SUN Haosu, LI Yunkai. Visualizing particles movement characteristics in drip irrigation emitters with digital particle image velocimetry[J]. Transactions of the CSAE, 2013, 29(13): 90-97. (in Chinese)

6 魏正英, 趙萬華, 唐一平, 等. 滴灌灌水器迷宮流道主航道抗堵設計方法研究[J]. 農業(yè)工程學報, 2005, 21(6): 1-7. WEI Zhengying, ZHAO Wanhua, TANG Yiping, et al. Anti-clogging design method for the labyrinth channels of drip irrigation emitters[J]. Transactions of the CSAE, 2005, 21(6): 1-7. (in Chinese)

7 蘆剛, 史玉升, 魏青松, 等. 基于兩相流模擬的高抗堵滴灌灌水器開發(fā)方法[J]. 華中科技大學學報：自然科學版, 2007, 35(7): 118-121. LU Gang, SHI Yusheng, WEI Qingsong, et al. A development method for high anti-clogged drip emitters based on two-phase flow simulations[J]. Journal of Huazhong University of Science and Technology: Natural Science Edition, 2007, 35(7): 118-121. (in Chinese)

8 喻黎明, 梅其勇. 迷宮流道灌水器抗堵塞設計與PIV試驗[J/OL]. 農業(yè)機械學報, 2014, 45(9): 155-160. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20140926&flag=1.DOI:10.6041/j.issn.1000-1298.2014.09.026. YU Liming, MEI Qiyong. Anti-clogging design and experimental investigation of PIV for labyrinth-channel emitters of drip irrigation emitters[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2014, 45(9): 155-160. (in Chinese)

9 喻黎明, 譚弘, 鄒小艷,等. 基于CFD-DEM耦合的迷宮流道水沙運動數(shù)值模擬[J/OL]. 農業(yè)機械學報, 2016, 47(8): 65-71. http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20160810&flag=1. DOI: 10.6041/j.issn.1000-1298.2016.08.010. YU Liming, TAN Hong, ZOU Xiaoyan, et al. Numerical simulation of water and sediment flow in labyrinth channel based on coupled CFD-DEM[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(8): 65-71. (in Chinese)

10 ZHANG J, ZHAO W, TANG Y, et al. Structural optimization of labyrinth-channel emitters based on hydraulic and anti-clogging performance[J]. Irrigation Science, 2011, 29(5): 351-357.

11 魏正英, 唐一平, 溫聚英, 等. 灌水器微細流道水沙兩相流分析和微PIV及抗堵實驗研究[J]. 農業(yè)工程學報, 2008, 24(6): 1-9． WEI Zhengying, TANG Yiping, WEN Juying, et al. Two-phase flow analysis and experimental investigation of micro-PIV and anti-clogging for micro-channels of emitter[J]. Transactions of the CSAE, 2008, 24(6): 1-9. (in Chinese)

12 牛文全, 吳普特, 喻黎明. 基于含沙量等值線的迷宮流道結構抗堵塞設計與模擬[J]. 農業(yè)工程學報, 2010, 26(5): 14-20. NIU Wenquan, WU Pute,YU Liming. Anti-clogging experimental investigation and optimized design of micro-channels of emitter based on isoline of sand content[J]. Transactions of the CSAE, 2010, 26(5): 14-20. (in Chinese)

13 ASSOULINE S, MOLLER M, COHEN S, et al. Soil-plant system response to pulsed drip irrigation and salinity[J]. Soil Science Society of America Journal, 2006, 70(5): 35-42．

14 ELMALOGLOU S, DIAMANTOPOULOS E. Wetting front advance patterns and water losses by deep percolation under the root zone as influenced by pulsed drip irrigation[J]. Agricultural Water Management, 2007, 90(1-2) : 160-163．

15 蘆剛. 工作水壓對滴灌灌水器水力性能影響規(guī)律的研究[D]. 武漢: 華中科技大學, 2010. LU Gang. Research on influence laws of the operating water pressure on hydraulic properties of drip emitters[D]. Wuhan: Huazhong University of Science and Technology, 2010. (in Chinese)

16 王聰. 波動水壓對灌水器水力性能影響的實驗研究[D]. 武漢: 華中科技大學, 2011. WANG Cong. Research on experiment of the effect of dynamic pressure on the hydraulic properties of drip Emitters[D]. Wuhan: Huazhong University of Science and Technology, 2011. (in Chinese)

17 劉潔, 王聰, 魏青松, 等. 波動水壓參數(shù)對灌水器水力性能影響試驗[J]. 河海大學學報: 自然科學版, 2014, 42(4): 361-366. LIU Jie, WANG Cong, WEI Qingsong, et al. Experimental study of effect of fluctuating water pressure factors on hydraulic properties of drip emitters[J]. Journal of Hohai University: Natural Sciences, 2014, 42(4): 361-366. (in Chinese)

18 GB/T 17187—2009 農業(yè)灌溉設備 滴頭和滴灌管 技術規(guī)范和試驗方法[S]. 2009.

19 鄭超, 吳普特, 張林, 等. 動態(tài)水壓下迷宮流道水流運動特性研究[J/OL]. 農業(yè)機械學報, 2015, 46(9): 167-172.http:∥www.j-csam.org/jcsam/ch/reader/view_abstract.aspx?file_no=20150924&flag=1.DOI:10.6041/j.issn.1000-1298.2015.09.024. ZHENG Chao, WU Pute, ZHANG Lin, et al. Flow characteristics in labyrinth channel under dynamic water pressure[J/OL]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(9): 167-172. (in Chinese)

20 錢寧, 萬兆惠. 泥沙運動力學[D]. 北京: 科學出版社，2003.

21 TAL S, ZUR B. Flow regime in helical long-path emitters[J]. Journal of the Irrigation & Drainage Division, 1980, 106(1): 27-35.

22 王文娥, 王福軍, 嚴海軍. 迷宮滴頭CFD分析方法研究[J]. 農業(yè)機械學報, 2006, 37(10): 70-73. WANG Wen’e, WANG Fujun, YAN Haijun. Study on CFD method for flow simulation in labyrinth emitter[J]. Transactions of the Chinese Society for Agricultural Machinery, 2006, 37(10): 70-73. (in Chinese)

23 金文, 張鴻雁. 灌水器內流道流場Micro-PIV試驗分析[J]. 農業(yè)工程學報, 2010, 26(2): 12-17. JIN Wen, ZHANG Hongyan. Micro-PIV analysis of flow fields in flow channel of emitter[J]. Transactions of the CSAE, 2010, 26(2): 12-17. (in Chinese)

24 禹明忠. PTV技術和顆粒三維運動規(guī)律的研究[D]. 北京: 清華大學, 2002. YU Mingzhong. Study on the PTV technique and 3D movement of particles[D]. Beijing: Tsinghua University, 2002. (in Chinese)

Particles Movement Characteristics in Labyrinth Channel under Different Dynamic Water Pressure Modes

ZHENG Chao1WU Pute1,2ZHANG Lin1,2ZHU Delan2ZHAO Xiao2AN Boda2

(1.InstituteofSoilandWaterConservation,NorthwestA&FUniversity,Yangling,Shaanxi712100,China2.InstituteofWater-savingAgricultureinAridAreasofChina,NorthwestA&FUniversity,Yangling,Shaanxi712100,China)

The movements of tracer particles in a labyrinth channel were observed by using the PIV technique. The pressure was in dynamic modes with trigonometric function, triangular, step and rectangular waveforms. The trajectories, velocity of single particle as well as the overall flow field in the channel were analyzed to reveal the anti-clogging mechanism of dynamic water pressure. The results showed that when compared with constant water pressure, flow rates in the channel under trigonometric function, triangular, step and rectangular waveforms dynamic pressure modes were not decreased significantly. Therefore, effect of dynamic water pressure on water delivery capacity of drip irrigation system was not obvious. Among the different dynamic water pressure modes mentioned above, the trigonometric function type demonstrated the shortest particles residence time in the flow stagnation region. The deposition probability of particles was the lowest under this mode, indicating a better particles transport capacity in the channel. Fluctuation effect generated by the water under pressure of trigonometric function form produced effect on the low-speed swirl in the stagnation region, and aggravated the turbulence of flow in this region. The turbulent flow was enhanced and the ability to carry the particles retained or even deposited in the stagnation zone was improved. As a result, the particles were suspended from the bottom and returned back to the main flow path. The transport capacity of the flow got strengthened, and the anti-clogging performance of the channel was improved.Therefore, the trigonometric function waveform dynamic water pressure model can most effectively improve the anti-clogging performance in the labyrinth channel.

dynamic water pressure; labyrinth channel; particles; velocity; trajectory

10.6041/j.issn.1000-1298.2017.03.037

2016-11-08

2016-12-28

中國科學院西部之光人才培養(yǎng)計劃項目、西北農林科技大學基本科研業(yè)務費專項資金項目(2014YB061)和“十二五”國家科技支撐計劃項目(2015BAD22B01-02)

鄭超(1989—)，男，博士生，主要從事節(jié)水灌溉理論研究，E-mail： zhengchao_xinong@163.com

吳普特(1963—)，男，研究員，博士生導師，主要從事水資源高效利用研究，E-mail： gizwpt@vip.sina.com

S275.6

A

1000-1298(2017)03-0294-08