張曉靜,王遠(yuǎn)成,潘鈺, 尉堯方, 高帥

(山東建筑大學(xué) 教育部可再生能源建筑利用技術(shù)實(shí)驗(yàn)室, 山東 濟(jì)南 250101)

?

靜態(tài)糧堆冷卻與干燥通風(fēng)溫濕度模擬研究

張曉靜,王遠(yuǎn)成*,潘鈺, 尉堯方, 高帥

(山東建筑大學(xué) 教育部可再生能源建筑利用技術(shù)實(shí)驗(yàn)室, 山東 濟(jì)南 250101)

小麥儲(chǔ)藏過(guò)程中，溫度和水分是安全儲(chǔ)糧的關(guān)鍵因素，通風(fēng)過(guò)程中小麥堆的溫度和水分的模擬研究對(duì)于安全儲(chǔ)糧具有指導(dǎo)意義。文章基于局部熱質(zhì)平衡原理采用Fortran程序建立了一套新的數(shù)學(xué)模型，圍繞靜態(tài)糧堆冷卻與干燥通風(fēng)過(guò)程，通過(guò)對(duì)不同通風(fēng)情況下糧堆溫度和水分變化的模擬研究，闡明了糧堆通風(fēng)過(guò)程中溫度和水分變化的一般規(guī)律。結(jié)果表明:對(duì)于小麥堆的干燥過(guò)程，小麥堆的溫度先升高達(dá)到峰值26 ℃后又逐漸趨近于進(jìn)風(fēng)溫度，隨著通風(fēng)時(shí)間的增加，小麥堆水分普遍降低；對(duì)于小麥的冷卻過(guò)程，小麥堆的溫度隨著進(jìn)風(fēng)狀態(tài)的改變而改變；通過(guò)模擬結(jié)果以及實(shí)驗(yàn)結(jié)果的對(duì)比分析，文中模型對(duì)溫度的模擬結(jié)果的誤差最大只有2 ℃，測(cè)量和預(yù)測(cè)水分的最大差值為1.19%。

Fortran程序；儲(chǔ)糧通風(fēng)；傳熱傳質(zhì)

0　引言

糧食儲(chǔ)藏過(guò)程中，水分和溫度是特別重要的兩個(gè)參數(shù)。環(huán)境溫度的季節(jié)性變化，會(huì)導(dǎo)致糧堆內(nèi)熱量傳遞和水分遷移，繼而引起糧堆局部的溫度和水分的升高，從而導(dǎo)致微生物和害蟲(chóng)的生長(zhǎng)，使得糧倉(cāng)中的糧食變得不安全[1]。糧堆的呼吸作用會(huì)消耗糧食的干物質(zhì)，同時(shí)呼吸過(guò)程會(huì)產(chǎn)生熱量和水分，如果糧堆通風(fēng)情況不好，有可能引起糧堆的“出汗”甚至“結(jié)頂”[2]。

經(jīng)過(guò)長(zhǎng)期的實(shí)踐和研究，人們發(fā)現(xiàn)當(dāng)小麥堆溫度低于 15 ℃(小麥安全儲(chǔ)藏的溫度)，小麥堆水分在 12%～12.5%(安全水分)時(shí)，可以有效地避免蟲(chóng)害的發(fā)生，抑制糧堆中生物體的生命活動(dòng)，延緩儲(chǔ)糧品質(zhì)的劣變[3]。通風(fēng)過(guò)程中糧堆內(nèi)部的流動(dòng)是一個(gè)非常復(fù)雜的過(guò)程，它與糧堆熱物性參數(shù)和糧食生物特性等多種因素有關(guān)，而且涉及到流體力學(xué)、傳熱學(xué)[4],工業(yè)通風(fēng)[5]、生物學(xué)以及多孔介質(zhì)內(nèi)部的流動(dòng)傳遞理論[6]。在機(jī)械通風(fēng)方面，我國(guó)近幾年一些技術(shù)力量較好的國(guó)儲(chǔ)庫(kù)在研究人員指導(dǎo)下進(jìn)行了研究和試驗(yàn)，提出原始水分在16.5%～18.0% 左右的糧食，利用機(jī)械通風(fēng)完全可均勻地將其干燥于安全水分之內(nèi)[7]。就倉(cāng)機(jī)械通風(fēng)在小麥儲(chǔ)藏中具有降溫效果顯著、費(fèi)用較低等特點(diǎn)，在確保儲(chǔ)糧安全方面，發(fā)揮著至關(guān)重要的作用。

