王遠(yuǎn)成 高 帥 邱化禹 魏 雷 趙會(huì)義 曹 陽(yáng)

橫向谷冷通風(fēng)過(guò)程的數(shù)值模擬研究

王遠(yuǎn)成1高 帥1邱化禹1魏 雷2趙會(huì)義2曹 陽(yáng)2

（山東建筑大學(xué)熱能工程學(xué)院1，濟(jì)南 250101）
（國(guó)家糧食局科學(xué)研究院2，北京 100037）

基于局部熱濕平衡原理和多孔介質(zhì)傳熱傳質(zhì)理論，建立了儲(chǔ)糧通風(fēng)過(guò)程中糧堆內(nèi)部流動(dòng)及熱濕耦合傳遞的數(shù)學(xué)模型。采用計(jì)算流體動(dòng)力學(xué)的方法，對(duì)橫向谷冷通風(fēng)時(shí)糧堆空氣內(nèi)部流動(dòng)、熱量傳遞和水分遷移過(guò)程進(jìn)行了數(shù)值分析。研究發(fā)現(xiàn)，橫向通風(fēng)約72 h，糧倉(cāng)進(jìn)風(fēng)口冷空氣平均溫度為17.5℃、相對(duì)濕度為85%，倉(cāng)內(nèi)糧堆一次降溫從32.2℃降低到23.6℃，降溫幅度為8.6℃。水分從12.2%降到12.0%，降水幅度為0.2%。相對(duì)于地上籠垂直通風(fēng)而言，橫向通風(fēng)時(shí)糧堆內(nèi)部速度分布均勻、溫度梯度較小，且具有降溫速度快，冷卻效率高的特點(diǎn)。

數(shù)值模擬 橫向通風(fēng) 傳熱傳質(zhì) 儲(chǔ)糧

糧倉(cāng)機(jī)械通風(fēng)在糧食儲(chǔ)藏中具有降溫效果顯著、費(fèi)用較低等特點(diǎn)，在確保儲(chǔ)糧安全方面，發(fā)揮著至關(guān)重要的作用［1］。目前國(guó)內(nèi)常用的機(jī)械通風(fēng)方式大多是地上籠垂直通風(fēng)，但這種通風(fēng)方式給糧食出入庫(kù)帶來(lái)極大的不便，加大了糧食進(jìn)出庫(kù)的工作量，增加了儲(chǔ)糧成本。橫向通風(fēng)是把2組通風(fēng)籠垂直安裝在糧倉(cāng)寬度方向的2個(gè)內(nèi)墻上，通過(guò)吸式或吹式的方式實(shí)現(xiàn)了沿著糧倉(cāng)寬度方向的橫向通風(fēng)。橫向通風(fēng)避免了地上通風(fēng)籠所帶來(lái)的糧食進(jìn)出倉(cāng)的不便，使得糧食的裝卸可以在糧倉(cāng)內(nèi)部完成，大大地提高了糧食進(jìn)出倉(cāng)的效率。

基于計(jì)算流體動(dòng)力學(xué)的數(shù)值模擬方法是國(guó)外近年發(fā)展起來(lái)的一種研究流動(dòng)、傳熱傳質(zhì)等現(xiàn)象的新方法，它可以形象地再現(xiàn)流動(dòng)、熱濕傳遞過(guò)程的情景，為解決儲(chǔ)糧通風(fēng)問(wèn)題也提供了一個(gè)良好的數(shù)值分析和優(yōu)化設(shè)計(jì)的工具［2］。

基于計(jì)算流體動(dòng)力學(xué)的方法，對(duì)橫向分布式谷冷通風(fēng)時(shí)糧堆內(nèi)部流動(dòng)和熱濕耦合傳遞過(guò)程進(jìn)行了數(shù)值模擬研究，探究了橫向通風(fēng)時(shí)糧堆內(nèi)部溫度和水分的變化規(guī)律，分析評(píng)價(jià)了橫向通風(fēng)時(shí)糧堆降溫的均勻性和降溫效果，并通過(guò)與試驗(yàn)結(jié)果的比較驗(yàn)證了數(shù)學(xué)模型。

