牛志娟，趙　奇，楊雪峰

（四川大學(xué)化學(xué)工程學(xué)院，四川成都610065）

?

環(huán)槽流量計(jì)的數(shù)值模擬與優(yōu)化

牛志娟，趙奇，楊雪峰

（四川大學(xué)化學(xué)工程學(xué)院，四川成都610065）

摘要：研究流量計(jì)內(nèi)部流場和結(jié)構(gòu)優(yōu)化，對(duì)改善流量計(jì)的測(cè)量性能和提高測(cè)量精度，具有重要的現(xiàn)實(shí)意義。將計(jì)算流體力學(xué)（CFD）仿真試驗(yàn)應(yīng)用于一種新型差壓流量計(jì)——環(huán)槽流量計(jì)，考查不同等效直徑比β、前端和尾部長度、等直徑段長度以及雷諾數(shù)對(duì)環(huán)槽流量計(jì)的流出系數(shù)和壓力損失的影響。結(jié)果表明：隨著雷諾數(shù)的增加，流出系數(shù)逐漸增大并達(dá)到穩(wěn)定值；隨著β增大，流出系數(shù)先增大后減?。磺岸思拔膊块L度對(duì)流出系數(shù)影響不大，但尾部長度越大，永久壓損越??；等直徑段長度越小，永久壓損越小。根據(jù)結(jié)果擬合出環(huán)槽流量計(jì)流出系數(shù)的公式，CFD數(shù)值模擬作為一種輔助設(shè)計(jì)和標(biāo)定手段，有助于指導(dǎo)環(huán)槽流量計(jì)的現(xiàn)場測(cè)試。

關(guān)鍵詞：環(huán)槽流量計(jì)；壓力損失；流出系數(shù)；數(shù)值模擬；測(cè)試

0　引言

差壓流量計(jì)是使用面最大最廣的一種流量計(jì)，大約占我國流量計(jì)使用量的60%。以標(biāo)準(zhǔn)孔板和噴嘴為代表的傳統(tǒng)的節(jié)流式差壓流量計(jì)有以下缺點(diǎn)[1]：范圍度窄，量程比一般僅為3∶1～4∶1；測(cè)量重復(fù)性和準(zhǔn)確度不高；現(xiàn)場安裝條件要求高，需要較長的直管段，難以滿足要求；易產(chǎn)生堵塞、信號(hào)失真；壓力損失大。2003年3月國際標(biāo)準(zhǔn)化組織（ISO）正式公布了差壓流量計(jì)的最新國際標(biāo)準(zhǔn)，對(duì)差壓流量計(jì)提出了更為苛刻的要求[2]。20世紀(jì)80年代中期由美國McCrometer公司設(shè)計(jì)的V錐流量計(jì)克服了以上缺點(diǎn)，V錐流量計(jì)采用邊壁節(jié)流[3]，具有一定的流動(dòng)調(diào)整能力，避免流動(dòng)分離，上游無需太長的直管段，具有自清掃能力，可用于臟污介質(zhì)[4]。

近年來，采用CFD對(duì)流量計(jì)的研究越來越多，魏燦等[5]利用FLUENT軟件對(duì)長腰內(nèi)錐式節(jié)流裝置進(jìn)行了優(yōu)化設(shè)計(jì)，陳家慶等[6]用CFD對(duì)標(biāo)準(zhǔn)孔板流量計(jì)的內(nèi)部流場進(jìn)行了數(shù)值模擬，黃珊芳等[7]對(duì)多孔孔板的流出系數(shù)進(jìn)行了數(shù)值模擬研究，Harvill等[8]在ASME會(huì)議上發(fā)表了關(guān)于內(nèi)錐流量計(jì)二維仿真的論文，論述了仿真方案，Manish等[9]對(duì)標(biāo)準(zhǔn)孔板流量計(jì)的內(nèi)部流場進(jìn)行CFD仿真，并探討了標(biāo)準(zhǔn)孔板流量計(jì)的特性，Bobovnik等[10]研究了流動(dòng)條件對(duì)科里奧利流量計(jì)的影響。

