張楊吳哲曹軍

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

木材干燥是改善木材物理力學(xué)性能、合理使用木材、減少木材降等損失、提高木材利用率的重要技術(shù)措施，也是保障家具及其他木制品質(zhì)量的重要工序。木材干燥過程的研究對有效地利用有限的木材資源、降低能源消耗，提高木材制品質(zhì)量顯得尤為重要。在實際生產(chǎn)中，由于干燥窯的特殊結(jié)構(gòu)，要保證窯內(nèi)的氣流分布均勻極為困難[1]。因此，為了提高木材干燥效率，應(yīng)該對干燥窯的內(nèi)部流場進行分析，以改進干燥窯的結(jié)構(gòu)。

利用計算流體力學(xué)(CFD)技術(shù)對木材干燥窯進行數(shù)值模擬能夠獲得干燥窯內(nèi)的速度場、溫度場和濕度場等，可以廣泛用于優(yōu)化木材干燥工程。文獻[2]利用計算流體力學(xué)軟件FLUENT對木材干燥窯建立三維穩(wěn)態(tài)的CFD模型，進行了數(shù)值模擬計算，預(yù)測了風(fēng)機給定不同氣流速度對干燥窯的氣流分布的影響。文獻[3]利用Fluent完成了6種不同干燥窯內(nèi)壁結(jié)構(gòu)下的干燥窯內(nèi)部流場的建模與仿真，對干燥窯壁進行了改進，獲得了較為均勻的干燥窯內(nèi)部流場。文獻[4]表明，在二維結(jié)構(gòu)下，導(dǎo)流擋板既可以提高窯內(nèi)不同木材堆間隙內(nèi)的速度分布均勻性，又可以顯著消除窯內(nèi)的回流低速區(qū)。從而可以較好地改善干燥窯內(nèi)的流場。而在真實干燥窯的三維結(jié)構(gòu)內(nèi)，多臺風(fēng)機布置在干燥窯頂部，木材堆進口氣流速度并不均勻，而且三維空間內(nèi)流場結(jié)構(gòu)更為復(fù)雜。

因此為了進一步研究導(dǎo)流擋板對改善干燥窯內(nèi)流場性能的作用，本研究進行了窯內(nèi)流場的三維數(shù)值模擬，詳細對比了不同截面的速度分布，為導(dǎo)流擋板結(jié)構(gòu)的應(yīng)用奠定基礎(chǔ)。

1 幾何模型和數(shù)學(xué)模型的建立

筆者研究的是實驗室用頂風(fēng)式木材干燥窯，干燥窯結(jié)構(gòu)為方形，長 2.10 m、寬 1.40 m、高 1.50 m。干燥窯頂部安裝3臺風(fēng)機，風(fēng)機直徑R=0.44 m。采用Gambit建立干燥窯三維計算模型，用混合網(wǎng)格，總網(wǎng)格數(shù)量160萬個左右。采用k－ε湍流模型，風(fēng)機進口采用速度進口，出口給定出風(fēng)速度。圖1分別給出了干燥窯的局部結(jié)構(gòu)和網(wǎng)格劃分，其中z為干燥窯長度方向，y為高度方向，x為寬度方向。

筆者對等溫等濕條件下裝有木材堆的干燥窯內(nèi)空氣的流動情況進行研究。由于送風(fēng)速度較低，將窯內(nèi)空氣的流動視為不可壓縮的湍流，密度為常數(shù)。采用隱式求解器，速度壓力耦合求解采用SIMPLE算法，湍流模型為標準k－ε模型[5]。流動控制方程包括連續(xù)方程、動量方程以及湍流動能輸運和耗散方程，具體如下:

①連續(xù)方程:

②動量方程:

式中:u為窯內(nèi)空氣流速(m·s－1);ρ為窯內(nèi)空氣密度(kg·m－3);μ為動力黏度(Pa·s);P為空氣靜壓(Pa);g 為重力加速度(kg·m－1·s－2)。

③湍動能k和耗散率ε的輸運方程分別為:

