孟楊　陳廣元

(東北林業(yè)大學(xué)，哈爾濱，150040)

?

木材干燥室強(qiáng)制進(jìn)氣孔口的氣流均勻性優(yōu)化1)

孟楊陳廣元

(東北林業(yè)大學(xué)，哈爾濱，150040)

摘要利用計(jì)算流體力學(xué)數(shù)值模擬、正交試驗(yàn)極差和方差分析等方法，仿真模擬木材干燥室進(jìn)氣道孔口位置、孔口直徑、進(jìn)氣風(fēng)機(jī)全壓與動(dòng)壓組合3個(gè)因素對(duì)進(jìn)氣孔口氣流均勻性的影響規(guī)律。結(jié)果表明：3個(gè)因素對(duì)進(jìn)氣孔口氣流速度、風(fēng)壓均勻性的影響程度從大到小依次為孔口直徑、孔口位置、進(jìn)氣風(fēng)機(jī)風(fēng)壓。在35～50 mm范圍內(nèi)，孔口直徑越小，氣流均勻性越好；孔口位置在靠近豎直平面中心偏上時(shí)，氣流均勻性更好；進(jìn)氣風(fēng)機(jī)全壓與動(dòng)壓組合對(duì)氣流均勻性的影響程度較低。根據(jù)分析結(jié)果和實(shí)際生產(chǎn)現(xiàn)狀，綜合得出氣流均勻性最優(yōu)時(shí)孔口直徑40 mm，孔口位置距豎直平面上端2/5處，進(jìn)氣風(fēng)機(jī)全壓690 Pa、動(dòng)壓90 Pa。按照這些參數(shù)設(shè)計(jì)的進(jìn)氣孔口氣流速度不均勻性系數(shù)0.035，最大速度偏差9%，平均氣流速度15.7 m/s，孔口氣流均勻性可滿足生產(chǎn)實(shí)際要求。

關(guān)鍵詞木材干燥室；進(jìn)氣孔口；氣流均勻性；數(shù)值模擬

木材干燥是保障木制品加工利用質(zhì)量的重要工序，木材干燥質(zhì)量的優(yōu)劣直接決定木材的利用率和產(chǎn)品的價(jià)值[1-3]。由于木材常規(guī)干燥方法比較靈活，易于操作，能夠滿足絕大多數(shù)木材干燥質(zhì)量的要求，在今后的木材加工和家具制造企業(yè)的干燥生產(chǎn)中，木材常規(guī)干燥仍將發(fā)揮主導(dǎo)作用[4]。木材干燥室強(qiáng)制集中換氣系統(tǒng)是改進(jìn)干燥生產(chǎn)中室內(nèi)外干濕空氣交換過程，有效降低生產(chǎn)能耗的一項(xiàng)新裝置。國(guó)內(nèi)現(xiàn)有常規(guī)木材干燥過程的干濕氣體交換是通過干燥室頂部或側(cè)壁上部的進(jìn)排氣道，利用干燥介質(zhì)循環(huán)風(fēng)機(jī)運(yùn)行中產(chǎn)生的正負(fù)壓直接將熱濕氣體排出室外，同時(shí)將新鮮干空氣吸入室內(nèi)的。在干濕空氣交換過程中，部分經(jīng)加熱器加熱后的熱空氣未經(jīng)過材堆直接排出室外，排出的熱濕空氣中熱能未被回收利用，造成熱能浪費(fèi)。相關(guān)研究表明，干燥室排出的熱濕空氣所含的熱能約占整個(gè)干燥過程中全部能耗的56%以上[5]。為降低干燥生產(chǎn)熱能消耗，陳廣元等設(shè)計(jì)了一種木材干燥裝置強(qiáng)制集中進(jìn)排氣系統(tǒng)[6]，已在企業(yè)生產(chǎn)中應(yīng)用。在木材常規(guī)干燥生產(chǎn)中，干燥室內(nèi)干燥介質(zhì)循環(huán)氣流均勻性對(duì)材堆干燥均勻度有重要影響[7-9]。強(qiáng)制集中換氣系統(tǒng)進(jìn)排氣道孔口氣流均勻性會(huì)直接影響沿材堆長(zhǎng)度方向干燥介質(zhì)循環(huán)均勻性，進(jìn)而影響材堆干燥均勻度。為進(jìn)一步優(yōu)化完善強(qiáng)制集中換氣系統(tǒng)的設(shè)計(jì)方法和理論，楊雪等研究了進(jìn)排氣孔口半徑對(duì)進(jìn)排氣均勻性的影響，得出了計(jì)算孔口半徑回歸方程[10]。筆者把進(jìn)氣道孔口作為研究對(duì)象，利用計(jì)算流體力學(xué)數(shù)值模擬、正交試驗(yàn)極差和方差分析等方法，探討進(jìn)氣道孔口面積、孔口位置、進(jìn)氣風(fēng)機(jī)參數(shù)等因素對(duì)孔口氣流均勻性的影響規(guī)律，優(yōu)化設(shè)計(jì)參數(shù)，為強(qiáng)制集中換氣系統(tǒng)在干燥生產(chǎn)實(shí)踐中的應(yīng)用提供一定的理論依據(jù)。

