用CFX軟件對超聲波燃?xì)獗須怏w流道仿真研究

翟義然，趙　勇，胡小川，劉　勛，劉　義，張　彬

（1.四川大學(xué)電子信息學(xué)院，四川 成都 610065；2.成都千嘉科技有限公司，四川 成都 610211）

用CFX軟件對超聲波燃?xì)獗須怏w流道仿真研究

翟義然1，2，趙勇2，胡小川1，劉勛2，劉義2，張彬1

（1.四川大學(xué)電子信息學(xué)院，四川 成都 610065；2.成都千嘉科技有限公司，四川 成都 610211）

針對超聲波燃?xì)獗淼脑瓍⒖荚O(shè)計流道和新改進(jìn)設(shè)計流道的氣體流動特性，利用流體動力學(xué)仿真軟件ANSYS中的CFX軟件進(jìn)行仿真計算和分析。獲得在9個不同的流量點上兩種超聲波燃?xì)獗淼牧鞯纼?nèi)氣體流動仿真數(shù)據(jù)，包括兩種流道的整體流道內(nèi)及超聲波測量流道部分斜截面的氣體流動速度分布。對仿真結(jié)果分析顯示，兩種流道在性能上基本一致，但新改進(jìn)流道在大流量和中流量上氣體流動速度分布性能比原參考設(shè)計流道性能更好；另外，與整體流道內(nèi)氣體流動速度分布相比，氣體流動速度在超聲波測量斜截面分布相對更均勻和穩(wěn)定。對采用新改進(jìn)流道設(shè)計制造的超聲波燃?xì)獗順訖C進(jìn)行測試，準(zhǔn)確度達(dá)到設(shè)計要求，且樣機的測量誤差穩(wěn)定。

超聲波；氣體流速；CFX仿真；燃?xì)獗?/p>

0　引　言

超聲波燃?xì)獗淼墓ぷ髟硎抢贸暡▊鞑ニ俣仁苋細(xì)饬鲃铀俣鹊挠绊憗碛嬃咳細(xì)獾捏w流量［1-2］?；诔暡y量的超聲波燃?xì)獗砭哂畜w積小、質(zhì)量輕、壓損小、無機械部分不受機械磨損影響、計量準(zhǔn)確度高，以及可實現(xiàn)階梯計價和方便數(shù)據(jù)遠(yuǎn)傳等優(yōu)點。因此，超聲波燃?xì)獗硖娲F(xiàn)在廣泛使用的膜式燃?xì)獗韺⒊蔀橐环N必然。目前，歐美和日本已經(jīng)研制出用于居民燃?xì)庥嬃康募矣贸暡ㄈ細(xì)獗?，并開始推廣應(yīng)用。我國的一些科研院所和相關(guān)企業(yè)也在開展超聲波燃?xì)獗淼难邪l(fā)工作［2］。

超聲波燃?xì)獗淼谋眢w結(jié)構(gòu)設(shè)計決定著表內(nèi)燃?xì)饬鲃犹匦?，而燃?xì)饬鲃拥姆€(wěn)定性直接影響到超聲波燃?xì)庥嬃繙?zhǔn)確度，實際工程設(shè)計應(yīng)用中對設(shè)計者經(jīng)驗和產(chǎn)品實驗有很大的依賴性，其設(shè)計周期長、成本高。近年來，計算流體動力學(xué)CFD軟件在實際工業(yè)設(shè)計中得到廣泛應(yīng)用，利用CFD各種流體模型在設(shè)計階段進(jìn)行仿真，研究結(jié)構(gòu)性能，可在結(jié)構(gòu)加工前使研究對象的速度場清晰直觀地顯示出來，因此，CFD仿真是一種理想的表體結(jié)構(gòu)設(shè)計分析手段，可以縮短設(shè)計驗證時間，提高工作效率［3-6］。ANSYS的CFX軟件是一款重要的CFD軟件，在過去20年被廣泛應(yīng)用于分析流體流動問題［7-8］。利用CFX軟件可進(jìn)行燃燒、多相流、化學(xué)反應(yīng)等復(fù)雜流場的仿真分析，廣泛用于航空航天、旋轉(zhuǎn)機械、能源、汽車、電子、生物技術(shù)、水處理、防火安全、冶金、環(huán)保等領(lǐng)域。

