王曉華

(中國石油遼河油田公司鉆采工藝研究院,遼寧盤錦 124000)

基于CFD的射流流量計的性能研究

王曉華

(中國石油遼河油田公司鉆采工藝研究院,遼寧盤錦 124000)

本文應(yīng)用SimpleC算法和Realizable湍流模型對垂直式天然氣射流表內(nèi)部結(jié)構(gòu)之一的震蕩器進行了壓力損失特性方面的數(shù)值模擬研究。分析了不同入口流量值對應(yīng)的震蕩器反饋通道上流體質(zhì)點的壓力脈動情況。發(fā)現(xiàn)在一定的入口流量范圍內(nèi),流速與震蕩器上下反饋通道的壓力震蕩頻率存在很明顯的線性關(guān)系。給出了該震蕩器某一截面上的流場速度矢量圖、靜壓力等值線圖和速度等值線圖,該數(shù)值模擬研究結(jié)果對垂直式射流表的優(yōu)化設(shè)計研究具有重要意義。

射流表 內(nèi)流場 數(shù)值模擬 振蕩器

射流技術(shù)是上世紀(jì)60年代發(fā)展起來的一門新興技術(shù)。射流流量計與其他流體振動式流量計相比一個突出的優(yōu)點是能在遠低于后者要求的雷諾數(shù)下穩(wěn)定起振。例如渦街流量計，當(dāng)管道雷諾數(shù)降至10000左右時便無法運行。但對于射流流量計，通過采取適當(dāng)?shù)慕Y(jié)構(gòu)，可以使其運行下限雷諾數(shù)降至120甚至100以下[5-10]。傳統(tǒng)的節(jié)流式和旋葉式流量計在微小流量計量的穩(wěn)定性和精度方面都存在一定的局限性。本課題研究的流量計，在整體上采用了立體結(jié)構(gòu)設(shè)計，在一定程度上消除流體振蕩引起的干擾信號，從而起到濾波的效果。對于射流元件，采用了共鳴腔設(shè)計，放大壓力脈動信號[11]，易于信號的檢測。同時采用了新的振蕩腔的設(shè)計，拓寬了流量計的量程比，提高了其工作的適應(yīng)性和穩(wěn)定性。因此，本課題研究的流量計具有量程較小、能測量微小流量且長期工作可靠，并且精度較高。

本文應(yīng)用Fluent流動分析軟件，計算垂直引流式流量計內(nèi)部流場，得到流量計內(nèi)部流場流動情況，分析射流振蕩過程，預(yù)測流量計整體能效，為進一步了解流量計的流動機理，改善流量計的氣動性能提供依據(jù)。

圖1 射流流量計三維模型

圖2 射流流量計的計算模型

1 射流流量計的幾何模型

本文利用SolidWorks軟件對射流流量計進行參數(shù)化建模。所建的射流流量計的三維造型如圖1所示。

圖3 射流流量計的壓分布圖

圖4 進口速度 v=0.5m/s 時震蕩層云圖

2 數(shù)值模擬

2.1 控制方程

對于氣體流量計，其氣動性能在很大程度上取決于流量計的結(jié)構(gòu)。由于流道形狀、哥氏力和粘性力的影響，流量計內(nèi)的氣體流動十分復(fù)雜[2]。一般認(rèn)為氣流在流量計內(nèi)的相對運動和在靜止元件內(nèi)的絕對運動為非定常流，而且震蕩器內(nèi)的氣體壓強變化不大，可忽略氣體的壓縮性。因此，流量計內(nèi)的流動可認(rèn)為是三維、非定常、不可壓縮流動。本文求解的流體動力學(xué)特性可以用流體力學(xué)基本方程描述：

連續(xù)方程：

動量方程：

圖5 不同速度下同一監(jiān)測點的頻率隨時間變化圖

圖6 不同流速下的震蕩頻率

式中， p是靜壓，ui，uj是流動速度分量，iF是質(zhì)量力，τij是應(yīng)力張量分量，定義為：

式中， μ是流體的湍流運動粘性系數(shù)。