1實(shí)驗(yàn)設(shè)計(jì)

1.1干燥室結(jié)構(gòu)

如圖1所示，該干燥室由強(qiáng)制集中進(jìn)排氣系統(tǒng)及常規(guī)干燥室主體兩部分構(gòu)成。排潮時(shí)，架設(shè)在干燥室外的進(jìn)排氣風(fēng)機(jī)強(qiáng)制將室外干冷空氣輸送入室內(nèi)，同時(shí)將室內(nèi)熱濕空氣排出室外。這個(gè)過程中，利用廢熱回收裝置進(jìn)行余熱回收，使排出的熱濕空氣溫度降低，同時(shí)進(jìn)入的干冷空氣溫度升高。干燥室內(nèi)部尺寸為4 600mm(深)×8 600mm(寬)×6 000mm(高)，實(shí)際裝材尺寸4 000mm(深)×6 600mm(寬)×4 200mm(高)，實(shí)際裝材容積約為60m3。該干燥室已在生產(chǎn)中應(yīng)用，圖1中A、B、C為應(yīng)用現(xiàn)場(chǎng)各部分的實(shí)例照片。研究對(duì)象為強(qiáng)制集中進(jìn)排氣系統(tǒng)的進(jìn)氣道部分，如圖1B所示，研究重點(diǎn)為進(jìn)氣道孔口的氣流均勻性。

A.廢熱回收裝置實(shí)例圖；B.進(jìn)排氣道實(shí)例圖；C.材堆實(shí)例圖；1.排氣道；2.進(jìn)氣道；3.加熱器；4.軸流循環(huán)風(fēng)機(jī)；5.擋風(fēng)板；6.大門；7.材堆。

圖2為干燥室進(jìn)氣道及孔口模型。利用CAD軟件進(jìn)行三維建模，孔口編號(hào)時(shí)都從靠近風(fēng)機(jī)的一端開始，模型大小和現(xiàn)有生產(chǎn)應(yīng)用中設(shè)備保持一致。進(jìn)氣道截面為三角形，兩直角邊分別為400、320mm，進(jìn)氣道結(jié)構(gòu)總長(zhǎng)度4 600mm；輸氣管路截面為正方形，長(zhǎng)度和寬度均為200mm，輸氣管路總長(zhǎng)度8 600mm；與進(jìn)氣風(fēng)機(jī)直接連接的輸氣管路截面為總進(jìn)氣口，進(jìn)氣道截面面積S=0.064m2。

圖2　干燥室進(jìn)氣道及孔口模型

1.2進(jìn)氣道孔口直徑、數(shù)量的理論計(jì)算

孔口所在平面高度與孔徑之比大于10時(shí)的孔口稱為小孔口，孔口斷面上各點(diǎn)的參數(shù)可以看作是常數(shù)，即忽略孔口處流體勢(shì)能的差別[11]。根據(jù)氣流原理，將孔口設(shè)計(jì)為小孔口，計(jì)算得到孔口最大允許直徑為50mm。再根據(jù)均勻氣流理論，要求孔口面積總和小于進(jìn)氣道截面面積，計(jì)算得到孔徑50mm時(shí)，允許的孔口最大數(shù)量為30個(gè)，此時(shí)孔口總面積占截面面積的91.99%。

