葉大鵬 沈碧河 張炳輝 黃俊煒 謝立敏

摘要：為保障密集烤房裝煙密度，降低烘烤的能源消耗，研發(fā)一套適用于密集烤房的變?nèi)莘e系統(tǒng)。在完成變?nèi)莘e裝置的設(shè)計后，基于CFD方法模擬分析裝置與煙葉的不同距離對烤房內(nèi)部氣體分布均勻性的影響。通過烘烤試驗獲取實際烘烤數(shù)據(jù)，對模擬值加以驗證。試驗結(jié)果表明：當(dāng)隔板與煙葉距離為0cm時，流速不均勻系數(shù)Kv為0.40，溫度不均勻系數(shù)Kt為0.41，距離為10cm時，Kv=0.41，Kt=0.43；距離為20cm時，Kv=0.42，Kt=0.49。裝煙區(qū)9個測量點的溫度模擬值與實測值基本吻合，誤差在6%以內(nèi)。變?nèi)莘e烤房在裝煙量為一半時，相比常規(guī)烤房的燃料消耗可節(jié)約13.4%。研究結(jié)果表明：當(dāng)隔板與煙葉距離為0cm時烤房內(nèi)部的氣體分布最均勻；CFD模型與數(shù)值模擬結(jié)果具有可靠性；變?nèi)莘e裝置具有較好的保溫效果，可保證裝煙密度，降低烤煙能耗。

關(guān)鍵詞：密集烤房；變?nèi)莘e；CFD方法；數(shù)值模擬

中圖分類號：S229+.2? 文獻(xiàn)標(biāo)識碼：A? 文章編號：2095-5553 （2024） 03-0090-07

CFD analysis and experimental research on variable volume bulk curing barn

Ye Dapeng1， 2，?Shen Bihe1，?Zhang Binghui3，?Huang Junwei1，?Xie Limin1， 2

（1. College of Mechanical and Electrical Engineering， Fujian Agricultural and Forestry University， Fuzhou，350100， China; 2. Key Laboratory of Agricultural Information Perception Technology in Fujian Province， Fuzhou，350100， China; 3. Tobacco Science Research Institute of Fujian Tobacco Monopoly Bureau， Fuzhou， 350003， China）

Abstract：

In order to ensure the density of smoke in the bulk curing barn and reduce the energy consumption of baking， a variable volume system suitable for the bulk curing barn was developed. After the design of the variable volume device was completed， the influence of different distances between the device and tobacco leaves on the uniformity of gas distribution inside the bulk cuing barn was simulated and analyzed based on the CFD method. The actual roasting data was obtained through the roasting test， and the simulated values were verified. The test results showed that when the distance between the separator and the tobacco leaf was 0 cm， the non-uniform coefficient of flow velocity Kv was 0.40， the non-uniform coefficient of temperature Kt was 0.41， Kv=0.41， Kt=0.43 when the distanc was 10cm， Kv=0.42， Kt=0.49 when the distanc was 20cm. The simulated temperature values of 9 measurement points in the cigarette loading area were basically consistent with the measured values， with an error of within 6%. When the capacity of the variable volume curing barn was half， the fuel consumption can be saved by 13.4% compared to the conventional curing room. The research results indicated that when the distance between the partition and the tobacco leaves was 0cm， the gas distribution inside the curing room was the most uniform. The CFD model and numerical simulation results had a certain degree of reliability. The variable volume device has good insulation effects， which can ensure the density of tobacco loading and reduce energy consumption of flue-cured tobacco.

Keywords：bulk curing barn; variable volume; CFD method; numerical simulation

0 引言