式中:ui、uj分別為 xi、xj方向的時均速度;xi為直角坐標系的3個軸坐標;μt為湍動黏度(Pa/s);ν為運動黏度(m2·s－1);E 為時均應(yīng)變張力，c1=1.44，c2=1.92，σk=1.0，σε=1.2。

圖1 裝有弧形導(dǎo)流板的方形干燥窯的局部結(jié)構(gòu)和網(wǎng)格劃分

2 數(shù)值模擬及結(jié)果分析

本研究主要針對干燥窯內(nèi)部速度場進行詳細分析，模擬結(jié)果為三維云圖，顏色的變化代表速度大小的變化，即分布趨勢。所以通過仿真模擬出的氣流分布云圖，可以直觀地觀察出干燥窯內(nèi)的氣流分布規(guī)律及均勻性情況。為了詳細對比兩種干燥窯內(nèi)的流場特性，取定進口速度為3m/s，選取x分別為－0.40、－0.25、0.25、0.40 m，y 分別為 0.1、0.6、1.0 m，z分別為 0.10、0.45、0.75、1.05 m，共 11 個截面進行詳細分析。各截面示意圖見圖2。

圖2 截面示意圖

2.1 x截面的速度分布圖

圖3給出了x=－0.4 m截面的速度分布。圖3a和圖3b中，由于無導(dǎo)流擋板，該截面頂部出風(fēng)口附近局部高速區(qū)比較集中，安裝導(dǎo)流擋板后，由于擋板對出口氣流的分流作用，頂部高速區(qū)消失，同時，擋板之間形成的流道內(nèi)部的氣流速度較為均勻。由圖4給出的x=－0.25 m截面的速度分布可知，無導(dǎo)流擋板時，干燥窯頂部高速區(qū)較為集中;而安裝導(dǎo)流擋板后，頂部氣流被明顯分成較為均勻的3股氣流，其中對應(yīng)出風(fēng)口附近的速度略大。圖5給出了x=0.25 m截面的速度分布?？傮w上看，由于頂部風(fēng)機的抽吸作用，在干燥窯進口處出現(xiàn)相應(yīng)的3個高速區(qū)。圖6給出了x=0.4 m截面的速度分布。由圖6a可知，由于無擋板結(jié)構(gòu)，風(fēng)機的抽吸作用使風(fēng)機口附近速度較大;在正對風(fēng)機口下方，氣流速度較低，而在風(fēng)機之間的下部位置，卻形成了一個較為明顯的高速區(qū)。圖6b顯示了擋板對氣流的分流作用，擋板分開的各流道間對應(yīng)風(fēng)機口位置的速度較大，同時風(fēng)機口下方的速度也較大;而且安裝導(dǎo)流擋板后風(fēng)機之間的下方位置也無圖6a中的低速區(qū)。

圖3 x=－0.4 m截面的速度分布

圖4 x=－0.25 m截面的速度分布

圖5 x=0.25 m截面的速度分布

圖6 x=0.4 m截面的速度分布

2.2 y截面的速度分布圖

為了研究干燥窯內(nèi)不同高度處木材堆間隙內(nèi)的流場均勻性，圖7至圖9分別給出了不同高度處的速度分布。由圖7可知，對于y=0.1 m截面，無導(dǎo)流擋板時，進風(fēng)口一側(cè)形成了高速區(qū)和低速區(qū)交替分布的趨勢;安裝導(dǎo)流擋板后，進風(fēng)口一側(cè)高速區(qū)和低速區(qū)交替分布的現(xiàn)象稍微緩和一點，同時，該截面上木材間隙內(nèi)的速度顯著提高。干燥窯中部y=0.6 m截面(圖8)以及頂部y=0.1 m截面(圖9)速度分布也表明了安裝導(dǎo)流擋板后該截面上木材間隙內(nèi)速度增加，出風(fēng)口附近擋板內(nèi)的速度均勻性也得到顯著改善。圖9表明，安裝導(dǎo)流擋板結(jié)構(gòu)后，木材堆頂部間隙內(nèi)的速度分布均勻性顯著提高。圖9a表明，y=1.0 m截面上中部存在明顯的低速區(qū)，而圖9b表明，安裝導(dǎo)流擋板后，該截面中部的低速區(qū)消失，而在兩側(cè)出現(xiàn)了兩個對稱分布的低速區(qū)。

