張修霖,楊開敏*,王遠成,石天玉

淺圓倉壓入式和吸出式通風效果對比研究

張修霖1,楊開敏1*,王遠成1,石天玉2

1. 山東建筑大學熱能工程學院, 山東 濟南 250101 2. 國家糧食與物資儲備局科學研究院, 北京 100037

基于CFD數(shù)值模擬技術，對淺圓倉儲糧不同通風方式降溫下溫度的數(shù)值模擬結果進行了分析，對比了通風降溫的效果。文章研究的糧種為不易儲藏的大豆，數(shù)值模擬了壓入式通風和吸出式通風方式下，糧堆內溫度分布及速度分布的變化規(guī)律，對比四種通風條件下糧堆內部的溫度變化，分析了其均勻性及降溫速率。研究發(fā)現(xiàn)四種通風方式下淺圓倉糧堆區(qū)域的通風效果接近，通風均勻性較好；但吸出式（中心集風管頂部密閉）通風時，降溫速率更高，糧溫均勻性更好，研究結果對淺圓倉倉儲大豆通風降溫方式的選擇提供了理論指導和依據(jù)。

糧倉; 通風效果; 對比分析

淺圓倉是一種倉內裝糧高度與直徑之比小于1.5倍的地上立筒倉[1]，淺圓倉自上世紀后期引入中國，經(jīng)過二十多年的發(fā)展和研究，淺圓倉徑向通風技術相關研究已比較成熟，并被大規(guī)模推廣應用。與傳統(tǒng)平房倉相比，淺圓倉具有占地面積小、堆糧高度高、機械化作業(yè)水平高，倉儲物流費用低等特點，在糧食倉儲行業(yè)應用更為廣泛[2]。

大豆不耐高溫，過高的溫度會引起大豆的主要成分發(fā)生物理、化學和生物變化，例如蛋白質變性、脂肪分解等[3]，因此溫度是大豆儲藏過程中至關重要的影響因素。淺圓倉由于儲量大、糧堆高、自動分級嚴重，存在通風阻力大、通風效率低、能耗高，通風后溫度水分分層明顯，通風過程中水分丟失嚴重的問題。特別是在儲存進口大豆過程中，因大豆發(fā)熱升溫速度較快，對通風效率要求極高，如果短時內解決不了發(fā)熱問題極易造成板結、結塊，甚至走油赤變，嚴重時可能會出現(xiàn)碳化等現(xiàn)象[4]。因此對大豆進行通風降溫，將溫度控制在合理范圍內就可以安全儲存。

使用數(shù)值模擬方法求解糧堆內傳熱傳質問題是近年來的一種新方法[5,6]。目前，對通風時糧堆內部溫度變化的模擬計算主要使用Thorpe建立的數(shù)學模型[7]，該模型可以較好地反映糧堆內部溫度分布。Gadton A等[8]通過建立糧堆熱濕耦合模型研究糧食內部溫度變化。王遠成等[9]和高帥等[10]建立了儲糧通風模型，對糧堆內部空氣流動及熱濕耦合規(guī)律進行模擬研究，重點研究了房式倉垂直和橫向通風過程中溫度前沿的變化。戚禹康等[11,12]運用CFD技術，模擬分析了淺圓倉模擬分析了淺圓倉橫向和垂直通風兩種狀態(tài)下糧堆溫度場和流場的分布及變化趨勢。

盡管國內外對淺圓倉機械通風時溫度變化的模擬計算和實倉測試研究已經(jīng)取得了一定進展，但缺乏對徑向壓入式通風和吸出式通風的對比研究，也大都沒有考慮糧食儲存中初始溫度分布不均勻的情況。為此，本研究以大豆為研究對象，探究在相同通風量條件下，壓入式和吸出式徑向通風時糧堆內溫度場和速度分布的情況及變化趨勢，并通過分析不同通風方式對糧堆溫度均勻性、降溫速率等指標的影響。得到的結果可以為淺圓倉儲糧最優(yōu)通風方式的選擇提供理論依據(jù)。

1 淺圓倉物理模型及數(shù)學模型建立