本文提出一種新型的節(jié)流式差壓流量計(jì)——環(huán)槽流量計(jì)，該流量計(jì)以流線型的幾何結(jié)構(gòu)為節(jié)流件，采用邊壁逐步收縮的方式對(duì)流體實(shí)現(xiàn)節(jié)流，利用Fluent對(duì)環(huán)槽流量計(jì)進(jìn)行仿真研究，不但可以詳細(xì)了解內(nèi)部流場，而且為其結(jié)構(gòu)優(yōu)化提供了指導(dǎo)方向，以提高流量計(jì)現(xiàn)場計(jì)量的準(zhǔn)確度，指導(dǎo)流量計(jì)的安裝和測(cè)試。

1　環(huán)槽流量計(jì)結(jié)構(gòu)與理論研究

1.1節(jié)流裝置結(jié)構(gòu)

環(huán)槽流量計(jì)的結(jié)構(gòu)如圖1所示，流線型節(jié)流件沿測(cè)量管的軸線安裝，形成一個(gè)由3部分組成的旋轉(zhuǎn)體，前端是半橢圓弧為母線的鏡像旋轉(zhuǎn)體，中部為一等直徑段，且與半橢圓弧相切，后端是以圓弧為母線的鏡像旋轉(zhuǎn)體，圓弧的圓心位于圓柱體的端面沿徑向向外延伸的平面上。節(jié)流件前端橢圓弧的長軸為h1，短軸為d/2，節(jié)流件等直徑段內(nèi)徑為d，等直徑段長度為L，尾部弧的長度為h2。

圖1　節(jié)流件示意圖（單位：m）

1.2測(cè)量原理

被測(cè)流體首先經(jīng)過節(jié)流件前端，流體被逐漸引到測(cè)量管的邊壁，對(duì)流體起到一定的整流作用，緊接著進(jìn)入等直徑段與測(cè)量管內(nèi)壁之間形成的環(huán)形通道，形成標(biāo)準(zhǔn)的槽道流動(dòng)[11]，節(jié)流件較長的尾部可以避免環(huán)形通道中高速流體因劇烈的減速增壓而產(chǎn)生流動(dòng)分離，使永久壓損大大降低。環(huán)槽流量計(jì)的高壓取自節(jié)流件前端，低壓取自環(huán)形槽道的中后部。根據(jù)流體連續(xù)性方程和伯努利方程，得出環(huán)槽流量計(jì)的流量公式為

則

式中：qν——體積流量，m3/h；

C——流出系數(shù)；

ε——可膨脹性系數(shù)；

D——管道內(nèi)徑，m；

d——節(jié)流件等直徑段的內(nèi)徑，m；

ΔP——差壓，Pa；

ρ1——上游流體密度，kg/m3；

ν1——上游流體平均速度，m/s。

2　環(huán)槽流量計(jì)模型及數(shù)值模擬結(jié)果

2.1流線型節(jié)流件建模與求解

因?yàn)樵摴?jié)流件為軸對(duì)稱結(jié)構(gòu)，因此采用圓柱坐標(biāo)，簡化為二維問題來模擬，上下游直管段長度各為10D，使用ICEM CFD軟件構(gòu)造四邊形結(jié)構(gòu)化網(wǎng)格，對(duì)節(jié)流件的位置進(jìn)行局部加密，網(wǎng)格數(shù)量為81171，用Fluent軟件進(jìn)行求解，求解方程組為不可壓縮二維軸對(duì)稱流動(dòng)方程組[12]。

本文采用標(biāo)準(zhǔn)的κ-ε湍流模型，壁面采用標(biāo)準(zhǔn)壁面法，流體為常溫下的水，不可壓縮流體，入口設(shè)為速度入口，湍流參數(shù)選擇湍流強(qiáng)度和水力直徑，出口設(shè)定為流出出口，節(jié)流件和管壁設(shè)為壁面，對(duì)稱軸設(shè)為axis，壓力-速度耦合項(xiàng)采用SIMPLE方案，然后進(jìn)行模擬，以各個(gè)方程的殘差達(dá)到10-3為止，保證計(jì)算結(jié)果充分收斂。

