劉兆宇　羅會(huì)龍　陳頤　鄒聰明

摘要：以標(biāo)準(zhǔn)密集烤房為對(duì)象，建立計(jì)算流體力學(xué)（CFD）模型，運(yùn)用大型模擬軟件（Fluent），對(duì)烘烤干筋期密集烤房?jī)?nèi)部的溫度場(chǎng)、濕度場(chǎng)及氣流組織進(jìn)行了系統(tǒng)的數(shù)值模擬?；跀?shù)值模擬結(jié)果，分析了密集烘烤典型送風(fēng)參數(shù)下溫度場(chǎng)、濕度場(chǎng)、氣流組織分布的規(guī)律及進(jìn)風(fēng)口風(fēng)速與煙葉間隙風(fēng)速的關(guān)系。綜合考慮風(fēng)機(jī)能耗、葉間風(fēng)速，進(jìn)風(fēng)口風(fēng)速宜選取2.2～3m/s。

關(guān)鍵詞：煙葉;干筋期;數(shù)值模擬;溫度場(chǎng);風(fēng)速場(chǎng);濕度場(chǎng);密集烤房

中圖分類號(hào)：TS43文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1002-1302（2021）01-0160-06

作者簡(jiǎn)介：劉兆宇（1996—），男，湖南邵東人，碩士研究生，主要從事熱泵烘烤及建筑節(jié)能相關(guān)研究。E-mail：zhaoyuliu123@163.com。

通信作者：羅會(huì)龍，博士，教授，主要從事熱泵烘烤及建筑節(jié)能相關(guān)研究。E-mail：huilongkm@126.com。

進(jìn)入21世紀(jì)后，隨著我國(guó)烤煙規(guī)模化生產(chǎn)的發(fā)展，密集烤房在我國(guó)已成為煙葉烘烤設(shè)備主要的發(fā)展方向[1-2]?？緹熀婵靖山钇谑菬熑~主脈水分變化的主要階段，同時(shí)也是煙葉香氣物質(zhì)含量和組分變化的關(guān)鍵階段[3-6]。在烘烤干筋期，密集烤房?jī)?nèi)合理的溫度、濕度可以有效提高煙葉烤后的外觀質(zhì)量，進(jìn)而提高煙葉烘烤品質(zhì)[7]，同時(shí)干筋期不同的葉間風(fēng)速對(duì)烤后煙葉的香氣品質(zhì)有直接影響[8-9]。因此，研究烤煙干筋期密集烤房?jī)?nèi)溫度、濕度分布及氣流組織分布有著重要意義。

近年來(lái)，隨著計(jì)算流體力學(xué)（CFD）的發(fā)展，其商業(yè)軟件Fluent已廣泛應(yīng)用于流體和熱量傳遞的數(shù)值模擬研究中[10-11]。高維潔等利用Fluent軟件對(duì)列管式換熱器性能校核的方法進(jìn)行了研究[12]。趙海霞等利用Fluent軟件對(duì)霧化器霧滴粒徑進(jìn)行了數(shù)值模擬[13]。方文康等利用Fluent軟件研究了多溫區(qū)保溫包裝箱不同條件對(duì)保溫箱內(nèi)部溫度場(chǎng)及保溫效果的影響[14]。隨著計(jì)算流體力學(xué)在農(nóng)業(yè)領(lǐng)域的發(fā)展，F(xiàn)luent軟件在農(nóng)業(yè)工程中的應(yīng)用也日益增多。徐陽(yáng)等利用Fluent軟件對(duì)水草收割臺(tái)進(jìn)行了結(jié)構(gòu)優(yōu)化[15]。王娟等利用Fluent對(duì)飼草料揉碎機(jī)的三維流場(chǎng)進(jìn)行了數(shù)值模擬[16]。段二亞利用Fluent對(duì)糧食干燥過(guò)程進(jìn)行了數(shù)值模擬研究[17]。

