惠 偉，王少飛

● （上海交通大學(xué) 動(dòng)力機(jī)械與工程教育部重點(diǎn)實(shí)驗(yàn)室，上海 200240）

閥門角度對蝶閥后雙彎管流場影響的數(shù)值模擬研究

惠 偉，王少飛

● （上海交通大學(xué) 動(dòng)力機(jī)械與工程教育部重點(diǎn)實(shí)驗(yàn)室，上海 200240）

本文采用計(jì)算流體力學(xué)（CFD）數(shù)值模擬的方法，研究了不同閥門角度下蝶閥后雙彎管模型中的復(fù)雜流動(dòng)現(xiàn)象，并將數(shù)值計(jì)算的速度云圖同流場實(shí)驗(yàn)測量結(jié)果進(jìn)行對比驗(yàn)證。分析結(jié)果表明：數(shù)值計(jì)算結(jié)果同粒子圖像速度場測量技術(shù)得到的速度云圖基本吻合；閥門角度對流場的影響較大，閥門角度越小，閥板迎背流面的流體擴(kuò)張也越不明顯，模型最大速度減??；閥板前駐點(diǎn)逐漸向閥板迎流面邊側(cè)移動(dòng)，一次分離區(qū)減小，而二次分離區(qū)先增大后減?。粡澒苤辛黧w質(zhì)點(diǎn)二次流、流動(dòng)分離及流動(dòng)剪切膨脹等是影響流場的重要因素。

計(jì)算流體力學(xué)；彎管；蝶閥；粒子圖像速度場測量；流動(dòng)分離

0 引言

由于結(jié)構(gòu)簡單、操作方便等原因，蝶閥作為流動(dòng)控制與調(diào)節(jié)的設(shè)備被廣泛應(yīng)用于動(dòng)力工程、城市供水、化工、船舶等工業(yè)場合[1-3]。然而，在調(diào)節(jié)蝶閥閥板角度以實(shí)施流動(dòng)控制時(shí)，蝶閥閥板下游常常有較大流動(dòng)分離，容易誘發(fā)噪聲輻射和結(jié)構(gòu)振動(dòng)等不利現(xiàn)象。布置在長直管內(nèi)部的蝶閥流動(dòng)分離現(xiàn)象已經(jīng)得到了非常廣泛的重視和研究[7-9,11]。然而，在某些特殊場合，由于空間布置限制等原因[3]，蝶閥常常需要和彎管就近配合使用，這使得蝶閥和彎管系統(tǒng)內(nèi)部的流動(dòng)現(xiàn)象更加復(fù)雜[4]。因此，很有必要針對蝶閥彎管系統(tǒng)內(nèi)部的復(fù)雜流場現(xiàn)象進(jìn)行系統(tǒng)的研究。

本文所選擇的蝶閥后雙彎管模型是對某電廠汽輪機(jī)系統(tǒng)中低壓缸間抽汽管道進(jìn)行的一定幾何相似縮比模型，如圖1所示。圖1中管道橫截面直徑d=100mm，閥門角度定義為來流與閥板迎流面夾角（銳角）θ。其中閥板板面由迎背流面雙凸型圓弧與前后緣扇形圓弧組成。本文采用計(jì)算流體力學(xué)手段對這一蝶閥后雙彎管系統(tǒng)中所存在的復(fù)雜流動(dòng)現(xiàn)象進(jìn)行了系統(tǒng)的數(shù)值模擬和分析，與流場測量結(jié)果進(jìn)行對比驗(yàn)證后，重點(diǎn)對比分析相同進(jìn)口質(zhì)量流量下不同閥門角度（0°、15°及30°）對模型復(fù)雜流場的影響。

圖1 模型示意圖

圖2 網(wǎng)格劃分策略

1 數(shù)值計(jì)算模型

本文采用商用計(jì)算流體力學(xué)軟件 ANSYS CFX 14.0實(shí)現(xiàn)流場穩(wěn)態(tài)計(jì)算。由于幾何模型結(jié)構(gòu)的對稱性，并基于相關(guān)比較分析及文獻(xiàn)[9]，實(shí)際計(jì)算區(qū)域僅取原始模型的一半，以避免計(jì)算網(wǎng)格數(shù)量過多帶來的諸多不利。采用ICEM CFD 14.0實(shí)施混合網(wǎng)格劃分，閥門角度為0°時(shí)網(wǎng)格劃分策略見圖2，其它兩種角度與之相似，出于篇幅考慮，本文不再附圖贅述。湍流模型采用該軟件中SST模型[6]，壁面函數(shù)采用automatic near-wall treatment（近壁面網(wǎng)格滿足Y+<2）。