家用超聲波燃?xì)獗淼脑瓍⒖荚O(shè)計中入口處流道過窄，不方便安裝開關(guān)閥門，而未來的家用超聲波燃?xì)獗頃驗檫\行安全控制、電源和預(yù)付費等原因，必須安裝開關(guān)閥門。因此，在原參考設(shè)計基礎(chǔ)上，改進(jìn)設(shè)計了一款新超聲波燃?xì)獗砹鞯溃饕窃黾恿舜怪绷鞯啦糠值目臻g［9］。因此，有必要對原參考流道和設(shè)計的新改進(jìn)流道進(jìn)行氣體流動性能仿真研究，以確保新改進(jìn)流道性能達(dá)到原參考流道性能要求。應(yīng)用ANSYS 15.0中的CFX軟件對燃?xì)獗砹鞯纼?nèi)部氣體流動特性進(jìn)行仿真分析［10］，模擬兩種流道在不同流量點下的氣體流動，并對兩種流道的整體流道內(nèi)氣體流動和超聲波測量斜截面上氣體流動進(jìn)行對比研究。

1　模型建立

1.1氣體流道模型建立

CFX軟件的內(nèi)部程序是利用計算機對流體的守恒控制偏微分方程組進(jìn)行求解，進(jìn)而得到流場連續(xù)域上的離散分布情況，并對其流動特性進(jìn)行模擬［7］。通過CFX仿真分析可以模擬氣體流場的變化，并在較短的時間內(nèi)準(zhǔn)確預(yù)測流場的性能，進(jìn)而通過改變流道結(jié)構(gòu)設(shè)計參數(shù)，獲得更佳的流動效果和測量性能。

通常采用Navier_Stokes（N-S）方程建立數(shù)學(xué)模型來模擬宏觀條件下的氣體流動特性［11］。N-S方程包含有3個偏微分方程：連續(xù)性方程、動量方程和能量方程［4］。這些方程在離散化之后可以求出數(shù)值近似解。模擬仿真以N-S方程作為數(shù)值計算的控制方程組，采用k-ε湍流模型，對超聲波燃?xì)獗韮?nèi)氣體的流動過程進(jìn)行仿真計算。

超聲波燃?xì)獗淼膬?nèi)部流道結(jié)構(gòu)采用CAD軟件生成，腔體內(nèi)的下方部位設(shè)計有一段相對狹窄的長方體多層流道結(jié)構(gòu)，可以使氣體流動更加穩(wěn)定，超聲波就在該段流道測量氣體流動速度。

1.2網(wǎng)格劃分

在Workbench中導(dǎo)入由CAD生成的超聲波燃?xì)獗須怏w流道結(jié)構(gòu)圖，采用自動網(wǎng)格劃分方法對整體流道進(jìn)行網(wǎng)格化。網(wǎng)格化后的原參考流道內(nèi)部結(jié)構(gòu)如圖1（a）所示，共得到136059個節(jié)點和724740個四面體。新改進(jìn)的流道內(nèi)部結(jié)構(gòu)如圖1（b）所示，共得到133 280個節(jié)點和710817個四面體，z軸方向向內(nèi)。從劃分的網(wǎng)格結(jié)構(gòu)圖中可以看出，在內(nèi)部結(jié)構(gòu)的轉(zhuǎn)角處和多層流道結(jié)構(gòu)部分網(wǎng)格劃分相對更細(xì)密，而多層流道結(jié)構(gòu)部分是本研究的重點，通過仿真計算可以得到超聲測量部位的氣體流動速度分布。

圖1　整體氣體流道網(wǎng)格劃分和斜截面示意圖

2　仿真分析

2.1模型建立和邊界條件設(shè)定

