鄭子喬，羅 星

（中國農(nóng)業(yè)大學(xué) 煙臺研究院，山東 煙臺 264670）

食用菌種植過程培養(yǎng)料發(fā)酵質(zhì)量與人們健康之間的關(guān)系十分密切，食用菌種植過程培養(yǎng)料發(fā)酵方式，或多或少存在著高溫氣體[1]，對食用菌的種植造成一定的影響。培養(yǎng)料發(fā)酵的目的主要是增加原料的腐熟程度，殺滅料內(nèi)細(xì)菌，雜菌，和蟲害。通過發(fā)酵后的培養(yǎng)料發(fā)酵對食用菌種植的吸收效果更加明顯，所以對培養(yǎng)料進(jìn)行發(fā)酵處理，是十分有必要的。由于培養(yǎng)料發(fā)酵中會釋放高溫氣體，導(dǎo)致食用菌種植過程培養(yǎng)料發(fā)酵氧氣含量得不到保障。這些高溫氣體中還存在致癌物質(zhì)，因此食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)設(shè)計變得尤為重要[2]。隨著人們對食用菌種植過程中氧氣含量的關(guān)注度越來越高，該領(lǐng)域相關(guān)學(xué)者提出了很多有關(guān)食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)設(shè)計的可行方案。

通風(fēng)可有效降低食用菌種植過程培養(yǎng)料發(fā)酵過程高溫氣體的含量，對食用菌種植過程種植人員的工作和生活十分重要[3]。針對上述食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)相關(guān)設(shè)計方法存在的問題，提出利用數(shù)值模擬方式將軸流風(fēng)機(jī)應(yīng)用在食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)設(shè)計中。

1 基于軸流風(fēng)機(jī)的食用菌種植過程培養(yǎng)料通風(fēng)發(fā)酵模擬

1.1 軸流風(fēng)機(jī)控制數(shù)學(xué)模型和邊界條件

根據(jù)經(jīng)驗可知，在食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)設(shè)計中，食用菌種植過程培養(yǎng)料發(fā)酵固體壁面充分發(fā)展下的湍流流動能夠分成湍流核心區(qū)域與近壁區(qū)域兩部分，由此湍流流動控制方程組會隨流動所處位置不同而不同[4]。

經(jīng)研究發(fā)現(xiàn)，RNGκ-ε模型能夠很好地處理高應(yīng)變率和流線彎曲程度比較大的流動情況。因此，在湍流中心位置使用RNGκ-ε模型。由單一介質(zhì)、不可壓縮和風(fēng)機(jī)運(yùn)行中無熱量交換假設(shè)可知，控制方程組無需考慮組分方程與能量方程，僅包含連續(xù)性方程和動量方程、κ方程以及ε方程即可。

綜合以上軸流風(fēng)機(jī)控制數(shù)學(xué)模型分析，設(shè)置軸流風(fēng)機(jī)食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)數(shù)值模擬邊界條件和初始參數(shù)。

針對軸流風(fēng)機(jī)在食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)設(shè)計中的應(yīng)用，構(gòu)建由集流器進(jìn)口至擴(kuò)散器出口的全場流幾何模型，再選擇適合的湍流模型。以上述內(nèi)容為基礎(chǔ)，對軸流風(fēng)機(jī)全場流中三維湍流流動通過CFD軟件實(shí)現(xiàn)數(shù)值計算。

選擇一個食用菌種植大棚，食用菌種植過程培養(yǎng)料發(fā)酵的空氣流速比較低，可當(dāng)作不可壓縮流體；食用菌種植過程培養(yǎng)料發(fā)酵過程含有的高溫氣體依照恒定速率揮發(fā)；空氣的物理性質(zhì)是定值，并且不會產(chǎn)生化學(xué)反應(yīng)；根據(jù)該前提假設(shè)，可有效提升研究結(jié)果與實(shí)際情況的擬合度。