密集烤房是專門用于卷煙煙葉烘烤的設(shè)備，具有容量大、裝煙密度高等特點［1］?？緹煏r密集烤房的裝煙密度對烤煙的經(jīng)濟(jì)效益有直接關(guān)系，裝煙密度過高或過低對于烤煙的經(jīng)濟(jì)效益都有不良的影響?？痉康哪茉纯傁呐c用電量隨著裝煙密度的提高而增加，但分?jǐn)偟絾挝恢亓扛蔁熒系暮婵境杀緞t相對降低。因此為了節(jié)省能耗，提高經(jīng)濟(jì)效益，密集烤房的裝煙密度不宜過低；但過高的裝煙密度又將對烤后的煙葉等級造成影響，導(dǎo)致均價降低，降低經(jīng)濟(jì)收益［2］。另外，不合適的裝煙密度會對煙葉內(nèi)在化學(xué)成分協(xié)調(diào)性和評吸質(zhì)量帶來不利影響［3］。相關(guān)研究表明，裝煙密度為70kg/m3時下烤后煙葉整體品質(zhì)最高，經(jīng)濟(jì)效益最好［4， 5］。

目前的卷煙生產(chǎn)過程中，鮮煙葉采用逐片采摘烘烤的方式。由于天氣或水肥差異等各方面因素的影響，煙葉成熟速度不同，在一個采收周期內(nèi)收獲的煙葉裝到烤房后會有所剩余或不足，大多數(shù)煙農(nóng)裝煙時只能改變裝煙密度將烤房裝滿，容易導(dǎo)致裝煙密度過高或過低，無法達(dá)到最優(yōu)裝煙密度，從而增加能耗，影響煙葉品質(zhì)，降低烤煙經(jīng)濟(jì)效益。因此，尋求一種密集烤房容積調(diào)節(jié)的方法，可根據(jù)裝載的煙葉數(shù)量調(diào)節(jié)烤房的實際烘烤區(qū)域，使煙葉保持在最佳裝煙密度下烘烤顯得極為必需。

近些年，CFD方法在研究各種農(nóng)業(yè)設(shè)施室內(nèi)的氣體流場分布中被廣泛應(yīng)用［68］。在各種畜禽舍的通風(fēng)換氣與濕熱傳遞分布等問題的研究中，多運用CFD方法進(jìn)行分析，如：運用計算流體動力學(xué)方法對豬舍冬季氣流環(huán)境進(jìn)行研究［9］、對保育豬舍垂直通風(fēng)進(jìn)行CFD模擬及擋風(fēng)板優(yōu)化設(shè)計［10］。一些屋舍結(jié)構(gòu)的優(yōu)化問題也可通過CFD方法得到較好的解決［11， 12］，對密集烤房的諸如通風(fēng)孔尺寸、位置之類的結(jié)構(gòu)優(yōu)化也多使用CFD模擬仿真方法［13， 14］，但有關(guān)烤房的容積調(diào)節(jié)研究卻鮮有報道。

基于以上討論，本文從劃分烤房區(qū)域，隔離裝煙區(qū)，減小加熱空間的角度研究設(shè)計了一套適用于密集烤房的容積調(diào)節(jié)系統(tǒng)。使每次烤煙總能以最佳裝煙密度進(jìn)行烘烤，通過烤煙對比實驗確定隔板厚度，利用CFD三維數(shù)值模擬方法分析隔板與煙葉的不同距離對裝煙區(qū)溫濕度分布均勻性的影響，探究隔板最佳安裝位置，提高烤煙均勻性，并降低能源消耗。

1 材料與方法

1.1 新型密集烤房基本結(jié)構(gòu)與工作原理

目前常用的密集烤房根據(jù)熱空氣流動方向的不同可分為氣流上升式與氣流下降式，二者除風(fēng)機(jī)的位置與朝向不同外結(jié)構(gòu)基本相同［15］。實驗所用烤房是一座位于福建省三明市永安市的氣流上升式密集烤房，烤房可分為燃燒室和烘烤室，烘烤室內(nèi)部尺寸為8.0m×2.7m×3.5m（長×寬×高）。烘烤室的天花板、左右墻面以及大門門板均由5cm厚的聚氨酯夾芯板構(gòu)成，內(nèi)部設(shè)有煙架，用于承載煙桿。烤房工作時燃燒室內(nèi)的氣體被加熱爐加熱，由風(fēng)機(jī)吹向出風(fēng)口，出風(fēng)口出來的熱空氣往大門方向移動的同時慢慢上升，與溫度較低的煙葉進(jìn)行熱交換后到達(dá)天花板，通過回風(fēng)口排出。烤房結(jié)構(gòu)與熱風(fēng)流動路線如圖1所示。

