陳匯龍，謝曉鳳，廖文言，桂鎧，謝昀彤，張杰，趙斌娟*

(1. 江蘇大學(xué)能源與動(dòng)力工程學(xué)院，江蘇 鎮(zhèn)江 212013； 2. 應(yīng)急管理部上海消防研究所，上海 200032)

水力消排煙機(jī)是一種新型的防爆型消防裝備，具有體積小、重量輕、移動(dòng)方便、消煙效果好等優(yōu)勢(shì)，適用于各類建筑特別是地下建筑、無(wú)窗建筑等區(qū)域的消防救援，以及化學(xué)品生產(chǎn)、加工、存放等區(qū)域的安全防爆[1-2]，也是備受關(guān)注的火場(chǎng)偵察和救災(zāi)的理想輔助器材.然而，現(xiàn)有水力消排煙機(jī)[3]效率只有40%左右，能耗高、排煙量小，不能滿足復(fù)雜或者大型火災(zāi)現(xiàn)場(chǎng)的排煙需求，已成為水力消排煙機(jī)推廣應(yīng)用的瓶頸，亟待深入研究.

目前，對(duì)水力消排煙機(jī)所做改進(jìn)和優(yōu)化研究，大多通過(guò)試驗(yàn)進(jìn)行.席慧敏[4]通過(guò)調(diào)研移動(dòng)排煙設(shè)備的配備和使用狀況，分析了移動(dòng)排煙設(shè)備在使用過(guò)程中存在的問(wèn)題以及遇到的主要困難，并提出了解決的具體方法和對(duì)策，為配備和使用移動(dòng)排煙設(shè)備提供了參考；李思成等[5]對(duì)移動(dòng)式火場(chǎng)排煙存在的問(wèn)題進(jìn)行了分析，主要包括排煙裝備的配備情況、火場(chǎng)排煙訓(xùn)練的開(kāi)展情況、影響指揮員排煙意識(shí)的因素、實(shí)戰(zhàn)中火場(chǎng)排煙的應(yīng)用情況等，并提出了提高移動(dòng)式火場(chǎng)排煙效果的對(duì)策；周娟利等[6]以水力排煙水輪葉片芯金屬模具逆向反求為例，設(shè)計(jì)了水輪葉片反求蠟形，介紹了曲面重構(gòu)的實(shí)施流程.然而，試驗(yàn)研究的工作量和經(jīng)費(fèi)投入大，過(guò)程復(fù)雜，周期長(zhǎng)且難以取得成效.

結(jié)合試驗(yàn)深入開(kāi)展理論研究是提升水力消排煙機(jī)性能的有效途徑，具有重要的學(xué)術(shù)價(jià)值和工程實(shí)際意義.水力消排煙機(jī)以消防水泵等的壓力水作為動(dòng)力，使水輪高速旋轉(zhuǎn)，從而驅(qū)動(dòng)風(fēng)機(jī)產(chǎn)生風(fēng)量，實(shí)現(xiàn)消排煙的功能.其中，水輪是水力消排煙機(jī)中的關(guān)鍵水力部件，其性能的優(yōu)劣直接決定了消排煙機(jī)的整機(jī)性能.因此，文中將單獨(dú)對(duì)移動(dòng)式水力消排煙機(jī)中的水輪進(jìn)行模擬計(jì)算，以揭示水輪內(nèi)部水力損失機(jī)理，并系統(tǒng)研究水輪幾何參數(shù)的變化對(duì)水輪性能的影響.

1 水輪主要結(jié)構(gòu)參數(shù)

水力消排煙機(jī)具有消煙和排煙功能，主要部件包括水輪組件(包含水輪固定管、水輪噴管和連接套筒)、風(fēng)機(jī)葉片、風(fēng)機(jī)網(wǎng)罩(包含前防護(hù)罩、風(fēng)機(jī)罩殼、后防護(hù)罩)和推車裝置4部分，如圖1所示.

圖1 水力消排煙機(jī)結(jié)構(gòu)示意圖Fig.1 Structural schematic of hydraulic fume exhauster

如圖1所示水力消排煙機(jī)以消防水泵等的壓力水作為動(dòng)力，通過(guò)消防水帶連接水輪固定管，壓力水流經(jīng)水輪噴管，從噴管上的噴口噴出，使水輪高速旋轉(zhuǎn)，并驅(qū)動(dòng)風(fēng)機(jī)葉片轉(zhuǎn)動(dòng)產(chǎn)生風(fēng)量，達(dá)到火場(chǎng)排煙效果；同時(shí)，壓力水經(jīng)水輪霧化噴頭噴出產(chǎn)生大量的細(xì)水霧，達(dá)到火場(chǎng)消煙效果.水輪噴管與水輪固定管通過(guò)連接套筒進(jìn)行連接.前防護(hù)罩、后防護(hù)罩和風(fēng)機(jī)罩殼起到保護(hù)水輪及風(fēng)機(jī)葉片的作用，且具有一定的導(dǎo)流作用，推車裝置便于消防員移動(dòng)水力消排煙機(jī).其中由水輪噴管、水輪固定管以及連接套筒組成的水輪組件是水力消排煙機(jī)整體結(jié)構(gòu)中關(guān)鍵的能量轉(zhuǎn)換部件，將水能轉(zhuǎn)換為機(jī)械能，驅(qū)動(dòng)同軸風(fēng)機(jī)工作，進(jìn)而起到消排煙的功能.