計算中采用Realizable k-ε湍流模型。Realizable k-ε湍流模型中采用了新的湍流粘度公式，方程是從渦量擾動量均方根的精確輸運方程推導(dǎo)出來的，為目前工程上使用最為廣泛的湍流模型之一，并且較其它湍流模型計算更容易收斂，在分離流計算和帶二次流的復(fù)雜流動計算中的研究標(biāo)明，Realizable k-ε模型是所有 k-ε模型中表現(xiàn)最出色的湍流模型。

該模型滿足對雷諾應(yīng)力的約束條件，因此可以在雷諾應(yīng)力上保持與真實湍流一致，可以更精確地模擬各種流動包括旋轉(zhuǎn)均勻剪切流、包含有射流和混合流的自由流動、管道內(nèi)部流動、邊界層流動和帶有分離的流動等。Realizable k-ε湍流模型是兩方程模型，需要求解的變量為湍動能 k與湍動能耗散率ε，湍動能輸出方程為：

其中，

其中，Gk為由于平均速度梯度導(dǎo)致的湍動能 k生成項，其表達式為：

而且，ν=μ/ ρ，系數(shù)C1的表達式為：

(4)式和(6)式中Realizable k-ε的模式常數(shù)為

2.2 邊界條件

進、出口處給定恒定的壓力條件。入口處為質(zhì)量入口，出口處為壓力出口，出口壓力與外界大氣壓相等。壁面為無滑移邊界條件，近壁面的流動采用Launder和Spalding提出的標(biāo)準(zhǔn)壁面函數(shù)處理。本文采用“凍結(jié)轉(zhuǎn)子”模型即多重參考坐標(biāo)系(MRF)模型，以得到穩(wěn)定的解，該模型假定震蕩腔的交界面是穩(wěn)態(tài)的，且相對位置保持不變。物理模型及劃分好的網(wǎng)格模型如圖2所示。

3 計算結(jié)果分析

3.1 壓力場分布

從壓力分布圖3中可以觀察到反饋通道壓力與蜂鳴腔壓力相同。監(jiān)測點設(shè)置在反饋通道直管處，并且從仿真圖中可知上下反饋通道中的壓力互為相反數(shù)。

3.2 震蕩層流場

從圖4中可以發(fā)現(xiàn)，由于存在附壁效應(yīng)，在震蕩器中將產(chǎn)生流體震蕩現(xiàn)象。計算結(jié)果表明，射流的切換過程總是伴隨著一個低壓區(qū)域沿同側(cè)的側(cè)壁從前向后的移動過程。反饋通道中的流體由控制噴嘴射出，與主射流相混合，在二者的相互作用下，在側(cè)壁面前段形成一個近壁漩渦。漩渦位于低壓中心，此漩渦對射流向另一個側(cè)壁面偏轉(zhuǎn)起著推動的作用。同時，在分流劈的凹部，由于分流劈的作用也會出現(xiàn)渦流。所不同的是，該渦流將使射流產(chǎn)生誘導(dǎo)速度，導(dǎo)致射流穩(wěn)定地向側(cè)壁附著，反而抑制了主射流的擴散和偏轉(zhuǎn)。當(dāng)分流劈凹部的渦流起主導(dǎo)作用時，射流仍可穩(wěn)定的附著側(cè)壁面而不被切換。隨著噴嘴流量漸漸增大，側(cè)壁面附近的漩渦強度增強，并向側(cè)壁后端移動，從而增強了對主射流的切換作用，而分流劈凹部的渦流逐漸變?nèi)?，便推動主射流偏轉(zhuǎn)的渦流最終占據(jù)主導(dǎo)地位，導(dǎo)致主射流被切換并依附另一個側(cè)壁面。此過程周期進行，便發(fā)展成為以某個固定頻率發(fā)生的周期性震蕩的流動現(xiàn)象，將誘導(dǎo)震蕩腔中的壓力發(fā)生周期性變化。流體的震蕩頻率由反饋通道中的壓力震蕩頻率決定。在振動階段，流場的壓力分布隨著時間變化而變化。該仿真結(jié)果如圖5所示。

