于洪仕，張 濤，趙珊珊，蔣 偉

(1. 天津大學電氣與自動化工程學院，天津 300072；2. 遼寧工程技術(shù)大學電子與信息工程學院，葫蘆島 125105)

多孔孔板流量計流場仿真

于洪仕1,2，張 濤1，趙珊珊1，蔣 偉1

(1. 天津大學電氣與自動化工程學院，天津 300072；2. 遼寧工程技術(shù)大學電子與信息工程學院，葫蘆島 125105)

為準確計算多孔孔板流量計流場，采用Standard-kω模型、SST(剪切應力傳輸) -kω模型、Standard k-ω+ SST k-ω(Standard k-ω模型的計算結(jié)果作為SST-kω模型仿真計算的初始值)組合模式分別對100,mm口徑、節(jié)流比為0.6的3種結(jié)構(gòu)多孔孔板流量計流場進行數(shù)值計算，并結(jié)合射流理論以及實流實驗結(jié)果對數(shù)值計算結(jié)果進行分析．結(jié)果表明：對于中心節(jié)流孔與環(huán)形排列孔之間距離較小的多孔孔板，SST-kω模型收斂性較好；反之，SST-kω模型計算結(jié)果收斂困難，Standard-kω+SST-kω組合模式在保證計算精度的前提下改善了收斂效果．相對于Standard-kω模型，SST-kω模型更適合計算多孔孔板流量計的流場，計算結(jié)果符合射流理論，能反映出不同多孔孔板流出系數(shù)線性度的差異，與實流實驗結(jié)果的最大誤差為4.2%．

多孔孔板；流出系數(shù)；湍流模型；仿真

多孔孔板流量計是一個對稱的多孔圓盤，是在標準孔板基礎(chǔ)上發(fā)展起來的非標準節(jié)流裝置．2006年該流量計被引入中國市場，開始應用于天然氣、化工、煉油等工業(yè)領(lǐng)域．從相關(guān)文獻[1-3]可以看出該流量計具有比標準孔板更為出色的計量性能，如測量精度高、量程范圍寬、壓力損失小、前后直管段要求低等優(yōu)點．多孔孔板流量計結(jié)構(gòu)參數(shù)多，如節(jié)流孔的大小、個數(shù)及排列方式等，優(yōu)化結(jié)構(gòu)參數(shù)是提高多孔孔板流量計性能的前提條件．實現(xiàn)這一研究有實流實驗和數(shù)值模擬2種方法．數(shù)值模擬方法是研究流量傳感器特性的有效手段之一，既可降低成本，又可提高效率．目前，關(guān)于對多孔孔板流量計流場仿真方法的研究國內(nèi)外尚鮮見文獻報道．

因此，筆者在對多孔孔板流量計的研究過程中引入該方法，一方面可以加速研究進程，另一方面通過選擇合適的計算模型提高多孔孔板流量計流場計算的準確性．

1 湍流模型的選擇

由于目前尚無對流場具有普適性的湍流模型，科研人員只能根據(jù)流場概況選擇相對合理的湍流模型．在近幾年的研究中，-kω湍流模型被廣泛應用，吳軍等[4]采用SST-kω模型對汽車繞流仿真；EIBehery等[5]采用Standard-kω模型對二維非對稱擴散器進行仿真；Garg等[6]利用SST-kω對跨音速渦輪葉片上熱量傳遞進行仿真．上述研究取得較好的效果，這說明雙方程形式的-kω模型在計算近壁區(qū)流場、含有尾渦及剪切層等流場具有較好的計算效果．

由于多孔孔板孔分布具有分散性，流體經(jīng)過多孔孔板后在管道中形成受限性多股射流．射流自孔口出射后與周圍靜止流體間形成速度不連續(xù)的間斷面，速度間斷面是不穩(wěn)定的，必定會產(chǎn)生波動，并發(fā)展成漩渦，從而引起紊動．這樣就把原來周圍處于靜止狀態(tài)的流體卷吸到射流中，形成射流的卷吸現(xiàn)象[7]．根據(jù)文獻[7]中的雙股射流理論，流體經(jīng)過多孔孔板后多股射流間形成會聚區(qū)，最終合而為一進入聯(lián)合區(qū)．由于卷吸現(xiàn)象的存在，會聚區(qū)內(nèi)形成射流間回流區(qū)，各股射流與壁面之間產(chǎn)生近壁面回流區(qū)，在壁面回流區(qū)和射流間回流區(qū)中有大量的漩渦存在，流場如圖1所示．

