超聲波流量計(jì)在測量過程中的(de)彎管誤差分析以(yi)🤩及修正研究

關鍵字：   超聲(shēng)波流量計   測(cè)量過程中   彎(wan)管誤差

一、本(běn)文引言 　 　 　

　超聲(shēng)波流量計 因爲具(jù)有非接觸測量 、計(jì)量準确度高、運行(háng)穩定、無壓力損失(shī)等諸多優點，目前(qian)怩在💜工業檢測領(ling)域有着廣泛的應(ying)用，市場對于相關(guān)産品的🈲需求💋十分(fen)地旺盛。伴随着上(shàng)個世紀 80年代(dai)電子技術和傳感(gan)器技術的迅猛發(fā)展，對于超聲波流(liu)量計的基礎研究(jiu)也在💜不斷地深👣入(ru)✏️，與此相關的各類(lèi)涉及到人們生産(chǎn)與生活的新産品(pin)也日新月異，不斷(duàn)出現。目前對于超(chāo)聲波流量計測量(liàng)精度的研究主要(yao)集中在 3個方(fāng)面：包括信号因素(sù)、硬件因素以及流(liú)場因素這三點。由(yóu)于超聲波流量計(jì)對流場狀态十分(fen)敏感，實際安裝現(xian)場的🔴流場不穩定(ding)會直接影響流量(liang)計的測量精度。對(dui)于超聲波流量計(jì)㊙️流場研究多采用(yòng)計算流體力學（ CFD）的方法，國内外(wài)諸多學者對超聲(shēng)波流量計在彎管(guǎn)流場情況下進行(hang)數值仿真，并進行(háng)了實驗驗證。以🔴往(wang)的💃🏻研究主要🐪是針(zhen)對規避安裝效應(ying)的影響。不過在一(yi)些中小口🧡徑超聲(shēng)波流量計的應🤩用(yòng)場合，因爲受到場(chang)地的限制，彎管下(xià)遊緩沖管道不足(zu)，流體在流經彎管(guǎn)後不能充分發展(zhan)，檢測精度受到彎(wan)管下💘遊徑向二次(cì)流分速度的極大(da)影響，安裝效應需(xū)要評估，并研究相(xiàng)應☔的補償方法。

　 　 　 　本研究采用 CFD仿真分析 90°單彎管下遊二次(ci)流誤差形成原因(yin)，并得出誤差的計(ji)算公💋式，定量地分(fen)析彎管下遊不同(tong)緩沖管道後，不同(tong)雷諾數下的二✉️次(cì)流誤差對測量精(jīng)度的影響，zui終得到(dào)誤差的修正規律(lü)。通過仿真發現，彎(wan)管出口處頂端和(he)💁底端的壓🐉力差與(yǔ)彎管二次流的強(qiáng)度有🈲關，提出在實(shí)🔴際測量中可通過(guò)測得此壓力差來(lái)對二次流誤差進(jin)行修正的方法。該(gai)研究可用于分析(xī)其他類型的超聲(shēng)波流量計的誤差(cha)🔞分析，對超聲波流(liu)量計的設計與安(an)裝具有重要意義(yi)。
二、測量原理(lǐ)與誤差形成
1.1 超聲波流量計測(ce)量原理
本研(yan)究針對一款雙探(tan)頭時差法超聲波(bo)流量計。時差法是(shì)利用🏃‍♂️聲脈沖波在(zài)流體中順向與逆(nì)向傳播的時間差(cha)來測量流體流速(sù)。雙探頭超聲波流(liú)量計原理圖如圖(tu) 1所示。
 

　 順(shùn)向和逆向的傳播(bo)時間爲 t1 和 t2 ，聲道線與管道(dào)壁面夾角爲 θ ，管道的橫截面積(ji)爲 S ，聲道線上(shang)的線平均流速 vl 和體積流量 Q 的表達式：

式中： L —超聲波(bō)流量計兩個探頭(tóu)之間的距離； D —管道直徑； vm —管(guan)道的面平均流速(sù)，流速修正系數 K 将聲道線上的(de)速度 vl 修正爲(wei)截面上流體的平(píng)均速度 vm 。
1.2 二次流誤差形成(chéng)原因
流體流(liu)經彎管，管内流體(ti)受到離心力和粘(zhan)性力相互作用，在(zài)管道徑向截面上(shàng)形成一對反向對(duì)稱渦旋如圖 2所示，稱爲彎管二(er)次流。有一無量綱(gang)數，迪恩數 Dn 可(ke)用來表示彎管二(er)次流的強度。當管(guan)道模型固定時，迪(di)恩數 Dn 隻與雷(lei)諾數 Re 有關。研(yan)究發現，流速越大(da)，産生的二次流強(qiáng)度越大，随着流動(dòng)的發展二次流逐(zhu)漸減弱。

