基于FLUENT的矩形液肥分配器出口流量均勻分配優(yōu)化設(shè)計(jì)

摘要：在液肥還田作業(yè)中，分配器各管路出口流量的均勻性是定量施肥的關(guān)鍵條件。雖然進(jìn)口的臥式分配器和立式分配器能夠?qū)崿F(xiàn)較高的輸肥精度，但成本高、維護(hù)難?，F(xiàn)有的矩形分配器成本低、結(jié)構(gòu)簡單，但各出肥管出口流量均勻性差，影響施肥精度與可控性。為解決這些問題，基于FLUENT流體仿真軟件開展矩形分配器的流體流動數(shù)值模擬，對矩形分配器結(jié)構(gòu)進(jìn)行改進(jìn)。結(jié)果表明，入肥口位置位于矩形分配器一側(cè)時，矩形分配器左側(cè)和右側(cè)流量差異達(dá)到10.27%，最大流量出口與最低流量出口質(zhì)量流率的差值為0.26 kg/s，各出口質(zhì)量流率的標(biāo)準(zhǔn)差為0.088 95，各出肥管口的流量差異性較大，這對液肥的均勻施用有較大的不利影響。通過改變?nèi)敕使苈凡季?，重新設(shè)計(jì)分配器結(jié)構(gòu)，將入肥口設(shè)置在分配器中部，在不增加制造難度和制造成本的基礎(chǔ)上能夠有效提升流量分配的均勻性。矩形分配器經(jīng)過改進(jìn)后，最大流量出口與最低流量出口質(zhì)量流率的差值由0.26 kg/s降為0.12 kg/s，左側(cè)和右側(cè)流量差異由10.27%降至0.18%，各出口質(zhì)量流率標(biāo)準(zhǔn)差由0.088 95降為0.038 7。改進(jìn)后的矩形分配器在流體分配的均勻性方面取得較大改善。

關(guān)鍵詞：液肥；計(jì)算流體力學(xué)；矩形分配器；均勻分配；數(shù)值模擬

中圖分類號：S224.21 文獻(xiàn)標(biāo)識碼：A 文章編號：2095?5553 （2024） 11?0234?07

Optimal design of uniform flow distribution of rectangular liquid fertilizer distributor

based on FLUENT

Liu Yi1， Song Pengxing1， Liu Rongguo1， Liu Jinghui2， Yang Yang2

（1. Shandong Agricultural Machinery Technology Extension Station， Jinan， 250000， China;

2. School of Engineering， Anhui Agricultural University， Hefei， 230036， China）

Abstract： In the operation of returning liquid fertilizer to the field， the uniformity of the flow rate of each pipeline outlet of the distributor is the key condition for quantitative application of fertilizer. However， although the imported horizontal distributor and vertical distributor can achieve high fertilizer delivery accuracy， they are costly and difficult to maintain. Although the existing rectangular distributor has low cost and simple structure， the uniformity of flow rate at the outlet of each fertilizer outlet pipe is poor， which affects the fertilizer application precision and controllability. In order to solve these problems， numerical simulation of fluid flow in the rectangular distributor was carried out based on FLUENT fluid simulation software， and the structure of the rectangular distributor was improved. The results show that when the position of the fertilizer inlet is located on one side of the rectangular distributor， the difference between the flow rate on the left side and the right side of the rectangular distributor reaches 10.27%， the difference between the mass flow rate of the maximum flow rate outlet and the minimum flow rate outlet is 0.26 kg/s， and the standard deviation of the mass flow rate of the outlets is 0.088 95， and the flow rate variability of the fertilizer outlet ports is large， which has a large negative impact on the uniform application of liquid fertilizer. By changing the layout of the fertilizer inlet pipeline， redesigning the distributor structure， and setting the fertilizer inlet port in the middle of the distributor， the uniformity of flow distribution can be effectively improved on the basis of not increasing the manufacturing difficulty and manufacturing cost. After the improvement of the rectangular distributor， the difference of mass flow rate between the outlet of the largest flow rate and the outlet of the lowest flow rate was reduced from 0.26 kg/s to 0.12 kg/s， the difference between the left and right flow rates was reduced from 10.27% to 0.18%， and the standard deviation of the mass flow rate of the outlets was reduced from 0.088 95 to 0.038 7. The improved rectangular distributor has achieved a great improvement of the homogeneity of the distribution of fluids.