2.2流場的數(shù)值模擬結(jié)果

2.2.1流線分布

圖2為Re=298 564.3，β=0.65時(shí)流線型節(jié)流件與孔板周圍的流線圖，從圖中可以看出流體穩(wěn)定地流過該節(jié)流件，沒有出現(xiàn)死區(qū)和漩渦，沒有發(fā)生流動(dòng)分離，可適用于臟污介質(zhì)，而在孔板的后方出現(xiàn)了較大的漩渦，出現(xiàn)流動(dòng)分離，導(dǎo)致該區(qū)域的壓力產(chǎn)生脈動(dòng)，造成較大的壓力損失，由此可以推斷環(huán)槽流量計(jì)不會(huì)造成太大的壓力損失。

圖2　流線型節(jié)流件與孔板周圍流線對(duì)比圖

2.2.2速度分布

圖3為流線型節(jié)流件前緣50mm、5mm及等直徑段中部的速度分布圖，圖中清楚地顯示出節(jié)流件對(duì)來流速度的調(diào)整過程。在節(jié)流件前緣的50mm處速度沒有受到影響，在前緣5mm處流體速度分布發(fā)生了明顯的改變，趨于對(duì)稱化，當(dāng)流體流經(jīng)節(jié)流件等直徑段的環(huán)形通道時(shí)，速度分布被調(diào)整為穩(wěn)定的槽道流動(dòng)，這也是在此處取低壓的原因，由此可見該節(jié)流件具有一定的流動(dòng)調(diào)整能力，可以縮短上下游的直管段長度。

2.2.3壓力分布

圖3　速度分布圖

圖4為壓力等值線與壓力分布圖，可以看出，在環(huán)形槽道的位置，壓力等值線沿軸向分布很均勻，且沿軸向線性下降，在節(jié)流件的尾部壓力很快恢復(fù)到穩(wěn)定值。

圖4　壓力等值線與壓力分布圖

2.3各參數(shù)變化對(duì)環(huán)槽流量計(jì)特性的影響

2.3.1節(jié)流比β的影響

以管道內(nèi)徑D=0.05 m，L=0.038 m，h1=0.038 m，h2=0.114m為不變參數(shù)，以0.45～0.85為變參數(shù)進(jìn)行數(shù)值模擬，并與相對(duì)應(yīng)的孔板流量計(jì)做比較，如圖5所示。

圖5　不同β時(shí)環(huán)槽流量計(jì)與對(duì)應(yīng)孔板流量計(jì)的流出系數(shù)圖

由圖可知，環(huán)槽流量計(jì)的流出系數(shù)隨著雷諾數(shù)的增加而增大，孔板流量計(jì)與之相反，但在雷諾數(shù)達(dá)到一定值的時(shí)候，兩者的流出系數(shù)都趨于常數(shù)，且環(huán)槽流量計(jì)的測(cè)量下限比孔板流量計(jì)小，說明環(huán)槽流量計(jì)的應(yīng)用范圍更寬。β為0.45，0.55，0.65的環(huán)槽流量計(jì)的流出系數(shù)明顯高于相應(yīng)的孔板流量計(jì)，分別相差62.1%，45.88%，45.35%。對(duì)于環(huán)槽流量計(jì)，流出系數(shù)隨著β的增大，先增大后減小，其中在β=0.45，0.55，0.65時(shí)很接近，在β=0.65時(shí)流出系數(shù)最大，不同β時(shí)，臨界雷諾數(shù)沒有很大的變化。

2.3.2橢圓弧長軸h1、尾部長度h2的影響

四只表演虎間隔2米，成東西方向一字排開，臥于表演場地北側(cè)中央，頭南尾北，所有表演者立于虎后，擊鼓人員位于東北方向。

圖6為不同h1、h2的管壁靜壓圖，可以看出，h1對(duì)差壓和永久壓力損失的影響很小，因而其對(duì)流出系數(shù)的影響也可忽略不計(jì)，圖中3條靜壓曲線在節(jié)流件前緣沒有重合，是因?yàn)楣?jié)流件的不同幾何尺寸造成的。改變h2并不影響低壓差壓信號(hào)，但對(duì)永久壓損還是有所影響，隨著h2的增大，永久壓損增大。綜合考慮各種因素，h2小一點(diǎn)比較好。