在圖1所示的烤房中搭建變?nèi)莘e裝置，該變?nèi)莘e裝置主要由保溫面板與鋁合金骨架構(gòu)成，六塊面板拼接后形成隔板，搭接在煙架上，將烤房內(nèi)部分隔為裝煙區(qū)與無煙區(qū)，如圖2所示。

為保證裝置的密閉性，隔板的保溫材料尺寸略大于烤房壁面之間的距離，面板互相擠壓減少空隙，因此，材料需要具備一定的可壓縮性和壓縮復(fù)原性；此外，所選取的保溫材料需有較好的隔熱能力以減少裝煙區(qū)的能量損耗。常見的保溫材料及其物理性質(zhì)如表1所示。

氣凝膠氈保溫性能最好但是復(fù)原性差且密度過大；聚氨酯保溫性能次之，但易出現(xiàn)鼓包與開裂；酚醛脆性較大，而聚苯乙烯則吸水性高，均無法滿足隔板的要求。珍珠棉保溫性能較好，壓縮性與壓縮復(fù)原性良好，密度較小但容易彎曲，可以通過骨架加固，使隔板烘烤時不會在風(fēng)壓的作用下發(fā)生變形。綜上，最終選定EPE高密度珍珠棉作為變?nèi)莘e裝置的保溫材料。

變?nèi)莘e裝置的隔板厚度對裝置的保溫性能影響巨大，隔板厚度越大，不同面板之間相互擠壓后的接觸面積越大，密閉性保溫性越好，但裝置的造價也越高。通過測定不同厚度下裝煙區(qū)與無煙區(qū)的溫度差確定面板厚度。

綜合考慮每千克干煙的平均烘烤成本與實際生產(chǎn)情況，研究的裝煙量為最佳裝煙密度條件下裝滿烤房容積一半，即在靠近燃燒室的一側(cè)，每層裝4m長的煙葉。

1.2 烤房建模與網(wǎng)格劃分

根據(jù)現(xiàn)有的烤房結(jié)構(gòu)，等比例建立烘烤室的幾何模型。忽略煙架，觀察窗等所占比重較小的結(jié)構(gòu)。模型由變?nèi)莘e裝置分隔成兩部分，右側(cè)裝煙區(qū)內(nèi)的三層煙葉簡化為三個多孔介質(zhì)區(qū)域，模型結(jié)構(gòu)簡單，在ICEM CFD 19.2中劃分六面體結(jié)構(gòu)化網(wǎng)格并進(jìn)行局部加密，共產(chǎn)生347734個節(jié)點，344667個單元，如圖3所示。

1.3 控制方程與模型選擇

煙葉烘烤過程中烤房內(nèi)部氣體保持流動，整個過程遵循質(zhì)量、動量和能量三大守恒定律［18］；

質(zhì)量守恒方程（連續(xù)性方程）如式（1）所示。

аρаt+а（ρu）аx+а（ρv）аy+а（ρw）аz=Sm（1）

式中：u、v、w——速度矢量在x、y、z方向分量，m/s；ρ——烤房內(nèi)流體密度，kg/m3；t——時間，h。動量守恒方程（N-S方程）如式（2）所示。

ааt（ρui）+ааxi（ρuiuj）=－аpаxi+аτijаxj+ρgi+Fi（2）

式中：p——靜壓，Pa；τij——應(yīng)力張量；gi、Fi——i方向上的重力體積力和外部體積力，N。

能量守恒方程如式（3）所示。

а（ρT）аt+аSymbolQC@（ρuT）=аSymbolQC@kcpgradT+ST（3）

式中：cp——比熱容，J/ （kg·K）；T——溫度，K；k——流體的傳熱系數(shù)；ST——黏性耗散項。

煙葉烘烤的空氣流動屬于湍流，選擇標(biāo)準(zhǔn)的k-ε模型進(jìn)行求解［16］，求解過程中將掛煙區(qū)域視為孔隙率為0.87%［17］的連續(xù)多孔介質(zhì)。