對(duì)于煙葉密集烘烤，通過(guò)試驗(yàn)來(lái)研究烤房?jī)?nèi)部熱濕環(huán)境不僅浪費(fèi)人力、物力，還大大地增加了生產(chǎn)成本，因此運(yùn)用Fluent軟件進(jìn)行數(shù)值模擬研究有著廣闊的應(yīng)用前景，可以較為精確地分析密集烤房?jī)?nèi)部溫度場(chǎng)、濕度場(chǎng)及氣流組織分布規(guī)律，為最佳的烘烤工藝提供理論依據(jù)。在此背景下，本研究運(yùn)用大型模擬軟件Fluent，對(duì)烘烤干筋期密集烤房?jī)?nèi)部的溫度場(chǎng)、濕度場(chǎng)及氣流組織進(jìn)行了系統(tǒng)的數(shù)值模擬分析。

1材料與方法

1.1基本假設(shè)

（1）烤房?jī)?nèi)的空氣為不可壓縮的理想氣體;（2）烤房的結(jié)構(gòu)密封性良好，沒(méi)有漏風(fēng)現(xiàn)象產(chǎn)生;（3）將烤房?jī)?nèi)的煙葉視為連續(xù)的多孔介質(zhì)，流體在煙葉區(qū)域定義為層流運(yùn)動(dòng)，忽略湍流運(yùn)動(dòng)的影響;（4）忽略煙葉骨架間的導(dǎo)熱和葉間的輻射傳熱作用，換熱過(guò)程中僅涉及對(duì)流換熱。

1.2裝煙室物理模型

根據(jù)國(guó)家煙草總局Q/GDYY019—2011《密集式烤房建設(shè)技術(shù)規(guī)范》建立模型，其尺寸為8200mm×2700mm×3500mm，裝煙室可裝3層煙葉，進(jìn)風(fēng)口尺寸為2700mm×400mm，出風(fēng)口尺寸為1400mm×400mm，裝煙室結(jié)構(gòu)圖如圖1所示。

1.3數(shù)學(xué)模型

任何流體的流動(dòng)都遵守質(zhì)量守恒定律與動(dòng)量守恒定律，由于裝煙室內(nèi)的熱空氣流不斷與煙葉中的水分發(fā)生熱質(zhì)交換，密集烤房?jī)?nèi)的熱濕傳遞及氣流組織的數(shù)學(xué)模型如下[18-20]：

1.4模型建立及參數(shù)設(shè)置

1.4.1網(wǎng)格劃分網(wǎng)格劃分見圖2。根據(jù)裝煙室原型，以出風(fēng)口底邊中點(diǎn)為原點(diǎn)建立坐標(biāo)軸，采用3D建模軟件SolidWorks建立烤房裝煙室模型，使用Fluent前處理軟件ICEM進(jìn)行網(wǎng)格劃分，為了加快運(yùn)算速度，使殘差更容易收斂，選用結(jié)構(gòu)化網(wǎng)格，網(wǎng)格數(shù)為162360，網(wǎng)格質(zhì)量為1，滿足模擬精度要求。

1.4.2邊界條件進(jìn)風(fēng)口設(shè)置為“velocity-inlet”，水蒸氣組分為0;出風(fēng)口設(shè)置為“pressure-outlet”，設(shè)置工作環(huán)境大氣壓為1個(gè)標(biāo)準(zhǔn)大氣壓，出風(fēng)口壓力為0;裝煙室壁面設(shè)置為“wall”，壁面條件設(shè)置為絕熱。

1.4.3Fluent參數(shù)的設(shè)置將網(wǎng)格模型導(dǎo)入Fluent中，F(xiàn)luent會(huì)對(duì)網(wǎng)格文件進(jìn)行初始化檢查，最小體積為1.228069×10-4m3，沒(méi)有出現(xiàn)負(fù)體積。根據(jù)烤房裝煙室模型，湍流模型選用K-epsilon（2eqn）方程，離散控制方程選用SIMPLE算法，設(shè)置絕對(duì)速度來(lái)計(jì)算空氣的流速。

2結(jié)果與討論

2.1裝煙室氣流場(chǎng)分布及分析