流動(dòng)控制方程組由可壓縮流動(dòng)雷諾方程、理想氣體狀態(tài)方程和各湍流模型方程組成[10,11]?？諝饨橘|(zhì)按理想氣體處理，參考壓力為97856Pa，并考慮粘性耗散效應(yīng)。離散方程采用有界高精度對流格式求解，并通過Rhie和Chow算法計(jì)算質(zhì)量流量，以保證壓力速度耦合。在299k和湍流度 5%下，進(jìn)口給定總壓（34.85Pa）；出口給定質(zhì)量流量（m=0.0463kg/s），且速度和溫度分布滿足局部單向化假設(shè)；閥板用丙烯黑色顏料涂黑處理，故設(shè)置為無滑移絕熱邊界條件；因?qū)嶒?yàn)管道材料（有機(jī)玻璃）導(dǎo)熱系數(shù)較小，故管道壁面設(shè)置為無滑移絕熱邊界條件；對稱面（圖1中平行紙面方向所取的中截面）設(shè)置為對稱邊界條件。分別在閥板迎背流面及兩彎管內(nèi)側(cè)各選定一個(gè)速度監(jiān)控點(diǎn)，殘差標(biāo)準(zhǔn)設(shè)為10-5。網(wǎng)格無關(guān)性驗(yàn)證指標(biāo)是閥板扭矩值是否隨網(wǎng)格數(shù)增加而變化，見圖4。此處扭矩表示流體對指定二維區(qū)域（這里指閥板面）作用力相對特定軸（這里指穿過閥板中心且垂直紙面的軸線）產(chǎn)生的力矩；對壁面邊界條件而言，作用力為氣體壓力和粘性切應(yīng)力之和在所選面上積分[7]；具體計(jì)算過程通過商業(yè)軟件中相關(guān)函數(shù)實(shí)現(xiàn)。

2 PIV試驗(yàn)驗(yàn)證

2.1 試驗(yàn)設(shè)備

為驗(yàn)證數(shù)值計(jì)算結(jié)果，本文首先進(jìn)行了流場測量實(shí)驗(yàn)，具體是通過高分辨率粒子圖像速度場測量技術(shù)（Particale Image Velocimetry,PIV）[8]來完成的，其可獲得閥門角度為30°時(shí)中截面閥板附近的速度云圖。該測量系統(tǒng)由高分辨率相機(jī)、激光器、同步控制機(jī)構(gòu)、圖像采集系統(tǒng)以及PIV圖像分析軟件組成?；驹砣缦拢杭す馄骱拖鄼C(jī)通過同步器控制，采集的數(shù)字圖像實(shí)時(shí)傳輸?shù)接?jì)算機(jī)硬盤中，經(jīng)過MicroVec粒子圖像分析軟件處理，得到瞬態(tài)流場整場速度分布。

實(shí)驗(yàn)中使用的激光發(fā)射裝置為英國 Litro公司 Nano系列YANG式脈沖激光器，能量135mJ/pluse，激光波長532nm，發(fā)射激光的脈沖寬度 8ns；高分辨率相機(jī)為為美國 IMPERX公司的高分辨率，黑白工業(yè)相機(jī)，型號為ICL-B4820M- KF000，分辨率可達(dá)4872×3248；同步控制器為MicroPulse725可編程同步時(shí)序控制器，集成有時(shí)間信號發(fā)生器，具有高精度獨(dú)立的7通道延時(shí)時(shí)間控制器，延時(shí)精度為0.25ns；示蹤粒子為癸二酸二辛酯（DEHS）；圖像采集和計(jì)算采用北京立方天地MicroVec V3.1.1軟件完成。

為了達(dá)到最好的拍攝效果以剔除背景噪聲，在成像前用丙烯黑色顏料對管道和蝶閥背面進(jìn)行涂黑處理，并在實(shí)驗(yàn)進(jìn)行過程中嚴(yán)禁除激光以外的光線干擾。另外在進(jìn)行數(shù)據(jù)處理之前，通過 matlab軟件對圖像進(jìn)行預(yù)處理，以進(jìn)一步剔除掉背景噪聲。