三維流場模型邊界條件為：流體是常溫下的空氣。把軸流風(fēng)機(jī)入口設(shè)定為質(zhì)量流量入口，在集流器和流線罩構(gòu)建的流體區(qū)域和第一級葉輪流道區(qū)域間、一級葉輪流道區(qū)域和二級葉輪流道間、二級葉輪流道和擴(kuò)散器構(gòu)成的流體區(qū)域間設(shè)定為混合面。軸流風(fēng)機(jī)前后兩級葉片和輪轂都為靜墻。

1.2 模型構(gòu)建

軸流風(fēng)機(jī)的內(nèi)部構(gòu)造與流道形狀十分復(fù)雜，以更加合理地反應(yīng)風(fēng)機(jī)內(nèi)部流場特征為目的，要對風(fēng)機(jī)內(nèi)部關(guān)鍵位置進(jìn)行網(wǎng)格加密[5]。結(jié)合食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)特性和通風(fēng)機(jī)內(nèi)部各個部件結(jié)構(gòu)特性，在構(gòu)建通風(fēng)機(jī)內(nèi)部全流場幾何模型過程中，采用分區(qū)分塊方式，將整個幾何模型劃分成集流器、前后級動葉輪、前后級導(dǎo)葉輪與擴(kuò)散器幾個區(qū)域。等幾個子區(qū)域流場模型生成之后，依據(jù)各子區(qū)域流場的幾何模型端面交界種類將各子區(qū)域沿著軸向進(jìn)行連接，以此生成全流場幾何模型。

全流場幾何模型選用專用模型制作軟件構(gòu)成。對構(gòu)建的全流場幾何模型進(jìn)行網(wǎng)格生成處理，進(jìn)而實(shí)現(xiàn)軸流風(fēng)機(jī)食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)設(shè)計的有限元模型構(gòu)建。

網(wǎng)格不僅是計算流體力學(xué)模型，即CFD的幾何表征形式，還是模擬和分析的載體。網(wǎng)格質(zhì)量對于流體力學(xué)模型的精準(zhǔn)計算十分重要。鑒于基于軸流風(fēng)機(jī)的食用菌種植過程培養(yǎng)料通風(fēng)發(fā)酵模擬中的幾何模型是十分復(fù)雜的3D幾何模型，因此，選用非結(jié)構(gòu)化網(wǎng)格，網(wǎng)格劃分中的網(wǎng)格單元選擇的是四面體單元。通過精細(xì)化的網(wǎng)格劃分，判斷軸流風(fēng)機(jī)疏散食用菌種植過程培養(yǎng)料發(fā)酵高溫氣體最佳參數(shù)，進(jìn)而最大程度降低食用菌種植過程培養(yǎng)料發(fā)酵高溫氣體。

網(wǎng)格是計算區(qū)域離散的基礎(chǔ)，也是食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)設(shè)計精細(xì)化的依據(jù)。合理的網(wǎng)格分布可以提升數(shù)值計算效率，節(jié)省計算時長，還能夠充分考慮流場中關(guān)鍵結(jié)構(gòu)對于流場參數(shù)大體分布的影響。

上述已經(jīng)將食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)整個流場幾何模型依據(jù)軸流風(fēng)機(jī)部件分為集流器等若干子區(qū)域。在網(wǎng)格生成過程中，也需要將軸流風(fēng)機(jī)全流場依照這些子區(qū)域分別生成，同時進(jìn)行不同程度的加密。

當(dāng)前針對復(fù)雜區(qū)域流場計算相關(guān)問題大多是經(jīng)專用網(wǎng)格生成軟件生成的。在此，也使用專用網(wǎng)格生成軟件生成軸流風(fēng)機(jī)有限元模型，如圖1所示。

圖1 軸流風(fēng)機(jī)網(wǎng)格劃分Fig.1 Axial fan meshing

2 實(shí)驗結(jié)果與分析

為驗證基于軸流風(fēng)機(jī)的食用菌種植過程培養(yǎng)料通風(fēng)發(fā)酵模擬有效性，在以下幾方面驗證所提方法。

1）軸流風(fēng)機(jī)數(shù)值模擬下的風(fēng)速分析

針對食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)具有的特殊性質(zhì)，軸流風(fēng)機(jī)的風(fēng)速不可太高，不然會影響食用菌種植過程種植人員舒適程度和正常活動，軸流風(fēng)機(jī)運(yùn)行工作頻率為50 Hz時的食用菌種植過程培養(yǎng)料發(fā)酵風(fēng)速矢量如圖2所示。

