1 引言
目前，在城市燃氣計量中應(yīng)用著基于各種測量原理的多種儀表。但是，經(jīng)過分析和比較后發(fā)現(xiàn)：這些燃氣流量測量儀表均存在著流量下限瓶頸問題，即在儀表流量下限以下流量不能實現(xiàn)準確計量或根本不能計量問題，城市燃氣小流量漏計問題突出。為了解決該問題，我們提出并設(shè)計研制了一種新型的城市燃氣流量計，能夠?qū)崿F(xiàn)對城市燃氣的寬量程計量。該流量計在燃氣大流量段時應(yīng)用流體振蕩測量原理進行測量，即采用渦街流量計的測量原理完成大流量測量[1]。
由于漩渦發(fā)生體的阻流作用，燃氣管道中的流場變得異常復(fù)雜，難以解析地求得流場分布情況，所以至今人們對漩渦發(fā)生體后漩渦的形成和脫落的認識大多依賴于經(jīng)驗和實驗。這對于優(yōu)化設(shè)計新型寬量程燃氣流量計都十分不利[2][3]。
本文從描述漩渦運動的經(jīng)典流體力學(xué)基本方程出發(fā)，以計算流體力學(xué)軟件FLUENT為平臺對寬量程燃氣流量計的方柱繞流流場進行數(shù)值模擬，并將模擬結(jié)果與實測數(shù)據(jù)比較，結(jié)果表明數(shù)值模擬結(jié)果與管道燃氣實際流動狀況具有很好的一致性。
2 數(shù)值模擬模型
2.1流體漩渦脫落現(xiàn)象的數(shù)值模擬方法
單鈍體繞流的漩渦脫落特性和受力情況，除采用實驗方法進行測定外，還可用數(shù)值模擬方法進行計算分析。
數(shù)值模擬方法主要有兩類，即經(jīng)驗的模型法和直接流場模擬法[4][5]。
經(jīng)驗的模型法，這類方法不考慮具體流體結(jié)構(gòu)，將流體及其中的振蕩物體視作一個整體系統(tǒng)。然后用一組適合的模型方程對其進行描述，其中包括采用一些已知結(jié)果的或用經(jīng)驗確定的系數(shù)，以便求解后可以較好地再現(xiàn)系統(tǒng)的運動特性，并從總體上可對現(xiàn)象本身的物理本質(zhì)上有一個較直觀的了解。目前得到廣泛使用的主要是尾流振蕩器模型和相關(guān)模型，可用于求解各種具體問題。但這類方法主要依賴于實驗結(jié)果，因而有一定局限性。
直接流場模擬法有直接求解Navier-Stokes方程式的各種差分方法，有限分析解方法，有限元方法，有限體積法、譜方法和基于邊界層方程的各種正、反解法以及粘性/無粘性干擾方法，離散渦方法或直接求解Navier-stokes方程等。其計算過程包括流動模型選擇(列出控制方程式)，方程式的離散、離散方程式的求解等。這些方法可以給出整個流場情況的詳細描述，但計算相當復(fù)雜，特別是流動Re數(shù)較高時，更為困難。在流場計算中應(yīng)用較多的為渦方法和有限差分和有限體積法。
鑒于本文的流體力學(xué)數(shù)值模擬使用的是FLUENT商業(yè)軟件，故將著重介紹非定常流體力學(xué)問題的有限體積法。
有限體積法(finite volume method)，簡稱FVM，就是在物理空間中選定的控制體積上把積分型守恒律直接離散的一類數(shù)值方法。離散一方面是指把計算區(qū)域分成網(wǎng)絡(luò)(或單元)，另一方面是指把積分守恒律離散成線性或非線性代數(shù)方程組。
有限體積法(FVM)將計算域劃分成若干規(guī)則或不規(guī)則形狀的單元或控制體。有限體積法是計算出通過每個控制體邊界沿法向輸入(出)的流量和動量通量后，對每個控制體分別進行流量和動量平衡計算，便得到計算時段末各控制體平均壓力和流速。
由此，F(xiàn)VM正是對于推導(dǎo)原始微分方程所用控制體途徑的回歸，其物理意義更直接明晰。如跨邊界通量的計算只使用時段初值，為顯式FVM；反之，當涉及時段始末的值時，則為隱式FVM。因為跨控制體間界面輸運的通量，對相鄰控制體來說大小相等，方向相反，故對整個計算域而言，沿所有內(nèi)部邊界的通量相互抵消。對由一個或多個控制體組成的任意區(qū)域，以至整個計算域，都嚴格滿足物理守恒律，不存在守恒誤差，并且能正確計算間斷。
2.2流體力學(xué)基本控制方程組
流體在寬量程流量計燃氣管道中的流動為時變渦流。由于在數(shù)值模擬的流量范圍內(nèi)流體都是在湍流狀態(tài)下，故需要對湍流進行處理。
現(xiàn)在工程計算中的湍流模型很多，選取不同的湍流模型會直接影響計算的結(jié)果。經(jīng)過多次試算，本文數(shù)值模擬中選擇FLUENT中湍流運動的RNG k- 方法對湍流進行處理。
RNG的基本思想是通過在任意空間尺度上的一系列連續(xù)的變換，對原本十分復(fù)雜的系統(tǒng)或過程實現(xiàn)粗分辨率的描述。采用RNG k- 湍流模型來模擬方柱繞流二維流場，將RNG方法用于雷諾平均N-S方程并引入湍能k及其耗散率 便可得到如下控制方程[6][7]：
連續(xù)方程：
方柱等鈍體繞流的數(shù)學(xué)模型都是建立在質(zhì)量守恒定律上的連續(xù)性方程、動量守恒定律上的運動方程和熱力學(xué)第一定律上的本構(gòu)方程基礎(chǔ)上的，決定方程解的因素是邊界條件的變化[8]。
2.3物理模型、初始條件及邊界條件
本文中將對管道燃氣計量中所用的方柱單鈍體燃氣繞流進行數(shù)值模擬研究，其物理模型如圖2-1所示，求解區(qū)域放大如圖2-2所示。
本例中為對比應(yīng)用DN200寬量程燃氣流量計，在模擬中設(shè)定方柱體迎風(fēng)面的寬度為0.05m，求解區(qū)域的寬度為0.20m，長度為0.50m。將漩渦發(fā)生體放置在離右側(cè)入口0.10m處的正中央。
為了對離散后的流體力學(xué)方程組進行數(shù)值計算，需要給出定解條件，包括邊界條件和初始條件。
①入口邊界
入口邊界設(shè)定速度入口，給定流體流速，為與實測值對比，速度分別取為：5.3 m/s，10.5 m/s， 20.2 m/s， 29.8 m/s，39.1 m/s, 50.6 m/s ，65.7m/s等7個流速情況。
②出口邊界
出口邊界條件設(shè)定為壓力出口，壓力出口的壓力為一個大氣壓，即表壓為零。
③壁面條件
包括流體流動管壁和方柱體壁面的條件，對這兩種壁面采用相同的處理方法。
④網(wǎng)格劃分
由于該問題的幾何邊界比較簡單，使用GAMBIT把求解區(qū)域按結(jié)構(gòu)化矩形網(wǎng)格進行劃分。
3 數(shù)值模擬結(jié)果
數(shù)值模擬中采用方柱鈍體形式，網(wǎng)格劃分如圖2-1和圖2-2所示，網(wǎng)格采用結(jié)構(gòu)化網(wǎng)格，入口邊界設(shè)定速度入口，給定流體平均流速。
圖3入口平均流速39.1 m/s時方柱體典型的一個漩渦脫落周期的速度等值線圖（天然氣）
圖3表示了基于RNG方法的湍流模型能夠捕獲方柱體燃氣繞流流場的不穩(wěn)定和劇烈分離的特征，描述了天然氣介質(zhì)中一個完整的漩渦脫落過程。
限于文章篇幅，這里只給出入口平均流速39.1 m/s時管道方柱天然氣繞流典型的一個漩渦脫落周期的速度場等值線圖。而動壓等值線圖、流線圖及其他流速和燃氣介質(zhì)的相關(guān)圖與此類似。
4 數(shù)值模擬與實測結(jié)果比較
為與實測結(jié)果進行對比，數(shù)值模擬時流體介質(zhì)同樣分別采用天然氣、人工煤氣等不同燃氣介質(zhì)的物性參數(shù)。溫度設(shè)定為實測溫度（294.6K），壓力為一個標準大氣壓（101325Pa）。
實測裝置采用由哈爾濱工業(yè)大學(xué)投資建立的燃氣流量寬量程計量標準裝置。該裝置于2004年8月獲得黑龍江省質(zhì)量技術(shù)監(jiān)督局授權(quán)，經(jīng)過2年多的檢定運行，并與國家原油大流量站（大慶）檢定裝置對比表明，該裝置具有良好的穩(wěn)定性、重復(fù)性，裝置精度0.5級。
數(shù)值模擬與實測結(jié)果比較見圖4所示。
圖4 數(shù)值模擬頻率與實測頻率比較
為與實測值對比，流體平均流速分別取為：
5.3 m/s，10.5 m/s， 20.2 m/s， 29.8 m/s，39.1 m/s, 50.6 m/s ，65.7m/s等7個流速情況。
5、結(jié)論
本文從描述漩渦運動的經(jīng)典流體力學(xué)基本方程出發(fā)，以計算流體力學(xué)軟件FLUENT為平臺對管道單鈍體燃氣繞流計量流場進行了數(shù)值模擬，并將模擬結(jié)果與實測數(shù)據(jù)進行了比較。結(jié)果表明數(shù)值模擬結(jié)果與管道燃氣實際流動狀況具有很好的一致性，計量流場中漩渦的脫離頻率不依賴于介質(zhì)的物性，在漩渦發(fā)生體尺寸一定的情況下，數(shù)值模擬和實測頻率均正比于被測介質(zhì)的流速。在新型燃氣流量計的研制中可以應(yīng)用數(shù)值模擬方法來指導(dǎo)和優(yōu)化流量計的結(jié)構(gòu)設(shè)計。
參考文獻：
[1] Li Chaohui,Dai Jingmin. A gas mass flowmeter with two measurement limits [J]. Measurement Techniques,2005,48(5): 487-491
[2] 夏泰淳.工程流體力學(xué)[M].上海:上海交通大學(xué)出版社,2006.
[3] 童秉剛，張炳暄,崔爾杰主編.非定常流與渦運動.北京:國防工業(yè)出版社，1993
[4] 孫志強,張宏建,黃詠梅. 渦街流量計流場特性的數(shù)值仿真研究[J].自動化儀表.2004,(25),5:10-13
[5] 王遠成,吳文權(quán).基于RNG k- 湍流模型鈍體繞流的數(shù)值模擬[J].上海理工大學(xué)學(xué)報.2004,(26),6:519-523
[6]YasutakaNagano,Yoshihiro Itazu.Renormalization group theory for turbulence:Assessment of the Yakhot Oiszag Smith theory [J]. Fluid Dynamics Research.1997,20(6):157-172
[7] Speziale G，Thangam S .Analysis of an RNG fased turbulence model for separated flows [J]. Int J Engrfg Sci.1992,30(2):1379-1384.
[8] 彭杰綱,傅新,陳鷹.雙鈍體渦街流量計流體振動特性實驗研究[J].機械工程學(xué)報. 2002，(26),8:519-523