1.4 邊界條件與數(shù)值計算

運用Ansys Fluent軟件的2022R1版本進(jìn)行數(shù)值求解。將烘烤室的熱空氣入口設(shè)置為速度入口，平均風(fēng)速3.5m/s，溫度42℃，湍流模型采用標(biāo)準(zhǔn)的k-ε模型，入口處設(shè)置湍流參數(shù)，湍流強(qiáng)度I=3.7%，水力直徑D=0.7m。出口設(shè)置為常壓狀態(tài)的壓力出口。烤房壁面與天花板設(shè)置為絕熱壁面，不考慮其與外界環(huán)境的換熱，即以壓力基進(jìn)行瞬態(tài)求解，壓力速度耦合方式為SIMPLE，壓力和動量采用二階迎風(fēng)離散格式。

2 試驗與分析

為驗證變?nèi)莘e裝置對烤煙燃料消耗的影響，選取同一烤房群內(nèi)的兩座相鄰烤房進(jìn)行對比試驗。一座為安裝了變?nèi)莘e裝置的試驗烤房，一座為未裝變?nèi)莘e裝置的對比烤房。在兩座烤房內(nèi)部搭建溫濕度監(jiān)測系統(tǒng)，檢測并記錄實時烘烤數(shù)據(jù)。在烤房兩側(cè)壁面與中間煙架三個平面上共布置18個溫濕度探頭，探頭型號為RS485型pvc探頭。在兩座烤房以相同裝煙量，相同開烤時間進(jìn)行烘烤，通過RS-XZJ-100-Y型監(jiān)控主機(jī)每隔10min記錄一次烘烤數(shù)據(jù)并發(fā)送至云平臺。同時記錄兩座烤房的燃料消耗情況。探頭的安裝位置如圖4所示。

烘烤試驗開始的時間為2023年7月4日，于福建省三明市永安市小陶鎮(zhèn)進(jìn)行，烘烤時間為250h。

2.1 隔板厚度的確定

變?nèi)莘e裝置隔板的密閉性和保溫性對于烤房的能源消耗有很大影響。隔板的保溫性越好，裝煙區(qū)的熱量往無煙區(qū)的傳遞越弱，無煙區(qū)的溫度也就越低。無煙區(qū)與外界溫差越低，烤房的熱量損耗越小，節(jié)能效果越明顯。在裝煙量為烤房容積一半的情況下測試隔板厚度為5cm，10cm時的變?nèi)莘e裝置保溫性以確定最佳隔板厚度。

選取以隔板為對稱平面的四個對稱觀測點P6、P7、P14和P16所采集的數(shù)據(jù)進(jìn)行對比。其中，P6、P14位于無煙區(qū)內(nèi)，分別距離隔板約90cm與200cm；P7、P16位于有煙區(qū)內(nèi)，距離隔板約90cm與200cm。烘烤過程中，厚度為5cm的隔板隔離下裝煙區(qū)與無煙區(qū)的溫度變化與溫差如圖5所示。

厚度為10cm的隔板隔離下裝煙區(qū)與無煙區(qū)的溫度變化與溫差如圖6所示。

由圖5和圖6可以看出，隨著烘烤時間的推移，位于有煙區(qū)的P7與P16位置的溫度逐漸上升，在烤房的溫度調(diào)控系統(tǒng)作用下，溫度曲線與烤煙烘烤工藝曲線吻合。位于無煙區(qū)的P6與P14位置的溫度曲線出現(xiàn)周期性波動，溫度下降處的記錄時間為夜晚，應(yīng)是夜晚溫度降低所致。P6與P7、P14與P16的溫度差值也出現(xiàn)周期變化，總體呈現(xiàn)上升趨勢，表明隨著烘烤溫度的上升，隔板兩側(cè)的溫度差呈增大趨勢。P6的溫度始終小于P7，P14始終小于P16，圖5與圖6中的監(jiān)測點溫度差值十分接近，表明10cm的隔板保溫性能與5cm相近。但10cm隔板的厚度增加一倍，設(shè)備成本提升，綜合考慮設(shè)備造價與節(jié)能效果，變?nèi)莘e裝置的隔板厚度定為5cm。

