郭文斌，李瑤，黃長華，杜建強(qiáng)*，錢珊珠，何澤民，高晶晶

（1.內(nèi)蒙古農(nóng)業(yè)大學(xué)機(jī)電工程學(xué)院，呼和浩特 010018；2.內(nèi)蒙古商貿(mào)職業(yè)學(xué)院，呼和浩特 010070）

新鮮草料收獲后含水率較高，不及時進(jìn)行加工處理容易造成養(yǎng)分損失，甚至腐爛、發(fā)酵。為解決農(nóng)業(yè)物料收獲后含水率高、易腐敗、不易儲藏的問題，科研人員嘗試了多種方法對其進(jìn)行干燥，如紅外線干燥、微波干燥、高壓靜電場干燥、太陽能熱風(fēng)干燥等。通過比較干燥效果，發(fā)現(xiàn)太陽能熱風(fēng)干燥不僅具有能耗小、效率高、成本低的特點(diǎn)[1]，而且還可以在干燥過程中減少牧草、果蔬等農(nóng)業(yè)物料葉片及嫩枝的脫落和營養(yǎng)物質(zhì)的損失，保證其干燥品質(zhì)[2]，因此被廣泛應(yīng)用于中藥材、牧草、糧食等的干燥。太陽能熱風(fēng)干燥以空氣作為傳熱介質(zhì)，通過空氣集熱器加熱后，由鼓風(fēng)機(jī)送入干燥箱內(nèi)，利用流動的熱空氣帶走物料水分，完成干燥作業(yè)。

在太陽能熱風(fēng)干燥過程中，熱空氣風(fēng)速流動往往受干燥設(shè)備內(nèi)部結(jié)構(gòu)影響，呈現(xiàn)不均勻分布，導(dǎo)致物料干燥效果不佳。干燥箱是太陽能干燥系統(tǒng)的主要組成部件，其內(nèi)部結(jié)構(gòu)是影響熱風(fēng)流動均勻性的關(guān)鍵因素之一，熱風(fēng)干燥箱結(jié)構(gòu)設(shè)計優(yōu)化[3-5]及相關(guān)農(nóng)業(yè)物料干燥性能的研究[6-8]成為關(guān)注的熱點(diǎn)。目前熱風(fēng)干燥箱的設(shè)計優(yōu)化主要針對箱體結(jié)構(gòu)，如內(nèi)部空間結(jié)構(gòu)、物料承載裝置、進(jìn)出風(fēng)口結(jié)構(gòu)、箱門及箱壁結(jié)構(gòu)等[9-13]，一定程度上提高了物料干燥的效率，而通過物理場模擬分析完成干燥箱結(jié)構(gòu)優(yōu)化和物料干燥均勻性的研究較少。本文針對小型熱風(fēng)干燥箱結(jié)構(gòu)，在驗證分析的基礎(chǔ)上，基于COMSOL Multiphysics 物理場數(shù)值模擬方法，提出結(jié)構(gòu)優(yōu)化方案，以提高干燥均勻性和改善熱風(fēng)干燥箱干燥性能。

1 材料與方法

1.1 試驗設(shè)備

為完成熱風(fēng)干燥箱內(nèi)物理場分布的測試與分析，搭建了太陽能熱風(fēng)干燥系統(tǒng)試驗平臺，由熱風(fēng)干燥箱、空氣集熱器、鼓風(fēng)機(jī)、TBQ-2 型太陽能輻射記錄儀、PC-2WS 型多通道溫濕度記錄儀及SHT10 溫濕度傳感器、FL-J 風(fēng)速傳感器等組成（圖1），傳感器及記錄儀均購自錦州陽光氣象科技有限公司。其中熱風(fēng)干燥箱由帶有保溫夾層的側(cè)壁和箱門組成，內(nèi)部設(shè)有進(jìn)風(fēng)口和出風(fēng)口，進(jìn)風(fēng)口位于干燥箱底部，通過通風(fēng)管與鼓風(fēng)機(jī)和空氣集熱器連接，出風(fēng)口位于箱體頂部與外界相通，進(jìn)、出風(fēng)口之間設(shè)置有承載物料的透風(fēng)托盤。試驗時，除進(jìn)、出風(fēng)口外應(yīng)保證干燥箱密閉，由鼓風(fēng)機(jī)將集熱器內(nèi)加熱的空氣從進(jìn)風(fēng)口送入干燥箱內(nèi)，試驗期間溫度為（40±1）℃。完成測試后，空氣經(jīng)出風(fēng)口流出。