1.3正交試驗(yàn)

利用正交試驗(yàn)法對(duì)多因素影響孔口氣流均勻性進(jìn)行研究。影響孔口氣流的主要因素有8個(gè)，由孔口直徑、數(shù)量的理論計(jì)算及生產(chǎn)應(yīng)用中設(shè)備的現(xiàn)狀，確定了孔口直徑、孔口位置、進(jìn)氣風(fēng)機(jī)全壓與動(dòng)壓組合3個(gè)因素(表1)。用孔口距進(jìn)氣道頂部所處的位置表示孔口位置，進(jìn)氣風(fēng)機(jī)提供的恒定風(fēng)壓表示進(jìn)氣全壓與動(dòng)壓組合。風(fēng)壓及氣流速度由標(biāo)準(zhǔn)畢托管風(fēng)速計(jì)(DP1000系列智能壓力風(fēng)速儀，風(fēng)壓測(cè)量范圍0～2 000Pa，氣流速度測(cè)量范圍0～57m/s)測(cè)得，分別對(duì)應(yīng)進(jìn)氣風(fēng)機(jī)用電頻率為50、40、30、20Hz運(yùn)行時(shí)的狀態(tài)。每個(gè)因素4個(gè)水平，其他因素都設(shè)為定量。孔口為圓形，孔口數(shù)量為30個(gè)，管道流向?yàn)檎?，進(jìn)氣道截面為三角形，進(jìn)氣道長(zhǎng)度4 600mm。采用L16(45)標(biāo)準(zhǔn)正交試驗(yàn)表[12]，3個(gè)因素占3個(gè)標(biāo)準(zhǔn)列，另設(shè)2個(gè)空白列。

表1干燥室強(qiáng)制進(jìn)氣道孔口氣流均勻性正交試驗(yàn)的因素與水平

水平孔口直徑(A)/mm孔口位置(B)進(jìn)氣風(fēng)機(jī)壓力(C)/Pa全壓動(dòng)壓1501/5處9601602452/5處690903403/5處360604354/5處16030

2數(shù)值模擬

2.1數(shù)值模擬原理

2.1.1基本方程

流體流動(dòng)受物理守恒定律的支配，包括質(zhì)量守恒定律、動(dòng)量守恒定律、能量守恒定律，這些定律直觀反映為控制方程。進(jìn)氣系統(tǒng)通過恒定動(dòng)力源供氣，可簡(jiǎn)化為三維不可壓流體的穩(wěn)態(tài)流動(dòng)，密度不隨時(shí)間變化；為方便研究，忽略流動(dòng)過程中的熱交換和重力影響?？刂品匠倘缡?1)、(2)所示。

連續(xù)方程：

(1)

u-動(dòng)量方程：

(2)

說明：將u-動(dòng)量方程中的u和x，替換成v和y得到v-動(dòng)量方程，替換成w和z得到w-動(dòng)量方程。式(1)和式(2)中：u、v、w為速度分量(m/s)；ρ為干空氣的密度，取1.20kg/m3；ν為流體的運(yùn)動(dòng)黏度，取1.57mm2/s；P為空氣靜壓(Pa)；U為速度矢量。

2.1.2湍流方程

湍流模型采用標(biāo)準(zhǔn)的k-ε[13]，可用于平面流動(dòng)、管道流動(dòng)、孔口氣流等方面研究，仿真度高。該模型包含湍動(dòng)能(k)和湍動(dòng)擴(kuò)散率(ε)兩個(gè)未知量，其中：

(3)。

式中：μ為流體的動(dòng)力黏度(Pa·s)；ρ為空氣密度(kg/m3)。

帶有經(jīng)驗(yàn)常數(shù)的湍動(dòng)黏度(μt)公式為

(4)

式中：ρ為空氣密度(kg/m3)；Cμ為經(jīng)驗(yàn)常數(shù)。

根據(jù)研究條件，三維穩(wěn)態(tài)不可壓流體的標(biāo)準(zhǔn)k-ε湍流控制方程為