圖1 多孔孔板流場示意Fig.1 Flow field of multi-hole orifice

由于射流與周圍靜止流體的卷吸與摻混，相應地產(chǎn)生了對射流的阻力，使射流邊緣部分流速降低，難以保持原來的初始流速．射流與周圍流體的摻混自邊緣部分向中心發(fā)展，經(jīng)過一定的距離發(fā)展到射流中心，自此以后射流的全斷面上都發(fā)展成湍流．由孔口邊界開始向內(nèi)外擴展的摻混區(qū)即為剪切層，因此，流體經(jīng)過多孔孔板形成的多股射流流場中存在較多的剪切層．綜上所述，多孔孔板流量計的流場情況較為復雜，這就要求湍流計算模型對含有大量漩渦及剪切層的流場具有較好的計算效果；由于多孔孔板流量計采用壁面取壓方式，該取壓方式要求湍流計算模型對近壁區(qū)域有較好的計算效果．

基于上述兩方面原因，本研究采用雙方程形式的Standard-kω模型、SST-kω模型以及Standard-kω+ SST-kω組合形式分別對10塊100,mm口徑、β= 0.6的多孔孔板進行了數(shù)值模擬與實流實驗，流速范圍為0.5～7.5,m/s．本文選擇了其中3塊具有代表性的多孔孔板對研究結(jié)果進行說明．

2 湍流模型

Standard-kω模型是一個通用雙方程湍流模型[8-9]，其中一個變量是湍動能k，另一個變量為耗散率ω．Standard-kω模型是基于Wilcox-kω模型，該模型對近壁區(qū)域的流動、尾流、射流、剪切層及低雷諾數(shù)流動有較好的預測效果．SST-kω模型是由Menter提出的雙方程湍流模型，該模型不但集成了Standard-kω模型特點與Standard -kω模型對高雷諾數(shù)流動具有較好計算效果的優(yōu)點，而且增加了橫向擴散導數(shù)項，在湍流黏度定義中考慮了湍流剪切應力的傳輸過程．其模型為

其中

3 建模網(wǎng)格剖分

3.1 多孔孔板流量計的幾何結(jié)構(gòu)

圖2為多孔孔板流量計結(jié)構(gòu)，其中圖2(a)為流量計的整體結(jié)構(gòu)，圖2(b)為多孔孔板的結(jié)構(gòu)及參數(shù)定義．圖2(b)中D為多孔孔板流量計的管徑；D1為中心節(jié)流孔直徑；D2為環(huán)狀排列孔直徑；D3為環(huán)狀排列孔的中心圓直徑；多孔孔板中心節(jié)流孔與環(huán)形排列孔之間的距離為d1，環(huán)形排列孔與管壁之間的距離為d2．圖3為多孔孔板實驗樣機，d1、d2的大小決定了射流間回流區(qū)及壁面回流區(qū)的尺寸，因此表1中給出了各樣機的d1、d2的具體數(shù)值．

圖2 多孔孔板流量計結(jié)構(gòu)Fig.2 Structure of multi-hole orifice flow meter

圖3 多孔孔板實驗樣機Fig.3 Experimental prototype of multi-hole orifice

3.2 網(wǎng)格剖分

按照流量計的實際尺寸在GAMBIT中建立三維計算模型，前直管段長度設(shè)置為10倍管徑，后直管段長度設(shè)置為30倍管徑．為了準確獲得多孔孔板附近的流場變化情況，多孔孔板附近采用size function函數(shù)進行加密處理，特別在多孔孔板的下游，加密區(qū)域更大，而在遠離多孔孔板的上下游直管段區(qū)域的網(wǎng)格逐漸變得稀疏，最密處網(wǎng)格尺寸與兩側(cè)稀疏處的比為1∶5．網(wǎng)格質(zhì)量為EquiSize Skew值為0.75，EquiAngle Skew值為0.80，Aspect Ratio值為1.0∶3.4．圖4為多孔孔板B仿真模型局部網(wǎng)格．

表1 實驗樣機的幾何參數(shù)Tab.1 Geometric parameters of experimental prototype