Keywords： liquid fertilizer; computational fluid dynamics; rectangular distributor; uniform distribution; numerical simulation

0 引言

隨著我國畜禽糞污資源化利用工作的加速推進(jìn)，當(dāng)前亟待解決的關(guān)鍵問題在于如何有效暢通糞污還田渠道，以克服糞污資源化利用“最后一公里”的障礙[1]。我國每年畜禽糞污產(chǎn)生量高達(dá)3.8×109 t，但綜合利用率卻不足60%[2]。如果不對這大量的畜禽糞污進(jìn)行合理處理和利用，將會造成嚴(yán)重的環(huán)境污染和資源浪費(fèi)。利用厭氧發(fā)酵技術(shù)[3]，可以有效地將畜禽糞便等有機(jī)物進(jìn)行轉(zhuǎn)化。除了發(fā)酵過程中產(chǎn)生的沼氣外，發(fā)酵后的產(chǎn)物，即沼液和沼渣，同樣具有豐富的利用價值[4]。然而，作為液態(tài)肥料施用的沼液的含水量較高，導(dǎo)致其施用的難度和成本較高。目前，國內(nèi)對于液態(tài)糞肥的施用相對比較粗放，雖然已經(jīng)成功解決了機(jī)械化液肥施用從無到有的問題，但是如何做到從有到優(yōu)亟需進(jìn)一步開展研究[5]。

在液態(tài)糞肥資源利用方面，液肥的施用方式有直接噴灑、滴流、淺/深施和拖管等方式[6]。在這些施肥作業(yè)中，液肥分配器是關(guān)鍵部件，其作用是將總管中的液肥均勻、穩(wěn)定地分配到各個分管中。由于液肥的儲存和運(yùn)輸空間成本較高，為了方便沼液的施用，人們通常會將液肥發(fā)酵場地設(shè)置在附近有足夠需求的田地上[7]。液體糞肥還田臍帶式裝備系統(tǒng)針對液肥輸運(yùn)成本問題，直接利用拖管將液肥由附近發(fā)酵場地運(yùn)輸至分配器分配后進(jìn)行施肥。但該系統(tǒng)液肥分配的均勻性不高，因此需要對該系統(tǒng)中的關(guān)鍵部件分配器進(jìn)行優(yōu)化設(shè)計(jì)。

流體分配的均勻性一直是流體力學(xué)研究的重要方向之一，已經(jīng)有研究人員對分配器的改進(jìn)進(jìn)行了研究[8]。國外在液體肥施用機(jī)械方面的研究較早，技術(shù)也比較成熟，但由于價格昂貴，在國內(nèi)的應(yīng)用并不廣泛[9]。李文哲等[10]設(shè)計(jì)了一種液態(tài)沼肥暗灌施肥機(jī)械，該機(jī)械采用分配器對糞液進(jìn)行分配，轉(zhuǎn)子刀具在分配器內(nèi)部旋轉(zhuǎn)實(shí)現(xiàn)高壓脈沖防堵、攪拌防堵、管口邊緣清掃防堵等功能。然而，這種類型的分配器對于還田臍帶式裝備系統(tǒng)并不適用。

在研究中，通過使用傳統(tǒng)方法來優(yōu)化分配器的結(jié)構(gòu)面臨多個挑戰(zhàn)，如操作復(fù)雜性高、成本昂貴以及難以獲取關(guān)鍵的內(nèi)部流動數(shù)據(jù)（如壓力、速度和流量）[11]。因此，本文采用計(jì)算流體動力學(xué)（CFD）方法，并借助FLUENT軟件平臺進(jìn)行數(shù)值模擬，特別選取Realizable k-ε湍流模型，以優(yōu)化模擬分配器內(nèi)部的流體流動過程。

1 矩形分配器機(jī)械機(jī)構(gòu)

矩形分配器的功能是將田邊儲液罐或沼液池中的有機(jī)液肥通過入肥管輸送到矩形分配器，通過該裝置的內(nèi)部壓力差以確保液肥能夠均勻地通過出肥管路施撒到田間。