(5)

(6)

式中：Gk為平均速度梯度引起的湍動(dòng)能的產(chǎn)生項(xiàng)；C1ε、C2ε為經(jīng)驗(yàn)常數(shù)；σε和σk分別為與湍動(dòng)能和耗散率對(duì)應(yīng)的Prandtl數(shù)。

2.2網(wǎng)格劃分

利用Gambit對(duì)進(jìn)氣系統(tǒng)模型進(jìn)行網(wǎng)格劃分，選擇Tet/Hybrid型三維不規(guī)則網(wǎng)格，網(wǎng)格總數(shù)96萬左右，部分三維網(wǎng)格如圖3所示。

2.3邊界條件

進(jìn)氣道入口是與進(jìn)氣風(fēng)機(jī)直接連接的方管，邊長(zhǎng)為200mm的正方形。進(jìn)氣道入口邊界條件的模擬參數(shù)為進(jìn)氣風(fēng)機(jī)出口的全壓與靜壓、雷諾數(shù)、湍流強(qiáng)度、湍流直徑；進(jìn)氣風(fēng)機(jī)出口為30個(gè)與干燥室內(nèi)部相通的孔口，采用壓力出口邊界條件模擬，保持與常壓環(huán)境一致，壓強(qiáng)為一個(gè)標(biāo)準(zhǔn)大氣壓。根據(jù)流體動(dòng)力學(xué)相關(guān)公式[14]，計(jì)算每組實(shí)驗(yàn)的入口邊界條件(表2)。

圖3　干燥室進(jìn)氣道網(wǎng)格劃分

試驗(yàn)組合入口邊界條件全壓/Pa靜壓/Pa雷諾數(shù)/10-5湍流強(qiáng)度/%湍流直徑/mA1B1C19608002.013.470.20A1B2C26906001.533.600.20A1B3C33603001.273.680.20A1B4C41601300.903.850.20A2B1C26906001.533.600.20A2B2C19608002.013.470.20A2B3C41601300.903.850.20A2B4C33603001.273.680.20A3B1C33603001.273.680.20A3B2C41601300.903.850.20A3B3C19608002.013.470.20A3B4C26906001.533.600.20A4B1C41601300.903.850.20A4B2C33603001.273.680.20A4B3C26906001.533.600.20

2.4模擬計(jì)算

利用CFD仿真軟件進(jìn)行數(shù)值模擬，選用分離式求解器，湍流模型選用標(biāo)準(zhǔn)k-ε，隱式求解，離散殘差收斂標(biāo)準(zhǔn)10-4，通過模擬計(jì)算得到不同實(shí)驗(yàn)組合的孔口氣流速度分布。

2.5模擬驗(yàn)證

為了保證模擬結(jié)果的真實(shí)性，將數(shù)值模擬得到的孔口氣流速度分布與生產(chǎn)實(shí)際中的換氣裝置進(jìn)氣孔口測(cè)量值進(jìn)行比較。模型大小與實(shí)際設(shè)備保持一致，此時(shí)孔口直徑50mm，孔口數(shù)量30個(gè)，孔口位置距進(jìn)氣道頂部1/5處均勻排列，風(fēng)機(jī)全壓960Pa、動(dòng)壓160Pa。風(fēng)壓及氣流速度取多次測(cè)量結(jié)果的平均值(見圖4)?？芍?，進(jìn)氣道靠近進(jìn)氣管路一側(cè)的幾個(gè)孔口，較遠(yuǎn)端孔口風(fēng)速有很大程度減小，風(fēng)速在最低點(diǎn)孔口4兩側(cè)對(duì)稱排列。這可能因?yàn)闅饬饔蛇M(jìn)氣管路進(jìn)入進(jìn)氣道的過程中，在兩個(gè)管路連接處發(fā)生擴(kuò)散流動(dòng)，且整體流線方向發(fā)生90°轉(zhuǎn)換。所以靠近進(jìn)氣管路孔口的氣流變化劇烈，整體流線首先要使流動(dòng)變得穩(wěn)定，然后再由孔口流出。對(duì)比模擬和測(cè)試結(jié)果，氣流速度整體趨勢(shì)一致，最大差值0.7m/s，均在允許的范圍內(nèi)，說明模擬仿真度高，驗(yàn)證了計(jì)算流體力學(xué)方法研究孔口氣流狀態(tài)的準(zhǔn)確性。