1 糧堆內(nèi)部流動(dòng)和熱濕耦合傳遞過(guò)程的數(shù)學(xué)模型

假設(shè)糧堆是連續(xù)性的、均勻分布的多孔介質(zhì)，糧堆內(nèi)部滿足局部熱濕平衡原理［3－5］，考慮糧食顆粒的吸濕和解吸濕特性，由于通風(fēng)時(shí)間相對(duì)較短，谷物和蟲霉呼吸作用所產(chǎn)生的熱量和水分可以忽略不計(jì)。糧堆內(nèi)部流動(dòng)及熱濕耦合傳遞的控制方程［6－8］：

連續(xù)性方程

動(dòng)量方程

能量方程

水分遷移方程

式中：ε為空隙率；ρa(bǔ)為空氣密度；ρb為糧堆的容重；dp為谷物顆粒的等效直徑；u為糧堆內(nèi)部空氣的表觀速度或達(dá)西速度；p為壓力；t為時(shí)間；?為微分算子；ca、cb分別為空氣和糧堆的比熱；T為糧堆絕對(duì)溫度；W為糧堆的水分；keff為糧堆的有效導(dǎo)熱系數(shù)；μ為空氣的動(dòng)力黏度；w為糧粒間空氣中的絕對(duì)含濕量；Deff為濕空氣在糧堆中的有效擴(kuò)散系數(shù)；hs為谷粒吸濕或解吸濕熱。

2 數(shù)值方法及初始和邊界條件

數(shù)值模擬的糧倉(cāng)為房式倉(cāng)，糧倉(cāng)長(zhǎng)度為60 m，寬度為21 m，裝糧高度為5.8 m，儲(chǔ)存5 700 t小麥。為了便于分析，本研究取出長(zhǎng)度方向一半的中截面上的溫度和水分變化進(jìn)行分析和討論，如圖1所示。

圖1 房式倉(cāng)橫向通風(fēng)系統(tǒng)和中截面示意圖

2.1 數(shù)值方法

對(duì)控制方程采用有限體積法進(jìn)行離散，離散格式為二階上風(fēng)差分格式。為了防止迭代過(guò)程的發(fā)散和數(shù)值不穩(wěn)定，對(duì)動(dòng)量方程、能量和標(biāo)量輸運(yùn)方程采用了欠松弛技術(shù)，壓力與速度耦合采用SIMPLE算法。

2.2 初始條件及熱物性參數(shù)

糧堆（小麥）平均初溫為32.2℃（絕對(duì)溫度為305.2 K），平均濕基水分為12.2%（干基水分為13.98%）。糧堆的容重為ρb＝ 805 kg／m3，糧堆的比熱容為cp＝1 790 J／kg·K，糧堆的導(dǎo)熱系數(shù)為λ＝0．159 W／m·K ，糧堆的孔隙度ε 取為0.4［6－8］。

2.3 邊界條件

出口邊界沿流線方向各流動(dòng)參數(shù)的一階導(dǎo)數(shù)取為零，在固體壁面和地面采用無(wú)滑移條件，并且為絕熱邊界，由于糧面覆膜，糧面假設(shè)不可滲透，但與上部空氣有熱量交換。由于采用的是橫向分布式谷冷通風(fēng)，即將谷冷機(jī)作為冷源放在糧倉(cāng)的一側(cè)，并通過(guò)糧倉(cāng)另一側(cè)的引風(fēng)機(jī)把谷冷機(jī)產(chǎn)生冷風(fēng)引入糧倉(cāng)內(nèi)部。為了避免通風(fēng)過(guò)程中糧堆水分損耗過(guò)大，一方面選取較小的通風(fēng)量；另一方面，調(diào)節(jié)谷冷機(jī)出口空氣的濕度，谷冷機(jī)出風(fēng)口的空氣保持在較高濕度。因此，在實(shí)際操作中，糧倉(cāng)進(jìn)風(fēng)口冷空氣平均溫度為17.5℃，糧倉(cāng)進(jìn)風(fēng)相對(duì)濕度為85%，單位通風(fēng)風(fēng)量為q ＝4.9 m3／（h·t）。