圖4 多孔孔板仿真模型局部網(wǎng)格Fig.4 Gridding shape in simulation of multi-hole orifice

4 計算結(jié)果分析

衡量湍流模型對節(jié)流式流量計數(shù)值計算效果優(yōu)劣標準如下.

(1)在同樣的流量范圍內(nèi)，比較數(shù)值計算得出的流出系數(shù)C與實流實驗結(jié)果是否具有一致性；

(2)通過對不同物理量的流場分析，判斷計算結(jié)果是否與相應流體力學理論一致.

4.1 流出系數(shù)C的計算結(jié)果與分析

節(jié)流式流量計測量不可壓縮流體的體積流量計算公式為

式中：qV為體積流量，m3/s；Δp為上下游取壓點測得的差壓值，Pa，在仿真實驗中，來自流場數(shù)值計算結(jié)束后壓力場數(shù)據(jù)的提取，在實流實驗中則直接來自差壓變送器的讀數(shù)；ρ為流體的密度，kg/m3；β與d分別是多孔孔板的等效直徑比和節(jié)流孔的等效直徑，在實驗中均為確定的幾何參數(shù)；C為節(jié)流式流量計的流出系數(shù)，該參數(shù)是從仿真計算或者是實流實驗中得出，因此節(jié)流式流量計的流出系數(shù)C是評價節(jié)流式儀表性能的最重要參數(shù)．

為了便于書寫，Standard-kω、SST-kω、Standard-kω+SST-kω分別采用如下縮寫形式：STD、SST、STD+SST．圖5～圖7是STD模型、SST模型及STD+SST組合形式在同一雷諾數(shù)范圍內(nèi)對不同結(jié)構(gòu)的多孔孔板流量計計算得出的流出系數(shù)C值和實流實驗值(EXP)的比較．每個湍流模型的8個仿真實驗點對應入口流速分別為0.5,m/s、1.0,m/s、2.0,m/s、3.0,m/s、4.0,m/s、5.0,m/s、6.0,m/s和7.5,m/s．

在數(shù)值計算過程中，對于多孔孔板A、B，SST模型在計算過程中發(fā)散．從圖5～圖7可以看出，在這3種數(shù)值計算方式中，SST模型或STD+SST模式計算得到的流出系數(shù)C在變化趨勢與實流實驗結(jié)果吻合得最好；STD模型計算得到的流出系數(shù)C的變化趨勢與實流實驗之間有輕微的差異，但總體趨勢一致．

圖5 多孔孔板A仿真結(jié)果與實驗結(jié)果比較Fig.5 Comparison of the results between simulations and experiments of multi-hole orifice A

圖6 多孔孔板B仿真結(jié)果與實驗結(jié)果比較Fig.6 Comparison of the results between simulations and experiments of multi-hole orifice B

圖7 多孔孔板C仿真結(jié)果與實驗結(jié)果比較Fig.7 Comparison of the results between simulations and experiments of multi-hole orifice C

表2和表3中定量地給出了采用各數(shù)值計算方法得出的計算結(jié)果．表2中定量地給出了采用各數(shù)值計算方法得到的流出系數(shù)平均值、實流實驗得出的流出系數(shù)平均值及其平均值相對誤差，該誤差定義為

表3中定量地給出了采用各數(shù)值計算方法得到流出系數(shù)線性度ELA以及實流實驗得出的流出系數(shù)線性度ELE，計算流出系數(shù)線性度的表達式為

式中：Cmax為所有流量點中流出系數(shù)最大值；Cmin為所有流量點中流出系數(shù)最小值．

表2 不同數(shù)值計算方法計算的多孔孔板流出系數(shù)結(jié)果比較Tab.2Comparisons of calculating discharge coefficient of multi-hole orifice using different models

表3 不同數(shù)值計算方法計算的多孔孔板流出系數(shù)線性度結(jié)果比較Tab.3Comparisons of calculating discharge coefficient linearity of multi-hole orifice using different models

從表2中可以看出，對于多孔孔板C，3種計算模式均收斂，STD模型計算結(jié)果的相對誤差為6.90%，SST模型與STD+SST模式計算結(jié)果的相對誤差較小，分別為4.30%與4.20%．對于多孔孔板A與B，STD模型與STD+SST模式計算結(jié)果的相對誤差均較小，其中STD+SST模式對多孔孔板計算結(jié)果的相對誤差隨著d2值的減小而減小．從表3中可以看出，利用STD+SST模式計算多孔孔板可以較好地反映出不同形式多孔孔板的流出系數(shù)線性度．