2.2 結(jié)果對比分析

PIV技術(shù)獲得的閥門角度為30°時(shí)的速度云圖如圖3中(a)圖所示。從圖可知閥板尾流區(qū)形成較大分離區(qū)（低速區(qū)），這是流體在閥板背流面流動(dòng)擴(kuò)張及在彎管向心運(yùn)動(dòng)造成的逆壓力梯度引起的。此時(shí)，流體經(jīng)過閥板迎流面與管道外側(cè)間喉部（最小通流面）后，由于喉部面積較小而速度較高，因而表現(xiàn)為高速（紅色）流體區(qū)域；但由于彎管外側(cè)對流動(dòng)阻礙，此處高速流體與尾流區(qū)低速流體的混合及流體在彎管中的向心運(yùn)動(dòng)等原因，該股高速（紅色）流動(dòng)區(qū)域的發(fā)展受到限制，表現(xiàn)為高速（紅色）區(qū)域面積很快消失。此時(shí)，經(jīng)過閥板背流面與管道內(nèi)側(cè)間喉部后，由于流動(dòng)空間擴(kuò)張較大及彎管內(nèi)側(cè)流體質(zhì)點(diǎn)的向心加速而形成較大的高速（紅色）區(qū)域；但進(jìn)一步的流動(dòng)擴(kuò)張作用及其與周圍低速流體的混合，此處流體速度逐漸降低。另外，由于背流面對流體向彎管外側(cè)的引導(dǎo)，彎管內(nèi)側(cè)流動(dòng)分離現(xiàn)象較為顯著。通過與數(shù)值計(jì)算結(jié)果(b)比較可知，該角度下實(shí)驗(yàn)測量和數(shù)值計(jì)算得到的速度云圖基本吻合。

圖3 中截面閥門角度為30°速度云圖對比

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

基于實(shí)驗(yàn)測量得到的閥門角度為30°時(shí)中截面閥板附近的速度云圖很好地驗(yàn)證了數(shù)值計(jì)算結(jié)果，本部分將數(shù)值計(jì)算得到的結(jié)果做進(jìn)一步分析，以了解各工況（閥門角度為0°，15°及30°時(shí)）中存在的復(fù)雜流動(dòng)現(xiàn)象。

本部分采用中截面速度矢量及云圖與五個(gè)剖面的速度矢量及云圖進(jìn)行三維流場分析。各剖面位置示意圖如圖4所示，其中α=arctan0.2，Slice-2和Slice-4分別位于彎管垂直紙面的對稱面上，其它截面位置如圖4所示。另外，本部分借助中截面流線圖進(jìn)行閥板周圍流動(dòng)分離現(xiàn)象分析。

圖4 五個(gè)剖面位置示意圖

3.1 閥門角度為30°工況

結(jié)合流線圖與中截面速度及矢量圖可知，此時(shí)閥板背流面形成兩個(gè)流動(dòng)分離區(qū)，第一分離區(qū)(L1)較大，這是因?yàn)閬砹鞴ソ禽^大而產(chǎn)生了較強(qiáng)的逆壓力梯度；第一分離區(qū)的存在有助于引導(dǎo)流體質(zhì)點(diǎn)在背流面中央到后緣的大范圍內(nèi)緊貼背流面流動(dòng)，因而造成第二分離區(qū)(L2)較??；同時(shí)，第二分離區(qū)的減小使閥板后緣扇形圓弧區(qū)的流動(dòng)分離大大削弱。

圖5為閥門角度為30°中截面速度矢量及云圖與流線圖，圖6為閥門角度為30°各剖面速度矢量及云圖。綜合Slice-1與圖5可知，閥板迎流面與管道外側(cè)間流體由于向彎管內(nèi)側(cè)擴(kuò)張導(dǎo)致剖面圖上部流體速度矢量向下；下部流體質(zhì)點(diǎn)以較高速向管道中間大面積擴(kuò)張導(dǎo)致速度矢量向上，且該部分流體質(zhì)點(diǎn)由于喉部通流面積較小因而速度較高，而喉部后通流面積擴(kuò)張較大，而形成面積較大的高速區(qū)，這與上部分速度云圖中的結(jié)論一致。同時(shí)，由于下部流體通流面積擴(kuò)張較大導(dǎo)致流體質(zhì)點(diǎn)速度較低、動(dòng)量較小，同上部流向相反的流團(tuán)碰撞后，會(huì)向管道兩旁流動(dòng)，因而形成一對低速漩渦（左側(cè)逆時(shí)針，右側(cè)時(shí)針），該低速流通同時(shí)與下部偏管道內(nèi)側(cè)的高速流團(tuán)發(fā)生流動(dòng)剪切，并接受高速流團(tuán)動(dòng)量和流量“補(bǔ)給”，該低速流團(tuán)也是蝶閥對流動(dòng)阻礙作用的體現(xiàn)。