2.2 隔板最佳位置確定

2.2.1 不同隔板位置的裝煙室氣體流速場分布

在模擬仿真的后處理階段，選取烤房內(nèi)距離地面0.35m、1.75m和2.85m的三個水平截面Z1、Z2、Z3與一個垂直地面與隔板的中心截面Y，展示烤房的氣體流速場分布情況。分別創(chuàng)建隔板與煙葉距離為0cm、10cm與20cm的不同模型，三個模型在Z1、Z2、Z3截面上的氣體流速分布如圖7所示。

圖7（a）所示的隔板與煙葉距離為0cm時，空氣從入口進(jìn)入后向上擴(kuò)散，Z1截面上方均為煙葉區(qū)域，截面上各處氣體擴(kuò)散速率相近，氣流撞到隔板后回流，隔板附近風(fēng)速出現(xiàn)周期性波動。Z2截面在隔板一側(cè)風(fēng)速高于燃燒室一側(cè)，高速區(qū)較分散。Z3氣流在往出口運動時速度逐漸增大。

圖7（b）所示的隔板與煙葉距離為10cm時，在Z2截面的氣流高速區(qū)集中的隔板一側(cè)，相比圖7（a）中的Z2截面更集中，應(yīng)是隔板與煙葉區(qū)域存在間隙，進(jìn)入烘烤室的氣流到達(dá)隔板后通過間隙上升至天花板，經(jīng)過Z2截面時往燃燒室運動時速度下降明顯。由于氣體在隔板處快速上升，使得Z3截面的氣流往出口運動的氣體流速增加，平均流速大于圖7（a）。

圖7（c）所示的隔板與煙葉距離為20cm時，隔板與煙葉區(qū)域間隙增大，間隙處的氣體流速增大，Z2截面高速區(qū)更集中，Z3截面平均流速大于圖7（b）。

2.2.2 不同隔板位置的裝煙室溫度場分布

三個模型的不同溫度場分布如圖8所示。

由圖8（a）可以看出，隔板與煙葉距離為0時，熱空氣由進(jìn)風(fēng)口進(jìn)入后，在到達(dá)隔板的過程中，從煙葉區(qū)域的最低端向上擴(kuò)散。在同一水平面各處擴(kuò)散速度相近，使得同一水平高度的煙葉區(qū)域受熱情況相近。圖8（b）中Z1、Z2截面溫度幾乎一致。Z3截面由于出口變窄，在出口附近的角落氣流難以到達(dá)，產(chǎn)生低溫區(qū)，但該區(qū)域無煙葉，影響較小。隨著時間的推移煙葉慢慢達(dá)到與進(jìn)口空氣相近的溫度。

圖8（c）所示隔板與煙葉距離為10cm時，熱空氣在到達(dá)隔板時，由于隔板與煙葉存在空隙，入口的空氣進(jìn)入烘烤室后很快便從空隙處排向出口，溫度傳遞呈“C”型向內(nèi)發(fā)展，煙葉區(qū)域慢慢達(dá)到均勻的溫度，但比入口處溫度略低。圖8（d）中Z1截面溫度均勻，Z2截面在隔板附近的溫度高于燃燒室一側(cè)，Z3截面在裝煙區(qū)的中心區(qū)域溫度較低。

圖8（e）所示隔板與煙葉距離為20cm時，隔板與煙葉空隙增大，熱空氣進(jìn)入烘烤室后更快從出口排出，“C”型溫度傳遞更明顯，由于熱空氣與煙葉接觸時間變短，出口下方的煙葉升溫更慢，出現(xiàn)低溫區(qū)。圖8（f）中Z1與Z3截面溫度均勻，Z2截面從隔板往燃燒室溫度遞減，兩端溫度差異較大。

2.2.3 最佳隔板位置

