太陽能熱泵干燥室氣流組織優(yōu)化研究

,,

(1.上海理工大學(xué) 食品科學(xué)與工程研究所,上海 200093；2.上海筑能環(huán)境科技有限公司,上海 201299)

0 引言

青海省格爾木地處青藏高原腹地，海拔高屬典型的大陸性氣候。晝夜溫差大、少雨、日照時間長、年日均太陽能輻照量達21.93 MJ/m2，光熱資源充足。這些有利因素決定了青海枸杞的產(chǎn)品品質(zhì)、營養(yǎng)成分及口感外觀的優(yōu)良，同時枸杞生產(chǎn)企業(yè)迅猛發(fā)展。枸杞種植規(guī)模由2007年的1 580 hm2增長到2015年的2.81萬hm2，2015年生產(chǎn)枸杞干果6.2萬t以上[1]。

枸杞干制多采用成本較低的自然曬干或熱風(fēng)烘干，自然曬干受天氣影響，枸杞的有效成分損失較大，不能保證干燥成品品質(zhì)。熱風(fēng)烘干是將空氣作為干燥介質(zhì)，加熱后送入干燥室，對枸杞進行干燥，該干燥方式具有對枸杞損傷小，裝載量大，干燥室投資成本低等優(yōu)點[2]。在此基礎(chǔ)上又出現(xiàn)了太陽能熱泵干燥機，干燥枸杞具有節(jié)能環(huán)保，系統(tǒng)運行效率高，物料有效成分損失少等特點。但為提高生產(chǎn)效率，增大干燥室，帶來了干燥不均勻等問題。

干燥室內(nèi)氣流組織的不均勻性造成設(shè)備的能源消耗增加、干燥效率低和物料干燥品質(zhì)差。翟大海[3]等對室內(nèi)空氣流場進行數(shù)值模擬分析，采用了四種不同送風(fēng)方式，分析了溫度場分布、速度場分布等參數(shù)，結(jié)果顯示上送下回的送風(fēng)方式氣流均勻性好。張仲欣[4]等對食品廠在橫流循環(huán)式烘干機進行氣流優(yōu)化，發(fā)現(xiàn)影響干燥均勻性的主要因素是物料通風(fēng)層內(nèi)氣流分布不均勻，出現(xiàn)風(fēng)速梯度，從而導(dǎo)致溫度差增加，產(chǎn)生干燥品質(zhì)不均勻現(xiàn)象。

Zhang[5]等人對蔬菜食品室內(nèi)采用強制對流換氣循環(huán)系統(tǒng)，設(shè)計具有三排空氣射流穿孔空氣管，并提供了四種不同的管道布置方式，改善了室內(nèi)氣流均勻性，解決了蔬菜生長不均勻現(xiàn)象。Majka[6]等人研究風(fēng)速對木材干燥均勻性的影響，研究表明：風(fēng)扇轉(zhuǎn)速的降低提高了木材干燥的均勻性；但風(fēng)扇速度過度減少，降低到標準水平以下會影響空氣擴散，同時也會導(dǎo)致干燥條件變差。

本文采用送風(fēng)均勻性較強的孔板送風(fēng)方式，研究風(fēng)速、孔徑、風(fēng)量因素對干燥室氣流不均勻性的影響，為確定各因素對氣流組織的影響，改善干燥室內(nèi)氣流組織不均勻，采用CFD數(shù)值模擬技術(shù)并結(jié)合正交實驗設(shè)計方法[7]，找出最佳的優(yōu)化方案，為干燥室氣流優(yōu)化提供參考。其中，物料形體大小在干燥過程中會影響氣流分布，本文主要分析孔板送風(fēng)方式對氣流均勻性的影響，忽略干燥物料形體對氣流的影響。

1 太陽能熱泵干燥系統(tǒng)

設(shè)計的太陽能熱泵干燥系統(tǒng)如圖1所示，分為太陽能集熱系統(tǒng)、熱泵系統(tǒng)、熱回收系統(tǒng)以及干燥室四部分。