水力消排煙機(jī)的水輪結(jié)構(gòu)如圖2所示(圖2a中的箭頭為高壓水流動(dòng)方向)，主要過(guò)流部件包括固定管、噴管、連接套筒，噴管與固定管之間采用螺紋連接.水輪是對(duì)稱結(jié)構(gòu)，其中固定管內(nèi)徑d1為45 mm、長(zhǎng)L1為250 mm；噴管內(nèi)徑d2為30 mm、長(zhǎng)L2為433 mm，根數(shù)N為2根，對(duì)稱均布；噴口是兩端直徑d3為4 mm、邊長(zhǎng)L3為20 mm的腰形孔口，噴口中心到噴管邊緣的距離L4為38 mm.

圖2 水輪結(jié)構(gòu)示意圖Fig.2 Structural schematic of hydro-turbine

2 數(shù)值模擬及結(jié)果驗(yàn)證

為了研究水輪的內(nèi)流場(chǎng)特性，獲得影響水輪效率的關(guān)鍵幾何參數(shù)，以試驗(yàn)所用水輪為研究對(duì)象進(jìn)行流道建模和內(nèi)流場(chǎng)數(shù)值模擬分析[7-8].

2.1 計(jì)算區(qū)域及網(wǎng)格劃分

結(jié)合圖1和2中的水輪結(jié)構(gòu)，建立如圖3所示的水輪流道模型.由圖可見(jiàn)按固體壁面邊界條件，整個(gè)計(jì)算區(qū)域分為2部分：一部分為固定管內(nèi)流道的固定域，另一部分為噴管內(nèi)流道的旋轉(zhuǎn)域，動(dòng)靜交界面為直徑45 mm的圓.為了使出口流動(dòng)充分發(fā)展，噴孔沿出流方向延長(zhǎng)10 mm.

圖3 水輪的計(jì)算區(qū)域Fig.3 Computational domain of hydro-turbine

水輪全流域采用四面體非結(jié)構(gòu)網(wǎng)格，以水輪效率η為檢驗(yàn)值進(jìn)行網(wǎng)格無(wú)關(guān)性檢驗(yàn)，結(jié)果如圖4所示.由圖可見(jiàn)，網(wǎng)格數(shù)量達(dá)到3.00×106后水輪效率基本不變，滿足網(wǎng)格無(wú)關(guān)性要求.此時(shí)，噴管水體網(wǎng)格數(shù)約為2.57×106，固定管水體網(wǎng)格數(shù)約為8.00×105.

圖4 網(wǎng)格無(wú)關(guān)性檢驗(yàn)Fig.4 Grid size independence verification

2.2 計(jì)算方法和邊界條件

水輪的工作介質(zhì)為水，雷諾數(shù)Re的數(shù)量級(jí)為104～105，說(shuō)明其內(nèi)流道流體流動(dòng)屬于湍流.針對(duì)水輪內(nèi)部流體流動(dòng)的特點(diǎn)假定：① 流動(dòng)過(guò)程中的流體不可壓縮；② 流動(dòng)參數(shù)的時(shí)均值不隨時(shí)間變化，采用RANS方法進(jìn)行水輪內(nèi)流場(chǎng)數(shù)值模擬.其中，控制方程組為定常、不可壓三維黏性流動(dòng)的連續(xù)性方程和動(dòng)量守恒方程(Navier-Stoke方程)；湍流模型為SSTk-ω湍流模型，該模型在標(biāo)準(zhǔn)k-ω模型的基礎(chǔ)上進(jìn)行修正得到，MEHRDAD等[9]對(duì)SST模型進(jìn)行分析，發(fā)現(xiàn)k-ω模型在近壁區(qū)計(jì)算中能較好地適用于各種壓力梯度下的邊界層問(wèn)題.

基于ANSYS CFX計(jì)算軟件，選取隱式分離求解器，旋轉(zhuǎn)域給定轉(zhuǎn)速n，壁面設(shè)為旋轉(zhuǎn)壁面且相對(duì)轉(zhuǎn)速為0；進(jìn)口采用靜壓入口方式，給定工作壓力；出口與大氣直接相通，故設(shè)為靜壓出口，相對(duì)壓力為0；動(dòng)靜交界面類型采用Genneral Connection形式；時(shí)間步長(zhǎng)為60/(2πn)，殘差取1.0×10-5.