以淺圓倉為研究對象，建立物理模型，直徑2 m，高3.3 m，裝糧高度為2.0 m，倉頂部有四個直徑口為0.1 m的風機出口。如圖1所示，淺圓倉徑向通風系統(tǒng)是由沿壁面上的16根支風道和中心集風管組成，支風道為緊貼倉壁的弧面，中心集風管為變徑開孔圓管。壓入式通風時，氣流由垂直支風道水平進入糧堆，并沿著半徑方向移動至糧倉中心的集風管，經(jīng)過集風管和糧面排到糧倉上部的空氣區(qū)域，然后通過倉頂?shù)念A留風機口排出倉外；吸出式通風期間氣流經(jīng)糧面及中心集風管進入糧堆，并沿圓周方向進入垂直支風道，最后匯入四周環(huán)形支風道，經(jīng)預留風機口排出倉外。

圖1 淺圓倉物理模型圖

圖2 淺圓倉徑向通風網(wǎng)格圖

1.1 淺圓倉壓入式通風和吸出式通風物理模型

本模型是根據(jù)試驗淺圓倉的實際尺寸通過ICEM建立了淺圓倉徑向壓入式通風和吸出式通風的數(shù)值模擬的物理模型。糧堆分為三個區(qū)域，即第一層為最下層，高度為1 m；第二層為中層，高度為0.5 m；第三層為最上層，高度為0.5 m，各層的初始溫度分布皆不相同。當壓入式通風時，支風道底部邊界條件為inlet，淺圓倉頂部四個離心機出口為outlet；當吸出式通風時，出入口邊界條件互換，即倉頂部離心機出口邊界條件為inlet，支風道底部為outlet。同時為了探究更優(yōu)通風方式，還建立了集風管頂部邊界條件為wall，其余條件和上述條件一致設立物理模型。然后對模型進行非結構網(wǎng)格劃分，如圖2所示。經(jīng)過網(wǎng)格獨立性驗證，最終確定網(wǎng)格數(shù)量約為37萬。

1.2 通風過程數(shù)學模型

1.2.1 連續(xù)性方程糧堆是一種典型的多孔介質結構，根據(jù)流體力學理論和質量守恒定律，糧堆上部空氣區(qū)域的連續(xù)性方程為：

式中一空氣密度，㎏·m-3；-時間，s；一氣體流速，m/s-1；▽一哈密頓微分算子。

在多孔介質結構模型中，根據(jù)谷物顆粒對氣流的影響，在流體的連續(xù)性方程中加入谷物顆粒的體積分數(shù)項[13]，即：

式中—氣體的體積分數(shù)項。

1.2.2 動量守恒方程空氣區(qū)域的流體的運動微分方程（N-S方程）為：

式中—氣體微元體上的壓強，Pa；—重力加速度，m·s-2；—氣體動力粘度，Pa·s。

因為給糧堆進行通風時，谷物顆粒會對氣流形成阻力，所以在流體對流流動及其阻力的動量方程中加入動量源項S，

式中動量源項S=-Rv-2。其中—氣體的動力粘度，Pa·s；是空氣的表觀速度或達西速度，m·s-1。

1.2.3 對流傳熱方程

1.3 邊界及初始條件

模擬研究對象為大豆，本次模擬噸糧通風量選用5.2 m3/(t·h)。送風溫度即進口溫度為15 ℃，糧堆初始平均溫度為22.75 ℃，糧堆區(qū)域溫度分層，最下層溫度為22 ℃、中層溫度為23 ℃、最上層溫度為24 ℃，空氣區(qū)域溫度為24 ℃，送風溫度與糧堆平均溫度溫差約為8 ℃。

2 模擬結果與分析

2.1 壓入式通風過程中糧堆溫度變化的模擬結果及分析

圖3和圖4為主風道頂部密閉及開孔壓入式通風時的氣流流線分布圖，從圖中可以看出兩種通風氣流方向十分相似，氣流由垂直支風道水平進入糧堆，并沿著半徑方向移動至糧倉中心的集風管，隨著高度升高，一部分氣流斜向上進入集風管，另一部分氣流垂直向上通過糧面排到倉體上部空氣區(qū)域，唯一區(qū)別是當集風管頂部密閉時，集風管內的氣流通過就近的糧面排至倉體上部空氣區(qū)域。在支風道及集風管中速度較大，糧堆內速度較小且均勻，無通風死角。

圖 3 壓入式通風（中心集風管頂部密閉）流線圖

圖4 壓入式通風（中心集風管頂部開孔）流線圖