為了掌握通風(fēng)過(guò)程中儲(chǔ)糧生態(tài)系統(tǒng)變化規(guī)律，儲(chǔ)藏工作者進(jìn)行了大量的實(shí)地測(cè)試工作，以研究特定的天氣條件下谷物通風(fēng)系統(tǒng)的性能和空氣流速以及糧堆內(nèi)部溫濕度變化規(guī)律。但是通風(fēng)實(shí)驗(yàn)需要投入較大的人力物力，成本較高，實(shí)驗(yàn)結(jié)果還不具有可重復(fù)性。為了減少工作量，Sharma開(kāi)發(fā)了仿真模型來(lái)描述在可變參數(shù)通風(fēng)下熱與質(zhì)量傳遞現(xiàn)象[8]。由于模型的限制，這些模型不能模擬谷物在整個(gè)通風(fēng)段的溫度?；跓崃亢唾|(zhì)量守衡定律，Jia等模擬了就倉(cāng)通風(fēng)過(guò)程中糧倉(cāng)內(nèi)儲(chǔ)藏小麥的溫度變化，但是沒(méi)有模擬小麥堆水分的變化[9]。王平等對(duì)平房倉(cāng)橫向通風(fēng)降溫進(jìn)行了模擬研究[10]，呂宗旺等以CFD技術(shù)對(duì)糧倉(cāng)進(jìn)行了降溫模擬[11]，王遠(yuǎn)成等對(duì)大型房式倉(cāng)地上籠糧堆溫度和水分變化規(guī)律[12]和圓筒倉(cāng)內(nèi)自然對(duì)流對(duì)糧堆熱濕傳遞的影響[13]以及倉(cāng)儲(chǔ)糧堆內(nèi)熱濕耦合傳遞[14]進(jìn)行了數(shù)值模擬研究。

以上發(fā)展的模型雖然能比較準(zhǔn)確地模擬通風(fēng)過(guò)程中溫度和水分的變化，但是一般都存在計(jì)算周期長(zhǎng)，占用計(jì)算機(jī)內(nèi)存大的缺陷。而且在設(shè)定初始條件的時(shí)候，大都需要知道所模擬糧種的導(dǎo)熱系數(shù)等參數(shù)，而糧種的這些熱物性參數(shù)受糧食水分以及糧溫的影響比較大，目前國(guó)內(nèi)對(duì)這方面的統(tǒng)計(jì)數(shù)據(jù)尚不完善。文章通過(guò)集總參數(shù)的方法結(jié)合Fortran程序設(shè)計(jì)了一套新的數(shù)學(xué)模型。該模型不針對(duì)某一特種糧食，而是普遍適用于各種儲(chǔ)糧品種，只需要在原程序的基礎(chǔ)上修改糧食的比熱容即可。文章對(duì)儲(chǔ)糧過(guò)程中小麥堆的溫度和水分變化進(jìn)行模擬，獲得在不同通風(fēng)條件下小麥堆的溫度和水分，通過(guò)與試驗(yàn)結(jié)果的比較也證明這個(gè)模型的準(zhǔn)確性與可行性。

1　模型的建立

1.1物理模型的建立

文章將糧堆劃分為相等厚度的谷物薄層，這些糧層間的空氣流動(dòng)的方向如圖1所示，其中，ΔH為空間增量，Ti為進(jìn)風(fēng)溫度，Wi為進(jìn)風(fēng)空氣濕度比率，To為出風(fēng)溫度，Tw為小麥堆溫度，Wo為出風(fēng)空氣濕度比率。建立仿真模型的目的是預(yù)測(cè)在規(guī)定的時(shí)間內(nèi)每一層糧層最終的溫度和水分含量。從上一層排出的空氣狀態(tài)被用作輸入下一層糧層的條件，此過(guò)程被重復(fù)，直到一段時(shí)間的增量上整個(gè)谷物柱上的模擬完成。然后該模型以這段時(shí)間增量末期的谷物狀態(tài)為初始條件，開(kāi)始模擬在進(jìn)風(fēng)狀態(tài)下下一個(gè)時(shí)間增量的第一層的小麥狀態(tài)，然后使用和第一次時(shí)間增量相同的程序，模擬了在第二個(gè)時(shí)間增量的整個(gè)小麥堆。不斷循環(huán)直到完成小麥的整個(gè)通風(fēng)過(guò)程。

