摘要:以DN50渦輪流量計為例，采用CFD方法模擬了方波和正弦波兩種脈動(dòng)流在渦輪流量計中的流動(dòng)。發(fā)現在不同波形、振幅、頻率下的脈動(dòng)流對渦輪流量計的影響具有差異性。各工況下差異性的來(lái)源，從流場(chǎng)速度和葉輪所受剪切力，分析渦輪流量計在不同工況的脈動(dòng)流下內部流場(chǎng)的變化。發(fā)現兩種脈動(dòng)流在同頻率同振幅的波峰時(shí)，對比正弦波脈動(dòng)流，方波脈動(dòng)流下葉輪周?chē)牧鲌?chǎng)速度更快，葉輪所受到的剪切應力更大，持續時(shí)間長(cháng)，力矩峰值更大，這是造成方波誤差大的根本原因。
0引言
天然氣作為清潔高效的綠色能源，已被廣泛應用于人民生活和工業(yè)生產(chǎn)中，天然氣的正確計量是利益相關(guān)方的關(guān)注重點(diǎn)。渦輪流量計在使用中具有精度高、重復性好、壓損小等優(yōu)點(diǎn)。但是，在天然氣大流量管輸、天然氣大流量檢定站中，由于流體管輸變化、動(dòng)力機械旋轉、調節閥的動(dòng)作等，都可能產(chǎn)生脈動(dòng)流，會(huì )對渦輪流量計現場(chǎng)應用產(chǎn)生的測量誤差。
目前都聚焦正弦波脈動(dòng)流對渦輪流量計的影響，關(guān)注脈動(dòng)頻率和振幅，但是正弦波脈動(dòng)流與實(shí)際脈動(dòng)流之間存在較大區別，為了更好地實(shí)現脈動(dòng)流修正，討論和分析脈動(dòng)流影響渦輪流量計的關(guān)鍵要素是必要的。
為了獲取脈動(dòng)流影響渦輪流量計的關(guān)鍵要素，利用CFD仿真，對比和分析正弦波以及方波脈動(dòng)流對渦輪流量計影響的差異性，從流場(chǎng)速度以及葉輪所受剪切力兩個(gè)方面來(lái)探究差異性的來(lái)源所在。
1計算流體力學(xué)仿真方法
1.1模型建立及網(wǎng)格劃分
以1.5級DN50氣體渦輪流量計為原型，并對其進(jìn)行仿真。其內部基本結構參數詳見(jiàn)表1。建立仿真模型，對渦輪量計進(jìn)行內部流場(chǎng)的提取。前后直管段同軸安裝，前直管段長(cháng)為10D，后直管段5D。前后直管段的網(wǎng)格尺寸為3mm,葉輪周?chē)木W(wǎng)格尺寸為1mm，導流片前后為1mm,葉輪邊沿的網(wǎng)格尺寸為0.2mm,對結構邊緣及狹縫進(jìn)行加密處理，網(wǎng)格總數為.262萬(wàn)。
1.2邊界條件
入口采用速度入口條件(Velocit-inlet)，出口采用壓力出口條件(Pressure-outlet)，選擇渦輪流.量計的分界流量點(diǎn)進(jìn)行數值仿真。該流量計的分界.流量點(diǎn)為20m3/h,換算平均入口流速為2.83m/s。出口壓力設置為0.5MPa。流體介質(zhì)為空氣，密度為1.225kg/m3，動(dòng)力黏度為1.79X10-5kg/(m.s)。選擇的湍流模型為Realizablek-e模型。運動(dòng)模型采用6DOF模型。
入口速度表達式如下所示：
V1=2.83+asin(2πƒt)(1)
V2=2.83+α(-1)[2ft](2)
式中，V1，,V2-----分別為正弦波、方波速度入口的瞬時(shí)速度，m/s。
1.3儀表系數獲取方法及結果
仿真結束后，在葉輪的輪轂上取一個(gè)點(diǎn)P，通過(guò)后處理軟件，獲得葉輪旋轉穩定后的平均線(xiàn)速度速度vp,通過(guò)公式(3)可計算葉輪的轉速。
式中。n一葉輪轉速，單位:r/s;rp一點(diǎn)P繞X軸旋轉的半徑。
再由轉速n計算渦輪流量計的儀表系數K。折算后的儀表系數結果如表2所示，脈動(dòng)流下的渦輪流量計存在著(zhù)正誤差，相同振幅、頻率下，不同波形脈動(dòng)流對渦輪流量計的影響依然存在差異性。
2流場(chǎng)速度分布
以5Hz不同波形脈動(dòng)流作用的渦輪流量計流場(chǎng)對比分析。圖1所示為方波脈動(dòng)流在四種振幅下對應穩定后一個(gè)周期T內葉輪周?chē)牧鲌?chǎng)速度分布圖，從左到右時(shí)間點(diǎn)依次為0T、0.25T、0.5T、0.75T。在方波工況下，0T、0.25T時(shí)的流場(chǎng)速度明顯比0.5T、0.75T時(shí)的流場(chǎng)速度更快，隨著(zhù)振幅的加大，0T.0.25T時(shí)的流場(chǎng)速度增大，而0.5T.0.75T時(shí)的流場(chǎng)速度則減慢。
對比兩種波形的速度分布圖，同振幅下的波峰時(shí)，方波工況下流場(chǎng)速度略大于正弦波工況，且方波的高流速持續時(shí)間比正弦波的高流速持續時(shí)間長(cháng)。
3葉輪剪切應力分布
圖3~圖4分別是5Hz時(shí)方波、正弦波脈動(dòng)流下葉輪，所受力矩。除了流量突變點(diǎn)外，葉輪所受力矩與流量波形近似。
以5Hz不同波形脈動(dòng)流作用的渦輪流量計葉輪所受剪切力對比分析。方波脈動(dòng)流在四種振幅下穩定后一個(gè)周期T內，分析葉輪所受到的剪切應力分布。葉輪的剪切應力在0T、0.25T達到最大，隨著(zhù)振幅增大而增大;在0.5T、0.75T時(shí)達到最小，隨著(zhù)振幅的增大而減小。正弦波脈動(dòng)流在四種振幅下穩定后一個(gè)周期T內，分析葉輪所受到的剪切應力分布。在0T、0.5T時(shí)，葉輪所受應力幾乎不隨振幅變化，在0.25T時(shí)即在波峰時(shí)的應力隨著(zhù)振幅加大而增加,在0.75T時(shí)隨著(zhù)振幅增大而減小。
4總結與討論
對口徑為DN50的氣體渦輪流量計仿真建模，進(jìn)行了一系列仿真實(shí)驗，以入口脈動(dòng)流的波形、頻率、振幅為變量條件，完成了相關(guān)仿真實(shí)驗。通過(guò)仿真數據分析可知:
(1)脈動(dòng)流作用下，渦輪流量計會(huì )出現正誤差。由流場(chǎng)速度分布、剪切應力分布及葉輪所受到的力矩可以看出，葉輪轉速相對穩定后，脈動(dòng)作用過(guò)程中葉輪加速和減速時(shí)間幾乎相同，而相同流量的脈動(dòng)流下的葉輪轉速比定常流下的高，表明加速過(guò)程響應更快。即葉輪在加速過(guò)程中更為敏感，在相同流量升降情況下，流體流速上升時(shí)，葉輪轉速立馬上升，而流速下降時(shí)，葉輪轉速下降較緩慢或者延遲下降，從而使得脈動(dòng)流下葉輪轉速高于定常流下的葉輪轉速，始終產(chǎn)生正誤差。
(2)相同頻率振幅下，不同波形的脈動(dòng)流產(chǎn)生的誤差大小不同。方波脈動(dòng)流對應的誤差最大，其次是正弦波。通過(guò)速度分布圖可以看出，在相同頻率振幅的波峰時(shí)，方波脈動(dòng)流下葉輪周?chē)牧鲌?chǎng)速度和應力都要比其他兩種波形的脈動(dòng)流更大。在方波作用下，一個(gè)周期內葉輪被加速的時(shí)間更長(cháng)，導致了相同頻率振幅下，方波工況下葉輪轉速比正弦波轉速更高。正弦波脈動(dòng)下，葉輪轉速是漸變的，從而相比方波脈動(dòng)的轉速略低。
(3)在方波脈動(dòng)流工況下，葉輪受的到?jīng)_擊力更強，使葉輪有遭受物理破壞的可能性。針對脈動(dòng)流對渦輪流量計的影響，應該還是盡量從振源避免脈動(dòng)流;在不能避免脈動(dòng)流時(shí)，應遠離振源，或者是增加流動(dòng)調整器，降低脈動(dòng)頻率、振幅，盡量消除類(lèi)似脈動(dòng)突增，保護渦輪流量計、控制測量誤差;智能渦輪流量計，可以考慮在流量計的積算模塊增加脈動(dòng)流修正公式，實(shí)現對渦輪流量計脈動(dòng)流誤差的修正補償，保障流量計在有脈動(dòng)流工況下的測量精度。
文章來(lái)源于網(wǎng)絡(luò ),如有侵權聯(lián)系即刪除！
