氣墊式熱風(fēng)箱及紙張傳熱性能模擬研究

張長飛　丁克強(qiáng)　李乾軍　張東平　吳功德

(南京工程學(xué)院環(huán)境工程系,江蘇南京,211167)

?

·氣墊式熱風(fēng)箱·

氣墊式熱風(fēng)箱及紙張傳熱性能模擬研究

張長飛丁克強(qiáng)李乾軍張東平吳功德

(南京工程學(xué)院環(huán)境工程系,江蘇南京,211167)

通過最小二乘法研究熱空氣在氣墊式熱風(fēng)箱內(nèi)的速度、靜壓、溫度分布,探討了入口風(fēng)速、高度比、入口溫度、開孔率和噴口方式對氣墊層不均勻指數(shù)、氣墊靜壓和表面?zhèn)鳠嵯禂?shù)的影響規(guī)律。通過量綱分析擬合Nu的函數(shù)關(guān)系式,數(shù)值模擬值與擬合值比較表明,最大相對誤差為9.53%,均方根誤差為9.5%,擬合方程可為熱風(fēng)箱氣墊層表面?zhèn)鳠嵯禂?shù)的優(yōu)化研究提供參考。

氣墊式；熱風(fēng)箱；傳熱特性；數(shù)值模擬

(*E-mail: dndx1984@126.com)

熱風(fēng)箱是形成空氣沖擊干燥的關(guān)鍵設(shè)備。空氣沖擊干燥是將高溫空氣通過噴嘴以極高的速度噴到紙幅的表面,高速氣流對紙幅表面產(chǎn)生極大的沖擊作用,消除了紙幅表面的滯流層進(jìn)而攜走水分的過程[1]。由于沖擊流產(chǎn)生大的湍流,因而具有較高的傳熱、傳質(zhì)系數(shù),可以有效地縮短干燥時間,具有極高的傳熱效率和蒸發(fā)速度。

研究表明，沖擊流的傳熱、傳質(zhì)系數(shù),隨流體相對噴嘴的位置而變化；與噴嘴的結(jié)構(gòu)、噴嘴的排列方式、流體出口截面距沖擊面的距離、介質(zhì)的溫度、被沖擊表面的幾何形狀、排氣孔的設(shè)置等因素有關(guān)；噴嘴的形式也嚴(yán)重影響傳遞過程,要達(dá)到相同的傳熱系數(shù),條形噴嘴的沖擊速度要比圓形噴嘴的低40%左右[2]。

鄭榮和[3]探討了氣墊式漿板機(jī)熱風(fēng)箱系統(tǒng)的氣體流動,首次建立氣墊式漿板機(jī)熱風(fēng)箱氣體流動實驗臺和熱風(fēng)箱氣體流動的數(shù)學(xué)模型。在不同的入口風(fēng)速和外形尺寸的工況下對熱風(fēng)箱內(nèi)的流場進(jìn)行了實驗研究,測量了熱風(fēng)箱內(nèi)的速度和靜壓分布,探討了入口風(fēng)速和外形尺寸對熱風(fēng)箱內(nèi)靜壓和速度的影響,分析了熱風(fēng)箱內(nèi)靜壓分布及氣墊層的均勻性。曾滿連[4]在此基礎(chǔ)上進(jìn)一步對氣墊式漿板機(jī)熱風(fēng)箱的流場、壓力場和傳熱性能進(jìn)行了研究,分析了入口風(fēng)速、溫度、開孔率等因素對氣墊式漿板機(jī)熱風(fēng)箱系統(tǒng)氣體流動特性和傳熱特性的影響。結(jié)果表明，傳熱系數(shù)與噴口傾角成反比,入口風(fēng)速、出風(fēng)板開孔率、高度比對傳熱系數(shù)有一定影響,均呈二次函數(shù)關(guān)系。