圖6　不同h1、h2的管壁靜壓圖

圖7　不同L時(shí)的管壁靜壓圖

圖8　流量與?的關(guān)系圖

2.3.3等直徑段長度L的影響

圖7為入口速度為4 m/s，等直徑段長度L=38，57，76mm時(shí)的管壁靜壓圖。永久壓損為節(jié)流件前緣與節(jié)流件后6D處的差壓值，由圖可以看出當(dāng)流體流經(jīng)節(jié)流件時(shí)靜壓急速減小，達(dá)到相同的壓力P1，此階段主要是由于節(jié)流件前端的節(jié)流作用，流體流經(jīng)等直徑段時(shí)壓力繼續(xù)減小，達(dá)到最小值P2，L愈大達(dá)到的壓力最小值最小，這是因?yàn)椴鄣懒鲃?dòng)時(shí)靜壓的降低是由于摩擦壓損造成的，摩擦壓損與長度成正比，所以L越大，造成的摩擦壓損越大，隨后壓力值不斷增大，恢復(fù)穩(wěn)定值，L愈大，永久壓損越大，因此應(yīng)該盡量減小L。

2.4量程比

2.5壓力損失

δp為流經(jīng)流量計(jì)的永久壓損，用δp/Δp來表征壓力損失占整個(gè)差壓的百分比，如圖9所示，可以看出，在相同β下，環(huán)槽流量計(jì)的壓力損失百分比小于孔板流量計(jì)，且隨著雷諾數(shù)的增大，壓損比越來越小，可達(dá)到24%，孔板流量計(jì)的壓損比基本維持在56%，但在低雷諾數(shù)時(shí)，環(huán)槽流量計(jì)的壓損比還是較大。

2.6流出系數(shù)公式擬合

流出系數(shù)C=f（ReD，節(jié)流件類型，D，β）[13]，利用牛頓法對(duì)模擬的數(shù)據(jù)進(jìn)行擬合得出流出系數(shù)的公式為

圖9　兩種流量計(jì)的壓損對(duì)比

3　結(jié)束語

穩(wěn)定的壓力信號(hào)可以保證流量計(jì)的準(zhǔn)確度和重復(fù)性，從流場、速度和壓力分布可以看出環(huán)槽流量計(jì)可以將來流快速地調(diào)整為環(huán)形槽道流動(dòng)，從而保證了靜壓的穩(wěn)定性。

由模擬的結(jié)果我們可以得出以下結(jié)論：

1）環(huán)槽流量計(jì)不會(huì)產(chǎn)生死區(qū)和漩渦，不會(huì)發(fā)生流動(dòng)分離，可適用于臟污介質(zhì)，不易磨損；具有流動(dòng)調(diào)整能力，與傳統(tǒng)差壓流量計(jì)比，不需太長的上下游直管段；量程比可達(dá)到7∶1，與孔板流量計(jì)相比應(yīng)用范圍更大，壓力損失更小。

2）隨著雷諾數(shù)的增加，流出系數(shù)逐漸增大，并趨于穩(wěn)定值，在β=0.65處可達(dá)到較大的流出系數(shù),約為0.89。節(jié)流件前端的長度對(duì)流出系數(shù)和壓力損失沒影響，h2對(duì)永久壓損影響較大，L愈小，永久壓損越小，綜合考慮h1=d、h2=d、L=d、β=0.65的節(jié)流件為最優(yōu)結(jié)構(gòu)。

3）根據(jù)模擬結(jié)果得出的流出系數(shù)公式可以給工程實(shí)際應(yīng)用提供參考。

CFD仿真模擬作為一種測(cè)試和輔助設(shè)計(jì)方法，為環(huán)槽流量計(jì)的研究提供了可靠的理論依據(jù)和基礎(chǔ)數(shù)據(jù)，用以指導(dǎo)流量計(jì)的安裝和使用。

參考文獻(xiàn)