圖1　糧層間空氣流動(dòng)方向示意圖

1.2數(shù)學(xué)模型的建立

1.2.1焓濕關(guān)系

Brooke提出焓濕關(guān)系式來(lái)計(jì)算空氣的屬性[15]。

當(dāng)干球溫度Tdb大于或等于0 ℃時(shí)，飽和蒸氣壓Ps1(g/cm2)由式(1)計(jì)算為

Ps1=exp[58.8858-12301.69/(491.69+1.8Tdb)-5.16923ln(491.69+1.8Tdb)]

(1)

當(dāng)干球溫度Tdb低于0 ℃時(shí)，飽和蒸氣壓Ps2(g/cm2)由式(2)計(jì)算為

Ps2=exp[27.6452-11286.1489/(491.69+1.8Tdb)-0.46057ln(491.69+1.8Tdb)]

(2)

在給定濕度下的情況下，氣體壓力為Pv=(RH)PS，其中，RH為進(jìn)風(fēng)的相對(duì)濕度。

空氣的濕度比率W=0.6219Pv/(Patm-Pv)，Patm是標(biāo)準(zhǔn)大氣壓，g/cm2。

空氣的比體積vsa=0.007573(491.69+1.8Tdb)(1+1.6055W)。

空氣的質(zhì)量流速ma=Qa/vsa；Qa是空氣的體積流量，m3/h?？諝獾谋葻醕pa=0.2405+0.44W，自由水分的蒸發(fā)潛熱ls=597.768-0.56983T， T是糧堆的溫度，℃

空氣焓h由ASHRAE基本手冊(cè)中的式(3)計(jì)算為[16]h=4.267+0.240Tdb+W(597.386+0.444Tdb)

(3)

1.2.2小麥基本屬性

Othmer等[17]推導(dǎo)出了小麥達(dá)到平衡時(shí)空氣的相對(duì)濕度RHe(小數(shù)) 由式(4)計(jì)算為RHe=exp[(l/ls)ln(2.04816Ps)+c]/2.04816Ps

(4)

式中：l是小麥中水分的蒸發(fā)潛熱，kcal/kg水;l/ls=1+23exp(-0.470 M);c=-3.34×104M-4.0;M是小麥的水分，干基百分比l/ls為1+23exp(-0.40M)。

Viravanichai[18]研究了熱量的以下關(guān)系，硬紅春小麥在溫度范圍為-33.5 至21.8 ℃的小麥的比熱cpw(kcal/kg·℃)值由式(5)～(9)表示為

cpw1=0.250+0.00743M(-33.5 ～-24.1 ℃)

(5)

cpw2=0.279+0.00580M(-21.4 ～-10.8 ℃)

(6)

cpw3=0.297+0.00708M(-10.8 ～3.6 ℃)

(7)

cpw4=0.245+0.01058M(0.6 ～8.9 ℃)

(8)

cpw5=0.273+0.00933M(8.9 ～21.8 ℃)

(9)

1.2.3熱濕傳遞公式

作為中國(guó)改革開(kāi)放的前沿，云南具有優(yōu)越的區(qū)位優(yōu)勢(shì)，是“一帶一路”、孟中印緬經(jīng)濟(jì)走廊、中國(guó)—中南半島經(jīng)濟(jì)走廊的重要節(jié)點(diǎn)，是中國(guó)唯一可以同時(shí)從陸上溝通東南亞、南亞的省，并可以通過(guò)中東聯(lián)接歐洲、非洲，開(kāi)拓西向貿(mào)易通道的省份。獨(dú)特的區(qū)位優(yōu)勢(shì)，決定了云南將在“一帶一路”建設(shè)中發(fā)揮不可替代的作用。

2.3.1熱平衡公式

(1) 根據(jù)局部熱平衡，可以得到式(10)計(jì)算為

maΔtcpa(Ti-To)+mgcpw(Tw-To)+(hi-ho)maΔt l/ls=0

(10)

