摘要:　通過對適合安裝于水平管道的特殊結(jié)構(gòu)的金屬管浮子流量計三維湍流流場的數(shù)值仿真及實驗研究提出一種基于計算流體力學的流量傳感器設(shè)計方法。流場仿真所需的模型采用GAMBIT來建立,通過FLUNT軟件進行仿真,仿真過程中利用受力平衡控制計算精度。數(shù)值仿真結(jié)果和實驗結(jié)果比較當浮子受力平衡度誤差為9. 5%時,流量誤差為0. 944%,證實了仿真結(jié)果的準確性,同時利用流場仿真信息對流量傳感器模型做了進一步的優(yōu)化。

 

　　1　引　言

 

　　金屬管式浮子流量計是一種傳統(tǒng)的差壓式流量計,為了適應(yīng)部分管道的特殊要求,本文設(shè)計研究的浮子流量計是左進右出型的,其測量原理與經(jīng)典的豎直型浮子流量計相同,但它是一種可以安裝于水平管道的特殊結(jié)構(gòu)的浮子流量計。

 

　　一般對浮子流量計的經(jīng)典研究是根據(jù)伯努利方程進行的[1],在推導浮子流量計流量測量公式時忽略了粘性應(yīng)力項,而該項的作用實際上是存在的;傳統(tǒng)流量計的設(shè)計要通過實驗來檢驗和修正設(shè)計圖紙,這樣不僅延長了設(shè)計周期而且增加了設(shè)計成本。

 

　　基于上述兩點原因,在設(shè)計水平式安裝浮子流量計時為了深入了解浮子流量傳感器的工作機理,引入了計算流體力學,即CFD[2]技術(shù),對傳感器流場進行數(shù)值模擬,通過對仿真及實驗結(jié)果進行分析來評價初樣設(shè)計,優(yōu)化流量傳感器的結(jié)構(gòu)參數(shù),使流量傳感器的設(shè)計更加精確,并提高了設(shè)計效率。

 

        2　水平式安裝金屬管浮子流量計的原理

 

　　2.1　檢測原理(圖1)

 

　　

 

　　子位置;———錐形管錐半角;Vf———浮子體積;

 

　　ρf———浮子材料密度;ρ———流體密度;Af———浮

 

　　子垂直于流向的最大截面積;D0———浮子最大迎流

 

　　面的直徑;Dh———浮子平衡在h高度時錐形管的直

 

　　徑;df———浮子最大直徑。

 

　　2.2　模型建立及其設(shè)計要求

 

　　浮子流量計傳統(tǒng)的設(shè)計方法是建立在式(1)的基礎(chǔ)之上,在該方程中流量系數(shù)α是一個受很多因素影響的變量。對于本文所研究的水平式安裝浮子流量計,測量介質(zhì)為20℃的水,設(shè)計要求流量測量范圍1~10m3/h,量程比為10∶1,行程50mm,其流量系數(shù)α的經(jīng)驗值為0.9~1.0。浮子位于41mm高處的傳感器三維流場模型如圖2所示。

 

　　2.3　計算精度的控制

 

　　利用浮子組件受力平衡來控制計算精度。在FLUNT的受力分析報告中會提供指定壁面所受到的凈壓力Fy↑和粘性摩擦力Fm↑以及這兩個力的合力Ff↑。這三個力遵循下面的公式:

 

　　

 

　　3　數(shù)值仿真

 

　　3.1　網(wǎng)格劃分及邊界設(shè)定

 

　　針對傳感器的流場模型,選擇三角形/四面體網(wǎng)格來進行網(wǎng)格劃分。如圖3所示為水平式浮子流量計浮子位于41mm高時的軸向網(wǎng)格剖分圖。

 

　　在進行邊界的設(shè)定過程中設(shè)定速度入口、壓力出口,并將導桿壁面設(shè)定為float.wall1,浮子壁面設(shè)定為float.wall2,除浮子組件和導向環(huán)外的空間設(shè)定為fluid。

 

　　

 

 

 

3.2　FLUNT計算條件

 

