竇秀敏,謝 東,*,查舒衛(wèi),王汝佳,王漢青
(1.南華大學(xué) 土木工程學(xué)院,湖南 衡陽 421001;2.南華大學(xué) 裝配式建筑節(jié)能技術(shù)湖南省重點(diǎn)實(shí)驗(yàn)室,湖南 衡陽421001;3.南華大學(xué) 建筑環(huán)境控制技術(shù)湖南省工程實(shí)驗(yàn)室,湖南 衡陽 421001;4.建筑環(huán)境氣載污染物治理與放射性防護(hù)國家地方聯(lián)合工程研究中心,湖南 衡陽 421001;5.中南林業(yè)科技大學(xué) 土木工程學(xué)院,湖南 長沙 410000)
工業(yè)廠房在生產(chǎn)過程中會(huì)產(chǎn)生大量的余熱、余濕,這嚴(yán)重影響了工人的健康。越來越多的研究表明長期處于高溫高濕條件下的工人會(huì)受到周圍熱輻射的影響,感到疲倦和不適,甚至?xí)a(chǎn)生中暑和暈倒等生理反應(yīng)[1]。然而以往的工業(yè)廠房設(shè)計(jì)主要考慮生產(chǎn)工藝的要求,而忽視了工人在廠房生產(chǎn)中的健康和舒適。因此,改善工業(yè)廠房的熱環(huán)境是必要且迫切的[2]。
通風(fēng)是改善室內(nèi)環(huán)境的重要手段,良好的通風(fēng)系統(tǒng)可以有效的改善室內(nèi)熱環(huán)境[3-4]。通風(fēng)方式可以分為全面通風(fēng)和局部通風(fēng)[5]。全面通風(fēng)是對(duì)整個(gè)房間進(jìn)行通風(fēng)換氣,其原理是,用清潔空氣稀釋室內(nèi)空間的有害氣體,同時(shí)不斷排出污染空氣。局部通風(fēng)是使用局部氣流,在局部區(qū)域形成良好的空氣環(huán)境。全面通風(fēng)包括自然通風(fēng)和機(jī)械通風(fēng),自然通風(fēng)在實(shí)際的工業(yè)建筑中可能由于進(jìn)風(fēng)口設(shè)計(jì)不合理或熱壓不足而不能達(dá)到較好的效果。機(jī)械通風(fēng)可以通過調(diào)節(jié)風(fēng)速彌補(bǔ)自然通風(fēng)的缺點(diǎn)。但是,較大的機(jī)械送風(fēng)速度可能會(huì)導(dǎo)致氣流短路,并降低通風(fēng)效率,增加能耗[6]。局部通風(fēng)包括局部送風(fēng)和局部排風(fēng),局部通風(fēng)所需通風(fēng)量小并且有針對(duì)性。然而局部通風(fēng)僅在較小范圍內(nèi)有效,并且對(duì)非封閉空間或傳統(tǒng)的銑削車間,局部排風(fēng)極其困難[7]。因此,通風(fēng)系統(tǒng)的合理設(shè)計(jì)對(duì)進(jìn)一步提高室內(nèi)空氣質(zhì)量并降低能源成本是至關(guān)重要的
辦公住宅類建筑的熱環(huán)境一直是人們關(guān)心的話題,因此針對(duì)辦公住宅建筑[8-11]通風(fēng)性能的研究已經(jīng)很多了。隨著人們對(duì)工人健康的關(guān)注,工業(yè)建筑通風(fēng)性能的研究也隨著開展。韓國浦項(xiàng)科技大學(xué)的J.H.Kang、S.J.Lee[12]提出在廠房進(jìn)風(fēng)口上部附加一定結(jié)構(gòu)和傾角的出風(fēng)口可以改善大型工業(yè)廠房內(nèi)部的空氣滯留問題。田力榮[13]對(duì)廠房內(nèi)排風(fēng)扇降低局部溫度的效果進(jìn)行研究發(fā)現(xiàn)由于結(jié)構(gòu)上的渦流作用,局部排風(fēng)扇反而使得局部的熱氣無法排除,增加了局部的平均溫度;較高的排風(fēng)速度,會(huì)使氣流短路加劇了局部溫升。王孟武[14]研究發(fā)現(xiàn)冗余的風(fēng)口設(shè)置會(huì)造成廠房內(nèi)溫度場(chǎng)均勻性變差,導(dǎo)致通風(fēng)效率及經(jīng)濟(jì)性的明顯下降。張大衛(wèi)[15]研究發(fā)現(xiàn)隨著進(jìn)排風(fēng)口中心高度差的增大,工作區(qū)溫度總體呈遞減趨勢(shì)。H.Q.Wang[16]在具有大功率熱源的煙廠中研究了置換通風(fēng)的去污效果,研究發(fā)現(xiàn)置換通風(fēng)形成向上的熱羽可以有效地去除顆粒物。S.J.Mei[17]分析了浮力驅(qū)動(dòng)的自然通風(fēng)和機(jī)械通風(fēng)條件下的熱環(huán)境改善效果,結(jié)果表明,在良好的氣流組織條件下,機(jī)械通風(fēng)可以有效的提高工作區(qū)域的換氣效率。從以上文獻(xiàn)綜述可知,送排風(fēng)參數(shù)(送風(fēng)和出風(fēng)口位置、送風(fēng)速度、溫度和角度)[18-21]、通風(fēng)方式[22-24]和室外參數(shù)[25]等是考慮較多的因素。針對(duì)舊廠房改造的研究相對(duì)較少,然而現(xiàn)存廠房中仍有大量廠房存在溫度高、溫度分布不均等情況,因此針對(duì)已有廠房的通風(fēng)改造研究也是十分急迫的。
本文以東莞某多熱源工業(yè)廠房為研究對(duì)象,提出增強(qiáng)全面通風(fēng)(全面通風(fēng)指對(duì)整個(gè)房間進(jìn)行通風(fēng)換氣)和局部送風(fēng)(局部送風(fēng)指向局部工作地點(diǎn)送風(fēng),使局部地帶造成良好的空氣環(huán)境)。兩種優(yōu)化方案,并運(yùn)用計(jì)算流體力學(xué)(computational fluid dynamics,CFD)方法對(duì)優(yōu)化前后廠房的熱環(huán)境進(jìn)行模擬。通過對(duì)比不同通風(fēng)方式下廠房內(nèi)熱環(huán)境的變化,以期探尋一種有效的改造方案,解決目前廠房存在的工作區(qū)溫度過高,熱環(huán)境惡劣的問題;并進(jìn)一步為工業(yè)廠房通風(fēng)系統(tǒng)的設(shè)計(jì)提供了依據(jù)。
本文以中國東莞市某紙漿模塑生產(chǎn)廠房為研究對(duì)象。東莞市位于南亞熱帶季風(fēng)氣候區(qū),一般最高溫度36 ℃,最低溫度1~3 ℃,平均室外空氣溫度為23.1 ℃,冬季盛行東風(fēng)和北風(fēng),夏季主要盛行南風(fēng)。廠房模型如圖1所示,廠房長78 m、寬30 m、高6 m。南、北側(cè)各有一扇門高2.2 m、寬1.4 m;西側(cè)墻體設(shè)有一扇門高3.2 m、長2.4 m,東側(cè)墻體設(shè)有一扇門高3.2 m、長3 m。東側(cè)墻體設(shè)有排風(fēng)扇長1.4 m、寬1.4 m,排風(fēng)扇距離地面3.6 m。在設(shè)計(jì)時(shí)考慮采光問題,在西側(cè)墻體設(shè)有高5 m的玻璃幕墻,距地面0.3 m。玻璃幕墻外為中庭,因此該側(cè)墻體受到的太陽輻射較少。該工廠包括9條生產(chǎn)線,每條生產(chǎn)線的熱源包括10臺(tái)設(shè)備和一根蒸汽管道,單個(gè)設(shè)備長1 m、寬0.8 m、高2 m,設(shè)備間距0.5 m,設(shè)備局部最高溫度可達(dá)200 ℃。設(shè)備朝西1 m范圍內(nèi)為工作區(qū)域,其他區(qū)域?yàn)檫^道。為提高網(wǎng)格離散效率及數(shù)值計(jì)算精確度,對(duì)廠房模型進(jìn)行簡化,人員熱量和蒸汽管道的熱量附加到設(shè)備上,照明燈具進(jìn)行忽略,房間內(nèi)其他面進(jìn)行封閉處理。目前該廠房的通風(fēng)系統(tǒng)由4個(gè)門和11個(gè)排風(fēng)扇組成,由于通風(fēng)口數(shù)量較少,氣流組織效果不佳,室內(nèi)熱環(huán)境較差。
圖1 廠房模型圖(Case 1)
本文基于該廠房多熱源的特點(diǎn),同時(shí)考慮可行性和經(jīng)濟(jì)性,提出了兩種廠房通風(fēng)系統(tǒng)的優(yōu)化方案:
1)增強(qiáng)全面通風(fēng);開啟廠房西面的10個(gè)窗戶,窗戶尺寸為3 000 mm×2 000 mm,增加室內(nèi)的自然通風(fēng),改造后模型如圖2所示。
圖2 優(yōu)化方案一廠房模型圖(Case 2)
2)增加局部送風(fēng);在人員活動(dòng)區(qū)增加送風(fēng),送風(fēng)口尺寸為350 mm×350 mm,間隔2 400 mm,優(yōu)化后模型如圖3所示。
圖3 優(yōu)化方案二廠房模型圖(Case 3)
工業(yè)廠房中的氣流流動(dòng)和溫度分布基于質(zhì)量、動(dòng)量和能量守恒定律。使用的數(shù)學(xué)模型包括連續(xù)性方程、動(dòng)量方程和能量方程一般形式表示如下:
(1)
其中,?為通用變量,可以表示諸如u,v,w,T等變量;Γ?為對(duì)應(yīng)?的擴(kuò)散系數(shù);S?為源項(xiàng)。
(2)
其中,a為吸收系數(shù),n為折射率,I為輻射強(qiáng)度,φ為相位函數(shù),Ω′為立體角,δs為散射系數(shù)。
本文研究的流場(chǎng)屬于湍流流場(chǎng),采用Realizable k-e湍流模型。壁面附近采用標(biāo)準(zhǔn)壁面函數(shù)法進(jìn)行處理;考慮熱源和壁面之間存在輻射換熱,采用DO輻射模型。并做了如下假設(shè):1)流體性質(zhì)滿足Boussinesq假設(shè),密度隨溫度線性變化,其他性質(zhì)恒定;2)流動(dòng)為單相穩(wěn)態(tài)流;3)建筑不漏風(fēng)。
根據(jù)機(jī)器的設(shè)計(jì)功率、蒸汽管和人員的散熱量計(jì)算得出熱源為850 W/m2的恒定熱流。由于門窗是工業(yè)建筑的自然通風(fēng)通風(fēng)口,風(fēng)速不確定,所以門窗是壓力入口的邊界。由于出風(fēng)口為排風(fēng)機(jī),出風(fēng)口為速度入口邊界,速度入口邊界為-2 m/s,為現(xiàn)場(chǎng)實(shí)測(cè)值,壁面均為無滑移壁面,邊界條件具體參數(shù)見表1。
表1 邊界條件
非結(jié)構(gòu)化網(wǎng)格是使用ANSYS FluentMeshing創(chuàng)建的。所使用的網(wǎng)格是不均勻的,對(duì)于廠房送、排風(fēng)口處速度變化梯度較大和且設(shè)備熱源局部溫度梯度變化大的區(qū)域進(jìn)行網(wǎng)格加密。并在熱源設(shè)備周圍增加了邊界層。圖4給出了局部的網(wǎng)格視圖。
圖4 局部網(wǎng)格視圖
為確保計(jì)算結(jié)果不受網(wǎng)格精度的影響,分別采用80萬、168萬、222萬共3種精度的網(wǎng)格對(duì)原物理模型進(jìn)行數(shù)值模擬。計(jì)算結(jié)果顯示,3種網(wǎng)格精度下的通風(fēng)量分別為51.691 kg/s,52.597 kg/s,52.604 kg/s,其中168萬與222萬網(wǎng)格數(shù)量的計(jì)算結(jié)果差別最小,僅為0.13%。因此權(quán)衡網(wǎng)格精度和計(jì)算效率之后,本文采用168萬的網(wǎng)格劃分方式。
本文采用SIMPLE算法完成壓力-速度耦合求解,壓力項(xiàng)取body-force-weighted格式,動(dòng)量、溫度等項(xiàng)取二階格式,湍流動(dòng)能及耗散率取一階格式。在具有兩個(gè)Intel Xeon-Sliver-4116處理器和192 GB RAM的工作站上使用并行處理來執(zhí)行計(jì)算。每次仿真大約需要30 h能完成。
在進(jìn)行數(shù)值模擬之前,有必要對(duì)數(shù)值模型進(jìn)行驗(yàn)證。為確保本文數(shù)值模擬結(jié)果的合理性與可靠性,本文將數(shù)值模擬結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)的數(shù)據(jù)進(jìn)行對(duì)比驗(yàn)證。
實(shí)驗(yàn)儀器和測(cè)量變量如表2所示。為減少太陽輻射的影響,將WatchDog溫度記錄儀放置在距離地面0.8 m的陰涼處。為監(jiān)測(cè)工作人員的工作環(huán)境,將Testo625溫度測(cè)量儀放置在距地面1.5 m的工作區(qū)域。實(shí)驗(yàn)測(cè)點(diǎn)分布如圖5所示。
圖5 現(xiàn)場(chǎng)測(cè)點(diǎn)布置圖
表2 實(shí)驗(yàn)儀器和測(cè)量變量
將實(shí)驗(yàn)測(cè)量的溫度值與數(shù)值模擬得到的溫度值進(jìn)行對(duì)比如圖6所示。由圖6可知,各測(cè)點(diǎn)模擬值與實(shí)測(cè)值有略微差異,模擬值與實(shí)測(cè)值的絕對(duì)誤差在0.1~1.6 ℃,最大相對(duì)誤差僅為3.9%。誤差原因在于本文對(duì)廠房模型進(jìn)行了簡化,同時(shí)實(shí)驗(yàn)儀器也存在一定測(cè)量誤差。兩者差異較小,說明所建CFD模型可用于進(jìn)一步的通風(fēng)研究。
圖6 實(shí)測(cè)溫度與模擬溫度對(duì)比圖
2.2.1 室內(nèi)溫度
圖7是Case1、Case2和Case3在水平方向(X-Y)Z=1.6 m截面的溫度云圖,平均溫度分別為38.9 ℃、36.3 ℃、37.1 ℃。由圖7可以看出,改造前廠房靠近北側(cè)的工作區(qū)溫度相對(duì)較高,最高溫度可達(dá)到43 ℃,主要原因是靠近北側(cè)的工作區(qū)僅有一個(gè)側(cè)門作為通風(fēng)口,室外冷空氣無法流過部分工作區(qū),熱量不能順利帶走;另外熱源相當(dāng)于一個(gè)阻礙物,對(duì)進(jìn)風(fēng)氣流產(chǎn)生了阻礙作用。增強(qiáng)全面通風(fēng)之后,廠房內(nèi)整體溫度均有所下降,熱環(huán)境得到改善,工作區(qū)局部高溫的現(xiàn)象有所緩解。靠近窗戶一側(cè)的工作區(qū)溫度基本保持在33~34 ℃。增加局部送風(fēng)之后,加快了工作區(qū)的空氣流動(dòng),帶走更多的熱量,廠房工作區(qū)的溫度有所下降,但局部區(qū)域還是存在高溫現(xiàn)象,最高溫度可達(dá)40 ℃。增強(qiáng)全面通風(fēng)之后不僅降低了工作區(qū)的溫度,整體廠房的溫度也有所降低,對(duì)改善工作區(qū)熱環(huán)境的作用更加明顯。
圖7 三種通風(fēng)系統(tǒng)水平(X-Y)面(Z=1.6 m)溫度的比較
圖8為Y=14 m、X=50 m處溫度隨高度變化圖,圖9為工作區(qū)測(cè)點(diǎn)對(duì)應(yīng)的溫度變化圖。由圖8可以看出,由于熱浮力的作用,溫度隨高度的升高而上升,出現(xiàn)分層現(xiàn)象。改造后廠房Z≤2 m區(qū)域的溫度明顯低于改造前廠房的溫度。其中Case2溫度在高度方向的變化最小,上下溫差將近3 ℃;Case3溫度在高度方向的變化最大,上下溫差將近4 ℃??梢钥闯鲈黾泳植克惋L(fēng)有效改善工作區(qū)的熱環(huán)境,而加強(qiáng)全面通風(fēng),使得整體廠房的溫度均有所降低,整個(gè)廠房熱環(huán)境得到改善。由圖9可以看出,改造前靠北側(cè)的測(cè)點(diǎn)(測(cè)點(diǎn)1、2、3、13、14)溫度相對(duì)較高,主要是由于該側(cè)僅有一個(gè)側(cè)門作為通風(fēng)口,氣流組織不均勻。增加局部送風(fēng)后,溫度有所降低。增強(qiáng)全面通風(fēng)后,通風(fēng)口數(shù)量增加,使得氣流組織的均勻性提高,整體廠房溫度降低。另外,增加窗戶后,靠近窗戶一側(cè)的測(cè)點(diǎn)(測(cè)點(diǎn)1-12)溫度相對(duì)較低,均在33 ℃左右。
圖8 高度方向的溫度分布圖
圖9 測(cè)點(diǎn)處溫度分布圖
2.2.2 室內(nèi)空氣流速
圖10為Case1、Case2和Case3在z=1.6 m截面的速度云圖,該截面處平均速度分別為0.45 m/s、0.3 m/s、0.33 m/s。由圖10可以看出,靠近通風(fēng)口的區(qū)域氣流速度相對(duì)較高,最高可達(dá)1.8 m/s,由于障礙物的阻礙作用,部分區(qū)域氣流流速較低。增強(qiáng)全面通風(fēng)后,廠房內(nèi)空氣流速均有降低,最高只達(dá)到0.7 m/s;主要是因?yàn)?,總進(jìn)風(fēng)量相同的情況下,進(jìn)風(fēng)口增多,每個(gè)進(jìn)風(fēng)口的空氣流速降低。增加局部送風(fēng)之后,工作區(qū)局部風(fēng)速增大,高達(dá)1.8 m/s,自然進(jìn)風(fēng)口的風(fēng)速有所降低。增加局部通風(fēng)和增強(qiáng)全面通風(fēng)后,速度分布更為均勻,這說明改造后,廠房氣流組織得到優(yōu)化,增加了空氣的流通,使得廠房熱環(huán)境得到改善。
圖10 三種通風(fēng)系統(tǒng)水平(X-Y)面(Z=1.6 m)速度的比較
2.2.3 工作區(qū)熱環(huán)境
長期處在高熱濕環(huán)境的工人對(duì)環(huán)境具有一定的熱適應(yīng)性[26-27],在具有失調(diào)危險(xiǎn)的熱環(huán)境中,以人體熱舒適為主的熱環(huán)境指標(biāo)不再適用[28]。許多學(xué)者建立熱應(yīng)力指標(biāo)與人體健康的關(guān)系,提出了多種不同的評(píng)價(jià)方法及指標(biāo)。本文選用1955年由Belding[29]提出熱應(yīng)力指數(shù)(heat stress index, HIS),其對(duì)環(huán)境的熱應(yīng)力評(píng)價(jià)反映了空氣流速、室內(nèi)溫度、濕球溫度和輻射溫度,還有新陳代謝產(chǎn)熱量的綜合作用。它是指為保持身體熱平衡所需的蒸發(fā)散熱與可以透過服裝散發(fā)到環(huán)境中去的最大蒸發(fā)散熱量的比值。熱應(yīng)力指數(shù)計(jì)算方法如下所示:
熱應(yīng)力指數(shù)Hsi:
(3)
允許暴露時(shí)間Tae:
(4)
最大蒸發(fā)熱損失Emax:
Emax=11.7×u0.6×(56-Pa)
(5)
所需排汗量Ereq:
Ereq=M-R-C
(6)
輻射熱損失R:
R=7.3×(35-TMRT)
(7)
對(duì)流熱損失C:
C=7.6×u0.6×(35-Tn)
(8)
其中,TMRT為室內(nèi)平均輻射溫度,℃;Tn為工作區(qū)平均溫度,℃;u為工作區(qū)平均氣流速度,m/s;Pa為空氣中的水蒸氣分壓力,Pa;M為新陳代謝產(chǎn)熱量,W/m2。
假定人體皮膚的溫度恒定在35 ℃,不計(jì)呼吸散熱。根據(jù)實(shí)測(cè)[28],水蒸氣分壓力取32 Pa。假定皮膚溫度恒定為35 ℃,并規(guī)定最大蒸發(fā)熱損失Emax的上限值為390 W/m2,相當(dāng)于一個(gè)典型男子的排汗量為1.0 L/h。用熱力指數(shù)100%作為可調(diào)節(jié)區(qū)域的上限。如果Ereq>Emax,人體就不能再保持熱平衡,體溫開始上升。在這樣的一種環(huán)境中長時(shí)間工作是不可能的,而只能在有限的時(shí)間內(nèi)進(jìn)行工作,必須在體內(nèi)蓄熱引起體內(nèi)溫升達(dá)到危險(xiǎn)程度之前即停止暴露。
表3為不同通風(fēng)方式下人員活動(dòng)區(qū)域(Z≤2 m)的熱應(yīng)力指數(shù)。由表3可知增強(qiáng)全面通風(fēng)對(duì)熱環(huán)境改善效果最好。人體的對(duì)流散熱量取決于工作區(qū)氣流溫度和流速,溫度為主要影響因素。改造后廠房工作區(qū)的平均風(fēng)速略有降低,但工作區(qū)平均溫度有所下降,導(dǎo)致人體對(duì)流換熱量有所降低。人體輻射得熱完全取決于工作區(qū)輻射溫度。改造后廠房工作區(qū)通風(fēng)效果增強(qiáng),帶走了更多的余熱量,使得維護(hù)結(jié)構(gòu)內(nèi)表面溫度降低,因此人體輻射得熱量下降。本文假定水蒸汽分壓力為定值,人體最大蒸發(fā)量Emax完全取決于空氣流速,并且與風(fēng)速成正比。改造前后,廠房工作區(qū)風(fēng)速變化較小,使得人體最大蒸發(fā)熱損失變化較小。由表3可以看出,人體輻射得熱量和所需排汗量的減小對(duì)于工作區(qū)熱應(yīng)力指數(shù)Hsi的降低起到了決定性作用。改造后熱應(yīng)力指數(shù)Hsi減小,人員的可允許工作時(shí)間增大。
表3 不同通風(fēng)情況下熱應(yīng)力指數(shù)(Hsi)
本文以東莞市某紙漿模塑生產(chǎn)廠房為研究對(duì)象,結(jié)合廠房實(shí)際情況提出增強(qiáng)全面通風(fēng)和增加局部送風(fēng)兩種通風(fēng)改造方式。并利用CFD方法對(duì)不同通風(fēng)方式下的廠房內(nèi)熱環(huán)境進(jìn)行模擬研究,確定最佳通風(fēng)方式。利用實(shí)測(cè)數(shù)據(jù)對(duì)數(shù)值模擬結(jié)果進(jìn)行了驗(yàn)證,模擬值與實(shí)測(cè)值的絕對(duì)誤差在0.1~1.6 ℃,最大相對(duì)誤差僅為3.9%,模擬值與實(shí)測(cè)值吻合良好。
分析模擬結(jié)果可知,原廠房工作區(qū)平均溫度為38.9 ℃,平均風(fēng)速為0.43 m/s,熱應(yīng)力指數(shù)為127.12%,允許暴露時(shí)間為53.2 min。增強(qiáng)全面通風(fēng)后工作區(qū)平均溫度為35.1 ℃,平均風(fēng)速為0.41 m/s,熱應(yīng)力指數(shù)為117.32%,允許暴露時(shí)間為85.7 min。增加局部送風(fēng)后工作區(qū)平均溫度為36.3 ℃,平均風(fēng)速為0.42 m/s,熱應(yīng)力指數(shù)為119.0%,允許暴露時(shí)間為77.1 min。兩種通風(fēng)改造方式均能有效的提高室內(nèi)熱環(huán)境。增加局部送風(fēng),主要是改善人員活動(dòng)區(qū)域的熱環(huán)境,提高工作人員的熱舒適。全面通風(fēng)則使室內(nèi)的整體熱環(huán)境得到改善,使室內(nèi)氣流和溫度分布更加均勻。從模擬結(jié)果可以看出,增強(qiáng)全面通風(fēng),對(duì)室內(nèi)熱環(huán)境的改善效果更為明顯,并且可以有效的節(jié)約能源。