圖5和圖6為壓入式通風時各層糧堆溫度及總糧堆平均溫度的變化曲線圖，可以看出各層糧堆溫度及總糧堆平均溫度都是隨著通風時間下降的，下層糧溫較低，上層糧溫較高，這是由于上層初始溫度就較下層初始溫度高，且低溫氣流首先通過最下層糧堆，會帶走一部分熱量，然后通過糧堆或集風管排出。通過對比圖5和圖6，可以看出兩種通風方式，各層糧堆及總糧堆平均溫度下降趨勢基本吻合，下邊兩層糧食及糧堆平均溫度在通風1~4 d內下降趨勢更為明顯，這是由于隨著通風時間增加，下邊兩層糧堆與低溫氣流間的溫差逐漸變小，降溫趨勢由前幾天的下降明顯也逐漸變緩；最上層糧堆溫度在通風3~6 d內下降趨勢更為明顯。

圖 5 壓入式通風（中心集風管頂部密閉）溫度變化

圖6 壓入式通風（中心集風管頂部開孔）溫度變化

圖7和圖8為主風道頂部密閉及開孔兩種情況下壓入式通風第5 d的溫度分布圖，由圖可以看出溫度分布與速度矢量相關，糧堆等溫線與風速方向幾乎成90°夾角，由倉壁沿著通風氣流方向糧溫逐漸增加，這是由于倉壁處糧溫先快速下降直至接近通風溫度。圖9和圖10為主風道頂部密閉及開孔2種情況下壓入式通風第11 d的溫度分布圖，表1為壓入式通風第11 d時各糧層溫度情況。

表1 壓入式通風第11 d時糧堆溫度

圖7 壓入式通風（頂部密閉）5 d溫度分布

圖8 壓入式通風（頂部開孔）5 d溫度分布

圖9 壓入式通風11 d（頂部密閉）溫度分布

圖10 壓入式通風11 d（頂部開孔）溫度分布

由圖和表中可以看出，通風11 d時，在上層糧堆靠近集風管處都會形成一小塊高溫區(qū)域，這是由于此處位于通風中氣流在糧堆最長流動路線的末端，同時該處阻力較大，風量較小，故溫度變化會有延遲。通過溫度分布圖可以看出，除去小塊高溫區(qū)域外，溫差都在1 ℃以內，糧堆的溫度也都均勻降低，說明通風均勻性較好。在通風第11 d時，上層糧堆最低糧溫接近15 ℃，已經(jīng)接近通風溫度，平均糧溫也降至16 ℃以下，但也在上層糧堆局部存在最高溫度16.37 ℃和16.25 ℃，從溫度變化趨勢圖亦可看出從第7~8 d糧溫下降趨勢已經(jīng)非常緩慢，這是由于低溫通風氣流和糧堆的溫度梯度已經(jīng)變得很小，而且上層糧堆大多處于通風氣流流動路線的末端。

2.2 吸出式通風過程中糧堆溫度變化的模擬結果及分析

圖11和圖12為主風道頂部密閉及開孔吸出式通風時的氣流流線分布圖，可以看出兩種通風方向十分相似，一部分氣流由糧堆表面和大部分氣流中心集風管進入糧堆，斜穿過糧層后進入倉壁四周的垂直支風道，然后排出倉外，區(qū)別是兩種通風方式產(chǎn)生渦流區(qū)域不同，且當集風管頂部密閉時，氣流會就近通過糧面進入中心集風管。通過流場可以看出，在支風道及集風管內速度較大，糧堆內速度較小且均勻，無通風死角。

圖13和圖14為吸出式通風時各層糧堆溫度及總糧堆平均溫度的變化曲線圖，可以看出各糧層溫度及糧堆平均溫度都是隨著通風時間下降的，但有部分糧堆平均溫度在通風初始會先升高然后再隨著通風時間下降，剛開始通風降溫時，下層糧溫較低，上層糧溫較高，這是由于上層糧食初始溫度就較下層糧食初始溫度高；隨著通風時間增加，上層糧溫較低，下層糧溫較高，這是由于低溫氣流從上層糧堆進入對糧堆進行降溫且會帶走一部分熱量經(jīng)過下層糧堆然后排出倉外。通過對比圖13和圖14，中層糧堆當中心集風管上部密閉時降溫趨勢比較明顯，與上層糧堆降溫趨勢趨于接近。上邊兩層糧食及糧堆平均溫度在通風1~5 d內下降趨勢更為明顯，這是由于隨著通風時間增加，上邊兩層糧食與低溫氣流溫度梯度逐漸變小，降溫趨勢由前幾天的下降明顯也逐漸變緩；最下層糧食溫度在通風降溫時下降趨勢比較緩慢，這是由于下層糧堆溫度與通過中上層糧堆的空氣溫度溫差較小，溫度梯度較小且通風氣流會攜帶熱量通過下層糧堆排出倉外。