目前,對氣墊干燥過程的研究較少,其他關(guān)于熱風(fēng)箱系統(tǒng)氣體流動特性和傳熱特性的研究寥寥無幾,沒有考慮到其他因素對氣體流動和傳熱特性的影響,如熱風(fēng)箱噴口形狀和布置方式等。熱風(fēng)箱噴口處氣墊層流場規(guī)律及傳熱特性非常復(fù)雜,不同的噴口形式和布置方式、入風(fēng)參數(shù)對流場規(guī)律和傳熱特性的影響很大[5]。目前，關(guān)于熱風(fēng)箱開孔形狀和布置方式對系統(tǒng)氣體流動特性和傳熱特性的影響未見報道。

本文將深入研究氣墊式熱風(fēng)箱在不同開孔形式和布置方式下風(fēng)速、開孔率、溫度等參數(shù)對氣體流動特性和傳熱特性的影響。

氣墊層傳熱系數(shù)的大小直接影響熱風(fēng)箱的干燥速率[6]。氣墊式熱風(fēng)箱干燥原理如圖1所示,熱空氣通過下吹箱上的噴孔將紙幅托起,使其懸浮在熱風(fēng)箱上,減少了紙幅在干燥器內(nèi)的摩擦阻力,紙幅斷紙、破損的機(jī)率大大減少,同時,熱空氣與紙幅間進(jìn)行強(qiáng)烈的傳熱傳質(zhì)過程,使紙幅上水分蒸發(fā)、擴(kuò)散至空氣中,紙幅不斷得以干燥。上吹箱把熱空氣吹向紙幅上部,進(jìn)一步增強(qiáng)紙幅水分蒸發(fā)效果[7-9]。上、下吹箱開孔方式不一樣,開孔作用也不同。上吹箱面上的開孔方式是圓形,而下吹箱面上的開孔方式是中間三角形,兩側(cè)圓形。熱空氣的壓力在熱風(fēng)箱中逐漸降低,為了使吹到紙幅上的風(fēng)力均勻,熱風(fēng)箱常設(shè)計成開口端向封閉端逐漸縮小的錐形,見圖2。

圖1　氣墊式熱風(fēng)箱原理圖

圖2　熱風(fēng)箱流場示意圖

根據(jù)相似原理,建立熱風(fēng)箱干燥實驗臺。相似原理的3個核心內(nèi)容：幾何相似、運動相似和動力相似。幾何相似是指模型與其原型形狀相同,但尺寸可以不同,而一切對應(yīng)的線性尺寸呈比例；運動相似是指對不同的流動現(xiàn)象,在流場中的所有對應(yīng)點處對應(yīng)的速度方向一致,且比值相等；動力相似是指作用在流體上相應(yīng)位置處各力組成的力多邊形幾何相似。在本研究中要保證實驗?zāi)Ｐ蛢?nèi)的介質(zhì)流動與熱風(fēng)箱原型中的氣體流動相似[10-11]。

符合物理現(xiàn)象相似的原則,應(yīng)做到以下幾方面：搭建實驗臺時,首先要保證氣體通道的幾何形狀相似,即實驗熱風(fēng)箱的尺寸與原型對應(yīng)呈比例；熱風(fēng)箱與原型的內(nèi)部流動為同類流動；進(jìn)入熱風(fēng)箱的氣體的物性參數(shù)相似；進(jìn)出界面的速度分布相似[12]。

關(guān)于特征數(shù)間的關(guān)系和現(xiàn)象的相似已進(jìn)行了大量的研究,結(jié)果表明，在黏性流體流動過程中,當(dāng)雷諾數(shù)Re大于某一值時,流動從層流狀態(tài)進(jìn)入湍流,這一數(shù)值稱為第一臨界雷諾值,Re小于第一臨界雷諾值的流動處于第一自模化區(qū),此時處在層流狀態(tài)。Re大于第一臨界雷諾值之后,隨著Re增大流體紊亂程度變化較大,當(dāng)其值增大到一定程度時這種影響便消失了,該Re稱第二臨界雷諾值。Re大于第二臨界雷諾值時流動處于第二自模化區(qū)[13]。

1　實驗?zāi)Ｐ?/h2>