在目前典型的送風(fēng)條件下（送風(fēng)溫度T進(jìn)=340K，進(jìn)風(fēng)口風(fēng)速v進(jìn)=3m/s，相對(duì)濕度=0），截取y=0m截面的速度分布云圖。由圖3-a可知，在裝煙室頂部和底部有明顯的速度梯度變化，熱氣流以3m/s的初速度進(jìn)入裝煙室后，迅速充滿裝煙室頂部區(qū)域，呈現(xiàn)出由裝煙室前部到后部遞減的錐形速度梯度，熱氣流到達(dá)裝煙室后部時(shí)，速度減至0.85m/s左右，而在裝煙室底部呈現(xiàn)出由裝煙室后部到前部遞增的錐形速度梯度，當(dāng)熱氣流到達(dá)出風(fēng)口附近時(shí)，氣流速度最大，為7.08m/s左右。在煙葉堆積區(qū)域，氣流速度最小，同時(shí)葉間風(fēng)速分布均勻，氣流速度為0.28m/s，但在上、中層煙葉和中、下層煙葉層間水平方向上出現(xiàn)了明顯的速度變化，上、中層間速度變化在0.28～0.56m/s，中、下層間速度變化為0.28～0.85m/s。為了更加直觀地觀察煙葉堆積區(qū)域的氣流組織分布，設(shè)置如圖3-b所示的速度矢量圖，矢量箭頭的大小表示氣流組織的強(qiáng)弱。在煙葉堆積區(qū)域，氣流組織整體分布較為均勻，氣流朝著出風(fēng)口方向運(yùn)動(dòng)，但裝煙室前部的氣流組織強(qiáng)度要明顯大于裝煙室后部，并且靠近大門底部的區(qū)域氣流組織十分微弱。

本研究分析認(rèn)為熱空氣進(jìn)入裝煙室后，由于上層煙葉的阻礙作用，氣流在裝煙室頂部水平流動(dòng)，隨著熱氣流動(dòng)能損失，流速逐漸減小，所以頂部呈現(xiàn)遞減的錐形速度梯度。熱空氣進(jìn)入煙葉堆積區(qū)域后進(jìn)行層流運(yùn)動(dòng)，同時(shí)受到煙葉的黏性阻力和慣性阻力作用，所以煙葉堆積區(qū)域風(fēng)速較均勻。在進(jìn)出口的壓差作用下，出風(fēng)口處出現(xiàn)遞增的錐形速度梯度，并且速度變化較大。

2.2裝煙室溫度場(chǎng)分布及分析

在目前典型的送風(fēng)條件下（送風(fēng)溫度T進(jìn)=340K，進(jìn)風(fēng)口風(fēng)速v進(jìn)=3m/s，相對(duì)濕度=0），截取Y=0m截面的溫度分布云圖，由圖4可知，裝煙室頂部和上層煙葉區(qū)域溫度分布均勻，溫度為339.96～340.00K，中層煙葉區(qū)域至下層煙葉區(qū)域溫度梯度變化較大，且溫度向裝煙室底部遞減，裝煙室底部溫度較低，最低溫度出現(xiàn)在裝煙室大門底部區(qū)域，溫度為339.30～339.40K。

為了更直觀地看出水平方向不同位置上的溫度分布情況，本研究截取z=0.8m、z=1.6m、z=2.4m3[CM（15]個(gè)截面上的溫度分布云圖。由圖5-a可知，在上層煙葉區(qū)域，水平方向上整體溫度分布較為均勻，但在前部出現(xiàn)了輕微的溫度分層，該部分溫度為339.93K左右，水平方向上最大溫差為0.03K。由圖5-b可知，在中層煙葉區(qū)域，前部小范圍區(qū)域出現(xiàn)了溫度分層現(xiàn)象，在中后部位置更加明顯，水平方向上最大溫差為0.1K左右。由圖5-c可知，在下層煙葉區(qū)域，水平方向上的溫度分布不均勻，且出現(xiàn)明顯的分層現(xiàn)象，且溫度由前部向后部遞減，最高溫度出現(xiàn)在最前端中間位置，溫度為339.80K左右，最低溫度出現(xiàn)在最前端左右2個(gè)角落處，最低溫度為339.30K左右，水平方向上最大溫差為0.50K左右。