假設(shè)超聲波燃?xì)獗韮?nèi)流動的為理想氣體。設(shè)置流道的入口和出口面，如圖1所示，其他全為固體壁面。仿真模型中的氣體從入口流入，從出口流出。針對本設(shè)計中超聲波燃?xì)獗頊y量的有效測量范圍，選擇9個典型流量點進(jìn)行氣體流動速度仿真，包括3個大流量、3個中流量和3個小流量，分別為6，4，2.5m3/h；1.2，0.8，0.5m3/h；0.04，0.025，0.016m3/h。仿真時根據(jù)某一個流量和出口或入口面積計算并設(shè)置相應(yīng)的氣體流動速度。

2.2求解控制

設(shè)置收斂條件為RMS=0.0001，最大收斂迭代步數(shù)為300。

3　結(jié)果分析

3.1整體流道氣體流動分析

首先選取流量點6m3/h，對原參考流道和新改進(jìn)流道分別獲得整個流道在z=0.0的xy軸向截面。原參考流道截面由11010個三角片組成，三角片上氣體流動速度的最大值（Max）為14.49m/s，最小值（Min）為0.04 m/s；新改進(jìn)流道的該截面由10 811個三角片組成，三角片上氣體流動速度的Max為14.35m/s，Min為0.06m/s，截面內(nèi)氣體流動速度分布分別如圖2所示?？梢娬w流道內(nèi)氣體流動速度分布變化范圍比較大，而參考流道和新改進(jìn)流道內(nèi)氣體整體流動速度分布情況類似。在多層流道部分氣體流動速度相對更快，但氣體流動速度分布較為均勻穩(wěn)定。根據(jù)超聲波測量理論可知，均勻穩(wěn)定的氣流有利于用超聲波準(zhǔn)確測量氣體流動速度。

圖2　原參考流道和新改進(jìn)流道整體氣體流道中xy軸向截面上氣體流動速度分布

進(jìn)一步選取由大到小的9個流量點對整體流道空間進(jìn)行仿真。分別獲取原參考流道和新改進(jìn)流道在不同氣體流動速度下的整體流道氣體流動速度的平均值（Mean）、標(biāo)準(zhǔn)方差（Std）、Max和Min，結(jié)果如表1所示。

表1　原參考流道和新改進(jìn)流道的整體流道在9個流量點的氣體流動速度仿真結(jié)果

將原參考流道和新改進(jìn)流道仿真結(jié)果中整體流道內(nèi)氣體流動速度的Mean繪制成曲線，并在具體流量點上繪制出對應(yīng)的Std。得到兩種流道的整體流道內(nèi)氣體流動速度Mean和Std曲線如圖3所示。

3.2超聲波測量斜截面的氣體流動速度分析

超聲波測量的多層流道內(nèi)氣體流動特性是本研究的重點，在圖1（b）中箭頭A的方向獲得多層流道部分的一個xz軸向截面，該截面如圖1（c）所示。兩個超聲波測量傳感器安裝在圖中所示的位置，因此，進(jìn)一步在多層流道上得到一個斜截面如圖1（d）所示，原參考流道和新改進(jìn)流道的斜截面分別被網(wǎng)格自動劃分為180個和183個三角片單元。

在9個流量點獲得的仿真結(jié)果如表2所示。將原參考流道和新改進(jìn)流道氣體流動速度Mean分別繪制曲線，并在每個流量點上繪制出對應(yīng)的Std，得到兩種流道超聲波測量斜截面氣體流動速度分布Mean和Std曲線如圖3所示。

圖3　原參考流道和新改進(jìn)流道在整體流道內(nèi)和超聲波測量斜截面上氣體流動速度Mean和Std曲線

表2　原參考流道和新改進(jìn)流道在超聲測量斜截面的9流量點的氣體流動速度仿真結(jié)果

在9個流量點分別仿真得到兩種流道的超聲波測量斜截面氣體流動速度分布，圖4（a）和圖4（b）分別為原參考流道和新改進(jìn)流道在流量點6，1.2，0.04m3/h時的各斜截面氣體流動速度分布圖?？梢钥闯觯暡y量斜截面內(nèi)氣體流動速度變化范圍不大，并且，新改進(jìn)流道斜截面氣體流動速度分布范圍變窄。