由圖2 a中可知，軸流風(fēng)機(jī)附近氣流速度比較大，高速氣流和扇葉以及風(fēng)機(jī)內(nèi)部壁面之間摩擦?xí)刹糠衷肼?。由圖2 b可知，當(dāng)軸流風(fēng)機(jī)運(yùn)行工作頻率為50 Hz時，食用菌種植過程培養(yǎng)料發(fā)酵氣流總體可以保持一個較為平穩(wěn)與低速的狀態(tài)，但風(fēng)機(jī)出口位置風(fēng)速相對比較大，這符合了風(fēng)機(jī)實(shí)際抽風(fēng)排風(fēng)特性。

圖2 風(fēng)機(jī)工作頻率為50 Hz時食用菌種植過程培養(yǎng)料發(fā)酵風(fēng)速矢量分布Fig.2 Vector distribution of fermenting wind speed of culture material in edible fungus culture process when the power frequency of the fan is 50 Hz

2）為了進(jìn)一步驗證所提方法可行性，分別在食用菌種植過程培養(yǎng)料發(fā)酵高溫氣體濃度去除率、模擬情況與實(shí)際值擬合度兩方面驗證基于軸流風(fēng)機(jī)的食用菌種植過程培養(yǎng)料通風(fēng)發(fā)酵模擬運(yùn)行效果。實(shí)驗結(jié)果如下。

圖3 不同方法食用菌種植過程培養(yǎng)料發(fā)酵高溫氣體濃度去除率對比Fig.3 Comparison of removal rates of high temperature gas in fermentation medium of different culture methods

分析圖3可知，參考文獻(xiàn)[4]方法和參考文獻(xiàn)[5]方法在不同實(shí)驗次數(shù)下，食用菌種植過程培養(yǎng)料發(fā)酵高溫氣體濃度去除率曲線波動比較大，且整體較低，可靠性較差。所提方法利用數(shù)值模擬的方式使用軸流風(fēng)機(jī)進(jìn)行食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)設(shè)計時，通過精細(xì)化的網(wǎng)格劃分，判斷軸流風(fēng)機(jī)疏散食用菌種植過程培養(yǎng)料發(fā)酵高溫氣體最佳參數(shù)，進(jìn)而最大程度降低了食用菌種植過程培養(yǎng)料發(fā)酵高溫氣體。

模擬情況與實(shí)際值擬合度表示研究方法考慮周全與否，該值越高，則表示提出方法魯棒性越強(qiáng)。由圖4可知，在進(jìn)行數(shù)值模擬前，所提方法設(shè)置了邊界條件，并且采用分區(qū)分塊方式，將整個幾何模型劃分成集流器、前后級動葉輪、前后級導(dǎo)葉輪與擴(kuò)散器幾個區(qū)域，考慮軸流風(fēng)機(jī)實(shí)際運(yùn)行狀況與食用菌種植過程培養(yǎng)料發(fā)酵環(huán)境，有效提升了模擬情況與實(shí)際值擬合度。

圖4 不同方法模擬情況與實(shí)際值擬合度對比Fig.4 Comparison of simulation results and actual values of different methods

3 結(jié)語

針對當(dāng)前食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)設(shè)計方法存在的問題，提出基于軸流風(fēng)機(jī)的食用菌種植過程培養(yǎng)料通風(fēng)發(fā)酵模擬方法。利用構(gòu)建軸流風(fēng)機(jī)控制數(shù)學(xué)模型和設(shè)置邊界條件為數(shù)值模擬奠定基礎(chǔ)，通過網(wǎng)格精細(xì)化劃分構(gòu)建軸流風(fēng)機(jī)有限元模型，完成軸流風(fēng)機(jī)食用菌種植過程培養(yǎng)料發(fā)酵通風(fēng)數(shù)值模擬。經(jīng)驗證，軸流風(fēng)機(jī)可有效去除食用菌種植過程培養(yǎng)料發(fā)酵高溫氣體，含氧量得到有效的補(bǔ)償，各方面性能良好。