此次我們對同一物理模型采用了12個不同入口邊界條件。從38.4 L/min-384 L/min的流量值開始，在流場中的震蕩層設(shè)置了5個監(jiān)測點，監(jiān)測流場中壓力、速度等物性參數(shù)隨時間的變化規(guī)律，共計算了12個算例，得到了不同流速下的震蕩頻率。為了簡單起見，圖5中僅示出了五種不同流速下的壓力變化，從圖中可以看出在反饋通道中的壓力變化類似于一條正弦曲線。通過采用FFT方法，我們得到了其變化頻率(見圖6)。

入口流速由1～10m/s變化時，在震蕩層反饋通道的監(jiān)測點上壓力的頻率隨著入口流速的增大而增大。在一定范圍內(nèi)，頻率與流速之間具有較好的線性關(guān)系。

4 結(jié)語

本文對射流表內(nèi)部流場進行了數(shù)值模擬，研究發(fā)現(xiàn)流體流經(jīng)射流流量計振蕩腔時產(chǎn)生自振現(xiàn)象，通過監(jiān)測反饋通道中的壓力得到振蕩頻率，在一定的流量范圍內(nèi)，振蕩頻率與流量呈線性關(guān)系。因此，采用數(shù)值仿真的方法可對流量計的結(jié)構(gòu)和測量特性進行研究，可為盡可能地在線性范圍內(nèi)拓寬射流流量計的量程比而進行結(jié)構(gòu)設(shè)計及參數(shù)設(shè)定提供一定的依據(jù)。

[1]王福軍.計算流體動力學(xué)分析—CFD軟件原理與應(yīng)用[M].北京：清華大學(xué)出版社,2007．

[2]梁國偉,滕武昌.流量測量技術(shù)及儀表[M].北京：機械工業(yè)出版社,2005：52-56.

[3]王瑞金,張凱,王剛.FLUENT技術(shù)基礎(chǔ)與應(yīng)用實例[M].北京：清華大學(xué)出版社,2007.

[4]景思睿,張鳴遠.流體力學(xué)[M].西安：西安交通大學(xué)出版社,2001.

[5]傅新,王馳宇,謝海波,等.射流流量計的仿真與試驗研究[J].機械工程學(xué)報,2006,42(7)：24-29.

[6]于曉洋,王雁,郭華,等.射流流量傳感器原理與設(shè)計[J].哈爾濱理工大學(xué)學(xué)報,1997,2(4)：34-36.

[7]方婷,謝代梁,梁國偉.射流流量計的數(shù)值仿真研究[J].中國計量學(xué)院學(xué)報,2008,19(2)：114-118.

[8]周云波,林佳端.射流流量計的原理及發(fā)展[J].中國儀器儀表,1997,1：3-5.

[9]王立勇,王興才.射流流量計的原理與應(yīng)用[J].工業(yè)計量,2002,3：44-45．

[10]趙漢臣,佘文晟.用途廣泛的射流流量計[J].自動化儀表,2001,22(9)：20-23．

[11]謝代梁,邢玉雷,梁國偉.大量程比射流流量計的仿真研究[J].傳感技術(shù)學(xué)報,2007,20(7)：1655-1658.

[12]白亞磊,明曉.附壁射流元件的仿真研究[J].南京航空航天大學(xué)學(xué)報,2008,40(1)：32-36.

“十二五”國家重大專項“碳酸鹽巖、火成巖及酸性氣藏高效安全鉆完井技術(shù)”(No.2010zx05021-002)和國家自然科學(xué)基金“基于分形理論的裂縫性儲層水力壓裂機理及三維動態(tài)模擬”(No.50974029)。

王曉華(1961-),女,工程碩士,高級工程師,主要從事油田測試技術(shù)研究及應(yīng)用方向;

李瑋,男,1979年2月生,2010年獲東北石油大學(xué)博士學(xué)位,現(xiàn)為東北石油大學(xué)副教授,主要從事油氣井工程力學(xué)及巖石力學(xué)方面的教學(xué)和研究工作。E-mail：cyyping@sina.com。