圖1 太陽能熱風(fēng)干燥系統(tǒng)試驗平臺Fig.1 Solar hot-air drying system test-bed

通過在透風(fēng)托盤上安裝風(fēng)速傳感器和溫濕度傳感器，測得干燥箱體內(nèi)各位置點(diǎn)處風(fēng)速及溫度，每組試驗測試時間為30 min，試驗數(shù)據(jù)采集頻率為每分鐘1次。衡量干燥箱內(nèi)風(fēng)速場是否均勻時，引入不均勻系數(shù)（M）作為試驗指標(biāo)[14]，M越大，各測點(diǎn)所在截面的風(fēng)速分布越不均勻，其計算公式如下。

式中，Vi為干燥箱某截面上各測點(diǎn)風(fēng)速，m·s-1；i為各個測點(diǎn)的數(shù)字編號；-V為同截面上各測點(diǎn)風(fēng)速的均值，m·s-1；n為測點(diǎn)數(shù)。

通過調(diào)控鼓風(fēng)機(jī)設(shè)置3 種干燥箱進(jìn)風(fēng)口處風(fēng)速：5、7、10 m·s-1；干燥箱內(nèi)部設(shè)置3 層透風(fēng)托盤，自下而上依次為第1層、第2層、第3層，距離進(jìn)風(fēng)口截面的層高依次為150、450、750 mm（圖2），以每層托盤所在平面為計算風(fēng)速不均勻系數(shù)M的截面，每層截面設(shè)計12 個測點(diǎn)，用于放置風(fēng)速傳感器和溫濕度傳感器，測點(diǎn)編號及位置如圖3所示。

圖2 熱風(fēng)干燥箱結(jié)構(gòu)Fig.2 Structural of hot-air drying box

圖3 干燥箱內(nèi)傳感器位置Fig.3 Positions of sensors drying box

干燥箱內(nèi)溫度是否分布均勻主要通過溫度分布不均勻系數(shù)（Tv）進(jìn)行描述，計算如公式（2）所示。Tv越小，各點(diǎn)的溫度值越接近于整體均值，當(dāng)Tv＜2%時，表明溫度場分布較均勻[15-16]。

式中，Ti為各測點(diǎn)溫度，℃為各測點(diǎn)溫度平均值，℃。

1.2 干燥箱物理場的COMSOL模擬

為改善熱風(fēng)干燥箱內(nèi)風(fēng)速場分布的均勻程度，基于計算流體力學(xué)完成其結(jié)構(gòu)優(yōu)化與改進(jìn)，在COMSOL 環(huán)境下對干燥箱內(nèi)風(fēng)速場、溫度場進(jìn)行模擬仿真，并與試驗所得物理場分布規(guī)律進(jìn)行驗證，確定合適的仿真分析條件與參數(shù)設(shè)置方法。

考慮到物理場分布計算的復(fù)雜性，為提高COMSOL 仿真環(huán)境下的計算效率和準(zhǔn)確程度，建立熱風(fēng)干燥箱仿真模型時簡化如下[14]：將干燥箱箱壁視為絕熱體，且壁面無滑移；除進(jìn)、出風(fēng)口外，干燥箱其他位置視為密閉；將空氣視為不可壓縮的連續(xù)介質(zhì)。

根據(jù)試驗時的條件和要求，對簡化后的干燥箱模型邊界條件參數(shù)進(jìn)行設(shè)置（表1）。

表1 模型及邊界條件參數(shù)設(shè)置Table 1 Parameters of boundary conditions setting

在設(shè)置干燥箱內(nèi)流場類型時，結(jié)合試驗所得熱風(fēng)流動情況，選擇較為普遍的湍流流場[17]，模型類型為RANS（reynolds-averaged navier-stokes equations），采用k-ε模型模擬干燥箱內(nèi)熱風(fēng)流動過程，通過湍動能與耗散率方程得到湍流動能k、耗散率ε，如公式（3）和（4）所示。

式中，ρ為流體密度，kg·m-3；σk和σε分別為湍流動能和耗散率的普朗特數(shù)；u為平均速度，m·s-1；Pk為平均速度改變引起的湍動能生成項；μ為動力黏度，μT為湍流黏度，Pa·s；Cε1、Cε2為常數(shù)項；ρ=1.132，μT=17.9×10-6，Cε1=1.44，Cε2=1.92。

完成參數(shù)設(shè)置后，以COMSOL 中物理場控制的網(wǎng)格劃分方式對干燥箱簡化模型進(jìn)行極粗化網(wǎng)格劃分，以完成其內(nèi)物理場的數(shù)值模擬分析。