圖4　干燥室進(jìn)氣孔口氣流速度模擬值與實(shí)測(cè)值對(duì)比

3結(jié)果與分析

3.1不同實(shí)驗(yàn)組合的氣流速度分布

由圖5可知，不同實(shí)驗(yàn)組合對(duì)氣流速度的分布都有不同程度的影響，有的組合從直觀來看波動(dòng)比

較大，表現(xiàn)為不均勻；有的波動(dòng)比較小，表現(xiàn)為均勻。為探討多組數(shù)據(jù)之間的均勻性關(guān)系，消除尺寸和量綱的影響，引入氣流速度不均勻系數(shù)(δ)對(duì)模擬結(jié)果比較分析，計(jì)算公式如式(7)所示。

(7)

3.2模擬結(jié)果的極差分析

表3是不同實(shí)驗(yàn)組合進(jìn)氣孔口的氣流速度不均勻系數(shù)。從模擬結(jié)果來看，實(shí)驗(yàn)組合A3B2C4的不均勻系數(shù)最小，為部分實(shí)驗(yàn)的最優(yōu)組合。對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行極差分析(見表4)，可知對(duì)進(jìn)氣孔氣流速度不均勻系數(shù)影響的順序從大到小依次為孔口直徑(A)、孔口位置(B)、全壓與動(dòng)壓組合(C)。

圖5　不同實(shí)驗(yàn)組合的干燥室進(jìn)氣孔氣流速度分布

孔口直徑直徑/mmδ孔口位置位置δ進(jìn)氣風(fēng)機(jī)全壓/Pa動(dòng)壓/Paδ500.1271/5處0.0769601600.073450.0662/5處0.057690900.081400.0563/5處0.080360600.074350.0544/5處0.090160300.075

正交試驗(yàn)設(shè)計(jì)力求以少數(shù)實(shí)驗(yàn)估計(jì)整體，最優(yōu)結(jié)果可能在沒有進(jìn)行實(shí)驗(yàn)的組合中出現(xiàn)。為找到可能存在的最優(yōu)結(jié)果，對(duì)模擬結(jié)果進(jìn)行位級(jí)趨勢(shì)分析?？芍卓谥睆皆叫?，氣流速度不均勻性系數(shù)越?。豢卓谖恢迷娇拷Q直平面中心偏上時(shí)，氣流速度不均勻性系數(shù)越小。全壓與動(dòng)壓組合趨勢(shì)變化范圍較小，對(duì)氣流均勻性影響較低。分析得到理論最優(yōu)組合為A4B2C1。綜合考慮平均風(fēng)速和進(jìn)氣質(zhì)量流量總量的關(guān)系，發(fā)現(xiàn)A3和A4影響近似；考慮到孔口面積變大，可使平均氣流變小，風(fēng)阻變小，進(jìn)氣質(zhì)量流量總量變大，選用孔口直徑40mm最合理；全壓與動(dòng)壓組合對(duì)氣流均勻性的影響不大，選取適中的全壓與動(dòng)壓組合C2，既能滿足均勻度要求，又有一定進(jìn)氣動(dòng)力。此時(shí)，組合A3B2C2為最終優(yōu)化組合，即孔口直徑40mm，孔口位置2/5處，全壓690Pa、動(dòng)壓90Pa。對(duì)該組合進(jìn)行數(shù)值模擬，得到氣流速度不均勻性系數(shù)為0.035，最大速度偏差9%(當(dāng)最大偏差小于10%，可認(rèn)為氣流均勻)[15]，平均氣流速度為15.7m/s，孔口氣流均勻性好。

表4　干燥室進(jìn)氣孔氣流速度不均勻系數(shù)的極差分析

3.3模擬結(jié)果的方差分析