為了研究在不同位置垂直方向上的溫度分布，截取x=2.0m、x=4.0m、x=6.0m這3個(gè)截面的溫度分布云圖。由圖6可知，在裝煙室內(nèi)不同的位置，上層煙葉區(qū)域的垂直溫度分布較為均勻，基本沒(méi)有溫度變化，溫度為339.96K左右。但在中層和下層煙葉區(qū)域，隨著位置的改變，在垂直方向上的溫度梯度有較為明顯的差異。由圖6-a可知，裝煙室前部（靠近出風(fēng)口一側(cè)）垂直方向上的溫度梯度較均勻，溫度由上到下均勻遞減，底部溫度為339.40K左右，最大溫差為0.56K左右。由圖6-b可知，裝煙室中部垂直方向上的溫度梯度在中層煙葉區(qū)域較為均勻，但在中、下層煙葉交界區(qū)域溫度梯度與裝煙室前部相比變化較大，并且底部低溫區(qū)域占比與裝煙室前部相比有所增加，底部最低溫度為339.38K左右，垂直方向上最大溫差為0.58K左右。由圖6-c可知，裝煙室后部垂直方向上的溫度梯度變化最大，并且底部的低溫區(qū)域占比較裝煙室前部和中部出現(xiàn)明顯增大，約占垂直平面的1/5左右，最低溫度為339.35K左右，垂直方向上最大溫差為0.61K左右。

本研究分析認(rèn)為，烘烤過(guò)程中在循環(huán)風(fēng)機(jī)的作用下，熱空氣流過(guò)煙葉區(qū)域，將熱量傳遞給煙葉，并帶走煙葉干燥過(guò)程中蒸發(fā)出的水分，使得煙葉區(qū)域的熱濕空氣被新的熱濕空氣替代;同時(shí)由于氣流組織分布不均勻，導(dǎo)致裝煙室內(nèi)水平方向和垂直方向出現(xiàn)溫度梯度，裝煙室前部的氣流組織明顯大于裝煙室后部;因此裝煙室前部區(qū)域內(nèi)的熱濕空氣能夠迅速被氣流組織帶走，使得該區(qū)域溫度分布更加均勻，由于裝煙室后部氣流組織十分微弱，滯留的熱空氣不能及時(shí)被新的熱空氣所替代，溫度逐漸降低，因此出現(xiàn)較大范圍的低溫區(qū)域。

2.3裝煙室濕度場(chǎng)分布及分析

在目前典型的送風(fēng)條件（送風(fēng)溫度T進(jìn)=340K，進(jìn)風(fēng)口風(fēng)速v進(jìn)=3m/s，相對(duì)濕度=0）下，截取x=2.0m、x=6.0m、y=0m、z=1.6m4個(gè)代表性的截面組成水蒸氣分布圖（圖7）。由圖7可知，在裝煙室頂部及上層煙葉區(qū)域水蒸氣含量極小，水蒸氣質(zhì)量分?jǐn)?shù)為0～0.0003，但在進(jìn)風(fēng)口下部

[FK（W36*3][TPLZY6.TIF][FK）]

區(qū)域含有少量水蒸氣，水蒸氣質(zhì)量分?jǐn)?shù)為0.0014左右。在中層煙葉區(qū)域，裝煙室前部水蒸氣含量略低于后部，但整體分布較為均勻，水蒸氣質(zhì)量分?jǐn)?shù)為0.0011左右。在下層煙葉區(qū)域，水蒸氣分布相對(duì)最不均勻，裝煙室前部水蒸氣質(zhì)量分?jǐn)?shù)為0.0014左右，裝煙室后部水蒸氣質(zhì)量分?jǐn)?shù)為0.0027左右。在裝煙室底部區(qū)域水蒸氣含量最高，并且由裝煙室出風(fēng)口至大門方向逐漸增大，在大門底部區(qū)域水蒸氣含量最大，水蒸氣質(zhì)量分?jǐn)?shù)為0.0038左右。

本研究分析認(rèn)為，裝煙室前部的氣流組織較強(qiáng)，因此熱空氣與煙葉熱交換后蒸發(fā)出的水蒸氣被氣流組織迅速帶走，因此裝煙室前部的水蒸氣含量少于后部。煙葉堆積區(qū)域蒸發(fā)的水蒸氣跟隨氣流組織運(yùn)動(dòng)，上層煙葉蒸發(fā)的水蒸氣跟隨氣流向下層運(yùn)動(dòng)，所以下層水蒸氣含量要高于上層，同時(shí)由于裝煙室前部到后部的氣流組織逐漸減弱，一部分水蒸氣受慣性作用影響滯留在裝煙室底部，在壓差作用下，底部出現(xiàn)遞增的速度梯度，靠近出風(fēng)口一側(cè)的滯留水蒸氣隨著氣流組織向出風(fēng)口運(yùn)動(dòng)，因此底部區(qū)域水蒸氣含量較高，并且該區(qū)域占比由出風(fēng)口一側(cè)向裝煙室后部逐漸增大。