式中(zhong)： d —管道直徑， R —彎管的曲率半(ban)徑。彎管下遊形成(cheng)的二次流在徑向(xiang)👅平面的流動，産生(shēng)了彎管二次流的(de)垂直誤差和水平(píng)誤🈲差。聲道線上二(èr)次流速度方向示(shi)意圖如圖 3所(suǒ)示。本研究在聲道(dao)線路徑上取兩個(ge)觀察面 A和 B，如圖 3（ a）所示；聲道線穿過(guò)這兩個二次流面(mian)的位置爲 a和(he) b，如圖 3（ b）所示。可見由于(yú)聲道線穿過截面(mian)上渦的位置不同(tong)，作用在聲道線上(shàng)的二次流速度方(fāng)向也不同，如圖 3（ c）所示。其中(zhong)，徑向平面二次流(liu)速度在水平方向(xiang)（ X 方向）上的分(fèn)速度，方向相反。

由于超聲波流(liu)量計的安裝，聲道(dào)線均在軸向平面(mian)，這導緻系統無法(fa)檢測到與軸向平(ping)面垂直的二次流(liú)垂直分速度（ Y 方向），産生了二次(cì)流的垂直誤差 Ea，得到 Ea 的計(jì)算公式如下：

式中： vf —聲道線(xian)在軸向平面上的(de)速度。
二次流(liu)水平速度（ X 方(fāng)向的分速度）直接(jie)影響了超聲波流(liú)量計的軸向檢測(cè)平面，對檢測造成(chéng)了非常大的影響(xiǎng)。聲道線在空間上(shang)先後收到方向相(xiang)反的二次流水平(ping)速度的作用，這🈲在(zài)很大程度上🤞削弱(ruo)了誤差。但反向速(sù)度并🔴不*相等，且超(chāo)聲波流量計是按(àn)固定角度進行速(su)度折算的，超聲波(bo)傳播速度 vs 對(duì)應地固定爲軸向(xiang)流速爲 vd ，而其(qí)真實流速爲 vf ，由此二次流徑向(xiang)兩個相反的水平(píng)速度，分别導緻了(le) Δv1（如圖 4（ a）所示）和 Δv2（如(rú)圖 4（ b）所示(shì)）兩個速度變化量(liàng)，其中 Δv1 導緻測(ce)得的流速偏大， Δv2 導緻測得的流(liú)速偏小，兩個誤差(chà)不能抵消，産生二(er)🥰次流的水平誤差(cha) Eb ：