以矩形液肥分配器為研究對象開展矩形分配器結(jié)構(gòu)設(shè)計(jì)優(yōu)化研究。矩形分配器結(jié)構(gòu)如圖1所示，主要由入肥管路、分配器殼體和出肥管路組成。現(xiàn)有的入肥管直徑為80 mm，各個出肥管直徑為30 mm。根據(jù)施肥系統(tǒng)對工作流量的需求，分配器入口的流體速度設(shè)為3 m/s，以實(shí)現(xiàn)合乎工作要求的液肥輸送和分布。

2 矩形分配器計(jì)算模型

為實(shí)現(xiàn)液態(tài)糞肥分配器內(nèi)流動行為的數(shù)值計(jì)算，需要建立精細(xì)的網(wǎng)格模型，細(xì)致設(shè)定符合實(shí)際施肥需求的邊界條件，并在關(guān)鍵位置設(shè)置監(jiān)測點(diǎn)定量評估液態(tài)糞肥分配特性。然后選取SIMPLE算法，SIMPLE算法是一個成熟且廣泛應(yīng)用的CFD求解器，適用于各種不可壓縮流動問題，包括層流和湍流流動，能夠提供良好的數(shù)值穩(wěn)定性。最后對仿真結(jié)果進(jìn)行分析。

2.1 網(wǎng)格劃分

為了進(jìn)行精確的數(shù)值模擬、分析液體糞肥的分配和流動效果，首先使用三維建模軟件SolidWorks創(chuàng)建矩形液肥分配器的三維模型，將建立好的三維模型被導(dǎo)入前處理軟件SpaceClaim以進(jìn)行幾何優(yōu)化和流體域提取。然后利用FLUENT軟件自帶的FLUENT Meshing網(wǎng)格劃分模塊對分配器模型進(jìn)行網(wǎng)格生成。矩形液肥分配器共計(jì)生成463 542個網(wǎng)格單元，網(wǎng)格細(xì)節(jié)如圖2所示。

高質(zhì)量的網(wǎng)格構(gòu)建是模擬的關(guān)鍵要素[12]，對于精確評估液肥分配器的分配均勻性以及準(zhǔn)確模擬其內(nèi)部流場狀況具有重要作用。

2.2 計(jì)算模型建立

經(jīng)過發(fā)酵處理的液態(tài)糞肥可直接用作農(nóng)田肥料，根據(jù)文獻(xiàn)[13]測定結(jié)果，其固形物含量（TS）為8%，密度為1 044.67 kg/m3。在流體動力學(xué)仿真中，流體流動通常被劃分為層流和湍流兩種狀態(tài)，根據(jù)雷諾數(shù)的大小確定其具體狀態(tài)，計(jì)算如式（1）所示。

[Re=ρvDμ] （1）

式中： Re——雷諾數(shù)；

ρ——流體密度，kg/m3；

v——流體流速，m/s；

μ——流體運(yùn)動黏度，Pa ? s；

D——入肥口管徑，m。

液肥的運(yùn)動黏度在0.02～0.1 Pa ? s之間[14]，本文設(shè)定液肥的運(yùn)動黏度[μ]為0.05 Pa ? s，入口流速[v]為3 m/s，入口直徑[D]為0.08 m，液肥密度[ρ]為1 044.67 kg/m3。根據(jù)式（1）可得雷諾數(shù)Re為5 014，雷諾數(shù)Re大于雷諾數(shù)臨界值2 300，因此，本文流體流動形態(tài)是湍流。為實(shí)現(xiàn)湍流渦流的精確模擬計(jì)算，選擇具有較高湍流渦流精度的k-ε模型[15]，k-ε模型的方程包括質(zhì)量守恒方程、動量方程以及k-ε方程。

質(zhì)量守恒方程

[?ρ?t+??xi（ρui）=0] （2）

式中： [xi]——流體中某一點(diǎn)在空間中的位置坐標(biāo)（[i=1，2，3]），m；

[ui]——流體在[xi]方向上的速度分量，m/s。

動量方程

[?（ρui）?t+?（ρuiuj）?xj=-?p?xi+??xj（μ?ui?xj）+?τij?xi] （3）

其中，[τij=-ρu'iu'j]。

式中： [p]——作用在流體微元體上的壓力，Pa；

[τij]——Reynolds應(yīng)力，Pa；