3進(jìn)風(fēng)口風(fēng)速與煙葉間隙風(fēng)速的關(guān)系

3.1進(jìn)風(fēng)口風(fēng)速對(duì)進(jìn)出口壓差的影響

裝煙室進(jìn)風(fēng)口與出風(fēng)口形成的壓差是裝煙室內(nèi)氣流組織流動(dòng)的動(dòng)力來(lái)源，為了研究?jī)烧叩穆?lián)系，在Fluent軟件中改變?nèi)肟谒俣群拖鄬?duì)應(yīng)的邊界條件參數(shù)進(jìn)行模擬。由圖8可知，隨著進(jìn)風(fēng)口風(fēng)速的增加，壓差呈指數(shù)上升趨勢(shì)。當(dāng)進(jìn)風(fēng)口風(fēng)速v≤3m/s時(shí)，隨著進(jìn)口風(fēng)速的增加，壓差緩慢上升，隨著壓差增大，裝煙室內(nèi)氣流組織增強(qiáng)，能有效提高裝煙室內(nèi)的溫度、濕度分布的均勻度。當(dāng)進(jìn)風(fēng)口風(fēng)速v>3m/s時(shí)，隨著進(jìn)風(fēng)口風(fēng)速增加，壓差急劇上升，雖然氣流組織明顯增強(qiáng)，但風(fēng)機(jī)能耗也將隨之上升。

3.2進(jìn)風(fēng)口風(fēng)速對(duì)煙葉間隙風(fēng)速的影響

煙葉間隙的氣流流速?zèng)Q定了水分從煙葉內(nèi)部散失的速度，流速過(guò)大會(huì)增加相同時(shí)間內(nèi)的能源消

耗，使烤后煙葉香氣淡，辛辣味重，刺激性大[6];而流速過(guò)慢則會(huì)影響煙葉的干燥速度，降低烘烤效率。采用“3.1”節(jié)中相同的方式進(jìn)行數(shù)值模擬。由圖9可知，隨著進(jìn)風(fēng)口風(fēng)速的增加，葉間風(fēng)速線性增加，并且葉間風(fēng)速變化幅度較小。宮長(zhǎng)榮等指出，煙葉烘烤干筋期葉間風(fēng)速以0.2～0.3m/s為宜[21]，根據(jù)模擬結(jié)果，其對(duì)應(yīng)的進(jìn)風(fēng)口風(fēng)速為2.2～3.2m/s，綜合考慮風(fēng)機(jī)能耗，進(jìn)風(fēng)口風(fēng)速宜在2.2～3m/s之間選取。

4結(jié)論

基于干筋期密集烤房?jī)?nèi)的熱濕傳遞及氣流組織數(shù)學(xué)模型，采用Fluent軟件對(duì)密集烤房裝煙室內(nèi)溫度場(chǎng)、速度場(chǎng)、濕度場(chǎng)進(jìn)行系統(tǒng)分析，模擬結(jié)果表明：

（1）在裝煙室內(nèi)溫度場(chǎng)、速度場(chǎng)和濕度場(chǎng)三者在同一個(gè)位置的分布是相互關(guān)聯(lián)的。

（2）進(jìn)風(fēng)口風(fēng)速與煙葉間隙風(fēng)速呈線性關(guān)系，并且煙葉風(fēng)速變化受進(jìn)風(fēng)口風(fēng)速影響較小。

（3）綜合考慮風(fēng)機(jī)能耗和葉間風(fēng)速，進(jìn)風(fēng)口風(fēng)速宜在2.2～3m/s之間選取。

參考文獻(xiàn)：

[1]鐵燕，和智君，羅會(huì)龍.煙葉烘烤密集烤房應(yīng)用現(xiàn)狀及展望[J].中國(guó)農(nóng)學(xué)通報(bào)，2009，25（13）：260-262.

