隨著我國(guó)工業(yè)化水平的快速發(fā)展，大多數(shù)傳統(tǒng)的農(nóng)業(yè)生產(chǎn)方式已被現(xiàn)代化機(jī)械所替代[1-2]。傳統(tǒng)的糧食收儲(chǔ)模式存在著收儲(chǔ)周期長(zhǎng)、產(chǎn)后損失嚴(yán)重以及糧食質(zhì)量較低等問(wèn)題，不能滿足現(xiàn)代化需求[3]。糧食烘干作為收儲(chǔ)過(guò)程的重要一環(huán)，對(duì)于糧食產(chǎn)后收儲(chǔ)減損有著至關(guān)重要的作用。于洋等[4]采用Fluent對(duì)烘干機(jī)烘干室的氣流分布進(jìn)行分析，解決了枸杞烘干過(guò)程中存在的粘黏問(wèn)題。于省元等[5]針對(duì)烘干機(jī)的排糧機(jī)構(gòu)進(jìn)行創(chuàng)新優(yōu)化，解決了烘干機(jī)排糧時(shí)存在的混糧混種現(xiàn)象。王曉明等[6]設(shè)計(jì)了塔式烘干機(jī)的氣力輸送系統(tǒng)，實(shí)現(xiàn)了谷物的定量供給和防擠壓破碎。本文對(duì)烘干塔干燥段進(jìn)行分析，為后續(xù)烘干裝置的優(yōu)化分析提供參考。

1 干燥段流場(chǎng)數(shù)值模擬設(shè)計(jì)方案

1.1 干燥段物理模型及仿真參數(shù)

烘干塔干燥段物理模型如圖1所示。熱風(fēng)從進(jìn)風(fēng)口以一定風(fēng)速充入進(jìn)風(fēng)空腔內(nèi)，熱空氣沿進(jìn)風(fēng)角狀盒逐漸充滿整個(gè)干燥段物料層，最后從出風(fēng)角狀盒到出風(fēng)口流出。

圖1 烘干塔干燥段物理模型圖

烘干塔的干燥效率與其內(nèi)溫度場(chǎng)和速度場(chǎng)均勻度有很大關(guān)系，為對(duì)干燥段內(nèi)溫度和速度分布進(jìn)行量化分析，并評(píng)價(jià)干燥段內(nèi)溫度場(chǎng)和速度場(chǎng)分布的均勻度，在每個(gè)角狀盒內(nèi)選擇從進(jìn)口處到出口處間隔200 mm的8個(gè)監(jiān)測(cè)點(diǎn)（計(jì)算域內(nèi)共88個(gè)監(jiān)測(cè)點(diǎn)），監(jiān)測(cè)點(diǎn)分布如圖2所示，其中in1～in11分別表示角狀盒進(jìn)氣口編號(hào)，in1距離干燥段進(jìn)氣口最近，in11距離干燥段進(jìn)氣口最遠(yuǎn)。

圖2 監(jiān)測(cè)點(diǎn)分布圖

1.2 模型構(gòu)建及網(wǎng)格劃分

利用SolidWorks三維建模軟件對(duì)干燥段外殼進(jìn)行建模并在Space Claim中抽取內(nèi)部流場(chǎng)空間，在FLUENT Meshing中進(jìn)行網(wǎng)格劃分后導(dǎo)入ANSYS FLUENT中進(jìn)行仿真模擬計(jì)算。干燥段尺寸和進(jìn)出風(fēng)角狀盒尺寸如圖3所示。干燥段內(nèi)流場(chǎng)網(wǎng)格的劃分如圖4所示。因模型較為復(fù)雜，結(jié)構(gòu)化網(wǎng)格效果欠佳，故采用擁有更高的質(zhì)量和計(jì)算精度的poly網(wǎng)格。

圖3 干燥段主視圖與進(jìn)出風(fēng)角狀盒尺寸圖