　　模型建好以后輸出.msh文件,利用FLUNT進行仿真,FLUNT中相應(yīng)計算條件如表1所示。其中流體介質(zhì)的屬性:密度998. 2 kg/m3,粘度0.001 003 kg/ms,定壓比熱4 182 J/kg·K,熱導率　　0.6W /m·k。速度入口采用的是平均速度,針對浮子位于41mm高的模型計算達到平衡時的入口條件,如表2所示。

 

　　

 

　　3.3　仿真過程

 

　　水平式浮子流量計三維流場的仿真過程如圖4所示。該過程需要解釋的幾點如下所示:

 

　　(1)因每個模型入口流速的準確值未知,是根據(jù)經(jīng)典流量公式計算的一個假設(shè)的流量,因此仿真最終結(jié)束的判斷依據(jù)為浮子受力平衡的程度,即通過檢查仿真結(jié)果,對浮子進行受力分析,距離受力平衡點誤差小于10%時,認為達到計算精度,仿真計算結(jié)束。當誤差大于10%,首先考慮改進該模型的網(wǎng)格精度,如圖4中的左側(cè)方案1;當網(wǎng)格精度改進到一定程度后誤差仍大于10%,可修正入口條件(主要指入口流速,其余條件可相應(yīng)計算調(diào)整),如圖4中的右側(cè)方案2,直到滿足計算精度。

 

　　(2)利用SIMPLE算法計算時,每次計算迭代次數(shù)為500次,當不足500次SIMPLE算法就已經(jīng)達到收斂精度(10-4)時,程序自動結(jié)束,此時可檢查計算結(jié)果;當?shù)螖?shù)大于500次仍未收斂時,停止計算,此時需重新檢查網(wǎng)格狀況和邊界設(shè)定,進行網(wǎng)格的合理剖分和邊界的合理設(shè)定。實踐證明,網(wǎng)格布置的恰當與否會直接影響收斂速度和收斂結(jié)果,不合理的網(wǎng)格布置將導致計算發(fā)散或者結(jié)果不正確。

 

　　

 

　　(3)迭代前首先打開監(jiān)視器,監(jiān)視X、Y、Z三個方向的流速以及k方程和ε方程的收斂狀況,實踐證明,即使未達到預(yù)計的迭代次數(shù),若在監(jiān)視器中已出現(xiàn)明顯的發(fā)散現(xiàn)象,可強行中止本次計算。

 

　　4　仿真結(jié)果及實驗結(jié)果分析

 

　　4.1　壓力場分析(圖5、圖6)

 

　　比較壓力的絕對值可以看到:浮子底部左右壓力不對稱,這種不對稱現(xiàn)象的存在使得流量比較大時浮子會出現(xiàn)抖動。

 

　　4.2　速度場分析(圖7、圖8)

 

　　據(jù)圖分析如下:

 

　　(1)據(jù)顏色分辨出環(huán)隙流通面積最小處及下游靠近錐管壁的流場速度最大,前者是流通面積減小導致速度增大,后者則是因為流場方向的改變而引起的,特別是此處可能產(chǎn)生旋渦,導致有效流通面積減小,流體被擠向管壁,使得此處速度增大。

 

　　(2)流場下游,外直管左下角速度較小,主要是因為流場的出口在右邊,由于出口壓力小,流體流動都趨向出口。

 

　　(3)浮子的最小截面處,流場速度存在較大的變化。

 

　　

 

 

 

4.3　浮子受力定量分析

 

　　從FLUNT的受力報告中可以得到如表3所示數(shù)據(jù),根據(jù)設(shè)計初樣給出的浮子材料及尺寸結(jié)構(gòu),可得浮子重力為5.995 146N。根據(jù)仿真結(jié)果,浮子在Z方向上的合力為5.425 312 7N。根據(jù)受力平衡度誤差分析公式可得,|Ef|=9.5%,小于設(shè)定值10%,認為浮子受力達到平衡。

 

　　

 

　　

 

　　

 

　　4.4　物理實驗及結(jié)果分析

 