在本實驗?zāi)Ｐ椭?研究對象為黏性不可壓縮定常流動熱空氣,熱風(fēng)箱系統(tǒng)示意圖見圖3。

圖3　熱風(fēng)箱系統(tǒng)模型示意圖

熱風(fēng)箱原型尺寸為：長度3720 mm、寬度400 mm、高度220 mm。

熱風(fēng)箱入口當(dāng)量直徑(de)：

熱空氣運動黏度為2.002×10-5m2/s,當(dāng)熱風(fēng)箱原型入口風(fēng)速為6 m/s時,雷諾數(shù)Re=93234.21,大于第二臨界雷諾值,流動處于第二自模化區(qū)。根據(jù)相似原理,熱風(fēng)箱模型按照原尺寸2∶5的比例制作,熱風(fēng)箱模型的外形尺寸分別為：長度1488 mm、寬度160 mm、進(jìn)口高度88 mm；熱風(fēng)箱入口當(dāng)量直徑(de)：

雷諾數(shù)Re表征的是流體的流動狀態(tài),Re=ρvl/η,其中v、ρ、η分別為流體的流速、密度與黏性系數(shù),l為熱風(fēng)箱長度。普朗特數(shù)代表了熱邊界層與流動邊界層的相對厚度,即流體中動量擴(kuò)散與熱量擴(kuò)散能力的對比。努塞爾數(shù)表征了流體對流換熱能力的大小。不同的換熱條件下,流體的換熱能力不相同。

表面?zhèn)鳠嵯禂?shù)計算公式為:

式中，A為紙張傳熱面積；A=1.488×0.16=0.2381(m2)；ta為紙張下底面不同位置測量溫度平均值，℃；tw為紙張溫度，℃；Q為熱空氣和紙張之間交換的熱量，J(包括水分汽化所吸收的熱量Q1,絕干紙張溫度升高所吸收熱量Q2,紙張中水分溫度升高所吸收熱量Q3,紙張對環(huán)境的散熱量Q4)。影響表面?zhèn)鳠嵯禂?shù)的主要因素見表1。

表1　影響表面?zhèn)鳠嵯禂?shù)的主要因素

實驗系統(tǒng)除熱風(fēng)箱模型主體外,還有動力系統(tǒng)、控制系統(tǒng)、加熱系統(tǒng)和測量系統(tǒng)。

動力系統(tǒng)：熱風(fēng)箱內(nèi)部空氣處于紊流狀態(tài),流動處在第二模化區(qū)。該實驗?zāi)Ｐ托枰L(fēng)量4500 m3/h,壓力4256 Pa、流量4792 m3/h的離心鼓風(fēng)機(jī),可完全滿足要求。

加熱系統(tǒng)：采用電加熱器把空氣加熱到一定的溫度,對紙幅進(jìn)行干燥。本實驗中采用的電加熱器的最大功率為33 kW,可將該流量的空氣加熱到150℃。

控制系統(tǒng)：采用XMT-3000系列智能專家PID工業(yè)控制/調(diào)節(jié)器控制溫度。調(diào)節(jié)器通過對電加熱器的通斷來控制溫度。當(dāng)熱電偶測量的空氣溫度高于設(shè)定值時,調(diào)節(jié)器可以把電加熱器暫時斷開,低于設(shè)定值時接通,通過控制使得空氣溫度恒定在一個設(shè)定值,其誤差小于1%。

測量系統(tǒng)：本實驗中主要測量的參數(shù)有紙張定量、熱風(fēng)箱和氣墊層靜壓、空氣流量、空氣和紙幅溫度。定量采用電子天平測量、風(fēng)速采用QFD-3熱球風(fēng)速儀測量,靜壓采用U型壓力計測量,溫度采用熱電阻測量。

2　量綱分析及擬合