干燥系統(tǒng)利用太陽能集熱系統(tǒng)收集太陽輻射轉(zhuǎn)化為熱量送入干燥室內(nèi)進行干燥；當太陽能滿足不了負荷時，由空氣源熱泵機組提供干燥熱負荷，保證系統(tǒng)的穩(wěn)定運行；該干燥系統(tǒng)還針對干燥室排出余熱高濕空氣，通過余熱回收器對新風(fēng)進行預(yù)熱[8]，從而避免能源浪費，提高系統(tǒng)運行效率。

空氣式太陽能集熱器由雙通螺旋真空管、蓄熱棒、聯(lián)箱組成，如圖2 所示。青海地區(qū)室外溫度在-15～12 ℃，由此在集熱器內(nèi)選用空氣為傳熱介質(zhì)，避免集熱器凍裂，延長集熱器壽命；在集熱器真空管內(nèi)裝設(shè)蓄熱棒，儲存太陽能輻射熱量，避免因溫度過高造成真空管爆裂，在沒有太陽輻射時，由蓄熱棒提供熱量，緩解了太陽能利用間斷性問題；雙通螺旋結(jié)構(gòu)，避免玻璃集熱管熱脹冷縮發(fā)生爆裂，同時增加換熱面積，提高了換熱效率。

利用空氣源熱泵系統(tǒng)作為輔助能源，在太陽能集熱系統(tǒng)不能滿足烘房熱量需求時，啟動空氣源熱泵系統(tǒng)，保證烘干過程連續(xù)性、穩(wěn)定性。在枸杞烘干過程中溫度超過60 ℃[9]，枸杞易發(fā)生褐變，影響產(chǎn)品品質(zhì)，此時通過空氣源熱泵系統(tǒng)對烘房溫度進行恒溫控制，可避免枸杞褐變，提高干燥品質(zhì)。

從干燥室內(nèi)排出的熱空氣通過余熱回收系統(tǒng)進行熱量回收，減少機組的熱負荷，起到了節(jié)能減排的效果。

干燥室為長方體結(jié)構(gòu)(長寬高2 400 mm×1 200 mm×2 000 mm)如圖3所示，外壁使用兩層成型金屬面板中間發(fā)泡、高分子隔熱內(nèi)芯，起到保溫效果及強度要求?？紤]冷、熱風(fēng)的下沉與上升現(xiàn)象，干燥室采用底部全面孔板送風(fēng)，左側(cè)上部設(shè)置長寬2 400 mm×200 mm出風(fēng)口。干燥過程中物料平鋪與干燥托盤，且托盤為孔板結(jié)構(gòu)，上下氣流能夠流體均勻，忽略氣流向上的運動阻力。由此，模擬過程不考慮干燥托盤對氣流分布的影響。

2 干燥室實驗設(shè)計及模擬

2.1 正交實驗設(shè)計

正交試驗設(shè)計和分析是處理多因素優(yōu)化問題的高效處理方法，采取具有代表性的試驗點來代替整體試驗，能有迅速找到最優(yōu)方案從而降低生產(chǎn)成本[10]。

影響孔板送風(fēng)的因素有風(fēng)速、孔徑以及風(fēng)量等，根據(jù)以上因素可計算出孔隙率以及空間距[11]，用于后期三維建模。

孔隙率計算公式如下

(1)

式中K——孔隙率；

Q——送風(fēng)風(fēng)量/m3·h-1；

V——氣流流速/m·s-1；

A——送風(fēng)板面積/m2。

孔間距計算如下

(2)

式中L——孔間距/mm；

d0——選取孔徑/mm。

為研究各因素對干燥室氣流組織的影響，并找出最優(yōu)設(shè)計參數(shù)，本文采用正交表格L9(34)進行正交實驗設(shè)計，實驗因素水平如表1所示。

表1實驗因素水平表