3 模擬結(jié)果與分析

模擬通風(fēng)時(shí)間為72 h，糧堆平均表觀風(fēng)速為0.023 m／s，由于流動(dòng)速度很小，糧堆內(nèi)的流動(dòng)主要是靜壓驅(qū)動(dòng)。圖2是橫向通風(fēng)時(shí)糧堆內(nèi)部的速度分布，其中圖2上部三角形區(qū)域?yàn)槟ど峡臻g的空氣自然對(duì)流運(yùn)動(dòng)，下部矩形區(qū)域?yàn)榧Z堆內(nèi)部流動(dòng)的流場(chǎng)。從圖2可以看出，橫向通風(fēng)時(shí)糧堆內(nèi)部的流場(chǎng)分布非常均勻，這給糧堆的快速均勻降溫奠定了基礎(chǔ)。

圖2 橫向通風(fēng)在糧堆內(nèi)的速度分布

圖3 通風(fēng)72 h時(shí)糧堆溫度分布

圖4 通風(fēng)72 h時(shí)糧堆水分分布

圖5 不考慮吸濕／解吸濕作用通風(fēng)72 h時(shí)糧堆溫度分布

圖3 和圖4是通風(fēng)72 h時(shí)糧堆內(nèi)部溫度（開爾文溫度）和水分分布。從圖3和圖4可以看出，由于通風(fēng)時(shí)進(jìn)入糧堆的空氣溫度低于糧堆的初始溫度，隨著冷風(fēng)前沿的推進(jìn)，糧粒與周圍空氣之間進(jìn)行對(duì)流換熱，糧堆與冷空氣進(jìn)行熱量交換而降低溫度。在糧堆與空氣進(jìn)行熱量交換的同時(shí)，糧粒還與周圍空氣進(jìn)行質(zhì)量交換，從而導(dǎo)致糧堆水分發(fā)生遷移。從圖3還可以看出，糧堆內(nèi)部溫度前沿（冷鋒面）是波浪式的，即糧堆內(nèi)部存在著溫度不同的多個(gè)冷鋒面，而不是只有1個(gè)冷鋒。為了說(shuō)明這個(gè)問(wèn)題，數(shù)值計(jì)算了不考慮水分變化，即只考慮通風(fēng)時(shí)糧堆與空氣之間的熱量交換，不考慮糧堆與空氣之間的質(zhì)量交換。也即僅僅求解流動(dòng)和傳熱，并且忽略方程3右邊第二項(xiàng)（吸濕或解吸濕熱）。此時(shí)糧堆內(nèi)部溫度隨著時(shí)間變化規(guī)律如圖5所示，比較圖3和圖5可以發(fā)現(xiàn)二者之間差別較大，前者降溫較快且存在多個(gè)冷鋒面，而后者的降溫較慢且只有1個(gè)冷鋒，主要原因是后者沒有考慮糧食吸濕性特征，忽略了糧堆的解吸濕作用，此時(shí)糧堆類同于沙堆，其內(nèi)部水分被認(rèn)為是恒定不變的，這顯然與事實(shí)不符。由于在糧堆水分蒸發(fā)（解吸濕）而降低時(shí)，水分蒸發(fā)潛熱沒有被考慮，而水分的蒸發(fā)必然消耗糧堆的熱力學(xué)能量，使得溫度降低加快。因此，在計(jì)算糧堆內(nèi)部熱濕傳遞過(guò)程時(shí)，必須要考慮糧堆內(nèi)部質(zhì)量交換，即糧粒與糧粒間空氣的水蒸氣交換，否則模擬計(jì)算的結(jié)果會(huì)產(chǎn)生較大的誤差，甚至產(chǎn)生謬誤。

圖6 糧堆平均溫度隨時(shí)間的變化規(guī)律