式中：ma為空氣的質(zhì)量流速，kg干空氣/h；mg為糧層中的小麥質(zhì)量，kg；hi、ho分別為進(jìn)風(fēng)的焓，kcal/kg干空氣。

對(duì)于干燥過(guò)程，可變形為

空氣損失的熱=糧堆獲得的熱+蒸發(fā)潛熱

(2) 水分平衡公式

根據(jù)質(zhì)量守恒原理，可以得到式(11)計(jì)算為

(Wi-Wo)100maΔt+(Mo-Mi)md=0

(11)

式中：md是糧層中的固體質(zhì)量，kg；Mo是在時(shí)間增量Δt內(nèi)模擬的小麥堆水分，干基百分比；Mi是小麥堆的初始水分，干基百分比；Δt是時(shí)間增量，h。

對(duì)于干燥過(guò)程，可變形為

空氣獲得的水分=糧堆失去的水分

文中模型以空氣逐時(shí)的溫度和相對(duì)濕度、糧堆的深度、糧層的厚度、空氣流速、谷物的初始水分和初始溫度為初始輸入數(shù)據(jù)，各種基礎(chǔ)屬性由上述公式計(jì)算得到。

比較Pvi和Pg、Tw和Ti的關(guān)系，存在下面四種可能：

加熱干燥小麥 Ti>Tw,Pvi

冷卻加濕小麥 Ti>Tw,Pvi>Pg

冷卻通風(fēng)小麥 Ti

加熱加濕小麥 Ti>Tw,Pvi>Pg

在以上四個(gè)過(guò)程中，都達(dá)到了熱質(zhì)平衡。為了發(fā)展數(shù)學(xué)模型，文中模型假設(shè)在空氣溫度下首先發(fā)生水分傳遞，然后發(fā)生熱量傳遞。假設(shè)達(dá)到平衡時(shí)，小麥堆的出風(fēng)溫度和小麥堆的溫度是相同的。同樣的，小麥堆的出風(fēng)空氣濕度比率和小麥堆的比濕達(dá)到平衡。而實(shí)際在小麥堆通風(fēng)過(guò)程中，溫度的變化對(duì)小麥堆的影響比水分變化快。因此，當(dāng)已經(jīng)達(dá)到溫度平衡時(shí)，水分平衡還沒(méi)有達(dá)到。鄰近通風(fēng)結(jié)束階段時(shí)，小麥堆的水分和溫度與空氣都達(dá)到平衡。

當(dāng)水分傳遞未達(dá)到平衡時(shí)，出風(fēng)空氣濕度比率可以由式(12)表示為

(12)

R(R以百分比表示)定義為：糧層內(nèi)實(shí)際傳遞的水分與在平衡狀態(tài)下傳遞水分之比。

干燥和調(diào)質(zhì)實(shí)驗(yàn)是在溫度為6.8～25.0 ℃的范圍內(nèi)進(jìn)行的，通過(guò)反復(fù)的實(shí)驗(yàn)，R值應(yīng)為80。在-19.0～0 ℃的溫度范圍調(diào)質(zhì)測(cè)試使用的R為經(jīng)驗(yàn)數(shù)值。

R值在小麥堆糧層厚度和空氣流速一定的情況下取決于空氣和小麥的狀態(tài)。

1.3數(shù)值方法和初始條件

文章模擬了Sharma[8]在實(shí)驗(yàn)室做的實(shí)驗(yàn)工況，得出了120cm深的糧堆干燥通風(fēng)和冷卻通風(fēng)時(shí)的溫度和濕度,實(shí)驗(yàn)條件見(jiàn)表1。

表1　用于通風(fēng)模型的變量

2　通風(fēng)實(shí)驗(yàn)

通風(fēng)實(shí)驗(yàn)的原理圖如圖2所示。實(shí)驗(yàn)設(shè)備可以被分成三部分：儲(chǔ)存糧食的垂直通風(fēng)柱，調(diào)節(jié)空氣的空氣室，調(diào)節(jié)穿過(guò)糧食的空氣的空氣分配室。下面是這個(gè)設(shè)備的簡(jiǎn)介，更多的細(xì)節(jié)Sharma已經(jīng)給出[8]。

圖2　通風(fēng)實(shí)驗(yàn)原理圖