圖5 閥門角度為30°中截面速度矢量及云圖與流線圖

圖6 閥門角度為30°各剖面速度矢量及云圖

從Slice-2剖面圖可知，在下部低速流團(tuán)碰撞剪切以及彎管對流體二次流的作用下，Slice-1中心線具有向下的速度矢量的上部流團(tuán)消失，只存在向管道兩側(cè)的運(yùn)動(dòng)，如該剖面中部兩側(cè)的紅色高速區(qū)域所示；而在上部高速流團(tuán)以及下部高速流團(tuán)的碰撞剪切下，低速流團(tuán)的漩渦核心向下運(yùn)動(dòng)并不斷縮小；同時(shí)可以看到，下側(cè)流團(tuán)的高速區(qū)域僅僅表現(xiàn)在中心線兩側(cè)周圍區(qū)域，這是slice-1剖面中低速漩渦流團(tuán)流動(dòng)剪切反作用的結(jié)果。

從中截面速度矢量及云圖可知，Slice-2到Slice-3間管道內(nèi)側(cè)與閥板背流面間的高速流團(tuán)由于二次流作用向外側(cè)運(yùn)動(dòng)，同時(shí)在直管內(nèi)側(cè)由于流動(dòng)分離作用出現(xiàn)低速區(qū)域。

管道內(nèi)側(cè)流動(dòng)分離造成的低速區(qū)體現(xiàn)為Slice-3剖面下方速度矢量向上的的藍(lán)色區(qū)域，而從中截面矢量及云圖中表現(xiàn)為向管道外側(cè)及后方的速度矢量；該剖面中部兩旁存在高速區(qū)域，這與Slice-2中在該區(qū)域表現(xiàn)相同，但渦團(tuán)核心位置有所下降，這是其二次流螺旋運(yùn)動(dòng)的結(jié)果。而Slice-2中心線上的高速流團(tuán)（管道內(nèi)側(cè)與閥板背流面間流體）由于流動(dòng)擴(kuò)張以及碰撞管道外側(cè)的能量損失（中截面矢量及云圖也可以看出），速度降低并向剖面左右兩側(cè)運(yùn)動(dòng)，其與剖面中部兩旁高速區(qū)域相混。

通過Slice-4剖面圖以及中截面速度-矢量剖面可知，由于Slice-3剖面圖左右兩側(cè)的高速流團(tuán)二次流螺旋運(yùn)動(dòng)到此處與原其管道內(nèi)側(cè)流體質(zhì)點(diǎn)混合，共同在第二管道內(nèi)側(cè)向心加速，形成Slice-4中下方的高速區(qū)域；而在高速區(qū)域上方存在一對低速流團(tuán)，這是管道內(nèi)側(cè)與閥板背流面間流團(tuán)降速后二次流運(yùn)動(dòng)到此處的表現(xiàn)。

通過剖面圖Slice-5可知，由于第二個(gè)彎管內(nèi)側(cè)的流動(dòng)分離，剖面下部流體質(zhì)點(diǎn)具有向上的速度矢量且大小較低；上部則是混合流體二次流作用的結(jié)果，且其速度高于下方。

3.2 閥門角度15°工況

