本文采用負(fù)壓噴嘴裝置對(duì)選取的DN80渦街流量計(jì)進(jìn)行重復(fù)性穩(wěn)定性試驗(yàn),該流量計(jì)流量范圍為70~700m3/h,參照檢定規(guī)程”選取Qmin、20%Qmax、40%Qmx、Qm四個(gè)流量點(diǎn)開展試驗(yàn)。LBM對(duì)DN80渦街流量計(jì)檢定流場的數(shù)值仿真計(jì)算,其速度根據(jù)試驗(yàn)中測得的流量點(diǎn)給定,數(shù)值模擬計(jì)算得到的渦街流量計(jì)流場信息,輸出升阻力變化信息。
通過流場后處理軟件Tecplot對(duì)其進(jìn)行分析,可以顯示渦街流量計(jì)數(shù)值仿真的流場信息，可以觀測到漩渦的脫落狀態(tài)分布以及流線速度特征，漩渦的脫落導(dǎo)致渦街流量計(jì)漩渦發(fā)生體表面的壓力分布和升/阻力系數(shù)的周期變化。圖3給出的是渦街流量計(jì)管道流量為640m3/h,Re為50050的流場流線分布及管向速度云圖;圖4給出的是升力系數(shù)在格子時(shí)間尺度下的周期變化。渦街流量計(jì)儀表的檢定過程中，關(guān)注的主要是儀表系數(shù),其正比于儀表脈沖系數(shù)(或漩渦脫落頻率),反比于儀表流量。于是,本文數(shù)值仿真關(guān)注的主要是渦街流量計(jì)的漩渦脫落頻率,該頻率與渦街漩渦發(fā)生體的升力系數(shù)變化周期相關(guān),可以通過對(duì)升力系數(shù)變化信息的分析得到漩渦脫落頻率。由于本文是基于LBM的數(shù)值模擬，,仿真結(jié)果也是立足于格子尺度,漩渦脫落頻率要經(jīng)過當(dāng)?shù)馗褡訒r(shí)間比例因子1,的變換獲得。利用Matlab傅里葉變換對(duì)升力系數(shù)變化周期進(jìn)行頻域變換分析,可以得到渦街流量計(jì)漩渦發(fā)生體的漩渦脫落周期,經(jīng)過比例因子l,的轉(zhuǎn)換獲得當(dāng)?shù)劁鰷u脫落周期,進(jìn)而求得當(dāng)?shù)劁鰷u脫落頻率。
表1是本文數(shù)值仿真得到的渦街漩渦脫落頻率與檢定裝置實(shí)測數(shù)據(jù)的對(duì)比,發(fā)現(xiàn)計(jì)算結(jié)果與實(shí)驗(yàn)
結(jié)果相比誤差絕對(duì)值在10%以內(nèi),考慮到本文沒有考慮流場三維效應(yīng)以及笛卡爾網(wǎng)格處理模型邊界分辨率問題,該方法作為基礎(chǔ)方法能夠用于渦街流量計(jì)仿真工程應(yīng)用。 |