式中： vx —聲道線線上 X 方向的分速度(du)即二次流水平速(sù)度， vz —Z 方向的分(fèn)速度即主流方向(xiang)分速度。
三、數(shù)值仿真
2.1 幾何(hé)模型
幾何模(mo)型采用的是管徑(jìng)爲 50 mm的管道，彎(wān)管流場幾何模型(xing)示意圖如圖 5所示。其由上遊緩(huǎn)沖管道、彎管、下遊(yóu)緩沖管道、測量管(guǎn)道、出口管道 5 部分構成。全美氣(qi)體聯合會（ AGA）發(fa)表的 GA-96建議，在(zài)彎管流場的下遊(yóu)保留 5倍管徑(jing)的直管作爲緩沖(chong)，但有研究表明這(zhè)個距離⁉️之🌏後二🎯次(cì)流的作用仍十分(fen)明顯。
據此，筆(bi)者設置流量計的(de) 3個典型安裝(zhuāng)位置來放置測量(liàng)管道，分别距上遊(yóu)彎道爲🎯 5D， 10D， 20D。本研究在彎(wān)管出口處頂部和(he)底部分别設置觀(guān)測點，測量兩點壓(yā)力，得到兩點的壓(ya)力差。
2.2 仿真與(yǔ)設定
在仿真(zhen)前，筆者先對幾何(hé)模型進行網格劃(huà)分。網格劃分采用(yòng) Gambit軟件，劃分時(shí)，順序是由線到面(mian)，由面到體。其中，爲(wèi)了得到更好的收(shou)斂性和精度，面網(wǎng)格如圖 6所示(shì)。其采用錢币畫法(fa)得到的矩形網格(ge)，體網格如圖 7所示。其在彎道處(chu)加深了密度。網格(ge)數量總計爲 1.53×106。畫好網格後，導入(ru) Fluent軟件進行計(ji)算，進口條件設爲(wei)速度進口，出口設(shè)爲🐆 outflow，介質爲空(kong)氣。研究結果表明(míng)，湍流模型采用 RSM時與真實測量(liang)zui接近［ 8］，故本研(yan)究選擇 RSM模型(xíng)。
爲了排除次(cì)要因素的幹擾，将(jiang)仿真更加合理化(huà)，本📞研究進行如下(xia)設定： ①幾何模(mo)型固定不變，聲波(bō)發射角度設置爲(wèi) 45°； ②結合流(liu)量計的實際量程(chéng)，将雷諾數（ Re）設(shè)置爲從 3000~50000，通過(guo)改變進口速度，來(lai)研究 Re 對測量(liàng)精度的影響； ③由于 Fluent是無法(fǎ)将聲波的傳播時(shí)間引入的，對于聲(shēng)道線上的速🍉度❤️，筆(bi)者采用提取聲道(dao)線每個節點上的(de)速度🧡，然後進行線(xian)積分的方法💛計算(suàn)。
四、仿真結果(guǒ)分析與讨論
3.1 誤差分析與讨論(lùn)
彎管下遊緩(huan)沖管道各典型位(wèi)置（ 5D， 10D， 20D）二次流垂直誤差(cha)如圖 8（ a）所(suǒ)示，當下遊緩沖管(guǎn)道爲 5D時，二次(cì)流垂直誤差基本(ben)可以分爲兩個階(jie)段，起初，誤差👉随着(zhe) Re 的增大而增(zēng)大，在 Re 值 13 000之前，增幅明顯，當(dāng) Re 值在 13 000~16 000時(shí)，增幅趨于平緩。在(zai)經過 Re 值 16 000這個後，誤差反而(er)随着 Re 值的增(zēng)大而減小。當下遊(yóu)緩沖管道爲 10D 時，誤差總體上随(sui)着 Re 的增大而(ér)增大，在 Re 值 14 000之前處于增幅(fu)明顯的上升趨勢(shì)，從 Re 值 14 000之(zhi)後增幅開始減小(xiǎo)。下遊緩沖管道爲(wei) 20D 時，誤差随 Re 值增大而增大(dà)，增幅緩慢，且并不(bu)十分穩定，這是由(you)于二次流✉️在流經(jīng) 20D時，已經發生(sheng)衰減，二次流狀态(tai)不是很穩定。二次(ci)流水平誤差如圖(tu) 8（ b）所示，其(qí)非常顯著的特點(diǎn)是誤差出現了正(zhèng)、負不同的👉情況， 10D 處由于 Δv1 比(bǐ) Δv2 要小，測得的(de)流速偏小，誤差值(zhi)變爲負，而在 5D 和 20D 處， Δv1和(hé) Δv2 的大小關系(xì)正好相反，流速偏(piān)大，誤差值爲正，這(zhe)表明二次流的水(shuǐ)平誤差跟安裝位(wei)置有很大關系，甚(shen)至出現了⭐誤差正(zhèng)、負🈲不同的情況。
對比不同下遊(you)緩沖管道，總體看(kàn)來，随着流動的發(fā)展，二💁次流強度減(jian)弱，誤差減小。但在(zài) Re 值 29 000之前(qián)， 5D 處的二次流(liú)垂直誤差比 10D 處大，在 Re 值 29 000之後，由于變化(huà)趨勢不同， 10D 處(chu)的誤差超過了 5D 處的誤差。可見(jiàn)，并不是距離上遊(yóu)彎管越近，誤差就(jiù)越大。對比兩種誤(wù)差可見，二次流的(de)垂直誤差總體🛀大(da)于二次流的水平(ping)誤差🙇‍♀️。
3.2 誤差修(xiū)正
實際測量(liàng)場合下，流量計本(běn)身就是測量流速(sù)的，所以事先并不(bú)知道彎管下遊的(de)二次流強度，這導(dao)緻研究人員在知(zhi)道誤🚶差規♋律的情(qíng)況下無法得知實(shi)際誤差。針對該🥰情(qíng)況，結合流體經過(guò)彎管後的特點，本(ben)研究在流體彎管(guǎn)出口處的頂端和(he)底端各設置一壓(ya)力測試㊙️點，得到其(qi)出口處的壓力差(chà)以反映二次流的(de)強度。雷諾數與彎(wān)管出口壓力如圖(tu) 9所示。由圖 9可見，壓力差随(suí)着雷諾數的增大(da)而增大，在實際安(ān)裝場合，管道模型(xing)固定，由此，壓力差(cha)可用來反映二次(cì)流的強度。将🔞雷諾(nuo)數用壓力差表示(shi)，得到壓力差跟二(er)次流的垂直誤差(cha)✨和水平誤差的關(guān)系。将兩種誤差結(jie)合，可得二次流的(de)總誤差 E總：
E總 =Ea Eb -Ea ×Eb （ 9）
壓力差與總誤(wu)差關系圖如圖 10所示。zui終通過壓(ya)力差來對彎管二(er)次流誤差進行修(xiū)正，得出壓力差與(yu)修正系數關系圖(tú)。