圖7為閥門角度為15°中截面速度矢量及云圖與流線圖，圖8為閥門角度為15°各剖面速度矢量及云圖。閥門角度15°時(shí)，綜合中截面流線圖和速度矢量及云圖可知，閥板背流面上第一分離區(qū)(L1)變小而第二分離區(qū)增大，且后者表現(xiàn)為一對逆向渦團(tuán)。第一分離區(qū)較小是由于來流攻角較小，流線型閥板對來流具有很好引導(dǎo)作用，因而產(chǎn)生較小的逆壓力梯度。第二分離區(qū)增大，則是由于第一分離區(qū)的存在以及彎管對流體質(zhì)點(diǎn)向內(nèi)側(cè)的引導(dǎo)，因而流線沿背流面運(yùn)動(dòng)至閥板中央附近即發(fā)生分離；之所以形成兩個(gè)渦團(tuán)，這是因?yàn)殚y板迎流面流體在閥板尾緣由于扇形圓弧存在導(dǎo)致流動(dòng)向內(nèi)側(cè)的擴(kuò)張，產(chǎn)生新渦團(tuán)，并與前述渦團(tuán)相互作用。

由于閥板迎背流面流體流出閥板時(shí)向閥板側(cè)的擴(kuò)張，因而在剖面圖Slice-1中上部速度矢量向下而下部向上，同時(shí)中部存由于閥板流動(dòng)分離導(dǎo)致的低速區(qū)域。在該分離區(qū)的下方兩側(cè)存在一個(gè)微小渦團(tuán)，這是閥板與管道間隙（見圖1中A-A向視圖，最小處僅為1mm）通過流體的表現(xiàn)，而在閥門角度30°時(shí)，由于閥板后緣造成的低速區(qū)域過大，該渦團(tuán)與其混合表現(xiàn)的不是很明顯。通過剖面圖Slice-2可知，閥板后分離區(qū)造成的藍(lán)色（低速）區(qū)域由于分離區(qū)影響變?nèi)醵鴾p小，剖面圖上方流體質(zhì)點(diǎn)向上運(yùn)動(dòng)則是其向管道外側(cè)離心運(yùn)動(dòng)的表現(xiàn)，而向下運(yùn)動(dòng)則是由于管道內(nèi)側(cè)流體質(zhì)點(diǎn)向心的表現(xiàn)，而導(dǎo)致向心離心主要與該處壓力梯度及所需的向心力的大小有關(guān)[10]；同時(shí)可知，由于管道內(nèi)側(cè)與閥板背流面間流體的擴(kuò)張，剖面下部高速區(qū)域速度降低。

圖7 閥門角度為15°中截面速度矢量及云圖與流線圖

圖8 閥門角度為15°各剖面速度矢量及云圖

從中截面速度矢量及云圖可知，Slice-2到Slice-3間存在流體質(zhì)點(diǎn)在管道內(nèi)側(cè)的流動(dòng)分離，其在Slice-3中表現(xiàn)為中部低速流團(tuán)。而分離形成的藍(lán)色（低速）區(qū)域下方，則存在具有較高速度（大約12m/s）的流體質(zhì)點(diǎn)，其是流體質(zhì)點(diǎn)二次流的結(jié)果。同時(shí)由于分離區(qū)的擠壓，Slice-2上部的低速區(qū)域加速并在Slice-3中表現(xiàn)為速度較高的向上流團(tuán)。而分離渦團(tuán)兩旁則的高速流圖則是Slice-2中下部高速區(qū)域二次流運(yùn)動(dòng)到此處的表現(xiàn)。

在剖面圖 Slice-4中，由于管道內(nèi)側(cè)向心作用而導(dǎo)致流體質(zhì)點(diǎn)加速，中間偏上部流體質(zhì)點(diǎn)的壓力梯度大于所需向心力，流體質(zhì)點(diǎn)向管道內(nèi)側(cè)運(yùn)動(dòng)并減速；上部則相反，速度矢量方向向上；可見由于彎管的存在，不同流體質(zhì)點(diǎn)的速度，所處位置曲率以及壓力梯度的不同，導(dǎo)致徑向壓力梯度和所需向心力的不同，繼而引起不同的流動(dòng)形態(tài)。

剖面圖Slice-4中的流體質(zhì)點(diǎn)在剖面圖Slice-5中發(fā)展并表現(xiàn)為兩對渦團(tuán)，這是二次流以及流動(dòng)碰撞剪切的結(jié)果，上部速度弱微的提高也是下方分離區(qū)流體質(zhì)點(diǎn)擠壓的表現(xiàn)。

3.3 閥門角度0°工況

