周良富，周立新，薛新宇，常 春，張學(xué)進(jìn)

(農(nóng)業(yè)部南京農(nóng)業(yè)機(jī)械化研究所，南京 210014)

0 引言

噴霧氣流的空間分布是影響霧滴輸運(yùn)的主要因素之一，氣流速度梯度是改變風(fēng)送霧滴動(dòng)力學(xué)特性的關(guān)鍵因子，因此開展噴霧氣流的速度空間分布及其在自由空間的衰減與擴(kuò)散規(guī)律研究，對(duì)深刻理解霧滴輸運(yùn)與沉積具有指導(dǎo)作用。

在風(fēng)送式果園噴霧技術(shù)研究中，國(guó)內(nèi)外學(xué)者一直關(guān)注著噴霧氣流場(chǎng)的研究。主要包括：①果園風(fēng)送噴霧機(jī)流道結(jié)構(gòu)對(duì)對(duì)空間氣流場(chǎng)分布影響研究。采用數(shù)值分析與試驗(yàn)相結(jié)合的方法[1]，研究噴霧機(jī)導(dǎo)流板角度[2]、風(fēng)機(jī)頻率[3-4]、導(dǎo)流片葉數(shù)[5]及風(fēng)機(jī)分布[6]等結(jié)構(gòu)參數(shù)對(duì)空間氣流場(chǎng)分布的影響。研究結(jié)果表明：風(fēng)機(jī)導(dǎo)流片有利于將旋轉(zhuǎn)氣流轉(zhuǎn)化為軸向氣流，但會(huì)造成壓力損失，普遍認(rèn)為導(dǎo)流葉片數(shù)以4～5 為宜；噴霧機(jī)外流場(chǎng)的氣流速度在垂直和水平兩個(gè)方向呈梯度減小趨勢(shì)；噴霧機(jī)出口氣流速度對(duì)氣流速度場(chǎng)的空間分布狀態(tài)無(wú)顯著影響,只是在數(shù)值上有差異。②噴霧氣流對(duì)果樹冠層的穿透、翻轉(zhuǎn)作業(yè)研究。Endalew 等建立了一個(gè)新的果園風(fēng)送噴霧 CFD 模型[7]，研究了噴霧機(jī)風(fēng)速與方向?qū)婌F氣流在冠層內(nèi)部分布的影響，證實(shí)氣流逆向線性衰減。同時(shí)，以該模型為基礎(chǔ)研究了3種(單個(gè)風(fēng)機(jī)、2個(gè)風(fēng)機(jī)和4個(gè)風(fēng)機(jī))不同噴霧機(jī)結(jié)構(gòu)[8]的噴霧氣流分布，結(jié)果顯示：該模型95%模型值的相對(duì)誤差小于20%。Gu等[9]測(cè)試了噴霧氣流在不同冠層結(jié)構(gòu)內(nèi)的速度分布規(guī)律，Hong等[10]在此數(shù)據(jù)基礎(chǔ)上仿真分析了冠層內(nèi)部氣流速度和壓力分布規(guī)律，模型的相對(duì)誤差小于30%。綜上文獻(xiàn)分析，采用試驗(yàn)與數(shù)值分析手段研究特定噴霧機(jī)結(jié)構(gòu)下噴霧氣流在自由空間及冠層內(nèi)分布規(guī)律已經(jīng)成為可能，但從自由射流理論角度出發(fā)，缺乏建立噴霧氣流在自由空間內(nèi)的衰減與擴(kuò)展規(guī)律模型的研究?；诖?，本文重點(diǎn)關(guān)注噴霧氣流衰減與擴(kuò)展規(guī)律理論模型建立及驗(yàn)證，為風(fēng)送噴霧機(jī)設(shè)計(jì)及噴霧技術(shù)參數(shù)選擇提供理論支持。

1 材料與方法

1.1 噴霧風(fēng)機(jī)及其氣流道

目前，果園風(fēng)送噴霧中的氣流道形式主要有環(huán)向出風(fēng)式、塔式、柔性管多頭式和獨(dú)立圓盤式等[11]，而環(huán)向出風(fēng)式氣流道是當(dāng)前普遍應(yīng)用的結(jié)構(gòu)，主要由旋轉(zhuǎn)葉輪、整流葉輪、分流板和流道板組成，如圖1所示。工作時(shí)，旋轉(zhuǎn)葉輪在驅(qū)動(dòng)力作用下旋轉(zhuǎn)，使氣流從軸向進(jìn)風(fēng)，流經(jīng)整流葉輪后氣流分布趨向均勻，在分流板和流道板作用下軸向氣流轉(zhuǎn)變?yōu)榄h(huán)向氣流，經(jīng)流道整流后的氣流環(huán)向分布趨于均勻，而后噴射向自由空間，并在自由空間衰減與擴(kuò)散。

1.2 測(cè)試項(xiàng)目及方法

1.2.1 氣流環(huán)向分布