　　為了進一步驗證傳感器流場仿真結(jié)果,需要進行物理實驗。按照設(shè)計圖紙加工設(shè)計模型,加工完后,配上流量顯示儀表,在標準裝置上進行標定。標定方法利用標準表法,標準表選擇渦輪流量計(精度0.5級)。結(jié)合仿真流量數(shù)據(jù)、物理實驗數(shù)據(jù)與根據(jù)浮子流量經(jīng)典測量公式得到的設(shè)計流量數(shù)據(jù)進行比較可以得到表4。

 

　　

 

　　5　大流量下流量傳感器結(jié)構(gòu)的優(yōu)化及改進結(jié)構(gòu)后的仿真

 

　　由上述對壓力場的分析可知浮子組件受力不平衡,物理實驗也表明在大流量下會出現(xiàn)流量計振動的現(xiàn)象,這是由于傳感器流場出現(xiàn)了變化。從流場的速度分布圖可以看出,浮子組件的右邊速度特別大,其原因有前流場引起的,也有后流場的因素,由于傳感器的出口在右邊,所以流體有向右邊流的趨勢。另外,由于浮子組件前直管段有個直角彎,容易產(chǎn)生二次流,對浮子組件的受力也有很大的影響。

 

　　所以,要減弱振動,解決的根本方法就是改變傳感器結(jié)構(gòu)參數(shù)、優(yōu)化流場、使浮子左右受力差盡量減小。根據(jù)上述分析下面對水平式流量傳感器的結(jié)構(gòu)提出幾點優(yōu)化方案:

 

　　(1)加入整流器,以消除或減小旋渦的產(chǎn)生,同時調(diào)整流速的分布狀況。

 

　　(2)將前流場的直管連接改為彎管連接,減少旋渦的產(chǎn)生,順滑流體的流動,使傳感器有比較平穩(wěn)的前流場。

 

　　(3)延長前直管段。這里提及的直管段指錐管前的垂直直管段,這也是為了使流體在通過整流器后有比較長的緩和段,使流場接近充分發(fā)展的流速分布。

 

　　改進結(jié)構(gòu)后的仿真結(jié)果如圖9、10所示,據(jù)圖分析如下:

 

　　

 

　　

 

　　(1)改進結(jié)構(gòu)后流場的壓力分布得到改善,浮子組件受力接近平衡,但是,由于整流器的引入,導致了整流器前后壓差增大,帶來比較大的壓損。

 

　　(2)改進結(jié)構(gòu)后流場的速度分布比較均勻,特別是使浮子組件周圍沒有太大的速度差,同樣由于整流器的使用,也使浮子組件的前流場更加復(fù)雜。

 

　　通過物理實驗也證實了這幾種優(yōu)化方案可以有效的減少浮子左右受力差,穩(wěn)定浮子。

 

　　6　結(jié)　論

 

　　由上述數(shù)據(jù)分析可知,對于浮子在41mm高處時的三維湍流流場進行仿真可得到設(shè)計要求的流量上限值。此位置處浮子受力平衡度誤差為9. 5%,傳感器物理實驗獲得的示值刻度流量與通過湍流數(shù)值模擬進行流場仿真實驗獲得的仿真流量值較為接近,仿真流量誤差為0. 944%。本文利用浮子受力平衡度誤差法確定仿真計算精度獲得了較為理想的效果,即仿真過程無需過分強調(diào)浮子受力平衡度誤差的減小,仿真流量誤差即可得到令人滿意的結(jié)果。

 

　　理論分析和實驗研究表明,這種設(shè)計方法不僅可以進一步地理解流體流動的機理和浮子流量計的測量原理,而且使流量傳感器的設(shè)計進一步得到優(yōu)化,使流量測量的靈敏度和精確度得到明顯的提高。

 

　　此外,對流場的數(shù)值仿真與實驗研究也是分析、解決流量計其它問題的一種有效方法。目前基于這種方法設(shè)計的水平式金屬管浮子流量計已經(jīng)投入市場,現(xiàn)場反饋這種流量計性能穩(wěn)定,精度可靠。