引入不均勻系數(shù)K［19］描述烤房內(nèi)部各個物理場的不均勻程度，評價隔板位置對烤房環(huán)境的影響。速度場不均勻系數(shù)Kv、溫度場不均勻系數(shù)Kt的公式如式（4）所示。

Kv，t=1m∑mj=11n∑ni=1（Ti－T-）2（4）

式中：j——水平截面數(shù)；i——該水平面上的測量點數(shù)；T-——該平面各個測量點的溫度平均值;m——截面總數(shù);n——同一截面的測量點總數(shù)。K值越小表明均勻程度越好。

在三個Z截面的裝煙區(qū)分別取裝煙區(qū)的對角線交點和兩條對角線上的四個四等分點計算不均勻系數(shù)，不同隔板位置的Kv與Kt見圖9所示。

由圖9可見，隨著隔板位置與煙葉的距離增大，氣體流速的不均勻系數(shù)Kv與溫度不均勻系數(shù)Kt均增大。隔板位置的變化與氣體流速的不均勻系數(shù)曲線較為平緩，位置變化與溫度不均勻系數(shù)曲線陡峭。綜上，隔板距離為0cm時烘烤室的氣體流速場與溫度場均勻程度最好。

2.3 模型可靠性驗證

記錄試驗烤房中裝煙區(qū)的9個監(jiān)測點的烘烤溫度數(shù)據(jù)，以烘烤時間t=68h時的溫度實測數(shù)據(jù)與仿真模擬值進(jìn)行對比，結(jié)果如表2所示。

2.4 燃料消耗分析

在試驗烤房與對比烤房內(nèi)裝煙194桿，占據(jù)烤房裝煙容積的一半，裝煙密度為70kg/m3。以顆粒狀生物質(zhì)燃料作為能量來源。記錄兩座烤房每日的生物質(zhì)燃料消耗，如圖10所示。

烘烤完畢后，兩座烤房的燃料料斗內(nèi)均無燃料剩余。由圖10可以看出，在烘烤過程中，試驗烤房每日的燃料消耗量均小于對比烤房，并且在烘烤后期，試驗烤房的燃料消耗明顯減少。整個烘烤實驗中，對比烤房共計消耗生物質(zhì)燃料820kg，加裝了變?nèi)莘e裝置的試驗烤房共計消耗生物質(zhì)燃料710kg。相較于對比烤房，試驗烤房共節(jié)約燃料110kg，約節(jié)省13.4%。

變?nèi)莘e裝置將烘烤室分隔為裝煙區(qū)與無煙區(qū)，縮小了烤房需要加熱的空間。在烤房沒有滿烤時，從燃燒室進(jìn)入烘烤室的熱空氣在遇到保溫隔板后更快上升，煙葉更快升溫，減少熱空氣在無煙區(qū)運動過程中的能量損耗。另外，變?nèi)莘e裝置使無煙區(qū)的溫度與外界溫差減少，降低了烤房無煙區(qū)與外界的熱交換產(chǎn)生的熱量損失，達(dá)到節(jié)約能源的作用。

3 結(jié)論

1） 本文設(shè)計了一套變?nèi)莘e裝置，并建立了變?nèi)莘e烤房的CFD仿真模型，分析該裝置在新型密集烤房內(nèi)的安裝位置對烤煙室氣流分布均勻性的影響。模擬結(jié)果顯示，當(dāng)隔板與煙葉距離為0cm時，流速不均勻系數(shù)Kv為0.40，溫度不均勻系數(shù)Kt為0.41，10cm時Kv=0.41，Kt=0.43，20cm時Kv=0.42，Kt=0.49，隔板安裝在緊靠煙葉的位置裝煙區(qū)流場分布最佳。

2） 對比烤煙室內(nèi)流場分布的模擬和實際烘烤數(shù)據(jù)測，裝煙區(qū)9個測量點的溫度實測值與模擬值基本吻合，相對誤差＜6%，表明所建立的CFD模型及模擬具有可靠性。

3） 變?nèi)莘e裝置的隔板厚度為5cm，材料為EPE保溫棉可保障裝置的密閉性與保溫性。烘烤試驗結(jié)果表明，烤房裝煙量為50%時，加裝了變?nèi)莘e裝置的烤房比常規(guī)烤房節(jié)約13.4%的燃料消耗。