圖9為閥門角度為0°中截面速度矢量及云圖與流線圖，圖10為閥門角度為0°下各剖面矢量圖。閥門角度為0°時(shí)，右側(cè)閥板與管道內(nèi)側(cè)間存在一個(gè)明顯的分離區(qū)(L2)，且該分離區(qū)較?。贿@是閥板型線導(dǎo)致流動(dòng)擴(kuò)壓較為緩慢導(dǎo)致逆壓力梯度較小，但由于彎管對流體的向心作用加大了該處的逆壓力梯度，因而在右側(cè)閥板擴(kuò)張段發(fā)生了微弱的流動(dòng)分離。而左側(cè)閥板擴(kuò)張段幾乎沒有發(fā)生分離，這是由于彎管存在導(dǎo)致該側(cè)流通面積變化不大，逆壓力梯度較小，因而流動(dòng)分離不明顯。但左側(cè)閥板與管道外側(cè)間通過流體在閥板后緣由于圓弧存在導(dǎo)致突然的擴(kuò)張，因而產(chǎn)生很小的渦團(tuán)，并與右側(cè)渦團(tuán)相互作用。

圖9 閥門角度為0°中截面速度矢量及云圖與流線圖

圖10 閥門角度為0°下各剖面矢量圖

閥門角度為0°時(shí)，剖面圖Slice-1速度矢量與15°時(shí)的不同之處在以下兩點(diǎn)：第一，流動(dòng)低速區(qū)（閥板阻流形成低速區(qū)）更小，這是來流攻角減小的表現(xiàn)；第二，中心線下部存在向下的速度矢量，這是彎管內(nèi)側(cè)向心作用較大的體現(xiàn)，15°時(shí)的微弱渦團(tuán)在此時(shí)與低速區(qū)流團(tuán)混合。其它同15°時(shí)。

剖面圖Slice-2中，由于彎管內(nèi)側(cè)流體質(zhì)點(diǎn)的向心作用，Slice-1中向下的速度矢量變?yōu)橄蛏希煌瑫r(shí)由于與上下高速流團(tuán)的混合，流動(dòng)低速區(qū)面積減小，剖面上部和下部的速度降低。

從中截面速度矢量及云圖可以看出，由于第一彎管內(nèi)側(cè)流動(dòng)分離以及二次流，分離區(qū)下方速度均具有向上的分量，因而剖面圖Slice-3中截面速度均向上；而該剖面圖中存在兩對渦團(tuán)，這是兩股流團(tuán)二次流的表現(xiàn)，高速流團(tuán)是右側(cè)閥板與管道內(nèi)側(cè)間流體二次流的表現(xiàn)，而低速流團(tuán)則是左側(cè)閥板與外側(cè)管道間流體二次流的表現(xiàn)。

剖面圖 Slice-4中速度矢量與閥門角度 15°時(shí)較為相似，此處不再贅言。而閥門角度0°時(shí)第二管道內(nèi)側(cè)的流體速度相對較高，則是因?yàn)檫@是第一管道內(nèi)側(cè)的流體在能量損失較小的情況下二次流運(yùn)動(dòng)到此處時(shí)，由于速度相對較高因而再次加速而獲得更高速度。

剖面圖 Slice-5中，由于第二管道內(nèi)側(cè)流動(dòng)分離，Slice-4中下部高速區(qū)域在Slice-5中上移，而在后者下部出現(xiàn)分離（低速）區(qū)域，并與上方高速區(qū)域混合而使其減速，這兩部分流團(tuán)速度矢量均向上；同時(shí)，在整個(gè)剖面中心線上部，由于下部流體質(zhì)點(diǎn)的擠壓，Slice-4中上部渦團(tuán)核心繼續(xù)上移并加速；而Slice-4中下部渦團(tuán)由于下方高速區(qū)域的流動(dòng)碰撞擠壓而在Slice-5中消失，但Slice-5中心線上由于兩股流團(tuán)流向相反，因而也會(huì)發(fā)生流動(dòng)碰撞，并在Slice-5分離區(qū)兩旁形成一對不是很明顯的小渦團(tuán)。

比較各閥門角度下中截流線圖可知，當(dāng)閥門角度從30°變?yōu)?°時(shí)，前駐點(diǎn)(點(diǎn)p)的位置向閥板迎流面邊側(cè)移動(dòng)。這與[9]中得到的結(jié)論是一致的。