圖4 網(wǎng)格劃分示意圖

1.3 仿真參數(shù)設(shè)計(jì)

為了探究干燥段內(nèi)溫度場(chǎng)和速度場(chǎng)分布的影響因素，針對(duì)不同進(jìn)風(fēng)風(fēng)速和不同模型結(jié)構(gòu)進(jìn)行仿真分析。在原始模型的基礎(chǔ)上，改變進(jìn)風(fēng)風(fēng)速為10 m·s-1、15 m·s-1、20 m·s-1、30 m·s-1和40 m·s-1進(jìn)行數(shù)值模擬分析。

2 原始模型仿真模擬結(jié)果及分析

2.1 干燥段內(nèi)溫度場(chǎng)和速度場(chǎng)的分布

如圖5（a）所示，在進(jìn)口風(fēng)速為20 m·s-1時(shí)，空氣在進(jìn)風(fēng)空腔和兩側(cè)進(jìn)氣角狀盒內(nèi)溫度高于周?chē)锪蠈硬糠郑译S著熱風(fēng)的充入，距離進(jìn)風(fēng)空腔越遠(yuǎn)的位置溫度越低。將溫度監(jiān)測(cè)點(diǎn)選取在進(jìn)風(fēng)口左側(cè)干燥段，比較各角狀盒內(nèi)溫度，如圖5（b）所示。由圖可知角狀盒in1內(nèi)的溫度低于其他角狀盒，且溫度不均勻度為4.66%，高于其他角狀盒的溫度不均勻度。

圖5 干燥段內(nèi)部溫度分布變化情況圖

（a）Z=820 mm剖面的溫度分布云圖

如圖6（a）所示，在進(jìn)口風(fēng)速為20 m·s-1時(shí)，進(jìn)風(fēng)口處和角狀盒內(nèi)的速度高于周?chē)锪蠈拥乃俣?，且距離進(jìn)風(fēng)空腔越遠(yuǎn)，角狀盒內(nèi)的風(fēng)速越低。圖6（b）給出不同角狀盒內(nèi)各監(jiān)測(cè)點(diǎn)的速度變化曲線圖。從圖中可以看出，對(duì)于除in1以外的其他角狀盒，監(jiān)測(cè)點(diǎn)7的速度均低于角狀盒內(nèi)的其他監(jiān)測(cè)點(diǎn)，且角狀盒in1的速度低于其他角狀盒。通過(guò)分析干燥段內(nèi)速度分布情況，可以解釋其溫度分布情況。因角狀盒in1內(nèi)的風(fēng)速較低，導(dǎo)致干燥段內(nèi)熱風(fēng)流動(dòng)不充分，溫度場(chǎng)和速度場(chǎng)不均勻，且不均勻位置位于角狀盒in1和監(jiān)測(cè)點(diǎn)7。

圖6 干燥段內(nèi)部速度分布變化情況圖

2.2 進(jìn)口風(fēng)速對(duì)溫度場(chǎng)和速度場(chǎng)的影響

不同風(fēng)速下的角狀盒in1內(nèi)各監(jiān)測(cè)點(diǎn)和監(jiān)測(cè)點(diǎn)7的溫度變化曲線圖如圖7、圖8所示。

圖7 不同風(fēng)速下角狀盒in1內(nèi)各監(jiān)測(cè)點(diǎn)溫度變化曲線圖

圖8 不同風(fēng)速下監(jiān)測(cè)點(diǎn)7溫度變化曲線圖

由圖7可知，在進(jìn)口風(fēng)速為10 m·s-1的情況下，角狀盒in1內(nèi)溫度最低，在進(jìn)口風(fēng)速為30 m·s-1的情況下，角狀盒in1內(nèi)的溫度升高，不均勻度從4.35%降低到4%。表明改變進(jìn)口風(fēng)速能改進(jìn)角狀盒in1內(nèi)溫度的不均勻性，在進(jìn)口風(fēng)速30 m·s-1時(shí)效果最好，溫度增幅量在3 ℃左右，不均勻度降幅0.35%。