出口風(fēng)速的測(cè)量數(shù)據(jù)是風(fēng)場(chǎng)數(shù)值分析邊界條件和初始條件的設(shè)置依據(jù)。噴霧機(jī)出口風(fēng)速受風(fēng)機(jī)轉(zhuǎn)向及流道影響，在不同出口位置氣流速度存在一定差異。試驗(yàn)將噴霧機(jī)出口均分8個(gè)測(cè)試點(diǎn)，分別編號(hào)從1到8。采用Kestrel 4500型風(fēng)速儀定點(diǎn)測(cè)定轉(zhuǎn)速為1 400、1 600、1 800、2 000r/min的各測(cè)試點(diǎn)的出口風(fēng)速，讀取5s測(cè)試時(shí)間內(nèi)的最大值，重復(fù)3次取平均值作為測(cè)試結(jié)果。

圖1 風(fēng)送噴霧機(jī)導(dǎo)流結(jié)構(gòu)

1.2.2 湍流特性

湍流強(qiáng)度是描述風(fēng)速隨時(shí)間和空間變化的程度，反映脈動(dòng)風(fēng)速的相對(duì)強(qiáng)度，是描述氣流場(chǎng)特性的重要特征量之一。本研究采用Kestrel 4500型風(fēng)速儀測(cè)試風(fēng)機(jī)轉(zhuǎn)速1 400r/min、噴霧機(jī)出口和距離出口1.5m時(shí)1min內(nèi)的氣流瞬態(tài)速度變化，即

(1)

1.2.3 氣流軸向衰減與擴(kuò)散

1)理論分析。設(shè)噴霧機(jī)氣流噴出射入靜止的自由空間內(nèi)，根據(jù)射流理論中射流主體段縱向時(shí)均速度分布的相似性(見圖2)，令

(2)

其中，u為某一斷面任一點(diǎn)縱向時(shí)均速度(m/s)；uc為斷面中心軸線上縱向最大時(shí)均速度(m/s)；y為縱向坐標(biāo)(m)；b為射流半擴(kuò)展厚度(m)。

圖2 自由射流理論

2)數(shù)值分析。數(shù)值分析選用Realizableκ-ε模型，其控制方程模型常數(shù)取值參考文獻(xiàn)[2]。采用Ansys15.0商業(yè)軟件中的Meshing模塊和Fluent模塊對(duì)計(jì)算域進(jìn)行網(wǎng)格劃分和流場(chǎng)仿真。為避免邊界對(duì)流場(chǎng)的影響，將模擬計(jì)算域擴(kuò)大到射流長(zhǎng)度為10m、擴(kuò)散末端為3m的錐扇形區(qū)域。計(jì)算域及網(wǎng)格劃分如圖3～圖5所示。其中，網(wǎng)格以六面體網(wǎng)格為主，網(wǎng)格節(jié)點(diǎn)217 700個(gè)，網(wǎng)格數(shù)量為204 012個(gè)。采用絕對(duì)壓力進(jìn)行計(jì)算，忽略地面對(duì)其風(fēng)場(chǎng)的影響，運(yùn)用二階隱式耦合求法，入口采用速度入口，3個(gè)出口均采用壓力出口邊界。采用壓力-速度耦合求解算法，壓力方程采用二階迎風(fēng)，其他方程采用QUICK法離散，設(shè)置計(jì)算殘差為1e-6，采用 Hybird Initialization 進(jìn)行初始化。

3)試驗(yàn)分析。在噴霧氣流射流軸向4.5m范圍內(nèi)，每0.5m設(shè)置測(cè)試線，在射流擴(kuò)散方向每0.25m設(shè)置測(cè)試線，在1m和2m高度上設(shè)置2個(gè)測(cè)試平面，如圖3所示。

圖3 風(fēng)速采樣點(diǎn)分布示意圖

在軸向主要采樣風(fēng)速儀讀取不同距離下的氣流速度值，而在射流擴(kuò)散方向主要結(jié)合卷尺和風(fēng)速儀測(cè)量出氣流速度為0.5m/s時(shí)的擴(kuò)展邊界值。

2 結(jié)果與分析

2.1 速度環(huán)向分布均勻性

不同轉(zhuǎn)速下8個(gè)出風(fēng)口測(cè)試點(diǎn)的氣流速度試驗(yàn)結(jié)果如表1所示。由表1可以看出：風(fēng)機(jī)轉(zhuǎn)速為1 400、1 600、1 800、2 000r/min時(shí)，噴霧機(jī)出口的平均風(fēng)速分別為13.56、18.88、20.46、22.69 m/s；出口風(fēng)速隨轉(zhuǎn)速增加而增大，同時(shí)相應(yīng)的標(biāo)準(zhǔn)差也隨轉(zhuǎn)速增加而增大，但不同轉(zhuǎn)速下各個(gè)測(cè)試點(diǎn)風(fēng)速變化的變異系數(shù)均小于14%。

表1 噴霧機(jī)出口處氣流速度 m/s