2 結(jié)果與分析

2.1 干燥箱風(fēng)速場和溫度場分析

進(jìn)風(fēng)口風(fēng)速為5 m·s-1時，3 層托盤上各測點(diǎn)風(fēng)速分布如圖4 所示?？梢钥闯觯捎谶M(jìn)、出風(fēng)口分別位于箱體底面和頂面的中心，且處于正對位置，使得各層中心處測點(diǎn)5、測點(diǎn)8 的風(fēng)速顯著高于周圍靠近箱壁處的其他測點(diǎn)。此外，受進(jìn)風(fēng)口下方聯(lián)接管位置、導(dǎo)向誤差的影響，箱壁附近其他測點(diǎn)的風(fēng)速大小也略有差異。

圖4 干燥箱內(nèi)各層風(fēng)速分布Fig.4 Wind velocity distribution at different layers of drying box

各層風(fēng)速場不均勻系數(shù)M如表2 所示，可以看出，不同進(jìn)風(fēng)口風(fēng)速條件下，M介于63.31%～86.32%，表明所測試干燥箱的內(nèi)部風(fēng)速場分布均勻性較差。

從進(jìn)風(fēng)口風(fēng)速為5、7、10 m·s-1時各層測點(diǎn)溫度（表2）可以看出，受環(huán)境溫度變化影響，3 種風(fēng)速下的溫度盡管有所不同，但同一風(fēng)速下，干燥箱內(nèi)各層溫度分布的不均勻系數(shù)Tv均小于2%，說明干燥箱工作時，箱內(nèi)各測點(diǎn)位置處溫度差異較小，溫度場分布比較均勻，并未受到風(fēng)速場分布不均勻的影響。

表2 干燥箱內(nèi)各層風(fēng)速、溫度分布及其不均勻系數(shù)Table 2 Distribution of wind velocity and temperature and their non-uniformity coefficients in different layers of drying box

上述結(jié)果表明，受干燥箱結(jié)構(gòu)影響，熱風(fēng)從底部進(jìn)入箱體內(nèi)部后流速并不均勻，風(fēng)速場呈現(xiàn)出中心高、四周低的分布規(guī)律，盡管該分布規(guī)律并未造成干燥箱內(nèi)溫度場的分布不均勻，但不同位置熱風(fēng)流速不同，單位時間內(nèi)從待干燥物料中帶走的水分也不同，從而對熱風(fēng)干燥箱的干燥性能產(chǎn)生影響，導(dǎo)致物料干燥效果較差。因此對干燥箱進(jìn)行結(jié)構(gòu)優(yōu)化時需以風(fēng)速場分布是否均勻為衡量指標(biāo)。

2.2 模擬與驗證結(jié)果分析

對干燥箱內(nèi)風(fēng)速場進(jìn)行COMSOL 模擬，得到各層托盤所在橫截面處的風(fēng)速分布如圖5 所示，可以看出，在距進(jìn)風(fēng)口150 mm 的層高處，由于熱風(fēng)剛進(jìn)入干燥箱內(nèi)，受進(jìn)風(fēng)口位置影響，流動較快的區(qū)域集中在靠近箱體中心的兩處對稱位置；而層高達(dá)到450 mm 時，與進(jìn)風(fēng)口位置對應(yīng)的2 處熱風(fēng)流速較快區(qū)域已逐漸從各自中心擴(kuò)散、合并；當(dāng)熱風(fēng)流至層高750 mm 處，風(fēng)速最高的區(qū)域由2 處合并為位于干燥箱中心的1 處。該分布規(guī)律表明，熱風(fēng)從干燥箱底部進(jìn)風(fēng)口進(jìn)入時風(fēng)速最高、最集中，隨后向上流動，風(fēng)速開始下降，并逐漸由箱體中心向四周緩慢擴(kuò)散，流至層高750 mm 處時，由2 個進(jìn)風(fēng)口進(jìn)入的熱風(fēng)已逐漸擴(kuò)散合并，到達(dá)出風(fēng)口時沿出風(fēng)口位置分離并流出。上述仿真風(fēng)速場分布特點(diǎn)與試驗所得風(fēng)速場分布特點(diǎn)相似，即現(xiàn)有干燥箱結(jié)構(gòu)內(nèi)部呈現(xiàn)中心風(fēng)速高、四周風(fēng)速低的規(guī)律。

圖5 干燥箱內(nèi)橫向截面風(fēng)速分布Fig.5 Distribution of wind velocity in the transverse section of drying box