[u'i]——[i]方向的速度平均值，m/s；

[u'j]——[j]方向的速度平均值，m/s。

k-ε方程

[?（ρk）?t+?（ρuik）?xi=??xjαkμeff?k?xj+Gk-ρε] （4）

[?（ρε）?t+?（ρuiε）?xi=??xjαεμeff?ε?xj+C1εεkGk-C2ερε2k-Rε] （5）

其中，[μeff=μ+μt]。

式中： k——湍流動能，[m2/s2]；

ε——湍流耗散率，[m3/s3]；

[αk]——湍動能的有效普朗特數(shù)的倒數(shù)；

[αε]——耗散率的有效普朗特數(shù)的倒數(shù)；

[μeff]——有效黏性系數(shù)；

[μt]——湍動黏度，Pa ? s；

[Gk]——湍動能的生成項(xiàng)；

[Rε]——附加項(xiàng)；

[C1ε、C2ε]——常數(shù)，分別為1.42、1.68。

以上述方程作為理論基礎(chǔ)，對于開展液體糞肥在分配器內(nèi)的流動行為的數(shù)值模擬，精確揭示液體糞肥流體的流動行為至關(guān)重要。

2.3 邊界條件設(shè)定與監(jiān)測點(diǎn)布置

為準(zhǔn)確模擬液態(tài)糞肥在分配器系統(tǒng)中的流動行為，根據(jù)設(shè)備運(yùn)行需求，恰當(dāng)?shù)卦O(shè)定邊界條件是關(guān)鍵。入口邊界采用速度入口條件，流速設(shè)為3 m/s，符合液肥施肥系統(tǒng)的工作要求。對于出口邊界，采用壓力出口條件，假設(shè)外部壓力等于標(biāo)準(zhǔn)大氣壓。同時，分配器內(nèi)壁面與流體接觸處采用無滑移邊界條件進(jìn)行模擬。

為全面分析液肥在分配器出口的流量分布和均勻性，鑒于矩形液肥分配器出肥管的對稱性，選擇在分配器一側(cè)的出口管道設(shè)置10個監(jiān)測點(diǎn)（P1～P10）進(jìn)行流量監(jiān)控，如圖3所示。通過獲取這些監(jiān)測點(diǎn)處的質(zhì)量流率數(shù)據(jù)，可以用來評估液肥分配的均勻性。

綜上，通過合理設(shè)置邊界條件和監(jiān)測點(diǎn)，能夠詳細(xì)了解液態(tài)糞肥在分配器中的流動特性和各出口的流量分布，這有助于評估和優(yōu)化液肥分配的效果。此外，這些數(shù)據(jù)也是進(jìn)一步優(yōu)化分配器結(jié)構(gòu)設(shè)計(jì)與參數(shù)的基礎(chǔ)，確保能更均勻有效地施用液肥。

2.4 計(jì)算方法

在對矩形分配器內(nèi)部流場進(jìn)行流體動力學(xué)仿真時，首先將在FLUENT Meshing軟件中生成的網(wǎng)格模型導(dǎo)入到流體仿真求解器中。在選擇求解器時，采用基于壓力的求解方法。為確保仿真結(jié)果的精確性，選擇Realizable k-ε模型作為描述流體黏性的模型，并輔以可擴(kuò)展的壁面函數(shù)，此組合有助于更準(zhǔn)確地模擬流場內(nèi)的湍流效應(yīng)。

在求解策略上，采用SIMPLE算法[16]來高效解決Navier-Stokes方程[17]。SIMPLE算法的一大優(yōu)點(diǎn)就是迭代過程的簡潔高效，通過反復(fù)修正壓力場，穩(wěn)步收斂到滿足連續(xù)性條件的速度和壓力場，巧妙地解決了速度場和壓力場的耦合問題。相比其他復(fù)雜的求解方法，SIMPLE算法不需要求解耦合方程系統(tǒng)，從而顯著降低了計(jì)算的復(fù)雜性和資源消耗，加之其迭代策略的穩(wěn)定性和可靠性，使得SIMPLE算法即便在網(wǎng)格不均勻或幾何結(jié)構(gòu)復(fù)雜的流場問題中也表現(xiàn)出色。仿真過程中，采用混合初始化方法設(shè)置初始條件，并對求解器進(jìn)行迭代計(jì)算，直至收斂以得到流場的數(shù)值模擬結(jié)果。這樣的流場數(shù)值模擬不僅有助于深入理解液體糞肥在分配器中的流動特性，而且為改進(jìn)設(shè)計(jì)提供依據(jù)，實(shí)現(xiàn)液肥施用的均勻性。

2.5 現(xiàn)有矩形分配器數(shù)值計(jì)算結(jié)果

現(xiàn)有矩形分配器數(shù)值計(jì)算如圖4所示。

液肥進(jìn)入矩形分配器時流速和壓力非常不均勻，特別在入肥管口分配器內(nèi)壁面位置。液肥流體動能因與管壁碰撞而產(chǎn)生較大損失，進(jìn)而導(dǎo)致跡線圖中出現(xiàn)大量渦流，造成分配器內(nèi)部壓力差異較大。由圖4（d）可知，在矩形分配器的左側(cè)，流體跡線十分雜亂，阻礙流體通過左側(cè)管口的流動。此外，由于入口處壓力存在，流體向右側(cè)流動的同時也會從兩側(cè)的出肥口流出一部分，然后到達(dá)右側(cè)壁面。由于壁面的阻礙，流體會發(fā)生回流，主要是通過右側(cè)的出肥管口流出。

矩形分配器各出口監(jiān)測點(diǎn)P1～P10質(zhì)量流率數(shù)值計(jì)算結(jié)果如圖5所示，分配器的左側(cè)質(zhì)量流率較低，而右側(cè)質(zhì)量流率較高。經(jīng)過統(tǒng)計(jì)P1～P5左側(cè)5個出口總的質(zhì)量流率為3.497 kg/s，P6～P10右側(cè)5個出口總的質(zhì)量流率為4.301 kg/s，分配器左側(cè)和右側(cè)流量差異達(dá)到10.27%，靠近入口一側(cè)的出口質(zhì)量流率較低；左側(cè)P2出口質(zhì)量流率最低，為0.635 kg/s，右側(cè)P9出口質(zhì)量流率最高，為0.895 kg/s，相差0.26 kg/s。

矩形分配器各出口監(jiān)測點(diǎn)P1～P10速度云圖如圖6所示，各出口截面平均流速如圖7所示，分配器的左側(cè)流速較低，右側(cè)流速較高。左側(cè)P2出口流速最低，為0.937 m/s，右側(cè)P9出口流速最高，為1.27 m/s，相差0.333 m/s。如圖6所示，左側(cè)出口的流速表現(xiàn)出較大的不均勻性。

由于入口流動的不均勻性和流體與壁面的相互作用，導(dǎo)致分配器內(nèi)流體流動的不均勻和渦流的產(chǎn)生。并且入肥口位置放置在矩形分配器左側(cè)這同樣對液體糞肥的分配效果產(chǎn)生了負(fù)面影響。因此，需要進(jìn)一步優(yōu)化分配器設(shè)計(jì)方案，改善流體的流動均勻性，減少渦流產(chǎn)生的負(fù)面影響，提高各出口流速的均勻性。

3 矩形分配器分配均勻性優(yōu)化

針對現(xiàn)有矩形分配器各出肥管流量不均勻問題，開展矩形分配器結(jié)構(gòu)優(yōu)化，通過改變?nèi)敕使苈凡季郑匦略O(shè)計(jì)分配器結(jié)構(gòu)，提高液體糞污施肥均勻性。