通風(fēng)測(cè)試使用的小麥放在由丙烯腈-丁二烯-苯乙烯(ABS)制造的直徑為10.16cm的觀察管內(nèi)。觀察管的高度為24.13cm，糧層的厚度為15.2cm。沿著觀察管設(shè)置了6個(gè)直徑為1.27cm的孔洞，通過(guò)孔洞里抽樣檢測(cè)小麥，沿著觀察管軸線(xiàn)分布有熱電偶。Brooks玻璃管流量計(jì)(型號(hào)1110)安裝在通風(fēng)柱和排風(fēng)支管之間的流動(dòng)管線(xiàn)上，用于測(cè)量空氣流速。系統(tǒng)采用機(jī)械通風(fēng)，處理過(guò)的空氣在風(fēng)機(jī)的作用下經(jīng)過(guò)空氣室穿過(guò)全孔板對(duì)糧食進(jìn)行通風(fēng)。

3　模擬結(jié)果分析

將文章利用程序計(jì)算出來(lái)的結(jié)果(程序值)與實(shí)驗(yàn)值和Sharma的模擬結(jié)果(模擬值)[8]作比較，得到以下圖形。

3.1工況1下模擬結(jié)果分析

工況1是在進(jìn)風(fēng)參數(shù)恒定的狀態(tài)下進(jìn)行模擬的，模擬結(jié)果如圖3、4所示。

由模擬結(jié)果可以發(fā)現(xiàn)，隨著通風(fēng)時(shí)間的增加，小麥堆溫度逐漸升高，達(dá)到峰值后逐漸趨于進(jìn)風(fēng)溫度(如圖3(a)所示)，在解吸濕過(guò)程中小麥釋放潛熱，使得底層小麥堆的溫度在通風(fēng)的前幾個(gè)小時(shí)里高于測(cè)量溫度而在通風(fēng)后期預(yù)測(cè)溫度又低于測(cè)量溫度，這個(gè)差值最大可達(dá)到2 ℃。同時(shí)由于局部熱平衡的作用，隨著糧層高度的增加，溫度變化越來(lái)越緩和(如圖3(b)所示)。通風(fēng)15h時(shí)，糧層內(nèi)水分出現(xiàn)一個(gè)峰值(如圖4(a)所示)，隨著通風(fēng)時(shí)間的增加，峰值小幅降低且峰值后水分含量普遍降低(如圖4(b)所示)，說(shuō)明當(dāng)溫度達(dá)到平衡時(shí)，水分尚未達(dá)到平衡。

3.2工況2下模擬結(jié)果分析

工況2是在進(jìn)風(fēng)參數(shù)隨時(shí)間不斷變化的狀態(tài)下進(jìn)行模擬的，模擬結(jié)果如圖5、6、7所示。

圖3　小麥堆各層溫度隨時(shí)間變化圖(a)小麥堆第1層；(b)小麥堆第5層

圖4　不同時(shí)刻下小麥堆各層水分變化圖(a)15 h ；(b)25 h

圖5　不同R值對(duì)小麥堆水分含量的影響圖(80 h)

低R值表征水分傳遞較少，即從空氣和小麥中帶走少量顯熱。而高R值表征水分傳遞較多，即從空氣和小麥中帶走大量顯熱。由模擬結(jié)果可知最終模擬R值等于80最切合實(shí)際(如圖5所示)。通風(fēng)前期，第一糧層的預(yù)測(cè)溫度比測(cè)量溫度要高。這個(gè)差值最大到2 ℃(如圖6(a)所示)。這是因?yàn)楫?dāng)空氣流過(guò)糧層時(shí)，它的濕度趨向于平衡濕度。因此，在較上面的糧層里，空氣的濕度與平衡時(shí)的濕度差值很小。這就會(huì)使得達(dá)到平衡時(shí)水分傳遞的實(shí)際速率比最大傳遞速率的80%大，即上層糧層的R值應(yīng)該比80大，這也是通風(fēng)前期預(yù)測(cè)水分低于測(cè)量水分(如圖7所示)的原因。模型預(yù)測(cè)的小麥堆溫度取決于在吸附過(guò)程中空氣潛熱的增加。在通風(fēng)的前幾個(gè)小時(shí)內(nèi)，最高可傳送水分的百分比是最低的。因此，R值應(yīng)該比80小。這可能是在通風(fēng)初期的幾個(gè)小時(shí)里預(yù)測(cè)溫度比測(cè)量溫度低(如圖6(b)所示)的原因。