為進(jìn)一步驗證干燥箱內(nèi)部風(fēng)速場模擬的準(zhǔn)確性，計算得到不同進(jìn)風(fēng)口風(fēng)速下各層檢測點(diǎn)的風(fēng)速不均勻系數(shù)M（表3），可以看出，模擬所得各層高度處的風(fēng)速分布不均勻系數(shù)M均高達(dá)55%以上，即仿真條件下干燥箱各層的風(fēng)速場分布均勻性均較差，與試驗所得結(jié)果一致。

表3 COMSOL環(huán)境下干燥箱各層風(fēng)速場不均勻系數(shù)Table 3 Non-uniformity coefficient of wind velocity field of drying box under the COMSOL model

表4 展示了不同進(jìn)風(fēng)口風(fēng)速條件下層高450 mm 處各測點(diǎn)風(fēng)速的試驗實(shí)測值與仿真模擬值，可以看出，模擬所得風(fēng)速分布規(guī)律與試驗所得結(jié)果相同，即中心測點(diǎn)5、8 處風(fēng)速較高，其余測點(diǎn)風(fēng)速相對較低。受箱體密封性的影響，測點(diǎn)3 的實(shí)測值與模擬值誤差相對較大，此外，其余測點(diǎn)的實(shí)測值與模擬值相對誤差均較?。ㄆ骄`差≤5.02%），因此模擬分析所得數(shù)據(jù)具有較高準(zhǔn)確性。對干燥箱內(nèi)溫度場分布進(jìn)行模擬后發(fā)現(xiàn)，各測點(diǎn)處的溫度基本相同，溫度場分布較試驗所得結(jié)果更為均勻理想(Tv＜0.001%)，并未受到風(fēng)速分布不均勻的影響，該結(jié)果與試驗測試結(jié)果一致。

表4 層高450 mm不同進(jìn)風(fēng)口風(fēng)速下各測點(diǎn)風(fēng)速實(shí)測值與模擬值對比Table 4 Comparison of measured and simulated wind velocities at different measurement points at 450 mm layer

上述模擬分析結(jié)果較準(zhǔn)確地印證了試驗所得風(fēng)速場、溫度場的分布規(guī)律。說明基于COMSOL仿真平臺進(jìn)行熱風(fēng)干燥箱內(nèi)部風(fēng)速場模擬能夠反映出干燥箱內(nèi)物理場的分布特點(diǎn)。因此，在驗證該模擬分析方法的基礎(chǔ)上，將其用于改進(jìn)結(jié)構(gòu)后干燥箱風(fēng)速場模擬，以降低干燥箱樣機(jī)生產(chǎn)和試驗測試帶來的研發(fā)成本、縮短設(shè)計周期。

2.3 干燥箱結(jié)構(gòu)改進(jìn)與風(fēng)速場分析

針對熱風(fēng)干燥箱作業(yè)時中心位置風(fēng)速高、靠近箱壁位置風(fēng)速低的特點(diǎn)，對箱體結(jié)構(gòu)進(jìn)行了重新設(shè)計與改進(jìn)，并依托COMSOL 仿真平臺對結(jié)構(gòu)改進(jìn)后的風(fēng)速場進(jìn)行模擬分析，以提升箱體內(nèi)風(fēng)速場分布的均勻程度，改善干燥箱的干燥性能。

2.3.1 干燥箱結(jié)構(gòu)改進(jìn)設(shè)計方案 對試驗及模擬驗證結(jié)果的分析表明，由于進(jìn)、出風(fēng)口均位于干燥箱頂面和底面的中心位置，直接導(dǎo)致熱風(fēng)進(jìn)入干燥箱后，中心風(fēng)速高、四周風(fēng)速低；而且由于進(jìn)風(fēng)口下方空間有限，使得與其連接的通風(fēng)管位置及導(dǎo)流方向容易對熱風(fēng)流向造成干擾。鑒于以上原因，在進(jìn)行干燥箱改進(jìn)時，首先將進(jìn)風(fēng)口位置調(diào)整至干燥箱側(cè)壁底部，使熱風(fēng)橫向進(jìn)入箱體內(nèi)部，進(jìn)風(fēng)口位置變化后箱體內(nèi)風(fēng)速場分布如圖6 所示，雖然干燥箱中心風(fēng)速快的情況有所緩解，但箱體內(nèi)風(fēng)速場分布仍不均勻，出現(xiàn)了熱風(fēng)沿著進(jìn)風(fēng)口對面箱壁向上流動較快的現(xiàn)象。

圖6 過進(jìn)風(fēng)口中心縱向截面風(fēng)速分布Fig.6 Wind velocity distribution in the vertical section through inlet center

在改變進(jìn)風(fēng)口位置基礎(chǔ)上，繼續(xù)在干燥箱底部透風(fēng)托盤下方正對進(jìn)風(fēng)口位置處增加3 層可調(diào)節(jié)帶孔擋流板，構(gòu)建隔斷式氣室，結(jié)構(gòu)如圖7 所示。在對擋流板間距、高度、開孔密度等逐一進(jìn)行仿真模擬后，最終采用了設(shè)置等高度、等間隔擋流板的設(shè)計方案，該方案中相鄰擋流板上開有位置高度不同的兩排圓孔（孔徑110 mm，各擋流板開孔高度隨熱風(fēng)流向依次增加30 mm），以便于熱風(fēng)進(jìn)入后能夠部分沿著擋流板向上流動、部分通過擋流板上圓孔進(jìn)入下一間隔區(qū)域，最終將熱風(fēng)分隔于各擋流板間區(qū)域，使其均勻向上流動，在干燥箱內(nèi)獲得分布較為均勻的風(fēng)速場。

圖7 含隔斷式氣室的干燥箱Fig.7 Drying box with a partition-type

2.3.2 改進(jìn)后仿真結(jié)果分析 對改進(jìn)設(shè)計后的具有隔斷式氣室的干燥箱進(jìn)行內(nèi)部風(fēng)速場仿真模擬結(jié)果如圖8 和9 所示（進(jìn)風(fēng)口風(fēng)速為5 m·s-1），熱風(fēng)經(jīng)干燥箱底部橫向進(jìn)入后，一部分受到擋流板的阻擋向上流動，另一部分則從開孔位置流向下一層擋流板處，隨后再次被分流，最終熱風(fēng)在擋流板作用下得到逐層分配，在箱體內(nèi)形成了較為均勻的風(fēng)速場。不同進(jìn)風(fēng)口風(fēng)速條件下，擋流板上方區(qū)域各層橫截面的風(fēng)速場分布不均勻系數(shù)如表5 所示，增加隔斷式氣室后干燥室內(nèi)的風(fēng)速場不均勻系數(shù)M較結(jié)構(gòu)改進(jìn)前（表3）明顯降低，且層高越高M(jìn)值越小，其中層高在1 100 mm 之上時，M值降至10%以下。由此可見，帶孔擋流板的設(shè)置起到了很好的勻風(fēng)作用，改善了干燥箱內(nèi)部風(fēng)速場的均勻性。

表5 不同進(jìn)風(fēng)口風(fēng)速下各層截面風(fēng)速場分布不均勻系數(shù)Table 5 Non-uniformity coefficients of wind velocity field under different inlet wind velocities

圖8 改進(jìn)后干燥箱縱向截面風(fēng)速分布Fig.8 Wind velocity distribution in the vertical section of improved drying box

3 討論

在試驗測試的基礎(chǔ)上，通過COMSOL 軟件完成了熱風(fēng)干燥箱內(nèi)風(fēng)速場和溫度場的數(shù)值模擬，得到了風(fēng)速和溫度的分布規(guī)律及不均勻系數(shù)，與試驗所得結(jié)果對比分析后，發(fā)現(xiàn)模擬所得結(jié)果與試驗一致，驗證了該數(shù)值模擬方法的準(zhǔn)確性和可行性。

基于COMSOL 數(shù)值模擬仿真分析方法，以降低風(fēng)速場分布不均勻系數(shù)為目標(biāo)，提出了干燥箱結(jié)構(gòu)優(yōu)化與改進(jìn)設(shè)計方案：改變進(jìn)風(fēng)口位置，實(shí)現(xiàn)側(cè)向進(jìn)風(fēng)；添加帶孔擋流板、構(gòu)建隔斷式氣室，實(shí)現(xiàn)干燥箱內(nèi)熱風(fēng)流向和流速的調(diào)節(jié)。模擬結(jié)果顯示改進(jìn)后干燥箱的風(fēng)速場分布有了明顯改善，不均勻系數(shù)下降至了10%以內(nèi)。該研究結(jié)果為熱風(fēng)干燥箱的結(jié)構(gòu)改進(jìn)優(yōu)化提供了新的技術(shù)手段和數(shù)據(jù)支撐，降低了新型干燥箱的研發(fā)成本。

圖9 改進(jìn)后干燥箱橫向截面風(fēng)速分布Fig.9 Wind velocity distribution in the transverse section of improved drying box