為分析極差結(jié)果的客觀性及各因素顯著性，采用P、Q、S的方法[16]對(duì)模擬結(jié)果進(jìn)行方差計(jì)算(見表5)。結(jié)果表明,對(duì)進(jìn)氣孔氣流速度不均勻系數(shù)影響的順序從大到小依次為孔口直徑(A)、孔口位置(B)、全壓與動(dòng)壓組合(C)，與極差分析結(jié)果一致，說明極差分析結(jié)果是正確的。

表5　干燥室進(jìn)氣孔氣流速度不均勻系數(shù)的方差分析

3.4模擬結(jié)果的顯著性分析

對(duì)模擬結(jié)果進(jìn)行顯著性分析。根據(jù)數(shù)理統(tǒng)計(jì)原理，p值為顯著性概率。當(dāng)p<0.05時(shí)，認(rèn)為該因素對(duì)實(shí)驗(yàn)結(jié)果有顯著影響。由表6可知，孔口直徑(A)是影響進(jìn)氣孔口氣流均勻性的重要因素，孔口位置(B)為次要因素，全壓與動(dòng)壓組合(C)相對(duì)于孔口直徑和孔口位置，影響不顯著。這與極差分析得到的優(yōu)化結(jié)果一致，說明本正交試驗(yàn)所取的因素和水平是合理的，誤差在一定的水平內(nèi)也具有合理性。

表6　干燥室進(jìn)氣孔氣流速度不均勻系數(shù)的顯著性分析

4結(jié)論

對(duì)木材干燥室強(qiáng)制集中換氣系統(tǒng)進(jìn)氣孔口氣流均勻性的研究結(jié)果表明：各因素對(duì)進(jìn)氣孔口氣流均勻性的影響順序從大到小依次為孔口直徑、孔口位置、進(jìn)氣風(fēng)機(jī)全壓與動(dòng)壓組合。其中孔口直徑為主要影響因素，在35～50mm范圍內(nèi)孔口直徑越小，進(jìn)氣均勻性越好；孔口位置為次要影響因素，孔口位置越靠近豎直平面中心偏上時(shí)，進(jìn)氣均勻性越好；進(jìn)氣風(fēng)機(jī)全壓與動(dòng)壓組合對(duì)進(jìn)氣均勻性的影響較低。結(jié)合實(shí)際生產(chǎn)現(xiàn)狀，綜合得到最優(yōu)設(shè)計(jì)參數(shù)應(yīng)為孔口直徑40mm，孔口位置2/5處，進(jìn)氣風(fēng)機(jī)全壓690Pa、動(dòng)壓90Pa。此時(shí)進(jìn)氣孔口氣流速度不均勻性系數(shù)0.035，最大速度偏差9%，平均氣流速度15.7m/s，孔口氣流均勻性滿足生產(chǎn)實(shí)際要求，可供強(qiáng)制集中換氣系統(tǒng)在干燥生產(chǎn)實(shí)踐中應(yīng)用參考。

參考文獻(xiàn)

[1]張璧光.我國(guó)木材干燥技術(shù)現(xiàn)狀與國(guó)內(nèi)外發(fā)展趨勢(shì)[J].北京林業(yè)大學(xué)學(xué)報(bào)，2002，24(5/6)：266-270.

[2]顧煉百.木材干燥理論在木材加工技術(shù)中的應(yīng)用分析[J].南京林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版)，2008，32(5)：27-31.

[3]李磊，陳廣元.改善頂風(fēng)機(jī)型木材干燥室空氣流動(dòng)特性研究[J].森林工程，2010，26(3)：26-28.

[4]艾沐野，李莉，閆佳玉，等.對(duì)改進(jìn)和完善常規(guī)木材干燥室設(shè)施的探討[J].林業(yè)機(jī)械與木工設(shè)備，2015，43(2)：4-7.

[5]馮立寧，李建榮，陳廣元，等.常規(guī)干燥室余熱回收裝置的研究[J].森林工程，2012，28(2)：22-25.

[6]陳廣元，楊雪，陳昊，等.木材干燥裝置強(qiáng)制集中進(jìn)排氣系統(tǒng)：201010614287.9[P].2012-11-14.

