|
新聞詳情
測量低粘度流體介質(zhì)水平安裝金屬浮子流量計的仿真研究 對于水平安裝金屬浮子流量計 是一種傳統(tǒng)的變截面流量計,具有結(jié)構(gòu)簡單、工作可靠、壓力損失小且穩(wěn)定、可測低速流體介質(zhì)等諸多優(yōu)點(diǎn),廣泛應(yīng)用于測量高溫、高壓及腐蝕性流體介質(zhì)。目前,金屬管浮子流量計設(shè)計采用經(jīng)典浮子流量公式計算,其中流量系數(shù)α受到浮子形狀和來流雷諾數(shù)等多種因素的影響,需要實(shí)驗(yàn)標(biāo)定來確定。該設(shè)計方法需進(jìn)行大量的實(shí)驗(yàn),費(fèi)用昂貴,周期很長,且獲取的數(shù)據(jù)有限。為深入研究金屬管浮子流量計的工作機(jī)理,筆者利用計算流體力學(xué)方法對水平安裝金屬浮子流量計進(jìn)行數(shù)值仿真,通過構(gòu)造仿真模型為研究金屬管浮子流量計的機(jī)理建立數(shù)值實(shí)驗(yàn)平臺,優(yōu)化傳感器的結(jié)構(gòu);同時為金屬管浮子流量計提供了低成本、短周期的設(shè)計方法。 1 基本原理 1.1 水平安裝金屬浮子流量計的工作原理 圖1為水平安裝金屬浮子流量計的工作原理示意。在垂直的錐形管中,放置一阻力件,即浮子。當(dāng)流體自下而上穿過錐管時,受到浮子迎流體積阻擋而產(chǎn)生一個壓差,從而對浮子形成向上的作用力,同時由于流體本身的黏性,對浮子產(chǎn)生黏性力,當(dāng)這兩個力的合力大于浮子本身的重力時,浮子就會向上升,同時浮子與錐形管問的環(huán)形流通面積增大,流速減低,此時浮子對流體阻力作用減小.當(dāng)浮子受到的力達(dá)到平衡時,浮子就會停留在某一高度。傳統(tǒng)的水平安裝金屬浮子流量計設(shè)計采用經(jīng)典浮子流量公式計算,即 式中:qv為水平安裝金屬浮子流量計的體積流量;α為流量系數(shù);D0為浮子非常大迎流面的直徑;h為浮子在錐管中的垂直位置;φ為錐形管夾角;Vf為浮子體積;ρf為浮子材料密度;ρ為流體密度;Af為浮子垂直于流向的非常大截面面積。 1.2 計算流體力學(xué)及其控制方程 計算流體力學(xué)(computational fluid dynamics,CFD)是一門利用計算機(jī)求解描述流體流動規(guī)律的控制方程組技術(shù),涉及到流體力學(xué)、計算方法和計算機(jī)圖形處理等技術(shù)。為了簡便,以不可壓縮湍流流動為例寫出筆者使用的k-ε模型通用形式的流體控制方程。在直角坐標(biāo)系中,流動可由連續(xù)性方程和雷諾時均N-S方程描述,即 水平安裝金屬浮子流量計 連續(xù)方程 雷諾方程 式中:ui為平均速度;P為平均壓力;ν和νt分別為水平安裝金屬浮子流量計分子黏性系數(shù)和渦黏性系數(shù)。 對高數(shù)湍流,渦黏性系數(shù)的計算公式為 νT=Cμk2/ ε (4) k和ε由水平安裝金屬浮子流量計各自的輸運(yùn)方程得到,對高Re數(shù)問題有K方程 ε方程 式中:Gk為湍動能產(chǎn)生項(xiàng) D?為平均應(yīng)變率張量。 筆者的仿真模型建立在已有仿真實(shí)驗(yàn)的基礎(chǔ)上,并針對以往模型中存在的入口速度剖面為等值面這一不合理設(shè)置,給出合理的圓管入口速度剖面;同時按照實(shí)際流量計的構(gòu)造,在流量計的入口與出口處加入導(dǎo)流架,從而提高了水平安裝金屬浮子流量計仿真結(jié)果的準(zhǔn)確性。 2 模型建立和數(shù)值計算 針對25mm口徑(DN25)www.shhzy3.cn?水平安裝金屬浮子流量計內(nèi)部機(jī)械結(jié)構(gòu)進(jìn)行流場計算和分析,金屬管浮子流量計結(jié)構(gòu)剖面如圖2(a)所示,金屬管浮子流量計采用1Cr18Ni9Ti材質(zhì)。流量計設(shè)計流量范圍為0.4~4.0 m3,量程比為10:1,為了檢測錐管中浮子的位置,在浮子中安裝有永久磁鐵,同時利用網(wǎng)格生成器建立了浮子流量計內(nèi)部結(jié)構(gòu)仿真模型,如圖2(b)所示,實(shí)驗(yàn)所用的低黏度流體介質(zhì)是在常溫下密度為998.0 kg/m3、黏度為0.001002 kg/m.s)的水。www.shhzy3.cn 水平安裝金屬浮子流量計內(nèi)部流場是高雷諾數(shù)完全發(fā)展湍流流動,故采用湍流模式理論中標(biāo)準(zhǔn)一模型來計算。 2.1 水平安裝金屬浮子流量計入口速度剖面設(shè)定 入口速度由經(jīng)驗(yàn)性光滑圓管湍流速度分布指數(shù)公式表示,即 式中:ur為入口截面半徑為r處的速度;R為圓管入口處半徑;m為常數(shù),由雷諾數(shù)的大小決定,見表1;umax為圓管入口中心處的非常大速度。 式中:q為入口流量;A為金屬管入口的計算面積,對本研究,A=2πR2m2/[(2m+1)(m+1)]。 表1 常數(shù)m與雷諾數(shù)的關(guān)系
為使構(gòu)建模型能滿足不同口徑水平安裝金屬浮子流量計的仿真結(jié)果,引入流量修正因子Cq,通過修正入口流量來適應(yīng)多種口徑水平安裝金屬浮子流量計流量的仿真模型。 流量的修正公式如下: qc=Cqq (9) 式中:q為實(shí)際流量;qc為修正流量;Cq為流量修正因子,對于25mm口徑的金屬管浮子流量計,Cq=1.13。 根據(jù)式(8)~式(11)編寫速度分布函數(shù)。圖3是q=4.8 m3/h時金屬管浮子流量計入口速度剖面的仿真。圖中色標(biāo)顏色由冷色調(diào)到暖色調(diào)表示速度由小到大,從圖中可以清楚地看到從邊壁到中心的速度是由小到大的非線性分布。 2.2 其他邊界條件的設(shè)定 1)水平安裝金屬浮子流量計壁面條件的設(shè)置 由于水平安裝金屬浮子流量計內(nèi)表面采用不銹鋼材料,根據(jù)壁面特征設(shè)定:粗糙高度Ks=0.04;粗糙常數(shù)Cs=1。 2)出入口的湍流參數(shù) 湍動能 耗散率 式中:uav為入口平均速度;Cμ為常數(shù);L為金屬管浮子流量計內(nèi)部特征尺度。 3 水平安裝金屬浮子流量計仿真結(jié)果及分析 圖4給出了不同入口流量下,仿真計算得到的水平安裝金屬浮子流量計內(nèi)部截面上的速度偽色圖。可以清楚地看到水平安裝金屬浮子流量計中流體在浮子周圍以及出入口的速度分布。隨著流量的增加,浮子在管中的位置升,浮子與管道之間環(huán)隙變大,流體在管中的速度分布也隨之發(fā)生明顯變化。通過觀察偽色圖中的速度分布,并根據(jù)流體力學(xué)基本原理,可以初步判斷出計算所得結(jié)果是合理的。 根據(jù)物理實(shí)驗(yàn)結(jié)果,仿真計算了水平安裝金屬浮子流量計在錐管中某一高度處流體對它表面的壓力、黏性力和浮子重力Gf 3項(xiàng)的合力,仿真結(jié)果中浮子所受合力結(jié)果見表2,合力向上為正方向。 表2 浮子所受合力Ff的仿真數(shù)據(jù)以及數(shù)值計算與物理實(shí)驗(yàn)所得流量的比較
由表2可知計算所得浮子在各位置的合力接近于零,趨于平衡,為驗(yàn)證仿真模型所能達(dá)到的仿真精度,根據(jù)浮子只有在所受合力為零時才能平衡,又計算了浮子在各個高度保持平衡,即合力為零時,所對應(yīng)的入口流量。表2還列出了數(shù)值計算與物理實(shí)驗(yàn)所得流量的比較,并計算出仿真實(shí)驗(yàn)的滿度誤差占,即 式中:qs為數(shù)值模擬流量,m3/h;qp為物理實(shí)驗(yàn)流量。 由表2可知,數(shù)值計算與物理實(shí)驗(yàn)所得流量的非常大滿度誤差為5.4695%,平均滿度誤差為2.4731%。利用本模型對50mm、80mm口徑水平安裝金屬浮子流量計仿真數(shù)據(jù)與物理實(shí)驗(yàn)數(shù)據(jù)的比較結(jié)果也相當(dāng)接近。表明筆者構(gòu)建的模型得到令人滿意的結(jié)果。 4 結(jié)語 通過計算流體力學(xué)方法對金屬管浮子流量計進(jìn)行仿真,并對流場進(jìn)行了定量分析,計算出浮子在金屬管中不同垂直位置的受力大小及受力平衡時所對應(yīng)流量。仿真實(shí)驗(yàn)結(jié)果與物理實(shí)驗(yàn)數(shù)據(jù)相比,非常大滿度誤差為5.4695%,平均誤差為2.4731%,表明筆者構(gòu)建的金屬管浮子流量計數(shù)值模型能夠滿足低黏度介質(zhì)金屬管浮子流量計設(shè)計要求,為金屬管浮子流量計傳感器結(jié)構(gòu)進(jìn)一步優(yōu)化提供了可靠的數(shù)值仿真平臺。 |