由圖8可知，增大進(jìn)口風(fēng)速在大部分角狀盒內(nèi)監(jiān)測(cè)點(diǎn)7位置都能提升溫度，最大溫度增幅達(dá)到4 ℃。在進(jìn)口風(fēng)速為40 m·s-1的情況下，除了in6、in7以外，監(jiān)測(cè)點(diǎn)7的溫度都高于其他低進(jìn)口風(fēng)速。因此，改變進(jìn)口風(fēng)速在一定程度上能改善監(jiān)測(cè)點(diǎn)7處溫度較低的情況，但對(duì)提升監(jiān)測(cè)點(diǎn)6的溫度效果不明顯。在進(jìn)口風(fēng)速為40 m·/s-1時(shí)，不均勻度降低至4%。

干燥段內(nèi)流場(chǎng)的平均溫度如圖9所示。由圖9可知，隨著進(jìn)口風(fēng)速的增加，干燥段內(nèi)流場(chǎng)的平均溫度逐漸上升，從進(jìn)口風(fēng)速10 m·s-1到40 m·s-1，平均溫度增幅達(dá)到3 ℃。這說(shuō)明，增大進(jìn)口風(fēng)速可提高干燥段內(nèi)的平均溫度。

結(jié)合圖7～圖9可知，增大進(jìn)口風(fēng)速能對(duì)干燥段內(nèi)的不均勻區(qū)域包括角狀盒in1和監(jiān)測(cè)點(diǎn)7進(jìn)行改善，在進(jìn)口風(fēng)速為40 m·s-1時(shí)，不均勻度降低至4%。

圖9 不同進(jìn)口風(fēng)速下干燥段內(nèi)流場(chǎng)的平均溫度圖

3 改進(jìn)模型仿真模擬結(jié)果及分析

3.1 改進(jìn)模型

本文主要采用3種改進(jìn)模型同原始模型N進(jìn)行對(duì)比，改進(jìn)模型分別為BISO模型、SIBO模型和CROS模型，分別表示增大角狀盒進(jìn)口減小角狀盒出口、減小角狀盒進(jìn)口增大角狀盒出口和大小角狀盒分布排列，3種改進(jìn)模型的形狀變化示意圖如圖10所示。

圖10 改進(jìn)模型結(jié)構(gòu)示意圖

通過(guò)前文對(duì)不同進(jìn)口風(fēng)速下干燥段內(nèi)部溫度場(chǎng)和速度場(chǎng)的分析可知，在進(jìn)口風(fēng)速為40 m·s-1時(shí)，干燥段內(nèi)不利區(qū)域的改進(jìn)效果最好，因此以上3種改進(jìn)模型均在進(jìn)口風(fēng)速為40 m·s-1的條件下進(jìn)行模擬仿真。通過(guò)在相同條件下的數(shù)值分析，得出改變角狀盒進(jìn)出口大小對(duì)干燥段內(nèi)部溫度場(chǎng)和速度場(chǎng)的影響，從而分析改變干燥段結(jié)構(gòu)對(duì)其內(nèi)部流場(chǎng)不均勻度的影響。

3.2 不同模型的溫度場(chǎng)和速度場(chǎng)變化

圖11給出不同改進(jìn)模型和原始模型N在角狀盒in1內(nèi)的溫度變化曲線圖。對(duì)比發(fā)現(xiàn)，3種改進(jìn)模型能夠提升角狀盒in1內(nèi)部流場(chǎng)大部分監(jiān)測(cè)點(diǎn)的溫度，溫度增幅最低達(dá)到1 ℃。其中SIBO模型溫度提升效果最好，不均勻度從N模型的4.66%下降到4%。但CROS模型在監(jiān)測(cè)點(diǎn)6的位置溫度反而下降，導(dǎo)致in1角狀盒內(nèi)不均勻度增大。