[2]段紹米，羅會(huì)龍，劉海鵬.煙葉密集烤房的研究應(yīng)用進(jìn)展與改進(jìn)設(shè)計(jì)[J].江蘇農(nóng)業(yè)科學(xué)，2019，47（1）：202-207.

[3]詹軍，武圣江，賀帆，等.密集烘烤干筋期溫濕度對(duì)上部煙葉外觀質(zhì)量和內(nèi)在品質(zhì)的影響[J].甘肅農(nóng)業(yè)大學(xué)學(xué)報(bào)，2011，46（6）：29-35.

[4]詹軍，宮長(zhǎng)榮，李偉，等.密集烘烤干筋期干球和濕球溫度對(duì)煙葉香氣質(zhì)量的影響[J].湖南農(nóng)業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2011，37（5）：484-489.

[5]宮長(zhǎng)榮，汪耀富，趙銘欽，等.烘烤過(guò)程中煙葉香氣成份變化的研究[J].煙草科技，1995（5）：31-33，27.

[6]白震.烤煙烘烤中干筋期的溫度與香吃味[J].煙草科技，1984（1）：56-60.

[7]沈少君，吳永茂，羅發(fā)健，等.密集烘烤干筋期溫濕度對(duì)煙葉內(nèi)在化學(xué)成分和評(píng)吸質(zhì)量的影響[J].現(xiàn)代農(nóng)業(yè)科技，2014（24）：170-172，181.

[8]王戰(zhàn)義，王剛，劉闖，等.定色期和干筋期不同葉間風(fēng)速對(duì)煙葉香氣質(zhì)量的影響[J].浙江農(nóng)業(yè)科學(xué)，2014（2）：242-244，248.

[9]黃山.烘烤溫濕度變化對(duì)煙葉氧化酶活性及其品質(zhì)的影響[D].長(zhǎng)沙：湖南農(nóng)業(yè)大學(xué)，2010.

[10]JohnDA.Computationalfluiddynamics[M].北京：清華大學(xué)出版社，2002：88-90.

[11]LimTG，KimYH.Analysisofairflowpatterninplantfactorywithdifferentinletandoutletlocationsusingcomputationalfluiddynamics[J].JournalofBiosystemsEngineering，2014，39（4）：310-317.

[12]高維浩，繩春晨，謝洪濤，等.基于Fluent的列管式換熱器性能校核方法研究[J].低溫與超導(dǎo)，2020，48（3）：70-74.

[13]趙海霞，張獻(xiàn)偉.基于Fluent的霧化器霧滴粒徑的研究[J].機(jī)械與電子，2020，38（3）：15-18.

[14]方文康，宋海燕，吳迪.基于Fluent的多溫區(qū)保溫包裝箱溫度場(chǎng)仿真分析[J].包裝工程，2020，41（3）：62-69.

[15]徐陽(yáng)，黃玲.基于CFD的含能材料3D打印機(jī)噴嘴[J].兵工自動(dòng)化，2019，38（6）：57-61.

[16]王娟，王春光，王芳.基于Fluent的9R-40型揉碎機(jī)三維流場(chǎng)數(shù)值模擬[J].農(nóng)業(yè)工程學(xué)報(bào)，2010，26（2）：165-169.

[17]段二亞.糧食干燥過(guò)程的數(shù)值模擬及實(shí)驗(yàn)研究[D].鄭州：河南工業(yè)大學(xué)，2014：40-74.

[18]沙汀鷗，楊潔，呂歡，等.木材干燥室內(nèi)部風(fēng)速場(chǎng)的數(shù)值模擬與優(yōu)化[J].南京林業(yè)大學(xué)學(xué)報(bào)（自然科學(xué)版），2016，40（1）：129-134.

[19]包亞峰.煙葉烘烤過(guò)程中的熱濕分析與優(yōu)化[D].重慶：重慶大學(xué)，2015：22-31.

[20]白志鵬.基于CFD的新型密集烤房流場(chǎng)模擬及其結(jié)構(gòu)優(yōu)化[D].楊凌：西北農(nóng)林科技大學(xué)，2017：27-35.

[21]宮長(zhǎng)榮，周義和，楊煥文.烤煙三段式烘烤導(dǎo)論[M].北京：科學(xué)出版社，2005：72-75，174.