圖4　原參考流道和新改進(jìn)流道在3個不同流量點的斜截面上氣體流動速度分布

從表1、表2中可以看出，兩種流道在仿真氣體流動速度的Mean、Std、Max和Min等參數(shù)上基本一致，能保證設(shè)計的新改進(jìn)流道與原來超聲波測量電子器件部分相匹配。進(jìn)一步對比表1中原參考流道和新改進(jìn)流道的整體流道仿真結(jié)果可以發(fā)現(xiàn)，在大流量和中流量時，即6～0.5m3/h 6個流量點，新改進(jìn)流道的整體流道內(nèi)氣體流動Mean在大于0.5 m3/h時全部下降，Std和Max全部下降，Min有升有降，但變化相對很小。而對于小流量點0.04～0.016m3/h，Mean、Std、Max和Min基本保持不變，或者變化很小。再對比表2中原參考流道和新改進(jìn)流道的超聲測量斜截面仿真結(jié)果同樣可以發(fā)現(xiàn)，在大流量和中流量情況下，即6～0.5m3/h 6個流量點，新改進(jìn)流道的超聲測量斜截面上Mean略有升降；Std除了在6.0 m3/h略有上升，其余部分均下降；Max除了在0.8 m3/h時略有上升外，其余部分均下降；而Min則全部有所上升。而對小流量點0.04～0.016 m3/h，新改進(jìn)流道與原參考流道相比，仿真結(jié)果Mean、Std、Max和Min基本保持不變或變化很小。上述對比結(jié)果說明，與原參考流道相比，在大流量和中流量點時，新改進(jìn)流道在整體流道和超聲測量部分流道氣體流動速度分布范圍基本上變得更小，說明氣體流動更為穩(wěn)定，這有助于超聲波在大流量和中流量時更準(zhǔn)確測量氣體流動速度；而小流量部分兩種流道仿真數(shù)據(jù)基本一致。

在各個流量點，超聲波測量的斜截面上的氣體流動速度與整體流道內(nèi)的氣體流動速度相比，平均氣體流動速度更高，而相對應(yīng)的氣體流動速度Max和Min更接近平均氣體流動速度，Std更小，即氣體流動速度變化范圍更小。該結(jié)論同樣也可以從圖3中整體流道內(nèi)和超聲波測量斜截面上氣體流動速度分布Mean和Std曲線得到。因此，超聲波測量的多層流道內(nèi)氣體流動相對整體流道更加平穩(wěn)，從而為超聲波的準(zhǔn)確測量提供了更好的條件。

3.3采用新改進(jìn)流道設(shè)計樣機實際測試結(jié)果

采用新改進(jìn)流道設(shè)計結(jié)構(gòu)做出一款超聲波燃?xì)獗順訖C，在環(huán)境溫度為20.5℃，濕度76%RH條件下進(jìn)行兩次測試，在9個流量點的測試結(jié)果如表3所示。實際樣機測試結(jié)果滿足大流量和中流量點誤差＜1.5%和小流量點誤差＜3%的1.5級超聲波燃?xì)獗淼脑O(shè)計要求。而且對比實驗結(jié)果可以看出超聲波燃?xì)獗順訖C的兩次測量誤差穩(wěn)定。

表3　用新改進(jìn)流道設(shè)計的超聲波燃?xì)獗順訖C測試結(jié)果

4　結(jié)束語