水平實驗因素風(fēng)速/m·s-1孔徑/mm風(fēng)量/m3·h-11261 00024102 00036163 000

按照正交表格L9(34)設(shè)計，并補充一個空列因素，實驗安排如表2所示。

表2L9(34)實驗方案

因素試驗號 風(fēng)速A/m·s-1孔徑B/mm風(fēng)量Cm3·h-1空列D實驗安排11111A1B1C121222A1B2C231333A1B3C342123A2B1C252231A2B2C362312A2B3C173132A3B1C383213A3B2C193321A3B3C2

2.2 干燥室模擬

2.2.1 網(wǎng)格劃分

實驗設(shè)計分為9組，以第三組為例，設(shè)置網(wǎng)格數(shù)為56萬、104萬、224萬三個算例進行數(shù)值求解，在迭代600次左右達到收斂條件。以干燥室高度方向的溫度和速度分布為參考，56萬網(wǎng)格組溫度和流速相差較大，104萬和224萬網(wǎng)格數(shù)據(jù)無明顯差別，如圖4所示。為簡化計算，縮短模擬時長，采用104萬網(wǎng)格處理，使用結(jié)構(gòu)化六面體網(wǎng)格，對于孔板出風(fēng)口處采用局部網(wǎng)格加密處理，以確保氣流、溫度場模擬的精確度。

2.2.2 邊界條件設(shè)置

采用Fluent分離隱式求解器進行穩(wěn)態(tài)計算，斷流模型采用Realizableκ-ε(2equ)模型，能量方程開啟，由于干燥室高2 m，則開啟重力加速度影響，壓力以及速度的耦合采用Simple算法。邊界條件如下：模型的維護結(jié)構(gòu)溫度恒定邊界，入口邊界采用velocity-inlet，回風(fēng)口為自由出流outflow，壁面邊界采用第一類邊界條件，速度情況按無滑移條件處理，具體模型的邊界條件參數(shù)如表3所示。

3 正交實驗結(jié)果分析

3.1 評價指標

采用正交實驗設(shè)計時必須明確實驗?zāi)康?，確定結(jié)果的考核指標。對干燥室內(nèi)氣流均勻性的評價，使用面均勻速度系數(shù)和面均勻溫度系數(shù)作為考核指標[12]，用于恒定截面上速度與溫度的均勻程度。

表3 CFD模型邊界條件匯總

流體:空氣干燥室空氣入口溫度:334 K干燥室壁面溫度:304 K重力加速度:9.8 m/s2項目性能參數(shù)能量模型打開粘性模型κ-ε模型入風(fēng)口速度口/m·s-12/4/6排風(fēng)口速度口Outflow聚氨酯保溫墻體導(dǎo)熱系數(shù)/W·m-1·K-10.052厚度/mm100密度/kg·m-3350

表4考核指標模擬結(jié)果

指標實驗組 面均勻速度系數(shù)面均勻溫度系數(shù)10.674 20.891 520.694 30.924 530.701 50.998 840.677 90.998 350.596 00.891 260.588 50.871 270.660 60.998 680.591 20.987 890.653 00.998 1

由表4的模擬結(jié)果可得，第三組A1B3C3實驗為最佳方案，干燥室內(nèi)氣流流速以及溫度分布均勻，有效避免了因流速過低造成干燥不足，過高造成物料過度失水；同時室內(nèi)溫度基本維持在進風(fēng)溫度334 K，避免不同區(qū)域干燥不均。

3.2 實驗結(jié)果極差分析

通過正交性實驗找出最優(yōu)方案，并進行實驗結(jié)果的極差分析，初步確定各因素對速度場和溫度場的影響程度。

由表5所示，通過比較各因素的極差大小，對干燥室速度的影響程度依次為：入口風(fēng)速A，風(fēng)量C，孔徑B；對干燥室溫度的影響程度依次為：入口風(fēng)速A，風(fēng)量C，孔徑B。繪制位級趨勢圖進一步分析各因素的影響趨勢，橫坐標為因素位級，縱坐標為考核指標系數(shù)。