圖11 吸出式通風（中心集風管頂部密閉）流線圖

圖12 吸出式通風（中心集風管頂部開孔）流線圖

圖13 吸出式通風（頂部密閉）溫度變化

圖14 吸出式通風（頂部開孔）溫度變化

圖15和圖16為主風道頂部密閉及開孔吸出式通風第5 d的溫度分布，由溫度分布圖可以看出溫度有明顯的分層現(xiàn)象，糧倉上部空氣區(qū)域及中心集風管附近溫度最低，水平方向糧堆溫度由中心集風管順著通風氣流方向逐漸增加，豎直方向遠離集風管靠近倉壁糧堆區(qū)域溫度變化不大，但糧堆頂部越靠近空氣區(qū)域溫度已經(jīng)有明顯降低。

圖17和圖18為主風道頂部密閉及開孔吸出通風第11 d的溫度分布圖，表2為由數(shù)值模擬導出的各層糧堆溫度。從中可以看出，通風11 d時，在下層糧堆靠近支風管處都會形成一小塊高溫區(qū)域，這是由于此處位于通風中氣流在糧堆最長流動路線的末端，同時該處阻力較大，風量較小，故溫度降低會有延遲。通過溫度分布圖可以看出，除了小塊高溫區(qū)域外，溫差都在1 ℃以內，糧堆的溫度也都均勻降低，說明通風均勻性較好。在通風第11 d時，下層糧堆最低糧溫為15 ℃，已經(jīng)接近通風溫度，平均糧溫也降至16 ℃以下，但也在下層糧堆存在最高溫度16.44 ℃和16.84 ℃，從溫度變化趨勢圖亦可看出從第7~8 d糧溫下降趨勢已經(jīng)非常緩慢，這是由于溫度梯度已經(jīng)變得很小，而且下層糧堆大多處于通風氣流流動路線的末端，溫度變化略有延遲。

表2 吸出式通風第11 d時糧堆溫度

圖 15 吸出式通風（頂部密閉）5 d溫度分布

圖 16 吸出式通風（頂部開孔）5 d溫度分布

圖 17 吸出式通風（頂部密閉）11 d溫度分布

圖 18 吸出式通風（頂部開孔）11 d溫度分布

2.3 壓入式和吸出式通風效果對比分析

由兩種工況計算結果的流線圖可以看出，4種通風方式均能夠達到近似的效果，在支風道及集風管速度較大，糧堆內速度較小且均勻，無通風死角。但通過對比4種通風方式對糧堆平均溫度變化的影響，由圖19可以看出當通風量相同時，隨著通風時間增加，糧堆平均溫度都呈相似的下降趨勢，都是在前5 d溫度下降明顯，5 d后下降趨勢逐漸平緩。通過對比不同工況下的糧堆平均溫度變化曲線，可以得出淺圓倉吸出式通風且中心集風管頂部密閉時的通風方式糧堆平均溫度下降趨勢最為明顯。

圖 19 四種通風方式糧堆平均溫度變化

3 結 論

當通風量為5.2 m3/(t·h)時，通過對淺圓倉徑向壓入式與吸出式通風進行數(shù)值模擬對比研究，得出以下結論。

（1）淺圓倉徑向壓入式與吸出式通風降溫的糧堆平均溫度變化趨勢基本一致，大致都在11天時降至到15 ℃；

（2）受空氣流動的影響，在通風降溫一段時間后，壓入式通風時糧堆上部存在小塊高溫區(qū)域，吸出式通風時在糧堆下部靠近支風道區(qū)域存在小塊高溫區(qū)域，由于其體積很小，對平均溫度基本沒有影響，且該區(qū)域與糧堆平均溫度不超過1 ℃，在允許的溫差范圍內，所以對糧食品質不會產(chǎn)生明顯影響；

（3）通過對4種通風降溫方式對糧堆降溫效果進行模擬，并進行對比分析，得出當通風量相同時，淺圓倉吸出式通風且中心集風管頂部密閉時的通風方式糧堆平均溫度比其余三種通風方式下降較快，通風效果更優(yōu)；