參 考 文 獻(xiàn)

［1］劉光輝， 聶榮邦. 我國烤房及烘烤技術(shù)研究進(jìn)展［J］. 作物研究， 2011， 25（1）： 76-80.

［2］黃建， 張海偉， 盧瑞杰. 密集烤房裝煙密度對煙葉質(zhì)量及烘烤效益的影響［J］. 湖北農(nóng)業(yè)科學(xué)， 2014， 53（7）： 1595-1598.Huang Jian， Zhang Haiwei， Lu Ruijie. Effects of loading densities on leaf qualities and benefits of bulk cured tobacco ［J］. Hubei Agricultural Science， 2014， 53（7）： 1595-1598.

［3］張?zhí)? 密集烤房裝煙密度對煙葉質(zhì)量的影響［D］. 長沙： 湖南農(nóng)業(yè)大學(xué)， 2013.Zhang Tao. The effect of tobacco density in dense curing on tobacco leaf quality ［D］. Changsha： Hunan Agricultural University， 2013.

［4］賀非， 皮本陽， 孫敬釗， 等. 裝煙密度對烤后煙葉品質(zhì)的影響［J］. 農(nóng)業(yè)工程技術(shù)， 2020， 40（26）： 32-33.

［5］詹軍， 李偉， 武圣江， 等. 密集烤房裝煙密度對烤煙上部煙葉質(zhì)量的影響［J］. 作物研究， 2011（5）： 473-476.Zhang Jun， Li Wei， Wu Shengjiang， et al. Effects of loading density of bulk curing barn on quality of upper leaves in flue-cured tobacco ［J］. Crop Research， 2011（5）： 473-476.

［6］吳照學(xué)， 梁偉， 鮑恩財， 等. 基于CFD的發(fā)酵床網(wǎng)養(yǎng)鴨舍冬季環(huán)境模擬與優(yōu)化［J］. 中國農(nóng)機(jī)化學(xué)報， 2023， 44（12）： 33-39.Wu Zhaoxue， liang Wei， Bao Encai， et al， Simulation and optimization of winter environment for duck house with fermentation bed net based on CFD ［J］. Journal of Chinese Agricultural Mechanization， 2023， 44（12）： 33-39.

［7］汪開英， 李開泰， 李王林娟， 等. 保育舍冬季濕熱環(huán)境與顆粒物CFD模擬研究［J］. 農(nóng)業(yè)機(jī)械學(xué)報， 2017， 48（9）： 270-278.Wang Kaiying， Li Kaitai， Li Wanglinjuan， et al. CFD simulation of indoor hygrothermal environment and particle matter of weaned pig building ［J］. Transactions of the Chinese Society for Agricultural Machinery. 2017， 48（9）： 270-278.

［8］吳志東， 房俊龍， 巴文革， 等. 冬季采暖保育豬舍送排風(fēng)管道組合換氣系統(tǒng)設(shè)計與評價［J］. 農(nóng)業(yè)工程學(xué)報， 2021， 37（10）： 152-158.Wu Zhidong， Fang Junlong， Ba Wenge， et al. Design and evaluation of combined ventilation system with air supply and exhaust ducts for a nursery heating piggery in winter ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2021， 37（10）： 152-158.

［9］周麗娜， 徐述胤， 馬精陽， 等. 基于計算流體力學(xué)的豬舍冬季氣流環(huán)境模擬與優(yōu)化研究［J］. 黑龍江畜牧獸醫(yī)， 2022（7）： 49-54， 59， 127-128.Zhou Lina， Xu Shuyin， Ma Jingyan， et al. Simulation and optimization of winter airflow environment in piggery based on CFD ［J］. Heilongjiang Animal Science and Veterinary Medicine， 2022（7）： 49-54， 59， 127-128.