圖11 不同模型在角狀盒in1內(nèi)的溫度變化曲線圖

圖12給出不同改進(jìn)模型和原始模型N之間的干燥段內(nèi)平均溫度變化曲線圖?？梢钥闯?，3種改進(jìn)模型的平均溫度均高于原始模型N的平均溫度，但溫度增幅在1 ℃以內(nèi)，說(shuō)明改進(jìn)模型對(duì)提高干燥段內(nèi)流場(chǎng)的平均溫度效果并不明顯。此外，由于改變角狀盒大小會(huì)造成干燥段內(nèi)糧食流動(dòng)發(fā)生變化，在角狀盒大的部分糧食流動(dòng)慢，角狀盒小的部分糧食流動(dòng)快，導(dǎo)致糧食在干燥段內(nèi)可能出現(xiàn)受熱不均勻現(xiàn)象，反而會(huì)降低干燥段的干燥效率。

圖12 不同模型干燥段內(nèi)流場(chǎng)的平均溫度圖

綜上可知，3種改進(jìn)模型在一定程度上確實(shí)能夠改善角狀盒in1的不利區(qū)域，提升其內(nèi)的溫度，降低其內(nèi)的溫度不均勻度。但3種改進(jìn)模型對(duì)于監(jiān)測(cè)點(diǎn)7和流場(chǎng)內(nèi)平均溫度的改善效果存在一定限制，需依據(jù)角狀盒的改變而導(dǎo)致的糧食流動(dòng)不均的問(wèn)題進(jìn)行進(jìn)一步的分析，使其達(dá)到最佳干燥效果。

4 結(jié)論

本文借助ANSYS FLUENT軟件建立烘干塔干燥段的簡(jiǎn)化模型，采用CFD的標(biāo)準(zhǔn)湍流模型和傳熱模型進(jìn)行數(shù)值模擬，對(duì)干燥段內(nèi)溫度場(chǎng)和速度場(chǎng)進(jìn)行分析，同時(shí)通過(guò)改變進(jìn)口熱風(fēng)的風(fēng)速和角狀盒的大小來(lái)判斷影響流場(chǎng)內(nèi)溫度和速度及其不均勻度的因素，結(jié)果如下。

（1）原始模型溫度場(chǎng)和速度場(chǎng)的不均勻區(qū)域主要集中在角狀盒in1和監(jiān)測(cè)點(diǎn)7位置。角狀盒in1和監(jiān)測(cè)點(diǎn)7處速度相對(duì)其他角狀盒較低，熱空氣流動(dòng)不充分，溫度場(chǎng)和速度場(chǎng)不均勻，導(dǎo)致周?chē)锪蠈邮軣岵怀浞郑档透稍锒蔚倪\(yùn)行效率。

（2）在原始模型的基礎(chǔ)上，通過(guò)改變進(jìn)口風(fēng)速的大小，得到增加進(jìn)口風(fēng)速能改善干燥段內(nèi)的不均勻區(qū)域。在進(jìn)口風(fēng)速為40 m·s-1的情況下，相較于10 m·s-1的情況，角狀盒in1和監(jiān)測(cè)點(diǎn)7的溫度增幅達(dá)到3 ℃左右，不均勻度從4.35%降低到4%。

（3）通過(guò)SIBO模型在提升角狀盒in1內(nèi)溫度來(lái)改善不均勻區(qū)域有一定的優(yōu)化效果，但結(jié)合監(jiān)測(cè)點(diǎn)7溫度變化和流場(chǎng)內(nèi)平均溫度和不均勻度變化，此方法存在一定的限制，因此后續(xù)可針對(duì)受熱均勻性、糧食流動(dòng)均勻性以及進(jìn)風(fēng)溫度之間的最佳配合比進(jìn)行深入的研究分析。