表5實驗指標極差分析

指標結(jié)果速度不均勻系數(shù)Ⅰ2.072.0131.854Ⅱ1.8621.8822.025Ⅲ1.9051.9431.958i0.690.6710.618ii0.6210.6270.675iii0.6350.6480.653極差R0.0690.0440.057溫度不均勻系數(shù)Ⅰ2.8152.8882.751Ⅱ2.7612.8042.921Ⅲ2.9852.8682.889i0.9380.9630.917ii0.9200.9350.974iii0.9950.9560.963極差R0.0750.0280.057

注：I、II、III分別是各列數(shù)字1、2、3相應(yīng)實驗指標之和，i、ii、iii是各列數(shù)字1、2、3相應(yīng)實驗指標的平均值，R為極差，是各列相應(yīng)位級下實驗指標平均值中最大值與最小值的差。

圖5(a)所示，干燥室內(nèi)空氣流速均勻性隨入口風(fēng)速的增加而降低，隨孔徑的增大而增加，風(fēng)量的變化對其影響相對較??；圖5(b)所示，干燥室內(nèi)空氣溫度均勻性有最佳入口風(fēng)速，隨孔徑的增大而增加，風(fēng)量變化對溫度均勻性影響較大。

3.3 實驗結(jié)果方差分析

以上通過極差分析對各因素的影響大小和影響趨勢做了簡單分析，并未考慮位級量和實驗誤差的影響，以下對實驗?zāi)M數(shù)據(jù)進行方差分析。

如表6所示，通過極差分析三個因素對速度均勻性和溫度均勻性的影響大小依次為：入口風(fēng)速A、風(fēng)量C、孔徑B，與前述極差分析結(jié)果相符，證明結(jié)果的可靠性。

表6實驗指標方差分析

指標結(jié)果速度不均勻系數(shù)Ⅰ 24.2854.0513.437Ⅱ23.4693.5404.101Ⅲ23.6283.7753.834P11.38211.36611.373Q3.7943.7893.791S0.0080.0030.005溫度不均勻系數(shù)Ⅰ 27.9238.3437.565Ⅱ27.6217.8608.532Ⅲ28.9078.2268.344P24.45224.42824.441Q8.1518.1438.147S0.0090.0010.005

4 結(jié)論

本文以實際工程需求，針對太陽能干燥室進行氣流組織的優(yōu)化研究，將數(shù)值模擬CFD與正交性實驗設(shè)計相結(jié)合，針對太陽能熱泵干燥室(2 400 mm×1 200 mm×2 200 mm)氣流不均勻性提出優(yōu)化方案。本文采用9組正交實驗減少了模擬次數(shù)，降低實驗成本，縮短實驗周期。主要分析了入口風(fēng)速、孔徑、送風(fēng)量對氣流均勻性的影響，并對模擬數(shù)據(jù)采用極差分析和方差分析確保結(jié)果的準確性。

通過實驗分析得出如下結(jié)論：

(1)流速均勻性隨入口風(fēng)速的增加而降低，隨孔徑的增大而增加，風(fēng)量的變化對其影響相對較??；

(2)溫度均勻性隨孔徑的增大而增加，風(fēng)量變化對溫度均勻性影響較大；

(3)得到最優(yōu)方案為：入口風(fēng)速2 m/s，孔徑16 mm，送風(fēng)量為3 000 m3/h，較低的風(fēng)速不易出現(xiàn)渦流，增加氣流均勻；合適的孔徑和送風(fēng)量有利于氣流擴散，避免出風(fēng)孔口之間的氣流干擾。

干燥室內(nèi)氣流組織均勻性的影響因素較多，本文分析了上述三種因素，另外孔隙率和孔板送風(fēng)的穩(wěn)壓層高度對氣流的影響也較大，還需進一步深入研究。