2.2 不同距離的湍流強(qiáng)度

噴霧機(jī)出口處和噴霧距離1.5m處的氣流速度瞬時(shí)值隨時(shí)間的變化如圖4所示。將試驗(yàn)數(shù)據(jù)按照式(2)計(jì)算湍流強(qiáng)度。試驗(yàn)結(jié)果顯示：噴霧機(jī)出口和距離出口1.5 m處的湍流強(qiáng)度分別為5.36%和11.5%。因此，在噴霧距離1.5m處的紊流強(qiáng)度高于噴霧機(jī)出口處，有利于氣流對(duì)果樹葉片的翻轉(zhuǎn)，便于霧滴在葉片背面沉積。

圖4 噴霧機(jī)氣流瞬態(tài)速度

2.3 氣流衰減與擴(kuò)散模型

2.3.1 基于量綱分析的理論模型

試驗(yàn)表明：噴霧機(jī)出口氣流雷諾數(shù)Re較大，所以在研究紊動(dòng)射流時(shí)忽略Re的影響。軸向最大速度僅與噴霧機(jī)出口氣流動(dòng)量、氣體密度ρ和軸向距離x有關(guān)[12]，即

uc=C1Mx1ρx2xx3

(3)

對(duì)式(3)進(jìn)行量綱分析如下：[uc]=LT-1；[M]=LML-3L2T-2=MT-2；[ρ]=ML-3；[x]=L。代入方程得：LT-1=(MT-2)x1(ML-3)x2Lx3。

對(duì)比量綱得：x1=1/2，x2=-1/2，x3=-1/2。代入方程(3)，得噴霧氣流軸向衰減規(guī)律為

(4)

同理，可令b=C1Mx4ρx5xx6，得x4=x5=0，x6=1。因此，噴霧氣流擴(kuò)展規(guī)律為

b=C1x

(5)

2.3.2 仿真結(jié)果分析

在初速度u為10、15、20m/s下的氣流速度分布數(shù)值分析結(jié)果中，提取距離噴霧出口0.5、1.5、2.5、3.5、4.5、5.5、7、8、9、10m處的軸向速度并繪制曲線，如圖5所示。由圖5可知：噴霧氣流在自由空間中快速衰減，噴霧氣流初速度為10、15、20m/s時(shí)，衰減到距離出口2.5m時(shí)的氣流速度僅為2.5、4.5、6m/s。

圖5 氣流軸向速度衰減規(guī)律

數(shù)值分析了初速度u0為10、15、20m/s下的氣流速度擴(kuò)散規(guī)律，試驗(yàn)結(jié)果如圖6所示。其中，圖6中所示的邊界為氣流速度大于0.5 m/s區(qū)域。

圖6 不同初速度下氣流擴(kuò)散規(guī)律

2.3.3 模型系數(shù)待定及驗(yàn)證

將已知的b0=0.13、u0=20m/s代入公式(4)，然后根據(jù)距離出口1.5、2.5、3.5、4.5、5.5、7、8、9、10m處軸向速度的測(cè)試值求出9個(gè)系數(shù)值后求平均值，作為模型系數(shù)C，得到該噴霧機(jī)的氣流軸向衰減規(guī)律模型式(6)。將不同距離的軸向氣流速度仿真值、模型值和試驗(yàn)值繪制曲線如圖7所示。試驗(yàn)結(jié)果顯示：試驗(yàn)值均略小于仿真值，將試驗(yàn)值與模型值作方差分析，顯示回歸方程顯著型檢驗(yàn)值F=1.634，遠(yuǎn)小于顯著值Fcrit=5.12。

(6)

按照?qǐng)D3的測(cè)試點(diǎn)布置方法，找出風(fēng)速為0.5m/s的邊界點(diǎn)后用卷尺測(cè)量處邊界點(diǎn)到射流軸線距離。圖8顯現(xiàn)的是距離出口4m內(nèi)的射流邊界值。以正半射流擴(kuò)展厚度線性擬合，結(jié)合公式(5)得到噴霧氣流擴(kuò)展規(guī)律為b=0.154 2x，模型決定系數(shù)為0.972 4。

圖7 氣流擴(kuò)展邊線

圖8 氣流衰減規(guī)律

3 結(jié)論

1)噴霧機(jī)環(huán)向出口風(fēng)速、風(fēng)速標(biāo)準(zhǔn)差隨轉(zhuǎn)速增加而增大，但其風(fēng)速變異系數(shù)均小于14%。

2)風(fēng)機(jī)轉(zhuǎn)速為1 400r/min時(shí)，噴霧機(jī)出口和距離出口1.5 m處的湍流強(qiáng)度分別為5.36%和11.5%。

3)噴霧氣流衰減模型uc=0.578u0x-0.5，模型顯著性檢驗(yàn)值F遠(yuǎn)小于顯著臨界值。

4)霧氣流擴(kuò)展模型為b=0.1542x，模型決定系數(shù)為0.972 4。