圖7 糧堆平均水分隨時(shí)間的變化規(guī)律

圖6 和圖7是糧堆平均溫度和水分的時(shí)間變化圖，從圖3和圖6看出，通風(fēng)約72 h，倉(cāng)內(nèi)糧堆溫度從32.2℃降低到23.6℃，降溫幅度為8.6℃。同時(shí)，從圖6還可以看出，溫度的數(shù)值模擬與試驗(yàn)測(cè)定結(jié)果相符，說(shuō)明本研究建立的數(shù)學(xué)模型可以用于儲(chǔ)糧通風(fēng)過(guò)程中糧堆內(nèi)部熱濕傳遞的模擬研究。從圖4和圖7看出，通風(fēng)72 h，糧堆平均水分從12.2%降到12.0%，降水幅度為0.2%，水分降低很小。這因?yàn)楸狙芯磕M的是分布式谷冷通風(fēng)，進(jìn)風(fēng)口的空氣濕度為85%，因此，通風(fēng)過(guò)程中糧堆的水分降低較少損耗較小。

圖8和圖9是寬度方向各個(gè)立面上和高度方向各個(gè)水平面上（見圖1所示）糧食溫度隨時(shí)間變化圖。從圖中可以看到，在糧堆寬度方向上，冷風(fēng)從進(jìn)口糧堆的北1立面推進(jìn)到冷風(fēng)出口的南1立面，冷鋒面逐步推進(jìn)，糧溫降低。在間距為5 m的相鄰2個(gè)立面之間的平均糧溫相差不大于2℃，北1和南1 2個(gè)立面的溫差為4.6℃。從圖8還可以看出，橫向谷冷通風(fēng)時(shí)糧堆的降溫速度快，平均每3.6 h冷風(fēng)向前推進(jìn)1 m（以通風(fēng)結(jié)束時(shí)刻南1立面上的平均溫度為準(zhǔn)），冷鋒面移動(dòng)速度為0.28 m／h。從圖9可以看出，在糧堆高度方向，除了頂層以外，相鄰的2層之間糧堆的溫差不超過(guò)0.2℃。頂層糧溫也有明顯降低，較空氣區(qū)域平均溫度低約0.5℃，頂層與底層之間的溫差為0.4℃。由此可以看出橫向谷冷通風(fēng)時(shí)糧堆高度方向溫度梯度較小，沒有出現(xiàn)地上籠垂直通風(fēng)造成的沿著高度方向溫度梯度過(guò)大的問(wèn)題。

圖8 寬度方向各立面上溫度變化規(guī)律

圖9 高度方向各水平面上溫度變化規(guī)律

4 結(jié)論

4.1 橫向谷冷通風(fēng)時(shí)糧堆內(nèi)部表觀風(fēng)速為0.023 m／s，冷鋒前沿移動(dòng)速度為0.28 m／h。橫向通風(fēng)約72 h，倉(cāng)內(nèi)糧堆一次降溫從32.2℃降低到23.6℃，降溫幅度為8.6℃。水分從12.2%降到12.0%，降水幅度為0.2%。橫向通風(fēng)時(shí)糧堆內(nèi)部速度分布均勻，且具有降溫速度快，冷卻效率高的特點(diǎn)。

4.2 橫向通風(fēng)時(shí)在糧堆寬度方向上，相鄰的間距為5 m的相鄰2個(gè)截面之間的平均糧溫變化不大于2℃；糧堆高度方向水平上，各層之間溫差不超過(guò)0.2℃。相對(duì)于地上籠垂直通風(fēng)而言，橫向谷冷通風(fēng)時(shí)糧堆內(nèi)部溫度梯度較小。

［1］王遠(yuǎn)成，魏雷，劉偉，等.儲(chǔ)糧保水降溫通風(fēng)關(guān)鍵技術(shù)研究［J］.中國(guó)糧油學(xué)報(bào)，2008，23（5）：141－145

