陶 凱

（句容市水利農(nóng)機局，江蘇 鎮(zhèn)江 212400）

滴灌作為國際公認的一種高效節(jié)水灌溉技術(shù)，具有灌水均勻、高效節(jié)水節(jié)肥、增產(chǎn)增收、便于管理等優(yōu)點，被世界各國廣泛應(yīng)用，并快速發(fā)展，在我國農(nóng)業(yè)可持續(xù)發(fā)展和新農(nóng)村建設(shè)中發(fā)揮著重要作用。灌水器是滴灌技術(shù)的核心部件，其成本占滴灌系統(tǒng)造價的25%～35%[1]，主要作用是使壓力水流通過內(nèi)部的流道消能，使水流穩(wěn)定、均勻的滴入土壤，其水力性能的優(yōu)劣對滴灌系統(tǒng)的造價、能耗、灌水質(zhì)量及使用壽命等具有決定性的作用。但灌水器流道尺寸狹小，流態(tài)復(fù)雜，滴灌水中的雜質(zhì)易在流道內(nèi)沉積、凝結(jié)與黏附，使灌水器發(fā)生堵塞現(xiàn)象，降低滴灌灌水質(zhì)量和使用壽命[2-5]。滴灌水肥一體化、低壓滴灌、地下滴灌的發(fā)展[6-9]，再生水、黃河水、微咸水等在滴灌技術(shù)中的應(yīng)用，增加了灌水器堵塞的風(fēng)險因素[10-14]。堵塞已成為滴灌技術(shù)的行業(yè)難題，是影響滴灌技術(shù)發(fā)展的關(guān)鍵問題之一。本文構(gòu)建一種三角環(huán)流滴灌灌水器，并采用Fluent軟件進行計算流體動力學(xué)（CFD）模擬分析，開展灌水器內(nèi)部流動、固-液-氣多相耦合運動的分析，預(yù)測灌水器水力性能及抗堵塞性能。

1 三角環(huán)流滴灌灌水器結(jié)構(gòu)

在對灌水器進行CFD數(shù)值模擬之前，需要通過CAD對灌水器進行造型，如圖1所示，構(gòu)建一種三角環(huán)流滴灌灌水器。三角環(huán)流滴灌灌水器由進口、出口與流道主體3個部分組成，流道主體由30個流道單元依次連接而成，流道單元為2個三角形扁柱體嵌套形成的三角環(huán)流通道。進口部分由3個長度為1.5 mm的長方體組成，每個長方體間隔0.5 mm，出口部是1個1.5 mm×0.6 mm×0.7 mm的長方體。

圖1 三角環(huán)流滴灌灌水器結(jié)構(gòu)

2 CFD數(shù)值模擬

灌水器流道中的水流是水與各種泥沙、雜質(zhì)等顆粒的混合流動，屬于液-固兩相流動，采用計算流體動力學(xué)（CFD）模擬分析，能夠便捷的獲取灌水器流道內(nèi)固體顆粒流動與分布情況，預(yù)測灌水器水力性能與抗堵塞性能。

2.1 網(wǎng)格劃分

利用Gambit進行網(wǎng)格劃分，基于分段劃分的思想，將滴頭分為進口段、流道段和出口段分別進行結(jié)構(gòu)網(wǎng)格劃分，采用0.08 mm基本尺寸的六面體網(wǎng)格，考慮到流道單元拐角處為三角形塊，選用TGrid網(wǎng)格結(jié)構(gòu)對其進行網(wǎng)格劃分。為了更準(zhǔn)確地模擬近壁區(qū)的流動狀況，在近壁區(qū)加密網(wǎng)格，最終生成18萬個左右的單元網(wǎng)格，流道單元網(wǎng)格如圖2所示。

圖2 流道網(wǎng)格局部放大圖

2.2 數(shù)值計算

求解器選擇Fluent5/6，流體選用不可壓縮的水，溫度為室溫，運動粘度為v=1.0×10-6。入口邊界類型選擇壓力入口，出口邊界類型選擇壓力出口，湍流模型采用k-ε模型，用分離隱式穩(wěn)態(tài)求解器求解各控制方程，速度和壓力的耦合采用SIMPLE算法處理，各參數(shù)的離散均采用二階精度的迎風(fēng)格式。采用歐拉-拉格朗日法離散相模型來進行模擬分析，單向藕合方法來計算顆粒的運動。假設(shè)顆粒的形狀為球形，密度為2500 kg/m3，顆粒直徑d分別取0.01 mm、0.05 mm、0.1 mm、0.12 mm、0.15 mm等5個水平，體積濃度取1%，進口壓力設(shè)置為50 kPa，出口壓力為0，其他壁面條件為默認值。

3 結(jié)果與分析

3.1 顆粒粒徑對抗堵塞性能的影響

不同粒徑模擬的流量與濃度結(jié)果如表1所示，不同粒徑的顆粒相濃度分布如圖3所示。

由表1可知，粒徑大小不同時，灌水器流量變化很小，流道內(nèi)顆粒濃度差異較大。粒徑為0.01 mm時，最大濃度為1.446%；當(dāng)粒徑增加到0.1 mm時，最大濃度達到39%；而當(dāng)顆粒粒徑為0.05 mm時，最大濃度只有1.582%，說明三角環(huán)流滴灌灌水器具有較好的抗堵塞能力。隨著粒徑的增大，顆粒濃度呈線性增大。雖然粒徑大小對瞬時流量影響不大，但當(dāng)?shù)喂嚅g歇灌水的停止期間，毛管內(nèi)殘留的顆粒會結(jié)塊，濃度大結(jié)塊就大，再灌水時易引起灌水器堵塞。因此，對于圖1所示的三角環(huán)流滴灌灌水器，當(dāng)灌溉水源中有超過0.1 mm的顆粒粒徑時，須采取有效水處理措施給予去除。

