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

關鍵字：   超聲(shēng)波流量計   測(ce)量過程中   彎(wān)管誤差

一、本(ben)文引言 　 　 　

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

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

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

式中： L —超聲波(bo)流量計兩個探頭(tou)之間的距離； D —管道直徑； vm —管(guǎn)道的面平均流速(sù)，流速修正系數 K 将聲道線上的(de)速度 vl 修正爲(wèi)截面上流體的平(ping)均速度 vm 。
1.2 二次流誤差形成(cheng)原因
流體流(liú)經彎管，管内流體(ti)受到離心力和粘(zhan)性力相互作用，在(zai)管道徑向截面上(shàng)形成一對反向對(dui)稱渦旋如圖 2所示，稱爲彎管二(er)次流。有一無量綱(gāng)數，迪恩數 Dn 可(ke)用來表示彎管二(èr)次流的強度。當管(guan)道模型固定時，迪(dí)🌈恩數 Dn 隻與雷(léi)諾數 Re 有關。研(yán)究發現，流速越大(dà)，産生的二次流強(qiang)度越大，随着流動(dong)的發展二次流逐(zhú)漸減弱。

式中(zhong)： d —管道直徑， R —彎管的曲率半(ban)徑。彎管下遊形成(cheng)的二次流在徑向(xiàng)平面的流動，産生(shēng)了彎管二次流的(de)垂直誤差和水平(ping)誤差。聲道線🥰上二(er)次流速度方向示(shi)意圖如圖 3所(suo)示。本研究在聲道(dao)線路徑上取兩個(gè)觀察面 A和 B，如圖 3（ a）所示；聲道線穿過(guò)這兩個二次流面(mian)的位置爲 a和(hé) b，如圖 3（ b）所示。可見由于(yu)聲道線穿過截面(miàn)上渦的位置不同(tóng)，作用在聲道線上(shang)的二次流速度方(fang)向也不同，如圖 3（ c）所示。其中(zhong)，徑向平面二次流(liú)速度在水平方向(xiang)（ X 方向）上的分(fen)速度，方向相反。

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

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

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