在熱風(fēng)箱的結(jié)構(gòu)尺寸參數(shù)相同時,熱風(fēng)箱入口風(fēng)速越大,氣墊層產(chǎn)生的靜壓值就越大,所以紙幅被浮起高度越大。紙幅的浮起高度直接影響氣墊層的傳熱性能。在本研究中,都是在氣墊層靜壓值和紙張定量相等時的表面?zhèn)鳠嵯禂?shù),該表面?zhèn)鳠嵯禂?shù)值才有研究意義。在模擬計算時,先模擬出靜壓值與紙張定量相等時的氣墊層高度尺寸,然后在此高度條件下求解表面?zhèn)鳠嵯禂?shù)。

在表1中,含有4個基本量綱,質(zhì)量的量綱M、時間的量綱T、長度的量綱L和溫度的量綱θ。根據(jù)π定理得到式(1)和式(2)。

(1)

[MT-3θ-1]a[θ-1]b[LT-1]c[L2T-2θ-1]d[L2T-1]e[MLT-3θ-1]f[L]g[L2T-1]h=1

(2)

由量綱分析得：

a+f=0

(3)

-3a-c-2d-e-3f-h=0

(4)

c+2d+2e+f+g+2h=0

(5)

b-a-d-f=0

(6)

10個因素變量中有4個基本量綱,因此需要存在6個獨立準(zhǔn)則項。整理得：

a=-f

(7)

b=d

(8)

c=-2d-e-h

(9)

g=-f-e-h

(10)

式(1)化簡得：

(11)

(12)

整理式(2)得：

(13)

式中，Pr為普朗特數(shù)；Re為雷諾數(shù)；Nu為努塞爾數(shù)。

將式(13)兩邊取對數(shù)化為線性方程見式(14)。

表2　最小二乘法擬合相關(guān)數(shù)據(jù)表

表3　不同風(fēng)速對Re和Nu的影響

(14)

綜合上述各因素的影響,選取穩(wěn)定的實驗點通過最小二乘法求解,相關(guān)數(shù)據(jù)見表2和表3。

代入擬合得到：k=-0.60,a=0.62,b=0.47,c=-1.43,d=4.64,e=8.26,帶回整理得：

(15)

為驗證擬合方程(15)的準(zhǔn)確性,將數(shù)值模擬結(jié)果與擬合值進(jìn)行比較,不同工況下Nu數(shù)的模擬值和擬合值對比如圖4所示。從圖4可以看出,風(fēng)速8m/s、開孔率7.42%、高度比0.5時，Nu數(shù)模擬值69.05，擬合值75.63，此時相對誤差最大，為9.53%，均方根誤差為9.5%，Nu數(shù)值模擬結(jié)果與擬合值吻合較好,表明擬合方程完全可用于熱風(fēng)箱氣墊層表面?zhèn)鳠嵯禂?shù)的研究。

圖4　不同開孔率Nu數(shù)的模擬值和擬合值對比

3　結(jié)　論

(1)實驗表明紙張干燥中對流換熱能力的大小與普朗特數(shù)(Pr)和雷諾數(shù)(Re)有關(guān)。

(2)通過量綱分析擬合努塞爾數(shù)(Nu)數(shù)與Re、Pr函數(shù)關(guān)系式為：

(3)Nu數(shù)的模擬值與擬合值比較表明,風(fēng)速8m/s、開孔率7.42%、高度比0.5時，Nu數(shù)模擬值69.05，擬合值75.63，最大相對誤差為9.53%,均方根誤差為9.5%，Nu數(shù)值模擬結(jié)果與擬合值吻合較好，表明擬合方程完全可用于熱風(fēng)箱氣墊層表面?zhèn)鳠嵯禂?shù)的研究。

[1]HOUShun-li.SuperHeatedSteamImpingementDringTechnology[J].ChinaPulp&Paper, 2008, 27(3): 63.

侯順利. 過熱蒸汽沖擊干燥技術(shù)[J]. 中國造紙， 2008， 27(3)： 62.

[2]Lin Zhi-zuo. The paper drying process modeling research[D]. Guangzhou: South China University of Technology, 2012.

林治作. 紙頁干燥過程建模研究[D]. 廣州： 華南理工大學(xué), 2012.

[3]Zhen Rong-he. Research on flow in hot air box system air cushion pulp board drier[D]. Nanjing: Southeast University, 2008.