[7]朱政賢，艾牧野，張曉峰，等.鋸材干燥設(shè)備性能檢測(cè)方法：GB/T17661-1999[S].北京：中國(guó)標(biāo)準(zhǔn)出版社，1999.

[8]李磊，陳廣元，李明君.多級(jí)弧形板導(dǎo)流裝置及其對(duì)木材干燥室介質(zhì)循環(huán)均勻性的改善[J].東北林業(yè)大學(xué)學(xué)報(bào)，2010，38(5)：128-130.

[9]ZADINV,KASEMGIH,VALDNAV,etal.Applicationofmultiphysicsandmultiscalesimulationstooptimizeindustrialwooddryingkilns[J].AppliedMathematicsandComputation,2015,267:465-475.

[10]楊雪，陳廣元，馮立寧，等.木材干燥室集中換氣時(shí)進(jìn)排氣孔氣流的均勻性模擬[J].北京林業(yè)大學(xué)學(xué)，2011，33(4)：113-117.

[11]陳卓如.工程流體力學(xué)[M].2版.北京：高等教育出版社，2004：282.

[12]任露泉.實(shí)驗(yàn)優(yōu)化設(shè)計(jì)與分析[M].2版.北京：高等教育出版社，2003：541.

[13]VERSTEEGH.Anintroductiontocomputationalfluiddynamics:Thefinitevolumemethod[M]. 2nded.Beijing：WorldPublishingCompany，2006：72-85.

[14]王福軍.計(jì)算流體動(dòng)力學(xué)分析：CFD軟件原理與應(yīng)用[M].北京：清華大學(xué)出版社，2008：215.

[15]楊雪.木材常規(guī)干燥設(shè)備進(jìn)排氣系統(tǒng)節(jié)能研究[D].哈爾濱：東北林業(yè)大學(xué)，2011：9-13.

[16]金莎，沈愜，孫大明，等.基于CFD方法的孔板送風(fēng)氣流組織優(yōu)化研究[J].低溫工程，2015(1)：23-28.

第一作者簡(jiǎn)介：孟楊，男，1990年9月生，東北林業(yè)大學(xué)材料科學(xué)與工程學(xué)院，碩士研究生。E-mail：lxf_mengyang@126.com。 通信作者：陳廣元，東北林業(yè)大學(xué)材料科學(xué)與工程學(xué)院，教授。E-mail：cgy32515@163.com。

收稿日期：2015年12月30日。

分類號(hào)S715.3

OptimizationofInletAirflowUniformityinWoodDryingKilnVentilationSystembyNumericalSimulation//

MengYang,ChenGuangyuan

(NortheastForestryUniversity,Harbin150040,P.R.China)//JournalofNortheastForestryUniversity，2016,44(8)：87-91.

Inordertomaketherunwithoutahitch,theComputationalFluidDynamic(CFD)andtheOrthogonalExperimentDesignwereutilizedtomodeltheeffectofthreeparameters.Therangeanalysisandthevarianceanalysiswerealsousedtodiscovertheobjectivelawsaboutvelocityuniformity.Theimpactfactorsrankingwasorificediameter,orificelocationandthecombinationoftotalanddynamicpressure.Orificediameterwasthemaininfluencefactorshowedinthefurtherstudy.Thesmallerorificediameterwas,themoreuniformlyinletflowed.Orificelocationwasthesecondaryinfluencefactor.Whenthelocationwasabovetheverticalplaneofthecenter,inletfloweduniformly.Fortheparametersoftheinletfan,ithadlittleimpacttouniformity.Bytheactualproductionsituationanalysis,theoptimalcombinationwasobtained,thatis,orificediameterwas40mm,locationwas2/5,andtheinlettotalpressureanddynamicpressureofthefanwere690and90Pa,respectively.Forthiscase,theorificeflowvelocityasymmetrycoefficientwas0.035,themaximumspeeddeviationwas9%,andaveragespeedwas15.7m/s.

KeywordsWood drying kiln； Inlet airflow； Air-flow uniformity； Numerical simulation

1)林業(yè)公益性行業(yè)科研專項(xiàng)(201304502)；林業(yè)科學(xué)技術(shù)推廣項(xiàng)目(2014-30)。

責(zé)任編輯：戴芳天。