3.1 改進(jìn)后矩形分配器機(jī)械結(jié)構(gòu)

在對現(xiàn)有矩形分配器內(nèi)部流場進(jìn)行數(shù)值模擬后，發(fā)現(xiàn)其各出肥管管路出口質(zhì)量流率差異較大。由圖5可知，現(xiàn)有矩形分配器靠近入肥口一側(cè)P1～P5出口的整體流量低于遠(yuǎn)離入肥口一側(cè)的P6～P10出口的整體流量，流量差異達(dá)到10.27%。同時根據(jù)圖4現(xiàn)有矩形分配器流體仿真的后處理結(jié)果可以發(fā)現(xiàn)，當(dāng)入肥口位于矩形分配器一側(cè)時，分配器內(nèi)部流體的速度矢量、壓力矢量都集中在分配器一側(cè)，不利于液肥均勻分配，同時跡線圖也說明流體在分配器內(nèi)部流動不均勻。為解決該問題，將入肥口的位置設(shè)置在矩形分配器的中間位置，出肥口的位置不變，在不增加制造難度和制造成本的基礎(chǔ)上，優(yōu)化后的矩形分配器機(jī)械結(jié)構(gòu)如圖8所示。并對改進(jìn)后的矩形分配器內(nèi)部流場進(jìn)行數(shù)值模擬，驗(yàn)證改進(jìn)后的分配器的性能能否得到有效提升。