2.3 水輪性能測(cè)試

為了給水輪內(nèi)流數(shù)值模擬提供必要的邊界條件并驗(yàn)證計(jì)算結(jié)果的準(zhǔn)確性，還需對(duì)水輪進(jìn)行性能測(cè)試.測(cè)試采用如圖5所示的試驗(yàn)裝置.試驗(yàn)時(shí)，通過(guò)調(diào)節(jié)消防水泵控制進(jìn)口壓力，高壓水通過(guò)消防水帶送入水輪固定管進(jìn)口處，先后流入水輪固定管、噴管，最后從噴管的噴口處高速噴出.此時(shí)，在壓力水的反作用下，噴管高速旋轉(zhuǎn)，同時(shí)帶動(dòng)同軸風(fēng)葉轉(zhuǎn)動(dòng)產(chǎn)生風(fēng)量，將噴管中噴出的水吹向前方.

圖5 試驗(yàn)裝置示意圖Fig.5 Schematic of experimental set-up

試驗(yàn)時(shí)，調(diào)節(jié)進(jìn)口靜壓分別為0.5，0.8和1.0 MPa，當(dāng)水輪運(yùn)行穩(wěn)定時(shí)，測(cè)得流量Q、扭矩T及轉(zhuǎn)速n，并根據(jù)式(1)計(jì)算水輪的效率η.

η=(Tn/9 550)/[Q(p1-p2)]，

(1)

式中：T為扭矩，N·m；Q為流量，kg/s；p1為進(jìn)口壓力，MPa；p2為出口壓力，MPa.

2.4 計(jì)算結(jié)果驗(yàn)證及討論

采用和性能試驗(yàn)一致的邊界條件，進(jìn)行水輪內(nèi)流場(chǎng)數(shù)值模擬，并對(duì)比數(shù)值模擬計(jì)算值與性能試驗(yàn)實(shí)測(cè)數(shù)據(jù)，結(jié)果見(jiàn)表1.可見(jiàn)流量Q及扭矩T的相對(duì)誤差σ均低于0.2%.因此所建立的計(jì)算模型是可靠的，計(jì)算結(jié)果能夠反映實(shí)際流動(dòng)情況.

表1 水輪計(jì)算及實(shí)測(cè)結(jié)果對(duì)比Tab.1 Comparison of simulated and tested hydro-turbine performance data

圖6為噴管速度云圖.由圖6a可見(jiàn)，水輪在相對(duì)較大的壓力及轉(zhuǎn)速條件下工作時(shí)，流道內(nèi)流速的均勻性更好，出現(xiàn)擾流的區(qū)域更少.

圖6 噴管速度云圖Fig.6 Velocity contour in nozzle

由圖6b可見(jiàn)，3個(gè)工況下的出口流速分布特征類似，均為近軸側(cè)速度梯度較大，高速區(qū)偏向出口斷面的近軸側(cè)區(qū)域，并且工作壓力和轉(zhuǎn)速提高，則出口斷面總體流速增大.由此可知，水輪內(nèi)流體擾動(dòng)集中在噴管位置，因而推測(cè)噴管幾何參數(shù)對(duì)水輪性能影響更大，而噴管主要幾何參數(shù)為噴管長(zhǎng)度、噴管直徑和噴口尺寸.

3 正交試驗(yàn)方案

實(shí)際工作中，通常要考慮3個(gè)或者3個(gè)以上試驗(yàn)因素的共同作用，正交試驗(yàn)設(shè)計(jì)的思路就是能夠安排多因素試驗(yàn)，分析各種因素的顯著性及敏感性程度[10-11].故以正交試驗(yàn)研究水輪幾何參數(shù)對(duì)水輪性能的影響.

水輪正交水力優(yōu)化設(shè)計(jì)目的是探求各結(jié)構(gòu)參數(shù)對(duì)水輪效率的影響程度，即目標(biāo)參數(shù)為水輪的效率.試驗(yàn)主要研究水輪噴口面積S、噴管直徑d2、噴管長(zhǎng)度L2等3個(gè)因素對(duì)水輪效率的影響，其中水輪噴管尺寸示意圖如圖2所示.

根據(jù)文獻(xiàn)[12-13]和經(jīng)驗(yàn)，每個(gè)因素分別取3個(gè)水平，見(jiàn)表2；選用 L9(33)正交表確定的正交設(shè)計(jì)方案見(jiàn)表3.原型水輪的因素組合即為A2B2C2.

表2 正交因素及水平表Tab.2 Factors and their levels in orthogonal experiments

表3 正交試驗(yàn)方案表Tab.3 Orthogonal experiment design cases

4 正交試驗(yàn)結(jié)果及分析