對比原來參考設(shè)計的和新改進(jìn)設(shè)計的兩種超聲波燃?xì)獗韮?nèi)部結(jié)構(gòu)，用ANSYS的CFX仿真軟件進(jìn)行了氣體流體力學(xué)仿真研究，對原參考流道和新改進(jìn)流道分別獲得了整體流道的和超聲波測量斜截面的氣體流動速度分布數(shù)據(jù)。數(shù)據(jù)分析結(jié)果顯示，新改進(jìn)流道與原參考流道的氣體流動性能基本一致，并且新改進(jìn)流道在大流量和中流量時能獲得更加平穩(wěn)的氣體流動速度。對采用新改進(jìn)流道設(shè)計制造的超聲波燃?xì)獗順訖C初步測試，結(jié)果顯示在9個流量點都能滿足設(shè)計準(zhǔn)確度要求，而且兩次測量的誤差值相差不大，測量結(jié)果穩(wěn)定。相關(guān)研究成果為改進(jìn)超聲波燃?xì)獗淼膬?nèi)部氣體流道結(jié)構(gòu)設(shè)計提供一定的理論依據(jù)，也可為提高超聲波氣體測量的穩(wěn)定度和準(zhǔn)確性提供一定的保障。

［1］LYNNWORTH L C，LIU Y.Ultrasonic flowmeters：halfcentury progress report，1955—2005［J］.Ultrasonics，2006（44）：1371-1378.

［2］于洋，宗光華，丁鳳林.超聲波流量測量中流速計算方法的對比［J］.北京航空航天大學(xué)學(xué)報，2013（1）：37-41.

［3］BERNARDOS，MORIM，PERESA P，et al.3-D computational fluid dynamics for gas and gas-particle flows in a cyclone with different inlet section angles［J］. Powder Technology，2006，162（3）：190-200.

［4］高井輝.基于CFX的離心式壓氣機內(nèi)部流場數(shù)值研究［D］.沈陽：大連理工大學(xué)，2011.

［5］姚美紅，張恒峰.基于CFX的氣體靜壓軸承數(shù)值研究［J］.中國機械工程，2012，23（22）：2681-2683.

［6］JIA H Y，KONG J Y，ZHOU S Z，et al.Study on the method for CFX numerical simulation of flow field in screw conveyer of sand fracturing blender［J］.Applied Mechanics and Materials，2013（318）：207-211.

［7］謝龍漢，趙新宇，張炯明.ANSYS CFX流體分析及仿真［M］.北京：電子工業(yè)出版社，2012：104-112.

［8］李麗丹，李聲.基于CFX和Workbench的數(shù)值仿真技術(shù)［J］.中國測試，2010，36（5）：79-83.

［9］梅景康裕，乾善紀(jì)，長岡行夫，等.超聲波流體測量裝置：200480004936.6［P］.2008-07-16.

［10］朱紅鈞.ANSYS 15.0幾何建模與網(wǎng)格劃分實戰(zhàn)指南［M］.北京：人民郵電出版社，2014：79-83.

［11］范晉偉，李海涌，陳東菊，等.基于速度滑移的氣膜流態(tài)建模與CFX仿真分析［J］.潤滑與密封，2012，37（12）：22-26.

（編輯：李妮）

Simulation study of gas flowing channel of ultrasonic gas meter based on CFX software

ZHAI Yiran1，2，ZHAO Yong2，HU Xiaochuan1，LIU Xun2，LIU Yi2，ZHANG Bin1
（1.College of Electronics and Information Engineering，Sichuan University，Chengdu 610065，China；2.Chengdu Qianjia Technology Co.，Ltd.，Chengdu 610211，China）

The gas fluid dynamics of an originally suggested flowing channel and an improved flowing channel were simulated by the CFX software of ANSYS.Gas flowing velocities at 9 typical flow points were calculated and analyzed，including the gas velocity distribution of both the whole flowing channel and the oblique section of ultrasonic measurement.The simulation results indicate that the two channels are almost the same in performance.But the improved channel is better in gas velocity distribution at maximal and middle flowing points.Moreover，the simulation results show that，compared with the gas velocity distribution inside the whole channel，the oblique section where ultrasonic wave passes through is more uniform and stable.The test accuracy of an ultrasonic gas meter，a prototype with the improved flowing channel meets the design requirements and the measuring error is stable.

ultrasonic；gas velocity；CFX simulation；gas meter

A

1674-5124（2016）05-0135-05

10.11857/j.issn.1674-5124.2016.05.028

2015-11-17；

2015-12-12

翟義然（1973-），男，河北承德市人，高級工程師，博士，主要從事信息與信號處理方面的研究。