摘要;超聲流量計(jì)測(cè)量過(guò)程探頭大小和結(jié)構(gòu)設(shè)計(jì)所中探頭對(duì)流場(chǎng)的干擾是流量計(jì)流聲耦合仿真,定量分析了引起的執(zhí)流效應(yīng)、計(jì)算實(shí)檢量計(jì)探頭擾流的系統(tǒng)偏差;了聲道速度分布、探失聲壓1系統(tǒng)偏差�！眎并利用分段加極平均的方式，這一少靠群進(jìn)更長(zhǎng)聲道長(zhǎng)度情況下的探共挽流系統(tǒng)偏差。
超聲流量計(jì)由于其無(wú)壓損、易安裝、精度高等優(yōu)點(diǎn),近年來(lái)得到了廣泛應(yīng)用。超聲流量計(jì)是通過(guò)測(cè)量超聲波在流體中順流和逆流的時(shí)間差,計(jì)算聲道上的平均流速,再對(duì)不同聲道高度的平均流速進(jìn)行積分,從而求得流量"。常見(jiàn)的探頭安裝方式如圖1所示，在探頭附近容易產(chǎn)生旋渦,影響了時(shí)差的測(cè)量。
利用CFD方法和實(shí)流.實(shí)驗(yàn)研究了不同探頭插人深度時(shí)超聲流量測(cè)量的偏差。為了分析系統(tǒng)偏差的來(lái)源,Loland等利用PIV、LDV和CFD研究了探頭空腔內(nèi)的局部流動(dòng)結(jié)構(gòu);對(duì)探頭空腔內(nèi)的流動(dòng)也進(jìn)行了細(xì)致的實(shí)驗(yàn)研究。兩人的研究關(guān)注點(diǎn)在于流場(chǎng),實(shí)際上流量計(jì).測(cè)到的聲波信號(hào)里體現(xiàn)了波束范圍內(nèi)流動(dòng)的影響和壁面反射對(duì)聲波信號(hào)的干擾，流場(chǎng)和聲場(chǎng)兩者耦合作用共同造成了流量測(cè)量的偏差。
為了研究超聲流量計(jì)探頭擾流影響的機(jī)理,合理修正探頭擾流影響造成的系統(tǒng)偏差,利用多物理場(chǎng)仿真軟件對(duì)圖1(a)中的管道模型進(jìn)行了流聲耦合仿真,分析了管道探頭模型中的流場(chǎng)細(xì)節(jié)和超聲波耦合傳播方式，并通過(guò)互相關(guān)算法計(jì)算時(shí)差，探討了探頭擾流和壁面反射作用對(duì)流量測(cè)量的影響。
1計(jì)算模型
仿真計(jì)算采用多物理場(chǎng)建模軟件COMSOL。首先進(jìn)行流場(chǎng)仿真,計(jì)算模型為帶有--對(duì)超聲探頭安裝孔的管道,長(zhǎng)度為300mm,直徑為70mm，探頭安裝孔直徑為14mm,按照45°聲道角分布于管道兩側(cè)，如圖2所示。流場(chǎng)仿真采用不可壓縮流動(dòng)k-&湍流模型來(lái)模擬管道中流場(chǎng)的流動(dòng)過(guò)程,并用PARDISO算法進(jìn)行穩(wěn)態(tài)求解，管道平均流速為3m/s。
式中ƒ0為振動(dòng)頻率,A為振動(dòng)幅值。假設(shè)理想介質(zhì)水域?yàn)檫B續(xù)介質(zhì),聲波在水域中的能量損耗為零,利用聲波在流體中的連續(xù)性方程,并通過(guò)MUMPS算法進(jìn)行瞬態(tài)求解,對(duì)聲波在水流中傳播的方式進(jìn)行仿真，
式中,P為聲壓;P0為流壓力;po為流密度;c0為聲速;V0為流速。本文中聲速C0設(shè)置為1481m/s。分別在探頭A探頭B端添加式(1)振動(dòng)速度u,為減少計(jì)算量,設(shè)置了較低的振動(dòng)頻率(0.2MHz)。
2仿真結(jié)果分析
2.1流場(chǎng)仿真結(jié)果
計(jì)算得到的探頭處流場(chǎng)如圖3所示,探頭附近存在旋渦。把探頭端面分別分為5個(gè)區(qū)域,以5個(gè)區(qū)域的中點(diǎn)M、U、D、L、R作為計(jì)算依據(jù),分別提取5條連線上的流速分布，比較不同區(qū)域的流速變化，如圖4所示。各個(gè)點(diǎn)與中心M的距離為3.9mm。
圖4中橫坐標(biāo)表示聲道方向探頭面與聲道中心的距離;縱坐標(biāo)表示聲道方向的流速,A至B方向流速為正;R區(qū)域與L區(qū)域流速分布相同;Ref是指參考位置即未受到探頭擾流處,壁面連線之間的區(qū)域。探頭A的D區(qū)域和探頭B的U區(qū)域流場(chǎng)有明顯的速度變化，這是因?yàn)樵谔筋^安裝孔處形成了旋渦,流速在這兩個(gè)區(qū)域內(nèi)變化最為劇烈,而且相對(duì)流場(chǎng)而言，安裝孔內(nèi)的旋渦方向和大小并不相同;M區(qū)域和L區(qū)域處的流場(chǎng)相類似,受旋渦影響較小。
2.2聲場(chǎng)仿真結(jié)果
圖5展示了探頭A發(fā)射超聲波時(shí),超聲波的傳播過(guò)程。在探頭A、探頭B的壁面處,超聲波發(fā)生了反射,反射信號(hào)和原信號(hào)相互疊加,造成了接收面聲壓的不對(duì)稱,進(jìn)而影響傳播時(shí)間的測(cè)量。探頭B接收聲壓的分布情況如圖6所示,聲壓在接收面上非均勻分布,在接收面上分布-一個(gè)低壓區(qū),低壓區(qū)中心位于中心偏下游的位置。
圖7為超聲傳播過(guò)程中不同位置聲壓振幅的分布情況,其中位置1~位置5已在圖5(a)中標(biāo)注,統(tǒng)計(jì)的是聲波經(jīng)過(guò)該位置的聲壓變化的振幅。在探頭A處的壁面反射造成了發(fā)射聲壓分布的畸變,在傳輸過(guò)程中逐漸減少了它的影響,所以探頭B所接收的聲壓主要受到B處壁面反射的影響,在探頭B附近低壓中心從上游逐漸向下游移動(dòng)。
2.3傳播時(shí)間及流速計(jì)算
由于旋渦和璧面反射的影響,靠近探頭邊緣區(qū)域的聲壓曲線存在一定畸變。探頭B不同區(qū)域接收的聲壓與平均聲壓的關(guān)系如圖8所示。其中實(shí)線代表平均聲壓,虛線代表M區(qū)域處的聲壓變化。聲波傳播過(guò)程中,受到不同聲波傳播路徑和壁面反射的影響,接收面不同位置,接收聲壓幅值與過(guò)零點(diǎn)有明顯的區(qū)別。M區(qū)域處聲壓曲線幅值略高于平均聲壓曲線、過(guò)零點(diǎn)與平均聲壓曲線接近;U區(qū)域和D區(qū)域處聲壓曲線過(guò)零點(diǎn)與平均曲線有較大差異。
由于探頭不同區(qū)域聲壓變化曲線的差異,采用平均聲壓曲線來(lái)計(jì)算時(shí)間差,平均聲壓的計(jì)算結(jié)果接近聲壓中心,而且有更好的穩(wěn)定性。探頭A、探頭B接收到的平均聲壓變化曲線如圖9所示,探頭A由于流體的減速作用收到波形略晚，兩個(gè)波形的相似度較高。利用互相關(guān)函數(shù)計(jì)算時(shí)差:
式中,y1(m)和y2(m)為探頭A、探頭B接收聲壓信號(hào);m為數(shù)據(jù)長(zhǎng)度,由互相關(guān)理論，當(dāng)互相關(guān)函數(shù)取得最大值的時(shí)間位移,對(duì)應(yīng)的是兩波形之間的時(shí)差。對(duì)R(m)進(jìn)行優(yōu)化求解,假設(shè)在m0點(diǎn)處取得最大值,可以求得時(shí)差△t:
實(shí)際流量計(jì)測(cè)量時(shí),通常是通過(guò)正逆向傳播時(shí)間T1、t2去和時(shí)差△t計(jì)算流速,由于流速遠(yuǎn)小于聲速u0<<C0,可以進(jìn)一步得到:
式中,L為聲道長(zhǎng)度。將△t代人到式(5)中,可以解得聲道方向的平均流速V'm。聲場(chǎng)仿真計(jì)算中,△t=1.986x10-7s,L=113.0mm,求得V'm=I.928m/s。
對(duì)比管道流場(chǎng)計(jì)算結(jié)果,對(duì)聲束范圍內(nèi)流速取平均值,求得V'm=1.934m/s,未受到探頭擾流區(qū)的聲道方向平均流速Vm=2.247m/s,求得聲場(chǎng)和流場(chǎng)計(jì)算的系統(tǒng)偏差E分別為一14.2%和-12.5%。兩者的差異體現(xiàn)了壁面反射對(duì)修正系數(shù)的影響。
3推論和討論
超聲探頭對(duì)流場(chǎng)的擾動(dòng)通常只發(fā)生在探頭附近--定范圍內(nèi),該范圍之外流場(chǎng)與上游充分發(fā)展的流場(chǎng)相同，因此可以用加權(quán)平均的方式將第2節(jié)中的計(jì)算結(jié)果向更長(zhǎng)的聲道進(jìn)行推論。在圖10所示的探頭安裝方式下,聲道角度為φ,探頭直徑為D。假設(shè)探頭在一定范圍內(nèi)影響流場(chǎng),上下游流場(chǎng)受影響的范圍為b,所以將流場(chǎng)沿聲道方向劃分成3個(gè)區(qū)域,分別為兩端的流場(chǎng)受影響速度區(qū)和中間的非影響速度區(qū)。
通過(guò)流場(chǎng)和聲場(chǎng)耦合仿真計(jì)算3個(gè)區(qū)域內(nèi)平均投影速度,然后用加權(quán)分析的方法計(jì)算凸出效應(yīng)造成的系統(tǒng)偏差,如式(7)所示。
式中,V1、V2分別為未受到探頭擾流處上下游受影響速度區(qū),聲道方向的平均流速;V'1、V'2分別為探頭擾流處,上下游受影響速度區(qū),聲道方向的平均流速;V為非影響速度區(qū),聲道方向的平均流速。其中幾何尺寸L、D、φ為固定值,V1、V2可以正確計(jì)算,所以確認(rèn)b的范圍和V'1、V'2的值是確定修正系數(shù)的關(guān)鍵。
通過(guò)對(duì)圖4分析可知,在這種安裝方式下,在管道中心位置附近,各方向的流速是相近的,流速差值小于1%,可以認(rèn)為管道中心附近為非影響速度區(qū);在管道中心兩側(cè),不同區(qū)域的流速變化情況不同,可以認(rèn)為受影響速度區(qū)的范圍b=3.535D。通過(guò)耦合仿真計(jì)算,將非影響區(qū)域的流速平均值V=2.458m/s和平均聲壓計(jì)算流速值V"m=1.928m/s帶人式(7)中,求得上下游影響區(qū)域內(nèi)的平均流速(V'1+V"2)/2為1.928m/s。再利用式(7)進(jìn)行加權(quán)計(jì)算,可以推算出更長(zhǎng)聲道時(shí)的探頭擾流影響,設(shè)定非影響區(qū)域的流速為1,求得不同管道口徑下,流速的系統(tǒng)偏差E如表1所示,其中聲道角度φ=45°,探頭直徑D=14mm。
4結(jié)論
超聲流量計(jì)探頭局部結(jié)構(gòu)帶來(lái)的擾流效應(yīng)造成了其流量測(cè)量的系統(tǒng)偏差,這一偏差通常利用實(shí)驗(yàn)室實(shí)流校準(zhǔn)來(lái)修正。為了更好地分析探頭擾流影響機(jī)理,利用多物理場(chǎng)仿真軟件對(duì)其進(jìn)行了流聲耦合分析,主要結(jié)論如下:
①探頭凹坑內(nèi)存在低速區(qū)且有漩渦,聲束范圍內(nèi)各區(qū)域的平均流速與探頭中心區(qū)域上的平均流速不同,再加上探頭附近的壁面聲波反射,造成探頭端面不同區(qū)域接收到的聲壓信號(hào)有差異,流量計(jì)測(cè)到的聲波傳播時(shí)間體現(xiàn)的是聲壓信號(hào)統(tǒng)計(jì)平均的結(jié)果。
②對(duì)于帶有直徑14mm的斜插縮進(jìn)式探頭的DN70流量計(jì),按照探頭收到的面平均聲壓信號(hào)計(jì)算時(shí)差,探頭擾流造成的系統(tǒng)偏差約為-14.2%。
③在仿真結(jié)果的基礎(chǔ)上,假設(shè)探頭擾流影響范圍只限于其附近一定范圍,利用分段加權(quán)平均的方式,推導(dǎo)了更長(zhǎng)的聲道長(zhǎng)度情況下的探頭擾流系統(tǒng)偏差,發(fā)現(xiàn)該偏差均為負(fù)偏差,其絕對(duì)值近似等于探頭縮進(jìn)比，隨著聲道長(zhǎng)度的增加而降低。
本文來(lái)源于網(wǎng)絡(luò),如有侵權(quán)聯(lián)系即刪除！