4.2 不同物理量流場分析

圖8 多孔孔板A的速度云圖(STD)Fig.8 Velocity nephogram of multi-hole orifice A(STD)

圖9 多孔孔板A的速度云圖(SST)Fig.9 Velocity nephogram of multi-hole orifice A(SST)

圖10 多孔孔板B的速度云圖(STD)Fig.10 Velocity nephogram of multi-hole orifice B(STD)

圖11 多孔孔板B的速度云圖(STD+SST)Fig.11 Velocity nephogram of multi-hole orifice B(STD+SST)

(1) 從上述分析可知，分別用STD湍流模型和STD+SST組合模式計算多孔孔板A、B得出的流出系數(shù)計算結(jié)果與實流實驗結(jié)果相對誤差均較小，但是速度場和湍流強度場卻有很大差別，如圖8～圖13所示．Standard k-ε 湍流模型對高雷諾數(shù)湍流及具有自由剪切層的湍流具有很好的計算效果，SST模型中集成了Standard k-ε 湍流模型的這一優(yōu)點，所以利用STD+SST模式仿真多孔孔板A得到的下游速度流場具有明顯的會聚趨勢，符合文獻[4]中的雙股理論，而利用STD仿真多孔孔板A得到的下游射流沒有明顯會聚趨勢．多孔孔板B的速度場云圖雖然符合射流理論，但是利用STD+SST模式計算的湍流強度場中湍流強度最大的位置在射流的剪切層中，與文獻[10]結(jié)論一致．因此可以看出SST湍流模型比STD湍流模型更適合計算受限性多股射流相互作用的流場．

(2) 從圖9、圖11和圖14中可以看出，相對于多孔孔板C，多孔孔板A、B的射流間回流區(qū)域較大，壁面回流區(qū)域較?。苯邮褂肧ST模型計算射流間回流區(qū)域較大多孔孔板時的收斂比較困難，而STD+SST組合模式不但克服了上述缺點并且計算效果較好．

圖12 多孔孔板B的湍流強度云圖(STD)Fig.12 Turbulence intensity nephogram of multi-hole orifice B(STD)

圖13 多孔孔板B的湍流強度云圖(STD+SST)Fig.13 Turbulence intensity nephogram of multi-hole orifice B(STD+SST)

圖14 多孔孔板C的速度云圖(SST)Fig.14 Velocity nephogram of multi-hole orifice C(SST)

(3)如前文所述，SST模型在近壁區(qū)以外及剪切層中集成了Standard k-ε 湍流模型的特點，而Standard k-ε 湍流模型本身存在缺陷，該模型在彎曲壁面、彎曲流線等情況下會產(chǎn)生失真．多孔孔板A、B、C的壁面回流區(qū)依次增大，所以采用壁面取壓方式時，計算得出流出系數(shù)平均值與實流實驗得出的流出系數(shù)平均值之間的相對誤差依次減?。?/p>

5 結(jié) 語

通過有限體積法數(shù)值求解Reynolds平均N-S方程，湍流模型分別用STD模型、SST模型及STD+ SST組合模式對3塊多孔孔板流量計進行了數(shù)值模擬．結(jié)果表明：對于中心節(jié)流孔與環(huán)形排列孔之間距離較小的多孔孔板，SST模型收斂性較好；對于中心節(jié)流孔與環(huán)形排列的小孔之間距離較大的多孔孔板，SST模型計算結(jié)果收斂困難，STD+SST組合模式在保證計算精度的前提下改善了收斂效果．相對STD模型，SST模型更加適合計算多孔孔板流量計的內(nèi)部流場，計算結(jié)果與射流力學中的雙股射流理論一致，與實流實驗結(jié)果誤差的最大值為4.2%，并且能反映出不同多孔孔板流出系數(shù)線性度的差異．因此，利用該方法計算多孔孔板流場對優(yōu)化多孔孔板結(jié)構(gòu)具有一定的指導意義，并且對其他具有射流性質(zhì)的流場仿真具有一定的參考價值．

［1］ 周 人. A+K平衡流量計的獨特性能和典型應用[J].醫(yī)藥工程設(shè)計，2011，32(3)：54-57.