4.1 計(jì)算結(jié)果驗(yàn)證及討論

通過(guò)數(shù)值模擬，得到了9組正交試驗(yàn)方案的效率數(shù)值模擬計(jì)算值，見(jiàn)表4.

表4 試驗(yàn)方案的效率計(jì)算值Tab.4 Predicted efficiency of experiment design cases

4.2 極差分析

為了分析幾何參數(shù)對(duì)水輪的性能影響權(quán)重，并找到影響效率的主要因素和優(yōu)化方案，需要對(duì)模擬結(jié)果進(jìn)行極差分析，分析結(jié)果見(jiàn)表5，表中Ki為相應(yīng)水平i的幾種模擬情況指標(biāo)之和，而Mi為相應(yīng)水平i的指標(biāo)平均值，R為各平均值的極差.極差主要反映了每個(gè)因素下所選取的水平對(duì)試驗(yàn)指標(biāo)影響權(quán)重的大小.極差越大，說(shuō)明該因素的水平變化對(duì)試驗(yàn)指標(biāo)的影響權(quán)重越大.

表5 數(shù)值模擬結(jié)果極差分析Tab.5 Range analysis of numerical results

為了更直觀地反映試驗(yàn)因素水平改變對(duì)試驗(yàn)結(jié)果的影響，以因素的變化水平為橫坐標(biāo)、指標(biāo)平均值η*為縱坐標(biāo)所做的折線圖如圖7所示.

圖7 指標(biāo)與各因素關(guān)系Fig.7 Relationships between turbine efficiency and factors

由圖7可以看出，水輪3個(gè)變量對(duì)其水力性能影響的程度不同.因此，通過(guò)比較極差大小，可以獲得噴口面積、噴管直徑、噴管長(zhǎng)度對(duì)水輪水力效率的影響程度.所選幾何參數(shù)對(duì)水輪效率影響的主次順序?yàn)镃，B，A.

對(duì)于單個(gè)因素，通過(guò)比較該因素下各個(gè)水平試驗(yàn)指標(biāo)數(shù)值的大小，可以得出該因素各個(gè)水平對(duì)該試驗(yàn)指標(biāo)影響的主次順序.以因素A各水平對(duì)效率試驗(yàn)指標(biāo)的影響為例，K3>K2>K1，因此因素A各個(gè)水平對(duì)效率影響的主次順序是A2，A1，A3.同理，其他2個(gè)因素的3個(gè)水平對(duì)效率影響的主次順序分別為B2，B3，B1以及C1，C2，C3.綜上分析，可得所選參數(shù)中最優(yōu)的因素組合為A2B2C1.

最佳方案的噴口面積S為92.57 mm2、噴管直徑d2為30 mm、噴管長(zhǎng)度L2為438 mm，與原型相比噴管長(zhǎng)度減少了10%.利用ANSYS CFX軟件對(duì)最佳方案進(jìn)行全流道的數(shù)值模擬計(jì)算，將優(yōu)化后的水輪效率與優(yōu)化前對(duì)比，結(jié)果見(jiàn)表6.由表6可見(jiàn)，優(yōu)化后水輪的流量減小0.76 kg/s、扭矩增大2.79 N·m、效率提高7.67%.

表6 優(yōu)化前后模擬結(jié)果Tab.6 Simulated results before and after optimization

圖8為優(yōu)化后噴管的中心截面速度云圖.由圖可見(jiàn)，優(yōu)化后水輪噴管高速區(qū)域增大、低速區(qū)域減少，速度分布更加均勻，流體流動(dòng)狀態(tài)更優(yōu)，因此水輪效率得到明顯提升.

圖8 優(yōu)化后噴管中心截面速度云圖Fig.8 Velocity contour in axial-section of nozzle after optimization

5 結(jié) 論

1) 數(shù)值模擬結(jié)果與試驗(yàn)測(cè)量結(jié)果基本吻合，說(shuō)明計(jì)算模型及計(jì)算方法的可靠性較好.

2) 噴口面積、噴管直徑和噴管長(zhǎng)度3個(gè)試驗(yàn)因素中，對(duì)試驗(yàn)指標(biāo)——效率的影響程度由小到大為噴口面積、噴管直徑、噴管長(zhǎng)度.

3) 正交優(yōu)化設(shè)計(jì)能夠獲得水輪的最佳設(shè)計(jì)方案.經(jīng)數(shù)值模擬對(duì)比發(fā)現(xiàn)優(yōu)化后水輪的效率提高了7.67%.

4) 經(jīng)分析數(shù)值模擬所得水輪內(nèi)流場(chǎng)，發(fā)現(xiàn)壓力水從噴口中直接噴出的方式水力損失較大.在后續(xù)研究中建議將噴嘴改型為流線型，以期進(jìn)一步減少水力損失、提高水輪的效率.