[1]吳小平.差壓流量計(jì)的發(fā)展現(xiàn)狀[J].常州工學(xué)院學(xué)報(bào)，2007，20（4）：49-51.

[2] Measurement of fluid flow by means of pressure differential devices inserted in circular cross-section conduits running full-Part 2：Orifice plates：ISO5167-2：2003[S]. [2003-03-01]

[3] STEVEN R. Horizontally installed cone differential pressure meter wet gas flow performance[J]. Flow Measurement and Instrumentation，2009（5）：152-167.

[4]劉躍.內(nèi)錐流量計(jì)的探討[J].自動(dòng)化儀表，2004，25（3）：38-42.

[5]魏燦.長腰內(nèi)錐節(jié)流裝置的數(shù)值仿真與實(shí)驗(yàn)研究[D].天津：天津大學(xué)，2011.

[6]陳家慶，王波.標(biāo)準(zhǔn)孔板流量計(jì)內(nèi)部流場的CFD數(shù)值模擬[J].實(shí)驗(yàn)流體力學(xué)，2008，22（2）：51-55.

[7] HUANGS F，MA T Y，WANG D，et al. Study on discharge coefficient of perforated orifices as a new kind of flowmeter[J]. Experimental Thermal and Fluid Science，2013（46）：74-83.

[8] HARVILL L R，HARVILL T L. Numerical study of flow through a v-cone[J]. American Society of Mechanical Engineers-Fluids Engineering Division（Publication）FED，1988（70）：23-26.

[9] MANISH S，JYESHTHARAJ S，JOSHI B，et al. Analysis of flow through an orificemeter：CFD simulation [J]. Chemical Engineering Science，2013（71）：300-309.

[10] BOBOVNIK G，KUTIN J，BAJSIC I. The effect of flow conditions on the sensitivity of the coriolis flow meter[J]. Flow Measurement and Instrumentation，2004（15）：69-76.

[11]鐘偉，明曉.紡錘體流量計(jì)的流場數(shù)值模擬[J].計(jì)量學(xué)報(bào)，2007，28（3）：262-265.

[12]程勇，汪軍，蔡小舒.低雷諾數(shù)的孔板計(jì)量數(shù)值模擬及其應(yīng)用[J].計(jì)量學(xué)報(bào)，2005，26（1）：57-59.

[13]孫淮清，王建中.流量測(cè)量節(jié)流裝置設(shè)計(jì)手冊(cè)[M].北京：化學(xué)工業(yè)出版社，2005：12-17.

（編輯：李妮）

Numerical simulation and optimization of ring groove flow meter

NIU Zhijuan，ZHAO Qi，YANG Xuefeng
（College of Chemical Engineering，Sichuan University，Chengdu 610065，China）

Abstract:Studies on the internal flow field and structural optimization of flow meters is important for improving the measurement precision of flow meters. The simulation test of computational fluid dynamics（CFD）was carried out on a new type of differential pressure flow meter or ring groove flow meter. The study is to explore how some key factors have influenced the discharge coefficient and pressure loss of the meter. These key factors include equivalent diameter ratio，fore-and-aft length，equal diameter length and Reynolds number. The experimental results show that the discharge coefficient increases gradually and then reaches a steady value while the Reynolds number keeps increasing. When the equivalent diameter ratio is increasing，the discharge coefficient increases as well at first and then decreases. The fore-and-aft length has little effect on the discharge coefficient. The permanent pressure loss declines as the tail becomes longer and the equal diameter section is shortened. An equation for discharge coefficient was fitted out according to the results above. As a complementary design and certification tool for flow meter，the CFD numerical simulation can help guide the field testing of ring groove flow meters.

Keywords:ring groove flow meter；pressure loss；discharge coefficient；numerical simulation；measurement

作者簡介：牛志娟（1990-），女，河南安陽市人，碩士研究生，專業(yè)方向?yàn)閭髻|(zhì)與分離。

收稿日期：2015-05-25；收到修改稿日期：2015-07-30

doi：10.11857/j.issn.1674-5124.2016.02.028

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1674-5124（2016）02-0122-05