Zhou Ren. Unique characteristic and typical application of A+K balanced flow meter[J]. Pharmaceutical and Engineering Design，2011，32(3)：54-57(in Chinese).

［2］ 馬太義，王 棟，張炳東，等. 多孔板流量測量的實驗研究[J]. 核動力工程，2010，31(2)：126-130.

Ma Taiyi，Wang Dong，Zhang Bingdong，et al. Experimental on metering characteristics of multi-hole orifice[J]. Nuclear Power Engineering，2010，31(2)：126-130(in Chinese).

［3］ Zhao Tianyi，Zhang Jili. Experimental investigation of key parameters pertinent to multi-hole orifice throttling characteristic[J]. Journal of Harbin Institute of Technology：New Series，2008，15(1)：9-12.

［4］ 吳 軍，谷正氣，鐘志華. SST湍流模型在汽車繞流仿真中的應用[J]. 汽車工程，2003，25(4)：326-329.

Wu Jun，Gu Zhengqi，Zhong Zhihua. The application of SST turbulence model in the aerodynamic simulation of the automobile[J]. Automobile Engineering，2003，25(4)：326-329(in Chinese).

［5］ EI-Behery S M，Hamed M H. A comparative study of turbulence models performance for separating flow in a planar asymmetric diffuser[J]. Computers and Fluids，2011，44(1)：248-257.

［6］ Gary V K，Ameri A A. Two-equation turbulence models for prediction of heat transfer on a transonic turbine blade[J]. Heat and Fluid Flow，2001，22(6)：593-602.

［7］ 董志勇. 射流力學[M]. 北京：科學出版社，2005.

Dong Zhiyong. Jet Mechanics[M]. Beijing：Science Press，2005(in Chinese).

［8］ Bakabel A，Hegab A M，Nasr M，et al. Assessment of turbulence modeling for gas flow in two dimensional convergent-divergent rocket nozzle[J]. Applied Mathematical Modelling，2011，35：3408-3422.

［9］ Qiu Yunming，Liu Zuyuan，Chen Xide，et al. Numerical calculation of maneuvering hydrodynamic forces of drift ship based on SST κ-ω turbulence model[C] // The 2nd International Conference on Computer Engineering and Technology. Chengdu，China，2010：285-288.

［10］ Fujisawa N，Nakamura K，Srinivas K. Interaction of two parallel plane jets of different velocities[J]. Journal of Visualization，2004，7(2)：135-142.

(責任編輯：孫立華)

Simulation of the Flow Field of Multi-Hole Orifice Flow Meter

Yu Hongshi1,2，Zhang Tao1，Zhao Shanshan1，Jiang Wei1
(1. School of Electrical Engineering and Automation，Tianjin University，Tianjin 300072，China；2. School of Electronic and Information Engineering，Liaoning Technical University，Huludao 125105，China)

In order to calculate the flow field of multi-hole orifice flow meter accurately，Standard k-ω turbulence model，SST(shear-stress transport)k-ω turbulence model and the composite pattern of Standard k-ω and SST kω (the calculation result of Standard k-ω as the original value of SST k-ω) were applied to simulate the flow field of three multi-hole orifice flow meters，each with the inside diameter of 100,mm and the equivalent diameter ratio β of 0.6，but with different structures，and the calculation result was analyzed according to the real physical experiment data and theory of jet. The result shows：for the multi-hole orifice with a short distance between the center hole and the holes arrayed on the ring，the convergence of SST k-ω model is better；oppositely，the convergence of SST kω model is difficult，and the composite pattern of Standard k-ω and SST k-ω enhance effect of convergence on the premise of calculation accuracy. Relative to Standard k-ω，the SST k-ω model is more accurate and performs better in calculating flow field of multi-hole orifice，the calculation result conforms to the theory of jet and shows the difference of linearity of different multi-hole orifices，and the maximum deviation between calculation result and real physical experiment result is 4.2%.

multi-hole orifice；discharge coefficient；turbulence model；simulation

TH814

A

0493-2137(2014)01-0061-06

10.11784/tdxbz201207058

2012-07-21；

2012-10-22.

于洪仕（1980— ），男，博士研究生，yuhongke_2004@163.com.

張 濤，zt50@tju.edu.cn.