3.2 改進(jìn)后矩形分配器數(shù)值計(jì)算結(jié)果

在改進(jìn)后的矩形分配器一側(cè)管路出口仍然設(shè)置10個監(jiān)測點(diǎn)（P1～P10），用于評估液體糞肥在分配器中的流動情況和分布情況，如圖9所示。

改進(jìn)后的矩形分配器內(nèi)部液體糞肥流動數(shù)值計(jì)算如圖10所示。通過觀察速度矢量圖、壓力矢量圖和壓力云圖可以發(fā)現(xiàn)，液肥進(jìn)入矩形分配器時流速和壓力集中在分配器中部位置，沖擊壁面仍造成了較大的動能損失。液肥在分配器內(nèi)部左右兩側(cè)的行程是對稱的，在進(jìn)入分配器后一部分液肥從中間管口流出，其余部分向分配器兩側(cè)流動，到達(dá)左右兩側(cè)壁面后，液肥回流，由于分配器中部入口壓力的存在，液肥從左右兩側(cè)管道流出，在一定程度上降低了入口渦流帶來的負(fù)面影響。

改進(jìn)后的矩形分配器各出口監(jiān)測點(diǎn)P1～P10質(zhì)量流率數(shù)值計(jì)算結(jié)果如圖11所示，分配器各出口的質(zhì)量流率柱狀圖呈W形，左右兩側(cè)的中段位置的出肥口質(zhì)量流率較低，其余位置的出肥口質(zhì)量流率較高。經(jīng)過統(tǒng)計(jì)P1～P5左側(cè)5個出口總的質(zhì)量流率為3.878 kg/s，P6～P10右側(cè)5個出口總的質(zhì)量流率為3.892 kg/s，分配器左側(cè)和右側(cè)流量差異降至0.18%，這是由于將入肥口設(shè)置在分配器中部，流體流向左右兩側(cè)出口的工況相同。改進(jìn)后的矩形分配器P4出口質(zhì)量流率最低，為0.72 kg/s，右側(cè)P9出口質(zhì)量流率最高，為0.84 kg/s，相差0.12 kg/s。

改進(jìn)后矩形分配器各出口監(jiān)測點(diǎn)P1～P10速度云圖如圖12所示，各出口截面平均流速如圖13所示，分配器的左右兩側(cè)和中間部位出口流速較高，其余部位出口流速較低。P4出口流速最低，為1.05 m/s，中間P6出口流速最高，為1.22 m/s，相差0.17 m/s。

通過比較現(xiàn)有矩形分配器圖4的跡線圖和改進(jìn)后矩形分配器圖10的跡線圖，發(fā)現(xiàn)改變?nèi)敕士诘奈恢煤?，流體的軌跡更加均勻。現(xiàn)有矩形分配器左側(cè)和右側(cè)流量差異為10.27%，最大流量出口與最低流量出口質(zhì)量流率的差值為0.26 kg/s。改進(jìn)后的矩形分配器左側(cè)和右側(cè)流量差異降至0.18%，最大流量出口與最低流量出口質(zhì)量流率的差值降為0.12 kg/s。

為了定量評估各出口質(zhì)量流率的離散程度，通過標(biāo)準(zhǔn)差來衡量數(shù)據(jù)的波動性，以反映數(shù)據(jù)與均值之間的距離。標(biāo)準(zhǔn)差越大，說明數(shù)據(jù)之間的差異越大；標(biāo)準(zhǔn)差越小，說明數(shù)據(jù)之間的差異越小。標(biāo)準(zhǔn)差計(jì)算如式（6）所示。

[δ=1Ni=1N（mi-μ）2] （6）

式中： [mi]——各監(jiān)測點(diǎn)的質(zhì)量流率；

[μ]——平均值（算數(shù)平均值）；

