摘要:為了研究孔板流量計(jì)在動(dòng)態(tài)非穩(wěn)定流態(tài)或振蕩流態(tài)下的瞬時(shí)壓力-流量特性，理論分析了孔板前后的旋渦城大小隨流速變化是引起孔板進(jìn)出口瞬時(shí)流量差的主要原因.借助CFD數(shù)值解析方法,建立孔板模型,并在模型入口加載某一頻率下的正弦流速，對(duì)孔板流量計(jì)在振蕩流態(tài)下的瞬時(shí)壓力-流量特性進(jìn)行分析.結(jié)果表明:當(dāng)孔板流量計(jì)處于低頻振蕩流動(dòng)狀態(tài)時(shí),孔板兩端差壓也處于周期振蕩狀態(tài),差壓與節(jié)流孔瞬時(shí)流量同頻不同相，差壓幅值隨入口流速振幅增大而線性增大，且線性增長(zhǎng)系數(shù)與振蕩頻率相關(guān);孔板的入口與出口存在周期波動(dòng)的瞬時(shí)流量差,振蕩頻率越大或入口流速峰值越小,瞬時(shí)流量差的波動(dòng)越小，由于相位滯后和瞬時(shí)流量差的存在，使孔板流量計(jì)的測(cè)量流量與實(shí)際出口流量之間存在偏差.振蕩頻率越大，偏差也越大.
孔板流量計(jì)因其結(jié)構(gòu)簡(jiǎn)單、耐用而成為目前國(guó)際上標(biāo)準(zhǔn)化程度最高、應(yīng)用最為廣泛的一種流量計(jì),因此研究節(jié)流孔的流量特性,對(duì)提高孔板流量計(jì)測(cè)量不確定度的認(rèn)識(shí)具有很重要的意義.孔板流量計(jì)通過(guò)測(cè)量壓差進(jìn)而獲得流量.當(dāng)液流經(jīng)過(guò)節(jié)流孔,流束縮小,流速變大并伴隨著較大的壓力降.流束的最小橫斷面出現(xiàn)在實(shí)際縮口的下游,稱(chēng)為縮流斷面.在縮流斷面處,壓力最低.壓降的產(chǎn)生是由于在孔板的兩側(cè)側(cè)面出現(xiàn)回流區(qū)及旋渦域,造成較大的內(nèi)部紊流和能量損耗的結(jié)果[1-2].旋渦域的大小取決于流動(dòng)雷諾數(shù),隨著雷諾數(shù)的增大,渦旋強(qiáng)度增加[3].
流體力學(xué)中對(duì)孔口恒定出流的描述為孔板結(jié)構(gòu)的設(shè)計(jì)提供了理論依據(jù).但實(shí)際應(yīng)用中，由于外界激勵(lì)引起的壓力波動(dòng),圓管內(nèi)流體常處于動(dòng)態(tài)非穩(wěn)定流態(tài)或振蕩流態(tài)[4],孔板流量計(jì)內(nèi)部流場(chǎng)結(jié)構(gòu)變化極為復(fù)雜，因此,計(jì)量孔板的瞬時(shí)流量特性往往與理論分析結(jié)果存在偏差因此,有必要對(duì)孔板在非穩(wěn)定流態(tài)下的流量特性進(jìn)行研究.
通過(guò)Fluent流體仿真程序,對(duì)不同節(jié)流孔直徑比的孔板,以水為介質(zhì)在振蕩流態(tài)下的流動(dòng)過(guò)程進(jìn)行仿真,對(duì)其瞬時(shí)壓力-流量特性進(jìn)行分析.
1理論分析
通常在特定測(cè)壓位置和特定流體參數(shù)情況下,根據(jù)流體流動(dòng)的連續(xù)性方程和伯努利方程可推導(dǎo)出孔板前后差壓△p與流經(jīng)節(jié)流孔的體積流量QY滿足以下函數(shù)關(guān)系[5],即
式中:C為流出系數(shù);ρ為流體密度;β為節(jié)流孔的直徑比(β為節(jié)流孔直徑d與圓管內(nèi)徑D的比值,即β=d/D);sign為符號(hào)函數(shù)*.
圓管進(jìn)口流量可計(jì)算公式為
孔板流量計(jì)通過(guò)測(cè)量節(jié)流孔兩端差壓進(jìn)而獲得節(jié)流孔流量QV.對(duì)于不可壓縮的定常流，圓管進(jìn)口流量Qin和出口流量Qout與節(jié)流孔流量QY相等,聯(lián)立以上方程可得節(jié)流孔兩端差壓與人口流速的關(guān)系表達(dá)式為
由式(3)可知,孔板兩端差壓也呈周期性波動(dòng),.其振蕩頻率與孔板人口流速振蕩頻率相同.
孔板前后存在旋渦域.旋渦域的大小占據(jù)圓管空間,液體在旋渦域停留,不流向下游管道.旋渦域增大,則流向圓管出口的液流減少.由于孔板前后游渦旋強(qiáng)度隨流動(dòng)雷諾數(shù)增大而增大,即隨流速增加旋渦域變大3]，.在振蕩流態(tài)下,旋渦域大小隨入口流速變化也表現(xiàn)為周期性變化狀態(tài),變化頻率與流速振蕩頻率相同.因此，在某一極短時(shí)間段內(nèi),旋渦域的體積變化量表現(xiàn)為圓管進(jìn)、出口的瞬時(shí)流量之差.對(duì)于不可壓縮的非定常流,人口瞬時(shí)流量Qin與出口瞬時(shí)流量Qout和節(jié)流孔瞬時(shí)流量QY滿足以下關(guān)系,即
式中:△Q表示圓管進(jìn)出口瞬時(shí)流量差.瞬時(shí)流量差的存在，使孔板流量計(jì)實(shí)際測(cè)量流量Qv與出口瞬時(shí)流量Qout之間不可避免存在偏差.
事實(shí)上,由于節(jié)流孔的壓降作用,當(dāng)孔板下游壓力低于液體飽和蒸氣壓以下，氣泡將在下游管道產(chǎn)生,形成閃蒸現(xiàn)象.當(dāng)壓力上升,氣泡破裂瞬間產(chǎn)生局部空穴,高壓液體重新流向這些空間.顯然,氣泡和空穴占據(jù)了下游管道空間,使進(jìn)、出口流
量Qm與Qout,存在差異,出現(xiàn)瞬時(shí)流量差實(shí)際孔板流量計(jì)使用過(guò)程中避免閃蒸和空穴現(xiàn)象的出現(xiàn),故文中對(duì)其影響不做表述.
為進(jìn)一步探究孔板的瞬時(shí)壓力流量特性,文中以上述理論分析為基礎(chǔ)，結(jié)合有限元分析思想,對(duì)孔板流量計(jì)在低頻微幅振蕩流態(tài)下的壓力流量特性進(jìn)行了分析.
2期修仿真
2.1控基文圖
文中選用RNGk-ε湍流模型對(duì)孔板的流量特性進(jìn)行模擬.該模型的控制方程分別為連續(xù)性方程
上述式中:xi,xj分別為縱向和橫向坐標(biāo);ui,uj分別為縱向和想象的速度分量;p為流體壓力;v為流體運(yùn)動(dòng)黏度;vt為流體渦流黏度,vt=Cμk²/Ɛ,其中k為湍動(dòng)能,Ɛ為湍動(dòng)耗散率,Cμ=0.085.
模型邊界條件包括速度人口、壓力出口、無(wú)滑移壁面邊界,在近壁面區(qū)域采用標(biāo)準(zhǔn)壁面麗數(shù)進(jìn)行處理.采用軸對(duì)稱(chēng)邊界,即模型對(duì)稱(chēng)軸的徑向速度為0.在求解離散方程組和壓力速度耦合時(shí)選擇了SIM-PLE算法,動(dòng)量和湍流動(dòng)能分別采用的是二階迎風(fēng)與一階迎風(fēng)差分格式.
2.2仿真文收
利用孔板模型的軸對(duì)稱(chēng)性的特征,在圓柱坐標(biāo)系下建立它們的1/2實(shí)體模型,取壓方式采用D-D/2取壓其計(jì)算域如圖1所示.孔板.上游直管段長(zhǎng)度為20D,充足的上游管長(zhǎng)能夠確保液流在孔板上游為充分發(fā)展的湍流流動(dòng).模型具體尺寸,其中D=12.3mm,β=0.247,Lu=246mm,Lt=494mm,t=2mm.
為了表現(xiàn)孔板前后的流場(chǎng)變化情況,首先在壁面附近劃分邊界層網(wǎng)格,邊界層第一次厚度為0.1mm,共10層,高度增長(zhǎng)因子1.1.其次,為了提高孔板附近的計(jì)算精度,對(duì)靠近孔板部分的網(wǎng)格進(jìn)行局部加密,離節(jié)流孔越遠(yuǎn),網(wǎng)格越稀疏最后,利用結(jié)構(gòu)化網(wǎng)格生成方式劃分其余部分網(wǎng)格.
文中所選用的流體介質(zhì)為常溫狀態(tài)下的水.人口流速設(shè)定為某一-頻率下的正弦流動(dòng)u=uarg+uamp·sin(2πƒt)，選擇不同平均流速u(mài)arg、流速振幅uamp和振蕩頻率ƒ參數(shù)作為節(jié)流孔的人口流速,具體參數(shù)見(jiàn)表1.利用UDF功能將該自定義速度函數(shù)加載在模型的速度人口.
3網(wǎng)出第日
3.1振蕩差絡(luò)
通過(guò)后處理后可以觀察到,當(dāng)人口流速為某一頻率下的正弦流動(dòng)時(shí),孔板兩端將出現(xiàn)與人口流速頻率相同的振蕩差壓.如圖2所示，節(jié)流孔瞬時(shí)流量與差壓振蕩頻率相等且具有固定的相位滯后.相位滯后意味著測(cè)量壓差不能反映當(dāng)時(shí)的流量情況.此外，由于壓差測(cè)量裝置的動(dòng)作時(shí)限,測(cè)量壓差滯后,不能及時(shí)反映瞬時(shí)壓差的變化因此，在振蕩流態(tài)下,孔板流量計(jì)對(duì)瞬時(shí)流量的測(cè)量存在不確定性.
圖3為人口流速振幅與差壓幅值的關(guān)系.對(duì)于同一振蕩頻率的入口流速,孔板兩端差壓幅值隨人口流速振幅增大而線性增大,但其線性增長(zhǎng)系數(shù)與振蕩頻率有關(guān).
從圖4中可以看出當(dāng)人口流速振幅一定時(shí),節(jié)流孔兩端差壓的振蕩幅值隨振蕩頻率的增大而增大.差壓幅值與振蕩頻率存在近似一次線性關(guān)系.
在孔口恒定出流情況下,測(cè)量流量與實(shí)際節(jié)流孔流量Qv相同.而在振蕩流態(tài)下,差壓幅值隨振蕩頻率線性增大,則測(cè)量流量幅值越大,與實(shí)際節(jié)流孔流量的偏差也越大.
從圖5中可以看出,平均入口流速的變化,對(duì)壓力幅值的影響幾乎可以忽略.
3.2瞬時(shí)流量差
在振蕩流態(tài)下,孔板前后回流區(qū)和旋渦域的大小隨人口流速變化不斷改變，導(dǎo)致進(jìn)出口流量存在.瞬時(shí)流量差△Q,如圖6所示.瞬時(shí)流量差表現(xiàn)為復(fù)雜的周期性波動(dòng),其波動(dòng)周期與差壓振蕩周期相同,相位介于瞬時(shí)流量和差壓兩者之間,且稍滯后于振蕩差壓.當(dāng)差壓增大至峰值點(diǎn)時(shí)，瞬時(shí)流量差趨向其波峰,并在到達(dá)峰值點(diǎn)后反向階躍.
為研究人口流速各參數(shù)對(duì)瞬時(shí)流量差的波動(dòng)特性影響,對(duì)仿真記錄的瞬時(shí)流量差數(shù)據(jù)作方差分析和極差分析,以此描述瞬時(shí)流量差的波動(dòng)情況.瞬時(shí)流量差的極差和方差與振蕩頻率關(guān)系如圖7所示.當(dāng)入口平均流速和流速振幅不變時(shí),瞬時(shí)流量差的極差和方差隨人口流速的振蕩頻率增大而減小也即人口流速頻率越大,瞬時(shí)流量差的波動(dòng)程度越小，同時(shí)波動(dòng)的峰值也越小.
圖8為入口平均流速與瞬時(shí)流量方差及極差的關(guān)系.當(dāng)振蕩頻率和流速振幅相同時(shí),人口平均流速越大,瞬時(shí)流量差的方差和極差越大.圖9為瞬時(shí)流量差方差和極差與流速振幅關(guān)系.從圖中可以看出，當(dāng)人口平均流速相同時(shí),對(duì)于給定的振蕩頻率,瞬時(shí)流量差的方差和極差隨著流速振幅增大而增大.因此,當(dāng)人口流速峰值越大,瞬時(shí)流量差波動(dòng)也越大,瞬時(shí)流量差就越不穩(wěn)定.
4結(jié)論
當(dāng)孔板流量計(jì)所計(jì)量不可壓縮流體為低頻振蕩流動(dòng)狀態(tài)時(shí)，通過(guò)前述CFD分析,得到如下結(jié)論:
1)孔板兩端差壓為周期振蕩狀態(tài),差壓與節(jié)流孔瞬時(shí)流量同頻不同相.差壓幅值隨人口流速振幅增大而線性增大，且線性增長(zhǎng)系數(shù)與振蕩頻率相關(guān).
2)圓管人口與出口存在周期波動(dòng)的瞬時(shí)流量差,振蕩頻率越大或人口流速峰值越小，則瞬時(shí)流量差的波動(dòng)也越小
3)在振蕩流態(tài)下,由于相位滯后和瞬時(shí)流量差的存在,使孔板流量計(jì)的測(cè)量流量與實(shí)際出口流量之間存在偏差振蕩頻率越大,偏差也越大..
4)孔板流量計(jì)作為的流量計(jì)量的常用元件,該分析結(jié)果對(duì)孔板的結(jié)構(gòu)設(shè)計(jì)及系統(tǒng)的整體動(dòng)態(tài)特性研究具有重要意義. 本文來(lái)源于網(wǎng)絡(luò),如有侵權(quán)聯(lián)系即刪除！