通過以上三種角度的對比分析可知，閥板對流場的影響較大，隨閥板角度減小，閥板迎背流面的流體擴(kuò)張也越不明顯，模型最大速度減??；閥板前駐點(diǎn)逐漸向閥板迎流面邊側(cè)移動(dòng)，一次分離區(qū)減小，而二次分離區(qū)先增大后減??；彎管中流體質(zhì)點(diǎn)二次流、流動(dòng)分離及流動(dòng)剪切碰撞等是影響流場的重要因素。

4 結(jié)論

本文采用計(jì)算流體力學(xué)數(shù)值模擬方法，研究了閥門角度變化(0°、15°及30°)對蝶閥后雙彎管模型中復(fù)雜流動(dòng)現(xiàn)象的影響，并將閥板角度為30°時(shí)的速度場云圖與實(shí)驗(yàn)測量結(jié)果進(jìn)行對比驗(yàn)證。發(fā)現(xiàn)閥門角度對流場的影響較大，且隨閥門角度的減?。?/p>

1）第一個(gè)彎管內(nèi)側(cè)的流動(dòng)分離現(xiàn)象越不明顯，閥板迎背流面的流體擴(kuò)張也越不明顯，模型最大速度減?。?/p>

2）閥板前駐點(diǎn)逐漸向閥板迎流面邊側(cè)移動(dòng)，一次分離區(qū)減小，而二次分離區(qū)先增大后減??；

3）彎管中流體質(zhì)點(diǎn)二次流運(yùn)動(dòng)、流動(dòng)分離及流動(dòng)剪切膨脹等是影響流場的重要因素。

[1]張新萍,鄭莉玲.長供水管道水泵出口閥門的選擇[J].電網(wǎng)與清潔能源,2009,25(4):71-73.

[2]熊衛(wèi),嚴(yán)啟明.淺析閥門在化工生產(chǎn)中的應(yīng)用[J].瀘天化科技,2009(3):266-270.

[3]魏志.閥體后 90°圓形彎管內(nèi)流場和噪聲的數(shù)值模擬[D].上海:上海交通大學(xué),2013.

[4]Danbon F,Solliec C.Aerodynamic torque of a butterfly valve:Influence of an elbow on the time-mean and instantaneous aerodynamic torque[J].Journal of fluids engineering,2000,122(2):337-344.

[5]Weerachai Chaiworapuek,et al.An Investigation of the Water Flow Past the Butterfly Valve.AIP Conf.Proc.1225,562 (2010); doi:10.1063/1.3464904.

[6]ANSYS,2011,“CFX User’s Guide”,Ver.14.0,ANSYS,Inc

[7]陶文銓.數(shù)值傳熱學(xué)(第二版)[M].西安:西安交通大學(xué)出版社,2001:2-4.

[8]萬津津.近壁方柱繞流湍流場非定常特性的 PIV實(shí)驗(yàn)研究[D].上海:上海交通大學(xué),2009.

[9]靳允立,胡俊.翼型失速及雷諾數(shù)變化對風(fēng)力機(jī)氣動(dòng)性能影響的數(shù)值研究[J].太陽能學(xué)報(bào),2009,30(9).

[10]湛含輝.二次流原理第1版[M].長沙:中南大學(xué)出版社,2006.

Numerical Study on the Influence of Valve Inclination Angle on Flow Field of Butterfly Valve Placed in a Dual-elbow Channel

HUI Wei,WANG Shao-fei
(Key Laboratory for Power Machinery and Engineering of Ministry of Education,Shanghai Jiao Tong University,Shanghai 200240,China)

By using computational fluid dynamics (CFD),the influence of different valve inclination angle on the complex flow phenomena is investigated in a butterfly valve placed in a dual-elbow channel model,and the velocity contour result is validated by experimental result.Numerical result fits better with data obtained by PIV.Valve inclination angle plays a significant role on the complex flow field.While it decreases,flow expansion in the upstream and downstream face become not obvious and the maximum velocity decreases.The front stagnation point moves gradually toward the upstream face edge.The first separation region decreases while the second one first increases and then decreases.In all,the second flow,flow separation and flow shear & collision are key factors in the complex flow formation.

computational fluid dynamics (CFD); elbow; butterfly valve; particle image velocimetry (PIV)

TP3

A

惠偉（1987-），男，碩士研究生。研究方向：閥門管道流動(dòng)傳熱數(shù)值模擬。