表1 不同粒徑條件下的流量與濃度

圖3 不同粒徑顆粒濃度分布

由圖3可知，在流道單元拐角處、流道交匯處濃度均有升高，其他區(qū)域相對較低。流道單元拐角處的顆粒濃度和入口處濃度，與顆粒大小無明顯關(guān)系；流道拐角處背水區(qū)的顆粒濃度隨著粒徑大小的增加而降低。

由表1可以看出，粒徑與最大濃度之間存在明顯的線性關(guān)系，通過擬合，得出粒徑與流道內(nèi)顆粒濃度間的關(guān)系，見式（1）。

式中，VF為流道內(nèi)的顆粒濃度，%；d為灌溉水中顆粒物粒徑，mm。公式擬合的相關(guān)度R2=0.9983，可見顆粒相的最大濃度與粒徑有非常高的相關(guān)度。根據(jù)式（1）能夠通過水中顆粒粒徑，預(yù)測三角環(huán)流滴灌灌水器流道內(nèi)的顆粒濃度，可以作為灌水器設(shè)計或滴灌工程中水處理設(shè)施配置的依據(jù)。

3.2 入口濃度對抗堵塞性能的影響

表2為顆粒直徑0.1 mm，灌水器入口（灌溉水）顆粒濃度0.1%～3%，入口濃度與灌水器流量、流道內(nèi)濃度分布的影響。圖4為不同入口濃度條件下流道內(nèi)的濃度分布圖。

由表2可知，在入口濃度增大時，灌水器的流量也受到一定影響，但變化不明顯，入口濃度對流量沒有顯著作用。隨著入口濃度的增大，流道內(nèi)的顆粒濃度按照一定規(guī)律逐漸增大。通過擬合，得出入口濃度與流道內(nèi)顆粒濃度之間的關(guān)系，見式（2）。

表2 入口濃度與灌水器流量、流道內(nèi)濃度

圖4 不同入流顆粒濃度的濃度分布圖

式中，VF為流道內(nèi)的顆粒濃度，%；x為灌溉水中顆粒濃度，%。公式擬合的相關(guān)系數(shù)R2=0.9761，較好的揭示了兩者間的關(guān)系。隨著入口濃度的增大，流道中顆粒濃度增大。根據(jù)式（2）可以通過水中顆粒濃度，預(yù)測三角環(huán)流滴灌灌水器流道內(nèi)的顆粒濃度，能夠作為灌水器設(shè)計的依據(jù)。

整體上看，流道內(nèi)顆粒濃度的分布隨著灌水器入口濃度的升高而升高，流道內(nèi)各個位置處的顆粒濃度變化規(guī)律略有不同。隨著灌水器入口濃度的增加，在單元進水口處的顆粒濃度有所升高。流道拐角處的顆粒濃度與入口濃度沒有明顯關(guān)系，顆粒粒徑一定時，此處的顆粒濃度差異不大。灌溉水中顆粒濃度達到3%時，流道內(nèi)的顆粒濃度只有5.25%，流道內(nèi)的顆粒濃度增加梯度小于灌溉水中顆粒濃度增加梯度，說明三角環(huán)流滴灌灌水器流道內(nèi)泥沙顆粒的累積效應(yīng)不顯著，具有較好的抗堵塞能力。

4 結(jié)論

通過對三角環(huán)流灌水器兩相流的數(shù)值模擬分析得出，隨著顆粒直徑和體積濃度的增大，灌水器的流量沒有顯著的變化，而顆粒相在流道中的濃度分布呈一定規(guī)律變化。當(dāng)灌溉水中顆粒粒徑為0.05 mm時，最大濃度只有1.582%；流道內(nèi)的顆粒濃度增加梯度小于灌溉水中顆粒濃度增加梯度，說明流道內(nèi)泥沙顆粒的累積效應(yīng)不顯著。因此，三角環(huán)流滴灌灌水器具有較好的抗堵塞能力。擬合了粒徑與流道內(nèi)濃度的關(guān)系式及入口濃度與流道內(nèi)濃度的關(guān)系式，可以作為灌水器設(shè)計及滴灌工程水處理的依據(jù)。

流道中懸浮顆粒物的運動與流道中水流速度密切相關(guān)，水流速度越大其泥沙輸運能力越強。三角環(huán)流流道結(jié)構(gòu)復(fù)雜，導(dǎo)致流道中的流場分布雜亂，流道中流場局部存在漩渦，單元兩側(cè)拐角處流速梯度較大，懸浮顆粒濃度在流道中的分布發(fā)生了變化，有利于增加灌水器流道的抗堵塞能力。當(dāng)流速一定時，流道內(nèi)水流的泥沙輸運能力也是有限的，粒徑越大、入口濃度越高，灌水器流道內(nèi)局部的顆粒濃度也就越高。停止灌水后，部分顆粒將會留在流道內(nèi)，多次使用后顆粒將在這些部位沉淀結(jié)塊，容易造成堵塞。因此，滴灌灌溉水源應(yīng)根據(jù)使用的灌水器流道尺寸及結(jié)構(gòu)，使用合適的過濾裝置，盡量降低懸浮顆粒的大小及濃度，延長灌水器的使用壽命。