[δ]——標(biāo)準(zhǔn)差。

利用式（6）計(jì)算傳統(tǒng)矩形分配器各出口質(zhì)量流率標(biāo)準(zhǔn)差[δ1=]0.088 95，改進(jìn)后矩形分配器各出口質(zhì)量流率標(biāo)準(zhǔn)差[δ2=]0.038 7，可以看出，[δ2]顯著小于[δ1]，說明改進(jìn)后的矩形分配器在流體分配的均勻性方面取得較大改善。通過定量分析方法，能夠更準(zhǔn)確地評估改進(jìn)后矩形分配器的性能，提高液肥施肥效率何質(zhì)量，為后續(xù)矩形分配器的優(yōu)化何改進(jìn)提供理論參考。

4 結(jié)論

1） 傳統(tǒng)矩形液肥分配器入肥口設(shè)置在矩形殼體左側(cè)，液肥進(jìn)入矩形分配器后，由于液肥流速較快，直對分配器內(nèi)壁面，形成較多渦流增大阻力從而導(dǎo)致左側(cè)出口質(zhì)量流率較低。由于入口壓力以及流體慣性的存在，右側(cè)出口質(zhì)量流率較高。左側(cè)5個出口總的質(zhì)量流率為3.497 kg/s，右側(cè)5個出口總的質(zhì)量流率為4.301 kg/s，分配器左側(cè)和右側(cè)流量差異達(dá)到10.27%，出口質(zhì)量流率最大值與最小值之間的差值達(dá)到0.26 kg/s。

2） 改進(jìn)后的矩形分配器中，入肥口直對中間部位的出肥管口，這導(dǎo)致中間部位的出肥口質(zhì)量流率較大。在流體向兩側(cè)流動后與壁面碰撞后發(fā)生回流，主要從兩側(cè)的出肥口流出，這些出肥口的質(zhì)量流率也相對較大，左右兩側(cè)的中段位置的出肥口質(zhì)量流率較低。左側(cè)5個出口總的質(zhì)量流率為3.878 kg/s，右側(cè)5個出口總的質(zhì)量流率為3.892 kg/s，分配器左側(cè)和右側(cè)流量差異僅為0.18%，出口質(zhì)量流率最大值與最小值之間的差值降至0.12 kg/s。通過改變?nèi)敕士谖恢?，矩形分配器左?cè)與右側(cè)的流量差異性顯著降低，最大流量出口與最低流量出口質(zhì)量流率的差值同樣顯著降低。由矩形分配器優(yōu)化前后各出口質(zhì)量流率的標(biāo)準(zhǔn)差[δ1=]0.088 95和[δ2=]0.038 7可知，將入肥口的位置設(shè)置在中間可以有效提升矩形分配器的工作性能。同時只是改變?nèi)敕士诘奈恢?，對分配器的整體結(jié)構(gòu)不進(jìn)行大的改變，可以發(fā)現(xiàn)，優(yōu)化后的矩形分配器在不增加制造難度和制造成本的基礎(chǔ)上能夠有效提升流量分配的均勻性。

參 考 文 獻(xiàn)

[ 1 ] 國家標(biāo)準(zhǔn)委 農(nóng)業(yè)農(nóng)村部 生態(tài)環(huán)境部關(guān)于推進(jìn)畜禽糞污資源化利用體系建設(shè)的指導(dǎo)意見[J]. 畜牧產(chǎn)業(yè)， 2023（12）： 5-13