［10］房俊龍， 吳爽， 吳志東， 等. 保育豬舍垂直通風(fēng)CFD模擬及擋風(fēng)板優(yōu)化設(shè)計［J］. 東北農(nóng)業(yè)大學(xué)學(xué)報， 2022， 53（5）： 59-68.Fang Junlong， Wu Shuang， Wu Zhidong， et al. CFD simulation of vertical ventilation and optimal design of wind shield in nursery pig house ［J］. Journal of Northeast Agricultural University， 2022， 53（5）： 59-68.

［11］李國建， 楊林， 許金釵， 等. 基于CFD的籠養(yǎng)雞舍冬季通風(fēng)方案模擬與驗證研究［J］. 黑龍江畜牧獸醫(yī)， 2021（21）： 57-62， 151-153.Li Guojian， Yang Lin， Xu Jinchai， et al .Simulation and verification studies of ventilation schemes in stacked-cage hen house based on CFD in winter ［J］. Heilongjiang Animal Science and Veterinary Medicine， 2021（21）： 57-62， 151-153.

［12］尚峰軍， 趙潔. 采用CFD技術(shù)分析和評估蛋雞舍溫?zé)岘h(huán)境質(zhì)量［J］. 中國家禽， 2021， 43（2）： 61-68.Shang Fengjun， Zhao Jie. Using CFD technology to analyze and evaluate the quality of the warm environment of layer house ［J］. Chinese Poultry， 2021， 43（2）： 61-68.

［13］Martínez-Martínez V， Baladrón C， Gomez-Gil J， et al. Temperature and relative humidity estimation and prediction in the tobacco drying process using artificial neural networks ［J］. Sensors， 2012， 12， 14004-14021.

［14］Bai Zhipeng， Guo Duoduo， Li Shoucang， et al. Analysis of temperature and humidity field in a new bulk tobacco curing barn based on CFD ［J］. Sensors， 2017， 17（2）： 279-279.

［15］詹軍， 樊軍輝， 宋朝鵬， 等. 密集烤房研究進(jìn)展與展望［J］. 南方農(nóng)業(yè)學(xué)報， 2011， 42（11）： 1406-1411.Zhang Jun， Fan Junhui， Song Chaopeng， et al. Comparative analysis of bulk curing barn and traditional flue-curing barn ［J］. Journal of Southern Agriculture， 2011， 42（11）： 1406-1411.

［16］劉兆宇， 羅會龍， 陳頤， 等. 煙葉烘烤干筋期密集烤房熱濕環(huán)境數(shù)值分析［J］. 江蘇農(nóng)業(yè)科學(xué)， 2021， 49（1）：160-165.

［17］孫力.基于水分遷移模型變孔隙率密集烤房煙葉干燥數(shù)值模擬及實驗分析［D］. 昆明： 昆明理工大學(xué)， 2020.Sun Li. Numerical simulation and experimental analysis of tobacco leaf drying in the closed bulk curing barn， based on variable porosity and moisture migration model ［D］. Kunming： Kunming University of Science and Technology， 2020.

［18］王福軍. 計算流體動力學(xué)分析—CFD軟件原理與應(yīng)用［M］. 北京： 清華大學(xué)出版社， 2004.

［19］鄧書輝， 施正香， 李保明， 等. 低屋面橫向通風(fēng)牛舍空氣流場CFD模擬［J］. 農(nóng)業(yè)工程學(xué)報， 2014， 30（6）： 139-146.Deng Shuhui， Shi Zheng Xiang， Li Baoming， et al. CFD simulation of airflow distribution in low profile cross ventilated dairy cattle barn ［J］. Transactions of the Chinese Society of Agricultural Engineering， 2014， 30（6）： 139-146.

基金項目：福建省農(nóng)業(yè)信息感知技術(shù)重點實驗室項目（2021ZDSYS0101）；中國煙草總公司福建省公司科技計劃項目（2022350000240084）

第一作者：葉大鵬，男，1971年生，福建寧德人，博士，教授；研究方向為農(nóng)業(yè)智能機(jī)械。E-mail： ydp@fafu.edu.cn

通訊作者：謝立敏，女，1985年生，福建南平人，博士，講師；研究方向為機(jī)器人運動控制。E-mail： lucy_min@163.com