圖7 y=0.1 m截面的速度分布

圖8 y=0.6 m截面的速度分布

圖9 y=1.0 m截面的速度分布

2.3 z截面的速度分布圖

圖10至圖13給出了不同z向截面的速度分布，其中z=0.1 m為靠近干燥窯側(cè)壁的截面，z=0.45 m為經(jīng)過靠近壁面的風(fēng)機軸線截面，z=0.75 m為兩風(fēng)機間的中間截面，而z=1.05 m為整個干燥窯z向?qū)ΨQ面，該截面過中間風(fēng)機的軸線。圖10至圖13中，風(fēng)機送風(fēng)口在右側(cè)上方，出風(fēng)口在左側(cè)上方。由圖10可知，擋板結(jié)構(gòu)顯著提高了z=0.1 m截面的速度均勻性，消除了出風(fēng)口附近的高速區(qū)。而在z=0.45 m和z=0.75 m截面上，擋板起到了較好的導(dǎo)流作用，同時消除了干燥窯側(cè)壁附近的高速區(qū)。從圖13給出的對稱面的速度分布可知，在方形干燥窯內(nèi)，風(fēng)機吹動形成的氣流沖擊著干燥窯的側(cè)壁(圖中右側(cè))，從而在側(cè)壁附近形成高速區(qū)，而在木材堆進口形成了低速區(qū)。安裝導(dǎo)流擋板后，方形干燥窯木材堆最上方速度有所降低，這主要是由干燥窯頂部轉(zhuǎn)折的結(jié)構(gòu)所致;而在安裝導(dǎo)流擋板的圓弧形干燥窯內(nèi)，整個流場速度分布較為合理，沒有出現(xiàn)低速區(qū)或者回流區(qū)。

圖10 z=0.1 m截面的速度分布

圖11 z=0.45 m截面的速度分布

圖12 z=0.75 m截面的速度分布

圖13 z=1.05 m截面的速度分布

3 結(jié)束語

用GAMBIT建立了實驗室用頂風(fēng)式木材干燥窯的三維幾何模型，并對該模型進行了網(wǎng)格劃分，然后利用FLUENT進行了干燥窯內(nèi)空氣流動的仿真計算。選取了 x分別為－0.40、－0.25、0.25、0.40 m，y分別為0.1、0.6、1.0 m，z分別為 0.10、0.45、0.75、1.05 m，共11個截面詳細分析了干燥窯內(nèi)的速度分布。結(jié)果顯示:無導(dǎo)流擋板時，送風(fēng)口附近速度較高，木材堆間隙內(nèi)沿橫向的速度均勻性較差，而安裝導(dǎo)流擋板后，從木材堆底部至頂部所有間隙內(nèi)沿橫向的速度均勻性顯著提高。由于擋板的導(dǎo)流作用，裝有弧形導(dǎo)流板的干燥窯內(nèi)流場比較均勻。分析證明該方法可以省時省力地模擬木材干燥窯內(nèi)任何時間任何空間內(nèi)的流場，預(yù)測風(fēng)速，進而提高木材干燥的質(zhì)量。

[1]朱政賢.木材干燥[M].北京:中國林業(yè)出版社，1992:102－103.

[2]張瑞雪，孫麗萍.計算流體力學(xué)在干燥窯風(fēng)速檢測中的應(yīng)用研究[J].機電產(chǎn)品開發(fā)與創(chuàng)新，2010，23(1):119－121.

[3]茹煜，賈志成，郁金.基于Fluent軟件的木材干燥窯內(nèi)部流場分析研究[J].木材加工機械，2010(4):12－15.

[4]張楊，吳哲，曹軍.木材干燥窯結(jié)構(gòu)對窯內(nèi)空氣流動特性的影響[J].東北林業(yè)大學(xué)學(xué)報，2012，40(9):116－119.

[5]王福軍.計算流體力學(xué)分析:CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社，2004.