圖6　小麥堆各層溫度隨時(shí)間變化圖(a)小麥堆第1層；(b)小麥堆第7層

圖7　80 h后小麥堆各層水分變化圖

在這些實(shí)驗(yàn)中測(cè)量和預(yù)測(cè)水分的最大差值為1.19%。產(chǎn)生差值的原因有可能和以下因素有關(guān)：(1)由于在每個(gè)實(shí)驗(yàn)中小麥和空氣的狀態(tài)不同，取相同的R值是不準(zhǔn)確的；(2)R值的變化預(yù)示著水分含量的變化；(3)糧層的邊緣效應(yīng)，進(jìn)風(fēng)相對(duì)濕度和小麥堆水分含量的測(cè)量誤差。

4　結(jié)論

通過(guò)上述研究可知：

(1) 儲(chǔ)糧通風(fēng)過(guò)程中溫度與濕度是相互耦合的。

(2) 對(duì)于小麥堆的干燥過(guò)程，小麥堆的溫度先升高達(dá)到峰值后(26 ℃)又逐漸降低并趨近于進(jìn)風(fēng)溫度(22 ℃)，隨著通風(fēng)時(shí)間的增加，小麥堆水分峰值降低且峰值后水分含量普遍降低；對(duì)于小麥的冷卻過(guò)程，小麥堆的溫度隨著進(jìn)風(fēng)狀態(tài)的改變而改變，水分含量與實(shí)驗(yàn)結(jié)果比較吻和。

(3) 通過(guò)與前人模擬結(jié)果以及實(shí)驗(yàn)結(jié)果的比較發(fā)現(xiàn)，文中模型具有更高的準(zhǔn)確性，溫度模擬結(jié)果的誤差最大只有2 ℃，測(cè)量和預(yù)測(cè)水分的最大差值為1.19%。文中模型在局部平衡的基礎(chǔ)上首次提出了R值，使水分傳遞更接近于實(shí)際情況。

(4) 通過(guò)與測(cè)量結(jié)果的比較，驗(yàn)證了該模型的準(zhǔn)確性與可行性。該模型的突出優(yōu)點(diǎn)就是占用內(nèi)存小，運(yùn)行速度快，操作簡(jiǎn)單，普適性強(qiáng)，與之前的數(shù)學(xué)模型相比，有更好的應(yīng)用價(jià)值。

[1]楊廣靖,任云虹,賈金，等. 我國(guó)糧食儲(chǔ)藏的現(xiàn)狀及發(fā)展趨勢(shì)[J].糧食加工，2012， 37(1):60-63.

[2]王遠(yuǎn)成, 武傳欣, 段海峰. 糧食儲(chǔ)藏過(guò)程中糧堆內(nèi)部熱濕耦合傳遞過(guò)程對(duì)糧食品質(zhì)變化的影響研究[C].中國(guó)糧油儲(chǔ)藏學(xué)會(huì)第六屆年會(huì)報(bào)告，德宏：中國(guó)糧油儲(chǔ)藏學(xué)會(huì),2009.

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

[4]陶文銓.數(shù)值傳熱學(xué)[M].西安:交通大學(xué)出版社,1998.

[5]孫一堅(jiān)，沈恒根. 工業(yè)通風(fēng)[M].北京：中國(guó)建筑工業(yè)出版社,2010.

[6]林瑞泰.多孔介質(zhì)傳熱傳質(zhì)引論[M].北京:科學(xué)出版社,1995.

[7]湯文良，狄育慧.關(guān)于我國(guó)糧倉(cāng)的通風(fēng)技術(shù)及低溫儲(chǔ)藏技術(shù)[J].潔凈與空調(diào)技術(shù),2013,3(1):44-46.

[8]SharmaS.C. .Simulationofheatandmasstransferduringaerationofwheatandrapeseedbu1ks[J].UnpublishedM.Sc.Thesis,UniversityofManitoba,Winnipeg, 1973(31):104-105.

[9]JiaC.C.，SunD.W.，CaoC.W..Computersimulationoftemperaturechangesinawheatstoragebin[J]．JournalofStoredProductsResearch. 2001, 37( 1) : 165 -177.