（4）通過進行各糧層及糧堆最高溫度的對比分析，可以得出淺圓倉壓入式通風時，糧堆除最上層糧堆還有部分高溫區(qū)域，其余糧層溫度都降至通風氣流溫度；淺圓倉吸出式通風時，糧堆除最下層糧堆還有部分高溫區(qū)域，其余糧層溫度都降至通風氣流溫度。通過計算，淺圓倉吸出式通風且中心集風管頂部密閉時的通風方式，糧溫下降較快且糧溫均勻性較好。綜上所述，當淺圓倉徑向通風對儲藏大豆進行降溫時，設立中心集風管頂部密閉吸出式通風時降溫效果更優(yōu)。

[1] 劉朝偉,呂建華,楊冰,等.淺圓倉負壓通風技術應用案例分析[J].糧食與食品工業(yè),2018,25(5):55-58

[2] 張來林,蔡育池,許國川,等.淺談淺圓倉的儲糧特點[J].糧油食品科技,2019,27(5):97-100

[3] 張永君,韓偉.大豆安全儲藏的應用技術[J].糧食與食品工業(yè),2010,17(5):48-50

[4] 白巖,曹毅,崔國華,等.大豆安全儲藏特性分析與工藝探討[J].糧食加工,2005(4):57-60

[5] 白忠權.吸濕性倉儲糧堆內熱濕耦合傳遞規(guī)律的研究[D].濟南:山東建筑大學,2013

[6] 吳銳.糧堆溫度場實測結果與數(shù)值模擬分析[D].鄭州:河南工業(yè)大學,2011

[7] Thorpe GR. Modelling ecosystems in ventilated conical bottomed farm grain silos [J]. Ecological Modelling, 1997,94(2-3):255-286

[8] Gadton A, Abalone R, Batosik RE,. Mathematical modelling of heat and moisture transfer of wheat stored in plastic bags ( silobags) [J]. Biosystems Engineering, 2009,104(1):72-85

[9] 王遠成,高帥,邱華宇,等.橫向谷冷通風過程的數(shù)值模擬研究[J].中國糧油學報,2016,31(7):103-106

[10] 高帥,王遠成,趙會義,等.平房倉橫向谷冷通風小麥糧堆傳熱傳質數(shù)值模擬[J].糧油食品科技,2015,23:15-19

[11] 戚禹康,王遠成,俞曉靜,等.淺圓倉不同裝糧高度時徑向與垂直通風模擬對比研究[J].河南工業(yè)大學學報(自然科學版),2019,40(4):108-113

[12] 戚禹康,王遠成,俞曉靜.淺圓倉橫向和垂直通風數(shù)值模擬對比研究[J].山東建筑大學學報,2019,34(3):50-56

[13] 郝世楊.稻谷通風阻力數(shù)值模擬與試驗研究[D].哈爾濱:東北農(nóng)業(yè)大學,2015

Research on Comparision between Ventilative Effects of Press-in and Suction in Shallow Round Grain Store

ZHANG Xiu-lin1, YANG Kai-min1*, WANG Yuan-cheng1, SHI Tian-yu2

1250101,2.100037,

In this paper, based on the CFD numerical simulation technology,the numerical simulation results of the temperature under different ventilation modes of shallow-circle storage grains were analyzed, and the effects of ventilation and cooling were compared. The grains studied in this article are soybeans that are not easy to store. Numerical simulations of the change law of temperature distribution and velocity distribution in the grain pile under the blowing ventilation and exhaust system of ventilation methods, compares the temperature change inside the grain pile under four ventilation conditions, and analyze its uniformity and cooling rate. The study found that the ventilation effect of the shallow-circle silo grain pile area under the four ventilation modes is similar, and the ventilation uniformity is better. However, the cooling rate is higher and the grain temperature uniformity is better when the suction type (the top of the central air collector is closed) is ventilated. The research results provide theoretical guidance and basis for the selection of ventilation and cooling methods for soybeans in shallow-circle silo storage.

Grain store; ventilative effect; comparative analysis

X954,X956

A

1000-2324(2021)01-0063-07

10.3969/j.issn.1000-2324.2021.01.011

2020-05-26

2020-06-05

國家重點研發(fā)項目(2016YDF0401002,2016YFD0400104)

張修霖(1996-),男,碩士研究生,主要研究方向為動力粘度. E-mail:1315659455@qq.com

Author for correspondence. E-mail:ykm_168@163.com