鄭榮和. 氣墊式漿板機(jī)熱風(fēng)箱系統(tǒng)氣體流動規(guī)律研究[D]. 南京： 東南大學(xué), 2008.

[4]Zeng Man-lian. Research on flow and heat transfer in hot air box system air cushion pulp board drier[D]. Nanjing: Southeast University, 2009.

曾滿連. 氣墊式漿板機(jī)熱風(fēng)箱系統(tǒng)氣體流動及傳熱特性研究[D]. 南京： 東南大學(xué), 2009.

[5]Zagar, Steve J. Computational fluid dynamics helps solve web stability problem[J]. Pulp & Paper Canada, 2006, 70(11): 110.

[6]Burch D M, Licata B A, Zarr R R. A Comparison of Two Test

Methods for Determining Transfer Function Coefficients for a Wall Using a Calibrated Hot Box[J]. Journal of Heat Transfer, 1990, 1(1): 112.

[7]Palfrey, Albert J. Thermal Performance of Low Admittance Building Sheathing[J]. Journal of Building Physics, 2006, 3(3): 129.

[8]Acharya S, Goldstein R J. Natural Convection in an Externally Heated Vertical or Inclined Square Box Containing Internal Energy Sources[J]. Journal of Heat Transfer, 2002, 53(53): 52.

[9]ZHANG Chang-fei, DING Ke-qiang, GE Shi-fu, et al. Research on Air Flow and Transfor Charactoristics of Hot Air Box in Air Cushion Drying System[J]. China Pulp & Paper, 2012, 31(9): 23.

張長飛, 丁克強(qiáng), 葛仕福, 等. 漿板機(jī)氣墊式熱風(fēng)箱氣體流動及傳熱特性模擬研究[J]. 中國造紙, 2012, 31(9): 23.

[10]Shrive A H, Mujumdar A S, Jujubes. Numerical study of heat transfer and fluid flow in multiple turbulent impinging jets[J]. Drying Technology, 2007, 15(5): 1343.

[11]Shi Yelling B Ray, A S Mujumdar. A parametric study of impingement heat transfer under a slot jet using CFD studies[C]//Proceedings of the 2nd Asian-Oceania Drying Conference, Batu Fearing i, Penang, M Malaysia, 2001.

[12]Devahastin S, Mujumdar A S. A Numerical Study of Flow and Mixing Characteristies of Laminar Confined Impinging Streams[J]. Chemical Engineering Javrnal, 2002, 85(85): 217.

[13]TANG Wei, ZHOU Yang, SU Ya-zhou. Optimization of Supply Air Temperature for the Paper Machine Dryer Section with Closed Steam Hood Based on Genetic Algorithm[J]. China Pulp & Paper, 2013, 32(11): 52.

(責(zé)任編輯：董鳳霞)

The Numerical Simulation Research on Heat Transfer in Hot Air Box System and Paper Sheet

ZHANG Chang-fei*DING Ke-qiangLI Qian-junZHANG Dong-pingWU Gong-De

(SchoolofEnvironmentalEngineering,NanjingInstituteofTechnology,Nanjing,JiangsuProvince, 211167)

The distribution of hot air speed, temperature and static pressure in the hot air box,and the influence of inlet velocity, height ratio, air temperature, percentage of opening and the style of jets on the uneven index, static pressure and surface heat transfer coefficient of the air cushion were discussed. With dimensional analysis, a functional relationship ofNuwas obtained. Comparing fitted values with numerical simulation values, the maximum relative error was 9.53%,root mean square error was 9.5%. The fitting formula could be used to study and calculate heat transfer coefficient in engineering cushion layer.

air cushion； hot air box； heat-transfer characteristic; numerical simulation

張長飛先生；講師；主要從事傳熱傳質(zhì)方面的研究。

2016- 02- 03(修改稿)

南京工程學(xué)院青年資金項目(QKJA2011007)。

TS733+.6

ADOI：10.11980/j.issn.0254- 508X.2016.07.011