[ 2 ] 吳浩瑋， 孫小淇， 梁博文， 等. 我國畜禽糞便污染現(xiàn)狀及處理與資源化利用分析[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報， 2020， 39（6）： 1168-1176.

Wu Haowei， Sun Xiaoqi， Liang Bowen， et al. Analysis of livestock and poultry manure pollution in China and its treatment and resource utilization [J]. Journal of Agro?Environment Science， 2020， 39（6）： 1168-1176.

[ 3 ] 李文哲， 徐名漢， 李晶宇. 畜禽養(yǎng)殖廢棄物資源化利用技術(shù)發(fā)展分析[J]. 農(nóng)業(yè)機(jī)械學(xué)報， 2013， 44（5）： 135-142.

Li Wenzhe， Xu Minghan， Li Jingyu. Prospect of resource utilization of animal faeces wastes [J]. Transactions of the Chinese Society for Agricultural Machinery， 2013， 44（5）： 135-142.

[ 4 ] 陳小華， 朱洪光. 農(nóng)作物秸稈產(chǎn)沼氣研究進(jìn)展與展望[J]. 農(nóng)業(yè)工程學(xué)報， 2007（3）： 279-283.

Chen Xiaohua， Zhu Hongguang. Research progress and prospect on producing biogas from crop straws [J]. Transactions of the Chinese Society of Agricultural Engineering， 2007（3）： 279-283.

（下轉(zhuǎn)第 252頁）

（上接第 239頁）

[ 5 ] 張克強(qiáng)， 杜連柱， 杜會英， 等. 國內(nèi)外畜禽養(yǎng)殖糞肥還田利用研究進(jìn)展[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報， 2021， 40（11）： 2472-2481， 2591.

Zhang Keqiang， Du Lianzhu， Du Huiying， et al. Application of livestock and poultry waste to agricultural land： A review [J]. Journal of Agro?Environment Science， 2021， 40（11）： 2472-2481， 2591.

[ 6 ] 何劍南. 液肥注射式工作部件的工作機(jī)理與試驗(yàn)研究[D]. 哈爾濱： 東北農(nóng)業(yè)大學(xué)， 2013.

[ 7 ] 王長軍. 沼氣池的選址與布局[J]. 吉林農(nóng)業(yè)， 2013（1）： 96.

[ 8 ] 唐永樂， 黃華， 胡偉東， 等. 數(shù)據(jù)中心用分流管流量均勻性及阻力特性的數(shù)值研究[J]. 制冷， 2022， 41（3）： 75-84.

[ 9 ] 白由路. 國內(nèi)外施肥機(jī)械的發(fā)展概況及需求分析[J]. 中國土壤與肥料， 2016（3）： 1-4.

[10] 李文哲， 袁虎， 劉宏新， 等. 沼液沼渣暗灌施肥機(jī)設(shè)計(jì)與試驗(yàn)[J]. 農(nóng)業(yè)機(jī)械學(xué)報， 2014， 45（11）： 75-80.

Li Wenzhe， Yuan Hu， Liu Hongxin， et al. Biogas slurry fertilizer applicator for dark irrigation [J]. Transactions of the Chinese Society for Agricultural Machinery， 2014， 45（11）： 75-80.

[11] 張文華， 李東來， 劉秀林， 等. 中心復(fù)合設(shè)計(jì)的水力旋流器結(jié)構(gòu)優(yōu)化與試驗(yàn)研究[J]. 機(jī)械科學(xué)與技術(shù)， 2023， 42（7）： 993-999.

[12] 李海峰， 吳冀川， 劉建波， 等. 有限元網(wǎng)格剖分與網(wǎng)格質(zhì)量判定指標(biāo)[J]. 中國機(jī)械工程， 2012， 23（3）： 368-377.

[13] Fu J， Chen Y， Xu B， et al. Investigation of distribution uniformity of distributor for biogas slurry application based on CFD analysis [J]. International Journal of Agricultural and Biological Engineering， 2023.

[14] 吳樹彪， 黎佳茜， 李偉， 等. 沼液回流對牛糞厭氧發(fā)酵產(chǎn)氣特性及其動力學(xué)的影響[J]. 農(nóng)業(yè)機(jī)械學(xué)報， 2015， 46（10）： 241-246.

Wu Shubiao， Li Jiaxi， Li Wei， et al. Effect of liquid digestate recirculation on biogas production and fermentation kinetics for anaerobic digestion of cattle manure [J]. Transactions of the Chinese Society for Agricultural Machinery， 2015， 46（10）： 241-246.

[15] 雷林， 王智祥， 孫鵬， 等. 計(jì)算流體力學(xué)k-ε二方程湍流模型應(yīng)用研究[J]. 船舶工程， 2010， 32（3）： 5-8.

[16] 任雄. 基于SIMPLE算法的三維湍流流場數(shù)值模擬方法的研究[D]. 西安： 西北工業(yè)大學(xué)， 2007.

[17] 肖志祥. 復(fù)雜流動Navier-Stokes方程數(shù)值模擬及湍流模型應(yīng)用研究[D]. 西安： 西北工業(yè)大學(xué)， 2003.