Wang Y C，Wei L，Liu W，et al.Study on key techniques of ventilation for cooling stored cereal grains ［J］.Journal of the Chinese Cereals and Oils Association，2008，23（5）：141 －145

［2］王遠(yuǎn)成，張忠杰，吳子丹，等.計(jì)算流體力學(xué)技術(shù)在糧食儲(chǔ)藏中的應(yīng)用［J］.中國(guó)糧油學(xué)報(bào)，2012，27（5）：86－91

Wang Y C，Zhang Z J，Wu Z D，et al.Application of computational fluid dynamics technology in grain storage systems［J］.Journal of the Chinese Cereals and Oils Association，2012，27（5）：86－91

［3］Thorpe G R，Whitaker S.Local mass and thermal equilibria in ventilated grain bulks.Part I：The development of heat and mass conservation equations［J］.Journal of Stored Products Research，1992，28：15 －27

［4］Thorpe G R，Whitaker S.Local mass and thermal equilibria in ventilated grain bulks.Part II：The developments of constraints［J］.Journal of Stored Products Research，1992b，28：29 －54

［5］林瑞泰.多孔介質(zhì)傳熱傳質(zhì)引論［M］.北京：科學(xué)出版社，1995：200－265

Lin R T.Introduction to heat and mass transfer in porous medium ［M］.Beijing，China：Science Press，1995，200 －265

［6］王遠(yuǎn)成，段海峰，張來(lái)林.就倉(cāng)通風(fēng)時(shí)糧堆內(nèi)部熱濕耦合傳遞過(guò)程的數(shù)值預(yù)測(cè)［J］.河南工業(yè)大學(xué)學(xué)報(bào)，2009，30（6）：75－79

Wang Y C，Duan H F，Zhang L L.Numerical simulation of heat and moisture transfer in stored grains at near－ambient drying ［J］.Journal of Henan University of Technology （Natural Science Edition），2009，30（6）：75 －79

［7］Wang Yuancheng，Duan Haifeng，Zhang Hao，et al.Modeling on heat and mass transfer in stored wheat during forced cooling ventilation［J］.Journal of Thermal Science，2010，19（2）：167 －172

［8］Wang Yuancheng，Zhang Congjv，Shi Shuo.Numerical simulation and validation of heat and moisture transfer in stored grains during ventilation［C］.Proceedings of the 7th International Symposium on Heating，Ventilating and Air Conditioning，Shanghai：2011：1006 －1011.

Numerical Study During the Horizontal Grain Aeration with Cooling Air

Wang Yuancheng1Gao Shuai1Qiu Huayu1Wei Lei2Zhao Huiyi2Cao Yang2
（College of Thermal Energy Engineering，Shandong Jianzhu University1，Jinan 250101）
（Academy of State Administration of Grains2，Beijing 100037）

A mathematical model of internal flow and heat and humidity coulped transport of grain bulk during the aeration of stored grains based on the principle of local heat and humidity balance and theory of heat and mass transfer in porous media is established.Based on the method of computational fluid dynamics，numerical analysis has done for internal air motion，heat transfer and water migration process of grain bulk during the horizontal grain aeration with cooling air.It’s found by study that in case of horizontal aeration for about 72 h，average temperature of cold air at the air intake of grain elevator，with relative humidity being 85%，temperature of grain bulk inside decreases from32.2℃to 23.6℃，with temperature drop range being 8.6℃.Moisture decreases from12.2%to 12.0%with dropping extent being 0.2%.Compared to ground cage vertical aeration，horizontal aeration is characterized by homogeneous distribution of internal speed of grain bulk，less temperature gradient，quick cooling rate and high cooling efficiency.

numerical simulation，horizontal aeration，heat and mass transfer，stored grains

TH432

A

1003－0174（2016）07－0103－04

國(guó)家自然科學(xué)基金（51276102），國(guó)家糧食公益專項(xiàng)（201313001，2015449 －001 －03）

2014－11－19

王遠(yuǎn)成，男，1963出生，教授，博士，復(fù)雜系統(tǒng)流動(dòng)和傳熱傳質(zhì)