[10]王平，周焰，曹陽(yáng)，等.平房倉(cāng)橫向通風(fēng)降溫技術(shù)研究[J]. 糧油倉(cāng)儲(chǔ)科技通訊2011( 2)：20.

[11]呂宗旺，馮黎明，孫福艷.糧倉(cāng)通風(fēng)CFD技術(shù)研究[J].河南工業(yè)大學(xué)學(xué)報(bào)，2013,34(6):105-109.

[12]王遠(yuǎn)成,張忠杰,吳子丹，等.計(jì)算流體力學(xué)技術(shù)在小麥儲(chǔ)藏中的應(yīng)用[J].中國(guó)糧油學(xué)報(bào)，2008,27(5):86-90.

[13]王遠(yuǎn)成, 亓偉, 張中濤,等. 圓筒倉(cāng)內(nèi)自然對(duì)流對(duì)糧堆熱濕傳遞的影響研究[J].水動(dòng)力學(xué)研究與進(jìn)展A,2015,25(4):487-496.

[14]王遠(yuǎn)成,白忠權(quán),張中濤，等.倉(cāng)儲(chǔ)糧堆內(nèi)熱濕耦合傳遞的數(shù)值模擬[J].中國(guó)糧油學(xué)報(bào),2015,3(11):64-66.

[15]BrookerD.B..Mathematicalmodelofthepychometrics[J].Chart.Trans.ASAE,1967,11( 10): 558-560.

[16]Ashrae.HandbookofFundamentals[M].NewYork:AmericanSocietyofHeating,RefrigeratingandAir-conditioningEngineer,1972.

[17]OthmerD.F,HuangH. .Correlatingvapourpressureandlatentheatdata[J],Ind.Eng.Chem. 1940,32: 841-846.

[18]ViravanichaiS..EffectofMoisturecontentandtemperatureonspecificheatofwheat[D].Winnipeg:Vniverstityofmanitoba, 1971.

(學(xué)科責(zé)編：吳芹)

Simulation of temperature and moisture in static grainstorageduringcoolanddryventilation

Zhang Xiaojing,Wang Yuancheng*, Pan Yu, et al.

(Key Laboratory of Renewable Energy Utilization Technologies in Building of the National Education Ministry, Shandong Jianzhu University, Jinan 250101, China)

Thetemperatureandmoisturecontentiscriticaltoachievingsafewheatstorage.Theresearchershavedesignedavarietyofmathematicalmodelstopredictthetemperatureandmoisturecontentofwheatduringventilation,whichgenerallytakesalotofmoneyandtime.TheproblemcanbesolvedeffectivelybyanewmathematicalmodelbasedonlocalheatandmassbalanceprincipleandFortranprogram.Itmakesthesimulationresearchofgraintemperatureandmoistureindifferentventilationconditionsandexploresthegenerallawofgraintemperatureandmoisturevariation.Theresultsshow:thetemperatureofwheatrisespeakfirstandthenapproachesthetemperatureofincomingairgraduallyduringthedryingventilation,andthemoistureofwheatdecreasesgenerallywiththeincreaseofventilationtime;thetemperatureofthewheatchangeswiththeconditionofincomingairduringcoolingventilation.Themodelhashigheraccuracycomparedwiththeprevioussimulationresultsandexperimentalresults,anditexists2 ℃errorintemperatureand1.19%errorinmoistercontentonly.Themodelisuniversalandpracticalfordifferentkindsofgrainanditisverysignificantfortheguidanceofwheatstorage.

Fortranprogramcalculation;ventilationofwheatstorage;heatandmasstransfer

2015-11-06

國(guó)家自然基金項(xiàng)目(51276102)；國(guó)家糧食公益專(zhuān)項(xiàng)項(xiàng)目(201313001)；國(guó)家糧食公益專(zhuān)項(xiàng)項(xiàng)目(2015449-001-03)

張曉靜(1990-)，女，在讀研究生，主要從事多孔介質(zhì)熱質(zhì)傳遞等方面的研究.E-mail:930065394@qq.com

*:王遠(yuǎn)成(1963-)，男,教授,博士，主要從事復(fù)雜系統(tǒng)的傳熱傳質(zhì)等方面的研究.E-mail:wycjn1@163.com

1673-7644(2016)01-0047-06

TS205

A